基于Marchenko理论一步法压制层间多次波

传统基于Marchenko理论自聚焦法压制层间多次波技术需要初始下行聚焦函数的估计,再进行多次波压制的自适应相减.论文基于传统方法的研究成果,提出了一种无需先验信息和自适应相减步骤的层间多次波压制方法.当聚焦点在实际反射层的反射点上时,上行格林函数为聚焦点上的脉冲源在表面所产生的上行波场,该反射层的反射将是上行格林函数...     

基于AVO正演模拟的自由表面多次波压制方法

针对SRME和Radon变换压制自由表面多次波的近道多次波残留和伤害有效波等问题,主要采用基于AVO属性分析的多次波剔除法对近道多次波进行压制,该方法在迭代更新过程中,受近道多次波影响容易出现更新误差,同时迭代更新降低计算效率.本文对基于AVO属性分析的多次波剔除法进行改进,利用中远道地震数据进行约束,通过AVO正演模...     

结合Marchenko聚焦的非线性最小平方基准面延拓方法

基准面延拓旨在通过模型驱动或数据驱动的方法消除上覆介质对地震信号造成的畸变或影响, 进而重构基准面上的虚拟观测数据.针对非线性波场效应(如层间多次波), 模型驱动的基准面延拓方法对地下速度模型精度提出严苛的要求, 实用性较差.如果检波器布设面与基准面恰好重合, 数据驱动的虚源法能够自适应干涉消除波场在上覆介质中的非线性传播效应.然而, 由于反射地震观测系统不够完备, 通常会在重构数据中引入串扰噪声.针对上覆介质具有强非线性波场效应、基准面位于观测面与目的层之间任意深度的一般情况, 本文将Marchenko方法和最小平方基准面延拓(LSR)相结合, 先估算地表到基准面的上、下行格林函数, 而后利用稀疏约束LSR方法克服观测系统和上覆地层的影响, 在压制层间多次波的同时高精度重构基准面上的虚拟观测数据.层状模型合成数据实验展示了非线性波路径的干涉抵消机制, 透镜体模型合成数据实验检验了该方法在压制上覆异常体多次散射效应、校正深部反射信号畸变的有效性.

     

基于Marchenko理论的成像域与数据域层间多次波消除方法对比

Marchenko成像可以理解为一种基于Marchenko理论的成像域层间多次波消除方法. 该方法利用Marchenko方程组将地表观测的地震数据准确的延拓到地下的成像点, 得到震源位于地下成像点、检波点位于地表的上下行格林函数, 进而采用多维反褶积获得该成像点不含多次波假象的像. 基于Marchenko方程组所推导出的数据域多次波消除方法可以准确的预测并消除所有阶数的层间多次波. 该方法不需要任何速度信息和匹配相减, 避免了传统多次波消除方法对一次波造成的损害. 基于数据域多次波消除之后的数据, 可以获得地下结构不含多次波假象的像. 本文对基于Marchenko理论的成像域及数据域层间多次波消除方法进行了对比分析, 给出了详细的理论及数值模型试算结果, 并对两种方法的优缺点进行了讨论.

     

基于稀疏反演三维表面多次波压制方法

三维表面多次波压制是海洋地震资料预处理中的重要研究课题,基于波动理论的三维表面多次波压制方法(3D SRME)是数据驱动的方法,理论上来说,可有效压制复杂构造地震数据表面多次波.但该方法因对原始地震数据采集要求高而很难在实际资料处理中广泛应用.本文基于贡献道集的概念,将稀疏反演方法引入到表面多次波压制中,应用稀疏反演代替横测线积分求和,无需对横测线进行大规模重建,进而完成三维表面多次波预测,这样可有效解决实际三维地震数据横测线方向稀疏的问题.基于纵测线多次波积分道集为抛物线的假设,为保证预测后三维表面多次波和全三维数据预测的多次波在运动学和动力学特征上基本一致,文中对预测数据实施基于稳相原理的相位校正.理论模型和实际数据的测试结果表明,本文基于稀疏反演三维表面多次波压制方法可在横测线稀疏的情况下,有效压制三维复杂介质地震资料中的表面多次波,从而更好地提高海洋地震资料的信噪比,为高分辨率地震成像提供可靠的预处理数据保障.     

基于数据一致性预测与压制自由表面多次波--理论研究与试处理

从声波波场理论出发,首先在不考虑自由表面反射的情况下研究了有效波波场传播及反射规律,继而考虑自由表面反射的情况,并详尽推导了压制自由表面多次波的正演和反演数学模型。从物理本质上对自由表面多次波的产生、预测及压制做了深入的剖析,讨论了基于数据一致性时空域褶积的多次波预测方法原理,实现了自由表面相关多次波自适应压制的迭代方法过程。复杂起伏海底的理论模型试算实例表明：在多次波产生系统非常复杂的情况下,本方法仍能非常有效地压制各种与自由表面相关的多次波,同时还能很好地保持有效波的能量,取得了很好的压制多次波效果。     

深水陡坡带绕射多次波压制方法研究

在深水陡坡带和复杂海底构造区域,即使采用长电缆,由绕射产生的自由表面相关多次波因为传播方向的异变,仍然难以避免在自由表面发生下行反射的反射点位置位于电缆长度范围之外,造成预测多次波的近偏移距信息不足或者负偏移距信息缺失等问题,以致不能准确地预测和压制陡坡带的绕射多次波.本文在常规反馈迭代自由表面多次波压制方法的基础上做出相应改进,将邻炮的炮记录依照炮点、检波点互换原理引入到反馈迭代循环当中,补充绕射多次波预测所需要的近偏移距和负偏移距信息.经过Sigsbee2b模型和实际资料测试,达到较好压制深水陡坡带的自由表面相关绕射多次波的目地.     

基于数据规则化和稀疏反演的三维表面多次波压制方法

表面多次波是海洋地震勘探中的主要问题.目前,二维数据驱动的表面多次波压制技术(SRME)已经比较成熟,并且已经成为工业界压制海洋表面多次波的主流方法.但是由于二维SRME算法没有考虑横测线方向上多次波的贡献,导致在处理实际三维海洋资料时存在比较大的误差.将二维SRME算法扩展到三维空间后可以得到三维SRME算法,但是由于目前实际采集的三维海洋资料的观测系统存在拖缆漂移,而且横测线方向采样过于稀疏,直接应用三维SRME算法无法准确预测表面多次波.本文提出的通过数据规则化配合稀疏反演的三维表面多次波压制方法能够解决这种实际资料和三维SRME算法之间的矛盾.本文通过研究数据规则化与反规则化技术,使得数据分布满足三维SRME的要求;通过研究稀疏反演技术,有效解决了横测线方向采样稀疏对于多次波预测的影响,三维实际海洋资料的应用结果验证了方法的有效性和可行性.     

基于反馈环的多次波压制方法研究

基于波动方程理论压制多次波,综合考虑了多次波传播的运动学特性和动力学特性,避免了非波动法解决问题所带来的局限性,已成为多次波去噪的发展趋势.首先对多次波压制国内外研究现状做了简单的阐述.在此基础上,给出了反馈环波动方程预测减去法的方法原理.介绍了抛物线拉东(Radon)变换外推近道数据、谱模拟法提取震源子波、均衡伪多道匹配和多次波减去法效果评价等四大关键技术的思路及主要计算公式.用JAVA语言自行开发了相应的模块,并应用射线法多次波模型和SEG/EAGEPluto模型,对以上方法技术进行了检验,去噪效果良好.对海上某多次波特发育实际资料进行了处理,取得了满意的效果.     

利用最小平方基准面延拓压制自由表面多次波Ⅰ:海底观测地震数据

海底多分量地震采集技术被广泛用于海洋地质调查、油气资源勘探以及地震与海啸等灾害监测领域.然而,自由表面效应导致地震观测数据中多次波干扰很强,影响一阶散射理论基础上的速度建模和偏移成像.本文基于数据驱动的地震干涉原理,在最小平方反演理论框架下提出消除自由表面和上覆介质影响的基准面延拓方法,结合海森算子与点扩散函数揭示其工作机理.针对海底分离的上、下行波记录,在无需进行多次波分离或分阶提取情况下,迭代反演获得在海底虚拟观测的叠前地震数据.理论模型合成数据实验结果表明,该方法能够有效压制源端自由表面与海水层的影响,提高有效信号的分辨率,进而改善地震成像质量.

     

基于广义的炮偏移方法实现地表多次波和一次波联合成像

多次波是地下反射层的多次波反射,也蕴含了地下反射界面的信息,因此并不是绝对地只能被当做噪音来处理.为实现对地下构造的准确成像,本文基于广义概念上的炮偏移成像算法,对常规一次波偏移方法从用于向下延拓的上、下行场以及成像条件方面进行了改进,将同时含有表层多次波的炮记录与脉冲震源之和作为下行延拓的震源波场,将同时含有表层多次波和一次波的炮记录作为接收波场,并应用2D反褶积成像条件进行成像,实现了表层多次波和一次波联合成像.在数值试验部分,通过对相对简单的五层模型和Sigsbee2b模型数据的数值测试表明,该方法可以很大程度上压制常规成像方法中由多次波带来的成像噪音.     

基于波动方程三维表面多次波预测方法研究

与传统的二维表面多次波预测算法相比,基于波动方程的全三维表面多次波预测方法无需对地下介质做简单近似,其更符合地震波在地下介质中传播的真实状况,是地震资料处理中解决多次波预测问题的强有力工具.本文从三维多次波预测的基本理论出发,给出了全三维多次波预测算法的预测矩阵表示、计算方法以及实现条件,采用GPU(图形处理器)加速全三维表面多次波预测,较传统的CPU串行计算,GPU并行预测表面多次波的计算效率约提高165倍.文中分别利用二维和三维表面多次波预测算法对理论模拟的含表面多次波的三维地震数据进行多次波预测计算,对比分析结果表明,相比于二维算法,文中所述的基于波动方程的全三维表面多次波预测效果明显改善,其计算精度更高,辅以合理有效的自适应相减算法,可获得高精度的地震勘探资料表面多次波压制数据.     

基于单程波偏移算子的地表相关多次波成像

在常规地震资料处理中,多次反射波被视为噪声并从地震数据中去除,以免在之后的地震资料解释中造成误解.而事实上,多次波也是地震信号,是照明波场的一部分,能够对地下构造成像的精度做出贡献.本文分析了多次波在传统单程波叠前深度偏移中产生构造假象的机制和表现,为实现基于单程波偏移算子的多次波成像,修改了单程波叠前深度偏移的边界条件,即将输入的震源波场用包含多次波的记录来替代,输入的记录波场用预测出的表层相关多次波来替代,实现了基于单程波偏移算子的地表相关多次波成像,并从理论上给出了其成像依据.通过基于二范式最小能量差原则求取的匹配因子,将多次波成像结果与一次波成像结果进行匹配叠加,应用多次波成像来弥补一次波成像的不足.简单模型验证了基于单程波偏移算子的多次波成像方法的有效性,最后对Sigsbee2B模型进行了一次波与多次波联合成像试算,盐边界高陡构造成像质量得到了明显改善.     

基于迭代反演的层间多次波压制方法

层间多次波与有效波在走时、频率和叠加速度上差异较小,因此层间多次波压制常难以获得理想效果.本文提出一种基于迭代反演的层间多次波压制方法（MSI,Multiple Suppression Inversion）,该方法以共聚焦点（CFP,Common Focus Point）层间多次波压制理论为基础,通过构建卷积因子,将层间多次波压制转变为迭代反演的问题,直接利用观测的地震数据进行迭代反演计算,进而完成多次波压制.MSI方法避免了共聚焦点方法中构建CFP道集的聚焦运算,大幅降低层间多次波预测的计算成本;同时该方法为全数据驱动方法,无需地下介质任何先验信息,算法容易实现.模型数据测试表明,本文提出的方法可有效压制层间多次波,而且对有效波也具有很好的保幅性.

     

基于泰勒多项式及曲面样条函数的地电场变化分析

选取江苏及邻区地电台网13个台站地电场观测分钟值,采用泰勒多项式和曲面样条函数模型建立地电日变化模型,使用模型曲线描述不同年份、同农历日期、同等磁情天气的实测地电日变化,结果表明:使用分钟值曲线建立的地电日变化模型,更能清晰反映日变化特点,并可反映一天当中台站观测环境、测量系统等变化引起的地电变化;使用同一区域地电台网中若干台站同方向测道的样本曲线,建立各测道电场变化模型,可反映其自身特征及共性特征。     

基于吸收衰减补偿的多分量高斯束逆时偏移

高斯束逆时偏移结合了射线类偏移的高计算效率和波动方程逆时偏移的高精度,能很好地处理焦散点、大倾角成像问题,并且具有面向目标成像的能力.多分量地震资料的偏移技术可以对地下复杂构造进行更准确的成像,由于实际地下介质具有黏滞性,研究黏弹性叠前逆时偏移具有一定的现实意义.本文采用高斯束逆时偏移方法对多分量地震数据进行吸收衰减补偿,首先分别给出纵波和转换波共炮域高斯束叠前逆时偏移方法原理,在此基础上推导补偿吸收衰减的表达式,校正Q引起的振幅衰减和相位畸变,实现基于吸收衰减补偿的多分量高斯束叠前逆时偏移.数值模型的测试结果显示,在考虑地下介质的黏滞性时,本文方法具有更高的成像分辨率.