首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 343 毫秒
1.
基于声学全波形反演的油气藏地震成像方法   总被引:4,自引:4,他引:0       下载免费PDF全文
岩性油气藏在我国天然气勘探开发中占有非常重要的位置,其分布区域的成像是合理布设井位,提高钻井成功率的关键之一.本文首先基于地下介质的声学近似和波场回传理论,利用频率域单程声波方程延拓计算地震波场,进行全波形反演,获得地层密度和体积模量的定量成像,并依据油气藏物性特征和流体饱和多孔介质岩石物理模型,简要讨论了孔隙度和饱和度与密度及体积模量的关系,明确了地震油气藏成像新概念.在此基础上,定义了基于流体体积模量和孔隙度的成像函数,进行油气藏成像.理论模型计算表明该方法是可行的.通过对西部地区某气田二维地震数据处理,实现了致密砂岩气藏成像,钻井结果证实了气藏区域成像位置的准确性和方法的有效性.  相似文献   

2.
基于贝叶斯理论的叠前多波联合反演弹性模量方法   总被引:8,自引:6,他引:2       下载免费PDF全文
AVO反演可以获得地层岩性和流体信息,而叠前反演问题都是高维的和非适定的,因此获得可靠稳定的解对叠前反演至关重要. 本文给出了一种基于贝叶斯理论的纵波和转换波联合反演密度比和模量比的方法. 鉴于剪切模量比、体积模量比可以更好地指示油气,基于岩石物理中速度比与模量比之间的关系,将此关系式代入Zoeppritz方程的近似形式Aki-Richards公式中,得到与模量比有关的反射系数近似公式. 联合纵波和转换波,利用最小二乘准则构建目标函数,最终反演出密度比、剪切模量比、体积模量比三个参数. 在反演过程中引入贝叶斯理论,假定先验信息服从高斯分布,待求参数服从改进的Cauchy分布,并去除待求参数之间的相关性. 利用模型数据和实际数据对本文方法进行测试,并与常规的单独利用纵波数据来反演方法进行比较,结果表明联合反演稳定性更好、精度更高、抗噪音能力更强,验证了本文方法的可行性和有效性.  相似文献   

3.
利用二维谱分析井径变化对声波测井响应的影响   总被引:1,自引:1,他引:0       下载免费PDF全文
本文简要介绍了井孔声场二维谱的计算方法,实现了对不同井径井孔中的声场二维谱和声场模拟计算,得到了不同井径井孔声场二维谱的分布图以及声压波形图.通过对不同井径井孔中声场二维谱的特征分析,讨论了不同井孔中的声波场中的纵波、横波以及各种模式波的特征.通过对比二维谱可以明显地看出纵波、横波以及各种模式波均受井径变化的影响.实践证明井孔声场二维谱能够直观地显示出不同井径中声场传播的特征.  相似文献   

4.
论文使用高分辨有限体积法模拟二维声波在多层介质中的传播现象.该方法需要求解双曲型偏微分方程的黎曼解,数值解在强速度变化的分界面上满足数值通量守恒,可以有效地计算数值解的不连续性.根据数值实验给出的波场快照和接收点上的波形变化,可以清晰地观察到声波在介质分界面上的多次反射波,说明该方法不但可以得到高分辨率的波场快照,而且...  相似文献   

5.
随着城市化的发展,城市地球物理日益成为地球物理研究的重要方向,地震成像是构建城市地下空间三维/四维图像的重要手段,但面临观测成本高的困难.近年来国际上新发展的分布式光纤声波传感器作为高密度地震观测系统已经在地震层析成像方面得到了应用,在提高成像分辨率的同时,又降低了观测成本.本研究使用国产分布式光纤声波传感器开展了观测实验,利用480m埋地光缆记录了13h背景噪声,计算得到噪声互相关函数,获得了高频Rayleigh面波信号.采用多道面波分析方法提取相速度频散曲线,其结果与传统检波器记录和主动源结果较为一致.采用遗传算法反演得到了研究区内二维S波速度剖面,获得了下方沉积物横向变化特征.通过本次实验,初步验证了国产设备开展地震背景噪声成像研究、构建地下浅层结构模型的可行性.  相似文献   

