地震 ›› 2020, Vol. 40 ›› Issue (3): 52-64.doi: 10.12196/j.issn.1000-3274.2020.03.005
肖卓1,2, 高原2
收稿日期:
2019-03-25
出版日期:
2020-07-31
发布日期:
2020-07-28
通讯作者:
高原,研究员。E-mail:qzgyseis@163.com
作者简介:
肖卓(1993-),男,湖南娄底人,在读博士研究生,主要从事地震成像研究。
基金资助:
XIAO Zhuo1,2, GAO Yuan2
Received:
2019-03-25
Online:
2020-07-31
Published:
2020-07-28
摘要: 基于高精度的三维波场模拟和伴随方法, 伴随成像成功实现了天然地震全波形成像在区域和全球尺度下的应用。伴随成像技术基于谱元法, 使用全三维、 多参数的初始模型对地震波场进行数值模拟, 通过正演波场和伴随波场的相互作用快速求取目标函数的梯度, 结合高性能计算技术, 实现大尺度全波场成像。相比于传统的地震层析成像方法, 伴随成像对地球内部结构异常体的描述更加精细和全面, 更适用于构造活跃地区的深部动力学研究。本文首先介绍伴随成像方法的基本原理, 随后阐述其具体实现流程, 并回顾该方法在地球深部结构研究中的应用, 最后对其未来发展做出展望。
中图分类号:
肖卓, 高原. 伴随成像及其在地球深部结构研究中的应用[J]. 地震, 2020, 40(3): 52-64.
XIAO Zhuo, GAO Yuan. A Review of Adjoint Tomography and Its Application in Earth’s Interior Imaging[J]. EARTHQUAKE, 2020, 40(3): 52-64.
[1] Aki K, Lee WHK. Determination of three-dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes: 1. A homogeneous initial model[J]. Journal of Geophysical Research, 1976, 81(23): 4381-4399. [2] Replumaz A, Capitanio F A, Guillot S, et al. The coupling of Indian subduction and Asian continental tectonics[J]. Gondwana Research, 2014, 26(2): 608-626. [3] 雷建设, 赵大鹏, 苏金蓉, 等. 龙门山断裂带地壳精细结构与汶川地震发震机理[J]. 地球物理学报, 2009, 52(2): 339-345. LEI Jian-she, ZHAO Da-peng, SU Jin-rong, et al. Fine seismic structure under the Longmenshan fault zone and the mechanism of the large Wenchuan earthquake[J]. Chinese Journal of Geophysics, 2009, 52(2): 339-345 (in Chinese). [4] Pei S, Zhang H, Su J, et al. Ductile Gap between the Wenchuan and Lushan Earthquakes Revealed from the Two-dimensional Pg Seismic Tomography[J]. Scientific Reports, 2015, 4: 6489. [5] Zhao D. Global tomographic images of mantle plumes and subducting slabs: Insight into deep Earth dynamics[J]. Physics of the Earth and Planetary Interiors, 2004, 146(1-2): 3-34. [6] Huang J, Zhao D. High-resolution mantle tomography of China and surrounding regions[J]. Journal of Geophysical Research, 2006, 111(B9): B09305. [7] Chen M, Niu F, Tromp J, et al. Lithospheric foundering and underthrusting imaged beneath Tibet[J]. Nature Communications, 2017, 8: 15659. [8] Rawlinson N, Fichtner A, Sambridge M, et al. Seismic tomography and the assessment of uncertainty[J]. Advances in Geophysics, 2014, 55: 1-76. [9] Liu Q, Gu Y J. Seismic imaging: From classical to adjoint tomography[J]. Tectonophysics, 2012, 566-567: 31-66. [10] 张丽娜, 罗艳, 陈智勇, 等. 福建及其邻区背景噪声面波群速度层析成像[J]. 地震, 2018, 38(3): 134-143. ZHANG Li-na, LUO Yan, CHEN Zhi-yong, et al. Surface wave group velocity tomography imaging for Fujian Area from ambient noise[J]. Earthquake, 2018, 38(3): 134-143 (in Chinese). [11] 张娜, 赵翠萍, 周连庆. 三峡水库区上地壳三维精细速度结构成像[J]. 地震, 2018, 38(4): 37-48. ZHANG Na, ZHAO Cui-ping, Zhou Lian-Qing. Upper crust velocity structure tomography of the Three Gorges Reservoir Area based on dense observation[J]. Earthquake, 2018, 38(4): 37-48 (in Chinese). [12] 杨峰. 华北北部地区地壳三维P波速度结构的双差地震层析成像[J]. 地震, 2019, 39(1): 58-71. YANG Feng. Double-difference seismic tomography of 3D P-wave velocity crustal structure under the northern part of north China[J]. Earthquake, 2019, 39(1): 58-71 (in Chinese). [13] Virieux J, Operto S. An overview of full-waveform inversion in exploration geophysics[J]. Geophysics, 2009, 74(6): WCC1-WCC26. [14] Wang J, Yang D, Jing H, et al. Full waveform inversion based on the ensemble Kalman filter method using uniform sampling without replacement[J]. Science Bulletin, 2019, 64(5): 321-330. [15] Komatitsch D, Tromp J. Introduction to the spectral element method for three-dimensional seismic wave propagation[J]. Geophysical Journal International, 1999, 139(3): 806-822. [16] Komatitsch D, Tromp J. Spectral-element simulations of global seismic wave propagation-I. Validation[J]. Geophysical Journal International, 2002, 149(2): 390-412. [17] Komatitsch D, Tromp J. Spectral-element simulations of global seismic wave propagation-II. Three-dimensional models, oceans, rotation and self-gravitation[J]. Geophysical Journal International, 2002, 150(1): 303-318. [18] Tape C, Liu Q, Alessia M, et al. Adjoint tomography of the southern California Crust[J]. Science, 2009, 325, 988. [19] Zhu H, Bozda E, Peter D, et al. Structure of the European upper mantle revealed by adjoint tomography[J]. Nature Geoscience, 2012, 5(7): 493-498. [20] Bozda E, Peter D, Lefebvre M, et al. Global adjoint tomography: First-generation model[J]. Geophysical Journal International, 2016, 207(3): 1739-1766. [21] Tarantola A. Linearized Inversion of seismic reflection data[J]. Geophysical Prospecting, 1984, 32(6): 998-1015. [22] Tromp J, Tape C, Liu Q. Seismic tomography, adjoint methods, time reversal and banana-doughnut kernels[J]. Geophysical Journal International, 2004, 160 (1): 195-216. [23] Fichtner A, Bunge H-P, Igel H. The adjoint method in seismology I. Theory[J]. Physics of the Earth and Planetary Interiors, 2006, 157(1-2): 86-104. [24] 谢春, 刘玉柱, 董良国, 等. 伴随状态法初至波走时层析[J]. 石油地球物理勘探, 2014, 49(5): 877-883. XIE Chun, LIU Yu-zhu, DONG Liang-guo, et al. First arrival traveltime tomography based on the adjoint state method. Oil Geophysical Prospecting, 2014, 49(5): 877-883 (in Chinese). [25] Aster RC, Borchers B, Thurber CH. Parameter Estimation and Inverse Problems[M]. Waltham: Elsevier Academic Press, 2005. [26] Metropolis N, Ulam S. The Monte Carlo Method[J]. Journal of the American Statistical Association, 1949, 44(247): 335-341. [27] Meier U, Curtis A, Trampert J. Global crustal thickness from neural network inversion of surface wave data[J]. Geophysical Journal International, 2007, 169(2): 706-722. [28] Epanomeritakis I, Akçelik V, Ghattas O, et al. A Newton-CG method for large-scale three-dimensional elastic full-waveform seismic inversion[J]. Inverse Problems, 2008, 24(3): 034015. [29] Nolet G. A Breviary of Seismic Tomography: Imaging the Interior of the Earthand Sun[M]. Cambridge: Cambridge University Press, 2008. [30] Aki K, Richards PG. Quantitative Seismology, 2nd edition[M]. Mill Valley: University Science Book, 2002. [31] Tape C, Liu Q, Tromp J. Finite-frequency tomography using adjoint methods-Methodology and examples using membrane surface waves[J]. Geophysical Journal International, 2007, 168(3): 1105-1129. [32] Nocedal J, Wright S J. Numerical Optimization, 2nd edition[M]. Evanston: Springer, 2006. [33] Fichtner A, Kennett BLN, Igel H, et al. Full seismic waveform tomography for upper-mantle structure in the Australasian region using adjoint methods[J]. Geophysical Journal International, 2009, 179(3): 1703-1725. [34] Hjörleifsdóttir V, Ekström G. Effects of three-dimensional Earth structure on CMT earthquake parameters[J]. Physics of the Earth and Planetary Interiors, 2010, 179(3-4): 178-190. [35] French SW, Romanowicz BA. Whole-mantle radially anisotropic shear velocity structure from spectral-element waveform tomography[J]. Geophysical Journal International, 2014, 199(3): 1303-1327. [36] Dziewonski A M, Anderson D L. Preliminary reference Earth model[J]. Physics of the Earth and Planetary Interiors, 1981, 25(4): 297-356. [37] Kennett BLN, Engdahl ER, Buland R. Constraints on seismic velocities in the Earth from traveltimes. Geophysical Journal International, 1995, 122(1): 108-124. [38] Laske G, Masters G, Ma Z, et al. Update on CRUST1.0 - A 1-degree global model of Earth’s crust[C]. EGU General Assembly Conference. EGU General Assembly Conference Abstracts, 2013. [39] Kustowski B, Ekström G, Dziewoński AM. Anisotropic shear-wave velocity structure of the Earth’s mantle: A global model[J]. Journal of Geophysical Research, 2008, 113(B6): B06306. [40] Tromp J, Komatitsch D, Hjórleifsdóttir V, et al. Near real-time simulations of global CMT earthquakes: Near real-time simulations of CMT earthquakes[J]. Geophysical Journal International, 2010, 183(1): 381-389. [41] Zhu H, Bozda E, Tromp J. Seismic structure of the European upper mantle based on adjoint tomography[J]. Geophysical Journal International, 2015, 201(1): 18-52. [42] Maggi A, Tape C, Chen M, et al. An automated time-window selection algorithm for seismic tomography[J]. Geophysical Journal International, 2009, 178(1): 257-281. [43] Inoue H, Fukao Y, Tanabe K, et al. Whole mantle P-wave travel time tomography[J]. Physics of the Earth and Planetary Interiors, 1990, 59(4): 294-328. [44] Zhu H. Crustal wave speed structure of North Texas and Oklahoma based on ambient noise cross-correlation functions and adjoint tomography[J]. Geophysical Journal International, 2018, 214(1): 716-730. [45] Fichtner A, Trampert J. Resolution analysis in full waveform inversion: Resolution in full waveform inversion[J]. Geophysical Journal International, 2011, 187(3): 1604-1624. [46] Rawlinson N, Fichtner A, Sambridge M, et al. Seismic tomography and the assessment of uncertainty[J]. Advances in Geophysics, 2014, 55: 1-76. [47] Chen M, Niu F, Liu Q, et al. Multiparameter adjoint tomography of the crust and upper mantle beneath East Asia: 1. Model construction and comparisons[J]. Journal of Geophysical Research: Solid Earth, 2015, 120(3): 1762-1786. [48] Zhu H, Tromp J. Mapping tectonic deformation in the crust and upper mantle beneath Europe and the North Atlantic Ocean[J]. Science, 2013, 341(6148): 871-875. [49] Chen M, Niu F, Liu Q, et al. Mantle-driven uplift of Hangai Dome: New seismic constraints from adjoint tomography[J]. Geophysical Research Letters, 2015, 42(17): 6967-6974. [50] Panning MP, Lekic′ V, Romanowicz BA. Importance of crustal corrections in the development of a new global model of radial anisotropy[J]. Journal of Geophysical Research, 2010, 115(B12): B12325. [51] Tromp J, Luo Y, Hanasoge S, et al. Noise cross-correlation sensitivity kernels: Noise cross-correlation sensitivity kernels[J]. Geophysical Journal International, 2010, 183(2): 791-819. [52] Chen M, Huang H, Yao H, et al. Low wave speed zones in the crust beneath SE Tibet revealed by ambient noise adjoint tomography[J]. Geophysical Research Letters, 2014, 41(2): 334-340. [53] Liu Y, Niu F, Chen M, et al. 3-D crustal and uppermost mantle structure beneath NE China revealed by ambient noise adjoint tomography[J]. Earth and Planetary Science Letters, 2017, 461: 20-29. [54] Zhang C, Yao H, Liu Q, et al. Linear array ambient noise adjoint tomography reveals intense crust-mantle interactions in North China Craton[J]. Journal of Geophysical Research: Solid Earth, 2018, 123(1): 368-383. [55] Bensen GD, Ritzwoller MH, Barmin MP, et al. Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements[J]. Geophysical Journal International, 2007, 169(3): 1239-1260. [56] Tong P, Chen C, Komatitsch D, et al. High-resolution seismic array imaging based on an SEM-FK hybrid method[J]. Geophysical Journal International, 2014, 197(1): 369-395. |
[1] | 曾宪伟, 龙锋, 任家琪, 蔡新华, 李文君. 2019年6月17日长宁MS6.0地震前后b值时空变化分析[J]. 地震, 2020, 40(3): 1-14. |
[2] | 孙凯, 孟国杰, 洪顺英, 黄星, 董彦芳. 基于Sentinel-1A数据的四川长宁MS6.0地震同震形变场分析及断层滑动分布反演[J]. 地震, 2020, 40(3): 15-27. |
[3] | 赵策, 左可桢, 赵翠萍. 2019年6月17日四川长宁6.0级地震序列活动特征[J]. 地震, 2020, 40(3): 28-40. |
[4] | 王玮铭, 谭大诚. 2019年长宁MS6.0地震邻区分钟采样地电场分析[J]. 地震, 2020, 40(3): 41-51. |
[5] | 蒋雨函, 高小其, 王阳洋, 张磊. 中国新疆北天山和台湾南部陆地泥火山研究进展[J]. 地震, 2020, 40(3): 65-82. |
[6] | 马禾青, 杨明芝, 罗国富. 地震活动频次场研究[J]. 地震, 2020, 40(3): 99-111. |
[7] | 戴盈磊, 万永革, 梁永朵, 张文静, 惠杨. 基于震源机制解资料的辽宁地区现今构造应力场[J]. 地震, 2020, 40(3): 112-130. |
[8] | 黄颂, 姚丽, 姜楚峰. 基于幅相法的地磁低点位移异常可靠性分析方法[J]. 地震, 2020, 40(3): 131-141. |
[9] | 王晨晖, 刘立申, 任佳, 袁颖, 王利兵, 陈凯男. 主成分分析法和遗传算法优化的支持向量机模型在地震伤亡人数预测中的应用[J]. 地震, 2020, 40(3): 142-152. |
[10] | 王洪体, 崔仁胜, 王宏远, 叶鹏, 李丽娟. 对齐误差对测试地震计自噪声的影响[J]. 地震, 2020, 40(3): 153-166. |
[11] | 罗毅, 田云锋, 张苏, 张景发. 利用InSAR技术研究青藏高原冻土形变特征[J]. 地震, 2020, 40(3): 179-188. |
[12] | 薛艳, 姜祥华, 刘桂萍. 华北地区强震活动状态研究[J]. 地震, 2020, 40(2): 1-17. |
[13] | 贾漯昭, 宋美琴, 谢健健, 邢康, 高家乙. 华北块体地震活动显著增强预测效能评估[J]. 地震, 2020, 40(2): 38-48. |
[14] | 万森林, 张军龙, 刘明军, 贺为民, 李海龙, 郭长宝, 李智敏. 岷山断块的发震构造与地震活动性分析[J]. 地震, 2020, 40(2): 49-70. |
[15] | 李瑞莎, 张希, 贾鹏, 白卓立. 青藏块体东北缘跨断层异常场地比指标及其在地震预测中的应用[J]. 地震, 2020, 40(2): 71-81. |
阅读次数 | ||||||
全文 |
|
|||||
摘要 |
|
|||||