海底电缆双检数据合并方法、装置、电子设备及介质制造方法及图纸

本申请公开一种海底电缆双检数据合并方法、装置、电子设备及介质，包括：采集陆检数据和水检数据；通过三维线性Radon将横波信号从所述陆检数据中分离，并通过三维高阶线性Radon保持有效波的能量；确定所述陆检数据和水检数据的振幅、相位的标定因子；根据标定因子合并陆检数据和水检数据。本发明专利技术通过利用三维高阶线性Radon，对变换过程中有效的保持了有效波的能量，获取水检数据和陆检数据的振幅、相位的标定因子，根据标定因子对水检数据和陆检数据进行合并，达到鬼波压制的目的。达到鬼波压制的目的。达到鬼波压制的目的。

【技术实现步骤摘要】
海底电缆双检数据合并方法、装置、电子设备及介质


[0001]本专利技术涉及油气物探工程
，更具体地，涉及一种海底电缆双检数据合并方法、装置、电子设备及介质。

技术介绍

[0002]双检合成技术作为海底电缆的配套处理技术最早出现于国外，国内于1997年引入海底电缆采集技术后，虽然也进行过相关的试验工作，但是由于处理技术的发展限制，双检四分量勘探系统未能得到很好的普及。直到2003年通过对此项技术进行攻克研究，拓展了四分量地震资料处理技术在国内的应用，使得海上勘探逐步活跃起来。国内很多地球物理专家也对双检资料合并技术进行了大量的探索工作，提出海底电缆(OBC)双检检波器采集技术可以提高地震资料的采集优良效果。
[0003]近几十年来，海上地震勘探技术也开始转向更难更深的领域，这对地震勘探的精度要求很高。尤其是我国南部和东部地区由于全球板块构造漂移，海底构造环境异常复杂，这些都给海上高精度地震勘探带来了巨大的挑战。海上地震勘探由于采集环境、检波器耦合情况等原因，噪声种类很多，比如海上涌浪，多次波，外源干扰等。多次波种类很多，有虚反射、微曲多次波、短程、全程多次波等，根据不同的传播路径和产生原因针对其特点应采用不同的方法进行压制。其中，鬼波和海上鸣震在海上地震资料中是比较常见的干扰，严重影响着有效信号的识别，对资料的成像和解释造成严重干扰。
[0004]因此，期待提出一种海底电缆双检数据合并方法，解决陆检数据的信噪比，水检数据和陆检数据合并过程中振幅、相位存在差异匹配性等的技术问题。
[0005]公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

技术实现思路

