Testing Gassmann fluid substitution: sonic logs versus ultrasonic core measurements

Although Gassmann fluid substitution is standard practice for time-lapse studies, its validity in the field environment rests upon a number of underlying assumptions. The impact of violation on the predictions of Gassmann equations can only ultimately be validated by in situ testing in real geological environments. In this paper we show a workflow that we developed to test Gassmann fluid substitution by comparing saturated P-wave moduli computed from dry core measurements against those obtained from sonic and density logs. The workflow has been tested on 43 samples taken from a 45 m turbidite reservoir from the Campos Basin, offshore Brazil. The results show good statistical agreement between the P-wave elastic moduli computed from cores using the Gassmann equation with the corresponding moduli computed from log data. This confirms that all the assumptions of the Gassmann are adequate within the measurement error and natural variability of elastic properties. These results provide further justification for using the Gassmann theory to interpret time-lapse effects in this sandstone reservoir and in similar geological formations.     

MODFLOW 2000 head uncertainty,a first-order second moment method

A computationally efficient method to estimate the variance and covariance in piezometric head results computed through MODFLOW 2000 using a first-order second moment (FOSM) approach is presented. This methodology employs a first-order Taylor series expansion to combine model sensitivity with uncertainty in geologic data. MODFLOW 2000 is used to calculate both the ground water head and the sensitivity of head to changes in input data. From a limited number of samples, geologic data are extrapolated and their associated uncertainties are computed through a conditional probability calculation. Combining the spatially related sensitivity and input uncertainty produces the variance-covariance matrix, the diagonal of which is used to yield the standard deviation in MODFLOW 2000 head. The variance in piezometric head can be used for calibrating the model, estimating confidence intervals, directing exploration, and evaluating the reliability of a design. A case study illustrates the approach, where aquifer transmissivity is the spatially related uncertain geologic input data. The FOSM methodology is shown to be applicable for calculating output uncertainty for (1) spatially related input and output data, and (2) multiple input parameters (transmissivity and recharge).     

A Taylor basis for kinematic nonlinear real-time simulations. Part II: The Taylor basis

Real-time simulations are used to a significant extent in many engineering fields. However, if nonlinearities are included, the real-time requirement significantly limits the size and complexity of numerical models. The present work constitutes the second of two papers where a general basis method to simulate kinematic nonlinear structures more efficiently is introduced. The advantage of the basis formulation is that it enables the number of basis vectors to be increased without increasing the number of unknown basis co-ordinates. This allows for larger numerical kinematically nonlinear models to run in real time. The basis is organized from a Taylor series that includes the system mode shapes and their complete first-order modal derivatives derived in Part I. The Taylor series predicts fixed linear relations between the modal co-ordinates of the system mode shapes and the modal derivatives, respectively. Thus, the full solution is known solely by determining the modal co-ordinates of the mode shapes, which significantly minimizes the computational costs. Furthermore, it is illustrated that the stability of the Taylor basis formulation is dependent on the mode shape frequencies only, allowing the applied time steps to be significantly larger than in standard nonlinear basis analysis. An example illustrates a case where the computational time can be decreased by one order of magnitude using a Taylor basis formulation compared with a standard basis formulation including identical basis vectors.     

用于声波方程数值模拟的时间-空间域有限差分系数确定新方法

声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.为精确求解空间偏导数,相继发展了高阶差分格式优化方法和伪谱方法.近期,为更好地缓解数值频散,提出了时间-空间域有限差分方法,该方法采用了泰勒展开近似方法来确定有限差分格式系数,因而只能保证在一定的小范围内很好的拟合波场传播规律.为进一步压制数值频散效应,本文引入了时间-空间域特定波数点满足频散关系的方法,根据震源、波速和网格间距确定波数范围,同时考虑了多个传播角度,然后建立方程确定了相应的有限差分格式系数,使得差分系数能在更大范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性.     

Elastic moduli of dry rocks containing spheroidal pores based on differential effective medium theory

