首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
We generalize the classical theory of acoustoelasticity to the porous case (one fluid and a solid frame) and finite deformations. A unified treatment of non‐linear acoustoelasticity of finite strains in fluid‐saturated porous rocks is developed on the basis of Biot’s theory. A strain‐energy function, formed with eleven terms, combined with Biot’s kinetic and dissipation energies, yields Lagrange’s equations and consequently the wave equation of the medium. The velocities and dissipation factors of the P‐ and S‐waves are obtained as a function of the 2nd‐ and 3rd‐order elastic constants for hydrostatic and uniaxial loading. The theory yields the limit to the classical theory if the fluid is replaced with a solid with the same properties of the frame. We consider sandstone and obtain results for open‐pore jacketed and closed‐pore jacketed ‘gedanken’ experiments. Finally, we compare the theoretical results with experimental data.  相似文献   

2.
Contrary to the traditional view, seismic attenuation in Biot's theory of fluid-saturated porous media is due to viscous damping of local (not global) pore-fluid motion. Since substantial inhomogeneities in fluid permeability of porous geological materials are to be expected, the regions of highest local permeability contribute most to the wave energy dissipation while those of lowest permeability dominate the fluid flow rate if they are uniformly distributed. This dichotomy can explain some of the observed discrepancies between computed and measured attenuation of compressional and shear waves in porous earth. One unfortunate consequence of this result is the fact that measured seismic wave attenuation in fluid-filled geological materials cannot be used directly as a diagnostic of the global fluid-flow permeability.  相似文献   

3.
Energy loss in porous media containing fluids is typically caused by a variety of dynamic mechanisms.In the Biot theory,energy loss only includes the frictional dissipation between the solid phase and the fluid phase,resulting in underestimation of the dispersion and attenuation of the waves in the low frequency range.To develop a dynamic model that can predict the high dispersion and strong attenuation of waves at the seismic band,we introduce viscoelasticity into the Biot model and use fractional derivatives to describe the viscoelastic mechanism,and finally propose a new wave propagation model.Unlike the Biot model,the proposed model includes the intrinsic dissipation of the solid frame.We investigate the effects of the fractional order parameters on the dispersion and attenuation of the P-and S-waves using several numerical experiments.Furthermore,we use several groups of experimental data from different fluid-saturated rocks to testify the validity of the new model.The results demonstrate that the new model provides more accurate predictions of high dispersion and strong attenuation of different waves in the low frequency range.  相似文献   

4.
Saturation of porous rocks with a mixture of two fluids (known as partial saturation) has a substantial effect on the seismic waves propagating through these rocks. In particular, partial saturation causes significant attenuation and dispersion of the propagating waves, due to wave-induced fluid flow. Such flow arises when a passing wave induces different fluid pressures in regions of rock saturated by different fluids. As partial fluid saturation can occur on different length scales, attenuation due to wave-induced fluid flow is ubiquitous. In particular, mesoscopic fluid flow due to heterogeneities occurring on a scale greater than porescale, but less than wavelength scale, is responsible for significant attenuation in the frequency range from 10 to 1000 Hz.Most models of attenuation and dispersion due to mesoscopic heterogeneities imply that fluid heterogeneities are distributed in a periodic/regular way. In 1D this corresponds to periodically alternating layering, in 3D as periodically distributed inclusions of a given shape (usually spheres). All these models yield very similar estimates of attenuation and dispersion.Experimental studies show that mesoscopic heterogeneities have less idealized distributions and that the distribution itself affects attenuation and dispersion. Therefore, theoretical models are required which would simulate the effect of more general and realistic fluid distributions.We have developed two theoretical models which simulate the effect of random distributions of mesoscopic fluid heterogeneities. The first model assumes that one fluid forms a random ensemble of spherical inclusions in a porous medium saturated by the other fluid. The attenuation and dispersion predicted by this model are very similar to those predicted for 3D periodic distribution. Attenuation (inverse quality factor) is proportional to ω at low frequencies for both distributions. This is in contrast to the 1D case, where random and periodically alternating layering shows different attenuation behaviour at low frequencies. The second model, which assumes a 3D continuous distribution of fluid heterogeneities, also predicts the same low-frequency asymptote of attenuation. However, the shapes of the frequency dependencies of attenuation are different. As the 3D continuous random approach assumes that there will be a distribution of different patch sizes, it is expected to be better suited to modelling experimental results. Further research is required in order to uncover how to relate the random functions to experimentally significant parameters.  相似文献   

5.
This article examines the conditions under which the pressure-work and viscous dissipation terms should be retained in the energy balance relation for single (liquid water or vapor) and two-phase (liquid water and vapor) fluid flow through porous media. It is shown that if one wishes to retain the pressure-work term, then one must also keep the viscous dissipation term in the energy balance. Consideration of steady non-isothermal radial flow demonstrates that both pressure-work and viscous dissipation are liable to have negligibly small effects in single phase liquid water and in two-phase liquid-vapor systems. This conclusion is, however, not generally valid for pure vapor systems; in this case, pressure-work and viscous dissipation can produce significant variations in the computed reservoir response.  相似文献   

6.
由于介观尺度的孔隙流体流动,弹性波传播过孔隙岩层时在地震频段表现出较强的频散和衰减。Johnson理论给出了在任意孔隙形状的条件下,部分气水饱和孔隙介质的理论相速度和品质因子的解析解。本文在Johnson模型的基础上,通过对Q值曲线的低频和高频近似,推导了Q值曲线的近似公式,以及基于孔隙介质基本地球物理参数和孔隙斑块几何形态参数T和比表面积S/V的最大衰减Qmin近似公式。通过与理论值的对比,对Qmin近似公式存在的线性误差进行改正,进一步提高了精度。复杂的斑块形态对最大衰减Qmin和过渡频率ftr的都产生一定影响,且对ftr影响更大。因为数值模拟直接求解介观尺度的Biot孔隙介质方程需要极大的计算量,我们使用Zener模型建立了等效粘弹模型,有效地模拟了地震频带内的衰减和频散现象。  相似文献   

7.
Wave‐induced fluid flow plays an important role in affecting the seismic dispersion and attenuation of fractured porous rocks. While numerous theoretical models have been proposed for the seismic dispersion and attenuation in fractured porous rocks, most of them neglect the wave‐induced fluid flow resulting from the background anisotropy (e.g. the interlayer fluid flow between different layers) that can be normal in real reservoirs. Here, according to the theories of poroelasticity, we present an approach to study the frequency‐dependent seismic properties of more realistic and complicated rocks, i.e. horizontally and periodically layered porous rock with horizontal and randomly orienting fractures, respectively, distributed in one of the two periodical layers. The approach accounts for the dual effects of the wave‐induced fluid flow between the fractures and the background pores and between different layers (the interlayer fluid flow). Because C33 (i.e., the modulus of the normally incident P‐wave) is directly related to the P‐wave velocity widely measured in the seismic exploration, and its comprehensive dispersion and attenuation are found to be most significant, we study mainly the effects of fracture properties and the stiffness contrast between the different layers on the seismic dispersion and attenuation of C33. The results show that the increasing stiffness contrast enhances the interlayer fluid flow of the layered porous rocks with both horizontal and randomly orienting fractures and weakens the wave‐induced fluid flow between the fractures and the background pores, especially for the layered porous rock with horizontal fractures. The modelling results also demonstrate that for the considered rock construction, the increasing fracture density reduces the interlayer fluid flow while improves the dispersion and attenuation in the fracture‐relevant frequency band. Increasing fracture aspect ratio is found to reduce the dispersion and attenuation in the fracture‐relevant frequency band only, especially for the layered porous rock with horizontal fractures.  相似文献   

8.
In sedimentary rocks attenuation/dispersion is dominated by fluid-rock interactions. Wave-induced fluid flow in the pores causes energy loss through several mechanisms, and as a result attenuation is strongly frequency dependent. However, the fluid motion process governing the frequency dependent attenuation and velocity remains unclear. We propose a new approach to obtain the analytical expressions of pore pressure, relative fluxes distribution and frame displacement within the double-layer porous media based on quasi-static poroelastic theory. The dispersion equation for a P-wave propagating in a porous medium permeated by aligned fractures is given by considering fractures as thin and highly compliant layers. The influence of mesoscopic fluid flow on phase velocity dispersion and attenuation is discussed under the condition of varying fracture weakness. In this model conversion of the compression wave energy into Biot slow wave diffusion at the facture surface can result in apparent attenuation and dispersion within the usual seismic frequency band. The magnitude of velocity dispersion and attenuation of P-wave increases with increasing fracture weakness, and the relaxation peak and maximum attenuation shift towards lower frequency. Because of its periodic structure, the fractured porous media can be considered as a phononic crystal with several pass and stop bands in the high frequency band. Therefore, the velocity and attenuation of the P-wave show an oscillatory behavior with increasing frequency when resonance occurs. The evolutions of the pore pressure and the relative fluxes as a function of frequency are presented, giving more physical insight into the behavior of P-wave velocity dispersion and the attenuation of fractured porous medium due to the wave-induced mesoscopic flow. We show that the specific behavior of attenuation as function of frequency is mainly controlled by the energy dissipated per wave cycle in the background layer.  相似文献   

9.
流体分布对松散介质中P波速度和衰减的影响   总被引:3,自引:0,他引:3       下载免费PDF全文
为了深入研究流体对岩石中弹性波速度和衰减的影响,必须考虑到流体的分布和粘性。引入气体包裹体模型来研究粘性流体的分布对松散介质中P波速度和衰减的影响,用气泡平均半径来描述流体分布的不均匀性,计算了不同气泡半径和频率下P波速度和衰减随饱和度变化的曲线,并与有效流体模型作了比较,由于流体喷流的存在会使Gassmann方程在高频下不适用,用干燥和饱和流体的P波、S波速度修正了理论曲线。测量了玻璃微珠中不同水饱和度下高频P波的速度和衰减,并尝试用峰值频率来计算衰减。此方法求出的Q和频谱比法求出的Q在干燥或饱和水时基本相同,随饱和度的变化规律也基本一致,但衰减峰的大小有差异。根据实测值得似事经修正的波速和衰减理论曲线从而估算出气泡平均半径,认为P波速度和衰减不仅与饱和度有关而且也与介质内部气体-液体压力平衡有关。  相似文献   

10.
We derived the velocity and attenuation of a generalized Stoneley wave being a symmetric trapped mode of a layer filled with a Newtonian fluid and embedded into either a poroelastic or a purely elastic rock. The dispersion relation corresponding to a linearized Navier–Stokes equation in a fracture coupling to either Biot or elasticity equations in the rock via proper boundary conditions was rigorously derived. A cubic equation for wavenumber was found that provides a rather precise analytical approximation of the full dispersion relation, in the frequency range of 10?3 Hz to 103 Hz and for layer width of less than 10 cm and fluid viscosity below 0.1 Pa· s [100 cP]. We compared our results to earlier results addressing viscous fluid in either porous rocks with a rigid matrix or in a purely elastic rock, and our formulae are found to better match the numerical solution, especially regarding attenuation. The computed attenuation was used to demonstrate detectability of fracture tip reflections at wellbore, for a range of fracture lengths and apertures, pulse frequencies, and fluid viscosity.  相似文献   

11.
岩石的粘弹性谐振Q模型   总被引:7,自引:2,他引:7       下载免费PDF全文
实验表明,在应变小于10-6范围内,砂岩对地震波的吸收主要由孔隙流体的局部运动引起,而且Q值随频率的变化出现谐振现象。据此,作者认为,地壳上部地震波能量的耗损可以表示为以下两种成分的线性组合:1.由滑动摩擦、热弹性驰豫、位错运动等引起的能量消散,它们主要与岩石的固相成分有关,可近似地用恒Q模型描述;2.由孔隙流体运动引起的能量消散,可以用谐振公式描述。根据这一认识,利用描述因果关系的Kramers-Krnig关系式可推导出表示地壳上部岩石粘弹性(复弹性模量,相速度频散和衰减函数)的公式,它们综合地描述了由各种机制引起的波的频散和吸收,并在谐振Q值等于参考常数Q值时退化为目前常用的Futterman模型。作为这种谐振Q模型的应用,介绍了它用于Q值测量结果外推和频散一吸收研究以及粘弹性介质中反射地震道合成的结果。  相似文献   

12.
Biot's theory is employed to study the reflection and transmission ofSH waves in a sandy layer lying over a fluid-saturated porous solid half-space. The entire medium is considered under constant initial stress. Effects of sandiness, initial stress, anelasticity and viscosity of the interstitial fluid on the partitioning of energy are studied. In the presence of initial stress the incident wave starts attenuating when incider beyond a certain angle (depending upon the amount of initial stress), even if the medium is perfectly clastic. Anelasticity of the solid layer results in the dissipation of energy during transmission. The direction of attenuation vector of incident wave affects the dissipation energy to a large extent. Effect on partitioning of energy reverse at incidence after the critical angle. A complete account of energy returmed back to the underlying half-space and that which is dissipated in the overlying layer has been discussed analytically as well as numerically.  相似文献   

