共查询到20条相似文献，搜索用时 512 毫秒
1.
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. 相似文献
2.
经典的微分等效介质(DEM)理论可用于确定多孔介质的弹性性质,但由于缺乏多重孔DEM方程,其估计的多重孔岩石的等效弹性模量依赖于包裹体(即不同孔隙纵横比的孔或缝)的添加顺序.本文首先从Kuster-Toksöz理论出发建立了Zimmermann和Norris两种形式的多重孔DEM方程.Norris形式的多重孔DEM方程预测的等效弹性模量总是位于Hashin-Shtrikman上下限内,而Zimmermann形式的多重孔DEM方程有时会越界.然后,通过使用干燥岩石模量比的解析近似式,对两个相互耦合的Norris形式DEM方程进行解耦得到干燥多重孔岩石的体积和剪切模量解析式.用全DEM方程的数值解对解析近似式的有效性进行了测试,解析公式的计算结果在整个孔隙度分布区间与数值解吻合良好.对实验室测量数据在假设岩石含有双重孔隙的情形下用双重孔DEM解析公式对岩石的弹性模量进行了预测,结果表明,解析式准确地预测了弹性模量随孔隙度的变化.双重孔(即软、硬孔)DEM解析模型可用来反演各孔隙类型的孔隙体积比,它可以通过实验室测量与理论预测之间的平方误差最小反演得到.砂岩样品的反演结果揭示,软孔的孔隙体积百分比与粘土含量没有明显的相关性. 相似文献
3.
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. 相似文献
4.
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. 相似文献
5.
Although there is no assumption of pore geometry in derivation of Gassmann's equation, the pore geometry is in close relation with hygroscopic water content and pore fluid communication between the micropores and the macropores. The hygroscopic water content in common reservoir rocks is small, and its effect on elastic properties is ignored in the Gassmann theory. However, the volume of hygroscopic water can be significant in shaly rocks or rocks made of fine particles; therefore, its effect on the elastic properties may be important. If the pore fluids in microspores cannot reach pressure equilibrium with the macropore system, assumption of the Gassmann theory is violated. Therefore, due to pore structure complexity, there may be a significant part of the pore fluids that do not satisfy the assumption of the Gassmann theory. We recommend that this part of pore fluids be accounted for within the solid rock frame and effective porosity be used in Gassmann's equation for fluid substitution. Integrated study of ultrasonic laboratory measurement data, petrographic data, mercury injection capillary pressure data, and nuclear magnetic resonance T2 data confirms rationality of using effective porosity for Gassmann fluid substitution. The effective porosity for Gassmann's equation should be frequency dependent. Knowing the pore geometry, if an empirical correlation between frequency and the threshold pore‐throat radius or nuclear magnetic resonance T2 could be set up, Gassmann's equation can be applicable to data measured at different frequencies. Without information of the pore geometry, the irreducible water saturation can be used to estimate the effective porosity. 相似文献
6.
Luděk Klimeš 《Studia Geophysica et Geodaetica》2018,62(2):255-260
For a given stiffness tensor (tensor of elastic moduli) of a generally anisotropic medium, we can estimate the extent to which the medium is transversely isotropic, and determine the direction of its reference symmetry axis. In this paper, we rotate the given stiffness tensor about this reference symmetry axis, and determine the reference transversely isotropic (uniaxial) stiffness tensor as the average of the rotated stiffness tensor over all angles of rotation. The obtained reference transversely isotropic (uniaxial) stiffness tensor represents an analytically differentiable approximation of the given generally anisotropic stiffness tensor. The proposed analytic method is compared with a previous numerical method in two numerical examples. 相似文献
7.
Luděk Klimeš 《Studia Geophysica et Geodaetica》2016,60(3):391-402
For a given stiffness tensor (tensor of elastic moduli) of a generally anisotropic medium, we estimate to what extent the medium is transversely isotropic (uniaxial) and determine the direction of its reference symmetry axis expressed in terms of the unit reference symmetry vector. If the medium is exactly transversely isotropic (exactly uniaxial), we obtain the direction of its symmetry axis. We can also calculate the first–order and second–order spatial derivatives of the reference symmetry vector which may be useful in tracing the reference rays for the coupling ray theory. The proposed method is tested using various transversely isotropic (uniaxial) and approximately transversely isotropic (approximately uniaxial) media. 相似文献
8.
《Journal of Applied Geophysics》2006,58(3):253-262
Attenuation data extracted from full waveform sonic logs is sensitive to vuggy and matrix porosities in a carbonate aquifer. This is consistent with the synthetic attenuation (1 / Q) as a function of depth at the borehole-sonic source-peak frequency of 10 kHz. We use velocity and densities versus porosity relationships based on core and well log data to determine the matrix, secondary, and effective bulk moduli. The attenuation model requires the bulk modulus of the primary and secondary porosities. We use a double porosity model that allows us to investigate attenuation at the mesoscopic scale. Thus, the secondary and primary porosities in the aquifer should respond with different changes in fluid pressure. The results show a high permeability region with a Q that varies from 25 to 50 and correlates with the stiffer part of the carbonate formation. This pore structure permits water to flow between the interconnected vugs and the matrix. In this region the double porosity model predicts a decrease in the attenuation at lower frequencies that is associated with fluid flowing from the more compliant high-pressure regions (interconnected vug space) to the relatively stiff, low-pressure regions (matrix). The chalky limestone with a low Q of 17 is formed by a muddy porous matrix with soft pores. This low permeability region correlates with the low matrix bulk modulus. A low Q of 18 characterizes the soft sandy carbonate rock above the vuggy carbonate.This paper demonstrates the use of attenuation logs for discriminating between lithology and provides information on the pore structure when integrated with cores and other well logs. In addition, the paper demonstrates the practical application of a new double porosity model to interpret the attenuation at sonic frequencies by achieving a good match between measured and modeled attenuation. 相似文献
9.
Although clay is composed of disconnected anisotropic clay platelets, many rock physics models treat the clay platelets in shale as interconnected. However, the clay matrix in shales can be modelled as anisotropic clay platelets embedded within a soft isotropic interplatelet region, allowing the influence of disconnected clay platelets on the elastic properties of the clay matrix to be analysed. In this model, properties of the interplatelet region are governed by its effective bulk and shear moduli, whereas the effective properties of the clay platelets are governed by their volume fraction, aspect ratio and elastic stiffness tensor. Together, these parameters implicitly account for variations in clay and fluid properties, as well as fluid saturation. Elastic stiffnesses of clay platelets are obtained from the literature, including both experimental measurements and first-principles calculations of the full anisotropic (monoclinic or triclinic) elastic stiffness tensors of layered silicates. These published elastic stiffness tensors are used to compile a database of equivalent transverse isotropic elastic stiffness tensors, and other physical properties, for eight common varieties of layered silicates. Clay matrix anisotropy is then investigated by examining the influence of these different elastic stiffnesses, and of varying model parameters, upon the effective transverse isotropic elastic stiffness tensor of the clay matrix. The relationship between the different clay minerals and their associated anisotropy parameters is studied, and their impact on the resulting anisotropy of the clay matrix is analysed. 相似文献
10.
Alexander Y. Rozhko 《Geophysical Prospecting》2020,68(2):631-656
It is evident from the laboratory experiments that shear moduli of different porous isotropic rocks may show softening behaviour upon saturation. The shear softening means that the shear modulus of dry samples is higher than of saturated samples. Shear softening was observed both at low (seismic) frequencies and high (ultrasonic) frequencies. Shear softening is stronger at seismic frequencies than at ultrasonic frequencies, where the softening is compensated by hardening due to unrelaxed squirt flow. It contradicts to Gassmann's theory suggesting that the relaxed shear modulus of isotropic rock should not depend upon fluid saturation, provided that no chemical reaction between the solid frame and the pore fluid. Several researchers demonstrated that the shear softening effect is reversible during re-saturation of rock samples, suggesting no permanent chemical reaction between the solid frame and the pore fluid. Therefore, it is extremely difficult to explain this fluid–rock interaction mechanism theoretically, because it does not contradict to the assumptions of Gassmann's theory, but contradicts to its conclusions. We argue that the observed shear softening of partially saturated rocks by different pore fluids is related to pore-scale interfacial phenomena effects, typically neglected by the rock physics models. These interface phenomena effects are dependent on surface tension between immiscible fluids, rock wettability, aperture distribution of microcracks, compressibility of microcracks, porosity of microcracks, elastic properties of rock mineral, fluid saturation, effective stress and wave amplitude. Derived equations allow to estimate effects of pore fluids and saturation on the shear modulus and mechanical strength of rocks. 相似文献
11.
Ida L. Fabricius Christian Høier Peter Japsen Uffe Korsbech 《Geophysical Prospecting》2007,55(4):487-506
In impure chalk, the elastic moduli are not only controlled by porosity but also by contact‐cementation, resulting in relatively large moduli for a given porosity, and by admixtures of clay and fine silica, which results in relatively small moduli for a given porosity. Based on a concept of solids suspended in pore fluids as well as composing the rock frame, we model P‐wave and S‐wave moduli of dry and wet plug samples by an effective‐medium Hashin–Shtrikman model, using chemical, mineralogical and textural input. For a given porosity, the elastic moduli correspond to a part of the solid (the iso‐frame value) forming the frame of an Upper Hashin–Shtrikman bound, whereas the remaining solid is modelled as suspended in the pore fluid. The iso‐frame model is thus a measure of the pore‐stiffness or degree of cementation of the chalk. The textural and mineralogical data may be assessed from logging data on spectral gamma radiation, density, sonic velocity and water saturation in a hydrocarbon zone, whereas the iso‐frame value of a chalk may be assessed from the density and acoustic P‐wave logs alone. The iso‐frame concept may thus be directly used in conventional log‐analysis and is a way of incorporating sonic‐logging data. The Rigs‐1 and Rigs‐2 wells in the South Arne field penetrate the chalk at the same depth but differ in porosity and in water saturation although almost the entire chalk interval has irreducible water saturation. Our model, combined with petrographic data, indicates that the difference in porosity is caused by a higher degree of pore‐filling cementation in Rigs‐1. Petrographic data indicate that the difference in water saturation is caused by a higher content of smectite in the pores of Rigs‐1. In both wells, we find submicron‐size diagenetic quartz. 相似文献
12.
在龙马溪页岩微观物性特征分析的基础上,综合利用测井解释、微观测试分析资料,建立了一种适用于龙马溪页岩的横观各向同性岩石物理模型,该模型建模过程：将各向异性SCA和DEM模型联合模拟得到的黏土和干酪根混合物作为背景介质;采用SCA模型对脆性矿物混合物进行模拟,利用各向异性DEM将脆性矿物混合物添加到背景介质;进一步将空孔隙添加到页岩基质,并利用Brown-Korringa模型进行各向异性条件下的流体替换,从而得到横观各向同性页岩岩石物理模型.通过对四川盆地A井龙马溪页岩进行岩石物理建模分析,计算了孔隙纵横比、纵横波速、各向异性系数和弹性参数,检验了模型的准确性.研究结果表明：矿物颗粒和孔隙纵横比是影响模型精度的关键参数,黏土和干酪根颗粒纵横比为0.05,图像识别获得的脆性矿物颗粒纵横比主要分布于0.45~1.0（集中分布于0.5~0.85）,横波波速反演获得的孔隙纵横比主要分布于0.1~0.3（平均值约为0.22）;模型预测和实测纵波波速之间误差为-2.40%~2.21%（平均绝对误差仅1.20%）,预测和实测横波波速之间误差为-1.93%~1.42%（平均绝对误差仅0.64%）,证实了本文模型的准确性和精度.本文模型能够准确计算页岩5个独立的刚度系数,为页岩弹性参数、声波波速、各向异性和脆性分析提供了有效手段,也为后续地球物理和工程地质参数分析提供了重要依据. 相似文献
13.
在Kuster-Toksöz（KT）模型与差分等效介质理论（DEM）结合过程中,通过用孔隙包含物逐渐替换基质的方法获得孔隙介质的岩石弹性模量。现有方法每次替换的孔隙体积是常量,而基质体积不断减小,实际参与替换的包含物体积与孔隙包含物计算体积是不同的。本文通过改进每次替换的孔隙体积计算公式,使得替换体积随迭代次数的增加而逐渐减小,保持替换体积相对基质体积的比率不变,在迭代次数足够大的条件下使得该比率足够小,满足了K-T计算公式的要求,计算结果更接近理论值。测试结果显示:随着孔隙度的增大,岩石等效弹性模量逐渐收敛于孔隙包含物的弹性模量,说明本方法与物理规律一致;与现有的KT迭代方法相比,采用本方法的计算结果与Xu-White模型更接近,本方法提出孔隙包含物的实际体积与计算体积的计算式更符合KT模型孔隙切分过程。 相似文献
14.
Maxim Lebedev Moyra E.J. Wilson Vassily Mikhaltsevitch 《Geophysical Prospecting》2014,62(6):1253-1265
Petrophysical properties of carbonate reservoirs are less predictable than that of siliciclastic reservoirs. One of the main reasons for this is the physical and chemical interactions of carbonate rocks with pore fluids. Such interactions can significantly change the elastic properties of the rock matrix and grains, making the applicability of Gassmann's fluid substitution procedure debatable. This study is an attempt to understand the mechanisms of fluid‐rock interactions and the influence of these interactions on elastic parameters of carbonates. We performed precise indentation tests on Savonnières limestone at a microscale level under dry, distilled water, and n‐Decane saturated conditions. Our experiments display softening of the rock matrix after water saturation. We have found that mainly the ooid cortices, peloid nuclei and prismatic intergranular cement are affected by water flooding. We also observed a shear modulus reduction in Savonnières limestone in an experiment performed at ultrasonic frequencies. One of the most important results obtained in our experimental study is that the Gassmann fluid substitution theory might not always be applicable to predict the elastic moduli of fluid‐saturated limestones. 相似文献
15.
Elastic properties of fluid saturated porous media with aligned fractures can be studied using the model of fractures as linear-slip interfaces in an isotropic porous background. Such a medium represents a particular case of a transversely isotropic (TI) porous medium, and as such can be analyzed with equations of anisotropic poroelasticity. This analysis allows the derivation of explicit analytical expressions for the low-frequency elastic constants and anisotropy parameters of the fractured porous medium saturated with a given fluid. The five elastic constants of the resultant TI medium are derived as a function of the properties of the dry (isotropic) background porous matrix, fracture properties (normal and shear excess compliances), and fluid bulk modulus. For the particular case of penny-shaped cracks, the expression for anisotropy parameter ε has the form similar to that of Thomsen [Geophys. Prospect. 43 (1995) 805]. However, contrary to the existing view, the compliance matrix of a fluid-saturated porous-fractured medium is not equivalent to the compliance matrix of any equivalent solid medium with a single set of parallel fractures. This unexpected result is caused by the wave-induced flow of fluids between pores and fractures. 相似文献
16.
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. 相似文献
17.
JAN DOUMA 《Geophysical Prospecting》1989,37(7):831-849
Media containing aligned cracks or ellipsoidal inclusions as well as media consisting of sequences of isotropic layers show transverse isotropy with respect to elastic wave propagation. However, the transversely isotropic media which are equivalent to media containing aligned inclusions do not necessarily have to be representable by sequences of stable isotropic layers. These transversely isotropic media can be modelled by such sequences if - and only if - several stability conditions are satisfied. Important parameters determining whether these conditions are satisfied are the aspect ratio of the inclusions and the material filling the inclusions, the‘fluid’. An analytical expression describing the range of aspect ratios for which the constraints are satisfied can be derived. This expression (which is a good approximation for several crack models) and numerical calculations show that media containing water-filled inclusions can be represented by sequences of stable isotropic layers if the inclusions have aspect ratios less than 0.1. The limiting aspect ratio decreases for a decreasing ratio of the bulk modulus of the fluid to the shear modulus of the matrix material. Finally, media containing dry inclusions of any aspect ratio cannot be modelled by thin isotropic layering. These results depend only weakly on the crack density and on the matrix material. The representation of crack-induced anisotropy by layer-induced anisotropy can be used to classify crack-induced anisotropy and might be useful in the separation of the cause of anisotropy and the determination of the nature of the fluid. 相似文献
18.
The shales of the Qiongzhusi Formation and Wufeng–Longmaxi Formations at Sichuan Basin and surrounding areas are presently the most important stratigraphic horizons for shale gas exploration and development in China. However, the regional characteristics of the seismic elastic properties need to be better determined. The ultrasonic velocities of shale samples were measured under dry conditions and the relations between elastic properties and petrology were systemically analyzed. The results suggest that 1) the effective porosity is positively correlated with clay content but negatively correlated with brittle minerals, 2) the dry shale matrix consists of clays, quartz, feldspars, and carbonates, and 3) organic matter and pyrite are in the pore spaces, weakly coupled with the shale matrix. Thus, by assuming that all connected pores are only present in the clay minerals and using the Gassmann substitution method to calculate the elastic effect of organic matter and pyrite in the pores, a relatively simple rock-physics model was constructed by combining the self-consistent approximation (SCA), the differential effective medium (DEM), and Gassmann’s equation. In addition, the effective pore aspect ratio was adopted from the sample averages or estimated from the carbonate content. The proposed model was used to predict the P-wave velocities and generally matched the ultrasonic measurements very well. 相似文献
19.
Yongyang Sun Boris Gurevich Maxim Lebedev Stanislav Glubokovskikh Vassili Mikhaltsevitch Junxin Guo 《Geophysical Prospecting》2019,67(4):888-899
Quantifying the effects of pore-filling materials on elastic properties of porous rocks is of considerable interest in geophysical practice. For rocks saturated with fluids, the Gassmann equation is proved effective in estimating the exact change in seismic velocity or rock moduli upon the changes in properties of pore infill. For solid substance or viscoelastic materials, however, the Gassmann theory is not applicable as the rigidity of the pore fill (either elastic or viscoelastic) prevents pressure communication in the pore space, which is a key assumption of the Gassmann equation. In this paper, we explored the elastic properties of a sandstone sample saturated with fluid and solid substance under different confining pressures. This sandstone sample is saturated with octadecane, which is a hydrocarbon with a melting point of 28°C, making it convenient to use in the lab in both solid and fluid forms. Ultrasonically measured velocities of the dry rock exhibit strong pressure dependency, which is largely reduced for the filling of solid octadecane. Predictions by the Gassmann theory for the elastic moduli of the sandstone saturated with liquid octadecane are consistent with ultrasonic measurements, but underestimate the elastic moduli of the sandstone saturated with solid octadecane. Our analysis shows that the difference between the elastic moduli of the dry and solid-octadecane-saturated sandstone is controlled by the squirt flow between stiff, compliant, and the so-called intermediate pores (with an aspect ratio larger than that of compliant pore but much less than that of stiff pores). Therefore, we developed a triple porosity model to quantify the combined squirt flow effects of compliant and intermediate pores saturated with solid or viscoelastic infill. Full saturation of remaining stiff pores with solid or viscoelastic materials is then considered by the lower embedded bound theory. The proposed model gave a reasonable fit to the ultrasonic measurements of the elastic moduli of the sandstone saturated with liquid or solid octadecane. Comparison of the predictions by the new model to other solid substitution schemes implied that accounting for the combined effects of compliant and intermediate pores is necessary to explain the solid squirt effects. 相似文献