Sedimentary environmental evolution for the past 30 ka of the northern continental slope of the South China Sea
-
摘要: 南海是西太地区最大的边缘海，汇集了周边陆地大量碎屑物质。这些陆源碎屑通过复杂的洋流系统经陆坡大量的输送到南海深海海盆中，使陆坡成为研究深海沉积物源汇体系不可缺少的重要环节。但陆坡区域水深变化大，洋流体系复杂，加之冰期间冰期海平面升降和季风的变化，使陆坡沉积环境一直成为研究的难点。为了研究陆坡沉积环境的演变过程，本文选用了南海北部陆坡中部和底部的两个重力柱开展元素地球化学方面的研究，探讨陆坡区域近三万年以来的沉积环境特征。研究发现海平面和季风是影响区域沉积环境的两个重要因素：（1）海平面变化是控制陆坡陆源物质/深海钙质碎屑变化的主要因素；（2）研究区域地层发育有“碳酸盐稀释事件”与东亚夏季风在全新世初期（11.5～8.5 kaBP）增强有关。Abstract: The South China Sea is the largest marginal sea in the Western Pacific, which receives massive sediments from the surrounding landmasses. Terrestrial sediments are transported into the deep basin via a complex ocean current system, with the continental slope as a key component of the source-to-sink process. However, past changes in sedimentary environment of the continental slope remain rarely understood, because of the drastic variations in water depth and complicated current systems in addition to sea level changes induced by glacial-interglacial cycles and the changes in trade wind. In this research, two gravity cores collected from the northern South China Sea are used to study the sedimentary environmental evolution of the continental slope. One is located in the middle of the slope, and the other at the lower part of the slope. According to the results of element geochemistry, it is revealed that: 1) the sea level change is the key factor which controls the changes in the ratio of terrigenous to biogenic components of the sediments; 2) a carbonate dilution event is found related to the intensification of the East Asian Summer Monsoon during the early Holocene from 11.5 to 8.5 kaBP.
-
1. 研究背景
在海洋资料中，多次波干扰非常发育并且种类也较多，有海水的鸣震、强海底尤其是崎岖海底产生的海底相关多次波、强反射界面产生的层间和长周期多次波等，这些多次波会造成地震记录中有效反射能量被压制，信噪比降低。因此，多次波的压制一直是海洋地震数据处理中的难点问题，也是海上资料处理的主要任务[1]。
深水海域地震资料数据处理是深水油气勘探的重要环节，其中多次波的压制又是重中之重，它直接影响到地震资料的品质，因此在偏移之前，尽可能地压制或衰减多次波。在深水海域，存在的多次波主要是自由表面多次波，该类多次波定义为地下介质反射的地震波到达自由表面后，至少发生一次下行反射，然后经一定传播路径后重新返回自由表面所接收的地震波[2-3]。可以说，在深水海域，如果能够压制自由表面多次波，也就压制了大部分的多次波干扰，因此自由表面多次波的压制是整个多次波压制的重点。针对此类多次波，学者们提出了很多压制的方法，有CMP叠加、f-k滤波法、Radon变换、聚束滤波法、预测反褶积和基于波动理论的多次波预测相减法等，其中目前最为广泛应用的是广义自由表面多次波预测技术（General-Surface Multiple Prediction，GSMP），相比于传统的二维自由表面多次波压制技术（Surface-Related Multiple Elimination，SRME），该技术预测的多次波模型更准确。同时，海上二维采集过程中电缆中—远偏移距难免受海流影响而偏离设计测线方向形成羽角，这是海上二维地震资料采集的固有特点。羽角的存在使共反射点发散无法满足SRME技术对规则化采集的要求，从而影响后续的多次波预测。因此，在本次多次波压制中，我们采用的是GSMP技术，但是在印度洋深水海域，海底相关多次波能量强，频带宽，常规的GSMP技术也不能得到很好的压制，因此，本文利用曲波变换，将多次波模型进一步优化，得到更加精确的多次波模型，从而使多次波的压制效果更好[4-9]。
2. 方法原理
2.1 广义自由表面多次波预测技术
广义自由表面多次波预测技术是近几年来逐渐兴起并广泛应用于海洋地震资料数据处理中的一项新技术。在理论上，该技术可以预测并衰减所有与地表相关的多次波，并且无需地下任何的先验信息，如速度、地层和构造等信息，是基于数据驱动的。广义自由表面多次波预测是通过模型建立和自适应减去法实现的，具体的实现途径为波动方程建模法，是在地表一致性褶积法的基础上进行改进的，通过波动方程外推来实现对多次波的模拟，该技术能适应任意观测系统，并且不受炮检点位置的约束。具体过程如下：首先对单炮数据进行时间反转，然后再向下外推，并与海底的反射系数进行褶积，再做向上的外推处理，最后完成整个单炮的多次波建模[6-8]。
2.2 曲波域多次波模型优化
广义自由表面多次波预测产生多次波模型，然后将地震数据和模型数据转换到曲波域，对多次波模型进一步优化，最后利用原始数据与多次波模型相减，对多次波进行压制。曲波变换使用的是第二代曲波变换，解决了第一代曲波变换大量数据冗余的问题，使曲波变换的实现更简单，运算效率更高。第二代曲波变换的公式为
$$ {\rm{c}}\left( {j,k,l} \right) = \left\langle {f,{{\rm{\varphi }}_{j,k,l}}} \right\rangle = \mathop \int \nolimits_{{R^2}}^{} f\left( x \right)\overline {{{\rm{\varphi }}_{j,k,{{l}}}}\left( x \right)} {\rm{d}}x $$ 其中，f（x）表示输入的原始地震信号或者多次波模型数据；φj,k,l为曲波函数，c（j,k,l）为曲波系数，其中j为尺度，l为方向，k为尺度j在l方向上的矩阵系数[10-13]。
具体的模型优化流程见图1，将地震数据和广义自由表面多次波预测产生的模型数据分为两部分，一部分是低频数据，一部分是高频数据，其中低频数据利用常规自适应减的方法得到低频多次波模型；高频数据动校后转换到曲波域，在曲波域中，比较不同尺度、不同角度的信号与多次波的振幅和相位差异（图2），具体的做法是：当信号与多次波的模型比较大于门槛值时，认为是信号，小于门槛值时，认为是多次波，依次来优化高频多次波模型，从而得到更加精确的多次波模型，再进行反动校（图3），最后用地震数据减去多次波模型，达到压制多次波的目的[14-17]。分高低频的主要原因是，在曲波域中，低频部分无法分角度和尺度对数据进行比较，见图4（分三个尺度）中Scale1，对低频模型无法进行优化，因此低频数据采用常规的自适应减，在高频数据中采用曲波变换对模型进行优化。高低频分界点的选取要稍大于Scale1的频率，低于Scale2的频率。
![]()
图 4 曲波变换示意图ω为频率，KN为空间奈奎斯特频率，N为尺度。Figure 4. Schematic diagram of curvelet transformω is the frequency, KN is the space Nyquist frequency, N is the scale.3. 实例分析
选取印度洋某深水海域的地震资料，该地区海底地形总体较为平坦，最大水深为5258 m。从原始炮集（图5）上可以看出，多次波主要是海底相关的多次波，图6是有效波与多次波频谱图的对比，其中红色是有效波频谱图，蓝色是多次波的频谱图，从图中可以看出，多次波能量强，频带宽，与有效波频谱基本一致。首先利用常规的广义自由表面多次波压制方法对其压制，图7是利用广义自由表面多次波压制方法得到的多次波模型，图8是压制后的炮集，可以看出多次波压制不干净，仍有较多残留。图9是利用本文方法，分4个尺度进行曲波变换，计算Scale1的频率为15.75 Hz，因此本文将原始数据和模型数据以20 Hz为界分为高频数据和低频数据，低频数据利用常规的自适应减的方法优化低频多次波模型，高频数据转到曲波域，在曲波域中根据不同尺度不同角度的信号与多次波的振幅和相位差异来优化高频多次波模型，然后将低频模型和高频模型相加得到优化后的多次波模型。为了更清晰地比较优化前后的多次波模型，将原始炮集的多次波与优化前后的多次波模型放大并进行比较，图10可以明显地看出，由浅至深，优化后的多次波模型与原始炮集的多次波更吻合，多次波模型的精确度更高。最后利用原始数据直接减去多次波模型，得到压制后的炮集，可以看出压制后炮集更干净，信噪比更高（图11）[18-21]。
![]()
图 6 有效波与多次波频谱图对比红色是有效波频谱，蓝色是多次波频谱。Figure 6. The spectrum of the effective wave compared with that of the multiple wavewhere red is the spectrum of the effective wave and blue is the spectrum of the multiple wave.![]()
图 10 多次波模型对比图从左到右依次为：原始数据多次波，常规方法得到的多次波模型，曲波域优化后的多次波模型。Figure 10. Multiples model comparison chartFrom left to right: multiples of raw data, multiples model obtained by conventional method, multiples model obtained by curvelet transform.![]()
图 11 利用优化后模型多次波压制效果Figure 11. Suppression of multiple waves using the optimized model下面从叠加剖面上看常规方法和本文方法的压制效果。选取印度洋该深水海域两条测线，图12是A测线原始剖面，图13是利用常规方法压制后的效果，可以看出压制效果不理想，多次波残留较为严重（图中箭头所指的地方）；图14 是利用本文方法压制后的效果，可以看出，压制效果较好，多次波去除的较为干净，剖面信噪比高，并且未损害有效信号，时间10.2 s的位置波组特征更加清晰，有利于后期地震资料的偏移和解释[22-25]。图15—17是B测线的原始剖面及利用常规方法和本文方法压制后的效果图，同样可以看出，利用本文方法压制多次波的效果更好，压制后的剖面信噪比更高，说明本文方法更适用于深水海域海底相关多次波的压制。
![]()
图 13 A测线常规方法压制后的叠加剖面Figure 13. The superimposed profile after conventional method of line A![]()
图 14 A测线利用曲波域优化模型压制的叠加剖面Figure 14. The stacked profile after optimization model in curved wave domain of line A![]()
图 17 B测线利用曲波域优化模型压制后的叠加剖面Figure 17. The stacked profile after optimization model in curved wave domain of line B![]()
图 16 B测线常规方法压制后的叠加剖面Figure 16. The superimposed profile after conventional method of line B4. 结论
本文通过在实际资料中的应用可以看出，多次波的压制效果较好，剖面的信噪比得到了较大的提高，同时压制后有效信号得到了凸显，波组特征更加清晰，有利于后期层位的识别和追踪。
该技术适用于海底地形较为平坦的深水海域，同时值得注意的是，本文方法在曲波域中对高频模型进行优化时，是根据信号和模型数据在不同尺度、不同角度上的振幅和相位差异，即当信号与多次波的模型比大于门槛值时，认为是信号，小于门槛值时，认为是多次波，因此门槛值的选择非常重要，直接决定优化后模型的精确度。门槛值的选择是选取有代表性的炮集，计算不同尺度、不同角度的振幅和相位差异，从而确定门槛值。
-
![]()
图 1 研究区域及周边地形图
a. 研究区域地形图，b.研究站位的三维地形图示，c. 研究站位地形剖面。图中地形数据根据参考文献[6]绘制，洋流数据根据参考文献[10-12]绘制。
Figure 1. Topographic map of the studied area
a. Topographic map of the area; b. three dimensional map of the studied region; c. topographic profile of studied sites; topographic data is derived from Ref.[6], current distribution citied from ref.[10-12].
![]()
![]()
![]()
![]()
图 5 SCSF39站位和SCSF41站位地球化学元素比值分布特征
Figure 5. Geochemical characteristics of the Core SCSF39 and SCSF41
![]()
图 6 SCSF39站位和SCSF41站位粒度分布特征
Figure 6. Grain size distributions of the Core SCSF39 and SCSF41
表 1 SCSF39站位和SCSF41站位基本信息
Table 1 Details of Core SCSF39 and SCSF 41
位置 水深/m 获取岩芯长度/cm 岩芯年龄/kaBP 平均沉积速/ （cm/ka） 经度/E 纬度/N SCSF39 114.97° 19.41° 1494 420 36.1 11.6 SCSF41 115.29° 18.61° 3717 460 36.7 12.5 
下载: 导出CSV
-
[1] Milliman J D, Farnsworth K L. River Discharge to the Coastal Ocean: A Global Synthesis[M]. Cambridge: Cambridge University Press, 2011: 1-382.
[2] Liu Z F, Zhao Y L, Colin C, et al. Source-to-sink transport processes of fluvial sediments in the South China Sea [J]. Earth-Science Reviews, 2016, 153: 238-273. doi: 10.1016/j.earscirev.2015.08.005
[3] Zhang Y W, Liu Z F, Zhao Y L, et al. Mesoscale eddies transport deep-sea sediments [J]. Scientific Reports, 2015, 4: 5937. doi: 10.1038/srep05937
[4] 舒业强, 王强, 俎婷婷. 南海北部陆架陆坡流系研究进展[J]. 中国科学: 地球科学, 2018, 61(5):560-571. [SHU Yeqiang, WANG Qiang, ZU Tinging. Progress on shelf and slope circulation in the northern South China Sea [J]. Science China Earth Sciences, 2018, 61(5): 560-571. doi: 10.1007/s11430-017-9152-y [5] 赵绍华, 刘志飞, 陈全, 等. 南海北部末次冰期以来深水沉积物组成及其堆积速率的时空变化特征[J]. 中国科学: 地球科学, 2017, 60(7):1368-1381. [ZHAO Shaohua, LIU Zhifei, CHEN Quan, et al. Spatiotemporal variations of deep-sea sediment components and their fluxes since the last glaciation in the northern South China Sea [J]. Science China Earth Sciences, 2017, 60(7): 1368-1381. doi: 10.1007/s11430-016-9058-6 [6] Becker J J, Sandwell D T, Smith W H, et al. Global bathymetry and elevation data at 30 arc seconds resolution: SRTM30_PLUS [J]. Marine Geodesy, 2009, 32(4): 355-371. doi: 10.1080/01490410903297766
[7] 姚伯初. 南海海盆新生代的构造演化史[J]. 海洋地质与第四纪地质, 1996, 16(2):1-13. [YAO Bochu. Tectonic evolution of the South China Sea in cenozoic [J]. Marine Geology & Quaternary Geology, 1996, 16(2): 1-13. [8] Wang P X, Li Q Y. The South China Sea: Paleoceanography and Sedimentology[M]. Netherlands: Springer, 2009: 1-497.
[9] 刘丛舒, 丁巍伟, 殷绍如, 等. 南海北部陆坡区海底峡谷地貌、沉积特征及控制因素[J]. 海洋学研究, 2019, 37(2):28-43. [LIU Congshu, DING Weiwei, YIN Shaoru, et al. Geomorphology, sedimentary characteristics and controlling factors of submarine canyons in the northern continental slope of the South China Sea [J]. Journal of Marine Sciences, 2019, 37(2): 28-43. [10] 俞慕耕, 刘金芳. 南海海流系统与环流形势[J]. 海洋预报, 1993, 10(2):13-17. [YU Mugeng, LIU Jinfang. Ocean current system and ocean current situation in the South China Sea [J]. Marine Forecasts, 1993, 10(2): 13-17. [11] 苏纪兰. 南海环流动力机制研究综述[J]. 海洋学报, 2005, 27(6):1-8. [SU Jilan. Overview of the South China Sea circulation and its dynamics [J]. Acta Oceanologica Sinica, 2005, 27(6): 1-8. [12] Liu Z Q, Gan J P. Three-dimensional pathways of water masses in the South China Sea: A modeling study [J]. Journal of Geophysical Research: Oceans, 2017, 122(7): 6039-6054. doi: 10.1002/2016JC012511
[13] Huang E Q, Tian J, Qiao P J, et al. Early interglacial carbonate-dilution events in the South China Sea: Implications for strengthened typhoon activities over subtropical East Asia [J]. Quaternary Science Reviews, 2015, 125: 61-77. doi: 10.1016/j.quascirev.2015.08.007
[14] Ge Q, Chu F Y, Xue Z, et al. Paleoenvironmental records from the northern South China Sea since the Last Glacial Maximum [J]. Acta Oceanologica Sinica, 2010, 29(3): 46-62. doi: 10.1007/s13131-010-0036-9
[15] Wang L J, Sarnthein M, Grootes P M, et al. Millennial reoccurrence of century-scale abrupt events of East Asian Monsoon: A possible heat conveyor for the global deglaciation [J]. Paleoceanography, 1999, 14(6): 725-731. doi: 10.1029/1999PA900028
[16] 范维佳, 陈荣华. 南海北部5万年来的表层海水盐度及东亚季风降水[J]. 第四纪研究, 2011, 31(2):227-235. [FAN Weijia, CHEN Ronghua. Sea Surface salinity and east Asian monsoon precipitation since the last 50000 years in the northern South China Sea [J]. Quaternary Sciences, 2011, 31(2): 227-235. doi: 10.3969/j.issn.1001-7410.2011.02.04 [17] Wang L, Sarnthein M, Erlenkeuser H, et al. East Asian monsoon climate during the Late Pleistocene: high-resolution sediment records from the South China Sea [J]. Marine Geology, 1999, 156(1-4): 245-284. doi: 10.1016/S0025-3227(98)00182-0
[18] 郑洪波, 陈国成, 谢昕, 等. 南海晚第四纪陆源沉积: 粒度组成、动力控制及反映的东亚季风演化[J]. 第四纪研究, 2008, 28(3):414-424. [ZHENG Hongbo, CHEN Guocheng, XIE Xin, et al. Grain size distribution and dynamic control of late Quaternary terrigenous sediments in the South China Sea and their implication for East Asian monsoon evolution [J]. Quaternary Sciences, 2008, 28(3): 414-424. doi: 10.3321/j.issn:1001-7410.2008.03.005 [19] Liu Z F, Li X J, Colin C, et al. A high-resolution clay mineralogical record in the northern South China Sea since the Last Glacial Maximum, and its time series provenance analysis [J]. Chinese Science Bulletin, 2010, 55(35): 4058-4068. doi: 10.1007/s11434-010-4149-5
[20] Folk R L, Andrews P B, Lewis D W. Detrital sedimentary rock classification and nomenclature for use in New Zealand [J]. New Zealand Journal of Geology & Geophysics, 1970, 13(4): 937-968.
[21] 汪品先. 十五万年来的南海: 南海晚第四纪古海洋学研究阶段报告[M]. 上海: 同济大学出版社, 1995: 1-172. WANG Pinxian. The Past 150 kyr of the South China Sea: A Periodic Report on the Late Quaternay Paleoceangraphy Research of the South China Sea[M]. Shanghai: Tongji University Press, 1995: 1-172.
[22] Zhao S H, Liu Z F, Colin C, et al. Responses of the East Asian Summer Monsoon in the low-latitude South China Sea to high-latitude millennial-scale climatic changes during the last glaciation: Evidence from a high-resolution clay mineralogical record [J]. Paleoceanography and Paleoclimatology, 2018, 33(7): 745-765. doi: 10.1029/2017PA003235
[23] 石学法, 刘升发, 乔淑卿, 等. 东海闽浙沿岸泥质区沉积特征与古环境记录[J]. 海洋地质与第四纪地质, 2010, 30(4):19-30. [SHI Xuefa, LIU Shengfa, QIAO Shuqing, et al. Depositional features and palaeoenvironmental records of the mud deposits in Min-Zhe Coastal Mud Area, East China Sea [J]. Marine Geology & Quaternary Geology, 2010, 30(4): 19-30. [24] 刘升发, 石学法, 刘焱光, 等. 东海内陆架泥质区表层沉积物常量元素地球化学及其地质意义[J]. 海洋科学进展, 2010, 38(1):80-86. [LIU Shengfa, SHI Xuefa, LIU Yanguang, et al. Geochemical characteristics and geological significance of major elements in the surface sediments from the inner shelf mud area of the East China Sea [J]. Advances in Marine Science, 2010, 38(1): 80-86. doi: 10.3969/j.issn.1671-6647.2010.01.011 [25] 肖尚斌, 陈木宏, 陆钧, 等. 南海北部陆架柱状沉积物记录的残留沉积[J]. 海洋地质与第四纪地质, 2006, 26(3):1-5. [XIAO Shangbin, CHEN Muhong, LU Jun, et al. New evidence for remnant deposits recorded by columnar sediments in the shelf of the Northern South China Sea [J]. Marine Geology & Quaternary Geology, 2006, 26(3): 1-5. [26] 田景春, 张翔. 沉积地球化学[M]. 北京: 地质出版社, 2016: 1-196. TIAN Jingchun, ZHANG Xiang. Sedimentary Geochemistry[M]. Beijing: Geological Publishing House, 2016: 1-196.
[27] Waelbroeck C, Labeyrie L, Michel E, et al. Sea-level and deep water temperature changes derived from benthic foraminifera isotopic records [J]. Quaternary Science Reviews, 2002, 21(1-3): 295-305. doi: 10.1016/S0277-3791(01)00101-9
-
期刊类型引用(5)
1. 刘鸿，徐华宁，刘欣欣，陈江欣，张菲菲，王小杰，颜中辉，杨佳佳，杨睿. 海洋地球物理数据处理现状及展望. 海洋地质与第四纪地质. 2024(03): 40-52 . 
本站查看
2. 王小杰，刘欣欣，颜中辉，刘鸿，杨佳佳. 基于曲波域模型优化的多次波压制方法在浅地层剖面的应用. 石油物探. 2024(06): 1155-1162 . 百度学术
3. 颜中辉，杨睿，冯京，刘欣欣，刘鸿，王小杰，姜春涛. 高分辨率小道距处理技术及属性分析在水合物识别中的应用. 海洋地质与第四纪地质. 2024(06): 46-59 . 本站查看
4. 冯晟，郑勋，王沛，宋峙潮. 基于预测反褶积理论的多次波压制技术研究. 黑龙江科学. 2023(04): 136-138 . 百度学术
5. 李森，梁杰，龚建明，廖晶，陈建文，张银国，王建强，杨传胜，袁勇，雷宝华. 巴基斯坦东部海域中-新生代沉积研究进展. 海洋地质前沿. 2022(02): 1-13 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 2158
- HTML全文浏览量: 270
- PDF下载量: 143
- 被引次数: 6
