Mogi模型在长白山天池火山区的应用.pdfVIP

维普资讯 第 29卷 第 1期 地 震 地 质 VO1．29．NO．1 Mar．．2007 2007年 3月 SEISMOLOGY AND GEOLOGY Mogi模型在长 白山天池火 山区的应用 胡亚轩 王庆 良 崔笃信 王文萍 李 克 陈红卫 1)中国地震局第二监钡!I中心 ，西安 710054 2)吉林省地震局，长春 130022 摘 要 文中首先介绍了 Mogi模型的阻尼最小二乘反演 ，得出应用形变资料联合反演压力源参数 的公式，再指 出用等效源半径反演时不同初始值对反演结果的影响，说 明选择合适 的初始值可以方便 地得到可靠 的反演结果 。然后用此方法反演长 白山天池火 山区压力源 的大小和位置参数 。采用垂直 和水平位移的联合反演 ，数据为 2002--2005年共 4期的水准和 GPS观测资料 ，最后结合其它观测资 料，分析近年来岩浆可能的活动特征。结果显示岩浆的位置在发生变化，体积增量逐年减小，表明火 山岩浆在 2002--2005年的活动逐年减弱 。 关键词 阻尼最小二乘法 联合反演 长 白山天池火 山 岩浆活动特征 中图分类号：P317 文献标识码 -A 文章编号：0253—4967(2007)01—0144—08 O 引言 Mogi模型是模拟研究火 山地表形变的最简单适用的模型。该模型已用来模拟和解释大量 的火山区地面变形 (Yokoyama，1971；Deckereta1．，1983；Berrinoeta1．，1984；Yangeta1．， 1992)。应用火 山区的形变测量数据来反演岩浆压力源参数对火山的危险性评价有很大 的意 义。反演方法很多，经典的最小二乘方法在很多模型中有很好的应用 ，但是该模型有一定 的局 限性。如果所求函数在某一初始值附近作泰勒展开形成的偏导矩阵为A，则当 A是奇异矩 阵或可近似看作奇异矩阵时出现发散困难。阻尼最小二乘方法可用来解决此难题 ，该方法 已在 其它模型中得到了成功的应用(王庆 良等，1998)。 2002年6月汪清地震后，地震 、地化和形变等资料均表明长 白山天池火山区经历了一次比 较强烈 的岩浆上涌活动过程 ，而且火山变形明显地出现在天池 旧火 山 口附近 ，水平位移接近放 射状 ，与 Mogi模型的理论位移形状相符 (Mogi，1958；Yangeta1．，1992)。由于岩浆 的相对变化 较小 ，AA是奇异矩阵。本文在理论证明阻尼最tJ~--乘方法可行性的前提下，应用此方法反演 了长 白山天池火 山区压力源的位置和等效半径等参数 ，分析了近年来岩浆可能的活动特征。 1 Mogi模型 Mogi模型是将岩浆压力源置于均匀弹性半空间中，在源 的半径远小于源的深度时，可得出 由爆炸性或收缩性压力源引起的弹性半空间地表的水平变形和垂直变形与源的参数之间的表 (收稿 日期】 2006—07—05收稿，2006一l1—1O改回。 【基金项 目) 国家 自然科学基金 和地震科学联合基金(605004)共同资助。 维普资讯 1期 胡亚轩等 ：Mogi模型在长 白山天池火 山区的应用 l45 达式 。假定地面直角坐标系为 (，Y，z)，岩浆压力源中心坐标为 ( ，Y。，一d)(图 1)，如果用等 效半径 R表示岩浆流入或流出的体积 ，则根据 Mogi模型，岩浆房膨胀所引起的地表位移表达式 为 分量位移 ，u 为 Y分量位移 。r= 为径 向半 径 ，d Fig ． i i