Differential effective medium (DEM) theory is applied to determine the elastic properties of dry rock with spheroidal pores. These pores are assumed to be randomly oriented. The ordinary differential equations for bulk and shear moduli are coupled and it is more difficult to obtain accurate analytical formulae about the moduli of dry porous rock. In this paper, we derive analytical solutions of the bulk and shear moduli for dry rock from the differential equations by applying an analytical approximation for dry-rock modulus ratio, in order to decouple these equations. Then, the validity of this analytical approximation is tested by integrating the full DEM equation numerically. The analytical formulae give good estimates of the numerical results over the whole porosity range. These analytical formulae can be further simplified under the assumption of small porosity. The simplified formulae for spherical pores (i.e., the pore aspect ratio is equal to 1) are the same as Mackenzie's equations. The analytical formulae are relatively easy to analyze the relationship between the elastic moduli and porosity or pore shapes, and can be used to invert some rock parameters such as porosity or pore aspect ratio. The predictions of the analytical formula for the sandstone experimental data show that the analytical formulae can accurately predict the variations of elastic moduli with porosity for dry sandstones. The effective elastic moduli of these sandstones can be reasonably well characterized by spheroidal pores, whose pore aspect ratio was determined by minimizing the error between theoretical predictions and experimental measurements. For sandstones the pore aspect ratio increases as porosity increases if the porosity is less than 0.15, whereas the pore aspect ratio remains relatively stable (about 0.14) if the porosity is more than 0.15.     

求取岩石基质体积模量的线形拟合方法(英文)

岩石基质的体积模量或其倒数--压缩系数,在进行油气预测的流体替换和孔隙度反演时,是重要的输入参数,但是利用现有方法很难准确求得。文中提出了一种求取该参数的线形拟合方法,该方法通过对Gassmann方程的合理简化并引入Eshelby-Walsh干燥岩石椭球包体近似公式,获得了计算岩石基质压缩系数的拟合公式,可方便地利用公式计算该参数。实际碳酸盐岩岩样的岩石物理测试分析显示:利用饱和岩样和干燥岩样测得的基质压缩系数的差异小于1%,说明所求参数是正确的,可靠的。     

Analytic expressions for the velocity sensitivity to the elastic moduli for the most general anisotropic media

For non‐linear kinematic inversion of elastic anisotropy parameters and related investigations of the sensitivity of seismic data, the derivatives of the wavespeed (phase velocity and group velocity) with respect to the individual elastic moduli are required. This paper presents two analytic methods, called the eigenvalue and eigenvector methods, to compute the derivatives of the wavespeeds for wave propagation in a general anisotropic medium, which may be defined by up to 21 density‐normalized elastic moduli. The first method employs a simple and compact form of the eigenvalue (phase velocity) and a general form of the group velocity, and directly yields general expressions of the derivatives for the three wave modes (qP, qS₁, qS₂). The second method applies simple eigenvector solutions of the three wave modes and leads to other general forms of the derivatives. These analytic formulae show that the derivatives are, in general, functions of the 21 elastic moduli as well as the wave propagation direction, and they reflect the sensitivity of the wavespeeds to the individual elastic moduli. Meanwhile, we give results of numerical investigations with some examples for particular simplified forms of anisotropy. They show that the eigenvalue method is suitable for the qP‐, qS₁‐ and qS₂‐wave computations and mitigates the singularity problem for the two quasi‐shear waves. The eigenvector method is preferable to the eigenvalue method for the group velocity and the derivative of the phase velocity because it involves simpler expressions and independent computations, but for the derivative of the group velocity the derivative of the eigenvector is required. Both methods tackle the singularity problem and are applicable to any degree of seismic anisotropy for all three wave modes.     

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.     

A semi-empirical velocity-porosity-clay model for petrophysical interpretation of P- and S-velocities

We design a velocity–porosity model for sand-shale environments with the emphasis on its application to petrophysical interpretation of compressional and shear velocities. In order to achieve this objective, we extend the velocity–porosity model proposed by Krief et al., to account for the effect of clay content in sandstones, using the published laboratory experiments on rocks and well log data in a wide range of porosities and clay contents. The model of Krief et al. works well for clean compacted rocks. It assumes that compressional and shear velocities in a porous fluid-saturated rock obey Gassmann formulae with the Biot compliance coefficient. In order to use this model for clay-rich rocks, we assume that the bulk and shear moduli of the grain material, and the dependence of the compliance on porosity, are functions of the clay content. Statistical analysis of published laboratory data shows that the moduli of the matrix grain material are best defined by low Hashin–Shtrikman bounds. The parameters of the model include the bulk and shear moduli of the sand and clay mineral components as well as coefficients which define the dependence of the bulk and shear compliance on porosity and clay content. The constants of the model are determined by a multivariate non-linear regression fit for P- and S-velocities as functions of porosity and clay content using the data acquired in the area of interest. In order to demonstrate the potential application of the proposed model to petrophysical interpretation, we design an inversion procedure, which allows us to estimate porosity, saturation and/or clay content from compressional and shear velocities. Testing of the model on laboratory data and a set of well logs from Carnarvon Basin, Australia, shows good agreement between predictions and measurements. This simple velocity-porosity-clay semi-empirical model could be used for more reliable petrophysical interpretation of compressional and shear velocities obtained from well logs or surface seismic data.     

