基于多t-x域联合的多震源地震数据混叠噪声分离方法

俞 岱，黄德智，孙 渊，杨文平

基于多t-x域联合的多震源地震数据混叠噪声分离方法

俞 岱1，黄德智2，孙 渊1，杨文平3

(1. 长安大学 地质工程与测绘学院，陕西 西安 710054；2. 黑龙江科技大学 矿业工程学院，黑龙江 哈尔滨 150022；3. 中石化石油工程地球物理有限公司华北分公司，河南 郑州 450000)

随着多震源同时激发采集方法在地震油气勘探领域的应用逐年增加，许多学者对其关键技术多震源地震数据混叠噪声分离进行了研究，就目前为止，仍需开展保真度更高的方法研究。基于以下设想：对多震源地震数据做正常时差校正(NMO)，可以加强反射信号的线性分布程度(尤其在共中心点CMP道集内)，并扩大反射信号和混叠噪声的分布差异；单域分离混叠噪声，仅针对混叠噪声和反射信号的某一分布特性差异进行分离，多t-x域联合分离可综合利用各t-x域内混叠噪声和反射信号在各域内的分布差异，其分离效果更好。首先，通过对多震源数据CMP道集做NMO以强化反射信号的线性程度、扩大随机分布混叠噪声的离散程度，利用中值滤波分离大部分混叠噪声，再在其他t-x域内(如共炮点CSG道集、共检波点CRG道集和共偏移距COG道集等)根据残余混叠噪声分布特征利用随机噪声衰减等方法进一步分离，并在整个分离处理中使用保幅处理方法。提出的多t-x域联合分离混叠噪声方法，经理论数据验证，相较于单域分离方法可以有效且保真分离多震源数据中的混叠噪声，并对其他噪声和横波有一定的压制作用；实际数据应用效果表明，该方法比单域方法分离效果较好，叠后成像效果品质较高。多t-x域联合分离混叠噪声方法，能够稳定、可靠和保真地分离多震源地震数据中的混叠噪声。

多震源地震采集；混叠噪声分离；多t-x域联合；噪声离散化

D. Silverman[1]在1979年首先提出多震源同时激发采集方法，并成功从多震源采集数据中分离出单炮数据；C. J. Beasley[2]引入了脉冲型震源的同时激发采集；C. Bagaini[3]对比讨论了各种不同可控震源同时激发采集方法。多震源采集技术在20世纪70年代末提出并在野外试验成功后，相较于传统地震采集技术，其高效性使其广泛应用于实际生产当中。但在采集过程中有多震源的混叠，使其存在严重的混叠噪声干扰问题，严重制约了地震信号的成像。

目前，混叠噪声分离主要有反演和滤波去噪2种分离方法[4]。反演法分离多震源数据是将数据分离过程视为数学物理反演问题；滤波法分离多震源数据是根据混叠噪声在不同域内与有效信号之间的分布特性差异来对混叠噪声进行分离。多位学者通过滤波法实现了混叠噪声分离，主要有：在伪分离道集的共偏移距道集进行滤波去噪[5-6]；在非炮域中使用矢量中值滤波方法分离[7]；根据反射信号和混叠噪声在不同时域的相干性差异使用迭代法分离[8]；通过调整阈值在共检波点域使用FK滤波方法进行分离[9-10]；通过联合Curvelet阈值迭代以及多级中值滤波进行混叠噪声分离[11]；利用动校正速度对局部时窗内同相轴斜率进行估算，据此拉平同相轴并通过中值滤波法分离[12]；采用局部正交化加权信号和噪声方法恢复去噪过程中损失的信号，并设计出将伪动校正、中值滤波和F-K滤波3种方法结合在一起的迭代算法[13-14]；利用自适应迭代多级中值滤波方法，在迭代过程中采用自适应减小中值滤波窗口方法，对海上多震源混合波场的分离效果好[15]；利用基于脉冲检测的混采波场分离方法，并与迭代的多级中值滤波方法作对比，时间延时范围较小时，脉冲检测法更有效[16]；联合中值滤波、动校正、复曲波域或seislet迭代阈值去噪方法，设计出改进型的混叠噪声压制技术流程，其综合了滤波方法计算效率快和反演类方法噪声压制效果精度高的优点[17-18]；利用基于奇异值分解约束迭代反演的混叠噪声压制方法，在共炮检距域或动校正后的共中心点域实现混叠噪声与有效信号的分离[19]。

