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

2.
裂缝发育导致的地震波场各向异性,是裂缝型油气藏的一个重要特征.然而针对该类油气藏的叠前弹性阻抗技术,一直没有通过各向异性弹性阻抗公式建立与裂缝参数的直接关系.本文在回顾封闭平行硬币状裂缝模型和线性滑动模型的基础上,将裂缝填充物性质、分布密度与围岩的横纵波速度比,引入裂缝型HTI介质弹性阻抗公式,并对相应的归一化弹性阻抗响应特征进行模拟分析.分析表明,弹性阻抗受以上三参数的影响规律存在明显差异,其中含气裂隙介质随裂缝密度的变化要明显高于对应的流体裂隙介质,流体填充裂隙介质随横、纵波速度比的变化幅度要高于对应的含气裂隙介.该结论初步为裂缝型油气藏流体识别工作提供了依据.  相似文献   

3.
刘财  兰慧田  郭智奇  冯晅  鹿琪 《地球物理学报》2013,56(10):3461-3473
改进BISQ(Biot-Squirt)机制在不引入特征喷流长度的情况下,将含流体孔隙介质中Biot流动和喷射流动两种重要的力学机制有机地结合起来,且各相关参数具有明确物理意义和可实现性.本文将改进BISQ机制一维孔隙流体压力公式推广到三维具有水平对称轴横向各向同性介质(HTI介质)情况,结合裂缝各向异性理论,给出了基于改进BISQ机制的双相HTI介质模型及其二维三分量波传播方程,采用伪谱法求解该方程,进行了不同相界、不同频率以及双层地质结构情况下该类介质中波场的数值模拟与特征分析.数值模拟结果表明:伪谱法模拟精度高,压制网格频散效果好,可以得到高精度的波场快照和合成记录;基于改进BISQ机制的双相HTI介质模型兼具裂缝各向异性特征和孔隙弹性特征,其为从双相各向异性理论角度深入研究裂缝性储层的地震响应奠定了理论基础.  相似文献   

4.
裂缝型储层流体识别方法   总被引:2,自引:0,他引:2       下载免费PDF全文
裂缝型储层的描述包括预测裂缝分布特征和识别裂隙充填物.依据等效介质理论计算的纵波速度随裂缝密度的增大而减小.正演地震记录显示,裂缝介质含气时反射振幅最大,且变化程度比含油或含水时大.叠前方位AVO反演所得的各向异性梯度Bani与裂缝密度成正比,可用于描述有效裂缝发育强度.对于不同的裂缝密度,各向异性梯度Bani与各向同性梯度Biso的比值I(1)fluid近似为常数,且对流体敏感.经裂缝纵横比和背景介质拉梅常数修正后,流体因子Ifluid既不随纵横比变化,又不受背景介质的影响,是裂缝型储层敏感的流体识别因子.在塔里木盆地塔北哈拉哈塘地区热瓦普区块碳酸盐岩储层裂缝发育区域,运用该参数在井点处的流体识别效果与钻井结果一致.  相似文献   

5.
二维横各向同性弹性随机介质中的波场特征   总被引:9,自引:4,他引:5  
本文通过交错网格有限差分正演.模拟了平面地震波在二维横各向同性弹性随机介质模型中的传播及其自激自收时间记录.为研究横各向同性弹性随机介质模型中的波场特征,我们在五个不同的时间区段上,分别计算剖面的三个统计特征(横向中心频率、纵向中心频率、波场能量相对值).这样,对应每一个横各向同性弹性随机介质模型.均可计算得到15个不同的波场特征量.我们通过在二维横各向同性弹性随机介质中的正演模拟.研究当自相关长度以及介质的各向异性系数变化时,对应的上述波场特征量的变化特点.证实了在随机介质模型中.各向异性系数的变化会引起波场记录上的某些统计特征的变化,归纳得出了若干结论.  相似文献   

