1989年以来华北地区中强震前地震矩加速释放空间扫描特征的研究*

李　霞，蒋长胜，闫德桥，张志慧，尹相好，陈传华，王　铮

(1.山东省地震局，山东 济南 250014；2.中国地震局地球物理研究所，北京 100081)

1989年以来华北地区中强震前地震矩加速释放空间扫描特征的研究*

李霞1，蒋长胜2，闫德桥1，张志慧1，尹相好1，陈传华1，王铮1

(1.山东省地震局，山东 济南 250014；2.中国地震局地球物理研究所，北京 100081)

摘要：利用地震矩释放程度m值对1989年以来华北地区12次MS≥5.0中强地震前的矩加速释放现象进行了研究。基于删除余震后的地震目录，提取了与矩加速释放相关的特征时—空尺度作为AMR现象在华北地区的扫描模型参数。采用震前1年逐3个月的扫描方式，得到了12次震例在震中附近与时间相依的m值的空间演变图像，并考察了震前震源区附近的Benioff应变释放情况。结果表明，华北地区MS5.0～5.5地震前AMR现象与实际震中和发震时刻的关联性不明显；而MS5.6～6.4地震前67%的震例在震中附近出现稳定和集中分布的矩释放加速区域，说明华北地区MS≥5.5中强震前AMR现象具有一定程度的普遍性，利用AMR现象对未来中强震发震地点的预测有一定的潜能和可挖掘性，但仍需其他方法相互佐证。

关键词：地震矩加速释放；Benioff应变释放；空间扫描；华北地区

0引言

根据临界点理论，地震的孕育过程是一种临界现象(Vere-Jones，1977)，这一临界现象的典型特征是大地震之前往往存在地震活动的相关长度增长及主震周围矩释放加速(AMR)现象。Bowman等(1998)研究发现，几次大震前震源区地震能量会出现明显幂率加速释放过程，符合临界理论的标度率，进而使用震前累积应变释放曲线的幂率拟合和线性拟合的均方根残差比最小来定量描述AMR的“临界”时、空尺度。随后人们对世界上其他地区的震例进行了类似的研究，Sykes和Jaume(1990)研究表明美国旧金山几次强震前都有5.0级左右中等地震活动的矩释放加速过程；Ellsworth 等(1981)指出，强震前小震发生率并没有明显变化，而中等地震的累积应变能呈指数增长。我国地震学家在长期预报研究实践中也将AMR现象作为强震预测的主要前兆之一(刘蒲雄，1982)。杨文政和马丽(1999)用类似方法对中国震例开展研究，发现三分之二以上的中强震在震前较长时间均出现了地震矩释放加速的现象；Jiang和Wu(2005)的研究指出，中国大陆有60%的MS≥6.0地震前存在稳定的地震矩释放加速特征，但在多时、空尺度的二维空间，一些地震前的矩释放表现出加、减速交替出现的复杂现象。更有研究表明，主震前临界区范围内地震活动呈现加速释放，而在较小的孕震区范围则出现地震活动的减弱态势(苗青壮，蒋海昆，2008；Papazachosetal.，2000，2006)。

华北地区是一个内部结构复杂和断裂活动强烈的地区。按其主要地震构造带和区域地震构造格局，可划分为两个一级活动地块：鄂尔多斯活动块体和太行冀鲁西活动块体，其分界带为山西地堑系，其中太行冀鲁西地块又可分为太行山块体、华北平原块体和胶辽块体，其分界带为太行山山前断裂带和郯庐断裂带，还有华北北部的阴山—燕山块体，其分界带为燕山—渤海断裂带(Xu，Gao，2000；Huang，Guo，1998)。华北地区的强烈地震几乎全部发生在这些活动块体的边界带上，其中山西地震带发生MS>7地震8次，华北平原带3次，郯庐地震带3次以及燕山—渤海地震带8次，给我们的国计民生带来了巨大的损失。因此，华北地区一直是我国地震监测防预的重点监视区和研究地震前兆机理、检验地震预测理论和探索地震数值预测方法的一个重要的“实验场地”。

本文在以往研究的基础上，以1989年以来华北地区12次MS≥5.0中强地震为例，利用地震矩释放程度的时—空扫描方法，系统验证华北地区中强震前地震能量的释放特征，进而讨论AMR现象在华北地区是否具有普遍性的前兆意义。

1地震矩加速释放现象和释放程度的时—空扫描

“地震破裂时间分析”法是利用给定区域的前兆性地震活动来预测未来地震的发震时间和强度的一种中长期预测方法，其函数定义一般可表示为(Bufe，Varnes，1993)

S(t)=A+B(tf-t)m.

(1)

式中，S(t)是[0,t]时段内地震活动的累积应变能，tf为预期地震的破裂时间，t为前兆性地震事件的发震时间，t1表征 Benioff应变为减速释放状态，相应的m越大减速释放特征越明显；m=1表示Benioff应变匀速释放。

对临界区域的估算依据是使区域内的样本数据能被幂率关系拟合得最好，“最优搜索办法”(Bowanetal.，1998)是以强震震中为圆心，用不同半径范围内的观测数据连续计算曲率参数C值(震前Benioff应变释放的幂律拟和均方差与线性拟和均方差之比)，参数C为最小且幂指数m<0.8的圆型区域被认为是主震的临界区域。基于上述理论，杨文政和马丽(1999)利用空间临界圆重叠率最高的方法来寻找未来强震的临界区域；也有研究认为这种“最优搜索”仅仅是对区域数据本身的拟合，显著存在的问题是相同区域构造和相近震级的中强震之间的临界时、空尺度差异较大*Kossobokovv.2007.Earthquake forecast/prediction：problem of verification，accuracy and limitations.ICTP.。因此蒋长胜和吴忠良(2009)仅利用m值作为描述震前矩加速释放“程度”的参量，在时—空—地震序列的完整性震级组成的三维空间内，使用固定时间尺度和多空间尺度系统研究了2008年3月21日新疆于田MS7.3强震前AMR现象存在的客观性。研究表明矩加速释放特征在时间上有逐渐逼近临界状态的趋势，空间上与实际震中位置也有较好的对应的关系。这种利用m值的时—空演化来研究地震孕育的物理过程和寻找未来中强震的临界区域的方法称为地震矩释放程度的时—空扫描。这种方法实质上是假定实际震中和发震时刻未知，使用固定的时—空尺度，在时—空坐标上扫描震前的矩释放程度，考察研究区中强震前AMR现象与实际震中和发震时刻之间的相关性。

2震例选取和所用资料

笔者以1989年以来华北地区12次MS≥5.0地震作为研究对象，选取(34°～42°N，108°～123°E)的空间范围为研究区域。研究区中强地震的活动表现出明显的分区、交替、成组发生的特征，在同一地震活跃期中发生的中强地震，往往时间相连，地点相近，构成一个典型的地震成组集中活动区，而与此同时其他地区几乎没有6级以上地震发生。1900年以来，华北地区中强地震的活动经历了5次周期性轮回过程(孙加林，章瑞，2007)，目前正处于第五活跃期的末期或第五个平静期(图1)。第五活跃期以1989年山西大同和阳高间MS5.9地震为起始，至今已发生了12次MS≥5.0的中强地震。如图2所示，这12次地震主要分布在燕山—渤海地震带和山西地震带，而在活动地块内部及郯庐地震带、华北平原地震带和黄海地震带仅有个别地震发生，在空间上分区成组分布的特征明显。

笔者使用中国地震台网中心提供的1970年以来《中国地震月报目录》，该目录由各区域台网产出的资料汇编而成。地震目录的完整性震级MC是进行区域地震活动性和地震危险性分析的重要基础，使用“最大曲率”法(MAXC)、分别为90%和95%的拟合优度检验法(GFT)计算研究区MC随时间的演化过程。通过设定优先级为GFT95%>GFT90%>MAXC，选取优先级高的MC结果作为最终结果。由图3可见，除1993年和2000年研究区最优MC分别为ML2.2和ML2.4，其它时段均小于或等于ML2.0，因此使用ML2.5作为计算Benioff应变释放的最小完整性震级。Benioff应变被定义为地震波能量的平方根，据经验关系可由震级直接得到应变能(Gutenberg，Richter，1956；Kanamori，1977)。再据式(1)即可计算出某一时段内地震矩释放参数A、B、m及未来中强地震的发震时刻。为更客观地考察华北地区中强震前的AMR现象，我们利用矩释放程度的时—空扫描方法对上述12次MS≥5.0地震进行回顾性检验，表1给出了12次震例的震源参数。

3矩加速释放时—空扫描模型参数的选取

本文在矩释放时—空扫描模型参数的选取中，结合前人给出的有关AMR现象中震级与“临界”时、空尺度的关系(Ben-Zion，Lyakhovsky，2002；彭克银，尹祥础，2003)分别计算了12次震例前的前兆性地震序列在时—空(t-r)二维空间内m值的分布特征，然后用归一化的方法给出华北地区中强震前具有普遍意义的AMR现象的“临界”时—空尺度。近年来在许多地震预测模型的研究中对所使用地震目录是否删除余震都采取了折衷的做法，同时给出删除余震前、后的两种预测结果。在以往对AMR现象的研究中，蒋长胜(2004)研究指出余震删除与否，对m值计算结果的影响并不明显；李宇彤和蒋长胜(2010)对东北地区中强震的研究表明，使用删除余震目录得到的时—空尺度与前人研究结果更接近。因此本文使用G-K法(Gardner，Knopoff，1974)删除余震后的地震目录，将(1)式中的tf固定为研究震例的实际发震时刻；t为震前2～14 a，以1 a为滑动窗长；r则是以实际震中为圆心，10 km为滑动步长，在30～250 km的范围内多尺度半径选取地震样本进行Benioff应变释放曲线的幂律拟和，由此得到m值及m值的拟合标准差Δm。在对累积Benioff应变释放曲线进行非线性拟合时，给定参数的初值，进行迭代计算函数极小点参数值。拟合参数的结果通常与初值的选取有关，因此选取多初值拟合，用拟合的最小残差值确定“拟合最佳解”的办法来提高拟合的精确度(蒋长胜，2004)。在计算结果中为确保m值的可靠性，删除了样本数少于5个、非线性拟合与线性拟合均方差残差比C>1、拟合标准差Δm>0.5的计算结果。因本文仅考虑地震矩加速释放模型(即m<1的情况)，图4给出研究时段内华北地区12次中强震前m值在t-r二维空间内克里金插值法绘出的等值线分布图，表1中t-r列给出各次震例前m最低值出现的时—空二维尺度。

表1　1989年以来华北地区MS≥5.0地震震源参数

注：-表示Benioff累积应变减速释放；+表示Benioff累积应变加速释放；0表示Benioff累积应变匀速释放.

由图4可见，12次中强震前均存在或强或弱的AMR现象，但m值空间分布较为复杂，低值出现的时—空尺度差异显著。如1998年河北张北MS6.2、1999年MS5.6和2006年河北文安MS5.1地震前m低值现象明显，均出现m<0.3的时—空分布；1992年黄海MS5.3和1991年山西大同MS5.8地震前出现大范围矩释放的加速与减速交替发生的复杂图像；1995年山东苍山MS5.2地震与1995年河北唐山MS5.0地震前累积应变能几乎为匀速释放状态；1991年山西忻州MS5.1、1989年山西大同MS5.7和1999年山西浑源MS5.6地震前，距离震中较远的区域出现应变能释放加速，而离震中较近的孕震区则表现出地震矩的匀速释放特征。为提取华北地区具有普遍应用价值的AMR时—空特征尺度，本文将12次中强震按震级大小分为两类，即MS5.0～5.5和MS5.6～6.4，同时定义了地震矩加速特征q值为

(2)

式中，n为震例个数，t和r是计算Benioff应变曲线的时、空参数，mi(t,r)为i第个震例的幂指数m值在二维空间的分布；RMSi-linear和RMSi-nonlinear分别为第i个震例在t-r二维空间中计算样本的线性拟合和幂律拟合的均方根残差，其中RMSi-linear/RMSi-nonlinear即为Bowman等(1998)定义的曲率参数，用于判断临界圆“最优搜索半径”的C值。这里的矩加速特征q值实际上是每个震例的特征值qi归一化后，多震例所反映的地震矩加速释放的平均时、空特征。图5为提取出的华北地区中强震前AMR“优势”时—空特征尺度，即MS5.0～5.5地震优势时间窗t=4 a、扫描半径r=110 km；MS5.6～6.4地震优势时间窗t=6 a、扫描半径r=180 km，此参数与前人研究结果较为接近，基本符合最佳搜索半径r与震级的线性相关。

4华北地区中强震前的矩加速释放时空扫描

按照地震的临界性概念，一次大地震的孕育初期，孕震系统并非一直处于临界状态，而是随着周围应力状态的不断演化逐渐接近并进入临界状态。在这一过程中，通常情况下孕震区中小地震活动的频度和强度不断增加，应变能释放加速，直到强震发生，临界状态被破坏(Ben-Zio，Lyakhovsky，2002)。如果主震的震中在临界点附近，主震发生前就可以通过对周围地震的频度或强度在一定时—空范围的扫描寻找未来震中的位置。本文利用上述提取出的华北地区矩加速释放时—空扫描模型参数对研究时段内12次MS≥5.0中强震进行时空扫描。首先将研究区划分为0.15°×0.15°的网格空间，网格节点作为未来可能的震中位置，地震可能的“破裂时间”tf假定为震前1年、9个月、6个月、3个月、1个月和实际发震时刻，分别计算震前矩释放程度m值的空间分布。扫描结果显示，6次MS5.0～5.5地震前m值空间分布均杂乱无序，没有明显的应力优势集中区，并且与实际震中和发震时刻之间没有关联性可言。图6为2006年7月4日河北文安MS5.1地震前AMR现象的扫描图像，震中附近没有出现明显的矩释放加速区域，而是随着实际发震时刻的逐渐逼近，AMR现象逐渐远离实际震中，外围远场的AMR演化图像也没有任何与实际震例时间相依的规律可言。对于6次MS5.6～6.4震例的AMR扫描结果则理想很多，6次中强震前都出现较长时间的AMR现象，且空间分布较为集中与实际震中有很好的对应关系。以1998年1月10日河北张北MS6.2地震为例(图7)，震前1年至6个月，在实际震中附近出现明显、稳定的矩加速释放区，随着周围应力场的调整，震前3个月至发震，震中附近的异常区域逐渐收缩，m<0.5的低值区域从不同方位向震中区逐渐逼近，时间上也表现出逐渐趋近于临界状态的趋势，这种时—空上的一致性说明地震的孕育过程与AMR现象在物理上有一定的相关性。与此同时在沂沭带南段也逐渐形成一个M<1的矩加速释放区，但在之后的时间里，沂沭带南段及周围地区并未发生MS≥5地震，也说明了AMR现象及“破裂时间分析”等方法在实际地震预测中的不确定性和复杂性。为进一步考察华北地区中强震前震中区附近的矩加速释放情况，分别以实际震中为圆心，采用与时—空扫描相同的模型参数(即MS5.0～5.5的地震半径r=110 km，时间窗t=4 a；MS5.6～6.4的地震半径r=180 km，时间窗t=6 a)，计算了12次MS≥5.0 地震震中区附近的Benioff累积应变能释放情况，相应的Benioff应变释放线性和非线性拟合情况(图8)。从6次MS5.0～5.5地震的Benioff应变释放曲线来看，除1992年黄海MS5.3地震外均表现为近匀速释放状态，与对m值的时—空扫描图像结果基本一致，即震前在孕震区没有明显的应变能加速释放或地震频次增加的现象。而6次MS5.6～6.4中强震虽然在时—空扫描中存在大面积、稳定的AMR现象，但对于破裂成核点附近的Benioff应变能释曲线仅4例地震存在矩加速释放现象，这说明AMR现象仅是中强震孕育发生的模式之一。从表1中可以清晰的看出，12次震例中仅5例震前在破裂成核点附近存在加速释放特征，3例震前表现出应变的减速释放特征，另外4例为匀速释放。值得一提的是，在MS5.0～5.5震例中仅1例在震前表现出应变能加速，占总数的17%；而MS5.6～6.4的震例中有4例在震前存在矩加速释放特征，占总数的67%。因此认为华北地区MS≥5.5中强震前的AMR现象更具有一定的普遍性，但却不是充要条件。

中小地震丛集现象往往被认为是强震发生的前兆或构造活动的外释信号，成为了地震预测的重要依据之一。2012年以来，华北地区最值得关注的震情活动是环渤海地区的中等地震及震群事件的频发，先后发生了辽宁灯塔M5.1、内蒙通辽M5.3、吉林前郭M5.8和河北滦县M4.8、山东莱州M4.6、乳山M4.3等中等地震，以及山东长岛、莱州、乳山、辽宁海城、盖州和河北赞皇、蔚县等一系列围绕渤海及郯庐断裂带北段的典型事件，这足以说明环渤海地区的应力背景处于高值水平。本文应用提取出的华北地区矩加速释放模型参数，采用删除余震并剔除干扰事件(M≥5.0)后MC2.5的地震目录，假定发震时刻为2014年1月26日，来检验目前环渤海地区的AMR现象及Benioff应变释放曲线情况，结果如图9所示。

