共查询到20条相似文献，搜索用时 15 毫秒
1.
Knowledge about the spatial distribution of the fracture density and the azimuthal fracture orientation can greatly help in optimizing production from fractured reservoirs. Frequency-dependent seismic velocity and attenuation anisotropy data contain information about the fractures present in the reservoir. In this study, we use the measurements of velocity and attenuation anisotropy data corresponding to different seismic frequencies and azimuths to infer information about the multiple fracture sets present in the reservoir. We consider a reservoir model with two sets of vertical fractures characterized by unknown azimuthal fracture orientations and fracture densities. Frequency-dependent seismic velocity and attenuation anisotropy data is computed using the effective viscoelastic stiffness tensor and solving the Christoffel equation. A Bayesian inversion method is then applied to measurements of velocity and attenuation anisotropy data corresponding to different seismic frequencies and azimuth to estimate the azimuthal fracture orientations and the fracture densities, as well as their uncertainties. Our numerical examples suggest that velocity anisotropy data alone cannot recover the unknown fracture parameters. However, an improved estimation of the unknown fracture parameters can be obtained by joint inversion of velocity and attenuation anisotropy data. 相似文献
2.
Gui-Hong Guo Jian-Ping Yan Zhi Zhang José Badal Jian-Wu Cheng Shuang-Hu Shi Ya-Wei Ma 《应用地球物理》2018,15(2):299-310
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. 相似文献
3.
裂缝广泛分布于地球介质中并且具有多尺度的特点,裂缝尺度对于油气勘探和开发有着重要的意义.本文制作了一组含不同长度裂缝的人工岩样,其中三块含裂缝岩样中的裂缝直径分别为2 mm、3 mm和4 mm,裂缝的厚度都约为0.06 mm,裂缝密度大致相同(分别为4.8%、4.86%和4.86%).在岩样含水的条件下测试不同方向上的纵横波速度,实验结果表明,虽然三块裂缝岩样中的裂缝密度大致相同,但是含不同直径裂缝岩样的纵横波速度存在差异.在各个方向上,含数量众多的小尺度裂缝的岩样中纵横波速度都明显低于含少量的大尺度裂缝的岩样中纵横波速度.尤其是对纵波速度和SV波速度,在不同尺度裂缝岩样中的差异更明显.在含数量多的小尺度裂缝的岩样中纵波各向异性和横波各向异性最高,而含少量的大尺度的裂缝的岩样中的纵波各向异性和横波各向异性较低.实验测量结果与Hudson理论模型预测结果进行了对比分析,结果发现Hudson理论考虑到了裂缝尺度对纵波速度和纵波各向异性的影响,但是忽略了其对横波速度和横波各向异性的影响. 相似文献
4.
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. 相似文献
5.
Measuring changes in fracture properties from temporal variations in anisotropic attenuation of microseismic waveforms 
下载免费PDF全文
下载免费PDF全文 We investigate fracture‐induced attenuation anisotropy in a cluster of events from a microseismic dataset acquired during hydraulic fracture stimulation. The dataset contains 888 events of magnitude ?3.0 to 0.0. We use a log‐spectral‐amplitude‐ratio method to estimate change in over a half‐hour time period where fluid is being injected and an increase in fracturing from S‐wave splitting analysis has been previously inferred. A Pearson's correlation analysis is used to assess whether or not changes in attenuation with time are statistically significant. P‐waves show no systematic change in during this time. In contrast, S‐waves polarised perpendicular to the fractures show a clear and statistically significant increase with time, whereas S‐waves polarised parallel to the fractures show a weak negative trend. We also compare between the two S‐waves, finding an increase in with time. A poroelastic rock physics model of fracture‐induced attenuation anisotropy is used to interpret the results. This model suggests that the observed changes in t* are related to an increase in fracture density of up to 0.04. This is much higher than previous estimates of 0.025 ± 0.002 based on S‐wave velocity anisotropy, but there is considerably more scatter in the attenuation measurements. This could be due to the added sensitivity of attenuation measurement to non‐aligned fractures, fracture shape, and fluid properties. Nevertheless, this pilot study shows that attenuation measurements are sensitive to fracture properties such as fracture density and aspect ratio. 相似文献
6.
本文首先由Christoffel方程推导出黏弹性EDA介质中均匀、 非均匀P波、 SV波和SH波的相速度表达式, 然后参照极端各向异性介质的相关计算方法, 推导出EDA介质中均匀、 非均匀地震波相衰减系数和群衰减系数的表达式, 并通过数值计算分析了相速度、 相衰减系数、 群衰减系数与裂隙方位的关系. 结果表明: 均匀介质中SH波的相速度和相衰减系数均可指示裂隙的走向; 非均匀介质中SH波相衰减系数随非均匀角的增大而增大, 且其对称轴与介质对称轴的夹角也相应增加; 由于地震波振幅的衰减随岩石物理性质的变化比地震波速度的变化更为灵敏, 而且携带了更多的岩石物理性质信息, 因此可用来探明裂隙走向、 密度及含水特性, 进而应用于预测、 预防地下工程地质灾害事故. 相似文献
7.
基于任意空间取向TI介质(简称ATI)中体波速度和偏振解析解,通过模型数值计算给出ATI介质中体波群速度和相速度的变化特征,说明TI空间取向与测线方位对速度的影响。研究表明,体波群速度图案和相速度图案相对TI对称轴固定,随TI对称轴倾角及其相对测线方位角的变化呈现出一定的对称性和重复性;可以针对ATI地区的地质情况,给出体波群速度和相速度变化图案,为进一步的理论研究提供便捷。此结果也可以直接用于VSP(垂直地震剖面)和井间地震资料的分析研究 相似文献
8.
非双曲（远偏移距）时差为各向异性介质中正、反演研究,特别是各向异性参数估计提供了重要信息.本文在任意空间取向TI（ATI）介质水平界面同类反射波四次时差系数(A4)精确解的基础上,进一步讨论我们推导得出的ATI介质中四次时差系数解析近似解,比较随CMP测线方位变化的近似解与精确解之间的差别,为利用近似解来解析研究ATI介质中非双曲时距以及参数反演提供有价值的信息．结合实际岩性资料的数值研究表明,ATI条件下四次时差系数近似解与精确解之间存在差别,不仅表现在A4数的大小及符号特征上,更突出地表现在A4系数随方位的变化特征上;在强各向异性条件下,近似解相比精确解存在较大误差．但在各向异性参数满足0<ε-δ<0.15、｜δ｜<0.20的情况下,对于TI对称轴的特殊倾角范围(75°～80°),近似解与精确解的差别很小,可用近似解进行各向异性观测解释及参数反演． 相似文献
9.
O.H. Al‐Harrasi A. Al‐Anboori A. Wüstefeld J.‐M. Kendall 《Geophysical Prospecting》2011,59(2):227-243
The study of seismic anisotropy in exploration seismology is gaining interest as it provides valuable information about reservoir properties and stress directions. In this study we estimate anisotropy in a petroleum field in Oman using observations of shear‐wave splitting from microseismic data. The data set was recorded by arrays of borehole geophones deployed in five wells. We analyse nearly 3400 microearthquakes, yielding around 8500 shear‐wave splitting measurements. Stringent quality control reduces the number of reliable measurements to 325. Shear‐wave splitting modelling in a range of rock models is then used to guide the interpretation. The difference between the fast and slow shear‐wave velocities along the raypath in the field ranges between 0–10% and it is controlled both by lithology and proximity to the NE‐SW trending graben fault system that cuts the field formations. The anisotropy is interpreted in terms of aligned fractures or cracks superimposed on an intrinsic vertical transversely isotropic (VTI) rock fabric. The highest magnitudes of anisotropy are within the highly fractured uppermost unit of the Natih carbonate reservoir. Anisotropy decreases with depth, with the lowest magnitudes found in the deep part of the Natih carbonate formation. Moderate amounts of anisotropy are found in the shale cap rock. Anisotropy also varies laterally with the highest anisotropy occurring either side of the south‐eastern graben fault. The predominant fracture strikes, inferred from the fast shear‐wave polarizations, are consistent with the trends of the main faults (NE‐SW and NW‐SE). The majority of observations indicate subvertical fracture dip (>70° ). Cumulatively, these observations show how studies of shear‐wave splitting using microseismic data can be used to characterize fractures, important information for the exploitation of many reservoirs. 相似文献
10.
P波通过裂缝地层传播时,其响应会随着地面观测方位的变化而变化,表现出明显的方位各向异性.这种随方位变化的裂缝响应信息(如振幅等)广泛存在于叠前方位一偏移距二维空间域中,检测此域中的裂缝响应,通过椭圆拟合的方法识别各向异性的存在,进而计算裂缝的方向——椭圆长轴的方向,和裂缝密度——椭圆长轴与短轴的比.本文研究出了一种基于常规叠前地震反射振幅来检测裂缝的三维叠前裂缝储层预测方法.它采用面元组合的方法扩大方位角和偏移距范围;采用多级数据控制技术获取所需的信息;采用稳妥的技术步骤来预测.通过大港油田实际工区的应用,结果表明本文研究出的方法技术应用效果好、精度高、经济实用. 相似文献
11.
利用全方位P波属性进行裂缝检测的地震方法 总被引：2，自引：0，他引：2
胜利油田存在着一种特殊的裂隙性油气藏,即泥岩裂隙油气藏。由于该类裂隙储层的孔隙度很小,岩石物性参数变化不灵敏,并表现出很强的各向异性,因此,其勘探难度很大。到目前为止,国内外还没有一套成熟的地质、物探、测井及钻井等资料综合分析的裂隙方位、分布、密度的识别方法。本文在研究国外裂缝检测方法的基础上,提出了波阻抗随方位变化(IPVA)的新方法,结合罗家地区的泥岩裂隙地震地质特征,研究了利用多方位 P 波资料检测定向垂直裂缝的地震采集、处理和识别方法,对不同共中心点 CMP 位置的全方位 P 波资料在速度随方位变化(VVA)及振幅随方位角变化(RVA)研究的基础上,进行波阻抗随方位变化(IPVA)的研究,用于识别裂缝的分布、走向及密度。经罗家地区实际资料的应用,见到了初步的效果,证明该方法是潜力较大的高角度裂隙的定量检测方法。 相似文献
12.
流体在断裂和岩石骨架间的交换被认为是影响岩石弹性参数各向异性的主要原因,理论研究表明断裂尺度同样对弹性参数的各向异性也有影响.为了说明两者对各向异性影响以实现多尺度断裂裂隙的识别,本文在等效介质模型的基础上,应用数值分析的方法研究速度和衰减(1/Q)随多尺度断裂、频率和流体因子变化规律.结果表明介质弹性参数是频率依赖的,并且参数中存在衰减项,而这种频率依赖性与介质物性参数中的断裂尺度及流体性质存在一定的联系;当断裂定向分布时,参数结果显示为各向异性;不同断裂尺度具有不同的波速频散特性,剪切波分裂程度依赖于频率,断裂尺度起着控制作用,高频时对小尺度的敏感,低频段对大尺度敏感.在地震频段Thomsen参数随着频率的增大而减小,随着断裂尺寸的增大而减小.因此地震数据可能区分断裂和微裂隙引起各向异性,从而可测量断裂尺度. 相似文献
13.
水力压裂是页岩气开发过程中的核心增产技术,微地震则广泛用于压裂分析、水驱前缘监测和储层描述.微地震反演过程中,用于反演的速度模型往往基于测井、地震或标定炮资料构建,忽略了压裂过程中裂缝及孔隙流体压力变化对地层速度的影响.本文首先基于物质守恒、渗流理论和断裂力学模拟三维水力压裂过程,得到地下裂缝发育特征和孔隙压力分布.继而根据Coates-Schoenberg方法和裂缝柔量参数计算裂缝和孔隙压力对速度场的影响,得到压裂过程中的实时速度模型.最后利用三维射线追踪方法正演微地震走时和方位信息,并采用常规微地震定位方法反演震源位置及进行误差分析.数值模拟结果表明,检波器空间分布影响定位精度,常规方法的定位误差随射线路径在压裂带中传播距离增加而变大,且不同压裂阶段的多点反演法与单点极化法精度相当. 相似文献
14.
方位角度域共成像点道集能够客观反映地下介质的速度、各向异性参数异常以及振幅随角度变化（AVA）和裂缝信息。传统Kirchhoff PSTM通常输出偏移距域共成像点道集,对于速度分析、各向异性分析、AVA分析、裂缝识别等均存在诸多不便。本文提出了基于走时梯度的Kirchhoff叠前时间偏移全方位角度集输出方法并提出工业上切实可行的实现方案。通过走时场梯度计算波场传播方向矢量,形成能够反映观测系统参数和波场传播情况的全方位角度域共成像点道集。为了在大规模地震数据Kirchhoff积分叠前时间偏移中输出全方位角度道集,本文给出基于输入道方式的偏移实现方法,采用逐条inline线进行线偏移成像,从而大大降低了全方位角度道集输出对计算机内存的压力,显著提高了Kirchhoff积分时间偏移输出全方位角度道集的可行性。三维盐丘模型测试和海上某区块三维实际资料试验证明了本文方法的正确性。 相似文献
15.
Wide-azimuth seismic data can be used to derive anisotropic parameters on the subsurface by observing variation in subsurface seismic response along different azimuths. Layer-based high-resolution estimates of components of the subsurface anisotropic elastic tensor can be reconstructed by using wide-azimuth P-wave data by combining the kinematic information derived from anisotropic velocity analysis with dynamic information obtained from amplitude versus angle and azimuth analysis of wide-azimuth seismic data. Interval P-impedance, S-impedance and anisotropic parameters associated with anisotropic fracture media are being reconstructed using linearized analysis assuming horizontal transverse anisotropy symmetry. In this paper it is shown how additional assumptions, such as the rock model, can be used to reduce the degrees of freedom in the estimation problem and recover all five anisotropic parameters. Because the use of a rock model is needed, the derived elastic parameters are consistent with the rock model and are used to infer fractured rock properties using stochastic rock physics inversion. The inversion is based on stochastic rock physics modelling and maximum a posteriori estimate of both porosity and crack density parameters associated with the observed elastic parameters derived from both velocity and amplitude versus angle and azimuth analysis. While the focus of this study is on the use of P-wave reflection data, we also show how additional information such as shear wave splitting and/or anisotropic well log data can reduce the assumptions needed to derive elastic parameter and rock properties. 相似文献
16.
各向异性ATI介质剪切位错点源P波远场辐射 总被引：2，自引：0，他引：2
蔡晓刚 《CT理论与应用研究》2011,20(4):433-452
假设震源区为ATI介质,传播介质为全空间各向同性介质情况下,给出了剪切位错源P波远场辐射解析表达式,讨论了震源区各向异性对远场P波振幅和震源球的影响。结果表明,震源区裂隙密度,裂隙方位等参数的动态变化导致P波远场辐射花样也随之动态变化,同时也会造成非双力偶百分含量的动态变化。研究认为,在井中对中小地震的非双力偶机制进行... 相似文献
17.
Although it is believed that natural fracture sets predominantly have near‐vertical orientation, oblique stresses and some other mechanisms may tilt fractures away from the vertical. Here, we examine an effective medium produced by a single system of obliquely dipping rotationally invariant fractures embedded in a transversely isotropic with a vertical symmetry axis (VTI) background rock. This model is monoclinic with a vertical symmetry plane that coincides with the dip plane of the fractures. Multicomponent seismic data acquired over such a medium possess several distinct features that make it possible to estimate the fracture orientation. For example, the vertically propagating fast shear wave (and the fast converted PS‐wave) is typically polarized in the direction of the fracture strike. The normal‐moveout (NMO) ellipses of horizontal reflection events are co‐orientated with the dip and strike directions of the fractures, which provides an independent estimate of the fracture azimuth. However, the polarization vector of the slow shear wave at vertical incidence does not lie in the horizontal plane – an unusual phenomenon that can be used to evaluate fracture dip. Also, for oblique fractures the shear‐wave splitting coefficient at vertical incidence becomes dependent on fracture infill (saturation). A complete medium‐characterization procedure includes estimating the fracture compliances and orientation (dip and azimuth), as well as the Thomsen parameters of the VTI background. We demonstrate that both the fracture and background parameters can be obtained from multicomponent wide‐azimuth data using the vertical velocities and NMO ellipses of PP‐waves and two split SS‐waves (or the traveltimes of PS‐waves) reflected from horizontal interfaces. Numerical tests corroborate the accuracy and stability of the inversion algorithm based on the exact expressions for the vertical and NMO velocities. 相似文献
18.
The hydrocarbon industry is moving increasingly towards tight sandstone and shale gas resources – reservoirs that require fractures to be produced economically. Therefore, techniques that can identify sets of aligned fractures are becoming more important. Fracture identification is also important in the areas of coal bed methane production, carbon capture and storage (CCS), geothermal energy, nuclear waste storage and mining. In all these settings, stress and pore pressure changes induced by engineering activity can generate or reactivate faults and fractures. P‐ and S‐waves are emitted by such microseismic events, which can be recorded on downhole geophones. The presence of aligned fracture sets generates seismic anisotropy, which can be identified by measuring the splitting of the S‐waves emitted by microseismic events. The raypaths of the S‐waves will have an arbitrary orientation, controlled by the event and geophone locations, meaning that the anisotropy system may only be partly illuminated by the available arrivals. Therefore to reliably interpret such splitting measurements it is necessary to construct models that compare splitting observations with modelled values, allowing the best fitting rock physics parameters to be determined. Commonly, splitting measurements are inverted for one fracture set and rock fabrics with a vertical axis of symmetry. In this paper we address the challenge of identifying multiple aligned fracture sets using splitting measured on microseismic events. We analyse data from the Weyburn CCS‐EOR reservoir, which is known to have multiple fracture sets, and from a hydraulic fracture stimulation, where it is believed that only one set is present. We make splitting measurements on microseismic data recorded on downhole geophone arrays. Our inversion technique successfully discriminates between the single and multiple fracture cases and in all cases accurately identifies the strikes of fracture sets previously imaged using independent methods (borehole image logs, core samples, microseismic event locations). We also generate a synthetic example to highlight the pitfalls that can be encountered if it is assumed that only one fracture set is present when splitting data are interpreted, when in fact more than one fracture set is contributing to the anisotropy. 相似文献
19.
利用新方法制作出含可控裂缝的双孔隙人工砂岩物理模型,具有与天然岩石更为接近的矿物成分、孔隙结构和胶结方式,其中裂缝密度、裂缝尺寸和裂缝张开度等裂缝参数可以控制以得到实验所需要的裂缝参数,岩样具有真实的孔隙和裂缝空间并可以在不同饱和流体状态下研究流体性质对于裂缝介质性质的影响.本次实验制作出一组具有不同裂缝密度的含裂缝人工岩样,对岩样利用SEM扫描电镜分析可以看到真实的孔隙结构和符合我们要求的裂缝参数,岩样被加工成八面棱柱以测量不同方向上弹性波传播的速度,用0.5 MHz的换能器使用透射法测量在饱和空气和饱和水条件下各个样品不同方向上的纵横波速度,并得出纵横波速度、横波分裂系数和纵横波各向异性强度受裂缝密度和饱和流体的影响.研究发现流体对于纵波速度和纵波各向异性强度的影响较强,而横波速度、横波分裂系数和横波各向异性强度受饱和流体的影响不大,但是对裂缝密度的变化更敏感. 相似文献