首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Elastic wave propagation and attenuation in porous rock layers with oriented sets of fractures, especially in carbonate reservoirs, are anisotropic owing to fracture sealing, fracture size, fracture density, filling fluid, and fracture strike orientation. To address this problem, we adopt the Chapman effective medium model and carry out numerical experiments to assess the variation in P-wave velocity and attenuation, and the shear-wave splitting anisotropy with the frequency and azimuth of the incident wave. The results suggest that velocity, attenuation, and anisotropy vary as function of azimuth and frequency. The azimuths of the minimum attenuation and maximum P-wave velocity are nearly coincident with the average strike of the two sets of open fractures. P-wave velocity is greater in sealed fractures than open fractures, whereas the attenuation of energy and anisotropy is stronger in open fractures than sealed fractures. For fractures of different sizes, the maximum velocity together with the minimum attenuation correspond to the average orientation of the fracture sets. Small fractures affect the wave propagation less. Azimuth-dependent anisotropy is low and varies more than the other attributes. Fracture density strongly affects the P-wave velocity, attenuation, and shear-wave anisotropy. The attenuation is more sensitive to the variation of fracture size than that of velocity and anisotropy. In the seismic frequency band, the effect of oil and gas saturation on attenuation is very different from that for brine saturation and varies weakly over azimuth. It is demonstrated that for two sets of fractures with the same density, the fast shear-wave polarization angle is almost linearly related with the orientation of one of the fracture sets.  相似文献   

2.
结合有限差分方法和等效介质理论,模拟了离散分布裂缝介质中地震波的传播. 基于等效介质理论,利用二维有限差分实现封闭裂缝的离散分布;裂缝可以处理成固体岩石中的高度柔性界面,并可以用线性滑动或者位移间断模型进行裂缝的物理描述. 对于含有多组裂隙的破裂固体,其有效柔度可以认为是固体骨架背景柔度和裂缝附加柔度之和. 在一阶近似条件下,固体骨架和裂缝参数可以通过有效各向异性系数联系起来,有效各向异性系数决定了各向异性（裂缝效应）对于地震波传播的影响. 通过与射线理论方法的对比检验,说明本文提出的模拟方法的有效性,并通过几个数值算例说明本方法可有效模拟不同的裂缝分布效应. 结果表明,即使在裂缝密度很小的情况下,具有相同裂缝密度的不同的空间分布可以产生不同的波场特征. 同时,也验证了不同裂缝尺度对波长的不同影响,以及裂缝尺度具有幂率分布（分形）时,尺度对波场的影响. 最后得出结论：在运用建立在等效介质理论基础上的地震各向异性概念来描述裂缝固体的特征时,要倍加小心,等效介质理论中尚未合理处理的裂缝尺度和空间分布对波的传播特征具有重要的影响.  相似文献   

3.
为了检测定向裂隙介质中横波分裂的方位属性特征,分析地震属性随裂隙密度和方位变化,采用人工吸收边界和反周期扩展边界,用伪谱法获得不同裂隙密度和不同方位地质模型三分量地面记录;应用时频分析和剪切波偏振分析研究由于裂隙方位和密度引起的横波分裂.结果显示,裂隙密度和方位决定着横波分裂的时差和偏振.快慢横波的延迟时间随裂隙密度增大而增加,不同方位相同裂隙密度的横波分裂时差有微小的变化.在45°方位检测时间延迟时间最大.通过时频分析,可以看到不同方位的瞬时主频有显著的变化,在横波分裂处瞬时主频有明显变化.因此,瞬时主频和快横波的偏振以及延迟时间可以作为裂隙方位和密度的指示.  相似文献   

