共查询到20条相似文献，搜索用时 46 毫秒
1.
A critical porosity model establishes the empirical relationship between a grain matrix and a dry rock by the concept of critical porosity. The model is simple and practical and widely used. But the critical porosity in the model is a fixed value that cannot relate to pore structure. The aim of this paper is to establish the theoretical relationship between critical porosity and pore structure by combining Kuster–Toksöz theory with the critical porosity model. Different from the traditional critical porosity model, critical porosity is not an empirical value but varied with pore shape and the ratio of bulk modulus versus shear modulus of the grain matrix. The substitution of the theoretical relationship into Kuster–Toksöz theory will generate the formulae for the bulk and shear moduli of multiple-porosity dry rocks, which is named the multiple-porosity variable critical porosity model. The new model has been used to predict elastic moduli for sandstone and carbonate rock. We compare the modelling results for P- and S-wave velocities and elastic moduli with the experimental data. The comparison shows that the new model can be used to describe the elastic properties for the rocks with multiple pore types. 相似文献
2.
We study the propagation of elastic waves that are generated in a fluid‐filled borehole surrounded by a cracked transversely isotropic medium. In the model studied the anisotropy and borehole axes coincide. To obtain the effective elastic moduli of a cracked medium we have applied Hudson's theory that enables the determination of the overall properties as a function of the crack orientation in relation to the symmetry axis of the anisotropic medium. This theory takes into account the hydrodynamic mechanism of the elastic‐wave attenuation caused by fluid filtration from the cracks into a porous matrix. We have simulated the full waveforms generated by an impulse source of finite length placed on the borehole axis. The kinematic and dynamic parameters of the compressional, shear and Stoneley waves as functions of the matrix permeability, crack orientation and porosity were studied. The modelling results demonstrated the influence of the crack‐system parameters (orientation and porosity) on the velocities and amplitudes of all wave types. The horizontally orientated cracks result in maximal decrease of the elastic‐wave parameters (velocities and amplitudes). Based on the fact that the shear‐ and Stoneley‐wave velocities in a transversely isotropic medium are determined by different shear moduli, we demonstrate the feasibility of the acoustic log to identify formations with close to horizontal crack orientations. 相似文献
3.
4.
The modelling of elastic waves in fractured media with an explicit finite‐difference scheme causes instability problems on a staggered grid when the medium possesses high‐contrast discontinuities (strong heterogeneities). For the present study we apply the rotated staggered grid. Using this modified grid it is possible to simulate the propagation of elastic waves in a 2D or 3D medium containing cracks, pores or free surfaces without hard‐coded boundary conditions. Therefore it allows an efficient and precise numerical study of effective velocities in fractured structures. We model the propagation of plane waves through a set of different, randomly cracked media. In these numerical experiments we vary the wavelength of the plane waves, the crack porosity and the crack density. The synthetic results are compared with several static theories that predict the effective P‐ and S‐wave velocities in fractured materials in the long wavelength limit. For randomly distributed and randomly orientated, rectilinear, non‐intersecting, thin, dry cracks, the numerical simulations of velocities of P‐, SV‐ and SH‐waves are in excellent agreement with the results of the modified (or differential) self‐consistent theory. On the other hand for intersecting cracks, the critical crack‐density (porosity) concept must be taken into account. To describe the wave velocities in media with intersecting cracks, we propose introducing the critical crack‐density concept into the modified self‐consistent theory. Numerical simulations show that this new formulation predicts effective elastic properties accurately for such a case. 相似文献
5.
Mechanical compaction and ultrasonic velocity of sands with different texture and mineralogical composition 总被引：1，自引：0，他引：1
Manzar Fawad Nazmul Haque Mondol Jens Jahren Knut Bjørlykke 《Geophysical Prospecting》2011,59(4):697-720
This study presents the results of experimental compaction while measuring ultrasonic velocities of sands with different grain size, shape, sorting and mineralogy. Uniaxial mechanical compaction tests up to a maximum of 50 MPa effective stress were performed on 29 dry sand aggregates derived from eight different sands to measure the rock properties. A good agreement was found between the Gassmann saturated bulk moduli of dry and brine saturated tests of selected sands. Sand samples with poor sorting showed low initial porosity while sands with high grain angularity had high initial porosity. The sand compaction tests showed that at a given stress well‐sorted, coarse‐grained sands were more compressible and had higher velocities (Vp and Vs) than fine‐grained sands when the mineralogy was similar. This can be attributed to grain crushing, where coarser grains lead to high compressibility and large grain‐to‐grain contact areas result in high velocities. At medium to high stresses the angular coarse to medium grained sands (both sorted sands and un‐sorted whole sands) showed high compaction and velocities (Vp and Vs). The small grain‐to‐grain contact areas promote higher deformation at grain contacts, more crushing and increased porosity loss resulting in high velocities. Compaction and velocities (Vp and Vs) increased with decreasing sorting in sands. However, at the same porosity, the velocities in whole sands were slightly lower than in the well‐sorted sands indicating the presence of loose smaller grains in‐between the framework grains. Quartz‐poor sands (containing less than 55% quartz) showed higher velocities (Vp and Vs) compared to that of quartz‐rich sands. This could be the result of sintering and enlargement of grain contacts of ductile mineral grains in the quartz‐poor sands increasing the effective bulk and shear stiffness. Tests both from wet measurements and Gassmann brine substitution showed a decreasing Vp/Vs ratio with increasing effective stress. The quartz‐rich sands separated out towards the higher side of the Vp/Vs range. The Gassmann brine substituted Vp and Vs plotted against effective stress provide a measure of the expected velocity range to be found in these and similar sands during mechanical compaction. Deviations of actual well log data from experimental data may indicate uplift, the presence of hydrocarbon, overpressure and/or cementation. Data from this study may help to model velocity‐depth trends and to improve the characterization of reservoir sands from well log data in a low temperature (<80–100o C) zone where compaction of sands is mostly mechanical. 相似文献
6.
Angus I. Best Lucy M. MacGregor Jeremy Sothcott Tim A. Minshull 《Geophysical Prospecting》2011,59(4):777-786
Improvements in the joint inversion of seismic and marine controlled source electromagnetic data sets will require better constrained models of the joint elastic‐electrical properties of reservoir rocks. Various effective medium models were compared to a novel laboratory data set of elastic velocity and electrical resistivity (obtained on 67 reservoir sandstone samples saturated with 35 g/l brine at a differential pressure of 8 MPa) with mixed results. Hence, we developed a new three‐phase effective medium model for sandstones with pore‐filling clay minerals based on the combined self‐consistent approximation and differential effective medium model. We found that using a critical porosity of 0.5 and an aspect ratio of 1 for all three components, the proposed model gave accurate model predictions of the observed magnitudes of P‐wave velocity and electrical resistivity and of the divergent trends of clean and clay‐rich sandstones at higher porosities. Using only a few well‐constrained input parameters, the new model offers a practical way to predict in situ porosity and clay content in brine saturated sandstones from co‐located P‐wave velocity and electrical resistivity data sets. 相似文献
7.
Estimates of depth, overpressure and amount of exhumation based on sonic data for a sedimentary formation rely on identification of a normal velocity–depth trend for the formation. Such trends describe how sonic velocity increases with depth in relatively homogeneous, brine‐saturated sedimentary formations as porosity is reduced during normal compaction (mechanical and chemical). Compaction is ‘normal’ when the fluid pressure is hydrostatic and the thickness of the overburden has not been reduced by exhumation. We suggest that normal porosity at the surface for a given lithology should be constrained by its critical porosity, i.e. the porosity limit above which a particular sediment exists only as a suspension. Consequently, normal velocity at the surface of unconsolidated sediments saturated with brine approaches the velocity of the sediment in suspension. Furthermore, porosity must approach zero at infinite depth, so the velocity approaches the matrix velocity of the rock and the velocity–depth gradient approaches zero. For sediments with initially good grain contact (when porosity is just below the critical porosity), the velocity gradient decreases with depth. By contrast, initially compliant sediments may have a maximum velocity gradient at some depth if we assume that porosity decreases exponentially with depth. We have used published velocity–porosity–depth relationships to formulate normal velocity–depth trends for consolidated sandstone with varying clay content and for marine shale dominated by smectite/illite. The first relationship is based on a modified Voigt trend (porosity scaled by critical porosity) and the second is based on a modified time‐average equation. Baselines for sandstone and shale in the North Sea agree with the established constraints and the shale trend can be applied to predict overpressure. A normal velocity–depth trend for a formation cannot be expressed from an arbitrary choice of mathematical functions and regression parameters, but should be considered as a physical model linked to the velocity–porosity transforms developed in rock physics. 相似文献
8.
成岩作用是影响砂岩声波速度的地质因素之一，定量表征压实和胶结作用的砂岩声波速度岩石物理模型具有重要的理论和实践应用意义.选取视压实率和视胶结率定量表征砂岩成岩作用，通过建立视压实率与颗粒配位数的关系将压实作用的影响引入修正的定量表征胶结作用的CCT模型，最终建立了一种能够定量表征压实和胶结作用对砂岩声波速度影响的岩石物理模型.理论考察发现，随胶结率的增大，岩石声波速度首先迅速增大，随后趋于稳定；随视压实率增大，岩石声波速度同样逐渐增大，当胶结率较大时声波速度变化更为明显.为了验证该声波速度模型，分别对人造砂岩和天然样品进行了声波速度实验观测，结果表明：实验结果与理论分析的趋势吻合良好.该模型易于使用，能够为应用地震和测井资料识别有利储层、定量评价孔隙度以及开展横波速度预测等应用提供理论基础. 相似文献
9.
Analytical approximations of bulk and shear moduli for dry rock based on the differential effective medium theory 总被引：2，自引：0，他引：2
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. 相似文献
10.
This paper is concerned with numerical tests of several rock physical relationships. The focus is on effective velocities and scattering attenuation in 3D fractured media. We apply the so‐called rotated staggered finite‐difference grid (RSG) technique for numerical experiments. Using this modified grid, it is possible to simulate the propagation of elastic waves in a 3D medium containing cracks, pores or free surfaces without applying explicit boundary conditions and without averaging the elastic moduli. We simulate the propagation of plane waves through a set of randomly cracked 3D media. In these numerical experiments we vary the number and the distribution of cracks. The synthetic results are compared with several (most popular) theories predicting the effective elastic properties of fractured materials. We find that, for randomly distributed and randomly orientated non‐intersecting thin penny‐shaped dry cracks, the numerical simulations of P‐ and S‐wave velocities are in good agreement with the predictions of the self‐consistent approximation. We observe similar results for fluid‐filled cracks. The standard Gassmann equation cannot be applied to our 3D fractured media, although we have very low porosity in our models. This is explained by the absence of a connected porosity. There is only a slight difference in effective velocities between the cases of intersecting and non‐intersecting cracks. This can be clearly demonstrated up to a crack density that is close to the connectivity percolation threshold. For crack densities beyond this threshold, we observe that the differential effective‐medium (DEM) theory gives the best fit with numerical results for intersecting cracks. Additionally, it is shown that the scattering attenuation coefficient (of the mean field) predicted by the classical Hudson approach is in excellent agreement with our numerical results. 相似文献
11.
12.
The frequency dependent mechanism of local fluid flow was found to be the decisive absorption and dispersion mechanism in fluid containing sandstones. In the ultrasonic frequency range local fluid flow and grain surface effects control the behaviour of highly porous and highly permeable rock if a pore fluid is present. Both mechanisms depend less on macroscopic rock parameters like porosity and permeability than essentially on microscopic parameters like crack size, crack density and grain contact properties. To demonstrate directly the important influence of the microstructure on the rock elastic and anelastic properties the microstructure of a sandstone was artificially changed by thermal cracking. The cracked rock exhibits a clearly changed behaviour at low uniaxial as well as at high hydrostatic pressure despite small changes of porosity and permeability. Fluid effects increase due to cracking. The experimental results are explained by means of a rock, model and local fluid flow. These results emphasize that it is the microstructure which controls the elastic and anelastic rock behaviour, even at high hydrostatic pressure. 相似文献
13.
Ismael Himar Falcon-Suarez Kelvin Amalokwu Jordi Delgado-Martin Ben Callow Katleen Robert Laurence North Sourav K. Sahoo Angus I. Best 《Geophysical Prospecting》2019,67(4):784-803
Synthetic rock samples can offer advantages over natural rock samples when used for laboratory rock physical properties studies, provided their success as natural analogues is well understood. The ability of synthetic rocks to mimic the natural stress dependency of elastic wave, electrical and fluid transport properties is of primary interest. Hence, we compare a consistent set of laboratory multi-physics measurements obtained on four quartz sandstone samples (porosity range 20–25%) comprising two synthetic and two natural (Berea and Corvio) samples, the latter used extensively as standards in rock physics research. We measured simultaneously ultrasonic (P- and S-wave) velocity and attenuation, electrical resistivity, permeability and axial and radial strains over a wide range of differential pressure (confining stress 15–50 MPa; pore pressure 5–10 MPa) on the four brine saturated samples. Despite some obvious physical discrepancies caused by the synthetic manufacturing process, such as silica cementation and anisotropy, the results show only small differences in stress dependency between the synthetic and natural sandstones for all measured parameters. Stress dependency analysis of the dry samples using an isotropic effective medium model of spheroidal pores and penny-shaped cracks, together with a granular cohesion model, provide evidence of crack closure mechanisms in the natural sandstones, seen to a much lesser extent in the synthetic sandstones. The smaller grain size, greater cement content, and cementation under oedometric conditions particularly affect the fluid transport properties of the synthetic sandstones, resulting in lower permeability and higher electrical resistivity for a similar porosity. The effective stress coefficients, determined for each parameter, are in agreement with data reported in the literature. Our results for the particular synthetic materials that were tested suggest that synthetic sandstones can serve as good proxies for natural sandstones for studies of elastic and mechanical properties, but should be used with care for transport properties studies. 相似文献
14.
15.
Osman Uyanık 《Journal of Applied Geophysics》2011,73(1):16-24
This study devises a new analytical relationship to determine the porosity of water-saturated soils at shallow depth using seismic compressional and shear wave velocities. Seismic refraction surveys together with soil sample collection were performed in selected areas containing water-saturated clay–silt, sand and gravely soils. Classification of clay–silt, sand and gravel dense soils provided the coefficient of experimental equation between the data sets, namely, Poisson's ratio, shear modulus and porosity values. This study presents a new analytical relationship between Poisson's ratio and shear modulus values, which are obtained from seismic velocities and porosity values of water-saturated material computed from water content and grain densities, which are determined by laboratory analysis of disturbed samples. The analytical relationship between data sets indicates that when the shear modulus of water-saturated loose soil increases, porosity decreases logarithmically. If shear modulus increases in dense or solid saturated soils, porosity decreases linearly. 相似文献
16.
孔隙度在岩石物理和地震勘探领域都是一个非常重要的参数,水饱和地层的孔隙度可以用于求海底地层的弹性模量进而求取纵横波速度、识别水合物以及估计水合物的储量等.然而在没有打钻的情况下要相对准确地了解海底未固结成岩地层的孔隙度变化趋势是困难的.历史上有学者提出了几个模型,都是基于孔隙度随深度的增加呈现指数下降的趋势,并且需要知道指数模型中的若干参数,因此就产生了两个问题,一是孔隙度是否随深度增加而呈现指数下降,二是指数模型中参数如何选取.本文从另外一个角度,即利用海底介质的压实性从表层开始往下逐层求取水饱和地层的孔隙度,此方法被称之为"体积模量法",它仅需要知道表层孔隙度和地层矿物组分,从而节省了勘探成本.本文将该方法应用于我国南海地区和美国ODP井数据,结果与测井数据以及地震反演数据吻合得较好. 相似文献
17.
18.
Yu-Qian Guo Hong-Da Ma Kai-Bo Shi Hong Cao Lu-Zhong Huang Feng-Chang Yao Tian-Yue Hu 《应用地球物理》2013,10(4):411-422
Most of the carbonates in the Tarim Basin in northwest China are low-porosity and low-permeability rocks. Owing to the complexity of porosity in carbonates, conventional rockphysics models do not describe the relation between velocity and porosity for the Tarim Basin carbonates well. We propose the porous-grain-upper-boundary (PGU) model for estimating the relation between velocity and porosity for low-porosity carbonates. In this model, the carbonate sediments are treated as packed media of porous elastic grains, and the carbonate pores are divided into isolated and connected pores The PGU model is modified from the porous-grain-stiff-sand (PGST) model by replacing the critical porosity with the more practical isolated porosity. In the implementation, the effective elastic constants of the porous grains are calculated by using the differential effective medium (DEM) model. Then, the elastic constants of connected porous grains in dry rocks are calculated by using the modified upper Hashin-Shtrikman bound. The application to the Tarim carbonates shows that relative to other conventional effective medium models the PGU model matches the well log data well. 相似文献
19.
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. 相似文献
20.
《Advances in water resources》2002,25(8-12):1105-1117
Macroscopic differential equations of mass and momentum balance for two immiscible fluids in a deformable porous medium are derived in an Eulerian framework using the continuum theory of mixtures. After inclusion of constitutive relationships, the resulting momentum balance equations feature terms characterizing the coupling among the fluid phases and the solid matrix caused by their relative accelerations. These terms, which imply a number of interesting phenomena, do not appear in current hydrologic models of subsurface multiphase flow. Our equations of momentum balance are shown to reduce to the Berryman–Thigpen–Chen model of bulk elastic wave propagation through unsaturated porous media after simplification (e.g., isothermal conditions, neglect of gravity, etc.) and under the assumption of constant volume fractions and material densities. When specialized to the case of a porous medium containing a single fluid and an elastic solid, our momentum balance equations reduce to the well-known Biot model of poroelasticity. We also show that mass balance alone is sufficient to derive the Biot model stress–strain relations, provided that a closure condition for porosity change suggested by de la Cruz and Spanos is invoked. Finally, a relation between elastic parameters and inertial coupling coefficients is derived that permits the partial differential equations of the Biot model to be decoupled into a telegraph equation and a wave equation whose respective dependent variables are two different linear combinations of the dilatations of the solid and the fluid. 相似文献