上述方法在混叠噪声分离方面取得了不错的效果，但针对混叠噪声离散程度方面的研究较少，多震源数据中反射信号在动校正(Normal Moveout Correction，NMO)后其线性分布程度会得到有效提高，尤其共中心点(Common Middle Point， CMP)道集内的反射信号经NMO后会变成线性水平分布，同时混叠噪声分布的离散程度会显著增大；在NMO前的道集中混叠噪声的分布形态类似于单炮数据，但在NMO后的非炮集时空(t-x)域地震数据中，这种分布特性被破坏，使混叠噪声和有效反射波之间的分布特性出现本质性差异。基于此，笔者提出多t-x域联合分离方法，通过联合NMO和多域混合去噪方法，对混叠噪声进行分离。该技术为多震源地震数据分离提供新的思路。

1 基本原理

1.1 混叠噪声离散化

CMP道集中反射信号在时空域内呈双曲线形态展布的二次函数，NMO后在时空域其会被拉成水平线性分布，如图1所示。

图1 CMP道集反射波

相同反射层处于不同偏移距的反射波，NMO可将相同反射层同相轴校至自激自收时间，式(1)为均匀水平层状介质下NMO时差公式：

由式(1)推导可知，对于相同非零偏移距位置，不同时刻的2个采样点、(t>t，h>h，t0>t0)，两点间NMO后的时差要大于它们之间的原有时差：

式中：t、t分别为采样点、的记录时间；t0、t0分别为采样点、在NMO后对应的自激自收时间； h、h分别为采样点、对应的地震反射深度；由几何关系推导可知不等式关系成立，如图1a所示。

由式(1)可知，同一反射深度的2个不同偏移距采样点、，偏移距x>x，大偏移距采样点的NMO时差大于小偏移距采样点的NMO时差：

综合式(2)、式(3)可知，对于非零偏移距位置的随机噪声而言，NMO后相同偏移距位置不同时刻的随机噪声之间的时差会变大；对于同一反射点上的随机噪声，其时差会随着偏移距的增大而增大。对于CMP道集中相同反射层位的有效反射信号而言， NMO后同相轴呈水平线性；对于其他t-x域地震数据中的反射信号而言，NMO后其线性分布也会得到不同程度提高。

在非NMO炮集中，混叠噪声的分布特性与正常单炮数据的分布特性非常相似，但在NMO炮集中，混叠噪声与正常单炮数据的分布特性会产生本质性差异。如图2所示，图2a为非NMO炮集数据，图2b为NMO后炮集数据，图2a中位于左侧的单炮数据为有效数据，位于右侧的单炮数据为本炮的混叠噪声。对比两图可以看出，经NMO后，图2a中的有效数据的反射波同相轴基本被拉成线性，而混叠噪声中的同相轴则变成反向双曲线形态(图2b)。故NMO后有效信号的分布特征与混叠噪声的分布特征会产生明显区别。

因此，对多震源数据中的反射信号进行NMO，可以提高道集中反射信号的线性化程度，扩大道集中随机分布混叠噪声的离散化程度以及多震源地震数据中混叠噪声与有效反射信号间的分布特征差异。

图2 多震源炮集NMO前后记录对比

1.2 去噪处理

1.2.1 中值滤波

中值滤波是一种基于排序统计理论抑制噪声的非线性信号处理技术。它将数字图像或数字序列中一点的值用该点一个邻域中各点值的中值代替，让该点的像素值接近真实值，从而消除孤立的噪声点。中值滤波对脉冲噪声有良好的滤除作用，特别是在滤除噪声的同时，能够保护信号的边缘，使之不被模糊[20]。

在混叠噪声随机分布的t-x域地震数据中，混叠噪声是以孤立脉冲噪声形式存在的。因此，中值滤波是一种比较适合消除该域中混叠噪声的方法。

1.2.2 振幅一致性处理(保幅处理)

地震波在传播过程中，因表层地质条件变化、激发震源能量差异、接收检波器与周围介质耦合程度的差异以及因传播距离而造成的能量损失等因素，会造成相同道集内的不同地震道之间出现振幅一致性差异。这种振幅非一致性问题会造成中值滤波后地震道的振幅属性改变，使得滤波后地震数据不能真实反映原数据所代表的真实地球物理信息。因此，在中值滤波前，要对道集内地震数据做振幅一致性处理，将道集内所有地震道的振幅调整到同一水平；滤波后再对所有地震道进行反振幅一致性处理，将滤波后地震道振幅水平恢复回原地震道振幅水平。这样会使得滤波前后地震数据保持相同振幅水平，从而达到保幅处理的目的。

对于一个包含个地震道的道集，首先统计出各道内n个非零采样点振幅绝对值(n为第道内非零采样点个数)的平均值，将其作为各道振幅一致性参考值vi；然后取各道振幅一致性参考值的平均值作为道集振幅一致性参考值v，求取各道振幅一致性参考值与道集振幅一致性参考值的比值作为各道的振幅一致性系数R，并将各道采样点振幅除以对应地震道振幅一致性系数，从而使地震道内各道能量处于同一振幅级别。在滤波后，各道内采样点振幅再乘以对应地震道振幅一致性系数，从而恢复各道原有振幅。

1.2.3 随机噪声衰减

根据f-x 域预测理论，若多个信号线性相干，利用前一道统计出空间预测算子可预测下一道。如果信号只有一个视速度，且各道信号的振幅谱相同，则相邻道的相位差对于固定的频率来说是常数。前一道的某个频率成分乘上适当的固定相移因子就等于后一道的同一频率成分[21-23]。

实际地震数据中，因噪声的存在使得相邻道的同一频率成分会不规则地偏离上述关系。采用多点预测算子，即使用多道来预测单道。通过这种方法，可以显著减小各道间的不规则程度，再将预测结果经过傅里叶反变换回时间域输出，即得到去噪后信号，从而达到衰减随机噪声的目的。

对于有个不同视速度的信号，预测算子长度应等于或大于个点。若一个信号的道间时差为，第道的频谱为：

第1道的频谱为：

因此，由傅里叶级数展开为：

式中：a为离散后的振幅，=1, 2, …,。

要求极小，对于同一圆频率预测算子是二维的，得到方程组：

对得到的各频率算子应用于实际数据，做f-x预测，再做傅里叶反变换后输出，得到去噪结果。

1.3 多t-x域联合分离

据前述可知，多震源地震数据经NMO后，CMP道集内的反射信号变成水平线性，并有效扩大随机分布混叠噪声的离散程度，使得NMO后CMP道集成为最有利于分离混叠噪声的道集。其他道集地震数据NMO后反射信号的线性程度也会得到有效提高，其中随机分布混叠噪声的离散程度也会得到扩大。

多t-x域联合分离方法实现步骤，如图3所示，CMP道集经过NMO后，反射信号成水平线性，通过中值滤波很容易将大部分随机分布的混叠噪声去除，在此基础上，根据残余混叠噪声在其他t-x域地震数据内的分布特征，利用随机噪声衰减等去噪方法进一步分离，以达到多t-x域联合分离的目的。

图3 多t-x域联合分离方法流程

2 理论数据验证

为验证本文所提方法，先使用理论数据合成多震源地震数据进行验证[24-25]，对理论数据进行混叠噪声分离，并将分离后地震数据与合成前数据进行对比。建立5层地质模型(图4a)，地层纵波速度分别为1 700、1 900、2 100、2 300、2 500 m/s，模型的纵横波速比为1.36，各向同性介质弹性波正演采用Tesseral软件，30 Hz雷克子波，固定排列接收，每炮251道，道间距20 m。图4是合成多震源地震数据示意图，图4d数据为图4b数据和图4c数据合成的多震源地震数据，假定图4b数据在合成数据中为混叠噪声。

理论数据验证包含两方面内容，一是多震源地震数据NMO后随机混叠噪声的离散化程度扩大验证；二是多t-x域联合分离混叠噪声能力与单域分离能力对比验证。

2.1 随机分布混叠噪声离散化扩大验证

混叠噪声在CMP道集、共检波点(Common Receiver Gather，CRG)道集和共偏移距(Common Offset Gather，COG)道集中是以随机噪声的形态分布的，如图5所示，除有效反射信号(红色箭头指示处)外，其余绝大部分信号成分横向不连续，呈随机分布。因此，多震源数据中随机混叠噪声的离散化扩大验证主要在这3种道集中进行，分别对这3种道集进行NMO，以验证NMO对混叠噪声的离散程度扩大效果[26-27]。