13.
Partially saturated reservoirs are one of the major sources of seismic wave attenuation, modulus defect and velocity dispersion in real seismic data. The main attenuation and dispersion phenomenon is wave induced fluid flow due to the heterogeneity in pore fluids or porous rock. The identification of pore fluid type, saturation and distribution pattern within the pore space is of great significance as several seismic and petrophysical properties of porous rocks are largely affected by fluid type, saturation and fluid distribution pattern. Based on Gassmann-Wood and Gassmann- Hill rock physics models modulus defect, velocity dispersion and attenuation in Jurassic siliclastic partially-saturated rocks are studied. For this purpose two saturation patterns - uniform and patchy - are considered within the pore spaces in two frequency regimes i.e., lower frequency and higher frequency. The results reveal that at low enough frequency where saturation of liquid and gas is uniform, the seismic velocity and bulk modulus are lower than at higher frequency where saturation of fluid mixture is in the form of patches. The velocity dispersion and attenuation is also modeled at different levels of gas saturation. It is found that the maximum attenuation and velocity dispersion is at low gas saturation. Therefore, the dispersion and attenuation can provide a potential way to predict gas saturation and can be used as a property to differentiate low from high gas saturation.  相似文献   

14.
Seismic wave attenuation in porous rocks consists of intrinsic or anelastic attenuation (the lost energy is converted into heat due to interaction between the waves and the rocks) and the extrinsic or geometric attenuation (the energy is lost due to beam spreading, transmission loss and scattering). The first is of great importance because it can give additional information on the petrophysical properties of rocks (permeability, degree of saturation, type of saturant, etc.). The most difficult problem in attenuation measurements is estimating or eliminating extrinsic attenuation, so that the intrinsic attenuation can be obtained. To date, in laboratory attenuation measurements using wave propagation, several methods have been used. The difficulties vary with the method. The coupling effect and the geometric divergence or beam spreading are the major problems. Papadakis’ diffraction corrections have been used extensively by Winkler and Plona in their modified pulse-echo high-pressure attenuation measurements. These corrections are computed for homogeneous liquid media and their failure to fit data for solid material implies that these corrections must be used with caution, especially for high Q values. Three new methods for laboratory ultrasonic attenuation measurements are presented. The first is the ‘ultrasonic lens’ method for attenuation measurements at atmospheric pressure, in which an ultrasonic lens placed between transmitter and sample transforms the initially oblique incident beam into normal incidence so that the geometric divergence is eliminated. The second method is the ‘panoramic receiver’, in which the beam spreading can be eliminated by integrating the ultrasonic energy over a large area. The third method is called 'self-spectral ratio’ and is applicable for all pressure conditions. Attenuation is estimated by comparing two signals recorded on the same rock but with two slightly different thicknesses under the same pressure conditions. Hence the extrinsic attenuation for both thicknesses is approximately the same. A comparison between the self-spectral ratio method and that of Winkler and Plona demonstrates a very good agreement for a broad band of frequencies. Hence the Winkler-Plona technique and Papadakis’ diffraction corrections can be accepted as reliable in any future work.  相似文献   

15.
The presence of fractures in fluid‐saturated porous rocks is usually associated with strong seismic P‐wave attenuation and velocity dispersion. This energy dissipation can be caused by oscillatory wave‐induced fluid pressure diffusion between the fractures and the host rock, an intrinsic attenuation mechanism generally referred to as wave‐induced fluid flow. Geological observations suggest that fracture surfaces are highly irregular at the millimetre and sub‐millimetre scale, which finds its expression in geometrical and mechanical complexities of the contact area between the fracture faces. It is well known that contact areas strongly affect the overall mechanical fracture properties. However, existing models for seismic attenuation and velocity dispersion in fractured rocks neglect this complexity. In this work, we explore the effects of fracture contact areas on seismic P‐wave attenuation and velocity dispersion using oscillatory relaxation simulations based on quasi‐static poroelastic equations. We verify that the geometrical and mechanical details of fracture contact areas have a strong impact on seismic signatures. In addition, our numerical approach allows us to quantify the vertical solid displacement jump across fractures, the key quantity in the linear slip theory. We find that the displacement jump is strongly affected by the geometrical details of the fracture contact area and, due to the oscillatory fluid pressure diffusion process, is complex‐valued and frequency‐dependent. By using laboratory measurements of stress‐induced changes in the fracture contact area, we relate seismic attenuation and dispersion to the effective stress. The corresponding results do indeed indicate that seismic attenuation and phase velocity may constitute useful attributes to constrain the effective stress. Alternatively, knowledge of the effective stress may help to identify the regions in which wave induced fluid flow is expected to be the dominant attenuation mechanism.  相似文献   

