地震数据POCS插值重建中的一种新型质量控制准则

王本锋,陈小宏,耿建华

(1.同济大学海洋地质国家重点实验室,海洋与地球科学学院,高等研究院,上海200092;2.中国石油大学(北京)油气资源与探测国家重点实验室,北京102249)

因障碍物、禁采区及海洋拖缆羽状漂移等各种因素的影响,采集的地震数据沿空间方向不规则,该不规则性会影响自由表面多次波去除、AVO分析、反演及偏移等处理结果的精度,因此有必要对采集到的地震数据进行插值重建[1]。插值重建方法可以大致归纳为4类[2-4]:基于稀疏变换的插值重建方法[5];基于预测滤波的插值重建方法[6];基于波动方程的插值方法[7]和基于降秩的插值方法[8-10]。

随着压缩感知理论的问世[11],基于稀疏变换的插值重建方法,特别是基于稀疏Curvelet变换的插值方法得到广泛关注[12-13]。高精度的凸集投影(POCS)方法[2-3,14-15]是基于稀疏变换的重要插值方法之一,该方法常用于图像重构[16],1991年,MENKE[17]首次将其引入地球物理领域的插值重建方面。目前,国内专家学者研究了基于压缩感知策略的地震数据插值重建问题[18-21]。大多数插值重建方法属于迭代方法的范畴[12-13,22-23],因此建立一个合理的质量控制准则(QCC)成为高效地震数据插值重建的关键,质量控制准则可在保证插值精度的同时,适时地结束迭代,达到提高插值重建计算效率的目的。GAO等[24]给出了两种不同的质量控制准则来监测插值迭代的收敛过程:第一种方法利用了完整的数据体,仅适用于插值重建理论方法研究;第二种方法利用了相邻的迭代更新解,但容易陷入局部极值,影响插值重建的计算精度[25]。

本文基于迭代硬阈值方法及投影算子推导了精确高效POCS方法,在此基础上,提出了一种新的质量控制准则,该质量控制准则仅利用了观测数据的信息,具有更好的适应性。基于稀疏Curvelet变换及POCS插值重建方法,利用两套模拟数据和实际数据验证了该质量控制准则的有效性及计算的高效性。

1 方法原理

1.1 POCS插值方法

空间方向不规则的观测地震数据dobs可以近似表征为:

(1)

式中:R为对角型采样算子;d0为需重建的完整数据体。由于受到观测地震数据频带有限及数据缺失等因素的影响,方程(1)的求解是不适定的,因此基于压缩感知理论,利用稀疏变换及稀疏约束构建目标泛函为:

(2)

式中:x为Curvelet系数向量;CT为Curvelet逆变换(C为Curvelet变换);λ为正则化因子,用于权衡Curvelet系数稀疏性与拟合残差。方程(2)可以由迭代硬阈值方法进行迭代求解[26-27]:

(3)

式中:xk为Curvelet域内第k次的更新解;Tλk为硬阈值函数;λk为阈值[3]。地震数据插值重建的目的是得到数据域重建结果,因此将变换域的更新结果转换到数据域得:

(4)

插值迭代收敛之前,在观测位置处dk+1的值偏离观测数据dobs,为了提高插值重建的效率,将dk+1投影到观测平面{d|dobs=Rd}上,得:

(5)

(6)

实际上POCS插值方法是一个两步法:稀疏变换域内更新解以及数据域内置入观测数据。该方法具有较高的插值重建精度及计算效率。

1.2 新型质量控制准则

由于POCS插值方法是一种迭代方法,有必要建立一个合理的质量控制准则来适时地终止迭代,进而提高插值重建的计算效率。GAO等[24]给出了两种质量控制准则:

(9)

图1 地震数据位置集合分解示意

2 数值实验

首先,利用两套模拟数据验证本文提出的质量控制准则的可行性及有效性;其次,利用一套实际资料进一步验证本文提出的质量控制准则能够在保持插值精度的同时,提高插值计算效率的有效性。

2.1 模拟算例

图2a为层状介质模型对应的模拟数据体(模拟数据1),包含201道,每道1001个采样点,空间采样率与时间采样率分别为12.5m和2ms。图2b为利用改进的Jittered欠采样方法进行50%地震道缺失得到的不规则采样数据。采用基于稀疏Curvelet变换的POCS插值方法进行插值重建,根据经验,最大迭代次数设置为50次,利用不同的质量控制准则J1,J2,J3,并进行归一化,得到的收敛曲线如图3所示。J1为理论的质量控制准则,可以作为参考解,合理且逐渐下降至收敛,其收敛性可以保证插值方法具有较高的精度。J2为利用相邻迭代更新解定义的质量控制准则,较大地偏离参考解,可能使插值迭代陷入局部极小值,进而影响插值重建结果的精度。然而,J3利用观测位置Ω1处的信息定义质量控制准则,其与参考解J1具有较好的一致性,具有更广的实际应用前景。图3验证了质量控制准则J3的可行性与有效性,即当J3满足预先设定的允许误差时,插值迭代可以适时地终止,在保证插值精度的同时提高插值计算效率。

图4a为利用基于稀疏Curvelet变换的POCS插值重建方法得到的插值结果,可以看出,其与原始数据具有较高的一致性。图4b为插值结果与原始数据之间的重构误差,可见,图上除了有微弱的边界效应外,残差较小。该边界效应是由边界处数据的不规则性产生,可以通过扩边的方式降低这种边界效应,但会增加计算量,具体量化分析留待后续研究。图3和图4验证了基于稀疏变换的POCS插值方法的有效性,且提出的质量控制准则与参考值的一致性较好。

图2 完整数据体(a)与不规则采样数据体(b)

图3 不同质量控制准则对应的归一化收敛曲线(模拟数据1)

为了进一步验证本文提出的质量控制准则的有效性,将其应用于含有干涉及弱反射的模拟数据(模拟数据2)中,图5a和图5b分别为完整数据体及不规则采样数据体,其中图5b缺失40%数据。利用不同的质量控制准则得到的归一化收敛曲线如图6所示,其中最大迭代次数为50次。图6与图3相似,均显示出本文提出的质量控制准则J3与理论解J1具有较好的一致性,但J2与理论解J1差别较大。图6进一步验证了本文提出的质量控制准则J3的可行性与有效性,其与理论值J1具有较好的一致性,同时也指出了J2的缺点。插值重建结果如图7a所示,其与原始完整数据体具有较好的一致性,且残差较小,可以忽略不计(图7b),插值重建结果与插值残差验证了基于稀疏变换的POCS插值方法的有效性。

上述两个模拟算例验证了本文提出的质量控制准则J3的有效性与可行性,其与理论值J1具有较好的一致性,但又与J1不同,J3仅利用了观测数据的信息,在实际数据插值重建中具有更好的灵活性,应用范围更广泛,在保证插值重建精度的同时,提高了插值重建计算效率。

2.2 实际数据算例

对不规则实际数据进行插值重建,理论质量控制准则J1不适用,因为J1利用了完整数据体的信息,但是本文提出的质量控制准则J3是一个较好的选择,而且与J1具有相似的性态。对于不规则缺失实际数据插值重建,下面仅比较质量控制准则J2和J3的应用效果,进一步验证质量控制准则J3的优势。

图8a为利用改进的Jittered欠采样方法得到的不规则实际数据,数据的不规则性会影响自由表面多次波去除、偏移与反演的精度。为了插值重建完整的数据体,利用基于稀疏Curvelet变换的POCS插值方法得到的重构数据体如图8b所示,数据插值中最大迭代次数为50。

图4 插值重建结果(a)与重建误差(b)

图5 完整数据体(a)与40%缺失不规则采样数据体(b)

图6 不同质量控制准则对应的归一化收敛曲线(模拟数据2)

由图8b可看出,插值重建后地震数据的横向连续性得到提高,该重建的地震数据可以作为后续地震数据处理的完整数据体。质量控制准则J2和J3对应的归一化收敛曲线如图8c所示,可以看出,J3对应的收敛曲线缓慢下降,具有一定的合理性,然而J2可能会使插值迭代陷入局部极值。为了利用J3进一步提高插值重建的计算效率,当满足预先设定的容许残差(J3<0.005)时,插值迭代被适时终止,插值结果如图8d所示。此时,迭代次数为30次(如图8c中纵向虚线所示),可以看出,相对50次迭代而言,计算效率提高了40%,但没有明显的降低插值重建精度。实际上,图8b与图8d具有较好的一致性,从而验证了本文提出的新型质量控制准则J3的有效性,在保证插值重建精度的同时,可以有效地提高插值重建的计算效率。从图8b和8d中还可以看出,插值结果中有方向各异的微弱假象,该假象由稀疏变换域系数截断引起,因此需研究更多的插值重建策略,进一步衰减假象、提高插值重建的精度。

图7 插值结果(a)与插值误差(b)

图8 实际数据插值重建结果a 不规则实际数据; b 迭代50次后插值重建结果; c 质量控制准则J2和J3对应的归一化收敛曲线; d 当J3<0.005时的插值重建结果

3 结论与认识

本文基于POCS插值重建方法推导过程,提出了一种新型质量控制准则J3来监测插值重建过程,可以适时地终止插值迭代过程,在保证插值重建精度的同时,进一步提高插值重建的计算效率。阐述并对比分析了不同类型的质量控制准则,提出的质量控制准则J3与理论值J1具有较好的一致性,而J3仅利用了观测数据的信息,在实际地震数据插值重建监测过程中具有更好的灵活性与适用性,应用范围更广。模拟与实际数据处理验证了新型质量控制准则J3的有效性。