图4 多震源地震数据合成

图5a、图5b、图5c分别为CMP道集数据、CRG道集数据和COG道集数据，图5d—图5f分别为对应道集数据NMO后结果。对比各道集数据NMO前后结果，在NMO后，3个道集内随机分布混叠噪声的离散化程度都得到有效扩大。但CMP道集和CRG道集中随机分布混叠噪声的离散化程度，明显大于常规分离方法中经常使用的COG道集(图5c)中的离散化程度。这是因为COG道集相对于其他道集，少了偏移距变化这个变量，在单个COG道集中，混叠噪声离散化的比例基本一致；但在不同的COG道集中，混叠噪声的离散比例是不同的，偏移距越小的COG道集，NMO后混叠噪声离散化越弱，反之则越强。CMP道集NMO后相同层位的有效反射波呈完全水平线性分布，CRG道集和COG道集NMO后，有效反射波同相轴分布的线性程度由其所对应地下反射地层的形态决定。

NMO后道集的随机混叠噪声分离能力较NMO前有很大提高。图6a数据是使用未做NMO的COG道集分离混叠噪声结果，图6b是对NMO后CMP道集分离混叠噪声结果。对比两图可见，NMO后CMP道集相较未做NMO的COG道集，其混叠噪声分离能力明显提高。

图5 NMO前后随机分布混叠噪声在多震源地震数据不同道集分布

图6 2种道集混叠噪声分离效果对比

2.2 多域联合分离验证

针对图4所示的理论数据进行多域联合分离验证，图7为对第50炮理论数据(图7a)进行多t-x域联合分离混叠噪声的过程。按照前述方法步骤，利用速度分析结果对CMP道集进行NMO处理，在此基础上，通过中值滤波分离大部分混叠噪声(图7b)，根据混叠噪声分离后结果在不同t-x域(图7c，以CSG道集为例)内分布特征，进一步利用随机噪声衰减等去噪方法分离残余混叠噪声(图7d)，最终得到多t-x域联合分离结果(图7e)。图7中的频谱是对红色时窗内反射波信号分析的结果，对比多t-x域联合分离前后频谱可见，主频均为30 Hz，时窗内最大振幅值均为3.7，证明该方法的保真度较高。

(a) 理论合成多震源数据；(b) NMO后CMP道集去噪结果；(c) 对b抽CSG道集；(d) CSG道集进一步去噪结果；(e) 多t-x域联合分离结果

图8为理论数据多t-x域联合与单域分离混叠噪声效果对比图。由图8c和图8d可看出，多t-x域联合方法分离混叠噪声能力明显高于单域方法(仅在未做NMO的COG道集中分离混叠噪声)，单域方法分离后数据中残留有严重的混叠噪声，而多t-x域联合方法分离后数据中混叠噪声得到了有效分离。

从图8a、图8c和图8d中提取第120道数据，绘制混叠噪声分离前后波形对比图(图9)，对比可见，多域分离结果中除完整保留4个主要反射信号波形(1 028、1 701、2 407、2 890 ms)外，其余混叠噪声的波形均被有效压制，单域分离结果中还残余一些混叠噪声的波形。对图中展示时间段内波形进行信噪比计算，信噪比计算公式为：

图8 理论数据多t-x域联合与单域分离混叠噪声效果对比

图9 理论数据分离混叠噪声前后第120道波形对比

分离前理论数据波形信噪比为–0.000 96，单域分离结果波形信噪比为0.083 177，多t-x域联合分离结果波形信噪比为0.098 696，从量化分析结果看多t-x域联合分离效果更好。对比图8a、图8b、图8d以及图9可以看出，多t-x域联合分离混叠噪声后，原数据中的有效反射波信号得到很好的还原，混叠噪声得到有效分离，且原数据中的横波反射也得到了很好的压制。因此，多t-x域联合分离混叠噪声方法，可以有效分离多震源地震数据中的混叠噪声，且对数据中的其他噪声和横波有一定的压制作用。

3 实际数据分离验证

为验证该方法对实际多震源地震数据的分离能力，选取了一段传统方法施工的二维地震数据，所用数据的野外采集参数为：181道接收、道间距30 m，排列固定不动，炮点滚动激发，炮间距60 m，共45炮，将其人工合成为多震源地震数据进行混叠噪声分离处理。