16.
Saturation of porous rocks with a mixture of two fluids has a substantial effect on seismic‐wave propagation. In particular, partial saturation causes significant attenuation and dispersion of the propagating waves due to the mechanism of wave‐induced fluid‐flow. Such flow arises when a passing wave induces different fluid pressures in regions of rock saturated by different fluids. Most models of attenuation and dispersion due to mesoscopic heterogeneities imply that fluid heterogeneities are distributed in a regular way. However, recent experimental studies show that mesoscopic heterogeneities have less idealized distributions and that the distribution itself affects attenuation and dispersion. Based on an approximation for the coherent wavefield in random porous media, we develop a model which assumes a continuous distribution of fluid heterogeneities. As this continuous random media approach assumes that there will be a distribution of different patch sizes, it is expected to be better suited to modelling experimental data. We also show how to relate the random functions to experimentally measurable parameters.  相似文献   

17.
A nonlinear model for single-phase fluid flow in slightly compressible porous media is presented and solved approximately. The model assumes state equations for density, porosity, viscosity and permeability that are exponential functions of the fluid (either gas or liquid) pressure. The governing equation is transformed into a nonlinear diffusion equation. It is solved for a semi-infinite domain for either constant pressure or constant flux boundary conditions at the surface. The solutions obtained, although approximate, are extremely accurate as demonstrated by comparisons with numerical results. Predictions for the surface pressure resulting from a constant flux into a porous medium are compared with published experimental data.  相似文献   

18.
地震波传播激发的不同尺度的流固相对运动(宏观、中观和微观)是许多沉积岩地层中地震波频散和衰减的主要原因,然而野外观测和试验测量都难以对非均匀多孔介质孔隙压力弛豫物理过程进行精细刻画.通过数字岩石物理技术,本文建立了三个典型的数字岩心分别用于表征孔隙结构、岩石骨架和斑状饱和流体引起的非均质性,利用动态应力应变模拟技术计算数字岩心的位移和孔隙流体增量图像.通过分析和比较三个数字岩心的位移和孔隙压力增量图像,细致刻画了发生于非均匀含流体多孔介质内的宏观、中观和微观尺度的流固相对运动:1)宏观尺度的波致孔隙流体流动导致波长尺度上数字岩心不同区域的孔隙压力和位移差异;2)中观尺度的流体流动发生在软层与硬层之间、气层与液层之间;3)微观尺度的流体流动发生在孔隙内部或相邻孔隙之间.数值模拟试验也证明基于数字岩心的动态应力应变模拟技术可以从微观尺度上更好的理解波致孔隙流体流动发生的物理机理,从而为建立岩石骨架、孔隙流体、孔隙结构非均质性和弹性波频散-衰减特征的映射关系奠定基础.  相似文献   

19.
Intrinsic wave attenuation at seismic frequencies is strongly dependent on rock permeability, fluid properties, and saturation. However, in order to use attenuation as an attribute to extract information on rock/fluid properties from seismic data, experimental studies on attenuation are necessary for a better understanding of physical mechanisms that are dominant at those frequencies. An appropriate laboratory methodology to measure attenuation at seismic frequencies is the forced oscillation method, but technical challenges kept this technique from being widely used. There is a need for the standardization of devices employing this method, and a comparison of existing setups is a step towards it. Here we summarize the apparatuses based on the forced oscillation method that were built in the last 30 years and were used to measure frequency‐dependent attenuation in fluid‐saturated and/or dry reservoir rocks under small strains (10?8–10?5). We list and discuss important technical aspects to be taken into account when working with these devices or in the course of designing a new one. We also present a summary of the attenuation measurements in reservoir rock samples performed with these apparatuses so far.  相似文献   

20.
Streaming potential is the result of coupling between a fluid flow and an electric current in porous rocks. The modified Helmholtz–Smoluchowski equation derived for capillary tubes is mostly used to determine the streaming potential coefficient of porous media. However, to the best of our knowledge, the fractal geometry theory is not yet applied to analyse the streaming potential in porous media. In this article, a fractal model for the streaming potential coefficient in porous media is developed based on the fractal theory of porous media and on the streaming potential in a capillary. The proposed model is expressed in terms of the zeta potential at the solid?liquid interface, the minimum and maximum pore/capillary radii, the fractal dimension, and the porosity of porous media. The model is also examined by using another capillary size distribution available in published articles. The results obtained from the model using two different capillary size distributions are in good agreement with each other. The model predictions are then compared with experimental data in the literature and those based on the modified Helmholtz–Smoluchowski equation. It is shown that the predictions from the proposed fractal model are in good agreement with experimental data. In addition, the proposed model is able to reproduce the same result as the Helmholtz–Smoluchowski equation, particularly for high fluid conductivity or large grain diameters. Other factors influencing the streaming potential coefficient in porous media are also analysed.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号