4.
准确模拟TTI介质中弹性波的传播是研究地震各向异性、AVO反演的基础. 在二维加权近似解析离散化(WNAD)算法的基础上, 本文发展的并行WNAD算法是一种研究三维横向各向同性(TI)介质中弹性波传播的、快速高效的数值模拟方法. 我们首先介绍三维WNAD方法的构造过程, 然后与经典的差分格式--交错网格(SG)算法进行了比较. 理论分析和数值算例表明, WNAD算法比交错网格算法更适合在高性能计算机上进行大规模弹性波场模拟. 同时, 本文利用并行的WNAD方法研究了弹性波在TTI介质中的传播规律, 观测了TI介质中弹性波传播的重要特征:横波分离、体波耦合和速度各向异性等. 在TTI介质分界面处, 弹性波产生更加复杂的折射、反射和波型转化, 使得波场非常复杂, 研究和辨别不同类型的波能够加深我们对由裂隙诱导的各向异性介质的认识.  相似文献   

5.
Measurements of seismic anisotropy in fractured rock are used at present to deduce information about the fracture orientation and the spatial distribution of fracture intensity. Analysis of the data is based upon equivalent-medium theories that describe the elastic response of a rock containing cracks or fractures in the long-wavelength limit. Conventional models assume frequency independence and cannot distinguish between microcracks and macrofractures. The latter, however, control the fluid flow in many subsurface reservoirs. Therefore, the fracture size is essential information for reservoir engineers. In this study we apply a new equivalent-medium theory that models frequency-dependent anisotropy and is sensitive to the length scale of fractures. The model considers velocity dispersion and attenuation due to a squirt-flow mechanism at two different scales: the grain scale (microcracks and equant matrix porosity) and formation-scale fractures. The theory is first tested and calibrated against published laboratory data. Then we present the analysis and modelling of frequency-dependent shear-wave splitting in multicomponent VSP data from a tight gas reservoir. We invert for fracture density and fracture size from the frequency dependence of the time delay between split shear waves. The derived fracture length matches independent observations from borehole data.  相似文献   

6.
We present numerical modeling of SH-wave propagation for the recently proposed whole Moon model and try to improve our understanding of lunar seismic wave propagation. We use a hybrid PSM/FDM method on staggered grids to solve the wave equations and implement the calculation on a parallel PC cluster to improve the computing efficiency. Features of global SH-wave propagation are firstly discussed for a 100-km shallow and 900-km deep moonquakes, respectively. Effects of frequency range and lateral variation of crust thickness are then investigated with various models. Our synthetic waveforms are finally compared with observed Apollo data to show the features of wave propagation that were produced by our model and those not reproduced by our models. Our numerical modeling show that the low-velocity upper crust plays significant role in the development of reverberating wave trains. Increasing frequency enhances the strength and duration of the reverberations. Surface multiples dominate wavefields for shallow event. Core–mantle reflections can be clearly identified for deep event at low frequency. The layered whole Moon model and the low-velocity upper crust produce the reverberating wave trains following each phases consistent with observation. However, more realistic Moon model should be considered in order to explain the strong and slow decay scattering between various phases shown on observation data.  相似文献   

7.
Natural fractures in hydrocarbon reservoirs can cause significant seismic attenuation and dispersion due to wave induced fluid flow between pores and fractures. We present two theoretical models explicitly based on the solution of Biot's equations of poroelasticity. The first model considers fractures as planes of weakness (or highly compliant and very thin layers) of infinite extent. In the second model fractures are modelled as thin penny-shaped voids of finite radius. In both models attenuation is a result of conversion of the incident compressional wave energy into the diffusive Biot slow wave at the fracture surface and exhibits a typical relaxation peak around a normalized frequency of about 1. This corresponds to a frequency where the fluid diffusion length is of the order of crack spacing for the first model and the crack diameter for the second. This is consistent with an intuitive understanding of the nature of attenuation: when fractures are closely and regularly spaced, the Biot's slow waves produced by cracks interfere with each other, with the interference pattern controlled by the fracture spacing. Conversely, if fractures are of finite length, which is smaller than spacing, then fractures act as independent scatterers and the attenuation resembles the pattern of scattering by isolated cracks. An approximate mathematical approach based on the use of a branching function gives a unified analytical framework for both models.  相似文献   