6.
地球介质的非线性弹性造成了岩石动态和静态弹性模量在数值上的差异。本文用超声波速度测量得到岩石的动模量,用含有小循环的加卸载实验得到小循环模量,用大循环的应力、应变曲线斜率得到静模量,用测点的应力、应变之比得到杨氏模量。通过改变小循环的应变振幅来研究模量与应变振幅的关系,进而比较用不同测量方法得到的模量值的差异。结果表明,动模量值最大,小循环模量值次之,大循环静模量值再次之,而杨氏模量最小。岩石非线性模量既是所受应变水平的函数,也是应变振幅的函数,模量随应变振幅的增加以指数形式下降,随应变水平增加以指数形式上升。最后简要介绍了P-M模型的基本概念,用岩石中微裂缝具有不同的分布密度,不同的打开和关闭压力来解释模量与应变振幅的关系。  相似文献   

7.
叠前逆时偏移是目前精度较高的成像方法,然而严重的低频噪声降低了逆时偏移的构造成像精度,偏移噪声的压制是逆时偏移必需要考虑的问题。分析了低频噪声的产生机理,并根据声波方程Poynting矢量的方向指示地震波场的传播方向的原理,分离出上、下、左、右行波,该方法计算量和存储量都远远小于常用的二维傅里叶变换分离方法,进而提出归一化的波场分离互相关成像条件,以压制逆时偏移低频噪声,提高成像精度。实现了Marmousi模型的试算,表明在波场延拓过程中利用Poynting矢量能够较好的分离上、下、左、右行波,与常规方法、拉普拉斯滤波、二维傅里叶变换波场分离的成像结果对比表明,成像时使用归一化的波场分离互相关成像条件能更好的压制偏移噪声,得到精度更高的逆时偏移成像结果。  相似文献   

8.
利用深海地层的声波测井数据进行远探测成像是海洋油气勘探的一个新的应用课题。声波远探测可以在一维井孔中对井周数十米范围内的复杂构造进行二维成像。本文基于偶极声源的指向发声特性和四分量数据,实现反射体的方位探测功能,从而对井周进行三维成像。我们首先对野外四分量数据进行矩阵旋转得到多方位数据;此基础上,对旋转后偶极数据采用一系列处理步骤得到反射体成像。根据以上横波成像原理,用中国南海深海斜井采集的正交偶极四分量数据对井周围50 m范围内构造进行反射横波成像,刻画了反射体的形态并确定反射体的走向。井旁地层和断层构造的横波成像结果与该井所在的正交地震剖面的产状相符,断层的产状也与区域应力和斜井应力分析一致。结合其他测井资料(井壁电成像),进一步确定了远探测成像中地层的倾向。本文的研究结果表明:将井中的声波远探测成像与地震成像相结合,可以精细刻画井旁构造的空间形态,并填补这两种测量尺度的空白。  相似文献   

9.
任意广角波动方程叠前逆时深度偏移   总被引:8,自引:0,他引:8       下载免费PDF全文
何兵寿  张会星  张晶 《地震学报》2008,30(5):491-499
从任意广角声波方程出发,在深度与时间方向采用二阶差分代替微分,在水平方向上采用高阶差分代替微分,得到任意广角声波方程逆时延拓的有限差分格式. 通过差分求解程函方程得到叠前逆时深度偏移的成像条件,在此基础上实现了二维各向同性介质中任意广角方程的叠前逆时深度偏移. 对理论模型的偏移处理表明,该方法能够有效地消除波场下行传播过程中的层间反射,减小偏移噪声. 同时,该方法还具有保真性高、 适用于大倾角地层的偏移等优点.   相似文献   

10.
本文将二维时间空间域和频率空间域声波全波形速度反演方法分别应用到Marmousi模型,进行数值试验.两种方法均采用相同的观测系统和其他的参数,理论模型的数值试验结果证实了:使用较多的计算集群的CPU进行二维频率空间域直接法声波全波形反演时,其加速有限(正演数值模拟的计算量主要用于稀疏矩阵的LU分解,炮点计算波场时为线性关系).二维时间空间域声波全波形反演计算时更灵活,多炮同时计算时,可以多倍提高其计算效率;二维声波全波形速度反演时,直接法求解频率空间域的计算速度远快于时间空间域,所需要的计算机内存也比时间空间域少.二维声波全波形速度反演时,相比较于时间空间域的方法,频率空间域直接法声波全波形反演具有计算速度快和节省计算机内存需求的优势.  相似文献   

