基于三维大地电磁AR-QN反演的长白山天池火山区电性结构  

ELECTRIC STRUCTRUE MODEL OF TIANCHI VOLCANO IN CHANGBAI MOUNTAINS BASED ON THREE-DIMENSIONAL AR-QN MAGNETOTELLURIC INVERSION

作　　者：阮帅 汤吉[1] 董泽义[1] 王立凤[1] 邓琰 韩冰[1] RUAN Shuai;TANG Ji;DONG Ze-yi;WANG Li-feng;DENG Yan;HAN Bing(State Key Laboratory of Earthquake Dynamics,Institute of Geology,China Earthquake Administration,Beijing 100029,China;Key Laboratory of Continental Collision and Plateau Uplift Institute of Tibetan plateau Research,Chinese Academy of Sciences,Beijing 100101,China)

机构地区：[1]中国地震局地质研究所,地震动力学国家重点实验室,北京100029 [2]中国科学院青藏高原研究所,大陆碰撞与高原隆升重点实验室,北京100101

出　　处：《地震地质》2020年第6期1282-1300,共19页Seismology and Geology

基　　金：国家自然科学基金(41704078,41674081)资助。

摘　　要：长白山天池火山是具有潜在喷发危险的巨型活火山,对其深部低阻岩浆囊的精确探测和准确定位是人们关注的一个重要问题。应用大地电磁方法研究深部电阻率结构是确定岩浆囊的最有效途径之一。文中采用自主研发的三维大地电磁自适应正则化拟牛顿反演程序,对长白山天池火山区大地电磁数据进行了带地形三维反演,获取了较合理的深部三维电性结构模型。反演除拟合大地电磁主阻抗外,同时增加了对横向不均匀性更灵敏的倾子数据,并通过定性分析对比、模型验证的方法确认了三维反演的电性结构可靠性。三维电性结构模型比以前的二维解释结果更加准确合理,在有限稀疏的数据集下尽量精细地勾画了壳内低阻岩浆囊及通道的三维空间分布。三维电性结构显示壳内岩浆囊位于天池东北部深10~30km处,整体上岩浆囊和岩浆通道系统呈"V"字形分布,天池火山口的岩浆通道在天池以北10km的EW向剖面图上清晰可见。The Tianchi volcano in Changbai Mountains is a giant active volcano with potential risk of eruption,so more detailed researches of deep magma chamber as well as its closely related threedimensional electric structure model are being concerned by many scholars.Based on adaptive regularized three-dimensional quasi-Newton inversion program developed by the author,this paper carries out three-dimensional magnetotelluric inversion including topography using the existing data observed decades ago.Our three-dimensional magnetotelluric adaptive regularization quasi-Newton inversion method improves the conventional quasi-Newton method by approximating Hessian matrix of the data misfit using LBFGS formula instead of total Hessian matrix which is approximated in conventional quasi-Newton inversion.This not only guarantees the precise regularization-term Hessian matrix,but also allows the regularization parameter to be adaptively updated on every inversion iteration without destroying the descent of total objective function.Thus,the regularization parameter’s adaptive updating strategy on every iteration can be established based on L2-norm ratio of the data misfit function and regularization function,and our AR-QN inversion algorithm was implemented.The model synthetic data inversion results indicate that AR-QN inversion has very stable iteration flow,and is weakly dependent on initial model selection,which tends to be more suitable for the threedimensional inversion task when MT survey sites are sparsely distributed on the surface.Besides the impedances data,the tipper data which is more sensitive to horizontal conductive discontinuity was fitted in this three-dimensional inversion as well.Through comparison of qualitative analysis to measured data and anomalous-body-erasing forward modeling test,our new threedimensional electric structure was proved to be a reliable model under the constraints of current magnetotelluric data.The three-dimensional electric structure is more accurate and rational compared with the electric

关 键 词：长白山天池火山 电性结构 大地电磁 拟牛顿反演 自适应正则化 

分 类 号：P317.2[天文地球—固体地球物理学]