图9a中显示，在环渤海及黄海地区集中出现了北北西向展布的大面积的低值区域，矩加速释放特征显著。Benioff应变释放曲线表明自2012年开始该区域处于应变能加速累积的阶段，曲率参数C值为0.62，低于华北地区1989年以来12次中强震前的任何一次(图9b)。

5结论和讨论

本文利用地震矩释放程度m值对1989年以来华北地区12次中强震前的AMR现象进行了回顾性震例检验。(1)通过对值在t-r二维空间的图像分析，提取出了研究区与AMR现象有关的特征时—空尺度，即MS5.0～5.5地震在半径r=110 km，时间窗长t=4 a，MS5.6～6.4 地震在半径r=180 km，时间窗长t=6 a的情况下，AMR现象更具有“优势”和普遍性。(2)利用得到的特征时—空尺度对每个震例进行震前1 a至震时，依次递减3个月在空间范围的扫描。结果显示，除1992年黄海MS5.3地震外的5次MS5.0～5.5地震前AMR现象没有形成明显的优势集中区，与实际震中和发震时刻没有关联性；6次MS5.6～6.4地震前在震中区附近均出现显著、稳定的距加速释放区域，值随着实际发震时刻的临近而逐渐减小或向震中位置逐渐收缩，AMR现象出现的空间区域与实际震中也有较好的对应关系。(3)考察了震前震中区附近的Benioff应变释放情况，12次震例中有5例在破裂成核点附近存在AMR现象，与空间扫描结果基本一致，说明华北地区MS≥5.5中强震前AMR现象存在一定的普遍性，与区域地震孕育的物理过程有一定的相关性。(4)对研究区目前地震矩释放程度时—空扫描的结果显示，环渤海地区存在显著的矩加速释放特征，Benioff应变释放曲线表明背景应力水平较高。(5)通过上述研究，不难发现地震能量的累积释放与时间的关系比简单的用AMR模型来描述要复杂很多，一些地震前没有AMR现象出现，出现AMR现象也不一定会有中强地震发生，因此用矩释放程度的时—空扫描方法对未来中强震发震地点的预测需要与其它地震学预测方法相互映证，综合考虑。

本研究在成文过程中得到山东省地震局刘希强研究员、山长仑高级工程师的帮助，审稿专家提出了诸多有益建议，在此深表感谢。

参考文献：

邓启东，张培震，冉勇康，等.2003.中国活动构造与地震活动[J].地学前缘，10(S1)：66-73.

蒋长胜，吴忠良.2009.强震前的“加速矩释放”(AMR)现象：对一个有争议的地震前兆的回溯性震例研究[J].地球物理学报，52(3)：692-702.

蒋长胜.2004.中国一些中强震前地震活动的矩加速释放现象[D].北京：中国地震局地球物理研究所.

蒋海昆，吴琼，董翔，等.2009.不同温压条件下声发射应变能释放特征——加速模型参数物理含义的初步讨论[J].地球物理学报，52(8)：2064-2073.

李宇彤，蒋长胜.2010.东北第五活跃期中强震前的矩释放空间扫描特征研究[J].地震学报，32(6)：680-6956.

刘蒲雄.1982.强震前地震活动增强特征[J].地震，2(2)：5-7.

苗青壮，蒋海昆.2008.地震临界现象及相关尺度增长定量检测方法研究综述[J].地震，28(2)：79-86.

彭克银，尹祥础.2003.用临界点理论讨论应变能加速释放现象和孕震区尺度[J].中国地震，19(4)：425-430.

孙加林，章瑞.2007.中国大陆20世纪以来第五活跃期的复杂性及其本质问题[J].地震，27(2)：31-40.

杨文政，马丽.1999.地震活动加速模型及其在中国的应用[J].地震学报，21(1)：32-41.

