留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

利用重力数据反演中国东海海域莫霍面深度

何慧优 方剑 陈铭 崔荣花

何慧优, 方剑, 陈铭, 崔荣花. 利用重力数据反演中国东海海域莫霍面深度[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 682-689. doi: 10.13203/j.whugis20170330
引用本文: 何慧优, 方剑, 陈铭, 崔荣花. 利用重力数据反演中国东海海域莫霍面深度[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 682-689. doi: 10.13203/j.whugis20170330
HE Huiyou, FANG Jian, CHEN Ming, CUI Ronghua. Moho Depth of the East China Sea Inversed Using Gravity Data[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 682-689. doi: 10.13203/j.whugis20170330
Citation: HE Huiyou, FANG Jian, CHEN Ming, CUI Ronghua. Moho Depth of the East China Sea Inversed Using Gravity Data[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 682-689. doi: 10.13203/j.whugis20170330

利用重力数据反演中国东海海域莫霍面深度

doi: 10.13203/j.whugis20170330
基金项目: 

国家重点研发计划 2016YFC0601101

国家自然科学基金 41621091

国家自然科学基金 41874096

中国科学院战略性先导科技专项 XDB18010304

详细信息
    作者简介:

    何慧优, 博士生, 主要从事地球重力场与磁场研究。hehuiyou@asch.whigg.ac.cn

  • 中图分类号: P312

Moho Depth of the East China Sea Inversed Using Gravity Data

Funds: 

The National Key Research and Development Program of China 2016YFC0601101

the National Natural Science Foundation of China 41621091

the National Natural Science Foundation of China 41874096

the Strategic Priority Research Program of Chinese Academy of Sciences XDB18010304

More Information
    Author Bio:

    HE Huiyou, PhD candidate, specializes in the earth?s gravity field and magnetic field. E-mail:hehuiyou@asch.whigg.ac.cn

  • 摘要: 莫霍面是地壳与上地幔的分界面,莫霍面的研究有助于认识地球深部构造及演化。在以往的研究中,由于重力反演方法、使用资料不同,中国东海海域莫霍面深度反演结果有较大不确定性。通过收集中国东海海域高精度地形和自由空气重力异常数据,计算了布格重力异常并分析了其构造特征。提取了中国福建省大田至中国台湾省花莲剖面和中国东海陆架盆地至菲律宾海盆剖面的重力异常数据,采用人机交互正演、反演拟合技术得到了这两条剖面的地壳密度与莫霍面深度。通过向上延拓的方法提取了莫霍面所引起的重力异常,在两条剖面反演结果和地震探测结果的约束下进行Parker-Oldenburg界面反演,得到了该地区的莫霍面深度。反演结果显示,中国东海海域莫霍面深度的变化范围为-34.4~-8.8 km,其结构呈两凹两凸的特征。
  • 图  1  中国东海海域地形及构造单元分布图

    Figure  1.  Terrain and Tectonic Map in East China Sea

    图  2  中国东海自由空气重力异常图

    Figure  2.  Free-Air Gravity Anomaly Map in East China Sea

    图  3  中国东海布格重力异常图

    Figure  3.  Bouguer Gravity Anomaly Map in East China Sea

    图  6  中国东海布格重力异常向上延拓不同高度的结果

    Figure  6.  Bouguer Gravity Anomalies of the Different Distances to Upward Continue in East China Sea

    图  7  中国东海海域莫霍面深度图

    Figure  7.  The Moho Depth Map in East China Sea

    图  8  中国东海海域莫霍面深度三维图(X方向为东向)

    Figure  8.  3D Map of Moho Depth in the East China Sea(the X Axis is the East)

    图  4  剖面1重力异常人机交互正反演拟合曲线与密度模型图

    Figure  4.  Fitting Curve of Gravity Anomaly and Density Model Map of the First Profile

    图  5  剖面2重力异常人机交互正反演拟合曲线与密度模型图

    Figure  5.  Fitting Curve of Gravity Anomaly and Density Model Map of the Second Profile

    表  1  莫霍面深度对比/km

    Table  1.   Comparison of the Moho Depths/km

    数据来源 中国浙闽隆起区 中国东海陆架盆地 中国钓鱼岛隆褶带 冲绳海槽盆地 琉球隆褶区 琉球海沟 菲律宾海盆
    文献[7] -31~-29 -28.5~-27 -27.5~-25 -25~-18
    文献[14] -34~-29 -29~-27 -27~-22 -22~-16 -25~-19 -18~-12 -13~-12
    文献[5] -28~-25 -28~-26 -26~-16 -32~-18
    文献[8] -33~-28 -24~-19 -30~-26 -26~-14 -18~-4
    本文 -34.4~-29 -28.5~-26 -26~-24 -24~-21 -26.4~-23 -22~-15 -14~-8.8
    注:文献[8]中,中国浙闽隆起区和中国东海陆架盆地作为一个整体区域进行研究
    下载: 导出CSV
  • [1] 蔡学林, 朱介寿, 曹家敏, 等.四川黑水-台湾花莲断面岩石圈与软流圈结构[J].成都理工大学学报(自然科学版), 2004, 31(5):441-451 doi:  10.3969/j.issn.1671-9727.2004.05.001

    Cai Xuelin, Zhu Jieshou, Cao Jiamin, et al. Structures of Lithosphere and Asthenosphere in Heishi of Sichuan-Hualian of Taiwan Profile [J]. Journal of Chengdu University of Technology (Science and Technology Edition), 2004, 31(5): 441-451 doi:  10.3969/j.issn.1671-9727.2004.05.001
    [2] 高德章, 赵金海, 薄玉玲, 等.东海重磁地震综合探测剖面研究[J].地球物理学报, 2004, 47(5): 853-861 doi:  10.3321/j.issn:0001-5733.2004.05.017

    Gao Dezhang, Zhao Jinhai, Bo Yuling, et al. A Profile Study of Gravitative-Magnetic and Seismic Comprehensive Survey in the East China Sea[J]. Chinese Journal of Geophysics, 2004, 47(5):853-861 doi:  10.3321/j.issn:0001-5733.2004.05.017
    [3] 江为为, 刘少华, 郝天珧, 等.应用重力资料估算东海冲绳海槽地壳厚度[J].地球物理学进展, 2002, 17(1): 35-41 doi:  10.3969/j.issn.1004-2903.2002.01.004

    Jiang Weiwei, Liu Shaohua, Hao Tianyao, et al. Using Gravity Data to Compute Crustal Thickness of East China Sea and Okinawa Trough[J]. Progress in Geophysics, 2002, 17(1): 35-41 doi:  10.3969/j.issn.1004-2903.2002.01.004
    [4] 韩波, 张训华, 裴建新, 等.东海及其邻域壳-幔结构与展布特征[J].地球物理学进展, 2007, 22(2):376-382 doi:  10.3969/j.issn.1004-2903.2007.02.007

    Han Bo, Zhang Xunhua, Pei Jianxin, et al. Characteristics of Crust-Mantle in East China Sea and Adjacent Regions[J]. Progress in Geophysics, 2007, 22(2):376-382 doi:  10.3969/j.issn.1004-2903.2007.02.007
    [5] 周志远, 高金耀, 吴招才, 等.东海莫霍面起伏与地壳减薄特征初步分析[J].海洋学研究, 2013, 31(1):16-25 doi:  10.3969/j.issn.1001-909X.2013.01.002

    Zhou Zhiyuan, Gao Jinyao, Wu Zhaocai, et al. Preliminary Analyses of the Characteristics of Moho Undulation and Crustal Thinning in East China Sea[J]. Journal of Marine Sciences, 2013, 31(1):16-25 doi:  10.3969/j.issn.1001-909X.2013.01.002
    [6] Wang T K, Lin S F, Liu C S, et al. Crustal Structure of the Southernmost Ryukyu Subduction Zone: OBS, MCS and Gravity Modelling[J]. Geophysical Journal International, 2004, 157(1): 147-163 doi:  10.1111/gji.2004.157.issue-1
    [7] 郝天珧, 徐亚, 胥颐, 等.对黄海-东海研究区深部结构的一些新认识[J].地球物理学报, 2006, 49(2):458-468 doi:  10.3321/j.issn:0001-5733.2006.02.019

    Hao Tianyao, Xu Ya, Xu Yi, et al. Some New Understandings on Deep Structure in Yellow Sea and East China Sea[J]. Chinese Journal of Geophysics, 2006, 49(2):458-468 doi:  10.3321/j.issn:0001-5733.2006.02.019
    [8] 胡珊.东亚大陆边缘俯冲带的重力异常特征及俯冲模式研究[D].青岛: 中国海洋大学, 2015 http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y2883929

    Hu Shan. The Study of Gravity Anomalies and Subduction Model in the East Asia Continental Margin Subduction Zones[D]. Qingdao: Ocean University of China, 2015 http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y2883929
    [9] Becker J J, Sandwell D T, Smith W H F, 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
    [10] 李家彪.东海区域地质[M].北京:海洋出版社, 2008

    Li Jiabiao. Regional Geology of the East China Sea[M]. Beijing: China Ocean Press, 2008
    [11] Sandwell D T, Müller R D, Smith W H F, et al. New Global Marine Gravity Model from CryoSat-2 and Jason-1 Reveals Buried Tectonic Structure[J]. Science, 2014, 346 (6205): 65-67 doi:  10.1126/science.1258213
    [12] 管志宁.地磁场与磁力勘探[M].北京:地质出版社, 2005

    Guan Zhining. Geomagnetic Field and Magnetic Exploration[M]. Beijing: Geological Publishing House, 2005
    [13] 曾华霖.重力场与重力勘探[M].北京:地质出版社, 2005

    Zeng Hualin. Gravity Field and Gravity Exploration[M]. Beijing: Geological Publishing House, 2005
    [14] 韩波.东海地球物理场及深部地质构造研究[D].北京: 中国科学院研究生院, 2008 http://cdmd.cnki.com.cn/Article/CDMD-80068-2008117235.htm

    Han Bo. Geophysical Field and Deep Tectonic Features of East China Sea[D]. Beijing: Graduate University of the Chinese Academy of Sciences, 2008 http://cdmd.cnki.com.cn/Article/CDMD-80068-2008117235.htm
    [15] 徐德琼, 李全兴, 蒋家祯.东海莫霍面及其地质意义[J].海洋通报, 1983, 2(5):35-41 http://www.cnki.com.cn/Article/CJFDTOTAL-HUTB198305006.htm

    Xu Deqiong, Li Quanxing, Jiang Jiazhen. The Mohorovicic Discontinuity in the East China Sea and Its Geologic Significance[J]. Marine Science Bulletin, 1983, 2(5):35-41 http://www.cnki.com.cn/Article/CJFDTOTAL-HUTB198305006.htm
    [16] 张力方, 徐杰, 彭艳菊, 等.东海地区新构造运动研究[J].地震地质, 2014, 36(3): 692-705 doi:  10.3969/j.issn.0253-4967.2014.03.012

    Zhang Lifang, Xu Jie, Peng Yanju, et al. A Study on Neotectonic Movement in the East China Sea[J]. Seismology and Geology, 2014, 36(3): 692-705 doi:  10.3969/j.issn.0253-4967.2014.03.012
  • [1] 袁华清, 吴云龙, 许闯, 张毅, 王振宇.  基于实测重力数据的秦岭地区均衡异常与垂向构造应力分析 . 武汉大学学报 ● 信息科学版, 2022, 47(6): 972-978. doi: 10.13203/j.whugis20220198
    [2] 王树良, 李大鹏, 赵柏翔, 耿晶, 张伟, 王海雷.  聊天机器人技术浅析 . 武汉大学学报 ● 信息科学版, 2021, 46(2): 296-302. doi: 10.13203/j.whugis20190177
    [3] 翟振和, 孙中苗, 肖云, 李芳.  自主海洋测高卫星串飞模式的设计与重力场反演精度分析 . 武汉大学学报 ● 信息科学版, 2018, 43(7): 1030-1035, 1128. doi: 10.13203/j.whugis20160234
    [4] 王伟, 李姗姗, 邢志斌, 李中全, 黄韵颖.  RAE并行滤波的重力异常匹配算法 . 武汉大学学报 ● 信息科学版, 2016, 41(12): 1664-1670. doi: 10.13203/j.whugis20140652
    [5] 刘金钊, 柳林涛, 梁星辉, 叶周润.  利用重力异常数据构建重力梯度场的方法比较 . 武汉大学学报 ● 信息科学版, 2015, 40(12): 1677-1682. doi: 10.13203/j.whugis20130797
    [6] 胡敏章, 李建成, 金涛勇.  顾及局部地形改正的GGM海底地形反演 . 武汉大学学报 ● 信息科学版, 2013, 38(1): 60-63.
    [7] 杜劲松, 陈超, 梁青, 张毅.  基于球冠滑动平均的球面重力异常分离方法 . 武汉大学学报 ● 信息科学版, 2012, 37(7): 864-868.
    [8] 胡敏章, 李建成, 金涛勇.  应用重力地质方法反演皇帝海山的海底地形 . 武汉大学学报 ● 信息科学版, 2012, 37(5): 610-612.
    [9] 李姗姗, 吴晓平, 田颜锋.  水下重力异常最近等值线迭代匹配算法的改进 . 武汉大学学报 ● 信息科学版, 2011, 36(2): 226-230.
    [10] 李振海, 罗志才, 汪海洪, 钟波.  给定内插重力异常精度时对重力数据的要求 . 武汉大学学报 ● 信息科学版, 2011, 36(11): 1328-1331.
    [11] 王瑞, 李厚朴.  基于逆Stokes公式的测高重力反演中央区效应计算 . 武汉大学学报 ● 信息科学版, 2010, 35(4): 467-471.
    [12] 王虎彪, 王勇, 陆洋, 周旭华.  联合多种测高数据确定中国海及其邻域1.5′×1.5′重力异常 . 武汉大学学报 ● 信息科学版, 2008, 33(12): 1292-1295.
    [13] 程芦颖, 许厚泽.  顾及测站点上重力场信息的大地水准面高的拟合方法 . 武汉大学学报 ● 信息科学版, 2008, 33(7): 701-705.
    [14] 杨元德, 鄂栋臣, 黄金维, 汪海洪.  Geosat/GM波形重跟踪反演中国沿海区域重力异常 . 武汉大学学报 ● 信息科学版, 2008, 33(12): 1288-1291.
    [15] 张赤军, 方剑, 边少锋, 高金跃.  重力均衡与地形波长——兼论均衡理论不适宜于珠峰重力值的推估 . 武汉大学学报 ● 信息科学版, 2007, 32(3): 229-233.
    [16] 黄谟涛, 王瑞, 翟国君, 欧阳永忠.  多代卫星测高数据联合平差及重力场反演 . 武汉大学学报 ● 信息科学版, 2007, 32(11): 988-993.
    [17] 史红岭, 陆洋, 王勇.  CHAMP重力卫星结果在中国海及邻近海域的初步分析 . 武汉大学学报 ● 信息科学版, 2003, 28(S1): 48-51.
    [18] 晁定波, 姚运生, 李建成, 徐菊生.  南海海盆测高重力异常特征及构造解释 . 武汉大学学报 ● 信息科学版, 2002, 27(4): 343-347.
    [19] 管泽霖, 鄂栋臣.  按克林索求和计算大地水准面差距垂线偏差及重力异常 . 武汉大学学报 ● 信息科学版, 1986, 11(4): 75-82.
    [20] 何慧优, 方剑.  重力异常频谱分析方法研究及其应用 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20200510
  • 加载中
图(8) / 表(1)
计量
  • 文章访问数:  854
  • HTML全文浏览量:  68
  • PDF下载量:  121
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-03-12
  • 刊出日期:  2019-05-05

利用重力数据反演中国东海海域莫霍面深度

doi: 10.13203/j.whugis20170330
    基金项目:

    国家重点研发计划 2016YFC0601101

    国家自然科学基金 41621091

    国家自然科学基金 41874096

    中国科学院战略性先导科技专项 XDB18010304

    作者简介:

    何慧优, 博士生, 主要从事地球重力场与磁场研究。hehuiyou@asch.whigg.ac.cn

  • 中图分类号: P312

摘要: 莫霍面是地壳与上地幔的分界面,莫霍面的研究有助于认识地球深部构造及演化。在以往的研究中,由于重力反演方法、使用资料不同,中国东海海域莫霍面深度反演结果有较大不确定性。通过收集中国东海海域高精度地形和自由空气重力异常数据,计算了布格重力异常并分析了其构造特征。提取了中国福建省大田至中国台湾省花莲剖面和中国东海陆架盆地至菲律宾海盆剖面的重力异常数据,采用人机交互正演、反演拟合技术得到了这两条剖面的地壳密度与莫霍面深度。通过向上延拓的方法提取了莫霍面所引起的重力异常,在两条剖面反演结果和地震探测结果的约束下进行Parker-Oldenburg界面反演,得到了该地区的莫霍面深度。反演结果显示,中国东海海域莫霍面深度的变化范围为-34.4~-8.8 km,其结构呈两凹两凸的特征。

English Abstract

何慧优, 方剑, 陈铭, 崔荣花. 利用重力数据反演中国东海海域莫霍面深度[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 682-689. doi: 10.13203/j.whugis20170330
引用本文: 何慧优, 方剑, 陈铭, 崔荣花. 利用重力数据反演中国东海海域莫霍面深度[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 682-689. doi: 10.13203/j.whugis20170330
HE Huiyou, FANG Jian, CHEN Ming, CUI Ronghua. Moho Depth of the East China Sea Inversed Using Gravity Data[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 682-689. doi: 10.13203/j.whugis20170330
Citation: HE Huiyou, FANG Jian, CHEN Ming, CUI Ronghua. Moho Depth of the East China Sea Inversed Using Gravity Data[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 682-689. doi: 10.13203/j.whugis20170330
  • 莫霍面是地壳与上地幔的分界面,上下界面物质的密度差异较大,是一个重要的密度界面。研究莫霍面结构特征对认识地球深部构造及地球动力学过程具有重要的意义。

    中国东海是西太平洋边缘海构造形成和演化的重要窗口。长期以来,关于中国东海海域莫霍面深度的研究一直是地球科学研究的热点问题之一。研究中国东海海域莫霍面深度对于该地区深部构造演化过程的认识和矿产资源的勘探具有重要意义。文献[1]运用现代构造解析理论和方法,对四川黑水-台湾花莲剖面地震数据进行了分析,结合地质学、地球化学及地球物理学研究得到了该剖面的岩石圈和软流圈结构模型。文献[2]采用重磁震综合反演解释系统,得到了中国东海陆架盆地至菲律宾海盆的重磁震综合探测剖面的地质解释模型。文献[3]利用重力资料,参考地震测深结果采用线性公式计算了冲绳海槽的莫霍面深度。文献[4]利用重力异常反演了中国东海及其邻域的莫霍面深度。文献[5]利用重力资料,结合地震约束反演得到了中国东海的莫霍面深度。文献[6]结合海底地震台与多道地震测量,通过建立密度模型进行重力数据拟合,得到了中国东海陆架盆地莫霍面深度。文献[7]采用小波变换与多尺度分析,利用调和级数法反演了中国东海海域莫霍面深度。文献[8]利用卫星重力资料和高精度水深数据,采用小波多尺度分析方法进行场分离,利用Parker-Oldenburg界面反演法计算了琉球海沟及邻区莫霍面深度。

    在以往的研究中,由于采用的反演方法、使用的资料不同,莫霍面深度反演结果有较大差异。本文收集了中国东海海域高精度地形和自由空气重力异常数据,以前人研究的中国福建省大田至中国台湾省花莲和中国东海陆架盆地至菲律宾海盆这两条剖面的地质模型为参考,进行重力异常人机交互正演、反演拟合(简称人机交互正反演拟合),得到了这两条剖面的地质-地球物理解释模型。采用向上延拓的方法提取了莫霍面引起的重力异常,以剖面反演得到的壳幔密度差和莫霍面深度及收集的地震资料揭示的莫霍面深度作为约束进行莫霍面界面反演,得到了中国东海海域莫霍面深度,并分析了莫霍面的结构特征。

    • 本文收集了中国东海海域高精度地形和自由空气重力异常数据,具体为:(1)地形数据,网格分辨率为15″×15″。该数据是文献[9]结合SRTM、ASTER和CryoSat-2冰盖地形得到的全球地形模型SRTM15。图 1是中国东海海域地形及构造单元分布图,显示了中国东海海域的构造单元从西向东依次为中国浙闽(浙江福建)隆起区、中国东海陆架盆地、中国钓鱼岛隆褶带、冲绳海槽盆地、琉球隆褶区、琉球海沟和菲律宾海盆[10]。(2)自由空气重力异常,网格分辨率为1′×1′,重力异常值精度为±2 mGal (见图 2)。该数据是文献[11]根据CryoSat-2和Jason-1雷达测高数据与现有数据相结合得到的高精度全球海洋重力模型。图 2显示,研究区域内的自由空气重力异常特征在陆地与海区存在较大差异。中国浙闽隆起区的自由空气重力异常在负异常背景上广泛分布着急剧升高的正异常圈闭,走向以北东、北西西向为主。中国东海陆架盆地自由空气重力异常变化平缓,走向基本为北东向。中国钓鱼岛隆褶带的自由空气重力异常主体走向为北东向,异常梯度较大。冲绳海槽盆地南部自由空气重力异常变化缓慢,北部则异常起伏较大。琉球隆褶区的自由空气重力异常存在高值正异常圈闭,北部异常呈北东向,南部异常呈近东西向。琉球海沟的自由空气重力异常存在规模大、幅值低的负异常带。菲律宾海盆的自由空气重力异常以变化平缓的正异常为主。本文中,图 1-3图 6-8为根据颜色等面积分布成图。

      图  1  中国东海海域地形及构造单元分布图

      Figure 1.  Terrain and Tectonic Map in East China Sea

      图  2  中国东海自由空气重力异常图

      Figure 2.  Free-Air Gravity Anomaly Map in East China Sea

      图  3  中国东海布格重力异常图

      Figure 3.  Bouguer Gravity Anomaly Map in East China Sea

      图  6  中国东海布格重力异常向上延拓不同高度的结果

      Figure 6.  Bouguer Gravity Anomalies of the Different Distances to Upward Continue in East China Sea

      图  7  中国东海海域莫霍面深度图

      Figure 7.  The Moho Depth Map in East China Sea

      图  8  中国东海海域莫霍面深度三维图(X方向为东向)

      Figure 8.  3D Map of Moho Depth in the East China Sea(the X Axis is the East)

      根据以上的地形和自由空气重力异常数据计算得到了布格重力异常(见图 3)。该异常主要反映了深部构造的演变趋势。从西北往东南异常幅值逐渐升高,说明莫霍面缓慢抬升。中国浙闽隆起区布格异常自西向东逐渐增大,由负异常逐渐变为正异常,负异常走向与海岸线一致,呈北东向,靠近海岸线与海区正异常连成一片,说明是大陆重力异常向海域的延伸。中国东海陆架盆地布格重力异常变化平缓,以正异常为主,呈北东向。中国钓鱼岛隆褶带布格重力异常呈北东向,中部异常值高。冲绳海槽布格重力异常呈北东向,异常值呈北部低、南部高的特点,表明地壳厚度从北向南逐渐减薄,莫霍面逐渐抬升。琉球隆褶区至菲律宾海盆,布格重力异常呈现高-低-高的趋势,反映了深部莫霍面先下降再抬升的变化。自由空气重力异常和布格重力异常与研究区域内的主要构造单元分布具有明显的相关性。

      图 3给出了本文用于参考的两条剖面的位置。其中,中国福建大田至中国台湾花莲剖面(剖面1)地质模型结合重力、地磁、地质、地震、大地热流等方面数据得到[1];中国东海陆架盆地至菲律宾海盆剖面(剖面2)地质模型依据详细的地震资料进行重磁震反演得到[2]。这两条剖面综合了多源数据,得到了较为详细的地层分布、莫霍面形态及物性特征,故将其作为本文人机交互正反演拟合的初始模型。

      在得到可靠的重力异常、初始剖面地质模型和地震揭示的莫霍面深度值后,采用频率域界面反演法对研究区域进行莫霍面反演,具体计算步骤为:

      1) 提取剖面1和剖面2的自由空气重力异常,将收集的两条剖面地质模型作为初始模型进行人机交互正反演拟合,建立地质-地球物理解释模型,得到剖面壳幔密度差和莫霍面深度。

      2) 通过对布格重力异常进行频率域向上延拓,分离出莫霍面引起的异常,在第1)步反演结果和地震资料的约束下进行界面反演,得到研究区域的莫霍面深度,并分析其特征。

      其中,频率域向上延拓方法和界面反演方法的原理分别如下。

      频率域向上延拓方法是假设场源位于z=H平面以下(H > 0),则重力场是z=H平面以上对xyz的连续函数[12]。若z=0,观测平面上的重力场Δg(x, y, 0)已知,则向上延拓公式为:

      $$\begin{array}{l} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\Delta g\left( {x, y, z} \right) = \\ \frac{{ - z}}{{2{\rm{ \mathsf{ π} }}}}\mathop {\mathop \smallint\limits_{ - \infty } }\limits^\infty \mathop {\mathop \smallint\limits_{ - \infty } }\limits^\infty \frac{{\Delta g\left( {\xi , \eta , 0} \right)}}{{{{\left[ {{{\left( {x - \xi } \right)}^2} + {{\left( {y - \eta } \right)}^2} + {z^2}} \right]}^{3/2}}}}{\rm{d}}\xi {\rm{d}}\eta \end{array}$$ (1)

      式中,ξη为积分变量。

      由褶积积分公式可知,式(1)即为Δg(x, y, 0)与$\frac{1}{{2{\rm{ \mathsf{ π} }}}} \cdot \frac{{ - z}}{{{{({x^2} + {y^2} + {z^2})}^{3/2}}}}$的褶积,空间域的褶积对应频率域的乘积。分别求Δg(x, y, 0)和$\frac{1}{{2{\rm{ \mathsf{ π} }}}} \cdot \frac{{ - z}}{{{{({x^2} + {y^2} + {z^2})}^{3/2}}}}$的傅里叶变换SΔg(u, v, 0)和e2π(u2+v2)1/2z。设Δg(x, y, z)的傅里叶变换为SΔg(u, v, z),频率域向上延拓的计算公式为:

      $${S_{\Delta g}}\left( {u, v, z} \right) = {S_{\Delta g}}\left( {u, v, 0} \right){{\rm{e}}^{2{\rm{ \mathsf{ π} }}{{({u^2} + {v^2})}^{1/2}} \cdot z}}$$ (2)

      式中,e2π(u2+v2)1/2z为频率域向上延拓滤波算子。将SΔg(u, v, z)进行傅里叶反变换即可得到向上延拓为z(z < 0)的重力场值。

      莫霍面界面反演则是采用Parker-Oldenburg界面反演方法[13],该方法采用快速傅里叶变换,具有高效的运算效率。

      假定在直角坐标系o-xz中,场源层的上部边界为z=0,下部边界为z=h(x),h(x)是相对于地面以下深度为z0的某个参考水平面度量的,定义其一维傅里叶变换为:

      $$F\left[ {h\left( x \right)} \right] = \mathop \smallint _{ - \infty }^\infty h\left( x \right){{\rm{e}}^{{\rm{i}}kx}}{\rm{d}}x$$ (3)

      式中,k是变换函数的波数。

      重力异常Δg(x)的一维傅里叶变换公式为:

      $$F\left[ {\Delta g\left( x \right)} \right] = - 2{\rm{ \mathsf{ π} }}G\sigma {{\rm{e}}^{ - \left| k \right|{z_0}}}\mathop {\mathop \sum\limits_{n = 1} }\limits^\infty \frac{{{{\left| k \right|}^{n - 1}}}}{{n!}}F\left[ {{h^n}\left( x \right)} \right]$$ (4)

      式中,σ为地壳与上地幔的界面密度差;G为重力常量;z0为参考面深度;n是变换点数。从式(4)的无限项中提出n=1,并重新排列得到:

      $$\begin{array}{l} F\left[ {h\left( x \right)} \right] = \frac{{F\left[ {\Delta g\left( x \right)} \right]{{\rm{e}}^{\left| k \right|{z_0}}}}}{{ - 2{\rm{ \mathsf{ π} }}G\sigma }} - \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\mathop {\mathop \sum\limits_{n = 2} }\limits^\infty \frac{{{{\left| k \right|}^{n - 1}}}}{{n!}}F\left[ {{h^n}\left( x \right)} \right] \end{array}$$ (5)

      假定σ已知,z0已知或给定,则可以应用式(5)进行计算。

    • 通过重力异常剖面反演来约束区域反演是减少区域反演多解性的有效方法。重力异常人机交互正反演拟合技术是重力异常剖面反演中常用的方法,它通过正演初始地质模型,对比重力异常的计算值与观测值,不断修正模型,直至拟合精度满足设定条件的技术要求。具体步骤为:分别提取剖面1和剖面2的自由空气重力异常作为观测值;在重磁正反演系统中输入初始地质模型,结合地震资料修正模型(地层位置和密度)来拟合模型计算值与观测值,直至拟合精度达到要求。其中,剖面1和剖面2的初始地质模型分别来自于文献[1]和文献[2]。模型的密度参数参考文献[14]中收集的中国东海多口钻井的岩芯密度值所总结的中国东海地壳分层密度参数。其中,海水密度为1.03 g/cm3,沉积层密度为2.2~2.4 g/cm3,侵入岩密度为2.68~2.85 g/cm3,上地壳密度为2.67 g/cm3,中地壳密度为2.75 g/cm3,下地壳密度为2.9 g/cm3,上地幔密度为3.3 g/cm3。由此可见,研究区域内下地壳和莫霍面上下物质存在较大的密度差异,是主要的密度界面。

      剖面1的反演结果如图 4所示。观测重力异常为本文收集的自由空气重力异常,计算重力异常为模型的重力正演值;模型中标注的数值为地层密度(下文同此),壳幔密度差为0.4 g/cm3。中国浙闽隆起区莫霍面埋深由西向东逐渐抬升,由-32 km慢慢抬升至-29 km。中国台湾海峡莫霍面深度为-29~-26 km。中国台湾地区莫霍面深度起伏较大,变化范围为-32~-22 km。沿剖面1向东至菲律宾海盆莫霍面深度抬升至-21 km,推测莫霍面的起伏变化与该地区菲律宾海板块向欧亚板块俯冲形成的构造有关。

      图  4  剖面1重力异常人机交互正反演拟合曲线与密度模型图

      Figure 4.  Fitting Curve of Gravity Anomaly and Density Model Map of the First Profile

      剖面2的反演结果如图 5所示,壳幔密度差为0.4 g/cm3。中国东海陆架盆地莫霍面深度为-28~-27 km,向东至中国钓鱼岛隆褶带莫霍面呈缓慢抬升的趋势。冲绳海槽盆地莫霍面抬升至-22.3 km。琉球隆褶区莫霍面慢慢下降至-25 km,向东至菲律宾海盆,莫霍面又抬升至-10 km。

      图  5  剖面2重力异常人机交互正反演拟合曲线与密度模型图

      Figure 5.  Fitting Curve of Gravity Anomaly and Density Model Map of the Second Profile

    • Parker-Oldenburg界面反演方法要求参与反演的异常是目标界面所产生的重力异常,故先从布格重力异常中提取莫霍面引起的重力异常。参考文献[14]将布格重力异常向上延拓多个高度的方法进行分析。向上延拓20 km的结果与延拓5 km和10 km的结果相比(图 6),清晰可见中国浙闽隆褶区、中国东海陆架盆地、中国钓鱼岛隆褶带、冲绳海槽盆地、琉球隆褶区、琉球海沟和菲律宾海盆等区域;向上延拓30 km的结果则损失了部分信息。故采用向上延拓20 km的布格重力异常作为莫霍面界面反演的重力异常较为合理。根据§2反演结果,取壳幔密度差为0.4 g/cm3,平均深度为-25 km进行Parker- Oldenburg界面反演,并绘制了中国东海海域莫霍面深度图(见图 7)。

      将本文反演结果、前人反演结果和地震探测的莫霍面深度值进行对比。文献[5]计算的莫霍面深度与地震探测的莫霍面深度值的均方差为2.4 km,文献[8]的均方差为1.81 km,而本文采用高精度的地形和自由空气重力异常数据反演得到的莫霍面深度值与地震探测值的均方差为1.37 km,吻合度更高。将本文反演的莫霍面深度与前人反演的结果(见表 1)进行对比,结果显示,在研究区域内的各个构造单元,本文结果与前人的研究结果较为接近。

      表 1  莫霍面深度对比/km

      Table 1.  Comparison of the Moho Depths/km

      数据来源 中国浙闽隆起区 中国东海陆架盆地 中国钓鱼岛隆褶带 冲绳海槽盆地 琉球隆褶区 琉球海沟 菲律宾海盆
      文献[7] -31~-29 -28.5~-27 -27.5~-25 -25~-18
      文献[14] -34~-29 -29~-27 -27~-22 -22~-16 -25~-19 -18~-12 -13~-12
      文献[5] -28~-25 -28~-26 -26~-16 -32~-18
      文献[8] -33~-28 -24~-19 -30~-26 -26~-14 -18~-4
      本文 -34.4~-29 -28.5~-26 -26~-24 -24~-21 -26.4~-23 -22~-15 -14~-8.8
      注:文献[8]中,中国浙闽隆起区和中国东海陆架盆地作为一个整体区域进行研究

      莫霍面深度三维图(见图 8)清晰地展现了中国东海海域莫霍面深度呈西北低、东南高的特征,由北西向南东依次呈现凹-凸-凹-凸的特征,地壳类型由陆壳向洋壳转变。中国浙闽隆起区莫霍面较深,该地区地壳很厚,是典型的大陆型地壳。中国东海陆架盆地的莫霍面变化非常平缓,呈现中间低、东北方向和西南方向偏高的特点。中国钓鱼岛隆褶带莫霍面逐渐抬升,是大陆莫霍面的延伸,变化平缓说明该地区地壳结构变化不大,基本处于均衡状态,是欧亚大陆前缘的板内盆地[15]。冲绳海槽盆地莫霍面抬升,在盆地中间形成一个明显的凸起,此处地壳发生了减薄,已经呈现出陆壳-过渡型地壳的特征。琉球海沟和菲律宾海盆的莫霍面变化剧烈,海沟西北侧莫霍面由低向高抬升,地壳具有过渡型地壳的特征。根据文献[16]分析的中国东海新构造运动,该地区菲律宾海板块沿琉球海沟往中国东海陆架边缘俯冲所引起的张裂、拉伸作用而形成的琉球岛弧隆褶带和弧后冲绳海槽盆地,具有典型的沟-弧-盆构造系统的特征。沿东南侧往菲律宾海盆方向,当莫霍面抬升至一定高度,莫霍面明显凸起,地壳明显减薄,该地区地壳具有大洋型地壳特征。

    • 本文以高精度的重力数据为基础,在剖面反演得到的地质-地球物理模型和地震资料约束下,采用Parker-Oldenburg界面反演方法得到了中国东海海域莫霍面深度,可为进一步认识该地区地壳深部结构和演化过程提供依据。

      中国东海海域莫霍面呈现两凹两凸的特征,与构造形态符合良好。整个区域莫霍面深度为-34.4~-8.8 km,且具有东西分带、南北分块的特点。浙闽隆起区莫霍面深度为-34.4~-29 km,走向呈北东、北北东向,属于典型的大陆型地壳。往东向中国东海陆架盆地方向的莫霍面抬升。中国东海陆架盆地莫霍面深度变化平缓,深度范围为-28.5~-26 km,且与中国浙闽隆起区莫霍面有良好的相关性。中国钓鱼岛隆褶带莫霍面逐渐抬升,莫霍面深度为-26~-24 km,至冲绳海槽盆地的莫霍面深度为-24~-21 km,走向呈北东向,此处地壳呈现出陆壳-过渡型地壳的特征。琉球隆褶区莫霍面由东侧的北东向,往西呈近东北向,莫霍面深度为-26.4~-23 km。琉球海沟莫霍面迅速抬升,东侧呈近北东向,往西呈近东北向,深度为-22~-15 km,至菲律宾海盆莫霍面抬升至-8.8 km,其地壳已完全属于大洋型地壳。

参考文献 (16)

目录

    /

    返回文章
    返回