一种深拖式高分辨率多道地震数据精细处理方法与流程

1.本发明涉及海洋地球物理勘探技术领域，具体涉及一种深拖式高分辨率多道地震数据精细处理方法，应用于深拖式高分辨率多道地震探测系统的资料处理环节，以获得高品质的深拖式高分辨率多道地震成像剖面。


背景技术：

2.海面拖曳的探测方式受大深度水体的影响，目标地层菲涅耳带非常大且震源高频能量大幅度衰减，很难实现真正意义上的高分辨率。现有技术中，采用将高主频深拖震源（主频约750 hz）和多道检波点阵列（48道）同时拖曳于距离海底约100 m的深拖式高分辨率多道地震探测方式，可极大地缩小菲涅耳带，减小海水吸收、海洋混响和避免多次波的影响，能够获得非常高的横向和纵向分辨率地震资料。然而，受大深度水体与拖船速度的影响，实现稳定、安全、高质量、高效的数据采集非常困难，深拖震源和多道检波点阵列在数据采集过程中是上下浮动的，其拖曳轨迹与常规海面拖曳式多道地震勘探方法也是不同的。在深拖式高分辨率多道地震数据的数据处理方面，由于深拖震源和多道检波点阵列无法使用gps实现高精度定位，数据采集过程中深拖震源-检波点相对位置关系时刻在变化，且拖曳轨迹与常规二维高分辨率多道地震探测方法具有差异性，常规海面拖曳式多道地震数据处理是以海平面为基准面的，其相应的数据处理技术与方法无法直接应用于深拖式高分辨率多道地震数据的处理。


技术实现要素：

