一种基于亚像素匹配的拖缆动态采集影响消除方法技术

一种基于亚像素匹配的拖缆动态采集影响消除方法，包括以下步骤：１、向计算机输入所需参数；２、逐条测线进行校正：将该测线的数据抽成Ｙ个共偏移距道集ｇ；叠加得到叠加剖面ｍ；对每个共偏移距道集ｇ横向上作道内插，时间方向上作线性插值，求叠加剖面ｍ中每个数据集中附近数据的最大互相关，并将其存为一个数据集ｄ；再对ｄ进行中值滤波；３、处理完所有共偏移距道集后，将数据抽成共成像点道集保存，即完成拖缆动态采集影响的消除。本发明专利技术采用的方法在共偏移距道集上进行匹配，从而在亚像素范围内校正数据，克服拖缆漂移带来的影响，能够使共成像点道集上的同相轴更平直，远近道的数据一致性更好，并扩展叠加剖面有效频带的范围，提高分辨率。

【技术实现步骤摘要】

本专利技术属于海洋地震信号处理领域，特别涉及一种消除拖缆动态采集影响的方 法。
技术介绍
拖缆采集是海上地震资料采集的主要方式。在采集时，航行在海面上的采集船拖 拽多条具有检波器的拖缆。采集船上装备有空气枪，并通过空气枪产生地震波，拖缆上检波 器完成地震数据的采集工作。在理想情况下，拖缆应与采集船航行的速度方向平行。但是， 由于海洋中浪涌、潮汐等自然现象的作用，使得拖缆偏离了原定位置，并成为与航行方向成 一定夹角的曲线。因此，位于拖缆上的检波点也偏移了原定位置，并较以前更加靠近震源， 且随着偏移距的增加，检波点与原定位置的差距也随着增加，这就造成通过叠加得到叠加 剖面时，原本应处在同一位置的各偏移距对应的成像点，会随偏移距增加而更靠近震源的 位置。因此在叠加后的剖面上，原本反映地层中一个点的数据，将变成来自于这个点附近的 一个区域，从而降低了叠加剖面的分辨率。另一方面，由于潮汐的作用，海平面的高度在多 次采集的过程中是随时间变化的，从而使拖缆采集得到的数据在时间方向发生扰动，这种 扰动也会破坏多次采集所得到的数据的一致性，降低叠加剖面分辨率。
技术实现思路
针对上述问题，本专利技术的目的是提供一种基于亚像素匹配的拖缆动态采集影响消 除方法。该方法在地震数据处理的过程中，能够以数据为驱动有效校正克服拖缆漂移带来 的影响，提高数据的可叠加性，使共成像点道集上同相轴的一致性更好，并抬升叠加剖面的 高频成分，使其具有更高的分辨率。为实现上述目的，本专利技术采取的技术方案是一种基于亚像素匹配的拖缆动态采 集影响消除方法，该方法包括以下步骤步骤(1)，向计算机输入以下设定量测线数L、每条测线上的检波点数X、炮数Y、采样点数T、横向插值大小Sx、纵向插 值大小St、匹配时所用窗口大小Wx和Wt ；步骤(2)，用计算机逐测线进行校正，对于每条测线，执行如下步骤步骤(2. 1)，将该测线的数据抽成共偏移距道集g ；步骤(2. 2)，将共偏移距道集g叠加，得到叠加剖面m ；步骤(2.3)，按下面步骤逐个处理每个共偏移距道集g 步骤(2.3. 1)，对当前的共偏移距道集g在横向上作道内插，时间方向上作线性插 值，插值后的共偏移距道集数据大小变为SxXXStT ；步骤(2. 3. 2)，对插值后的共偏移距道集进行等间距采样，从而获得N个大小为 XXT 的数据集 gi，i = 1，…，N，N = SxXSt ；步骤(2. 3. 3)，对应原共偏移距道集上任意一个点(X，t)，求叠加剖面m中的相同位置点附近数据和每个数据集中(X，t)附近数据的互相关Cn，即 从N个Cn中找出最大的一个，并将其编号为i，存为一个数U=X-W V-t-W/据集d;步骤(2. 3. 4)，上一步的结果得到一个大小为XXT的数据集d，其中记录的内容为 最大互相关所对应的编号i，再对d进行中值滤波以减小噪声；步骤(2. 3. 5)，用i = d(x,t)所对应编号的数据集中该位置处的数据gi(x，t)替 换原共偏移距道集中数据g ；步骤(2. 4)，按步骤(2. 3. 1) 步骤(2. 3. 5)处理完一条测线上的所有共偏移距道 集后，再将数据抽成共成像点道集，并保存。步骤(3)，按步骤(2)各子步骤处理完所有测线后，即完成拖缆动态采集影响的消 除。由上述采用的技术方案可以看出，本专利技术在数据处理阶段，以采集得到的数据为 驱动，完成数据校正和动态采集影响的消除，不对采集过程提出额外的要求，因此，方法具 有较强的适应性，对于由多种拖缆采集方式采集的地震数据，均可以用本方法进行校正。附图说明图1是本专利技术方法的流程示意图；图2a是以第50个共成像点道集为例，200-1600ms，处理前共成像点道集图；图2b是以第50个共成像点道集为例，200-1600ms，处理后共成像点道集图；图2c是以第50个共成像点道集为例，1600-3000ms,处理前共成像点道集图；图2d是以第50个共成像点道集为例，1600-3000ms,处理后共成像点道集图；图3a是处理前的第20个共成像点道集的局部放大图；图3b是处理后的第20个共成像点道集的局部放大图；图4a是叠加剖面处理前示意图；图4b是叠加剖面处理后示意图；图5a是叠加剖面处理前局部放大示意图；图5b是叠加剖面处理后局部放大示意图；图5c是经过110Hz-120Hz滤波的叠加剖面处理前局部放大示意图；图5d是经过110Hz-120Hz滤波的局部叠加剖面处理后局部放大示意图；图6是某叠加道的振幅谱。具体实施例方式针对上述方案，具体来说，对于一个实际海上采集的地震数据，本专利技术的方法是逐 测线进行处理的。将测线上的所有共偏移距道集叠加得到叠加剖面，并将叠加剖面作为匹 配时所用的模板。每条测线对应一个模板，当处理下条测线时，应重新叠加产生新的模板。对每个共偏移距道集插值时，为了保证同相轴的连续性，需要采用道内插技术进 行水平方向的插值，即插值出新的地震道。有多种道内插算法可供使用，如f-χ域等道距内 插法、f-k域内插法以及倾角叠加法等。我们采用斯来蒙.司贝兹(Simon Spitz)提出的5f_x域的内插方法进行插值。该方法能产生2的整数次幂道数的数据。利用线性同相轴在 f"X域的可预测性，在f_x域中应用最小平方原理求取横向预测算子，继而通过预测算子与 内插算子之间的对偶关系，求得内插算子，然后应用最小平方原理求取内插道。为了减少系统内存的使用量并提高计算速度，我们采用划分重叠区域的方法对共 偏移距道集进行处理。若划分区域的大小为LxXLt,则第一个区域的范围是，那么水平方向上下一个区域的范围是，时间方向上下一个区 域的范围是，以此类推。因此，同一个点会被分配到4个不同的区 域中，在每个区域中分别得到4个不同的值后，进行加权求和得到最终的值。其中加权系数 为2维汉宁窗系数。本方法是在亚像素范围内对数据进行匹配和校正。在进行匹配时查找最大相关系 数的搜索范围不宜过大，否则会受到其他同相轴的干扰。在叠加剖面m中，我们一般在对应 叠加剖面m中与每个数据集相同位置点附近10个点范围内进行搜索，用这10个采样点计 算互相关值Cn。下面是对本专利技术实施例的一个具体实验验证我们采集的实际数据为一条测线，其中包含201个共成像点道集，每个道集有117 道，每道1501个采样点，采样率为2ms。我们采用的插值大小为Sx = St = 4，求互相关时采 用窗口大小为wx = 8, wt = 3。选取第50个共成像点道集和第20个共像点道集分析处理结果。图2a d给出了 第50个共成像点道集处理前后的比较结果，图3a和图3b是第20个共成像点处理前后的 局部放大结果，从图中的对比区域可以看出，经过校正后的共成像点道集中成像点扰动减 小了，同相轴的同相性得到明显改善。将所有共偏移距道集叠加后获得叠加剖面，图4a、图 4b是处理前后的叠加剖面图，处理后的叠加剖面细节更丰富，同相轴的连续性得以提高。图 5a d是对处理前后的叠加剖面选取了同一个局部区域进行超高频滤波的结果，图5a、5b 是对图4a、4b进行局部放大的结果，图5c、5d是对图5a、5b分别进行相同的超高频处理后 的结果，我们采用的是110Hz-120Hz的高通滤波。从图5c、5d中可以看出本文档来自技高网...

