一种井控区深度域陆上层间多次波直接压制方法及装置与流程

本公开涉及地震资料处理，尤其涉及一种井控区深度域陆上层间多次波直接压制方法及装置。

背景技术：

1、地震资料处理中提高成果保真度是首要目标，陆上地震资料通常存在层间多次波，层间多次波严重降低处理成果的保真度甚至影响地震预测结果。去除多次波的方法很多，最常用的是应用拉冬变换在τ-p域压制多次波，拉冬变换由奥地利数学家radon在1917年提出。claerbout率先将τ-p变换引入到地震勘探中来；hampson首次提出抛物线拉冬变换在cmp道集上压制多次波；kostov将抛物线拉冬变换最小二乘化，使其变换算子具有toeplitz结构，可快速求解；wang应用radon变换域自适应滤波提高了拉冬变换的分辨率；ng利用高斯迭代法计算相似系数加权，对时间轴和偏移距轴两个方向同时稀疏，显著地提高了拉冬变换的分辨率；sacchi和ulrych在最小二乘抛物线拉冬变换的算子中引入非线性规则化对角矩阵，提出的高分辨拉冬变换解决拉冬域数据发散问题，提高拉冬域多次波与一次波的分辨能力；mauricio等提出共轭梯度法实现高分辨率拉冬变换的快速求解；herrmann等提出去假频高分辨率抛物线拉冬变换，解决分辨率问题和空间假频问题；熊登等将以上理论应用于多次波的衰减中，取得良好效果；刘喜武等采用最小二乘反演方法研究抛物线radon变换和双曲线radon变换，给出稀疏约束预条件的高分辨率radon变换的算法，压制多次波效果良好。上述方法中大部分基于拉东变换的压制方法都与地震频率有关，尤其是高精度拉东变换是在频率域逐一频率进行的拉东变换与逆变换。

2、在大量陆上实际地震资料的处理中，我们发现叠前深度偏移比叠前时间偏移在复杂构造刻画能力方面更强，叠前深度偏移成果剖面和深度域crp道集保真度都更高，但上述压制多次波的方法都是在时间域实现的，深度偏移的成果要通过速度模型进行时深转换到时间域，在时深转换过程中就会出现这样或那样的畸变，导致时间域道集对陆上层间多次波的刻画没有深度域道集那么精确。而如果粗暴的将深度偏移成果的单位m直接替换为ms后进行拉东变换，会出现频率失真和逆变换产生大量畸变的问题，所以不能简单地将深度偏移的成果利用时间域的处理办法进行处理。综上所述，急需找到一种能够在深度域压制陆上层间多次波的方法。

技术实现思路

1、本公开提出了一种井控区深度域陆上层间多次波直接压制方法及装置，以解决以往通过拉冬变换进行压制层间多次波时，需要通过时深转换将深度域偏移处理成果转换到时间域，在时深转换过程中会出现畸变，导致得到的时间域道集对陆上层间多次波的刻画不够精确，不利于后续压制处理的问题。

2、根据本公开的一方面，提供了一种井控区深度域陆上层间多次波直接压制方法，包括：

3、获取工区测井数据及叠前深度偏移地震成果叠加剖面，其中，测井数据深度单位及叠前深度偏移地震成果叠加剖面数据的单位都为m；

4、确定工区内的典型井，根据所述典型井的测井曲线，利用考虑传播效应正演公式，得到深度域单位为m的去除一次反射波信息的叠前crp道集；

5、对所述去除一次反射波信息的叠前crp道集上的多次波的位置进行标定，并将标定后的多次波叠前crp道集的深度单位m替换为ms；

6、对深度单位为ms的标定后的多次波crp道集，通过速度分析进行nmo拉平及切除远道拉伸部分处理，形成多次波叠加剖面；

7、将多次波叠加剖面的深度单位的ms替换为m，将所述叠前深度偏移地震成果叠加剖面与所述深度单位为m的多次波叠加剖面相减，得到相减后的叠加剖面，并将相减后的叠加剖面的深度单位的m替换为ms；

8、对所述深度单位为ms的相减后的叠加剖面进行脉冲反褶积处理，得到井控区深度域陆上层间多次波直接压制后的叠前深度偏移成果。

9、优选地，所述确定工区内的典型井的方法，包括：

10、确定工区内所有井中，垂直深度最深的井为所述典型井。

11、优选地，所述根据所述典型井的测井曲线，利用考虑传播效应正演公式，得到深度域单位为m的去除一次反射波信息的叠前crp道集的方法，包括：

12、将所述典型井的测井曲线的纵横波速度代入考虑传播效应正演公式，得到单位为m的深度域波动方程正演全波场crp道集；

13、将所述典型井声波曲线代入考虑传播效应正演公式，忽略除一次反射以外的全部传播效应，正演得到单位为m的深度域一次反射波crp道集；