6.
介质的弹性常数为三维四阶张量的分量,共有81个,由于应力张量和应变张量的对称性及能量密度是应变的二次函数,一般各向异常性介质的独立弹性常数可减为21个,如果介质具有较高的对称性,独立弹性常数的数目会更少。 对于地壳和上地幔,具有5个独立弹性常数的横向各向同性介质是一个非常好的近似,本研究中横向各向同性介质的对称轴方向可以是任意的(即对称轴可以不平等于铅直方向),在此情况下,需要进行坐标变换,如果已知介质在某一坐标系(其坐标轴平行或垂直于介质的对称轴)中的弹性常数,我们能够容易地利用变换公式得到变换后新坐标系中的弹性常数。 本文提出了一种方案,利用伪谱法既能模拟横向各向同性介质中的平面波,也能模拟点源激发的波场。在勘探地球物理和地震学中,模拟横向各向同性介拮中传播的平面波及区域源产生的波是最重要的研究课题之一。然而在一般各向异性介质中,很难或不可能确定弹性波的相速度和偏振方向,但在横向各向同性介质中,则可以通过坐标变换来实现,这里我们所提出的方法可以用于横向各向同性介质中弹性波的模拟。  相似文献   

7.
长波长假设条件下,各向同性背景地层中发育一组平行排列的垂直裂缝可等效为具有水平对称轴的横向各向同性(HTI)介质.基于不同观测方位的岩石地震响应特征变化,宽方位地震数据不仅可实现裂缝岩石弹性参数与各向异性参数的预测,同时也蕴含着丰富的孔隙度等储层物性参数信息.本文结合实际地震资料提出了贝叶斯框架下岩石物理驱动的储层裂缝参数与物性参数概率地震联合反演方法,首先基于AVAZ反演裂缝岩石的弹性参数与各向异性参数,并在此基础上通过统计岩石物理模型表征孔隙度、裂缝密度等各向异性介质储层参数与裂缝岩石参数的相互关联,并采用马尔科夫链蒙特卡洛(MCMC)抽样方法进行大量样本的随机模拟,使用期望最大化(EM)算法估计后验条件概率分布,最终寻找最大后验条件概率对应的孔隙度、裂缝密度等HTI裂缝介质储层参数即为反演结果.测井及实际地震数据处理表明,该方法能够稳定合理地从方位地震资料中获取裂缝岩石弹性参数与各向异性参数,并提供了一种较为可靠的孔隙度、裂缝密度等裂缝介质储层参数概率地震反演方法.  相似文献   

8.
储层中的裂隙会导致介质表现出各向异性.传统等效介质方法假定储层可以被等效成为横向各向同性介质,不适用于复杂裂缝储层的弹性等效建模和各向异性特征分析.本文在前人研究基础上,利用最小二乘方法开展3D复杂裂隙弹性等效数值模拟研究,并与常规等效介质理论方法进行比较,验证数值方法的精确性.数值方法不仅能有效避免Hudson理论和Eshelby-Cheng理论等解析方法在大裂隙密度条件下的参数估计误差,而且兼顾裂隙扰动作用对等效介质参数的影响.数值研究表明,数值方法能够较好地描述复杂裂缝介质的各向异性特征,为复杂裂缝储层地震波各向异性参数反演提供依据.  相似文献   

9.
裂缝预测是页岩气甜点预测的重要内容.页岩气储层高角度裂缝发育,具有明显的HTI介质各向异性特征,基于各向异性参数与裂缝密度的相关性,提出了一种裂缝密度反演的新方法.首先推导了一种能反映裂缝密度的方位AVO近似公式及其方位弹性阻抗方程(EVAAz方程);然后应用贝叶斯理论,形成了组稀疏方位弹性阻抗反演目标泛函,在此基础上...  相似文献   

10.
秦海旭  吴国忱 《地震学报》2014,36(6):1062-1074
裂缝密度是评价裂缝性储层的重要参数. 本文提出一种利用各向异性梯度计算裂缝密度的新方法,即运用裂缝介质等效理论将高陡倾角裂缝介质等效为横向各向同性介质,得到裂缝介质AVO响应特征,从而得到裂缝介质的各向异性梯度,之后由该梯度与裂缝密度的关系得到裂缝介质的裂缝密度. 实际处理中可以根据地震数据直接求取裂缝介质的各向异性梯度,最后求取裂缝介质的裂缝密度. 模型试算表明, 该方法能得到裂缝介质准确的裂缝密度. 反演结果与初始模型基本一致验证了该方法的正确性; 另外地震记录添加一定的信噪比后也能反演得到准确的裂缝密度,证明了该方法的稳定性.   相似文献   