8.
In fractured reservoirs, seismic wave velocity and amplitude depend on frequency and incidence angle. Frequency dependence is believed to be principally caused by the wave‐induced flow of pore fluid at the mesoscopic scale. In recent years, two particular phenomena, i.e., patchy saturation and flow between fractures and pores, have been identified as significant mechanisms of wave‐induced flow. However, these two phenomena are studied separately. Recently, a unified model has been proposed for a porous rock with a set of aligned fractures, with pores and fractures filled with two different fluids. Existing models treat waves propagating perpendicular to the fractures. In this paper, we extend the model to all propagation angles by assuming that the flow direction is perpendicular to the layering plane and is independent of the loading direction. We first consider the limiting cases through poroelastic Backus averaging, and then we obtain the five complex and frequency‐dependent stiffness values of the equivalent transversely isotropic medium as a function of the frequency. The numerical results show that, when the bulk modulus of the fracture‐filling fluid is relatively large, the dispersion and attenuation of P‐waves are mainly caused by fractures, and the values decrease as angles increase, almost vanishing when the incidence angle is 90° (propagation parallel to the fracture plane). While the bulk modulus of fluid in fractures is much smaller than that of matrix pores, the attenuation due to the “partial saturation” mechanism makes the fluid flow from pores into fractures, which is almost independent of the incidence angle.  相似文献   

9.
We present the analysis of a multi-azimuth vertical seismic profiling data set that has been acquired in a tight gas field with the objective of characterizing fracture distributions using seismic anisotropy. We investigate different measurements of anisotropy, which are shear-wave splitting, P-wave traveltime anisotropy and azimuthal amplitude variation with offset. We find that for our field case shear-wave splitting is the most robust measure of azimuthal anisotropy, which is clearly observed over two distinct intervals in the target. We compare the results of the vertical seismic profiling analysis with other borehole data from the same well. Cross-dipole sonic and Formation MicroImager data from the reservoir section suggest that no open fractures intersect the well or are present within half a metre of the borehole wall. Furthermore, a detailed dispersion analysis of the sonic scanner data provides no indication of stress-induced seismic anisotropy along the logged borehole section. We therefore explain the azimuthal anisotropy measured in the vertical seismic profiling data with a model that contains discrete fracture corridors, which do not intersect the well itself but lie within the vertical seismic profiling investigation radius. We show that such a model can reproduce some basic characteristics of azimuthal anisotropy observed in the vertical seismic profiling data. The model is also consistent with well test data that suggest the presence of a fracture corridor away from the well. With this study we demonstrate the necessity of integrating different data types that investigate different scales of rock volume and can provide complementary information for understanding the characteristics of fracture networks in the subsurface.  相似文献   

10.
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.  相似文献   

11.
This paper aims to improve current understanding of the subsurface fracture system in the Coso geothermal field, located in east-central California. The Coso reservoir is in active economic development, so that knowledge of the subsurface fracture system is of vital importance for an accurate evaluation of its geothermal potential and day-to-day production. To detect the geometry and density of fracture systems we applied the shear-wave splitting technique to a large number of high-quality seismograms from local microearthquakes recorded by a permanent, 16-station, down-hole, 3-component seismic array running at 480 samples/s. The analysis of shear-wave splitting (seismic birefringence) provides parameters directly related to the strike of the subsurface fractures and their density (number of cracks per unit volume), and, consequently, is an important technique to outline zones of high permeability. Three major fracture directions N10–30W, N0–20E, and N40–50E, of which the first and the second are the most prominent, were identified from the seismograms recorded by the 16-station down-hole array. All orientations are consistent with the known strike of local sets of faults and fractures in local wells and at the surface, as well as with previous analyses of seismic anisotropy in the region. The high quality of the recordings has allowed us to launch an unprecedented investigation into the characteristics of the temporal variations in crack polarization and crack density in a producing geothermal environment. Preliminary results point to significant temporal changes in shear-wave time delays, probably influenced by temporal changes in crack density within a period of 5 years (1996–2000). They are tentatively interpreted as due to a local 3% increase in shear-wave velocity in the southwestern part of the field during 1999.  相似文献   

