CN103308944B - 一种海洋地震资料处理中鬼波衰减的方法及装置 - Google Patents
一种海洋地震资料处理中鬼波衰减的方法及装置 Download PDFInfo
- Publication number
- CN103308944B CN103308944B CN201310187398.XA CN201310187398A CN103308944B CN 103308944 B CN103308944 B CN 103308944B CN 201310187398 A CN201310187398 A CN 201310187398A CN 103308944 B CN103308944 B CN 103308944B
- Authority
- CN
- China
- Prior art keywords
- ripple
- terrible
- sided
- reflection
- bilateral
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
本发明提供了一种海洋地震资料处理中鬼波衰减的方法及装置，包括以下步骤：将原始地震记录中的直达波和海面反射波切除，得到记录A；在所述原始地震记录中提取出海面反射波M和反射层一次反射波P；将所述海面反射波M和所述反射层一次反射波P褶积，得到该反射层相关的震源侧的单侧鬼波G1；将所述海面反射波M与所述震源侧的单侧鬼波G1褶积，得到震源和检波器的双侧鬼波G2；将所述单侧鬼波G1、双侧鬼波G2和记录A进行能量匹配，得到匹配后的单侧鬼波G1m和双侧鬼波G2m；在所述记录A中减去所述匹配后的单侧鬼波G1m和双侧鬼波G2m，得到衰减鬼波后的记录Sdghost；由于该方案通过纯数据驱动进行鬼波预测和衰减，因而不依赖与速度场的精度，提高了海洋地震资料的处理质量。
Description
技术领域
本发明涉及地震资料处理技术领域，特别涉及一种海洋地震资料处理中鬼波衰减的方法及装置。
背景技术
随着油气勘探的深入，海洋油气勘探受到了越来越多的重视，海洋地震资料处理成果是海洋油气勘探的重要依据，而多次波是影响海洋地震资料处理的最为关键因素之一。由于海平面和海底都存在强速度差的反射界面，反射系数很大，所以表面相关多次波非常发育，严重的影响了地震资料处理的可靠性，干扰了地震资料解释人员对构造的判断。海洋地震资勘探中，激发震源空气枪与接收电缆都在海面以下某个深度，震源波场传播至海面后又发生反射被检波器接收到的波场称为鬼波，鬼波是一种普遍存在于海洋地震资料中的表面相关多次波，难以在地震资料处理中衰减。
为了消除鬼波对地震资料的影响，一种思路是从资料采集中入手，海底电缆OBC技术将检波器置于海底，通过两种检波器记录数据的结合来压制鬼波。另一种思路是在处理中进行衰减，一般是通过子波整形类的方法，但是当震源与比较深时，误差会变大。当前工业界普遍应用的表面相关多次波压制方法(SRME)，通过记录本身的自褶积来预测自由表面相关多次波。鬼波分为震源一侧鬼波、检波点一侧鬼波和震源与检波点双侧鬼波，传统SRME方法不能直接应用到这种鬼波压制中，而通过采集方式的改变来压制鬼波会较大的增加工业成本。Weglein(1997)等人提出完全数据驱动型的逆散射去层间多次波方法，该方法可从逆散射子序列中一次预测所有的层间多次波，但运行成本非常高，在工业界至今还没有应用。Jakubowicz(1998)提出了层相关层间多次波预测方法，但是该方法主要是针对地下强反射界面层间多次波的预测。
鬼波的衰减方法是海洋地震资料处理普遍存在的难题，随着油气勘探开发的不断深入，该项技术的研究有着重要的应用潜力，对提高海洋地震资料处理品质有着重要意义，因而亟需一种行之有效的鬼波预测和衰减的方法。
发明内容
为克服现有技术中存在的问题，本发明提供一种海洋地震资料处理中鬼波衰减的方法及装置，以实现纯数据驱动进行鬼波预测和衰减，从而不依赖与速度场的精度，提高海洋地震资料的处理质量。
为达到上述目的，本发明一方面提供了一种海洋地震资料处理中鬼波衰减的方法，包括：
将原始地震记录中的直达波和海面反射波切除，得到记录A；
在所述原始地震记录中提取出海面反射波M和反射层一次反射波P；
将所述海面反射波M和所述反射层一次反射波P褶积，得到该反射层相关的震源侧的单侧鬼波G1；
将所述海面反射波M与所述震源侧的单侧鬼波G1褶积，得到震源和检波器的双侧鬼波G2；
将所述单侧鬼波G1、双侧鬼波G2和记录A进行能量匹配，得到匹配后的单侧鬼波G1m和双侧鬼波G2m；
在所述记录A中减去所述匹配后的单侧鬼波G1m和双侧鬼波G2m，得到衰减鬼波后的记录Sdghost。
本发明另一方面提供了一种海洋地震资料处理中鬼波衰减的装置，包括：
记录A提取单元，用于将原始地震记录中的直达波和海面反射波切除，得到记录A；
海面反射波提取单元，用于在所述原始地震记录中提取出海面反射波M；
一次反射波提取单元，用于在所述原始地震记录中提取出反射层一次反射波P；
单侧鬼波预测单元，与所述海面反射波提取单元机一次反射波提取单元相连接，用于将所述海面反射波M和所述反射层一次反射波P褶积，得到该反射层相关的震源侧的单侧鬼波G1；
双侧鬼波预测单元，与所述海面反射波提取单元及单侧鬼波预测单元相连接，用于将所述海面反射波M与所述震源侧的单侧鬼波G1褶积，得到震源和检波器的双侧鬼波G2；
能量匹配单元，与所述记录A提取单元、单侧鬼波预测单元及双侧鬼波预测单元相连接，用于将所述单侧鬼波G1、双侧鬼波G2和记录A进行能量匹配，得到匹配后的单侧鬼波G1m和双侧鬼波G2m；
鬼波衰减单元，与所述记录A提取单元及能量匹配单元相连接，用于在所述记录A中减去所述匹配后的单侧鬼波G1m和双侧鬼波G2m，得到衰减鬼波后的记录Sdghost。
本发明能够达到的有益效果：通过纯数据驱动型的鬼波预测和衰减，能够不依赖与速度场的精度，提高了海洋地震资料的处理质量。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1是本发明的原理示意图；
图2是具体实施方式的简单模型速度场；
图3是本发明实施例一种海洋地震资料处理中鬼波衰减的方法的流程图；
图4是本发明实施例一种海洋地震资料处理中鬼波衰减的装置的结构图；
图5中，a是激发与接收在深度100米处正演模拟得到的原始单炮记录，b是切除直达波与海面反射波的记录；
图6中，a是近偏移距部分的海面反射波，b是拾取出的海底一次反射波；
图7中，a是预测出的单侧鬼波，b是预测出的双侧鬼波；
图8中，a是切除直达波与海面反射波的记录，b是衰减完鬼波后的记录。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
如图1所示，是本发明的原理示意图，双侧鬼波传播路径是ABCDEFG，虚线处为激发与接收面，ABCDE与CDEFG是单侧鬼波，因此单侧鬼波的传播可以认为是海面反射波ABC与一次反射波CDE的褶积，双侧鬼波可以认为是单侧鬼波ABCDE与海面反射波EFG的褶积。
图2是是测试模型的速度场，观测面位于海平面以下100米深度处，海水速度为1500m/s，海水下层的速度是2500m/s；
先从原始地震记录中提取出海面反射波与地下反射层的一次反射波，通过将海面反射波与地下反射层的一次反射波褶积得到震源一侧或检波一侧的单侧鬼波，然后在通过海面反射波与预测出的单侧鬼波褶积预测出震源侧和检波器侧的双侧鬼波，最后应用匹配相减的方法将预测出的鬼波从地震记录中衰减掉。
具体通过以下实施例进行说明：
实施例一
如图3所示，为本实施例一种海洋地震资料处理中鬼波衰减的方法的流程图，包括以下步骤：
步骤301，将原始地震记录中的直达波和海面反射波切除，得到记录A；
如图5所示，
图5a是正演模拟得到的单炮记录，存在直达波、海面反射波、一次反射波、单侧鬼波和双侧鬼波，即原始地震记录；
图5b是切除直达波和海面反射波的记录波场，即记录A。
步骤302，在所述原始地震记录中提取出海面反射波M和反射层一次反射波P；
如图6所示，
图6a是提取出的海面反射波M，由于在大偏移距的地方，海面反射波与直达波重叠在一块，所以只提取处双曲规律部分，也就是近偏移距部分；
图6b是提取出的地下反射层一次反射波P。
步骤303，将所述海面反射波M和所述反射层一次反射波P褶积，得到该反射层相关的震源侧的单侧鬼波G1；
褶积的具体算法如下：
G1(xs,wr,f)＝M(xs,ws,f)·P(ws,wr,f)；
其中，M(xs,ws,f)是xs位置震源，ws位置接收到的海面反射波；
P(ws,wr,f)是ws位置震源，wr位置接收到的反射层一次反射波；
G1(xs,wr,f)是xs位置震源，wr位置接收到的单侧鬼波；
f代表频率域。
步骤304，将所述海面反射波M与所述震源侧的单侧鬼波G1褶积，得到震源侧和检波器侧的双侧鬼波G2；
褶积的具体算法如下：
G2(xs,xr,f)＝G1(xs,wr,f)·M(wr,xr,f)；
其中，M(wr,xr,f)是wr位置震源，xr位置接收到的海面反射波；
G2(xs,xr,f)为xs位置震源，xr位置接收到的震源与检波器的双侧鬼波。
如图7所示，
图7a是通过海面反射波与一次记录的褶积得到的单侧鬼波G1，由于海面反射波在大偏移距处被切除，所以预测出的单侧鬼波在大偏移距的地方会有一定的误差；
图7b是单侧鬼波与海面反射波的褶积得到的双侧鬼波G2。
步骤305，将所述单侧鬼波G1、双侧鬼波G2和记录A进行能量匹配，得到匹配后的单侧鬼波G1m和双侧鬼波G2m；
匹配的具体算法如下：
G1m＝F1·G1；
G2m＝F2·G2；
其中，F1与F2是基于最小二乘法在记录A中得到的匹配滤波因子。
最小二乘法匹配法是通过求取一匹配因子使得输入的两个信号在同一能量级，对于F1来说分别将单侧鬼波G1的自相关和单侧鬼波G1与记录A的互相关带入Toeplitz矩阵，通过Levinson算法求解该矩阵得到匹配因子F1。对于F2来说，分别将单侧鬼波G2的自相关和单侧鬼波G2与记录A的互相关带入Toeplitz矩阵，通过Levinson算法求解该矩阵得到匹配因子F2。
步骤306，在所述记录A中减去所述匹配后的单侧鬼波G1m和双侧鬼波G2m，得到衰减鬼波后的记录Sdghost；具体算法如下：
Sdghost＝A-G1m-G2m；
其中，Sdghost是衰减鬼波后的记录。
如图8所示，
图8a是切除直达波与海面反射波的记录，
图8b衰减完鬼波后得到的有效反射波记录，单侧与双侧鬼波都得到了较好的预测与衰减。
本实施例可以达到以下有益效果：通过纯数据驱动型的鬼波预测和衰减，能够不依赖与速度场的精度，提高了海洋地震资料的处理质量。
实施例二
如图4所示，为本发明实施例一种海洋地震资料处理中鬼波衰减的装置的结构图，包括：记录A提取单元401、海面反射波提取单元402、一次反射波提取单元403、单侧鬼波预测单元404、双侧鬼波预测单元405、能量匹配单元406以及鬼波衰减单元407。
该装置工作时，需要事先将原始地震记录的资料送入记录A提取单元401、海面反射波提取单元402和一次反射波提取单元403。其中，
记录A提取单元401，用于将原始地震记录中的直达波和海面反射波切除，得到记录A；
海面反射波提取单元402，用于在所述原始地震记录中提取出海面反射波M；
一次反射波提取单元403，用于在所述原始地震记录中提取出反射层一次反射波P；
单侧鬼波预测单元404，与所述海面反射波提取单元机一次反射波提取单元相连接，用于将所述海面反射波M和所述反射层一次反射波P褶积，得到该反射层相关的震源侧的单侧鬼波G1；具体计算方法如下：
G1(xs,wr,f)＝M(xs,ws,f)·P(ws,wr,f)；
其中，M(xs,ws,f)是xs位置震源，ws位置接收到的海面反射波；所述单侧鬼波预测单元根据以下公式进行计算得到单侧鬼波：
P(ws,wr,f)是ws位置震源，wr位置接收到的反射层一次反射波；
G1(xs,wr,f)是xs位置震源，wr位置接收到的单侧鬼波；
f代表频率域。
双侧鬼波预测单元405，与所述海面反射波提取单元及单侧鬼波预测单元相连接，用于将所述海面反射波M与所述震源侧的单侧鬼波G1褶积，得到震源和检波器的双侧鬼波G2；具体计算方法如下：
G2(xs,xr,f)＝G1(xs,wr,f)·M(wr,xr,f)；
其中，M(wr,xr,f)是wr位置震源，xr位置接收到的海面反射波；
G2(xs,xr,f)为xs位置震源，xr位置接收到的震源与检波器的双侧鬼波。
能量匹配单元406，与所述记录A提取单元、单侧鬼波预测单元及双侧鬼波预测单元相连接，用于将所述单侧鬼波G1、双侧鬼波G2和记录A进行能量匹配，得到匹配后的单侧鬼波G1m和双侧鬼波G2m；具体计算方法如下：
G1m＝F1·G1；
G2m＝F2·G2；
其中，F1与F2是基于最小二乘法在记录A中得到的匹配滤波因子。
最小二乘法匹配法是通过求取一匹配因子使得输入的两个信号在同一能量级，对于F1来说分别将单侧鬼波G1的自相关和单侧鬼波G1与记录A的互相关带入Toeplitz矩阵，通过Levinson算法求解该矩阵得到匹配因子F1。对于F2来说，分别将单侧鬼波G2的自相关和单侧鬼波G2与记录A的互相关带入Toeplitz矩阵，通过Levinson算法求解该矩阵得到匹配因子F2。
鬼波衰减单元407，与所述记录A提取单元及能量匹配单元相连接，用于在所述记录A中减去所述匹配后的单侧鬼波G1m和双侧鬼波G2m，得到衰减鬼波后的记录Sdghost；具体计算方法如下：
Sdghost＝A-G1m-G2m；
其中，Sdghost是衰减鬼波后的记录。
上述实施例能够达到的有益效果：通过纯数据驱动型的鬼波预测和衰减，能够不依赖与速度场的精度，提高了海洋地震资料的处理质量。
本领域一般技术人员在此设计思想之下所做任何不具有创造性的改造，均应视为在本发明的保护范围之内。
Claims (10)
1.一种海洋地震资料处理中鬼波衰减的方法，其特征在于，包括：
将原始地震记录中的直达波和海面反射波切除，得到记录A；
在所述原始地震记录中提取出海面反射波M和反射层一次反射波P；
将所述海面反射波M和所述反射层一次反射波P褶积，得到该反射层相关的震源侧的单侧鬼波G1；
将所述海面反射波M与所述震源侧的单侧鬼波G1褶积，得到震源侧和检波器侧的双侧鬼波G2；
将所述单侧鬼波G1、双侧鬼波G2和记录A进行能量匹配，得到匹配后的单侧鬼波G1m和双侧鬼波G2m；
在所述记录A中减去所述匹配后的单侧鬼波G1m和双侧鬼波G2m，得到衰减鬼波后的记录Sdghost。
2.如权利要求1所述海洋地震资料处理中鬼波衰减的方法，其特征在于，所述将所述海面反射波M和所述反射层一次反射波P褶积，得到该反射层相关的震源侧的单侧鬼波G1，具体根据以下公式进行：
G1(xs,wr,f)＝M(xs,ws,f)·P(ws,wr,f)；
其中，M(xs,ws,f)是xs位置震源，ws位置接收到的海面反射波；
P(ws,wr,f)是ws位置震源，wr位置接收到的反射层一次反射波；
G1(xs,wr,f)是xs位置震源，wr位置接收到的单侧鬼波；
f代表频率域。
3.如权利要求2所述海洋地震资料处理中鬼波衰减的方法，其特征在于，所述将所述海面反射波M与所述震源侧的单侧鬼波G1褶积，得到震源侧和检波器侧的双侧鬼波G2具体根据以下公式进行：
G2(xs,xr,f)＝G1(xs,wr,f)·M(wr,xr,f)；
其中，M(wr,xr,f)是wr位置震源，xr位置接收到的海面反射波；
G2(xs,xr,f)为xs位置震源，xr位置接收到的震源侧和检波器侧的双侧鬼波。
4.如权利要求3所述海洋地震资料处理中鬼波衰减的方法，其特征在于，所述将所述单侧鬼波G1、双侧鬼波G2和记录A进行能量匹配，得到匹配后的单侧鬼波G1m和双侧鬼波G2m具体根据以下公式进行：
G1m＝F1·G1；
G2m＝F2·G2；
其中，F1与F2是基于最小二乘法在记录A中得到的匹配滤波因子。
5.如权利要求4所述海洋地震资料处理中鬼波衰减的方法，其特征在于，所述在所述记录A中减去所述匹配后的单侧鬼波G1m和双侧鬼波G2m，得到衰减鬼波后的记录Sdghost，根据以下公式进行：
Sdghost＝A-G1m-G2m；
其中，Sdghost是衰减鬼波后的记录。
6.一种海洋地震资料处理中鬼波衰减的装置，其特征在于，包括：
记录A提取单元，用于将原始地震记录中的直达波和海面反射波切除，得到记录A；
海面反射波提取单元，用于在所述原始地震记录中提取出海面反射波M；
一次反射波提取单元，用于在所述原始地震记录中提取出反射层一次反射波P；
单侧鬼波预测单元，与所述海面反射波提取单元及一次反射波提取单元相连接，用于将所述海面反射波M和所述反射层一次反射波P褶积，得到该反射层相关的震源侧的单侧鬼波G1；
双侧鬼波预测单元，与所述海面反射波提取单元及单侧鬼波预测单元相连接，用于将所述海面反射波M与所述震源侧的单侧鬼波G1褶积，得到震源侧和检波器侧的双侧鬼波G2；
能量匹配单元，与所述记录A提取单元、单侧鬼波预测单元及双侧鬼波预测单元相连接，用于将所述单侧鬼波G1、双侧鬼波G2和记录A进行能量匹配，得到匹配后的单侧鬼波G1m和双侧鬼波G2m；
鬼波衰减单元，与所述记录A提取单元及能量匹配单元相连接，用于在所述记录A中减去所述匹配后的单侧鬼波G1m和双侧鬼波G2m，得到衰减鬼波后的记录Sdghost。
7.如权利要求6所述海洋地震资料处理中鬼波衰减的装置，其特征在于，所述单侧鬼波预测单元根据以下公式进行计算得到单侧鬼波：
G1(xs,wr,f)＝M(xs,ws,f)·P(ws,wr,f)；
其中，M(xs,ws,f)是xs位置震源，ws位置接收到的海面反射波；
P(ws,wr,f)是ws位置震源，wr位置接收到的反射层一次反射波；
G1(xs,wr,f)是xs位置震源，wr位置接收到的单侧鬼波；
f代表频率域。
8.如权利要求7所述海洋地震资料处理中鬼波衰减的装置，其特征在于，所述双侧鬼波预测单元根据以下公式进行计算得到双侧鬼波：
G2(xs,xr,f)＝G1(xs,wr,f)·M(wr,xr,f)；
其中，M(wr,xr,f)是wr位置震源，xr位置接收到的海面反射波；
G2(xs,xr,f)为xs位置震源，xr位置接收到的震源侧和检波器侧的双侧鬼波。
9.如权利要求8所述海洋地震资料处理中鬼波衰减的装置，其特征在于，所述能量匹配单元根据以下公式进行能量匹配：
G1m＝F1·G1；
G2m＝F2·G2；
其中，F1与F2是基于最小二乘法在记录A中得到的匹配滤波因子。
10.如权利要求9所述海洋地震资料处理中鬼波衰减的装置，其特征在于，所述鬼波衰减单元根据以下公式进行鬼波衰减：
Sdghost＝A-G1m-G2m；
其中，Sdghost是衰减鬼波后的记录。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310187398.XA CN103308944B (zh) | 2013-05-20 | 2013-05-20 | 一种海洋地震资料处理中鬼波衰减的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310187398.XA CN103308944B (zh) | 2013-05-20 | 2013-05-20 | 一种海洋地震资料处理中鬼波衰减的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103308944A CN103308944A (zh) | 2013-09-18 |
| CN103308944B true CN103308944B (zh) | 2015-10-14 |
Family
ID=49134364
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310187398.XA Active CN103308944B (zh) | 2013-05-20 | 2013-05-20 | 一种海洋地震资料处理中鬼波衰减的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103308944B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104360394B (zh) * | 2014-12-01 | 2017-02-22 | 中国海洋石油总公司 | 一种海上变深度缆地震数据鬼波压制的方法 |
| CN104345343B (zh) * | 2014-12-01 | 2017-02-22 | 中国海洋石油总公司 | 一种复杂海底相关的层间多次波预测方法 |
| CN104536045A (zh) * | 2015-01-16 | 2015-04-22 | 中国海洋石油总公司 | 一种基于子波处理的鬼波压制方法 |
| CN106353817A (zh) * | 2016-08-11 | 2017-01-25 | 北京软岛科技有限公司 | 一种地震勘探数据多次波自动搜索和压制方法 |
| CN106526677B (zh) * | 2016-10-26 | 2018-12-21 | 中海石油(中国)有限公司 | 一种海上自适应压制鬼波的宽频逆时偏移成像方法 |
| CN107193045B (zh) * | 2017-07-07 | 2019-02-15 | 中国石油天然气集团公司 | 一种地震数据处理方法及装置 |
| CN108919357B (zh) * | 2018-05-16 | 2019-10-11 | 中国海洋石油集团有限公司 | 一种基于频谱重构的鬼波压制方法 |
| CN109975868B (zh) * | 2019-03-26 | 2020-03-03 | 中国科学院地质与地球物理研究所 | 一种基于泰勒展开的波动方程鬼波压制方法 |
| CN112578453A (zh) * | 2019-09-27 | 2021-03-30 | 中国石油化工股份有限公司 | 多次波压制方法及系统 |
| CN116774280A (zh) * | 2023-06-25 | 2023-09-19 | 中海石油(中国)有限公司深圳分公司 | 一种鬼波压制定量质控方法、装置、电子设备及存储介质 |
| CN117148443B (zh) * | 2023-10-27 | 2024-03-19 | 胜利信科(山东)勘察测绘有限公司 | 基于鬼波提取与转换的浅剖数据信噪比增强方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2405473B (en) * | 2003-08-23 | 2005-10-05 | Westerngeco Ltd | Multiple attenuation method |
| US9110191B2 (en) * | 2009-03-30 | 2015-08-18 | Westerngeco L.L.C. | Multiple attenuation for ocean-bottom seismic data |
| US8600680B2 (en) * | 2011-04-05 | 2013-12-03 | Pgs Geophysical As | Method for eliminating spectral constraints of acquisition system and earth filtering effects |
-
2013
- 2013-05-20 CN CN201310187398.XA patent/CN103308944B/zh active Active
Non-Patent Citations (4)
| Title |
|---|
| 《Deconvolving the ghost effect of the water surface in marine seismics》;Santi Kumar Ghosh;《GEOPHYSICS》;20001231;第65卷(第6期);第1831-1836页 * |
| 《Partial Surface Related Multiple Elimination》;Pierre Hugonnet et al.;《SEG International Exposition and 72nd Annual Meeting》;20021231;第1-4页 * |
| 《均衡拟多道匹配滤波法在波动方程法压制多次波中的应用》;李鹏等;《地球物理学报》;20071130;第50卷(第6期);第1844-1853页 * |
| 《多次波衰减技术》;郭全仕;《油气地球物理技术新进展—第71、72届SEG年会论文概要》;20021006;第205-215页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103308944A (zh) | 2013-09-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103308944B (zh) | 一种海洋地震资料处理中鬼波衰减的方法及装置 | |
| US10444393B2 (en) | Coherent noise attenuation using statistical methods | |
| EP2641109B1 (en) | Separating interfering signals in seismic data | |
| EP1660914B1 (en) | Method for attenuating water layer multiples in a seismic survey | |
| US9395456B2 (en) | System and method for narrow beam scanning microseismic monitoring | |
| CN103308943B (zh) | 一种海洋地震资料处理中层间多次波衰减的方法及装置 | |
| CN102985850B (zh) | 使用多能级传感器阵列的无源地震数据获取和处理 | |
| CN101937099B (zh) | 用于三维表面相关的多次波消除的动态孔径确定的方法 | |
| CN101620276A (zh) | 用于在地震数据中衰减多次反射的方法 | |
| CN104345343B (zh) | 一种复杂海底相关的层间多次波预测方法 | |
| CN104536045A (zh) | 一种基于子波处理的鬼波压制方法 | |
| WO2015126889A2 (en) | Methods and systems for using known source events in seismic data processing | |
| Kobayashi et al. | A field experiment of walkaway distributed acoustic sensing vertical seismic profile in a deep and deviated onshore well in Japan using a fibre optic cable deployed inside coiled tubing | |
| US10444387B2 (en) | Monitoring microseismic events | |
| CN102590863A (zh) | 确定反射系数频散、反射能量及识别气层的方法及装置 | |
| US8441890B2 (en) | Removing ground-roll from geophysical data using source-receiver interferometry | |
| EP2917857B1 (en) | Fracture characterization from refraction travel time data | |
| Ma et al. | Single borehole radar for well logging in a limestone formation: Experiments and simulations | |
| Cao et al. | High‐resolution and super stacking of time‐reversal mirrors in locating seismic sources | |
| Ringrose et al. | Developing optimised and cost-effective solutions for monitoring CO2 injection from subsea wells | |
| Grandi et al. | Efficient containment monitoring with distributed acoustic sensing: feasibility studies for the former Peterhead CCS Project | |
| EP3387466B1 (en) | Efficient internal multiple prediction methods | |
| CN105607128A (zh) | 一种软硬互层地质条件下的地质断层探测方法 | |
| Forgues et al. | Continuous high-resolution seismic monitoring of SAGD | |
| Melosh et al. | Seismic reflection data and conceptual models for geothermal development in Nevada |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |