共偏移距可分表示法叠前深度偏移

基于双平方根方程的共偏移距可分表示法叠前深度偏移用于复杂介质成像, 该方法在中点-偏移距坐标中同时向下延拓炮点和检波点波场, 实现方式采用正反傅立叶变换.构造的双平方根方程波场延拓算子能够使波数域变量与空间(速度) 域变量分离, 波数域内进行相移计算, 在空间域对因介质横向变速引起的时移作修正.地震数据偏移不需要逐炮计算, 具有较高的效率.对Marmousi模型数据的偏移成像结果显示, 该方法较好地成像强横向变速介质中的复杂构造.      

波动方程角度域共成像道集

针对基于波场深度递推的波动方程叠前深度偏移成像方法难以准确地给出如Kirchhoff积分叠前深度偏移成像方法的偏移距域共成像道集的不足,在分析研究现有的各种波动方程偏移成像共成像道集方法的基础上,利用波场外推的单平方根算子和波场的窗口Fourier框架展开与重构方法,提出一种局部角度域共成像道集方法——成像点处反射角共成像道集方法.把这种角度域共成像道集方法应用于国际标准的Marmousi模型数据和一条实际二维地震数据都取得了理想的结果.提出的波动方程反射角度域共成像道集方法可为进一步的叠前偏移数据振幅随角度变化分析和偏移速度分析工作提供基础.     

泌阳凹陷南部陡坡带地震资料的炮域波动方程叠前深度偏移

泌阳凹陷南部陡坡带是河南油田的重点勘探老区,其构造产状不仅陡,断层也多,叠后时间偏移难以得到较好的地震成像效果,基于Kirchhoff积分偏移方法具有运算速度快、效率高的特点,以及炮域波动方程方法能保持地震波的动力学特征,应用Kirchhoff积分偏移方法进行了偏移速度分析,建立了精细的速度模型.然后基于该模型在深度域进行了剩余速度分析以及层析成像,进一步提高了速度模型的精度. 通过分析偏移孔径、去假频参数以及延拓步长等成像处理参数对成像效果的影响,确定出最佳偏移处理参数,最后利用炮域波动方程方法对工区的三维地震数据进行了叠前偏移成像.成像结果表明,该方法能够使该区陡坡带反射层位得到较好的偏移,信噪比和横向分辨率都得到提高.     

复杂地表地震资料叠前深度偏移成像

分析了复杂地表地区复杂地震波场的形成原因,探讨了基于高角度时空域单程波动方程适用于起伏地表之下构造成像的逆时叠前深度偏移方法。从声波方程有限差分法合成起伏地表上的炮集记录知,当激发点位于低速层中,共炮点道集中低速层表面和高速层表面接收到的信号强度相差较大:低速层表面的信号很强,而高速层表面的信号很弱,低速层中产生了较强的槽波。将该逆时叠前深度偏移方法应用于起伏地表地震波场的偏移处理。虽然在偏移成像前既没有压制与低速层有关的槽波,又没有压制随机噪声,但偏移剖面上界面清晰、位置正确,断层面也得到了很好的聚焦。     

单程波法地震波衰减特征数值模拟分析

在非均匀粘弹性介质中,地震波衰减研究能够有效揭示地层吸收特性对地震波的影响。这里从单程波方程出发,结合粘弹性介质中的复波数理论,考虑介质的横向变化,采用分步傅立叶波场延拓技术,推导出适用于横向变速介质的带Q地震记录计算公式,最后将此方法应用于四个不同含Q地层模型。试验表明：单程波法计算速度快,波场信息简单清晰,易于深入分析和反演算法的建立。并且,Q值对地震波波形、能量、频率等都产生了较为明显的影响,对实际资料能量补偿、地震资料分辨率的提高,以及利用衰减特征进行储层识别有一定指导价值。     

基于保幅波动方程的广义高阶屏地震偏移成像方法

原始地震数据中饱含丰富的反映相位的走时信息和反映反射率的振幅信息。基于波场延拓理论的波动方程保幅地震偏移成像,是在给出正确构造成像位置的同时也给出真实反射振幅的有效完善。基于全声波方程,利用严格的解耦理论进行单程波动保幅分解,得到一个由波场传播项与振幅补偿项构成的,在走时与振幅上满足全声波方程对应的程函方程与输运方程的保幅单程波动方程;利用摄动理论进行单平方根算子渐进展开,推导出基于保幅波动方程的广义高阶屏地震偏移算子方程。模型测试和实际资料处理表明,该方法不但可以凭借更准确的相位归位和散射能量聚焦提高构造成像精度,而且输出了能更正确反映地下反射属性的能量信息,从而可以为更深层次的勘探开发,提供地球物理技术支撑。     

炮检分离浮动基准面处理方法在复杂山地地震成像中的应用——以准噶尔盆地南缘为例