12.
Kirchhoff–Helmholtz (KH) theory is extended to synthesize two-way elastic wave propagation in 3D laterally heterogeneous, anisotropic media. I have developed and tested numerically a specialized algorithm for the generation of three-component synthetic seismograms in multi-layered isotropic and transversely isotropic (TI) media with dipping interfaces and tilted axes of symmetry. This algorithm can be applied to vertical seismic profile (VSP) geometries and works well when the receiver is located near the reflector interface. It is superior to ray methods in predicting elliptical polarization effects observed on radial and transverse components. The algorithm is used to study converted-wave propagation for determining fracture-related shear-wave anisotropy in realistic reservoir models. Results show that all wavefront attributes are strongly affected by the anisotropy. However, it is necessary to resolve a trade-off between the effects of fractures and formation dip prior to converted-wave interpretation. These results provide some assurance that the present scheme is sufficiently versatile to handle shear wave behaviour due to various generalized rays propagating in complex geological models.  相似文献   

13.
We present a finite difference (FD) method for the simulation of seismic wave fields in fractured medium with an irregular (non-flat) free surface which is beneficial for interpreting exploration data acquired in mountainous regions. Fractures are introduced through the Coates-Schoenberg approach into the FD scheme which leads to local anisotropic properties of the media where fractures are embedded. To implement surface topography, we take advantage of the boundary-conforming grid and map a rectangular grid onto a curved one. We use a stable and explicit second-order accurate finite difference scheme to discretize the elastic wave equations (in a curvilinear coordinate system) in a 2D heterogeneous transversely isotropic medium with a horizontal axis of symmetry (HTI). Efficiency tests performed by different numerical experiments clearly illustrate the influence of an irregular free surface on seismic wave propagation in fractured media which may be significant to mountain seismic exploration. The tests also illustrate that the scattered waves induced by the tips of the fracture are re-scattered by the features of the free surface topography. The scattered waves provoked by the topography are re-scattered by the fractures, especially Rayleigh wave scattering whose amplitudes are much larger than others and making it very difficult to identify effective information from the fractures.  相似文献   

14.
Certain crack-influence parameters of Sayers and Kachanov are shown to be directly related to Thomsen's weak-anisotropy seismic parameters for fractured reservoirs when the crack/fracture density is small enough. These results are then applied to the problem of seismic wave propagation in polar reservoirs, i.e., those anisotropic reservoirs having two axes that are equivalent but distinct from the third axis), especially for horizontal transversely isotropic seismic wave symmetry due to the presence of aligned vertical fractures and resulting in azimuthal seismic wave symmetry at the Earth's surface. The approach presented suggests one method of inverting for fracture density from wave speed data. A significant fraction of the technical effort in the paper is devoted to showing how to predict the angular location of the true peak (or trough) of the quasi-SV-wave for polar media and especially how this peak is related to another angle that is very easy to compute. The axis of symmetry is always treated here as the x 3-axis for either vertical transversely isotropic symmetry (due, for example, to horizontal cracks), or horizontal transversely isotropic symmetry (due to aligned vertical cracks). Then, the meaning of the stiffnesses is derived from the fracture analysis in the same way for vertical transversely isotropic and horizontal transversely isotropic media, but for horizontal transverse isotropy the wave speeds relative to the Earth's surface are shifted by  90o  in the plane perpendicular to the aligned vertical fractures. Skempton's poroelastic coefficient B is used as a general means of quantifying the effects of fluids inside the fractures. Explicit Biot-Gassmann-consistent formulas for Thomsen's parameters are also obtained for either drained or undrained fractures resulting in either vertical transversely isotropic or horizontal transversely isotropic symmetry of the reservoir.  相似文献   

