45°方程叠后深度偏移的方向导数法

时间偏移由于忽略了速度的横向变化而导致成像偏离正确位置。作者从45°深度偏移方程出发，引出方向导数，并将速度横向变化所引起的折射效应项纳入此方向导数内，使所构造的深度偏移方法在形式上与相应的时间偏移方法一致。因此，只需对目前生产中应用的时间偏移模块略作修改，即增加非整数节点插值部分，就可实现深度偏移。这样不仅节省了计算工作量，提高了计算效率，而且也缩短了编程的时间。数值试验结果表明，方法是可靠有效的。

双复杂条件下非倾斜叠加精确束偏移方法及应用Ⅰ——声波方程

近年来,油气勘探的重心正转向具有复杂地表和复杂地质体的双复杂区域.本文发展了一种精确的双复杂条件下基于地表倾角信息的非倾斜叠加束偏移方法,相对于传统束成像方法无需进行三方面处理:(1)高程静校正;(2)相位校正;(3)束中心与接收点之间关于速度和束出射角的近似替换,因而具有更高的成像精度.通过加拿大逆掩断层模型、中原油田断层模型及实际资料的偏移试算,并与传统束偏移及波动方程偏移成像结果对比可知:本文非近似束偏移方法在近地表、高陡倾等构造处的成像精度、反射界面成像振幅等方面优于传统的偏移方法,以此验证了本文非倾斜叠加精确束偏移方法的正确性、优越性及适应性.

VSP的双程无反射波动方程逆时偏移

VSP 资料通常为共炮点道集,因所含数据量较少,采用常规偏移方法,难以获得良好的偏移结果。而用双程无反射波动方程作逆时偏移,可以剔除内部界面反射波和层间多次波,同时使用 Clayton 的吸收边界条件,还可以消除VSP 中的边界反射波。可见采用这种逆时偏移方法能够较好地改善偏移效果。其成像条件仍由 Claerbout 成像原理推广,求解镜像方程得到,此法是一种适用性很强的波动方程偏移方法,既适用叠前 VSP 记录,又适用常规多偏移距的炮点道集,还可用于 VSP 和地面资料的联合偏移,并能适应各种速度的变化。此法用于实际 VSP 资料处理时,要预先切除初至波,剪接边道或设置边道阻尼,并要作空间插值补道、频率滤波等项处理,否则将影响 VSP 偏移剖面的质量。

波动方程深度偏移的频率相关变步长延拓方法

发展了波动方程深度延拓的频率相关变步长深度延拓方法和表驱动的单点波场插值技术.前者通过减少深度延拓的次数减少了波动方程深度偏移的计算量,而后者用很少的计算量实现了等间距、理想采样的深度成像.就同一偏移方法,采用频率相关变步长深度延拓加单点插值,其计算量大约是常规的等间距采样延拓方法的三分之一,但两者的成像效果基本相同.文中以最优分裂Fourier方法为例,用二维理论数据(Marmousi模型)和三维实际地震资料验证了这一方法,但这一方法可适用于各类频率域波动方程深度偏移方法.

波动方程叠前深度偏移角度道集求取方法

针对常规直接转换法和直射线叠前时间偏移(PSTM)法提取的角度道集无法满足地质构造复杂地区的AVA分析与反演精度的问题,研究基于共炮检距道集波动方程保幅叠前深度偏移方法,提取角度域共成像点道集。首先,对地震数据进行叠前深度偏移得到炮检距域共成像点道集;其次,在频率—波数域对炮检距域共成像点道集利用快速插值映射法抽取角度域共成像点道集。通过与直接转换法和直射线PSTM法提取的角度道集比较,分析不同方法求得的角度道集信息的准确性和AVA特征。数值试算和实际资料处理结果表明:共炮检距道集波动方程保幅叠前深度偏移方法求取的角度道集的振幅相对保真,角度范围更广,更适用于AVA分析。

波动方程地震偏移成像方法的现状与进展

综述了波动方程地震偏移成像方法研究的现状 ,指出了各种算法的问题和解决方案 ;给出了相空间小波分析步进算法偏移成像以及保幅保结构李群算法偏移成像的最新进展 ;对波动方程偏移成像发展趋势进行展望 .

波动方程偏移速度建模:一种直接反演方法

基于炮域波动方程叠前深度偏移所产出的角道集,发展了不需要多次迭代的波动方程偏移速度建模方法.文中分析了炮域波动方程偏移生成角道集的方法,给出了非均匀介质中角道集同相轴曲率与速度误差的定量关系,发展了基于同相轴曲率的速度模型直接更新算法.在初始模型较合理的情况下,应用一次炮域波动方程偏移计算,即可得到较准确的速度模型.这一方法可用于解决由走时层析成像等方法得到的速度模型与波动偏移方法不匹配的问题,避免了直接由波动方程偏移进行速度建模所需要的多次偏移计算.数值结果证明了本文方法的有效性.

单程波方程偏移算法的相位问题研究

利用稳相原理对叠后单程波方程偏移、二维和三维叠前偏移算法中的相位进行了分析和研究,从理论上论证了叠后偏移算法可以保持输入子波的相位特征,而叠前偏移算法会改变输入子波的相位特征,叠前单炮深度偏移结果相位与输入数据子波相位之间相差了一个因子。上述理论通过基本单位脉冲响应和一个水平反射界面模型的数值模拟得到了验证。从而从理论上解决了波动方程偏移与Kirchhoff积分偏移结果相位不一致问题,这对于正确标定反射层的振幅和深度具有实用意义。

基于GPU计算平台的三维波动方程叠前深度偏移

GPU/CPU协同并行计算是一项全新的、具有革命性的技术,将对地震资料处理行业产生深远的影响。本文在充分发挥GPU存储带宽宽、多寄存器和多处理器结构特点的基础上,通过改进地震成像并行计算模块的核心函数,形成一套波动方程地震成像GPU高效并行算法,并取得理想的加速比和可扩展性。所形成的一套GPU/CPU平台下波动方程叠前深度偏移软件,已投入地震资料处理实际应用。理论模型试算和实际资料的处理验证了本方法的有效性和实用性。

VSP数据波动方程叠前深度偏移成像及立体地震成像

,徐蔚亚 ,王建民 ,马在田 .VSP数据波动方程叠前深度偏移成像及立体地震成像 .石油地球物理勘探 ,2 0 0 1,36 (5 ) :5 17～ 5 2 5 本文从三维声波方程出发 ,通过引入参考速度 ,推导出折射项和绕射项方程 ,并对绕射项进行优化展开 ,得到适应任意变速情况下的 VSP数据波动方程叠前深度偏移方程及其差分格式 ;讨论了有限差分法的误差补偿、波场成像条件 ;基于波场传播的线性叠加原理 ,在一定成像条件下 ,提出了把地面地震记录、VSP记录及井间地震记录偏移到同一成像空间的完整的地震成像方法——立体地震成像。立体地震成像的实质是 :地面地震数据、VSP数据及井中反射地震数据的成像过程可以统一成一个波动方程叠前深度偏移过程 ,同时 ,每种数据的成像结果也可以单独输出。模型试算结果表明 ,在立体地震成像过程中 ,每一种观测方式对目标地质体的照明范围不同 ,它们对目标地质体的成像互为补充与加强 。

基于波动方程偏移的地震数据地层成像

地层成像是对地下地层物性空间变化的近似反演,对于岩性油气藏和非常规油气藏的勘探开发具有重要价值。文中根据波动方程偏移的思想提出了一种有关反射体波阻抗相对扰动的地层成像方法。基于地震波运动学特征准确的地下介质光滑模型,应用地震波传播理论推导了基于反射体波阻抗相对扰动的地震数据线性正演表达式;然后结合线性反演理论,利用地震数据近似反演波阻抗相对扰动,构建了基于波动方程偏移的以波阻抗变化为主要目标的地震数据地层成像方法。在地层成像中,如果对波场做角度分解,则可得到地层成像的角度域共成像点道集,如果对波场不做角度分解,则可快捷地得到角度域平均的地层成像结果。对于声波方程描述的地震数据,利用该地层成像方法可获得地层的声波阻抗相对扰动成像;对于标量波方程描述的地震数据,利用该成像方法可获得地层的速度相对扰动成像。该地层成像方法与地震数据逆时偏移方法的计算过程相同、计算量相当,合成地震数据试验取得了理想的地层成像效果。

双复杂条件下的波动方程叠前深度偏移

复杂地表条件下的地震勘探越来越被人们所关注.双复杂条件下的叠前深度偏移方法是解决复杂地表条件和复杂地质构造成像的有效手段."波场上延"法能实现由非水平观测界面开始的偏移过程,解决复杂地表对地下构造成像的影响.复杂理论模型的试算以及实际资料处理表明,"波场上延"方法较好地克服了起伏地形对地下构造成像的影响,取得了令人满意的效果,实现了波动方程基准面校正和深度成像的有机结合.

大别造山带宽角度反射射线方程偏移研究

二维介质中的宽角度反射射线方程偏移是一种基尔霍夫型叠前深度偏移算法．基于射线理论，这一算法对二维不均匀介质具有广泛的适用性，在宽角度反射研究中能够显示其优越性．通过计算射线场，插值出震源和接收器的波场，对旅行时间ｔＰ＝ｔＳＰ＋ｔＰＲ的深度点Ｐ，可以在地震图中找到与该点对应的强度（或振幅）ＡＰ．将ＡＰ按照一定的关系进行转化，移至该深度点，对多个震源和接收点进行同样的计算，并将结果进行叠加，就完成了射线方程偏移．数值模拟计算取得了令人满意的结果；处理了通过大别造山带的庄墓—岳西—黄梅—张公渡人工地震测深剖面，获得了沿测线的莫霍界面构造形态．结果表明：南大别岩块高压变质岩带与这里莫霍界面局部抬升有关．北大别岩块下面有“山根”存在，与重力均衡原理相一致．晓天—磨子潭断裂作为华北板块和扬子板块的碰撞缝合带穿透了莫霍界面．在六安附近有一深大断裂，它是洛南—明港断裂的东延部分，在这里得到了进一步的证实．

盐丘模型VSP数据的波动方程偏移

本文对一个盐丘侧翼模型的VSP数据成像进行了理论试算和分析,试算的结果表明了VSP数据波动方程偏移方法的可行性,展示了该方法的效果和潜力。

基于Xeon Phi平台的波动方程叠前深度偏移

波动方程叠前深度偏移适用于强横向变速介质,是一种高精度成像方法,但其巨大的计算量阻碍了该技术的应用。Xeon Phi是一种全新的高性能计算设备,为波动方程叠前深度偏移方法的推广应用提供了新的技术支持。以裂步傅里叶算子为例,介绍了面向Xeon Phi平台的偏移算法移植和优化方法,即采用offload模式将计算核函数加载到Xeon Phi设备上,在Xeon Phi协处理器上采用多线程方式,并且调整程序结构,充分利用SIMD矢量引擎提高向量化处理效率。扩展负载动态均衡的并行框架,形成了一套适用于大规模异构系统、基于Xeon Phi平台的波动方程叠前深度偏移软件。实际数据测试表明Xeon Phi平台可以极大地提高地震偏移处理效率,具有良好的可扩展性。

波动方程偏移技术在泌阳凹陷资料处理中的应用

利用美国CoreLab公司的叠前深度偏移VIEWS软件系统,基于PC集群并行运算系统,对泌阳凹陷南部陡坡带进行叠前深度偏移处理。系统介绍了针对该区的深度域速度建模的基本思想和叠前深度偏移的原理。在叠前深度偏移剖面上,各主要目的层成像清楚,尤其是边界大断面及盆地基底较清晰,位置可靠;层间接触关系—地质层位与断层的边界清楚,小断层断点干脆,便于地质人员分析解释。

VSP资料的叠前波动方程逆时偏移

从波动方程出发,利用有限差分法导出波场外推公式,由此可以进行波场正演计算和偏移。用 P 波方程编制的偏移程序 VSP MIG 既可用于地面资料的偏移,也可用于井中资料的偏移;既可用于 P 波偏移,也可用于 P-SV 波偏移,但需相应的变化速度模型。对于 P 波和 P-SV 波,偏移的成像时间也是完全一样的。分别偏移分离后的上行 P 波和上行 P-SV 波,可以得到纵波和横波两张偏移剖面,这样比纵、横波资料的联合一次偏移可以得到更多的信息。

波动方程叠前深度偏移的GPU技术

偏移成像是地震资料处理流程中重要的一环,目前工业界所采用的方法是克希霍夫时间偏移,该方法的优点在于计算效率高,模型易于获取。但因其基于高频近似,难以适应多次到达,欠照明影响及横向变速剧烈的地质情况。波动方程叠前深度偏移则能较好地解决这些问题,但其计算成本远高于克希霍夫偏移方法,这也成为其在实际生产应用中的瓶颈之一。这里将利用GPU技术,研究波动方程叠前深度偏移的高性能计算方法,使之能够满足生产中海量数据处理的需要。

基于波动方程偏移的宽方位三维角道集提取

基于三维波动方程偏移的宽方位三维角度道集不仅包含了局部反射角度的信息,还包含反射平面在地下的方位角信息,它充分反映了宽方位采集观测系统下勘探目标不同方位的照明情况。这一特点在复杂构造的三维偏移速度分析中是非常重要的。笔者从角度分解的思想出发,讨论了射线参数、反射角以及局部方位角的关系,实现了宽方位情况下局部偏移距域共成像点道集向角度域共成像点道集的映射方法,并给出具体流程。模型测试和实际资料应用结果表明,用该方法提取的角度道集可以真正反映局部反射处的不同方位角的反射信息,可以为全三维勘探地震成像中的深度偏移速度分析研究提供道集依据,也可以用于三维的振幅随角度关系分析。

基于超大内存节点的波动方程逆时偏移高效实现

随着地震数据规模不断增大,进行波动方程逆时偏移计算时,采用的重复正传计算与波场数据硬盘缓存的波场数据存储策略是制约逆时偏移计算效率的一个主要因素。为此提出了利用混合内存技术构建高达1.5 TB的超大内存计算节点,以满足大数据量波场数据无损压缩后内存存储需求,实现最小计算量的全波场数据内存存储策略,解决大数据量波场数据存储问题,提升偏移计算效率。实际生产数据测试结果表明,对于大规模数据,节点配置超大内存可使高密度炮数据偏移计算效率提升46.1%以上,计算效率提升显著,可缩短大规模地震数据逆时偏移计算时间,降低计算能源消耗。该技术还可应用于其它大内存需求的地震数据处理,因此具有很好的实用性和经济效益。