Estimation of Dry Fracture Weakness,Porosity, and Fluid Modulus Using Observable Seismic Reflection Data in a Gas-Bearing Reservoir

Fracture detection and fluid identification are important tasks for a fractured reservoir characterization. Our goal is to demonstrate a direct approach to utilize azimuthal seismic data to estimate fluid bulk modulus, porosity, and dry fracture weaknesses, which decreases the uncertainty of fluid identification. Combining Gassmann’s (Vier. der Natur. Gesellschaft Zürich 96:1–23, 1951) equations and linear-slip model, we first establish new simplified expressions of stiffness parameters for a gas-bearing saturated fractured rock with low porosity and small fracture density, and then we derive a novel PP-wave reflection coefficient in terms of dry background rock properties (P-wave and S-wave moduli, and density), fracture (dry fracture weaknesses), porosity, and fluid (fluid bulk modulus). A Bayesian Markov chain Monte Carlo nonlinear inversion method is proposed to estimate fluid bulk modulus, porosity, and fracture weaknesses directly from azimuthal seismic data. The inversion method yields reasonable estimates in the case of synthetic data containing a moderate noise and stable results on real data.     

The influence of dependence in characterizing multi-variable uncertainty for climate change impact assessments

Few approaches exist that explicitly use the uncertainty associated with the spread of climate model simulations in assessing climate change impacts. An approach that does so is second-order approximation (SOA). This incorporates quantification of uncertainty to ascertain its impact on the derived response using a Taylor series expansion of the model. This study uses SOA in a statistical downscaling model of monthly streamflow, with a focus on the influence of dependence in the uncertainty of multiple atmospheric variables. Uncertainty is quantified using the square root error variance concept with a new extension that allows the inter-dependence terms among the atmospheric variable uncertainty to be specified. Applying the model to selected point locations in Australia, it is noted that the downscaling results differ considerably from downscaling that ignores uncertainty. However, when the effects of dependence in uncertainty are incorporated, the results differ according to the regional variations in dependence structure.     

Fractional conservation of mass

The traditional conservation of mass equation is derived using a first-order Taylor series to represent flux change in a control volume, which is valid strictly for cases of linear changes in flux through the control volume. We show that using higher-order Taylor series approximations for the mass flux results in mass conservation equations that are intractable. We then show that a fractional Taylor series has the advantage of being able to exactly represent non-linear flux in a control volume with only two terms, analogous to using a first-order traditional Taylor series. We replace the integer-order Taylor series approximation for flux with the fractional-order Taylor series approximation, and remove the restriction that the flux has to be linear, or piece-wise linear, and remove the restriction that the control volume must be infinitesimal. As long as the flux can be approximated by a power-law function, the fractional-order conservation of mass equation will be exact when the fractional order of differentiation matches the flux power-law. There are two important distinctions between the traditional mass conservation, and its fractional equivalent. The first is that the divergence term in the fractional mass conservation equation is the fractional divergence, and the second is the appearance of a scaling term in the fractional conservation of mass equation that may eliminate scale effects in parameters (e.g., hydraulic conductivity) that should be scale-invariant.     

Computation of strain and rotation tensor as well as their uncertainties for small arrays in spherical coordinate system

Based on Taylor series expansion and strain components expressions of elastic mechanics, we derive formulae of strain and rotation tensor for small arrays in spherical coordinates system. By linearization process of the formulae, we also derive expressions of strain components and Euler vector uncertainties respectively for subnets using the law of error propagation. Taking GPS velocity field in Sichuan-Yunnan area as an example, we compute dilation rate and maximum shear strain rate field using the above procedure, and their characteristics are preliminarily car- ried on. Limits of the strain model for small array are also discussed. We make detailed explanations on small array method and the choice of small arrays. How to set weights of GPS observations are further discussed. Moreover relationship between strain and radius of GPS subnets is also analyzed.     

Analytical approximations of bulk and shear moduli for dry rock based on the differential effective medium theory