11.
A modern approach to migration is to perform wavefield extrapolation, subject to an imaging condition. Correct wavefield extrapolation requires that the boundary conditions at the array of geophones satisfy the wave equation. A sufficient condition is to perform the survey with a single stationary source. Contrary to this condition, many VSPs are conducted in deviated wells, where the source is maintained vertically above the down-hole geophone at each well station. Such a survey fails to provide the boundary conditions theoretically necessary for wave-equation migration. A recently published inversion scheme, referred to as acoustic generalized Radon transform migration (GRT migration), was developed to handle any configuration of sources and geophones, including moving-source deviated-well VSP surveys. GRT migration may be viewed as a weighted version of the generalized Kirchhoff migration, derived in this paper from the exploding-reflector model. When a VSP-survey geometry has been specified, GRT migration can be expressed in terms of array parameters, and compared with the equivalent expression for Kirchhoff (wave-equation) migration. The differences between the two integrals are significant and their effect is demonstrated on VSP data.  相似文献   

12.
波场延拓深度滤波方法   总被引:5,自引:1,他引:4       下载免费PDF全文
消除面波是地震数据处理中的一个重要内容.本文提出了基于15°波动方程的深度滤波方法.由于面波与有效反射波具有不同的传播深度,可利用波场向下延拓方法将二者进行波场分离.把向下延拓后的波场中集中在地表附近的面波能量切除后,再将波场重新延拓回原始的观测面,达到去除干扰的目的.实际资料处理显示:方法计算稳定,消除面波能力强,能更好地保持波场的有效成份和幅值,符合波场的实际传播状态,表明该方法正确可行.  相似文献   

13.
复杂地表的单程波动方程地震叠前正演   总被引:4,自引:0,他引:4  
作者基于数学检波器和等时叠加原理,实现了复杂地表的单程波动方程地震叠前正演模拟。该方法采用虚拟的数学检波器接收地下的反射地震信号,灵活地将接收点布置在地表的任何地方,从而满足地表起伏的要求。此外,根据等时叠加原理, 该方法采用单程波动方程进行波场延拓和成像,计算简单快速。通过复杂正断层的数值模拟,得到了高信噪比的共炮集地震记录,并采用适用于起伏地形的深度偏移方法对该共炮集地震记录进行了叠前深度偏移,较好地实现了地震波的偏移归位,从而证明了这里提出的起伏地表的单程波动方程地震叠前正演方法是正确和有效的。  相似文献   

14.
A transmission + reflection wave-equation traveltime and waveform inversion method is presented that inverts the seismic data for the anisotropic parameters in a vertical transverse isotropic medium. The simultaneous inversion of anisotropic parameters and ε is initially performed using transmission wave-equation traveltime inversion method. Transmission wave-equation traveltime only provides the low-intermediate wavenumbers for the shallow part of the anisotropic model; in contrast, reflection wave-equation traveltime estimates the anisotropic parameters in the deeper section of the model. By incorporating a layer-stripping method with reflection wave-equation traveltime, the ambiguity between the background-velocity model and the depths of reflectors can be greatly mitigated. In the final step, multi-scale full-waveform inversion is performed to recover the high-wavenumber component of the model.  We use a synthetic model to illustrate the local minima problem of full-waveform inversion and how transmission and reflection wave-equation traveltime can mitigate this problem. We demonstrate the efficacy of our new method using field data from the Gulf of Mexico.  相似文献   

15.
弹性波方程正演的粗粒度并行算法   总被引:3,自引:3,他引:0       下载免费PDF全文
波动方程正演在地震资料采集、处理、解释与反演中均能发挥重要作用,但现有的基于求解地震波动方程的正演算法由于受庞大计算量的制约而难于大规模应用于工业生产.本文从二维弹性波动方程出发,研究了利用有限差分法并行求解该方程的基本思路与方法,给出了适于并行求解的计算空间划分方法与通信方案,分析了不同参数条件下并行程序的运行时间、加速比与效率.引入消息传递接口(MPI)实现了弹性波方程的并行求解,极大地提高了数值求解弹性波方程的计算效率.  相似文献   