Ben-Zion Y.，Lyakhovsky V..2002.Accelerated seismic release and related aspects of seismicity patterns on earthquake faults[J].Pure and Applied Geophysics，159(10)：2385-2412.

Bowman D.D.，Ouillon G.，Sammis C.G.，etal..1998.An observational test of the critical earthquake concept[J].J.Geophys.Res.，103(B10)：24359-24372.

Bufe，V.，Varnes D.J.1993.Predictive modeling of seismic cycling of the greater San Francisco Bay region[J].J.Geophys.Res.，98(B6)：9871-9883.

Ellsworth W.L.，Lindh A.G.，Prescott W.H.，etal..1981.The 1906 San Francisco earthquake and the seismic cycle[J].In：Simpson，D.W.and Richards，P.G.(eds.)，Earthquake Prediction：An International Review，21-27.

Gardner J.J.，Knopoff L..1974.Is the sequence of earthquakes in Southern California，with aftershocks removed，Poissonian[J].Bull.Seism.Soc.Am.，5(64)：1363-1367.

Gutenberg B.，Richter C.F..1956.Magnitude and energy of earthquake[J].Annali di Geofisica，9(1)：1-15.

Huang L.R.，Guo L.Q..1998.Crustal horizontal motion feature and its mode in North China[J].Crustal Deformation and Earthquake.18(2)：20-27.

Jiang C.S.，Wu Z.L..2005.Accelerating strain release before strong earthquakes：more complex in the real world[J].Journal of the Graduate School of the Chinese Academy of Sciences，22(3)：286-291.

Kanamori H..1977.The energy release in Great earthquakes[J].J.Geophys.Res.，82(B201)：2981-2987.

Papazachos B.，Papazachos C..2000.Accelerated pre-shock deformation of broad regions in the Aegean Area[J].Pure and Applied Geophysics，157(10)：1663-1681.

Papazachos C.B.，Karakaisis G.F.，Scordilis E.M.，etal..2006.New observational information on the precursory accelerating and decelerating strain energy release[J].Tectonophysics，423(1-4)：83-96.

Sykes L.R.，Jaume S.C..1990.Seismic activity on neighboring fault s as a long-term precursor to large earth-quakes in the San Francisco Bay area[J].Nature，348(6302)：595-599.

Vere-Jones D..1977.Statistical theories of crack propagation[J].Math.Geol，9：455-481.

Xu J.，Gao Z.W..2000.The structura1 characters of the piedmont fault zone of Taihang mountain[J].Seismology and Geology，22(2)：111-122.

Firstly，we study the accelerating moment release(AMR)before the 12MS≥5.0 moderate-strong earthquakes in North China region since 1989 by using the accelerating moment releasem-value.Secondly，based on the catalogue of deleted aftershock，we extracted a temporal-spatial window related with AMR as the scanning model parameter in North China.Using the three monthly spatial scanning modes in one year before the earthquake，we obtained the temporal evolution variation ofm-value with time，and investigated the Benioff train release around the epicenter area of 12 earthquakes.The result indicates that there was no obvious correlation between the AMR appeared beforeMS5.6～6.4 earthquakes and origin time.There appeared stable and concentrated AMR area before the 67% earthquakes out of theMS5.6～6.4 studied events around the epicenter area，which showed that the AMR phenomenon appeared beforeMS≥5.0 moderate-strong earthquakes had a certain universality in North China.Using AMR to forecast the site of future moderate-strong earthquakes has some potential and develop ability，but it still needed other ways to evident.

Key words：seismic accelerating moment release；Benioff strain release；spatial scanning；the north China region

*收稿日期：2014-02-27. 基金项目：2012年度山东省自然科学基金青年项目(ZR2012DQ009)和“中国大陆未来10年M≥7.0地震危险区预测”工作专项联合资助.

中图分类号：P315.5

文献标识码：A

文章编号：1000-0666(2015)03-0359-11

Research on Spatial Scanning Characteristic of Accelerating Moment Release before Moderate-strong Earthquakes in North China Region since 1989

LI Xia1，JIANG Chang-sheng2，YAN De-qiao1，ZHANG Zhi-hui1，YIN Xiang-hao1，CHEN Chuan-hua1，WANG Zheng1

(1.Earthquake Administration of Shandong Province，Jinan 250014，Shandong，China) (2.Institute of Geophysics，CEA，Beijing 100081，China)

Abstract