14、利用深度域波动方程正演全波场crp道集，减去深度域一次反射波crp道集，得到去除一次反射波信息的叠前crp道集。

15、优选地，所述对去除一次反射波信息的叠前crp道集上的多次波的位置进行多次波标定的方法，包括：

16、若所述去除一次反射波信息的叠前crp道集上的同相轴发生弯曲现象，则将该同相轴对应偏移距为0的深度位置，标记为多次波存在位置，从而得到标定后的多次波叠前crp道集。

17、优选地，所述对深度单位为ms的标定后的多次波crp道集，通过速度分析进行nmo拉平及切除远道拉伸部分处理，形成多次波叠加剖面的方法，包括：

18、根据所述标定后的多次波crp道集，确定加权叠加地震道；

19、对所述加权叠加地震道进行速度分析，找到nmo速度；

20、利用所述nmo速度，对所述标定后的多次波叠前crp道集进行拉平，得到拉平后的crp道集；

21、对拉平后的crp道集进行切除远道拉伸部分处理，得到切除后的crp道集；

22、对所述切除后的crp道集进行水平叠加，形成多次波叠加剖面。

23、优选地，所述将拉平后的crp道集进行切除远道拉伸部分处理，得到切除后的crp道集的方法，包括：

24、对拉平后的crp道集，切除其偏移距的前第一预定百分数和后第二预定百分数，得到多次波叠加剖面。

25、优选地，所述将相减后的叠加剖面的深度单位的m替换为ms之前，还包括：对所述相减后的叠加剖面进行质控，其方法，包括：

26、对所述相减后的叠加剖面与多次波叠加剖面进行互相关计算，判断其互相关结果是否小于第三预定百分数，若为是，则质控通过，否则，质控不通过。

27、优选地，所述脉冲反褶积处理的反滤波算子的长度m为：若工区地震记录总道长小于或等于10000ms，则反滤波算子长度为400ms，若工区地震记录总道长大于10000ms，则反滤波算子长度为600ms。

28、优选地，还包括：对得到的所述叠前深度偏移成果进行质控，其方法，包括：

29、将所述叠前深度偏移成果的深度单位的ms替换为m；

30、对深度单位为m的所述叠前深度偏移成果与多次波叠加剖面进行互相关计算，判断其互相关结果是否小于第四预定百分数，若为是，则质控通过，否则，对所述相减后的叠加剖面进行新一轮迭代的脉冲反褶积处理，直至得到的叠前深度偏移成果和多次波叠加剖面的互相关结果小于第四预定百分数。

31、根据本公开的一方面，提供了一种井控区深度域陆上层间多次波直接压制装置，包括：

32、获取单元，用于获取工区测井数据及叠前深度偏移地震成果叠加剖面，其中，测井数据深度单位及叠前深度偏移地震成果叠加剖面数据的单位都为m；

33、多次波叠前crp道集确定单元，用于确定工区内的典型井，根据所述典型井的测井曲线，利用考虑传播效应正演公式，得到深度域单位为m的去除一次反射波信息的叠前crp道集；

34、多次波标定单元，用于对所述去除一次反射波信息的叠前crp道集上的多次波的位置进行多次波标定，并将标定后的多次波叠前crp道集的深度单位m替换为ms；

35、多次波叠加剖面确定单元，用于对深度单位为ms的标定后的多次波crp道集，通过速度分析进行nmo拉平及切除远道拉伸部分处理，形成多次波叠加剖面；

36、一次反射波叠加剖面确定单元，用于将多次波叠加剖面的深度单位的ms替换为m，将所述叠前深度偏移地震成果叠加剖面与所述深度单位为m的多次波叠加剖面相减，得到相减后的叠加剖面，并将相减后的叠加剖面的深度单位的m替换为ms；

37、叠前深度偏移成果确定单元，用于对所述深度单位为ms的相减后的叠加剖面进行脉冲反褶积处理，得到井控区深度域陆上层间多次波直接压制后的叠前深度偏移成果。

38、本发明至少具有如下有益效果：

39、本公开提出了一种井控区深度域陆上层间多次波直接压制方法及装置，利用测井数据通过波动方程正演分离出多次波叠前crp道集，通过对多次波crp道集进行速度分析后找到nmo速度，利用nmo速度拉平并切除远道拉伸部分，叠加后形成多次波叠加剖面，从原始psdm剖面减去多次波剖面后再进行脉冲反褶积处理。本公开全程不进行时深转换，而是直接对多次波刻画相对更清晰的叠前深度偏移地震成果叠加剖面的多次波进行压制，可以避免深度域和时间域来回转换时产生的畸变，从而得到保真度更高的多次波压制后叠前深度偏移成果。