准噶尔盆地南缘山地地表起伏剧烈、地下构造复杂且速度场纵横向变化大,地震数据原始单炮一致性差、信噪比低,这对山前带起伏地表地震成像造成了巨大挑战。为完善复杂条件下低信噪比山地资料的起伏地表处理技术,以提高山地资料地震成像精度,本文提出了保持炮检点一致性的浮动基准面处理方法。通过配套的时间域处理和深度域速度建模技术,实现了炮检点分离的浮动地表处理思路。基于南缘四棵树工区连片深度偏移生产中的实际应用,证明了该处理思路能够有效提高复杂构造成像质量,同时大大缩短处理项目周期,现已成为双复杂地区地震成像的主力技术手段,可以广泛推广到其它地区山前带起伏地表地震成像中。     

碳酸盐岩喀斯特溶洞和裂缝系统的地震模拟与预测

提出了有助于确定碳酸盐岩裂缝带和喀斯特溶洞系统空间分布的地震正演模拟和预测方法。地震正演方法包含两个关键步骤，一是利用弹性波动方程计算含缝洞介质的有效弹性参数，另一步骤是根据有效弹性参数构成的等效地质模型，用单程声学波动方程进行波场延拓。数值模拟和实际地震数据的例子说明，将地震正演剖面与地震偏移剖面相结合进行对比分析有可能识别和预测碳酸盐岩地层中的缝洞发育带。     

地震波传播的褶积微分算子法数值模拟

为了解决地震波场数值模拟中的速度与精度的匹配及局部信息与全局信息的耦合等问题, 该文借助新推出的基于Forsyte广义正交多项式的褶积微分算子, 将计算数学中的Forsyte多项式褶积微分算子应用到地震波传播的数值模拟中.复杂非均匀介质模型数值模拟结果说明了该方法的可行性和优越性.该方法同时具有广义正交多项式褶积微分算子的高精度和有限差分短算子算法的高速度.通过对算子长度的调节及算子系数的优化, 可同时兼顾波场解的全局信息与局部信息.      

基于等效偏移距的偏移成像方法

等效偏移距偏移(EOM)成像方法将克希霍夫时间偏移的双平方根方程转化为一个单平方根方程,使输入采样在没有发生时移的情况下映射到以等效偏移距为变量的共散射点道集上,使散射能量按照双曲线规律分布。然后通过克希霍夫动校正和适用在共散射点道集上的叠加来完成偏移成像。基于等效偏移距叠前地震偏移(EOM)的基本理论,对CSP道集的抽取过程进行模拟实现,并对模型进行速度分析,实现偏移成像。     

基于单程波真振幅分步傅里叶叠前深度偏移方法

这里将单程波真振幅方程与分步傅里叶算子（SSF）相结合,同时还结合了保幅算法和分步傅里叶算法的优点,因此该方法具有计算量小,占内存少,能处理横向变化的速度等优点。并且克服了傅里叶有限差分方法偏移后的振幅都有很大的偏差的不足。与目前广泛应用的常规的分步有限差分叠前深度偏移相比,具有成像精度高,保持地震波动力学特征等优点。在Marmousi模型上成功地进行了真振幅分步傅里叠前深度偏移处理,取得了理想的成像效果。     

偏移速度分析与偏移成像

偏移速度分析和偏移成像是地震资料处理的两个重要组成部分.目前时间偏移技术比较成熟,而深度偏移技术也在逐步完善,时间域偏移成像主要推崇叠前时间偏移法.采用沿层叠加速度分析技术可获得层面上准确的叠加速度,并通过倾角校正、叠前时间偏移和CRP反偏移速度分析逐步优化速度,得到一个符合地质规律、准确的均方根速度场;通过深度偏移方法的研究,总结了建立精确偏移速度场的方法,并提出了一种地震资料处理的思路,即基于射线追踪的Kirchhoff偏移和基于波场延拓的波动方程偏移的结合,使偏移速度分析和偏移成像在应用效果和效率上得到了很大的提高.     

基于各向异性理论的深水区地震资料叠前处理技术

南海北部深水区白云凹陷海底地质构造复杂、速度变化剧烈,各向异性现象普遍存在,多次波极为发育,地震成像难度大。为了更好地提高该深水区地震资料的成像质量,笔者引入了基于各向异性介质的叠前地震资料处理技术,并将其应用到深水白云6-1工区。该技术系列主要由三维广义多次波衰减、各向异性速度分析和偏移成像等技术组成。深水资料应用结果表明：GSMP建立的多次波模型,能有效地压制多次波,使剖面更显著;基于各向异性的高密度速度分析可以更好地描述地质信息;剖面的断层能量在基于各向异性叠前时间偏移处理后得到聚焦,最终该技术系列提高了地震资料的成像质量。     

VTI介质多波速度与各向异性系数求取及应用

针对具有垂直对称轴的横向各向同性(VTI)介质,研究了多波速度与地层各向异性系数之间的关系。采用适应大偏移距以及强各向异性介质的双平方根方程,进行速度和各向异性系数分析,并将批量计算与交互解释相结合,实现了高精度各向异性介质多波速度与各向异性系数提取,并形成了一套利用多波地震资料求取地震波纵、横波速度及各向异性系数的方法和软件。所求取的参数可用于多波地震资料处理,提高多波资料的成像质量,也可为地层地质解释提供参考。实际资料的应用效果证明了该方法和软件的正确性及有效性。     

变参考慢度Born近似傅氏偏移

针对常规Born近似傅氏偏移方法对于剧烈横向变速介质不能精确成像的状况而提出了变参考慢度Born近似傅氏偏移，理论上解决了任意速度变化地质模型的偏移成像问题。此外，为进一步提高复杂地层的成像精度和波场延拓算子的稳定性，对散射波场的计算公式作了改进。将改进的方法运用于盐丘模型的正演和偏移试验，并与常规Born近似偏移方法相比较，可明显看出变参考慢度Born近似傅氏偏移方法在效果上要优于后者，其处理速度横向变化的能力大大增强。     

几种单程波动方程解法的对比

在均匀介质条件下, 比较了ω-k域内Stolt算法和Gazdag相移法、t-x域内部分分式算法、混合域内SSF算法和FFD算法的计算流程, 零炮间距叠后时间偏移的单位脉冲响应、脉冲振幅变化、中心地震道数据的子波振幅变化和脉冲波场频谱特征等体现了波场动力学特性的信息。分析表明, 在均匀介质条件下, ω-k域方法在脉冲响应、振幅变化和频谱特征方面能得到较好的结果, 但该方法的子波分辨率相对其他方法较低, 无法应用于横向速度变化介质问题, 所以复杂地质问题应采用FFD法进行处理。同时以Stolt算法为例具体分析了子波振幅、子波频谱在偏移过程中的变化。     

Seismic depth migration with pseudo-screen propagator

The generalized screen propagator (GSP) is a fast and accurate method for computing wave propagation in laterally varying velocity media. The method solves the acoustic one-way wave equation by shuffling FFTs between frequency-space and frequency-wavenumber domain. The velocity model is decomposed into a vertically varying reference velocity and a laterally varying perturbation. Downward continuation of the wavefield is performed with a phase shift in a background medium followed by a correction that accommodates the lateral velocity variations. All frequencies of interest are migrated and summed to produce an image at each depth level. The screen propagators are applicable to both poststack and prestack depth migrations. We illustrate the method with the SEG-EAGE salt model and Marmousi data set. The image results show that the pseudo-screen propagator is a big improvement over the phase-screen propagator for steeply dipping events in strong-contrast heterogeneous media.     

地震波混合阶褶积算法模拟

基于正反傅立叶变换,提出了地震波模拟的混合阶褶积算法。该方法原理简单、易于实现,结合了有限差分与伪谱法的优点,具有较高的精度和计算效率,适用于地震波场的正演计算,同时还给出了二维地震波场的理论计算实例。计算结果表明,此算法模拟结果正确、精度高、速度快、能适应较为复杂地质模型,并且易于推广到各向异性介质中去。     

地震分析中人工边界处理与地震动输入方法研究

基于柱面波波动方程,推导建立了适用于土-结构地震动力相互作用分析的地震动输入和人工边界的处理方法.其中,地震动的输入是通过在人工边界上施加等效节点力来实现的,等效节点力的大小与入射地震波波速成正比;而人工边界的处理方法使得人工边界条件不仅在时间上是局部的,而且在空间上也是局部的.这种处理方法简单、有效,物理意义清晰,且很容易在有限元法中实现,结合Newmark时间积分是无条件稳定的.为了验证方法的有效性和精度,给出了两个算例,分别用于检验人工边界条件的性能以及地震动输入方法的正确性.算例分析结果表明,所提出的方法是十分有效的.     

黏滞声波高斯束叠前深度偏移

高斯束偏移不仅具有接近于波动方程偏移的成像精度,而且保留了Kirchhoff 积分法高效、灵活的优点,可以对复杂介质准确成像。由于实际地下介质具有黏滞性,因此研究黏滞声波叠前深度偏移具有一定的现实意义。笔者采用高斯束偏移方法对地震数据进行吸收衰减补偿。首先给出共炮域高斯束叠前深度偏移原理;然后在此基础上推导补偿吸收衰减的表达式,校正品质因子Q引起的振幅衰减和相位畸变,实现基于吸收衰减补偿的高斯束叠前深度偏移;最后用两层模型和气云模型对偏移方法进行了测试。结果表明,在考虑地下介质的黏滞性时,黏滞声波高斯束叠前深度偏移比声波高斯束叠前深度偏移具有更高的成像分辨率。