Differential effective medium theory has been applied to determine the elastic properties of porous media. The ordinary differential equations for bulk and shear moduli are coupled and it is more difficult to obtain accurate analytical formulae about the moduli of dry porous rock. In this paper, in order to decouple these equations we first substitute an analytical approximation for the dry‐rock modulus ratio into the differential equation and derive analytical solutions of the bulk and shear moduli for dry rock with three specific pore shapes: spherical pores, needle‐shaped pores and penny‐shaped cracks. Then, the validity of the analytical approximations is tested by integrating the full differential effective medium equation numerically. The analytical formulae give good estimates of the numerical results over the whole porosity range for the cases of the three given pore shapes. These analytical formulae can be further simplified under the assumption of small porosity. The simplified formulae for spherical pores are the same as Mackenzie's equations. The analytical formulae are relatively easy to analyse the relationship between the elastic moduli and porosity or pore shapes and can be used to invert some rock parameters such as porosity or pore aspect ratio. The predictions of the analytical formulae for experimental data show that the formulae for penny‐shaped cracks are suitable to estimate the elastic properties of micro‐crack rock such as granite, they can be used to estimate the crack aspect ratio while the crack porosity is known and also to estimate the crack porosity evolution with pressure if the crack aspect ratio is given.     

Seismic wave velocities of rare gas solids through elastic properties in Earth’s lower mantle

The expressions for second (SOE) and third order elastic (TOE) constants for rare gas solids are derived for comparative study of elastic behavior within the framework of many body potentials including the effect of pressure. The derived expressions are used to obtain the relations for pressure derivatives of bulk and shear moduli of RGS solids. The values of SOE, TOE constants and pressure derivative of bulk and shear modulus for Ne up to 100 GPa, Ar up to 75 GPa, for Kr up to 136 GPa and Xe up to 53.4 GPa...     

Gassmann Modeling of Acoustic Properties of Sand-clay Mixtures

—The feasibility of modeling elastic properties of a fluid-saturated sand-clay mixture rock is analyzed by assuming that the rock is composed of macroscopic regions of sand and clay. The elastic properties of such a composite rock are computed using two alternative schemes.¶The first scheme, which we call the composite Gassmann (CG) scheme, uses Gassmann equations to compute elastic moduli of the saturated sand and clay from their respective dry moduli. The effective elastic moduli of the fluid-saturated composite rock are then computed by applying one of the mixing laws commonly used to estimate elastic properties of composite materials.¶In the second scheme which we call the Berryman-Milton scheme, the elastic moduli of the dry composite rock matrix are computed from the moduli of dry sand and clay matrices using the same composite mixing law used in the first scheme. Next, the saturated composite rock moduli are computed using the equations of Brown and Korringa, which, together with the expressions for the coefficients derived by Berryman and Milton, provide an extension of Gassmann equations to rocks with a heterogeneous solid matrix.¶For both schemes, the moduli of the dry homogeneous sand and clay matrices are assumed to obey the Krief’s velocity-porosity relationship. As a mixing law we use the self-consistent coherent potential approximation proposed by Berryman.¶The calculated dependence of compressional and shear velocities on porosity and clay content for a given set of parameters using the two schemes depends on the distribution of total porosity between the sand and clay regions. If the distribution of total porosity between sand and clay is relatively uniform, the predictions of the two schemes in the porosity range up to 0.3 are very similar to each other. For higher porosities and medium-to-large clay content the elastic moduli predicted by CG scheme are significantly higher than those predicted by the BM scheme.¶This difference is explained by the fact that the BM model predicts the fully relaxed moduli, wherein the fluid can move freely between sand and clay regions. In contrast, the CG scheme predicts the no-flow or unrelaxed moduli. Our analysis reveals that due to the extremely low permeability of clays, at seismic and higher frequencies the fluid has no time to move between sand and clay regions. Thus, the CG scheme is more appropriate for clay-rich rocks.     

Sensitivity and moment analyses of head in variably saturated regimes

A numerical approach for approximating statistical moments of hydraulic heads of variably saturated flows in multi-dimensional porous media is developed. The approximation relies on a first-order Taylor series expansion of a finite element flow model and an adjoint state numerical method for variably saturated flows to evaluate sensitivities. This approach can be employed to analyze uncertainties associated with predictions of head of steady-state or transient flows in variably saturated porous media, with any type of boundary and initial conditions. Limitations of stochastic analytical methods such as spectral/perturbation approaches and the time-consuming Monte Carlo simulation technique are thus alleviated. An example is given to demonstrate the utility of the approach and to investigate the temporal evolution of head variances in a variably saturated flow regime. Results show that the fluctuation of the water table can have significant impacts on the propagation of the head variance.     

Effect of micro‐inhomogeneity on the effective stress coefficients and undrained bulk modulus of a poroelastic medium: a double spherical shell model