[0006]本专利技术提出了一种海底电缆双检数据合并方法、装置、电子设备及介质，至少解决陆检数据的信噪比，水检数据和陆检数据合并过程中振幅、相位存在差异匹配性等的技术问题的技术问题。
[0007]第一方面，本公开实施例提供了一种海底电缆双检数据合并方法，包括：
[0008]采集陆检数据和水检数据；
[0009]通过三维线性Radon将横波信号从所述陆检数据中分离，并通过三维高阶线性Radon保持有效波的能量；
[0010]确定所述陆检数据和水检数据的振幅、相位的标定因子；
[0011]根据所述标定因子合并陆检数据和水检数据。
[0012]作为本公开实施例的一种具体实现方式，所述三维线性Radon的变换公式如下：
[0013][0014][0015]其中，d为三维地震数据；p
x
,p
y
为纵、横向测线的路径曲率参数；
[0016]所述三维高阶线性Radon的变换公式为：
[0017][0018][0019][0020][0021]其中，d(x,y,t)为时间域三维地震数据，q
x
,q
y
为纵、横测线的曲率，p
j
(x
i
)表示第j阶多项式。
[0022]作为本公开实施例的一种具体实现方式，所述确定标定因子的方法包括：通过对所述水检数据与所述陆检数据进行相位、振幅等匹配，通过采集数据计算标定因子。
[0023]作为本公开实施例的一种具体实现方式，所述陆检数据的相位匹配公式如下：
[0024][0025]其中，φ为相位匹配因子；g[n]是陆检数据序列；g
H
[n]是陆检数据的希尔伯特地震道序列，表示为g
H
[n]＝g[n]*h
w
[n]，*表示褶积；
[0026]相位匹配因子公式如下：
[0027][0028]相位匹配因子特征方程系数的公式为：
[0029][0030]式中，A[i]为水检数据和陆检数据的互相关；B[i]为水检数据和陆检数据希尔伯特地震道的互相关函数，表示成：
[0031][0032][0033]在这个式子中，e[k]是水检数据；N是地震数据序列的长度；i＝0,1,2,L,L，L是互相关函数的长度；根据匹配因子φ，带入式(2.1)实现陆检数据相位匹配。
[0034]作为本公开实施例的一种具体实现方式，所述水检数据和陆检数据合并公式如下：
[0035][0036]其中，S
i,j
是水检数据和陆检数据合并数据；H
i,j
是水检数据序列；G
i,j
是陆检数据序列的长度；i是时窗参数的顺序号，i＝0,1,2,L,L，L是时窗参数的总道数，j＝0,1,2,L,N，N是时窗参数时间样点数，j是时窗参数时间样点的顺序号，α是标定因子；是水检数据和陆检数据的均方根振幅比其计算公式为：
[0037][0038]为计算标定因子，使用最大方差模(Varimax)准则，其计算公式为：
[0039][0040]其中，P
n
为方程分子的系数，Q
n
为方程分母的系数；可得到标定因子α的15阶特征方程：
[0041]T0+T1α+T2α2+
…
+T
15
α
15
＝0
ꢀꢀ
(2.9)
[0042]式中T
n
为特征方程的系数。
[0043][0044][0045]对上式进行求解，得到标定因子α方程的15个根α
n
，n＝1,2,3,...15，将这些根带入式(2.3)，即可计算出15个对应的最大方差模Varm(α
n
)；
[0046]α
best
＝Max{Varm(α1),Varm(α2),Varm(α3),...,Varm(α
15
)}
[0047]其中最大值即为标定因子。
[0048]第二方面，本公开实施例还提供了一种海底电缆双检数据合并装置，包括：
[0049]采集模块，采集陆检数据和水检数据；
[0050]去噪模块，通过三维线性Radon将横波信号从所述陆检数据中分离，并通过三维高
阶线性Radon保持有效波的能量；
[0051]获取模块，确定所述陆检数据和水检数据的振幅、相位的标定因子；
[0052]合并模块，根据所述标定因子合并陆检数据和水检数据。
[0053]作为本公开实施例的一种具体实现方式，所述三维线性Radon的变换公式如下：
[0054][0055][0056]其中，d为三维地震数据；p
x
,p
y
为纵、横向测线的路径曲率参数；
[0057]所述三维高阶线性Radon的变换公式为：
[0058][0059][0060][0061][0062]其中，d(x,y,t)为时间域三维地震数据，q
x
,q
y
为纵、横测线的曲率本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种海底电缆双检数据合并方法，其特征在于，包括：采集陆检数据和水检数据；通过三维线性Radon将横波信号从所述陆检数据中分离，并通过三维高阶线性Radon保持有效波的能量；确定所述陆检数据和水检数据的振幅、相位的标定因子；根据所述标定因子合并陆检数据和水检数据。2.根据权利要求1所述的海底电缆双检数据合并方法，其特征在于，所述三维线性Radon的变换公式如下：所述三维线性Radon的变换公式如下：其中，d为三维地震数据；p
x
,p
y
为纵、横向测线的路径曲率参数；所述三维高阶线性Radon的变换公式为：所述三维高阶线性Radon的变换公式为：所述三维高阶线性Radon的变换公式为：所述三维高阶线性Radon的变换公式为：其中，d(x,y,t)为时间域三维地震数据，q
x
,q
y
为纵、横测线的曲率，p
j
(x
i
)表示第j阶多项式。3.根据权利要求1所述的海底电缆双检数据合并方法，其特征在于，所述确定标定因子的方法包括：通过对所述水检数据与所述陆检数据进行相位、振幅等匹配，通过采集数据计算标定因子。4.根据权利要求3所述的海底电缆双检数据合并方法，其特征在于，所述陆检数据的相位匹配公式如下：其中，φ为相位匹配因子；g[n]是陆检数据序列；g
H
[n]是陆检数据的希尔伯特地震道序列，表示为g
H
[n]＝g[n]*h
w
[n]，*表示褶积。5.根据权利要求4所述的海底电缆双检数据合并方法，其特征在于，所述水检数据和陆检数据合并公式如下：其中，S
i,j
是水检数据和陆检数据合并数据；H
i,j
是水检数据序列；G
i,j
是陆检数据序列的长度；i是时窗参数的顺序号，i＝0,1,2,L,L，L是时窗参数的总道数，j＝0,1,2,L,N，N是
时窗参数时间样点数，j是时窗参数时间样点的顺序号，α是标定因子；是水检数据和陆检数据的均方根振幅比为计算标定因子，使用最大方差模(Varimax)准则，其计算公式为：其中，P
n
为方程分子的系数，Q
n
为方程分母的系数；可得到标定因子α的15阶特征方程：T0+T1α+T2α2+
…
+T
15
α
15
＝0式中T...

【专利技术属性】
技术研发人员：于鑫，郭恺，
申请(专利权)人：中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：