图10为实际数据分离混叠噪声效果对比图。图10c为图10a数据与图10b数据合成的多震源地震数据，图10d为单域分离方法(与前同)分离后结果，图10e为多t-x域联合分离方法分离后结果。对比图10d与图10e可看出单域分离方法得到结果中残存有严重的混叠噪声，多域分离方法中混叠噪声得到了很好的分离。图10b为图10c中的有效炮部分，对比图10b与图10e可以看出，图10b中其他噪声也得到有效压制。

在前述2种分离混叠噪声方法得到的单炮记录基础之上，进一步完成后续成像处理，如图11所示，为实际数据分离混叠噪声后叠加剖面对比图，对比两剖面可以看出，多t-x域联合分离方法得到的叠加剖面反射同相轴清晰，且浅部被混叠噪声压制的弱能量同相轴也突显出来，剖面成像品质更高，证明该分离方法稳定、可靠。

图10 实际数据分离混叠噪声效果对比

图11 实际数据混叠噪声分离效果剖面对比

通过以上两方面对比，充分说明本方法能可靠、稳定分离多震源地震数据中的混叠噪声。

4 结论

a. 通过理论与实际数据验证，提出的多t-x域联合分离混叠噪声方法能有效分离多震源数据中的混叠噪声；相较于单域(在未做NMO的COG道集内分离)方法，分离结果信噪比更高。

b. 该方法的优点在于，多震源数据经NMO后，可以有效改善反射波在各t-x域地震数据内的线性分布，扩大各种t-x域道集内混叠噪声分布的离散化程度。该方法应用中值滤波，在滤波前需对道集内各地震数据做振幅一致性处理，滤波后再进行反一致性处理，以达到保留原始地震单炮道间振幅关系的目的，为后续利用传统地震处理方法提供数据保障。

c.本次未对多震源数据采集的混合度对多t-x域内混叠噪声随机分布特征的影响进行讨论，在今后工作中应针对其做进一步的研究。

[1] SILVERMAN D. Method of three dimensional seismic prospecting：US4159463[P].1979-06-26.

[2] BEASLEY C J. A new look at marine simultaneous source[J]. The Leading Edge，2008，27(7)：914–917.

[3] BAGAINI C. Overview of simultaneous vibroseis acquisition methods[C]//SEG Technical Program Expanded Abstracts 2006. New Orleans：Society of Exploration Geophysicists，2006：70–74.

[4] 邬达理. 多源地震资料重叠炮信息消除方法研究[J].石油物探，2013，52(5)：519–523.

WU Dali. Elimination of duplicate shots information of multi-sources seismic data[J]. Geophysical Prospecting for Petroleum，2013，52(5)：519–523.

[5] MOORE I，DRAGOSET B，OMMUNDSEN T，et al. Simultaneous source separation using dithered sources[C]//SEG Technical Program Expanded Abstracts Society of Exploration Geophysicists. 2008. Las Vegas，2008：2806–2810.

[6] AKERBERG P，HAMPSON G，RICKETT J，et al. Simultaneous source separation by sparse radon transform[C]//SEG Technical Program Expanded Abstracts 2008. Las Vegas：Society of Exploration Geophysicists，2008：2801–2805.

[7] HUO Shoudong，LUO Yi，KELAMIS P G. Simultaneous sources separation via multi-directional vector-median filtering[C]//SEG Technical Program Expanded Abstracts 2009. Houston：Society of Exploration Geophysicists，2009：31–35.

[8] MAHDAD A，BLACQUIÈRE G. Iterative method for the seperation of blended encoded shot records[C]//72ndEAGE Conference and Exhibition incorporating SPE EUROPEC 2010. Barcelona：European Association of Geoscientists & Engineers，2010：1–4.

[9] DOULGERIS P，MAHDAD A，BLACQUIÈRE G. Separation of blended data by iterative estimation and subtraction of interference noise[C]//DOULGERIS P. SEG Technical Program Expanded Abstracts 2010. Denver：Society of Exploration Geophysicists，2010：3514–3518.

[10] MAHDAD A，DOULGERIS P，BLACQUIRE G. Separation of blended data by iterative estimation and subtraction of blending interference noise[J]. Geophysics，2011，76(3)：9–17.

[11] 韩立国，谭尘青，吕庆田，等. 基于迭代去噪的多源地震混合采集数据分离[J]. 地球物理学报，2013，56(7)：2402–2412.