11.
12.
洪泽地区由于沉积的特点,储层横向变化快,油藏受构造、岩性、油源多因素控制。在对该区三维AVO属性体解释中,利用多元回归方法求取了横波曲线,分岩性和含油气性统计了纵、横波、泊松比参数分布规律,建立了本区的含油砂岩的正演模型,从而降低了AVO解释的多解性。通过井-震结合对四种AVO属性数据体进行了标定,并确定了各属性体应用范围,进而进行了储层和含油气检测。实践表明,该方法能有效地利用AVO属性数据体进行储层预测及油气检测,具有一定的推广价值。  相似文献   

13.
This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, temporal-and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial resuits for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.  相似文献   

14.
We analyze the characteristics of different floating datums for static corrections and discuss the methods for determining them. The effect of different floating datum corrections was studied using theoretical model experiments, resulting in the conclusion that the velocity obtained after the floating datum correction with the minimum static correction errors depends on the velocity of the layer below the low velocity layer (LVL) lower boundary and is not related to topographic relief and LVL structure. For the real data processing case, wave equation numerical model experiments were conducted which resulted in a new method for calculating objective functions based on the waveform and modifications to the calculation equation for minimum static correction errors to make the method suitable for real data static correction processing using inhomogeneous velocity models with lower velocity boundary relief. Real data processing results demonstrate the method's superiority.  相似文献   

15.
Tilt and strainmeter data collected in the Friuli seismic area between 1977 and 1985 have been analyzed to study the stress-strain field and/or the elastic parameter variations in the area after the strong 1976 seismic event.We concentrated on three different frequency bands, respectively linked with the tidal waves, the seasonal effect and the secular one. The results from the first two bands point at an increase of the local rigidity modulus.Relevant stresses induced by the seasonal groundwater level variations have been found, thus supporting the hypothesis that they can act as triggering effects. The secular trend reveals significant tectonic stress modifications. It appears that the post-seismic tectonic stresses faded during 1981–1983, and were followed by a new loading phase.On leave from the Institute of Geodesy and Geophysics of the Chinese Academy of Sciences. Ph. D. student at the Institute of Geodesy and Geophysics, Trieste University.  相似文献   

16.
The working area is located in the industrially developed region of Rongshengpu-Qianjin, where a surface water system is developed, surface-layer lithology is complicated, and various kinds of hydrocarbon traps are buried at depth. The seismic data acquired previously couldn’t be interpreted due to the complex surface and geological conditions. Taking secondary 3D seismic from the Rongshengpu-Qianjin area as an example, this paper describes a set of techniques designed to overcome these difficulties and improve the quality of seismic data. The applied techniques included flexible acquisition geometry, low-noise receiver conditions, quantitative quality control, and so on. First Author Tian Qingbo, Engineer, graduated in Applied Geophysics Department from Chengdu Geology College in 1991. He engaged in acquisition and processing of seismic data, research of acquisition method for a long time. Now, he works in the department of Project and Technique in Geophysical Exploration Company of Liaohe Petroleum Exploration Bureau.  相似文献   

17.
非常规油气藏(如致密性地层及蕴藏油气的页岩地层)的重要特征是低孔、低渗,但裂隙或裂缝比较发育.为满足非常规勘探的需求,本文将孔、裂隙介质弹性波传播理论应用于多极子声波测井的井孔声场模拟,重点研究了致密介质中裂隙发育时多极子声波的传播机理以及衰减特征.井孔声场的数值计算结果表明裂隙的存在明显改变了弹性波和井孔模式波的频散、衰减和激发强度,尤其是井壁临界折射纵波的激发谱的峰值随着频率的增加逐渐降低,这与应用经典的Biot理论下的计算结果相反,且裂隙的存在也使得饱含水和饱含气时临界折射纵波激发强度的差异变大.井孔模式波的衰减与地层横波衰减和井壁流体交换有关,井壁开孔边界下致密地层裂隙发育还使得井孔斯通利波和艾里相附近的弯曲波对孔隙流体的敏感性增强,在井壁闭孔边界条件下引起井孔模式波衰减的主要因素是裂隙引起的地层横波衰减造成的,且在截止频率附近弯曲波的衰减与地层的横波衰减一致.数值计算结果为解释非常规油气地层的声学响应特征提供了参考.  相似文献   