【技术保护点】
一种基于亚像素匹配的拖缆动态采集影响消除方法，其包括以下步骤：步骤（１），向计算机输入以下设定量：测线数Ｌ、每条测线上的检波点数Ｘ、炮数Ｙ、采样点数Ｔ、横向插值大小Ｓ↓［ｘ］、纵向插值大小Ｓ↓［ｔ］、匹配时所用窗口大小Ｗ↓［ｘ］和Ｗ↓［ｔ］；步骤（２），用计算机逐测线进行校正，对于每条测线，执行如下步骤：步骤（２．１），将该测线的数据抽成共偏移距道集ｇ；步骤（２．２），将共偏移距道集ｇ叠加，得到叠加剖面ｍ；步骤（２．３），按下面步骤逐个处理每个共偏移距道集ｇ：步骤（２．３．１），对当前的共偏移距道集ｇ在横向上作道内插，时间方向上作线性插值，插值后的共偏移距道集数据大小变为Ｓ↓［ｘ］Ｘ×Ｓ↓［ｔ］Ｔ；步骤（２．３．２），对插值后的共偏移距道集进行等间距采样，从而获得Ｎ个大小为Ｘ×Ｔ的数据集ｇ↓［ｉ］，ｉ＝１，…，Ｎ，Ｎ＝Ｓ↓［ｘ］×Ｓ↓［ｔ］；步骤（２．３．３），对应原共偏移距道集上任意一个点（ｘ，ｔ），求叠加剖面ｍ中的相同位置点附近数据和每个数据集中（ｘ，ｔ）附近数据的互相关Ｃ↓［ｎ］，即Ｃ↓［ｉ］（ｘ，ｔ）＝＊＊ｇ↓［ｉ］（ｕ，ｖ）ｍ（ｕ，ｖ）；从Ｎ个Ｃ↓［ｎ］中找出最大的一个，并将其编号为ｉ，存为一个数据集ｄ；步骤（２．３．４），上一步的结果得到一个大小为Ｘ×Ｔ的数据集ｄ，其中记录的内容为最大互相关所对应的编号ｉ，再对ｄ进行中值滤波以减小噪声；步骤（２．３．５），用ｉ＝ｄ（ｘ，ｔ）所对应编号的数据集中该位置处的数据ｇ↓［ｉ］（ｘ，ｔ）替换原共偏移距道集中数据ｇ；步骤（２．４），按步骤（２．３．１）～步骤（２．３．５）处理完一条测线上的所有共偏移距道集后，再将数据抽成共成像点道集，并保存。步骤（３），按步骤（２）各子步骤处理完所有测线后，即完成拖缆动态采集影响的消除。...

【技术特征摘要】

【专利技术属性】
技术研发人员：赵伟，陆文凯，王小六，郝振江，张云鹏，刘永江，杨俊，朱振宇，薛东川，王季，
申请(专利权)人：中国海洋石油总公司，中海石油研究中心，清华大学，
类型：发明
国别省市：11[中国|北京]