基于波场分解的一次波和多次波分离及同时成像方法技术

本发明专利技术公开基于波场分解的一次波和多次波分离及同时成像方法，涉及地震信号处理技术领域。本发明专利技术通过改进双程波方程波场分解方案，实现高效的分离一次波和不同阶次的多次波，消除一次波和不同阶多次波成像中的串扰噪音，并同时实现一次波和多次波的高效成像，本发明专利技术提供的方法为：在采集面将上行波场分解为一次波和不同阶次的多次波，将下行波场分解为不同阶次的多次波虚震源；基于这种广义的上/下行波场分解基础，使用改进的双程波方程波场深度延拓方案进行一次波和不同阶次的多次波的同时深度延拓和成像，仅需计算一次，就可实现多种波场类型高效的成像。现多种波场类型高效的成像。现多种波场类型高效的成像。

【技术实现步骤摘要】
基于波场分解的一次波和多次波分离及同时成像方法


[0001]本专利技术涉及地震信号处理
，涉及基于波场分解的一次波和多次波分离及同时成像方法。

技术介绍

[0002]针对地震信号处理技术中的一次波和多次波分离及成像：
[0003]现有技术中：
[0004]CN105334537A用交替分裂Bregman迭代算法求解对一次波施加稀疏约束的优化问题，实现3D匹配滤波器的估计，并利用估计的3D匹配滤波器自适应分离3D数据窗口中的一次波和多次波；
[0005]CN112083492A通过输入速度场、Q补偿参数场及实际观测炮记录，建立观测系统，输入海底崎岖地表高程和观测系统文件，生成正交贴体网格，并将速度场和Q补偿参数场变换到曲坐标系，基于深海环境特征计算崎岖海底条件下全路径Q补偿波场的正演模拟算子、一次波与多次波联合伴随算子与反偏移算子，构建全路径Q补偿一次波与多次波联合成像的目标泛函并求取梯度，计算曲坐标系下的全路径补偿一次波与多次波联合成像结果并变换至笛卡尔坐标系下，输出成像结果；
[0006]CN106324669B结合聚焦变换的加权互相关思想和SRME的加权褶积思想，先用SRME实现一次波与各阶表层多次波整体的分离，再用正聚焦变换实现各阶表层多次波的降阶，接着用非平稳匹配滤波提取准地震记录，然后用SRME提取准地震记录中的准一次波，最后用反聚焦变换实现准一次波的升阶返回原始数据域，获取分离的各阶表层多次波；
[0007]综上所述，现有技术对于一次波和多次波分离及成像方法，通常情况下，当分别单独偏移不同阶数的多次波时，由于需要进行多次偏移计算，会导致计算成本显著增加，这尤其对基于双程波波动方程的深度偏移提出了挑战，针对多次波偏移中通常面临的两个挑战：串扰伪影干扰和计算成本增加；本专利技术基于双程波方程波场延拓实现高效的分离一次波和不同阶次的多次波及成像。

技术实现思路