HAN Liguo，TAN Chenqing，LYU Qingtian，et al. Separation of multi-source blended seismic acquisition data by iterative denoising[J]. Chinese Journal of Geophysics(in Chinese)，2013，56(7)：2402–2412.

[12] GAN Shuwei，WANG Shoudong，CHEN Xiaohong，et al. Deblending using a structural-oriented median filter[C]// GAN Shuwei. SEG Technical Program Expanded Abstracts 2015. New Orleans：Society of Exploration Geophysicists，2015：59–64.

[13] CHEN Yangkang，FOMEL S. Random noise attenuation using local signal-and-noise orthogonalization[J]. Geophysics，2015，80(6)：WD1–WD9.

[14] CHEN Yangkang，JIN Zhaoyu，GAN Shuwei，et al. Iterative deblending using shaping regularization with a combined PNMO-MF-FK coherency filter[J]. Journal of Applied Geophysics，2015，122(5)：18–27.

[15] 周丽，庄众，成景旺，等. 利用自适应迭代多级中值滤波法分离海上多震源混合波场[J]. 石油地球物理勘探，2016，51(3)：434–443.

ZHOU Li，ZHUANG Zhong，CHENG Jingwang，et al. Multi-source blended wavefield separation for marine seismic based on an adaptive iterative multi-level median filtering[J]. Oil Geophysical Prospecting，2016，51(3)：434–443.

[16] 魏亚杰，韩立国，单刚义，等. 基于脉冲检测的混合震源数据分离[J]. 地球物理学报，2018，61(3)：1157–1168.

WEI Yajie，HAN Liguo，SHAN Gangyi，et al. Separation of mixed source seismic data based on the impulse detection[J]. Chinese Journal of Geophysics(in Chinese)，2018，61(3)：1157–1168.

[17] 董烈乾，张慕刚，周大同，等. 一种优化的复曲波变换压制混叠噪声方法[J]. 地球物理学进展，2019，34(2)：517–522.

DONG Lieqian，ZHANG Mugang，ZHOU Datong，et al. Blended noise suppression using an optimized complex curvelet transform approach[J]. Progress in Geophysics(in Chinese)，2019，34(2)：517–522.

[18] 董烈乾，周恒，郭善力，等. 一种改进型seislet域迭代阈值压制混叠噪声方法[J]. 物探与化探，2020，44(3)：568–572.

DONG Lieqian，ZHOU Heng，GUO Shanli，et al. An optimized blending noise suppression based on seislet domain iterative threshold denoising approach[J]. Geophysical and Geochemical Exploration，2020，44(3)：568–572.

[19] 董烈乾，张慕刚，骆飞，等. 应用SVD约束的迭代反演混叠噪声压制方法[J]. 石油地球物理勘探，2021，56(1)：57–61.

DONG Lieqian，ZHANG Mugang，LUO Fei，et al. Blended noise suppression based on SVD constrained iterative inversion[J]. Oil Geophysical Prospecting，2021，56(1)：57–61.

[20] 刘财，李红星，陶春辉，等.模糊嵌套多级中值滤波方法及其在地震数据处理中的应用[J]. 地球物理学报，2007，50(5)：1534–1542.

LIU Cai，LI Hongxing，TAO Chunhui，et al. A new fuzzy nesting multilevel median filter and its application to seismic data processing[J]. Chinese Journal of Geophysics(in Chinese)，2007，50(5)：1534–1542.

[21] 刘伟，曹思远，王征，等. 随机噪声的多尺度多方向域压制方法研究[J]. 石油物探，2011，50(3)：301–306.

LIU Wei，CAO Siyuan，WANG Zheng，et al. Multi-scale & multi-direction suppressing method for random noise[J]. Geophysical Prospecting for Petroleum，2011，50(3)：301–306.

[22] 康冶，于承业，贾卧，等. f–x域去噪方法研究[J]. 石油地球物理勘探，2003，38(2)：136–138.

KANG Ye，YU Chengye，JIA Wo，et al. A study on noise- suppression method in f-x domain[J]. Oil Geophysical Prospecting，2003，38(2)：136–138.

[23] 田新琦，周彤，王志明，等. 滑动扫描可控震源地震数据谐波干扰的消除方法[J]. 石油物探，2011，50(6)：565–574.