Although most rocks are complex multi‐mineralic aggregates, quantitative interpretation workflows usually ignore this complexity and employ Gassmann equation and effective stress laws that assume a micro‐homogeneous (mono‐mineralic) rock. Even though the Gassmann theory and effective stress concepts have been generalized to micro‐inhomogeneous rocks, they are seldom if at all used in practice because they require a greater number of parameters, which are difficult to measure or infer from data. Furthermore, the magnitude of the effect of micro‐heterogeneity on fluid substitution and on effective stress coefficients is poorly understood. In particular, it is an open question whether deviations of the experimentally measurements of the effective stress coefficients for drained and undrained elastic moduli from theoretical predictions can be explained by the effect of micro‐heterogeneity. In an attempt to bridge this gap, we consider an idealized model of a micro‐inhomogeneous medium: a Hashin assemblage of double spherical shells. Each shell consists of a spherical pore surrounded by two concentric spherical layers of two different isotropic minerals. By analyzing the exact solution of this problem, we show that the results are exactly consistent with the equations of Brown and Korringa (which represent an extension of Gassmann's equation to micro‐inhomogeneous media). We also show that the effective stress coefficients for bulk volume α, for porosity n_? and for drained and undrained moduli are quite sensitive to the degree of heterogeneity (contrast between the moduli of the two mineral components). For instance, while for micro‐homogeneous rocks the theory gives n_? = 1, for strongly micro‐inhomogenous rocks, n_? may span a range of values from –∞ to ∞ (depending on the contrast between moduli of inner and outer shells). Furthermore, the effective stress coefficient for pore volume (Biot–Willis coefficient) α can be smaller than the porosity ?. Further studies are required to understand the applicability of the results to realistic rock geometries.     

Latin hypercube approach to estimate uncertainty in ground water vulnerability

A methodology is proposed to quantify prediction uncertainty associated with ground water vulnerability models that were developed through an approach that coupled multivariate logistic regression with a geographic information system (GIS). This method uses Latin hypercube sampling (LHS) to illustrate the propagation of input error and estimate uncertainty associated with the logistic regression predictions of ground water vulnerability. Central to the proposed method is the assumption that prediction uncertainty in ground water vulnerability models is a function of input error propagation from uncertainty in the estimated logistic regression model coefficients (model error) and the values of explanatory variables represented in the GIS (data error). Input probability distributions that represent both model and data error sources of uncertainty were simultaneously sampled using a Latin hypercube approach with logistic regression calculations of probability of elevated nonpoint source contaminants in ground water. The resulting probability distribution represents the prediction intervals and associated uncertainty of the ground water vulnerability predictions. The method is illustrated through a ground water vulnerability assessment of the High Plains regional aquifer. Results of the LHS simulations reveal significant prediction uncertainties that vary spatially across the regional aquifer. Additionally, the proposed method enables a spatial deconstruction of the prediction uncertainty that can lead to improved prediction of ground water vulnerability.     

Velocity dispersion and fluid substitution in sandstone under partially saturated conditions

The elastic moduli of four sandstone samples are measured at seismic (2?2000 Hz) and ultrasonic (1 MHz) frequencies and water- and glycerin-saturated conditions. We observe that the high-permeability samples under partially water-saturated conditions and the low-permeability samples under partially glycerin-saturated conditions show little dispersion at low frequencies (2?2000 Hz). However, the high-permeability samples under partially glycerin-saturated conditions and the low-permeability samples under partially water-saturated conditions produce strong dispersion in the same frequency range (2?2000 Hz). This suggests that fluid mobility largely controls the pore-fluid movement and pore pressure in a porous medium. High fluid mobility facilitates pore-pressure equilibration either between pores or between heterogeneous regions, resulting in a low-frequency domain where the Gassmann equations are valid. In contrast, low fluid mobility produces pressure gradients even at seismic frequencies, and thus dispersion. The latter shows a systematic shift to lower frequencies with decreasing mobility. Sandstone samples showed variations in V_p as a function of fluid saturation. We explore the applicability of the Gassmann model on sandstone rocks. Two theoretical bounds for the P-velocity are known, the Gassmann–Wood and Gassmann–Hill limits. The observations confirm the effect of wave-induced flow on the transition from the Gassmann–Wood to the Gassmann–Hill limit. With decreasing fluid mobility, the P-velocity at 2–2000 Hz moves from the Gassmann–Wood boundary to the Gassmann–Hill boundary. In addition,, we investigate the mechanisms responsible for this transition.