16.
针对处理起伏地表(或含地下不规则波阻抗界面)条件下发展起来的地震波场数值模拟算法的模拟结果与解析解(大多数情形下无法得到)无法进行对比,且其有效性和正确性难以验证的情况,本文提出了一种可以相互验证波场数值模拟结果与射线追踪数值模拟结果的正确性和有效性的佐证方法,验证了参考射线追踪法.其中,波场数值模拟中采用曲线网格DRP/opt MacCormack有限差分法,射线追踪模拟则采用分区多步三角网格最短路径算法.通过系统对比上述两种方法得到的波场快照、单炮地震记录,以及合成理论地震图的结果显示,本方法相互作证了两种方法所得结果的正确性和有效性.双层和三层起伏层状模型的对比分析结果表明,这种方法不但可以加深理解地震波在复杂介质中的传播规律,同时射线法的引入为清晰识别和标定地震波场数值模拟中各种不同震相提供了一种便捷的途径.   相似文献   

17.
基于波动方程有限差分算法的接收函数正演与偏移   总被引:20,自引:6,他引:14       下载免费PDF全文
针对接收函数正演与偏移, 本文采用波动方程有限差分算法. 借鉴成熟的勘探地震学方法, 引入等效速度概念, 建立接收函数转换波与地震勘探反射波的等效走时方程, 实现了基于波动方程有限差分算法的接收函数正演与偏移. 数值计算表明, 波动方程有限差分叠后偏移方法可以对点绕射和穹隆构造模型实现高精度成像. 本文利用数值计算讨论了波动方程有限差分叠后偏移与Kirchhoff叠后偏移对于接收函数偏移的适用性, 还对偏移过程中速度模型的误差进行了分析.  相似文献   

18.
A new source location method using wave-equation based traveltime inversion is proposed to locate microseismic events accurately. With a sourceindependent strategy, microseismic events can be located independently regardless of the accuracy of the source signature and the origin time. The traveltime-residuals-based misfit function has robust performance when the velocity model is inaccurate. The new Fréchet derivatives of the misfit function with respect to source location are derived directly based on the acoustic wave equation, accounting for the influence of geometrical perturbation and spatial velocity variation. Unlike the mostly used traveltime inversion methods, no traveltime picking or ray tracing is needed.Additionally, the improved scattering-integral method is applied to reduce the computational cost. Numerical tests show the validity of the proposed method.  相似文献   

19.
Imaging in geological challenging environments has led to new developments, including the idea of generating reflection responses by means of interferometric redatuming at a given target datum in the subsurface, when the target datum lies beneath a complex overburden. One way to perform this redatuming is via conventional model-based wave-equation techniques. But those techniques can be computationally expensive for large-scale seismic problems since the number of wave-equation solves is equal to two times the number of sources involved during seismic data acquisition. Also conventional shot-profile techniques require lots of memory to save full subsurface extended image volumes. Therefore, we can only form subsurface image volumes in either horizontal or vertical directions. To exploit the information hidden in full subsurface extended image volumes, we now present a randomized singular value decomposition-based approach built upon the matrix probing scheme, which takes advantage of the algebraic structure of the extended imaging system. This low-rank representation enables us to overcome both the computational cost associated with the number of wave-equation solutions and memory usage due to explicit storage of full subsurface extended image volumes employed by conventional migration methods. Experimental results on complex geological models demonstrate the efficacy of the proposed methodology and allow practical reflection-based extended imaging for large-scale five-dimensional seismic data.  相似文献   

20.
裂缝模型声波衰减系数的数值模拟   总被引:3,自引:0,他引:3       下载免费PDF全文
超声波数据是研究碳酸盐岩储层孔隙结构的重要资料,超声波响应可用来反演孔隙结构.本文从声波波动理论出发,以裂缝模型为研究对象,开展超声波数值模拟研究.结果表明:随着裂缝孔隙度的增加,岩石声波衰减系数总体呈上升趋势,但孔隙度增长方式不同,衰减系数的递增函数规律亦不同;随着裂缝延伸方向与声波激发方向的夹角增大,岩石声波的衰减系数线性增加;衰减系数随着裂缝密度增大呈幂函数规律递增;利用数值模拟的方法能避免人为误差、节约实验成本,为裂缝结构的岩石声波实验研究拓展了思路.  相似文献   

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

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