CN104777015B - 一种海上地震资料中侧反射波的压制方法和装置 - Google Patents
一种海上地震资料中侧反射波的压制方法和装置 Download PDFInfo
- Publication number
- CN104777015B CN104777015B CN201410014290.5A CN201410014290A CN104777015B CN 104777015 B CN104777015 B CN 104777015B CN 201410014290 A CN201410014290 A CN 201410014290A CN 104777015 B CN104777015 B CN 104777015B
- Authority
- CN
- China
- Prior art keywords
- wave
- side reflection
- reflection
- source
- marine seismic
- 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
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000011410 subtraction method Methods 0.000 claims abstract description 13
- 239000013535 sea water Substances 0.000 claims description 37
- 238000001514 detection method Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000004088 simulation Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Abstract
本发明公开了一种海上地震资料中侧反射波的压制方法，属于海上地震资料处理领域。该方法包括：解析海上地震资料，得到海上地震总反射波信息；根据上述信息，计算侧反射源的地理坐标；并通过绕射波的波场方程模拟该地理坐标处的侧反射波的波场；在海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场。本发明通过计算侧反射源的地理坐标，并模拟该地理坐标处的侧反射波的波场，对侧反射源进行精确定位，真实全面地得到侧反射波信息。通过匹配相减法去除模拟的侧反射波的波场。该方法不受限于侧反射波与有效反射波的频谱形态，不仅能对侧反射波进行有效压制，还能使有效反射波进行保护并增强，提高海上地震资料的信噪比。
Description
技术领域
本发明涉及海上地震资料处理领域，特别涉及一种海上地震资料中侧反射波的压制方法和装置。
背景技术
海上地震资料是通过海上地震反射波勘探采集得到的、反映地下构造情况的海上地震数据，其对于海上开采地下油气资源有着重要的意义。在海上地震勘探过程中，炮源激发的反射波受海底地形高陡凸起、过往船只、岛礁等影响，会产生呈线性或双曲线性的侧反射波，侧反射波的存在会降低海上地震资料的信噪比，甚至会完全掩盖地下的有效反射波，其能量强，分布广，严重影响海上地震资料的质量。因此，有必要对海上地震资料中的侧反射波进行压制。
现有技术通过采用炮源-侧反射源-检波器探测系统来采集由侧反射波和有效反射波组成的海上地震总反射波。基于有效反射波和侧反射波在频域可以区分的前提条件下，采用分频的方式将侧反射波从海上地震总反射波中滤波出来，从而达到对侧反射波进行压制的目的。
在实现本发明的过程中，发明人发现现有技术至少存在以下问题：
在侧反射波与有效反射波的频谱形态几乎完全重叠，且侧反射波与有效反射波的能量相当的情况下，难以将侧反射波从海上地震总反射波中分离出来，无法有效地对侧反射波的进行压制。
发明内容
为了解决现有技术无法有效地对侧反射波的进行压制的问题，本发明实施例提供了一种海上地震资料中侧反射波的压制方法及装置。所述技术方案如下：
一方面，本发明实施例提供了一种海上地震资料中侧反射波的压制方法，所述方法为：
解析海上地震资料，得到海上地震总反射波信息，所述海上地震总反射波信息至少包括：所述侧反射波的旅行时、所述侧反射波在海水中的传播速度、所述侧反射波的入射角、偏移距、有效反射波在海水中的传播速度和炮源的地理坐标；
根据所述海上地震总反射波信息，计算所述海上地震资料中侧反射源的地理坐标；
根据所述海上地震总反射波的信息，通过绕射波的波场方程模拟所述地理坐标处的侧反射波的波场；
在所述海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场，得到处理后的海上地震反射波。
具体地，作为优选，所述根据所述海上地震总反射波信息，计算所述海上地震资料中侧反射源的地理坐标具体为包括：
根据所述侧反射波的旅行时、所述侧反射波在海水中的传播速度、所述侧反射波的入射角和所述偏移距，通过所述侧反射波的时距曲线方程确定所述炮源到所述侧反射源的距离；
根据所述炮源到所述侧反射源的距离和所述炮源的地理坐标，通过所述侧反射源地理坐标的求解方程，确定所述侧反射源的地理坐标。
具体地，所述侧反射波的时距曲线方程为：
其中，O'O为所述炮源到该侧反射源的距离、t为所述侧反射波的旅行时、x为所述偏移距、v为所述侧反射波在海水中的传播速度、α为所述侧反射波的入射角。
具体地，所述侧反射源地理坐标的求解方程为：
x、y分别为该侧反射源地理坐标的横、纵坐标点；O'O1、O'O2为不同炮源到该侧反射源的距离；P1x、P2x为不同炮源的地理坐标的横坐标；P1y、P2y为不同炮源的地理坐标的纵坐标。
具体地，所述绕射波的波场方程为：
其中，u1(rS,t)为所述侧反射波的波场、rS为所述侧反射源的地理坐标、u0(r,t)为有效反射波的波场、r为有效反射波的波场中任一点的地理坐标、t为侧反射波的旅行时、v1(r)为侧反射波在海水中的传播速度、v0(r)为有效反射波在海水中的传播速度、V为波场。
另一方面，本发明实施例还提供了一种海上地震资料中侧反射波的压制装置，所述装置包括：
解析模块，用于解析海上地震资料，得到海上地震总反射波信息，所述海上地震总反射波信息至少包括：所述侧反射波的旅行时、所述侧反射波在海水中的传播速度、所述侧反射波的入射角、偏移距、有效反射波在海水中的传播速度和炮源的地理坐标；
计算模块，用于根据所述海上地震总反射波信息，计算所述海上地震资料中侧反射源的地理坐标；
模拟模块，用于根据所述海上地震总反射波信息，通过绕射波的波场方程模拟所述地理坐标处的侧反射波的波场；
处理模块，用于在所述海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场，得到处理后的海上地震反射波。
具体地，所述计算模块用于根据所述侧反射波的旅行时、所述侧反射波在海水中的传播速度、所述侧反射波的入射角和该偏移距，通过所述侧反射波的时距曲线方程确定炮源到所述侧反射源的距离；根据所述炮源到所述该侧反射源的距离和所述炮源的地理坐标，通过所述侧反射源地理坐标的求解方程，确定所述侧反射源的地理坐标。
具体地，所述侧反射波的时距曲线方程为：
其中，O'O为所述炮源到所述侧反射源的距离、t为所述侧反射波的旅行时、x为所述偏移距、v为所述侧反射波在海水中的传播速度、α为所述侧反射波的入射角。
具体地，该侧反射源地理坐标的求解方程为：
其中，x、y分别为该侧反射源地理坐标的横、纵坐标点；O'O1、O'O2为不同炮源到该侧反射源的距离；P1x、P2x为不同炮源的地理坐标的横坐标；P1y、P2y为不同炮源的地理坐标的纵坐标。
具体地，所述绕射波的波场方程为：
其中，u1(rS,t)为所述侧反射波的波场、rS为所述侧反射源的地理坐标、u0(r,t)为有效反射波的波场、r为有效反射波的波场中任一点的地理坐标、t为侧反射波的旅行时、v1(r)为侧反射波在海水中的传播速度、v0(r)为有效反射波在海水中的传播速度、V为波场。
本发明实施例提供的技术方案带来的有益效果是：
本发明实施例提供的海上地震资料中侧反射波的压制方法及装置，通过海上地震总反射波信息，计算得到侧反射源的地理坐标，并根据绕射波的波场方程模拟所述地理坐标处的侧反射波的波场，对侧反射源进行精确定位，真实全面地得到侧反射波信息。然后在所述海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场，有效压制了侧反射波。本发明实施例所述海上地震资料中侧反射波的压制方法及装置，不受限于侧反射波与有效反射波的频谱形态，不仅能够对侧反射波进行有效压制，而且能够对有效反射波进行保护并增强，使得海上地震资料的信噪比得到提高。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的海上地震资料中侧反射波的压制方法流程图；
图2是本发明又一实施例提供的海上地震资料中侧反射波的压制方法流程图；
图3是炮源-侧反射源-检波器探测系统模拟示意图；
图4a是本发明又一实施例提供的侧反射波压制前的海上地震总反射波的示意图；
图4b是本发明又一实施例提供的对侧反射波压制后的海上地震总反射波的示意图；
图4c是本发明又一实施例提供的海上地震总反射波中侧反射波的示意图。
其中，O表示炮源，
O′表示侧反射源，
B、D均表示检波器中接收点,
1表示海上地震总反射波，
11表示有效反射波，
12表示侧反射波。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
实施例1
如附图1所示，本发明实施例提供了一种海上地震资料中侧反射波的压制方法，该方法为：
101、解析海上地震资料，得到海上地震总反射波信息，该海上地震总反射波信息至少包括：该侧反射波的旅行时、该侧反射波在海水中的传播速度、该侧反射波的入射角、偏移距、有效反射波在海水中的传播速度和炮源的地理坐标。
102、根据该海上地震总反射波信息，计算该海上地震资料中侧反射源的地理坐标。
103、根据该海上地震总反射波的信息，通过绕射波的波场方程模拟该地理坐标处的侧反射波的波场。
104、在该海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场，得到处理后的海上地震反射波。
本发明实施例提供的海上地震资料中侧反射波的压制方法，通过海上地震总反射波信息，计算得到侧反射源的地理坐标，并根据绕射波的波场方程模拟该地理坐标处的侧反射波的波场，对侧反射源进行精确定位，真实全面地得到侧反射波信息。然后在该海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场，有效压制了侧反射波。本发明实施例该海上地震资料中侧反射波的压制方法，不受限于侧反射波与有效反射波的频谱形态，不仅能够对侧反射波进行有效压制，而且能够对有效反射波进行保护并增强，使得海上地震资料的信噪比得到提高。
实施例2
图2是本发明实施例提供的一种海上地震资料中侧反射波的压制方法的流程图。参见图2，该方法为：
201、解析海上地震资料，得到海上地震总反射波信息，该海上地震总反射波信息至少包括：该侧反射波的旅行时、该侧反射波在海水中的传播速度、该侧反射波的入射角、偏移距、有效反射波在海水中的传播速度和炮源的地理坐标。
其中，海上地震资料为实际采集到的实际地震道集数据，该实际地震道集数据包括由侧反射波和有效反射波组成的海上地震总反射波信息。侧反射波的旅行时为侧反射波由炮源经侧反射源到检波点所用的时间。侧反射波的入射角为侧反射波在侧反射源处反射方向与其法线方向之间的夹角。偏移距为炮源到检波点的距离。
本发明实施例通过使用炮源-侧反射源-检波器探测系统来采集海上地震资料，并对海上地震资料进行解析得到海上地震总反射波的信息，其具体过程如下：通过在海上炮源处人工激发地震波，该地震波在海水中传播，并穿过海上侧反射源进行反射后由检波器中的接收点（以下简称：检波点）进行接收，得到海上地震资料。对上述接收到的海上地震资料进行解析，得到海上地震总反射波信息，为后续模拟侧反射波的波场提供真实的数据支持。
由于采集海上地震总反射波信息时，不会受限于侧反射波与有效反射波的频谱形态，所以即使在侧反射波与有效反射波的频谱形态几乎完全重叠、且两者能量相当的情况下，也能够准确而全面地得到海上地震总反射波信息。可见，步骤201通过解析海上地震资料，得到海上地震总反射波信息，具有更好的保真性。
202、根据该侧反射波的旅行时、该侧反射波在海水中的传播速度、该侧反射波的入射角和该偏移距，通过该侧反射波的时距曲线方程确定该炮源到该侧反射源的距离。
可选地，该时距曲线方程为：
其中，O'O为该炮源到该侧反射源的距离、t为该侧反射波的旅行时、x为该偏移距、v为该侧反射波在海水中的传播速度、α为该侧反射波的入射角。
通过侧反射波的时距曲线方程能够确定侧反射波在海上地震资料中的分布范围，即可确定侧反射源到震源的距离，为下一步计算侧反射源的地理坐标奠定基础。
该侧反射波的时距曲线方程通过如下步骤确定：如附图3所示的炮源-侧反射源-检波器探测系统模拟示意图中，O表示炮源，O′表示侧反射源，B、D均表示检波器中接收点。O'B及O'D均为侧反射源到检波点的距离，两者均可根据所采集的炮源到检波点的距离参数来确定。其中D点满足侧反射源到检波点的距离最短。tD及tB均为侧反射波的旅行时，且tD及tB分别通过公式（1）-（2）来计算得到：
侧反射波的入射角α通过公式（3）来确定：
根据余弦定理，基于公式（1）-（3），即可推导得到侧反射波的时距曲线方程。
203、根据该炮源到该侧反射源的距离和该炮源的地理坐标，通过该侧反射源地理坐标的求解方程，确定该侧反射源的地理坐标。
该侧反射源地理坐标的求解方程为：
x、y分别为该侧反射源地理坐标的横、纵坐标点；O'O1、O'O2为不同炮源到该侧反射源的距离；P1x、P2x为不同炮源的地理坐标的横坐标；P1y、P2y为不同炮源的地理坐标的纵坐标。
步骤203中，确定该侧反射源的地理坐标是为了能够对侧反射源进行精确定位，避免侧反射源的遗漏，以便真实、全面地对侧反射源处发出的侧反射波的波场进行模拟。
优选地，该侧反射源的数量至少为一个。所以步骤203可计算得到至少一个反射源的地理坐标，从而更加准确地模拟侧反射波的波场。
204、根据该海上地震总反射波的信息，通过绕射波的波场方程模拟该地理坐标处的侧反射波的波场。
该侧反射波的波场方程为：
其中，u1(rS,t)为该侧反射波的波场、rS为该侧反射源的地理坐标、u0(r,t)为有效反射波的波场、r为有效反射波的波场中任一点的地理坐标、t为侧反射波的旅行时、v1(r)为侧反射波在海水中的传播速度、v0(r)为有效反射波在海水中的传播速度、V为波场。
地震波的波场是指有地震波传播的空间。在该空间的每一点上，一定时刻都有一定的地震波通过，且该地震波的能量也按—定的规律传播。可见，根据地震波的波场就能够确定在空间内传播的地震波的分布。所以，步骤204通过绕射波的波场方程模拟得到该地理坐标处的侧反射波的波场，即可确定该侧反射波在海上地震资料中的分布。
由于侧反射波为一种绕射波，为了真实地模拟和保留上述海上地震资料中的侧反射波，本发明实施例以绕射理论为基础，利用绕射波波场的运动学和动力学特征，来模拟该侧反射波的波场。
根据绕射理论，在炮源-侧反射源-检波器探测系统中，入射声波在某种参考介质中的波场称为参考波场（该参考波场可理解为本发明实施例该有效反射波的波场）。如果参考波场中的参考介质发生变化或者扰动，此时观测到的波场称为总波场（该参考波场可理解为本发明实施例该海上地震总反射波的波场）。该总波场与参考波场之差称为散射波场（该散射波场可理解为本发明实施例该侧反射波的波场）。所以，在声波绕射过程中，将一定能量的声源放在参考介质C0(x)中，可接收到参考波场u0(x,t)；把同样能量的声源放在发生变化或者扰动的参考介质C(x)中，可接收到总波场为u(x,t)；则散射波场为u1(x,t)=u(x,t)-u0(x,t)。其中，x为波场内任一点的地理坐标，单位为米；t为波场内声波的传播时间，单位为秒。
设从侧反射源rs(x,y,z)处发出的海上地震总反射波d(t)被检波点rG(x,y,z)处的检波器接收，此时该海上地震总反射波的三维波动方程如公式（4）-（5）所示，在时间域内表示为：
其中r为波场空间一点，其地理坐标可表示为(x，y，z)；v(r)为海上地震总反射波的传播速度，单位为米/秒。
根据绕射理论，实际的海上地震总反射波的传播速度v(r)可以表示成有效反射波的传播速度v0(r)和侧反射波的传播速度ν1(r)的叠加，具体叠加方程如公式（6）所示：
基于公式（6）与公式（4）及公式（5），得到海上地震总反射波的波动方程，如公式（7）所示：
对于微扰动系统有ν1(r)→0，u1→0，所以：
ν1(r)u1(r,t)<<ν1(r)u0(r,t) (8)
最后根据Born近似（Born Approximation），即可得到侧反射波的波场方程。由侧反射波的波场方程可以看出侧反射波场与侧反射波海水中的传播速度、有效反射波在海水中的传播速度以及有效反射波的波场直接相关，所以根据炮源-侧反射源-检波器探测系统所采集到的海上地震总反射波信息，能够准确地模拟侧反射波的波场。
205、在该海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场，得到处理后的海上地震反射波。
对于步骤205，可将采集到的实际海上地震总反射波的波场和模拟的侧反射波的波场转换为二维的频率波或时间波，然后将转换后的频率波或时间波进行匹配相减即可达到压制侧反射波的目的。由于模拟的侧反射波的波场与实际采集到的海上地震总反射波中所记录的侧反射波的波场不可能完全一致（例如，它们在相位、振幅以及波的到达时间上都有可能存在一定的差异）。因此，本发明实施例采用匹配相减法，将模拟的侧反射波的波场从实际海上地震总反射波的波场中去除，能够更加精确地保证侧反射波得到有效压制。
如附图4a所示，使用本发明实施例提供的方法对含有侧反射波的海上地震总反射波1进行处理前，在海上地震总反射波1的有效信息区域能清晰地分辨出呈倒抛物线性的侧反射波12。通过采用本发明实施例提供的方法对该侧反射波进行模拟，并在实际海上地震总反射波上匹配相减后，能够清晰地分辨出侧反射波得到压制后的有效反射波11（如附图4b所示），以及从海上地震总反射波1上去除的侧反射波12（如附图4c所示）。
可见，通过采用本发明实施例提供的海上地震资料中侧反射波的压制方法，侧反射波得到了有效地压制，并使海上地震资料达到了保真处理和提高信噪比的目的。由于海上地震资料中侧反射波得到了有效压制，使得其成像和地质解释更加精确，显著提高了海上地震资料的处理效率，对于海上开采地下油气资源具有重要的贡献。
实施例3
本发明实施例提供了一种海上地震资料中侧反射波的压制装置，该装置包括：
解析模块，用于解析海上地震资料，得到海上地震总反射波信息，该海上地震总反射波信息至少包括：该侧反射波的旅行时、该侧反射波在海水中的传播速度、该侧反射波的入射角、偏移距、有效反射波在海水中的传播速度和炮源的地理坐标。
计算模块，用于根据该海上地震总反射波信息，计算该海上地震资料中侧反射源的地理坐标。
模拟模块，用于根据该海上地震总反射波信息，通过绕射波的波场方程模拟该地理坐标处的侧反射波的波场。
处理模块，用于在该海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场，得到处理后的海上地震反射波。
本发明实施例提供的海上地震资料中侧反射波的压制装置，通过海上地震总反射波信息，计算得到侧反射源的地理坐标，并根据绕射波的波场方程模拟所述地理坐标处的侧反射波的波场，对侧反射源进行精确定位，真实全面地得到侧反射波信息。然后在所述海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场，有效压制了侧反射波。本发明实施例中该海上地震资料中侧反射波的压制装置，不受限于侧反射波与有效反射波的频谱形态，不仅能够对侧反射波进行有效压制，而且能够对有效反射波进行保护并增强，使得海上地震资料的信噪比得到提高。
实施例4
本发明实施例提供了一种海上地震资料中侧反射波的压制装置，该装置包括：
解析模块，用于解析海上地震资料，得到海上地震总反射波信息，该海上地震总反射波信息至少包括：该侧反射波的旅行时、该侧反射波在海水中的传播速度、该侧反射波的入射角、偏移距、有效反射波在海水中的传播速度和炮源的地理坐标。
计算模块，用于根据该侧反射波的旅行时、该侧反射波在海水中的传播速度、该侧反射波的入射角和该偏移距，通过该侧反射波的时距曲线方程确定炮源到所述侧反射源的距离；根据该炮源到该侧反射源的距离和该炮源的地理坐标，通过该侧反射源地理坐标的求解方程，确定该侧反射源的地理坐标。
其中，该侧反射波的时距曲线方程为：
其中，O'O为该炮源到所述侧反射源的距离、t为该侧反射波的旅行时、x为该偏移距、v为该侧反射波在海水中的传播速度、α为该侧反射波的入射角。
该侧反射源地理坐标的求解方程为：
其中，x、y分别为该侧反射源地理坐标的横、纵坐标点；O'O1、O'O2为不同炮源到该侧反射源的距离；P1x、P2x为不同炮源的地理坐标的横坐标；P1y、P2y为不同炮源的地理坐标的纵坐标。
模拟模块，用于根据该海上地震总反射波信息，通过绕射波的波场方程模拟该地理坐标处的侧反射波的波场。
其中，该绕射波的波场方程为：
其中，u1(rS,t)为该侧反射波的波场、rS为该侧反射源的地理坐标、u0(r,t)为有效反射波的波场、r为有效反射波的波场中任一点的地理坐标、t为侧反射波的旅行时、v1(r)为侧反射波在海水中的传播速度、v0(r)为有效反射波在海水中的传播速度、V为波场。
处理模块，用于在该海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场，得到处理后的海上地震反射波。
本发明实施例提供的海上地震资料中侧反射波的压制装置，侧反射波得到了有效地压制，并使海上地震资料达到了保真处理和提高信噪比的目的。由于海上地震资料中侧反射波得到了有效压制，使得其成像和地质解释更加精确，显著提高了海上地震资料的处理效率，对于海上开采地下油气资源具有重要的贡献
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (2)
1.一种海上地震资料中侧反射波的压制方法，其特征在于，所述方法为：
解析海上地震资料，得到海上地震总反射波信息，所述海上地震总反射波信息至少包括：所述侧反射波的旅行时、所述侧反射波在海水中的传播速度、所述侧反射波的入射角、偏移距、有效反射波在海水中的传播速度和炮源的地理坐标；
根据所述海上地震总反射波信息，计算所述海上地震资料中侧反射源的地理坐标；
根据所述海上地震总反射波信息，通过绕射波的波场方程模拟所述地理坐标处的侧反射波的波场；
在所述海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场，得到处理后的海上地震反射波；
所述海上地震资料通过使用炮源-侧反射源-检波器探测系统采集得到；
所述绕射波的波场方程为：
<mrow>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>S</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mi>V</mi>
<mn>0</mn>
</msubsup>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>v</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<mi>u</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>V</mi>
</mrow>
其中，u1(rS,t)为所述侧反射波的波场、rS为所述侧反射源的地理坐标、u0(r,t)为有效反射波的波场、r为有效反射波的波场中任一点的地理坐标、t为侧反射波的旅行时、v1(r)为侧反射波在海水中的传播速度、v0(r)为有效反射波在海水中的传播速度、V为波场；
所述根据所述海上地震总反射波信息，计算所述海上地震资料中侧反射源的地理坐标包括：
根据所述侧反射波的旅行时、所述侧反射波在海水中的传播速度、所述侧反射波的入射角和偏移距，通过侧反射波的时距曲线方程确定炮源到所述侧反射源的距离；
根据所述炮源到所述该侧反射源的距离和所述炮源的地理坐标，通过侧反射源地理坐标的求解方程，确定所述侧反射源的地理坐标；
所述侧反射波的时距曲线方程为：
<mrow>
<mi>t</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<mi>O</mi>
</mrow>
<mi>v</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<mi>O</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<mi>O</mi>
<mo>&times;</mo>
<mi>x</mi>
<mo>&times;</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mi>v</mi>
</mfrac>
</mrow>
其中，O'O为所述炮源到所述侧反射源的距离、t为所述侧反射波的旅行时、x为所述偏移距、v为所述侧反射波在海水中的传播速度、α为所述侧反射波的入射角；
所述侧反射源地理坐标的求解方程为：
<mrow>
<msubsup>
<mi>P</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mi>x</mi>
<mo>+</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mi>x</mi>
<mo>&times;</mo>
<mi>x</mi>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mi>y</mi>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<mi>y</mi>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mi>y</mi>
<mo>=</mo>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<msubsup>
<mi>O</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>P</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mi>x</mi>
<mo>+</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<mi>x</mi>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mi>x</mi>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mi>y</mi>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<mi>y</mi>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mi>y</mi>
<mo>=</mo>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<msubsup>
<mi>O</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
其中，x、y分别为该侧反射源地理坐标的横、纵坐标点；O'O1、O'O2为不同炮源到该侧反射源的距离；P1x、P2x为不同炮源的地理坐标的横坐标；P1y、P2y为不同炮源的地理坐标的纵坐标。
2.一种用于权利要求1所述海上地震资料中侧反射波的压制方法的海上地震资料中侧反射波的压制装置，其特征在于，所述装置包括：
解析模块，用于解析海上地震资料，得到海上地震总反射波信息，所述海上地震总反射波信息至少包括：所述侧反射波的旅行时、所述侧反射波在海水中的传播速度、所述侧反射波的入射角、偏移距、有效反射波在海水中的传播速度和炮源的地理坐标；
计算模块，用于根据所述海上地震总反射波信息，计算所述海上地震资料中侧反射源的地理坐标；
模拟模块，用于根据所述海上地震总反射波信息，通过绕射波的波场方程模拟所述地理坐标处的侧反射波的波场；
处理模块，用于在所述海上地震总反射波的波场上，运用匹配相减法去除模拟的侧反射波的波场，得到处理后的海上地震反射波；
所述海上地震资料通过使用炮源-侧反射源-检波器探测系统采集得到；
所述绕射波的波场方程为：
<mrow>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>S</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mi>V</mi>
<mn>0</mn>
</msubsup>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>v</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<mi>u</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>V</mi>
</mrow>
其中，u1(rS,t)为所述侧反射波的波场、rS为所述侧反射源的地理坐标、u0(r,t)为有效反射波的波场、r为有效反射波的波场中任一点的地理坐标、t为侧反射波的旅行时、v1(r)为侧反射波在海水中的传播速度、v0(r)为有效反射波在海水中的传播速度、V为波场；
所述计算模块用于根据所述侧反射波的旅行时、所述侧反射波在海水中的传播速度、所述侧反射波的入射角和该偏移距，通过侧反射波的时距曲线方程确定炮源到所述侧反射源的距离；根据所述炮源到所述该侧反射源的距离和所述炮源的地理坐标，通过侧反射源地理坐标的求解方程，确定所述侧反射源的地理坐标；
所述侧反射波的时距曲线方程为：
<mrow>
<mi>t</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<mi>O</mi>
</mrow>
<mi>v</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<mi>O</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<mi>O</mi>
<mo>&times;</mo>
<mi>x</mi>
<mo>&times;</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mi>v</mi>
</mfrac>
</mrow>
其中，O'O为所述炮源到所述侧反射源的距离、t为所述侧反射波的旅行时、x为所述偏移距、v为所述侧反射波在海水中的传播速度、α为所述侧反射波的入射角；
侧反射源地理坐标的求解方程为：
<mrow>
<msubsup>
<mi>P</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mi>x</mi>
<mo>+</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mi>x</mi>
<mo>&times;</mo>
<mi>x</mi>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mi>y</mi>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<mi>y</mi>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mi>y</mi>
<mo>=</mo>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<msubsup>
<mi>O</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>P</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mi>x</mi>
<mo>+</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<mi>x</mi>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mi>x</mi>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mi>y</mi>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<mi>y</mi>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mi>y</mi>
<mo>=</mo>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<msubsup>
<mi>O</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
其中，x、y分别为该侧反射源地理坐标的横、纵坐标点；O'O1、O'O2为不同炮源到该侧反射源的距离；P1x、P2x为不同炮源的地理坐标的横坐标；P1y、P2y为不同炮源的地理坐标的纵坐标。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410014290.5A CN104777015B (zh) | 2014-01-13 | 一种海上地震资料中侧反射波的压制方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410014290.5A CN104777015B (zh) | 2014-01-13 | 一种海上地震资料中侧反射波的压制方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104777015A CN104777015A (zh) | 2015-07-15 |
| CN104777015B true CN104777015B (zh) | 2018-06-01 |
Family
ID=
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102890291A (zh) * | 2012-10-09 | 2013-01-23 | 中国海洋石油总公司 | 一种地震资料中外源干扰的自适应压制方法 |
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102890291A (zh) * | 2012-10-09 | 2013-01-23 | 中国海洋石油总公司 | 一种地震资料中外源干扰的自适应压制方法 |
Non-Patent Citations (2)
| Title |
|---|
| 海上地震数据中侧面干扰波的特征与压制方法;蔡希玲等;《中国海上油气》;20050425(第02期) * |
| 海上船噪声干扰的压制技术;张雅勤等;《物探与化探》;20010820(第04期) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3227725B1 (en) | Seismic acquisition method | |
| CN102183435B (zh) | 一种基于多路径反射理论的海底密度和声速测量方法 | |
| US10466378B2 (en) | Impact assessment of marine seismic surveys | |
| US20020118602A1 (en) | Angle dependent surface multiple attenuation for two-component marine bottom sensor data | |
| US10073183B2 (en) | Methods and systems that attenuate noise in seismic data | |
| Hui et al. | Underwater acoustic channel | |
| Kargl et al. | Synthetic aperture sonar imaging of simple finite targets | |
| Marsset et al. | Deep-towed high resolution seismic imaging II: Determination of P-wave velocity distribution | |
| Crone et al. | Estimating shallow water sound power levels and mitigation radii for the R/VM arcus G. L angseth using an 8 km long MCS streamer | |
| CN106610507B (zh) | 一种去除虚反射干扰的方法和装置 | |
| Matthews et al. | Comparing modeled and measured sound levels from a seismic survey in the Canadian Beaufort Sea | |
| US9791580B2 (en) | Methods and systems to separate wavefields using pressure wavefield data | |
| JP2008014830A (ja) | ハイドレートの存在領域探査方法及び探査システム | |
| CN104777015B (zh) | 一种海上地震资料中侧反射波的压制方法和装置 | |
| US20120195168A1 (en) | Method for detection of a plume of a first fluid within a second fluid | |
| CN115236592B (zh) | 一种基于单阵元时频曲线匹配的冰声定位方法 | |
| CN114442152B (zh) | 地震数据的海底多次波预测方法及装置 | |
| Skinner et al. | Development of underwater beacon for Arctic through-ice communication via satellite | |
| Hines et al. | Time-of-flight measurements of acoustic wave speed in a sandy sediment at 0.6–20 kHz | |
| GB2375606A (en) | Angle dependent surface multiple attenuation for two - component marine bottom sensor data | |
| Culver et al. | Sonar signal processing using probabilistic signal and ocean environmental models | |
| RU2292569C1 (ru) | Способ определения предвестника цунами | |
| RU2300781C1 (ru) | Устройство гидрометеорологоакустических наблюдений за акваторией морского полигона | |
| Caley | Development of a dynamic underwater acoustic communication channel simulator with configurable sea surface parameters to explore time-varying signal distortion | |
| WO2016011581A1 (zh) | 一种海上地震资料侧反射压制方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |