METHODS FOR DELINEATING THE RESOURCE PERSPECTIVE AREA OF COBALT-RICH CRUSTS
-
摘要: 提出了4种不同的富钴结壳资源远景区圈定方法,即水深法、工程间距外推法、网格微分法和综合法,对比了它们的优缺点,对同一研究区运用不同的方法圈定了不同的资源远景区并进行了资源评价,结果表明:这几种方法是层层递进、相辅相成的,水深法和工程间距外推法可以给出研究区资源远景的大致范围,网格微分法是圈定较为精确资源远景区的基础,但是,加入地质规律的综合法可以提供更接近真实情况的资源远景区。不同的方法适合不同的勘探程度和不同研究目的,并且有不同的精度。Abstract: The paper brings forward four different methods for delineating the resource perspective areas of cobalt-rich crusts, and they are as follows:water depth method, exploration engineering distance extrapolation method,gridding differential coefficient method and comprehensive method. The advantage and disadvantage are also compared with each other. The different resource perspective areas of cobalt-rich crusts are delineated using different methods, and their resources are also evaluated. The results show that the four methods are related to one another, and the water depth method and exploration engineering distance extrapolation method can give the approximate variation range of resource future, and gridding differential coefficient method is the base to delineate better resource perspective area, but the comprehensive method considers the geological rules, so it can give one resource perspective area which is closer to reality. Each of the four methods fits its own investigate degree and research aim, and has its own precision.
-
随着复杂油气藏勘探精度的提高，精细的地震属性反演、油气储层的定量描述都需要高分辨率的地震资料。而地层的吸收衰减作用是影响地震资料分辨率的一个主要因素[1]，从而导致了有效反射波频带变窄，相位产生畸变，分辨率和信噪比降低，从而影响了地震资料的成像精度。Q补偿作为一种确定性的处理技术[2]，不仅可以补偿地震反射信号的频率和振幅衰减，同时能改善地震子波的相位特性，使中深层的高频部分得到加强，频谱被拓宽，从而有效提高地震资料的质量。
反Q滤波补偿方面也有许多学者进行了研究，基本可以分为3类：用级数展开作近似高频补偿的反Q滤波、基于波场延拓的反Q 滤波和其他反Q滤波方法。最早Hale提出Futterman模型进行反Q滤波的方法[3]，该方法基于预测误差滤波分析迭代求出Q值，然后进行Q补偿，其缺点是计算量大。Bickel和Natarajan根据复函数分析技术提出一种反Q滤波方法，此方法基于平面波假设，通过复数对平面波传播进行描述，从而得到时变的Q补偿[4]。Hargreaves和Calvert根据常Q模型得到相位Q补偿方法[5]，其做法是通过波场外推来计算常Q模型，从而对相位进行补偿，缺点是未考虑振幅的影响。Bano利用常Q模型发展到层常Q值模型的相位Q补偿方法，优势在于补偿结果能稳定地改善相位的畸变，而且准确性比较高，但同样忽略了振幅的影响[6]。WANG基于等间隔形式的Q模型提出一种稳定有效的反Q滤波方法，此方法以波场延拓理论为基础，并对Q值求取的抗噪性和稳定性进行深入讨论，并不断完善，形成了最常用的反Q滤波算法[7-9]。严红勇和刘洋在WANG方法的基础上，发展成沿射线路径进行波场延拓从而进行Q补偿的方法，主要是针对多分量资料进行处理[10]。刘财等提出基于频率域分时窗吸收补偿的反Q滤波方法[11]；Zhang和Ulrych在反Q滤波过程中引入了最小二乘方法和Bayasian原理，从而弥补了Q求取的不稳定性，并在实际数据中得到了很好的效果[12]；王本锋等结合反演的思想正则化方法提出了一种新的补偿方法，该方法基于波场延拓的正Q滤波方程，对有效频带内的频率分量进行分析，提高了计算效率，补偿结果稳定高效[13]。叶秋焱利用广义S变换求取时频域表层Q值的方法，通过表层Q值与井控深层Q值建立综合Q体，从而对地震数据进行Q补偿[14]。
虽然地层Q值补偿方法的理论和算法比较成熟，但是在实际地震资料处理应用中一直存在其稳定性和抗噪性不足的问题[15-19]，尤其是地层Q值的准确性。本文在HHT方法的基础上，基于EMD自适应分解和HHT时频谱局部刻画能力[20-23]，通过将地震资料分成不同的频率段，在HHT时频域计算不同目标层的地层Q值，从而进行Q值补偿。确保了Q值提取精度和Q值补偿的准确性。
1. 基于HHT的地层Q值提取
HHT方法的核心内容是固有模态分解（EMD）[24-25]，其主要思路是将具有复杂尺度成分的非平稳、非线性信号分解成一组较优Hilbert变换性能的固有模态函数（IMF）之和。我们给定复杂信号X(t)，经过EMD自适应分解可表示为：
$$ X\left( t \right) = \sum\nolimits_{i = 1}^n {{c_i}\left( t \right) + {r_n}\left( t \right)} $$ (1) 式中，ci(t)是固有模态函数，rn(t)称为余量。EMD分解其优势在于信号分解过程是自适应的。得到的各个不同尺度IMF都是平稳窄带信号，从而使得波形轮廓对称更好。实际上，EMD方法主要是从频率特征尺度进行分解，第一步得到的是信号中频率特征尺度最小的IMF分量；第二步得到频率特征时间尺度较大的IMF分量；最后分离出来的IMF具有频率特征时间尺度最大的特征。
原信号X(t)分解成多个不同尺度IMF后，基于复数道构建技术，可以构成解析信号Z(t)，
$$ Z\left( t \right) = X\left( t \right) + iY\left( t \right) $$ (2) $$ Z\left( t \right)\text{可写为:}\;\;\;\;\;\;Z\left( t \right) = a\left( t \right){{\rm{e}}^{i\theta \left( t \right)}} $$ (3) Y(t)为X(t)经过Hilbert变换的信号，每一个IMF分量经过Hilbert变换后，构造解析函数可以表示为：
$$ H\left( t \right) = \sum\nolimits_{j + 1}^n {{a_j}\left( t \right)\exp \left[ {i\int {{w_j}\left( t \right){\rm{d}}t} } \right]} $$ (4) 其中，aj(t)为分解的第j阶IMF分量，wj(t)为特定的瞬时频率，对照上式看出，HHT变换实际是傅里叶变换的扩展，优势在于容易反映局部振幅和频率的信息。去除残余函数，展开后得到Hilbert幅值谱，记作：
$$ H\left( {w,t} \right) = Re\sum\nolimits_{j + 1}^n {{a_j}\left( t \right){{\rm{e}}^{i\int {{w_j}\left( t \right){\rm{d}}t} }}} $$ (5) 基于HHT分频技术的地层Q值提取方法如下：首先基于HHT时频分析技术对目的层进行分析，选择适合目的层的优势频率范围；然后利用EMD分解技术对地震数据根据频率尺度进行分解，选择只含有目的层优势频率范围的地震数据重构；最后在HHT时频域内通过谱比法计算地层Q值[26]。本文方法在实际地震数据对地层Q补偿应用中的具体步骤为：将实际地震数据分成浅层、中深层，分别对这两套数据进行HHT时频分析，提取局部的频率属性特征，选择适合不同层的频率范围，基于EMD分解技术对这两套数据进行分解、选择、重构处理[27]，得到具有不同频率尺度的地震资料；基于复数道技术，在HHT时频域内通过谱比法计算Q值。此做法的优势在于不同频率成分地震信号估算的Q值避免了造成相互影响，同时局部的Q值异常求得更加准确，补偿效果更佳。HHT时频域内谱比法是假定地震子波为零相位子波，其估算地层Q 值的公式为：
$$ \ln \frac{{W\left( {a,{t_1}} \right)}}{{W\left( {a,t} \right)}} = - \alpha \left( {\frac{{{t_1} - t}}{Q}} \right) $$ (6) 其中，
$\alpha = {\rm{\pi }}\left( {\dfrac{{{f_1} + \tau f}}{{\tau + 1}}} \right)$ 式中，τ表示能量衰减率，为一常数，f1表示零相位子波的视频率；f表示平面波频率；t和t1表示平面波从地面开始传播到目的层上下界面所用的时间。W(a，t)表示HHT时频谱。通过公式（6）进行线性回归可以求得地层Q值[28]。
2. 地层Q值补偿
地震波在地下传播时，由于信号的中心频率要向较低的频率方向移动，因此找出Q和中心频率的关系至关重要。以此通过Q补偿可以消除地震波在地下介质中传播的频散吸收，提高地震数据的垂向分辨率[29-31]。
在均匀介质中传播的平面波频域表达式为：
$$ U\left( f \right) = {U_0}\left( f \right)\exp \left[ {{{i}}\left( {2{\rm{\pi }}ft - kx} \right)} \right] $$ (7) 式中，U0(f)表示震源脉冲响应；k为波数；x为波前面沿着射线路径从震源传播到检波器的距离。假如介质是粘弹性，则k为复数，表达式如下：
$$ k\left( f \right) = \frac{{2{\rm{\pi }}f}}{{c\left( f \right)}} = \frac{{2{\rm{\pi }}f}}{{v\left( f \right)}} - i\alpha \left( f \right) $$ (8) 其中，c(f)为复速度，v(f)为相速度，α(f)为衰减系数。
在Kolsky-Futterman模型中，相速度和衰减系数对应的关系满足：
$$ \alpha \left( f \right) = \frac{{\rm{\pi }}}{{v\left( f \right)Q}}f $$ (9) $$ \frac{1}{{v\left( f \right)}} = \frac{1}{{{v_r}}}\left( {1 - \frac{1}{{{\rm{\pi }}Q}}\ln \left| {\frac{f}{{{f_r}}}} \right|} \right) $$ (10) 式中，fr代表参考频率；vr代表参考相速度。
于是将公式（8）至公式（10）代入公式（7）中，展开后对应表达式：
$$ U\left( f \right) = {U_0}\left( f \right)\exp \left[ { - \left( {\frac{{{\rm{\pi }}ft}}{Q} + {{i}}\frac{{2ft}}{Q}\ln \left| {\frac{f}{{{f_r}}}} \right|} \right)} \right] $$ (11) 式中，
$t = \dfrac{x}{{{v_r}}}$ 。定义
$G\left( {f,t} \right) = \exp \left[ { - \left( {\dfrac{{{\rm{\pi }}ft}}{Q} + {{i}}\dfrac{{2ft}}{Q}\ln \left| {\dfrac{f}{{{f_r}}}} \right|} \right)} \right]$ 大地滤波算子，则反Q滤波算子为：$$ {G^{ - 1}}\left( {f,t} \right) = A\left( {f,t} \right)P\left( {f,t} \right) = {\rm{exp}}\left( {\frac{{{\rm{\pi }}ft}}{Q} + {{i}}\frac{{2ft}}{Q}{\rm{ln}}\left| {\frac{f}{{{f_r}}}} \right|} \right) $$ (12) 式中，
$P\left( {f,t} \right) = {\rm{exp}}\left( {i\dfrac{{2ft}}{Q}{\rm{ln}}\left| {\dfrac{f}{{{f_r}}}} \right|} \right)$ 是相位补偿算子，$A\left( {f,t} \right) = {\rm{exp}}\dfrac{{{\rm{\pi }}ft}}{Q}$ 是振幅补偿算子。由此推导出地震波衰减补偿在时频域的表达式：
$$ {U_0}\left( {f,t} \right) = U\left( {f,t} \right)A\left( {f,t} \right)P\left( {f,t} \right) $$ (13) 最后进行HHT反变换，得到补偿后的时间域地震信号。
3. 模型试验
基于上述理论基础，首先进行模型参数的试验，通过层状介质的衰减模拟分析，以及Q值求取与补偿研究对理论进行验证。模型参数如下：对应层的速度依次为1 500、2 000、3 000 m/s，Q值依次为120、60、120。正演模拟采用基于单程波动方程的分步傅里叶波场延拓方法，采用的地震子波为雷克子波，频率为30 Hz，时间采样间隔1 ms。图1a为模拟得到的地震记录，图2a为放大的波形显示。由于地下介质吸收衰减作用，体现出振幅减小和波形畸变的现象，从放大的波形分析，子波波形不再是零相位，即分辨率逐渐降低。
![]()
图 1 模型数据Q补偿前后的对比a. 补偿前的地震数据，b. 补偿后的地震数据。Figure 1. Comparison of model seismic data before and after Q compensationa. Model data before compensation，b. Model data after compensation.![]()
图 2 模型数据Q补偿前后的波形显示对比a. 补偿前的波形显示，b. 补偿后的波形显示。Figure 2. Waveform comparison of model seismic data before and after Q compensationa. Waveform before compensation，b. Waveform after compensation.图3是基于小波分频技术提取的地层Q值曲线，与模型参数基本一致。图1b、图2b是经过本文方法进行补偿后的结果。对比图1a和图2a可以看出，通过Q补偿处理的不同深度地层反射波形均恢复到了原始的雷克子波形态，对应每一层地震信号的振幅和相位都得到了恢复。图4是经过Q补偿前后的振幅谱对比，蓝色为补偿前的频谱，红色为补偿后的频谱。由图可以看出，补偿后的频谱得到展宽。由此可以说明，通过本文方法对地层进行Q值补偿，能够达到使地震子波波形恢复原状，主频提高，频带变宽，提高分辨率的作用。
![]()
图 4 模型数据Q补偿前后的频谱对比Figure 4. Spectrum comparison before and after model data Q compensation为了验证本文Q补偿方法的抗噪性，通过对模型数据加入8%的随机噪声，同样的步骤对加噪的数据进行Q补偿测试。图5为补偿前后的数据对比，图6为补偿前后的频谱对比。从图中的对比分析，补偿后的地震记录有效信息得到加强恢复，噪声的部分基本没变。也就是说本文方法在信号补偿的同时，也能确保信噪比。
![]()
图 5 加噪的模型数据Q补偿前后的对比a. 补偿前的地震数据，b. 补偿后的地震数据。Figure 5. Comparison of model seismic data with noise before and after Q compensationa. Model data before compensation，b. Model data after compensation.![]()
图 6 加噪的模型数据Q补偿前后的频谱对比Figure 6. Spectrum comparison of model seismic data with noise before and after Q compensation4. 实际数据应用
本文选择某海域地震资料实际数据对本文方法进行应用效果分析。选取某一段信噪比较好的资料，对其浅层、中深层进行Hilbert分频处理，选择适合不同层的频率区间在HHT时频域计算地层的Q值，进而进行反Q滤波，实现对地层振幅和相位的补偿。图7对应的是根据本文方法求得的Q值剖面，图8a是原始的叠后剖面，从图中可看出Q值剖面跟地层吻合较好。图8b是整个剖面进行补偿后的效果，可以看出分辨率有了显著提高，同时剖面的信噪比也能得到很好的保持。图9为浅层经过补偿后的对比图，补偿后同相轴的连续性有了明显的改善。图10是中深层补偿后的对比图，同样的，补偿后的剖面分辨率更高，地层信息更加丰富。
![]()
图 8 实际地震数据Q补偿结果对比a. 补偿前的地震数据，b. 补偿后的地震数据。Figure 8. Comparison of actual seismic data before and after Q compensationa. Seismic data before compensation, b. Seismic data after compensation.![]()
图 9 浅层数据Q补偿结果对比a. 补偿前的地震数据，b. 补偿后的地震数据。Figure 9. Comparison of shallow data before and after Q compensationa. Seismic data before compensation，b. Seismic data after compensation.![]()
图 10 中深层数据Q补偿结果对比a. 补偿前的地震数据，b. 补偿后的地震数据。Figure 10. Comparison of deep data before and after Q compensation resultsa.Seismic data before compensation，b. Seismic data after compensation.图11是整个补偿前后振幅谱对比图，图12 为浅层补偿前后频谱对比图。对比可知，Q补偿后的数据频谱有了较大的改善，从浅层到深层表现出高频能量提高，数据主频范围向高频移动，且频谱的形态基本没发生变换，说明补偿是有效的，同时频带变宽，对应的剖面高频成分也有很明显的补偿效果，分辨率有所提高。尤其对中深层的补偿，从图13分析可知，原始资料的高频部信息缺失严重，通过Q补偿后，高频的有效信息被恢复。由于实际数据为海上资料，补偿后高频部分振幅能量得到加强，相应的会导致海上虚反射引起的陷波点效应加剧，海上常规资料处理的方法是通过鬼波压制解决。
![]()
图 12 浅层数据Q补偿前后频谱对比Figure 12. Spectrum comparison of shallow data before and after Q compensation![]()
图 13 中深层数据 Q补偿前后频谱对比Figure 13. Spectrum compensation before and after Medium-deep data Q-compensation为了验证本文方法在提高分辨率的同时，不会损失信噪比，提取噪音部分进行频谱分析，如图14所示，Q补偿后的频谱基本没变换，也就是说噪音的成分没有加强。从补偿的剖面来看，存在信噪比降低的表象，其实只是整个资料主频向高频移动的表现。
![]()
图 14 噪音成分Q补偿前后频谱对比Figure 14. Spectrum comparison before and after noise component Q-compensation5. 结论
（1）HHT变换具有很高的局部分析能力，因此，基于HHT时频谱的地层Q值提取方法能得到更加准确的地层Q值，同时EMD分解算法是基于频率尺度进行分解，可针对不同目标区频率区间进行地震数据的Q补偿，从而达到提高资料分辨率的目的。
（2）通过本文方法进行Q值补偿，在有效提高地震资料主频，展宽频带，提高分辨率的同时，确保了资料的信噪比。从模拟数据到实际地震数据的测试表明本文方法是可行且有效的。
（3）由于EMD分解的过程存在模态混叠效应，同时Q值提取的正确性对资料信噪比有一定的要求。因此，在处理之前应对地震数据进行去噪预处理，以达到更佳的补偿效果。
-
[1] James Hein. Cobalt-rich ferromanganese crusts:global distribution,composition,origin and research activities[J].International Seabed Authority,2003:36-89.
[2] 武光海,马维林,刘捷红,等.海山富钴结壳的圈矿方法[J].海洋学研究,2005,23(4):15-19. [WU Guang-hai,MA Wei-lin,LIU Jie-hong,et al. A method for finding the ore block of Co-rich crusts on seamounts[J].Joural of Marine Science, 2005,23(4):15-19.]
[3] 张国祯.富钴结壳矿区申请需要解决的两个重要问题[J].海洋地质动态,2005,21(1):19-22. [ZHANG Guo-zhen. Two important questions to solve before China applying for survey of cobalt-rich crust in mining areas[J].Marine Geology Development, 2005,21(1):19-22.]
[4] 初凤友,孙国胜,李晓敏,等.中太平洋海山富钴结壳生长习性及控制因素[J].吉林大学学报,2005,25(3):320-325. [CHU Feng-you,SUN Guo-sheng,LI Xiao-min,et al.The growth habit and controlling factors of the cobalt-rich crusts in seamount of Central Pacific[J].Journal of Jilin University:Earth Science Edition, 2005,25(3):320-325.]
[5] 曹瑜,胡光道,杨志峰,等. GIS环境下地质变量自动提取与地质异常的圈定[J]. 计算机工程与应用,2003(14):81-85.[CAO Yu,HU Guang-dao,YANG Zhi-feng,et al.Geo-variables extraction and geo-anomaly circle based on GIS[J].Computer Engineering and Application, 2003 (14):81-85.]
[6] 周林立,胡光道.富钴结壳资源量评估方法[J].海洋地质动态,2005,21(6):29-31. [ZHOU Lin-li,HU Guang-dao. Methods to evaluate cobalt-rich crust resources[J].Marine Gelogody Letters, 2005,21(6):29-31.]
[7] 马维林,初凤友,金翔龙.富钴结壳资源评价和圈矿方法探讨[J].海洋学报,2007,29(2):67-73. [MA Wei-lin,CHU Feng-you,JIN Xiang-long. Method approach of resource assessment and ore delineation for cobalt-rich crust[J].Acta Oceanologica Sinica, 2007,29(2):67-73.]
[8] 何高文,梁东红,宋成兵,等.浅地层剖面测量和海底摄像联合应用确定平顶海山富钴结壳分布界线[J].地球科学,2005,30(4):509-512. [HE Gao-wen,LIANG Dong-hong,SONG Cheng-bing,et al.Deterring the distribution boundary of cobalt-rich crusts of guyot by synchronous application of sub-bottom profiling and deep-sea video recording[J]. Earth-Science-Journal of China University of Geosciences, 2005, 30(4):509-512.]
-
期刊类型引用(1)
1. 刘鸿，徐华宁，刘欣欣，陈江欣，张菲菲，王小杰，颜中辉，杨佳佳，杨睿. 海洋地球物理数据处理现状及展望. 海洋地质与第四纪地质. 2024(03): 40-52 . 
本站查看
其他类型引用(2)
下载:
计量
- 文章访问数: 1607
- HTML全文浏览量: 114
- PDF下载量: 8
- 被引次数: 3