TIAN Xinqi，ZHOU Tong，WANG Zhiming，et al. Method of eliminating harmonic noise for vibroseis data in slip sweep technique[J]. Geophysical Prospecting for Petroleum，2011，50(6)：565–574.

[24] 刘强，韩立国，李洪建. 混采数据分离中插值与去噪的同步处理[J]. 地球物理学报，2014，57(5)：1647–1654.

LIU Qiang，HAN Liguo，LI Hongjian. Synchronous interpolation and denoising in simultaneous-source data separation[J]. Chinese Journal Geophysics(in Chinese)，2014，57(5)：1647–1654.

[25] 佘德平. 多震源地震正演数值模拟技术[J]. 地球物理学进展，2012，27(4)：1533–1540.

SHE Deping. A numerical modeling of simultaneous sources[J]. Progress in Geophysics(in Chinese)，2012，27(4)：1533–1540.

[26] 朱立华. 基于稀疏反演的多震源地震混合采集数据分离技术[J]. 石油物探，2018，57(2)：208–212.

ZHU Lihua. Deblending simultaneous-source seismic data using sparse inversion[J]. Geophysical Prospecting for Petroleum，2018，57(2)：208–212.

[27] ZU Shaohuan，ZHOU Hui，LI Qingqing，et al. Shot-domain deblending using least-squares inversion[J]. Geophysics，2017，82(4)：1–96.

Blending interference noise separation method of simultaneous source seismic data based on multi-t-x domain combination

YU Dai1, HUANG Dezhi2, SUN Yuan1, YANG Wenping3

(1. College of Geology Engineering and Geomatics, Chang’an University, Xi’an 710054, China; 2. College of Mining Engineering, Heilongjiang University of Science & Technology, Harbin 150022, China; 3. North China Branch of SINOPEC Geophysical Corporation, Zhengzhou 450000, China)

With the increasing application of simultaneous source seismic data acquisition method in seismic oil and gas exploration, many scholars have studied the key technology of simultaneous source seismic data blending interference noise separation. So far, research on methods with higher fidelity is still needed. This article is based on the following assumptions: Normal Moveout Correction(NMO) in simultaneous source seismic data can enhance the linear distribution degree of reflection(especially in Common Middle Point(CMP) gather), and expand the distribution difference of reflection and blending interference noise; single domain separation of blending interference noise only aims at a certain distribution characteristic difference of blending interference noise and reflection; multi-t-x domain combined separation can comprehensively utilize the distribution difference of blending interference noise and reflection in each t-x domain, and its separation effect is better. We first use NMO for simultaneous source seismic data CMP gather to enhance the linearity of the reflection, and expand the dispersion of randomly distributed blending interference noise. Secondly, we use median filtering to separate most of the blending interference noise. Finally, based on the residual blending interference noise distribution in other t-x domains(such as Common Shot Gather(CSG), Common Receiver Gather(CRG) and Common Offset Gather(COG) etc.), the random noise attenuation is used for further separation, and the amplitude-preserving method is used in the whole separation process. The multi-t-x domain combination separation method is verified by theoretical data, and compared with the single domain separation method, it can effectively separate blending interference noise from simultaneous source seismic data with high fidelity, and suppress other noises and shear waves to a certain extent. The actual data results show that the method has better suppression effect and better seismic profile results.

simultaneous source seismic acquisition; blending interference noise separation; multi-t-x domain combination; noise discretization

移动阅读

语音讲解

P631

A

1001-1986(2021)05-0209-11

2021-05-31；

2021-08-06

陕西省自然科学基金项目(2021-JQ-588)

俞岱，1986年生，男，陕西西安人，博士，从事井中地震资料处理和偏移成像研究工作. E-mail：tablemax47@126.com

俞岱，黄德智，孙渊，等. 基于多t-x域联合的多震源地震数据混叠噪声分离方法[J]. 煤田地质与勘探，2021，49(5)：209–219. doi: 10.3969/j.issn.1001-1986.2021.05.023

YU Dai，HUAN Dezhi，SUN Yuan，et al. Blending interference noise separation method of simultaneous source seismic data based on multi-t-x domain combination[J]. Coal Geology & Exploration，2021，49(5)：209–219. doi: 10.3969/j. issn.1001-1986.2021.05.023

(责任编辑 聂爱兰)