3.针对现有技术存在的问题，本发明提供一种深拖式高分辨率多道地震数据精细处理方法，其能够解决深拖式高分辨率多道地震探测系统的拖曳轨迹和深拖震源-检波点相对位置关系时刻变化带来的数据处理问题，提高深拖式高分辨率多道地震成像剖面的信噪比和分辨率，使其能够提供更为详实的海底沉积和构造细节。
4.本发明采取以下技术方案：一种深拖式高分辨率多道地震数据精细处理方法，所述方法具体包括：步骤s1：将深拖震源和各道检波点置于一个起伏较为平缓的浮动基准面上，获取所述多道地震数据浮动基准面校正量，对cmp道集进行浮动基准面校正；步骤s2：对cmp道集进行动校正，获取动校正后cmp道集内剩余时差并对cmp道集进行剩余时差校正；步骤s3：将cmp道集内剩余时差校正后的海底反射波走时与其对应的自激自收时间计算自激自收时差，进行自激自收时差校正；步骤s4：利用相位替代法对自激自收时差校正后的cmp道集进行非一致性时差校正；步骤s5：对非一致性时差校正后的叠加剖面进行海底平滑时差校正。
5.进一步地，所述步骤s1具体为：
步骤s1.1：计算基于固定基准面的常规静校正量，计算公式如下：其中，表示第j个炮点或检波点对应的常规静校正量，表示第j个炮点或检波点的离底高度，表示海水声速，为固定基准面的高程，为第j个炮点或检波点的入水深度，为替换速度；步骤s1.2：计算浮动基准面和固定基准面之间的静校正量，计算公式如下：其中，表示第k个cmp道集对应的静校正量，为固定基准面的高程，为第k个cmp道集的浮动基准面高程，为替换速度；步骤s1.3：将式（1）和式（2）计算得到的静校正量应用到cmp道集内从而完成浮动基准面校正。
6.进一步地，所述步骤s2具体为：对cmp道集进行动校正，拾取动校正后cmp道集内海底反射波走时，统计cmp道集内其他道相对某道的时差，并选取时差总和最小的道为参考道，计算其他道与参考道的海底反射波走时时差，在cmp道集内应用所述海底反射波走时时差进行剩余时差校正。
7.进一步地，所述步骤s3具体为：利用浮动基准面高程计算自激自收时间，并将cmp道集内剩余时差校正后的海底反射波走时与其对应的自激自收时间计算自激自收时差，进行自激自收时差校正。
8.进一步地，所述自激自收时间的计算公式如下：其中，表示利用浮动基准面高程计算的第k个cmp道集的自激自收时间，为第k个cmp道集的浮动基准面高程，为海水平均声速。
9.进一步地，所述步骤s4具体为：对cmp道集中的每一道使用快速傅里叶变换，求其相位谱和振幅谱，然后利用模型道的相位谱替换cmp道集中每一道的相位谱后再进行反傅里叶变换。
10.进一步地，所述模型道为做过自激自收时差校正的叠加道。
11.进一步地，所述步骤s4中进行海底平滑时差校正的海底平滑时差校正量的计算公式如下：其中，表示进行海底平滑时差校正的海底平滑时差校正量，表示利用浮动基准面高程计算的第k个cmp道集的自激自收时间，为第k个cmp位置对应海底位置的海底深度，为海水平均声速。
12.进一步地，和来自深拖式高分辨率多道地震接收阵列形态重建反演结果。
13.由于采取以上技术方案，本发明与现有技术相比实现了以下有益效果：（1）在深拖式高分辨率多道地震探测系统工作过程中，首先深拖震源和检波点是连续移动的，且深拖震源的位置间隔取决于船的速度变化；其次深拖震源和检波点的深度
是实时变化的；基于地表一致性假设或常规非地表一致性剩余静校正方法均不适用于深拖式高分辨率多道地震探测的高频数据。本发明通过综合引入浮动基准面校正、cmp道集内剩余时差校正、自激自收时差校正、利用相位替代法的非一致性时差校正及海底平滑时差校正等，充分消除各类时差，实现海底以下地层构造形态的高分辨率、高信噪比、高保真度成像。
14.（2）本发明通过浮动基准面校正聚焦了cmp道集的速度谱能量团，有利于精细速度分析。
15.（3）本发明通过剩余时差校正消除了cmp道集剩余时差，提高数据处理精度。
16.（4）本发明通过自激自收时差校正提高了同相轴连续性。
17.（5）本发明通过非一致性时差校正消除了cmp道集的非一致性时差，提高了叠加剖面的信噪比和分辨率。
附图说明
18.图1为本发明实施例浮动基准面设置示意图；图2为本发明实施例浮动基准面校正前后的速度谱对比图；图3为本发明实施例剩余时差校正前后的cmp道集波形对比图；图4为本发明实施例自激自收时差校正前后的叠加剖面对比图；图5为本发明实施例非一致性时差校正前后的cmp道集波形对比图；图6为本发明实施例非一致性时差校正前后的叠加剖面对比图；图7为本发明实施例海底平滑时差校正前后的叠加剖面对比图。
具体实施方式
19.为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
20.本发明提供一种深拖式高分辨率多道地震数据精细处理方法，在对深拖式高分辨率多道地震数据进行精细处理前，首先需要通过深拖式高分辨率多道地震探测系统采集得到深拖式高分辨率多道地震数据，所述方法包括以下步骤：步骤s1：将深拖震源和各道检波点置于一个起伏较为平缓的浮动基准面上，获取所述多道地震数据浮动基准面校正量，对cmp道集进行浮动基准面校正。
21.本实施例中，深拖式高分辨率多道地震探测系统的浮动基准面设置情况如图1所示，其中，c为固定基准面，s为深拖震源拖曳轨迹，f为浮动基准面，f1、f2、f3、f4为炮点，浮动基准面高程采用平滑炮点拖曳轨迹来确定，平滑半径为排列长度158 m。具体到本实施例中，在炮点f1、f2、f3、f4处，以排列长度158 m为平滑半径，获得浮动基准面。
22.所述步骤s1具体为：步骤s1.1：计算基于固定基准面的常规静校正量，计算公式如下：其中，表示第j个炮点或检波点对应的常规静校正量，表示第j个炮点或检波
点的离底高度，表示海水声速，为固定基准面的高程，为第j个炮点或检波点的入水深度，为替换速度。在本实施例中和来自深拖式高分辨率多道地震接收阵列形态重建反演结果，海水声速取1481.9 m/s，替换速度取海底以下地层速度的平均值1610 m/s。由式（1）计算得到炮点和检波点相对固定基准面具有地表一致性特点的静校正量。
23.步骤s1.2：计算浮动基准面和固定基准面之间的静校正量，计算公式如下：其中，表示第k个cmp道集对应的静校正量，为固定基准面的高程，为第k个cmp道集的浮动基准面高程，为替换速度。在本实施例中，替换速度取海底以下地层速度的平均值1610 m/s。
24.步骤s1.3：将式（1）和式（2）计算得到的静校正量应用到cmp道集内从而完成浮动基准面校正。
25.图2左图为浮动基准面校正前的速度谱，图2右图为浮动基准面校正后的速度谱，经对比可以看出，经过浮动基准面校正后，速度谱能量团聚焦效果非常好，有利于精细速度分析。
26.步骤s2：对cmp道集进行动校正，获取动校正后cmp道集内剩余时差并对cmp道集进行剩余时差校正。
27.所述步骤s2具体为：对cmp道集进行动校正，拾取动校正后cmp道集内海底反射波走时，统计cmp道集内其他道相对某道的时差，并选取时差总和最小的道为参考道，计算其他道的海底反射波走时与参考道的海底反射波走时之差，得到其他道与参考道的海底反射波走时时差，在cmp道集内应用所述海底反射波走时时差进行剩余时差校正。
28.步骤s3：将cmp道集内剩余时差校正后的海底反射波走时与其对应的自激自收时间计算自激自收时差，进行自激自收时差校正。
29.所述步骤s3具体为：利用浮动基准面高程计算自激自收时间，并将cmp道集内剩余时差校正后的海底反射波走时与其对应的自激自收时间计算自激自收时差，进行自激自收时差校正，避免了同相轴不连续的问题。
30.进一步的，所述自激自收时间的计算公式如下：其中，表示利用浮动基准面高程计算的第k个cmp道集的自激自收时间，为第k个cmp道集的浮动基准面高程，为海水平均声速。在本实施例中，海水平均声速取1481.9 m/s。
31.利用上述步骤对本实施例经过浮动基准面校正后的cmp道集进行处理，图3左图为经过动校正后的cmp道集，可以发现存在一定的时差，图3右图为经过剩余时差校正后的cmp道集，时差基本上得到了消除；图4左图为自激自收时差校正前的叠加剖面效果图，图4右图为自激自收时差校正后的叠加剖面效果图，经对比可以发现在同相轴出现不连续的位置，其时差基本上得到了消除，同相轴连续性明显变好。
32.步骤s4：利用相位替代法对自激自收时差校正后的cmp道集进行非一致性时差校正。
33.本实施例在经过上述自激自收时差校正后的cmp道集，利用相位替代法校正非一致性时差，非一致性时差是指在不同层位存在的不同走时残差。
34.所述步骤s4具体为：对cmp道集中的每一道使用快速傅里叶变换，求其相位谱和振幅谱，然后利用模型道的相位谱替换cmp道集中每一道的相位谱后再进行反傅里叶变换。
35.优选的，本实施例中选取做过自激自收时差校正的叠加道作为模型道，即图4右图中的叠加剖面各道集。图5左图为本实施例非一致性时差校正前的cmp道集波形图，图5右图为本实施例非一致性时差校正后的cmp道集波形图，经对比可见cmp道集的非一致性时差基本上得到了消除。如图6所示，相较于非一致性时差校正处理前（参见图6左图），经过非一致性时差校正处理后（参见图6右图）叠加剖面的信噪比和分辨率也得到了有效的提高。
36.步骤s5：对非一致性时差校正后的叠加剖面进行海底平滑时差校正。
37.进一步的，进行海底平滑时差校正的海底平滑时差校正量的计算公式如下：其中，表示进行海底平滑时差校正的海底平滑时差校正量，表示利用浮动基准面高程计算的第k个cmp道集的自激自收时间，为第k个cmp位置对应海底位置的海底深度，为海水平均声速。在本实施例中，海水平均声速取1481.9 m/s。
38.进一步的，为深度计和高度计测量值之和，经低通滤波后得到。
39.通常情况下，最终通过将叠加或偏移剖面从浮动基准面校正到水平基准面上，使得成像结果与地下地层真实构造相吻合。然而深拖式高分辨率多道地震数据所建立的浮动基准面并不一定与海底地形相关，只是深拖震源拖曳轨迹的平滑，所以采取在浮动基准面上进行cmp道集叠加成像。经过非一致性时差校正后的叠加剖面不再校正到固定基准面上，而是利用深度计和高度计测量值之和计算海底位置的海底深度，获得测线内海底地形，对计算结果进行低通滤波后与海水的平均声速计算海底自激自收反射双程时间，进行海底平滑校正。如图7所示，相较于海底平滑时差校正处理前（参见图7左图），经过海底平滑时差校正处理后（参见图7右图）的最终叠加剖面更符合海底以下地层的构造形态。
40.综上所述，本发明通过综合引入浮动基准面校正、cmp道集内剩余时差校正、自激自收时差校正、利用相位替代法的非一致性时差校正及海底平滑时差校正等，避免了基于地表一致性假设或常规非地表一致性剩余静校正方法均不适用于深拖式高分辨率多道地震探测的高频数据的问题，可以充分消除各类时差，能够实现海底以下地层构造形态的高分辨率、高信噪比、高保真度成像。
41.应当理解的是，除非本文中有明确的说明，本发明方法各个步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
42.以上对本公开的实施例进行了描述。但是，这些实施例仅仅是为了说明的目的，而并非为了限制本公开的范围。本公开的范围由所附权利要求及其等价物限定。不脱离本公开的范围，本领域技术人员可以做出多种替代和修改，这些替代和修改都应落在本公开的
范围之内。