[0008]本专利技术的目的在于，提供基于波场分解的一次波和多次波分离及同时成像方法，通过改进双程波方程波场深度延拓方案，能够以非常有限的计算成本和高效的方式同时完成一次波和不同阶数的多次波的成像。
[0009]为实现上述技术目的，达到上述技术效果，本专利技术是通过以下技术方案实现：
[0010]本专利技术的第一目的在于，提供一种双程波方程波场深度延拓方法，所述方法如下：
[0011]所述双程波方程波场深度延拓方法如下：
[0012]在频率空间域中，二维声波方程可表示为：
[0013][0014]其中，表示压力波场，ω表示角频率x和z分别表示水平和垂直空间坐标；v(x,z)表示二维速度模型；
[0015]定义并施加两个边界条件，这些条件在初始深度z＝z0处被定义；通常把及其导数作为方程(1)的边界条件；从z＝z0开始的双程波方程深度延拓方案可以表示为以下矩阵
‑
向量形式：
[0016][0017]其中，表示垂直波数，z为深度，Δz为延拓步长；
[0018]通过引入压力波场来替换方程(2)中的压力波场导数；参数的定义为：
[0019][0020]其中，I是虚数单位；
[0021]此外，在实际多分量数据采集中，通常可记录到速度垂直分量v
z
(x,z,ω)，基于该速度垂直分量数据，也可以使用以下公式来计算波场
[0022][0023]其中，I是虚数单位；
[0024]方程(2)中的矩阵
‑
向量形式可以改写为如下形式：
[0025][0026]与方程(2)相比，方程(5)中改写后的矩阵
‑
向量形式使得上行波场和下行波场的分离更加简单，可以通过加减法计算高效完成：
[0027][0028]其中，和分别代表上行波场和下行波场；上下波场分离是使用一阶泰勒公式并忽略高阶项所得，为保持能量守恒，加1/2；
[0029]因此，双压力波场和可以通过与方程(6)类似的计算高效地完成，具体如下：
[0030][0031]本专利技术的另一目的在于，提供基于双程波方程波场分解的广义一次波和多次波分
离方法，以解决串扰伪影问题，所述方法为：
[0032]根据一次波和多次波传播的几何模式，总波场包含两个分量：向下传播的震源波场和向上传播的检波器波场即总下行波场和总上行波场其表达式分别表示为：
[0033][0034][0035]下行波场作为一个广义源，包含两种类型的源：真震源(或S)和多次波虚震源D；上行波场包括一次波和多次波M；分别表示第i阶下行和上行波场，其中多次波序数i＝1,2,...,N；N为多次波的阶数；
[0036]从方程(8)和方程(9)可得：
[0037][0038][0039]其中，
[0040][0041][0042]1其中，S即真震源D是上界面的反射下来的多次波，M是下界面反射上去的多次波；
[0043]如果直接将相互成像条件应用于上行波场和下行波场，得到的成像不仅包含所需要的真正偏移信息，还包含不需要的串扰噪音；
[0044][0045]其中，右边方程的第一项和第二项分别对应真正偏移成像结果和串扰噪音成像结果；
[0046]针对常规互相关成像中的串扰噪音问题，本专利技术在采集面上将上行波场分解为一次波和不同阶次的多次波，将下行波场分解为一次波和不同阶次的多次波虚震源；
[0047]基于上行波场和下行波场分离基础，使用双程波方程波场深度延拓方法进行一次波和不同阶次的多次波的同时深度延拓和成像，避免成像中的串扰噪音，仅需一次计算，即可实现高效的一次波和多次波同时成像。
[0048]进一步的，将上行波场和下行波场分解为一次波和不同阶次的多次波，假设震源波场已知，并遵循上下分离的方法，应用两个压力参数和或等效的上行波场
和下行波场
[0049]其中，上行波场的分离，从震源、地球响应和反射之间的关系式开始：
[0050][0051]其中：
[0052][0053]M＝G
·
D (17)
[0054]从方程式(12)、(13)、(17)可以得到：
[0055][0056]G是地球响应，从方程式15
‑
18可知，G由以下四个方程式之一表示出来：
[0057][0058]使用方程式(19)的第一个公式和方程(16)，我们可以推到得到如下的和之间的关系：
[0059][0060][0061]其中：
[0062][0063][0064]上述D
′
代表反褶积后的下行波，上述A代表频域中震源特征；
[0065]为了将上行波场分解成一次波和不同阶次的多次波，将方程式(21)右边按级数展开：
[0066][0067][0068]其中，N为多次波的阶数；
[0069]方程式(24)、(25)表达的和之间的关系，可以被证明与使用只有上行波场导出的其它方法是等效的；本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于双程波方程波场深度延拓方法的广义一次波和多次波分离方法，其特征在于：采集面上将上行波场分解为一次波和不同阶次的多次波，将下行波场分解为一次波和不同阶次的多次波虚震源；基于上行波场和下行波场分离基础，使用双程波方程波场深度延拓方法进行一次波和不同阶次的多次波的同时深度延拓和成像，避免成像中的串扰噪音。2.如权利要求1所述的基于双程波方程波场深度延拓方法的广义一次波和多次波分离方法，其特征在于：所述双程波方程波场深度延拓方法如下：在频率空间域中，二维声波方程可表示为：其中，表示压力波场，ω表示角频率x和z分别表示水平和垂直空间坐标；v(x,z)表示二维速度模型；定义并施加两个边界条件，这些条件在初始深度z＝z0处被定义；通常把及其导数作为方程(1)的边界条件；从z＝z0开始的双程波方程深度延拓方案可以表示为以下矩阵
‑
向量形式：其中，表示垂直波数，z为深度，Δz为延拓步长；通过引入压力波场来替换方程(2)中的压力波场导数；参数的定义为：方程式(3)中，I为虚数单位；此外，在实际多分量数据采集中，通常可记录到速度垂直分量v
z
(x,z,ω)，基于该速度垂直分量数据，也可以使用以下公式来计算波场(x,z,ω)，基于该速度垂直分量数据，也可以使用以下公式来计算波场其中，ρ为介质密度；方程(2)中的矩阵
‑
向量形式可以改写为如下形式：I为虚数单位；与方程(2)相比，方程(5)中改写后的矩阵
‑
向量形式使得上行波场和下行波场的分离更加简单，可以通过加减法计算高效完成：
其中，和分别代表上行波场和下行波场；上下波场分离是使用一阶泰勒公式并忽略高阶项所得，为保持能量守恒，加1/2；因此，双压力波场和可以通过与方程(6)类似的计算高效地完成，具体如下：3.如权利要求2所述的基于双程波方程波场深度延拓方法的广义一次波和多次波分离方法，其特征在于：所述方法为：根据一次波和多次波传播的几何模式，总波场包含两个分量：向下传播的震源波场和向上传播的检波器波场即总下行波场和总上行波场其表达式分别表示为：其表达式分别表示为：下行波场作为一个广义源，包含两种类型的源：真震源(或S)和多次波虚震源D；上行波场包括一次波和多次波M；分别表示第i阶下行和上行波场，其中多次波序数i＝1,2,...,N；N为多次波的阶数；从方程(8)和方程(9)可得：从方程(8)和方程(9)可得：其中，其中，其中，S即真震源D是上界面的反射下来的多次波，M是下界面反射上去的多次波；如果直接将相互成像条件应用于上行波场和下行波场，得到的成像不仅包含所需要的真正偏移信息，还包含不需要的串扰噪音；其中，右边方程的第一项和第二项分别对应真正偏移成像结果和串扰噪音成像结果。4.如权利要求3所述的基于双程波方程波场深度延拓方法的广义一次波和多次波分离
方法，其特征在于：将上行波场和下行波场分解为一次波和不同阶次的多次波，假设震源波场已知，并遵循上下分离的方法，应用两个压力参数和或等效的上行波场和下行波场其中，上行波场的分离，从震源S、地球的脉冲响应G和上界面的反射多次波D之间的关系式开始：其中：M＝G
·
D (17)从方程式(12)、(13)、(17)可以得到：G是地球响应，从方程式15
‑
18可知，G由以下四个方程式之一表示出来：使用方程式(19)的第一个公式和方程(16)，我们可以推到得到如下的和之间的关系：系：其中：其中：上述D
′
代表反褶积后的下行波，A代表频域中震源特征。5.如权利要求4所述的基于双程波方程波场深度延拓方法的广义一次波和多次波分离方法，其特征在于：为了将上行波场分解成一次波和不同阶次的多次...

【专利技术属性】
技术研发人员：尤加春，任强，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：