18.
李向阳  Jianxin  Yuan 《应用地球物理》2005,2(3):153-167,i0001,F0003
在具有垂直对称轴横向各向同性介质中,利用四种参数来确定中间至远偏移距转换波(C-波)动校正。它们是C-波叠加速度Vc2,垂直速度比和有效速度比γ0和γeff以及各向异性参数χeff。我们将这四种参数作为C波叠加速度模型。C-波速度分析的目的就是确定这种叠加速度模型。C-波叠加速度模型Vc2,γ0,γeff,和χeff可以由P-波和C-波反射动校正资料获得。然而错误的传播是C-波反射动校正反演中的严重问题。当前短排列叠加速度由于是从双曲线动校正推算而得,因而其精度不足以为各向异性参数提供有意义的反演值。中间偏移非双曲线动校正不再被人们所勿略,而是可以用一个背景γ加以量化。非双曲线分析通过中间偏移距的γ校正量可以产生Vc2,若数据不含燥音,其误差小于1%。方法稳健,允许γ启始假定值的误差达20%。该方法也适用垂直非均匀各向异性介质。精度的提高使能够用4分量地震资料计算各向异性参数。为此提出了两种工作流程:双扫描和单扫描流程。理论数据和实际数据的应用表明这两种流程得出的结果其精度相似,但是单扫描流程比双扫描更有效。  相似文献   

19.
In transversely isotropic media with a vertical symmetry axis (VTI), the converted-wave (C-wave) moveout over intermediate-to-far offsets is determined by four parameters. These are the C-wave stacking velocity V C2, the vertical and effective velocity ratios γ 0and γ eff, and the anisotropic parameter X eff. We refer to the four parameters as the C-wave stacking velocity model. The purpose of C-wave velocity analysis is to determine this stacking velocity model. The C-wave stacking velocity model V C2, γ 0, γ geff, and X eff can be determined from P- and C-wave reflection moveout data. However, error propagation is a severe problem in C-wave reflection-moveout inversion. The current short-spread stacking velocity as deduced from hyperbolic moveout does not provide sufficient accuracy to yield meaningful inverted values for the anisotropic parameters. The non-hyperbolic moveout over intermediate-offsets (x/z from 1.0 to 1.5) is no longer negligible and can be quantified using a background γ. Non-hyperbolic analysis with a γ correction over the intermediate offsets can yield V C2 with errors less than 1% for noise free data. The procedure is very robust, allowing initial guesses of γ with up to 20% errors. It is also applicable for vertically inhomogeneous anisotropic media. This improved accuracy makes it possible to estimate anisotropic parameters using 4C seismic data. Two practical work flows are presented for this purpose: the double-scanning flow and the single-scanning flow. Applications to synthetic and real data show that the two flows yield results with similar accuracy but the single-scanning flow is more efficient than the double-scanning flow. This work is funded by the Edinburgh Anisotropy Project of the British Geological Survey. First Author Li Xiangyang, he is currently a professorial research seismologist (Grade 6) and technical director of the Edinburgh Anisotropy Project in the British Geological Survey. He also holds a honorary professorship multicomponent seismology at the School of Geosciences, University of Edinburgh. He received his BSc(1982) in Geophysics from Changchun Geological Institute, China, an MSc (1984) in applied geophysics from East China Petroleum Institute (now known as the China University of Petroleum), and a PhD (1992) in seismology from the University of Edinburgh. During 1984–1987, he worked as a lecturer with the East China Petroleum Institute. Since 1991, he has been employed by the British Geological Survey. His research interests include seismic anisotropy and multicomponent seismology.  相似文献   

20.
岩石孔隙结构是控制砂岩和碳酸盐岩的地震波速度和渗透率的重要参数之一。如果两种类似的岩石其孔隙度一定,而渗透率不同,那么它们的声波速度相差2km/s,而渗透率两者则可能相差近6个数量级,即从0.01mD到20mO。在本文中我们总结了由一个广义孔隙弹性理论简化的一个双参数弹性速度模型,以描述孔隙结构对弹性波的影响。由于矿物和储层的流体是给定的,我们用孔隙度和骨架的柔性(挠性)因子来确定速度模型,这个模型可以用于地震反演和储层表征,已改善孔隙度和储量的计算骨架柔性因子可用于岩石结构(PST)类型的定量分类,并可以利用叠前、叠后的两种地震资料将其与孔隙的连通性和渗透率联系起来。本项研究同时也有助于说明为什么振幅与偏移距(AVO)分析用于流体检测在某些情况下失败。这是由于孔隙结构对地震波的影响能够掩盖所有流体效应,特别是在碳酸岩中。  相似文献   

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

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