15.
Average elastic properties of a fluid‐saturated fractured rock are discussed in association with the extremely slow and dispersive Krauklis wave propagation within individual fractures. The presence of the Krauklis wave increases P‐wave velocity dispersion and attenuation with decreasing frequency. Different laws (exponential, power, fractal, and gamma laws) of distribution of the fracture length within the rock show more velocity dispersion and attenuation of the P‐wave for greater fracture density, particularly at low seismic frequencies. The results exhibit a remarkable difference in the P‐wave reflection coefficient for frequency and angular dependency from the fractured layer in comparison with the homogeneous layer. The biggest variation in behaviour of the reflection coefficient versus incident angle is observed at low seismic frequencies. The proposed approach and results of calculations allow an interpretation of abnormal velocity dispersion, high attenuation, and special behaviour of reflection coefficients versus frequency and angle of incidence as the indicators of fractures.  相似文献   

16.
Many investigations of the propagation elastic waves within the earth require a technique for producing synthetic seismograms which is capable of modelling 3D propagation effects. Ray methods are an excellent option for these problems, because they can be made fully 3D and allow a relatively quick and flexible computation of synthetic seismograms. However, the two point problem of finding the ray which connects exactly a specific source and receiver, may still be difficult and time consuming. Therefore, application of the paraxial method, which allows extrapolation of the information on a given ray to nearby receiver locations, is very valuable. With this approach, great savings in computation time and significant simplification of computer codes are possible. We investigate the application of the paraxial ray method to two problems in which the effects of 3D seismic wave propagation are important. The first is a model of a reef structure. In this case, we consider synthetic seismograms for a VSP experimental configuration. When the SV source and well are located along the axis of the reef, only 2D propagation effects are observed. If the source-receiver plane is located to the side of the reef, however, the 3D shape of the reef causes significant amplitudes to be predicted for shear-wave arrivals on the transverse component of the synthetic seismograms. The second example is a ID, layered earth model, but it includes two layers which are azimuthally anisotropic due to the presence of aligned, vertical fractures. This anisotropy leads to 3D raypaths. Synthetic seismographs are presented for a cross-hole geometry both for an equivalent isotropic model and for the direction parallel to the fractures and at an angle of 45° to the fractures in the anisotropic case. These synthetics show that the differences between the isotropic case and the case for source and receivers aligned with the plane of the fractures are small and subtle. On the other hand, the predictions for the direction at 45° to the cracks show shear-wave splitting and significant transverse component signal. These results have important implications for both modelling and for applications such as tomography. It is clear that for some of the cases considered, a 2D algorithm will lead to errors in interpretation of data. In addition, ray-based tomographic techniques will have great difficulty in obtaining a well-defined 2D planar image when the signals are propagating in regions outside the image plane.  相似文献   

17.
裂缝广泛分布于地球介质中并且具有多尺度的特点,裂缝尺度对于油气勘探和开发有着重要的意义.本文制作了一组含不同长度裂缝的人工岩样,其中三块含裂缝岩样中的裂缝直径分别为2 mm、3 mm和4 mm,裂缝的厚度都约为0.06 mm,裂缝密度大致相同(分别为4.8%、4.86%和4.86%).在岩样含水的条件下测试不同方向上的纵横波速度,实验结果表明,虽然三块裂缝岩样中的裂缝密度大致相同,但是含不同直径裂缝岩样的纵横波速度存在差异.在各个方向上,含数量众多的小尺度裂缝的岩样中纵横波速度都明显低于含少量的大尺度裂缝的岩样中纵横波速度.尤其是对纵波速度和SV波速度,在不同尺度裂缝岩样中的差异更明显.在含数量多的小尺度裂缝的岩样中纵波各向异性和横波各向异性最高,而含少量的大尺度的裂缝的岩样中的纵波各向异性和横波各向异性较低.实验测量结果与Hudson理论模型预测结果进行了对比分析,结果发现Hudson理论考虑到了裂缝尺度对纵波速度和纵波各向异性的影响,但是忽略了其对横波速度和横波各向异性的影响.  相似文献   

18.
Seismic attenuation mechanisms receive increasing attention for the characterization of fractured formations because of their inherent sensitivity to the hydraulic and elastic properties of the probed media. Attenuation has been successfully inferred from seismic data in the past, but linking these estimates to intrinsic rock physical properties remains challenging. A reason for these difficulties in fluid-saturated fractured porous media is that several mechanisms can cause attenuation and may interfere with each other. These mechanisms notably comprise pressure diffusion phenomena and dynamic effects, such as scattering, as well as Biot's so-called intrinsic attenuation mechanism. Understanding the interplay between these mechanisms is therefore an essential step for estimating fracture properties from seismic measurements. In order to do this, we perform a comparative study involving wave propagation modelling in a transmission set-up based on Biot's low-frequency dynamic equations and numerical upscaling based on Biot's consolidation equations. The former captures all aforementioned attenuation mechanisms and their interference, whereas the latter only accounts for pressure diffusion phenomena. A comparison of the results from both methods therefore allows to distinguish between dynamic and pressure diffusion phenomena and to shed light on their interference. To this end, we consider a range of canonical models with randomly distributed vertical and/or horizontal fractures. We observe that scattering attenuation strongly interferes with pressure diffusion phenomena, since the latter affect the elastic contrasts between fractures and their embedding background. Our results also demonstrate that it is essential to account for amplitude reductions due to transmission losses to allow for an adequate estimation of the intrinsic attenuation of fractured media. The effects of Biot's intrinsic mechanism are rather small for the models considered in this study.  相似文献   

19.
利用新方法制作出含可控裂缝的双孔隙人工砂岩物理模型,具有与天然岩石更为接近的矿物成分、孔隙结构和胶结方式,其中裂缝密度、裂缝尺寸和裂缝张开度等裂缝参数可以控制以得到实验所需要的裂缝参数,岩样具有真实的孔隙和裂缝空间并可以在不同饱和流体状态下研究流体性质对于裂缝介质性质的影响.本次实验制作出一组具有不同裂缝密度的含裂缝人工岩样,对岩样利用SEM扫描电镜分析可以看到真实的孔隙结构和符合我们要求的裂缝参数,岩样被加工成八面棱柱以测量不同方向上弹性波传播的速度,用0.5 MHz的换能器使用透射法测量在饱和空气和饱和水条件下各个样品不同方向上的纵横波速度,并得出纵横波速度、横波分裂系数和纵横波各向异性强度受裂缝密度和饱和流体的影响.研究发现流体对于纵波速度和纵波各向异性强度的影响较强,而横波速度、横波分裂系数和横波各向异性强度受饱和流体的影响不大,但是对裂缝密度的变化更敏感.  相似文献   

20.
建立包含震源、沉积盆地和表层低速沉积层的二维模型,采用交错网格有限差分/伪谱混合方法求解地震波传播,讨论沉积层厚度和速度对地震地面运动的作用。结果表明:沉积层内产生的地震波的多重反射以及转换会引起地面运动持续时间的延长,它们的相干叠加会造成地面运动峰值的放大;随着沉积层速度的增加,多重反射与转换波的能量减小,地面运动持续时间减小,但是不同速度或者不同厚度的低速层模型均显示出一致的地面运动峰值放大特征。结果说明,在包含震源、沉积盆地和沉积层的模型中,沉积层对地面运动的作用机理更复杂。在实际应用中有必要同时考虑这些因素的综合作用。  相似文献   

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

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