祝意青趙云峰李鐵明梁偉鋒徐云馬郭樹松

1）中國地震局第二監(jiān)測中心，西安 710054

2）中國科學院測量與地球物理研究所，大地測量與地球動力學國家重點實驗室，武漢 430077

3）中國地震局地質(zhì)研究所，北京 100029

2013年甘肅岷縣漳縣6.6級地震前后重力場動態(tài)變化

祝意青1，2）趙云峰1）李鐵明3）梁偉鋒1）徐云馬1）郭樹松1）

1）中國地震局第二監(jiān)測中心，西安 710054

2）中國科學院測量與地球物理研究所，大地測量與地球動力學國家重點實驗室，武漢 430077

3）中國地震局地質(zhì)研究所，北京 100029

利用青藏高原東北緣2011—2013年期間的流動重力觀測資料，系統(tǒng)分析了區(qū)域重力場動態(tài)變化及其與2013年7月22日岷縣漳縣6.6級地震發(fā)生的關(guān)系。結(jié)果表明：1）測區(qū)內(nèi)重力場異常變化與主要斷裂帶在空間上關(guān)系密切，反映沿主要斷裂帶（段）在2011—2013年期間發(fā)生了引起地表重力變化效應(yīng)的構(gòu)造活動或變形。2）岷縣漳縣6.6級地震前，測區(qū)內(nèi)先出現(xiàn)了較大空間范圍的區(qū)域性重力異常，后在震源區(qū)附近產(chǎn)生了局部重力異常及重力變化高梯度帶，其中，甘肅臨夏與岷縣一帶重力差異變化達150×10-8ms-2以上；這些可能反映岷縣漳縣地震前，區(qū)域及震源區(qū)附近均產(chǎn)生與該地震孕育、發(fā)生有關(guān)的構(gòu)造運動或應(yīng)力增強作用。3）重力場1a尺度動態(tài)變化圖像和差分動態(tài)演化圖像均反映岷縣漳縣6.6級地震孕育過程的最后2a出現(xiàn)較顯著的流動重力異常變化，地震發(fā)生在NE向重力變化高梯度帶上、重力變化零值線附近和等值線的拐彎部位。4）基于流動重力異常變化在岷縣漳縣6.6級地震前做過一定程度的中期預(yù)測，尤其是地點判定。

青藏高原東北緣 岷縣漳縣6.6級地震 重力觀測 重力變化 構(gòu)造活動

0 引言

2013年7月22日，在甘肅岷縣漳縣（34.5°N，104.2°E）發(fā)生了6.6級地震，震中位于甘肅東南部地區(qū)的臨潭-宕昌斷裂的中東段，該斷裂是這次地震的發(fā)震構(gòu)造（鄭文俊等，2013）。這次地震發(fā)生在中國地震重點監(jiān)視區(qū)的青藏高原東北緣地區(qū)，中國地震局系統(tǒng)自20世紀80年代以來就在該地區(qū)開展流動重力觀測工作（祝意青等，2012a），2008年汶川8.0級地震后觀測到甘東南地區(qū)存在重力異常變化，并對本次強震連續(xù)多年進行了年度震情跟蹤預(yù)測，尤其是強震地點的判定。中國地震局第二監(jiān)測中心在2011—2012年度地震趨勢研究報告中指出①，甘東南地區(qū)及其附近重力場在時間演化上變化強烈，空間上差異變化劇烈，存在強震危險背景，未來1a左右有發(fā)生6級以上強震的可能；2013年地震趨勢研究報告進一步指出①，重力變化沿區(qū)內(nèi)主要構(gòu)造斷裂帶出現(xiàn)重力變化高梯度帶，重力升、降差異達100×10-8ms-2以上，為顯著的強震背景異常。中國地震局年度趨勢預(yù)測匯總組編寫的《2013年度全國地震重點危險區(qū)匯總研究報告》，對這次地震也進行了時、空、強3要素的預(yù)測，2013年7月22日岷縣漳縣6.6級地震發(fā)生在重力異常判定的危險區(qū)內(nèi)。在此之前的相關(guān)研究也表明，基于流動重力觀測資料曾對2008年四川汶川8.0級和2013年四川蘆山7.0級地震做出了一定程度的中期預(yù)測（Zhu et al．，2010；祝意青等，2009，2013），預(yù)測危險區(qū)的中心距離中國地震臺網(wǎng)測定的汶川地震震中和蘆山地震震中均＜80km。這表明，依據(jù)區(qū)域重力場時-空變化的分析，可以開展強震中期預(yù)測的探索，尤其是強震可能發(fā)生地點的判定，這在地震3要素預(yù)測中尤為重要。因此，系統(tǒng)深入分析青藏高原東北緣重力觀測資料，進而對2013年甘肅岷縣漳縣MS6.6地震前后重力場的時空動態(tài)演化特征與規(guī)律、地震孕育發(fā)生過程中的前兆表現(xiàn)有一個較為客觀的認識，這無論是從地震預(yù)報角度，還是從地球動力學的角度來看，都是有實際意義的。

1 重力觀測與資料處理

青藏高原東北緣是中國大陸地殼運動最強烈、地震活動頻度最高、強度最大的地區(qū)之一。為了監(jiān)測青藏高原東北緣地區(qū)地殼運動變化、發(fā)現(xiàn)可能的地震中短期前兆，中國地震局第二監(jiān)測中心、甘肅省地震局和寧夏回族自治區(qū)地震局3個單位，自20世紀80年代分別在該地區(qū)建立了地震重力監(jiān)測網(wǎng)，通過不斷的優(yōu)化與整合，形成了目前的青藏高原東北緣整體重力監(jiān)測網(wǎng)（圖1），這有利于系統(tǒng)分析青藏高原東北緣重力場時空動態(tài)變化，探討地面重力變化與地殼運動、構(gòu)造活動和強震活動的關(guān)系（祝意青等，2012a，b）。2008年以前，每年對青藏高原東北緣進行1期重力觀測，汶川地震后，中國地震局加強了青藏高原東北緣的流動重力觀測，2009年起每年進行2期重力觀測。

圖1 青藏高原東北緣重力測量路線及活動斷裂略圖Fig.1 Map of the gravity surveying routes and main active faults of the northeastern margin of Qinghai-Tibet plateau.

對資料的處理：1）對青藏高原東北緣的流動重力觀測資料，以該地區(qū)準同期的絕對重力測量數(shù)據(jù)作為控制，實測標定重力儀的1次項系數(shù)，以消去儀器格值系數(shù)變化帶來的誤差（李輝等，2000；梁偉鋒等，2013）。2）平差計算時，先對多期重力觀測資料計算結(jié)果進行整體分析，初步了解各臺儀器的觀測精度后，合理確定各臺儀器的先驗方差，再重新平差計算，以得到最佳合理解算結(jié)果（祝意青等，2004，2005）。3）利用統(tǒng)一基準解算的重力平差計算結(jié)果較好，2011—2013年的點值平均精度在（8.6～11.8）×10-8ms-2之間。4）用最小二乘配置對重力觀測數(shù)據(jù)進行擬合推估，以便突出顯示構(gòu)造因素的重力效應(yīng)。

2 岷縣漳縣6.6級地震前后重力場動態(tài)變化

筆者曾對青藏高原東北緣重力場時變做過大量的研究（祝意青等，2004，2005），依據(jù)2005—2011年青藏高原東北緣的重力變化曾指出（祝意青等，2012a），汶川地震的發(fā)生使青藏高原深部物質(zhì)向NE運移，且部分通道受鄂爾多斯地塊的阻隔，青藏高原東北緣區(qū)域重力場出現(xiàn)大空間尺度的趨勢性顯著重力變化。汶川地震后，甘寧陜交界的六盤山斷裂帶與西秦嶺北緣斷裂帶及其之間地區(qū)重力場出現(xiàn)顯著的非均勻變化，該地區(qū)受汶川大地震的影響顯著。歷史上該地區(qū)及其附近曾發(fā)生過6級以上的強震，值得對該地區(qū)進行深入研究和加強監(jiān)視。2013年7月22日岷縣漳縣6.6級地震發(fā)生在西秦嶺北緣斷裂帶附近的臨潭-宕昌斷裂上，重力變化對這次強震有較好的反映。

本文主要分析2011年以來的重力變化，即岷縣漳縣6.6級地震前后的重力場動態(tài)變化（圖2，3）。

圖2 青藏高原東北緣1a尺度的重力場動態(tài)變化圖像（單位：10-8ms-2）Fig.2 Maps of gravity change patterns of the 1 year scale in the northeastern margin of Qinghai-Tibet plateau. a 2011-05—2012-05；b 2012-05—2013-05

2.1 1a尺度的重力場動態(tài)變化圖像

青藏高原東北緣每年進行2期重力測量，首先進行1a尺度的年變化分析，以了解其年尺度動態(tài)變化特征。圖2是研究區(qū)1a尺度的區(qū)域重力場的動態(tài)變化圖像，從圖2可以看出以下現(xiàn)象：

圖3 岷縣漳縣地震前后重力場差分動態(tài)變化圖像（單位：10-8ms-2）Fig.3 Maps of gravity change patterns between two adjacent observation campaigns before and after the Minxian-Zhangxian MS6.6 earthquake.a 2012-05—2012-10；b 2012-10—2013-05；c 2013-05—2013-10

（1）2011-05—2012-05期間，重力變化總體趨勢是自SW向NE出現(xiàn)由負向正逐漸增加的變化，重力變化等值線走向總體與祁連-海原大斷裂帶的走向基本一致，并沿斷裂構(gòu)造線出現(xiàn)重力變化高梯度帶，較好地反映了重力變化受區(qū)域應(yīng)力場作用和深大斷裂活動的控制。重力變化在自南向北由負向正的變化過程中，在武威地區(qū)及其附近出現(xiàn)局部重力異常區(qū)及等值線的拐彎，異常區(qū)中心位于莊浪河斷裂與毛毛山斷裂、金強河斷裂及冷龍嶺斷裂交會區(qū)附近，等值線拐彎的走向與NNW向的莊浪河斷裂走向基本一致。

（2）2012-05—2013-05期間，重力變化出現(xiàn)了新的態(tài)勢，以103°E為界重力出現(xiàn)分區(qū)變化。其一，測區(qū)東部的甘東南地區(qū)出現(xiàn)顯著重力差異運動，臨夏地區(qū)重力正值變化，岷縣、天水地區(qū)重力負值變化，重力差異變化達150×10-8ms-2，并沿岷縣—隆德形成重力變化高梯度帶，2013年7月22日甘肅岷縣漳縣6.6級地震發(fā)生在重力變化高梯度帶上。其二，測區(qū)西部的河西走廊地區(qū)在門源和祁連附近出現(xiàn)2個局部重力異常區(qū)，并沿祁連山構(gòu)造帶出現(xiàn)重力變化梯度帶，2013年9月20日甘肅肅南、青海門源交界（37.7°N，101.5°E）5.1級地震發(fā)生在門源和天祝2個局部異常區(qū)的過渡地帶，與門源異常有關(guān)的重力變化梯度帶上。

2.2 岷縣漳縣6.6級地震前后重力場差分動態(tài)變化圖像

為了分析岷縣漳縣6.6級地震前后區(qū)域重力場的短期前兆變化特征，以地震前后相鄰2期的觀測資料為時間基準，分別繪制了相鄰2期（0.5a尺度）區(qū)域重力場的差分圖像（圖3），分析圖3可以發(fā)現(xiàn)以下現(xiàn)象：

（1）2012-05—2012-10期間，重力變化具有分區(qū)特征，其一，甘東南地區(qū)重力變化較為復(fù)雜，總體為正負相間的波動變化，岷縣地區(qū)存在+40×10-8ms-2的局部重力異常變化；其二，寧夏隆德至銀川一帶自南向北出現(xiàn)由負向正的趨勢性變化，并在固原、隆德地區(qū)出現(xiàn)-60×10-8ms-2向北變化的局部重力異常區(qū)及伴生的重力變化梯度帶；其三，河西走廊地區(qū)重力變化較為平緩，重力變化在±30×10-8ms-2內(nèi)。

（2）2012-10—2013-05期間，整個測區(qū)重力變化非常劇烈，具有以下3個特征：1）甘東南地區(qū)出現(xiàn)顯著的重力差異運動，臨夏地區(qū)出現(xiàn)+50×10-8ms-2重力異常變化，岷縣、天水地區(qū)出現(xiàn)-100×10-8ms-2以上的重力變化，重力變化最大差異運動達150×10-8ms-2以上，并沿岷縣—隆德形成重力變化高梯度帶，岷縣漳縣6.6級地震發(fā)生在岷縣重力變化異常區(qū)伴生的重力變化高梯度帶上；2）寧夏隆德至銀川一帶自南向北出現(xiàn)由正向負的趨勢性變化，表現(xiàn)出與上期反向的變化特征，并沿固原、隆德出現(xiàn)局部重力異常區(qū)及重力變化梯度帶；3）河西走廊地區(qū)重力變化仍較平緩，重力變化在±30×10-8ms-2內(nèi)。

（3）2013-05—2013-10期間，重力變化總體與上期變化相反，表現(xiàn)為岷縣漳縣6.6級地震后的反向恢復(fù)變化。1）甘東南的臨夏東部由上期的+50×10-8ms-2轉(zhuǎn)為-30×10-8ms-2的重力異常變化，天水地區(qū)由上期的-100×10-8ms-2轉(zhuǎn)為+100×10-8ms-2的重力異常變化，仍沿岷縣—隆德形成重力變化高梯度帶，岷縣漳縣6.6級地震發(fā)生在重力變化梯度帶上、重力變化零值線附近和等值線的拐彎部位；2）寧夏固原至銀川一帶重力變化較為平緩；3）河西走廊地區(qū)重力變化與上2期變化相反，在門源和祁連附近出現(xiàn)-30和-40×10-8ms-2兩個局部重力異常區(qū)，2013年9月甘肅肅南—青海門源交界5.1級地震發(fā)生在門源和天祝2個局部異常區(qū)過渡地帶的重力變化梯度帶上。

3 重力變化與岷縣漳縣6.6級地震

強烈地震常常以其猝不及防的突發(fā)性和巨大的破壞力給社會經(jīng)濟發(fā)展、人類生存安全和社會穩(wěn)定以及社會功能帶來嚴重的危害。這在當今人的生命第一、以人為本地建設(shè)更高水平的小康社會的形式下更突現(xiàn)了防震減災(zāi)的重要性、緊迫性及其特殊地位。地震預(yù)報是減輕地震災(zāi)害最有效、最直接的方法之一，因此分析總結(jié)強震前后的觀測資料就顯得特別重要（馬瑾等，1995，2000）。

分析1a尺度的重力場動態(tài)圖像（圖2）可以發(fā)現(xiàn)，2011—2012年自西南向東北出現(xiàn)由負向正逐漸增加的趨勢性重力變化及沿構(gòu)造斷裂出現(xiàn)重力變化高梯度，而且區(qū)域重力場異常變化形態(tài)與其布格重力異常的空間分布具有很大程度的相關(guān)性。布格重力異常自西南向東北逐漸增加，由青海西寧的-340×10-5ms-2急劇上升至甘肅民勤的-240×10-5ms-2（祝意青等，2012b），重力變化也表現(xiàn)為自西南向東北逐漸增加，由青海西寧的-30×10-8ms-2上升至甘肅民勤的70× 10-8ms-2，祁連山斷裂帶所在的NW向布格重力梯度帶，其形態(tài)與2011—2012年的重力變化梯度帶基本一致。已有研究表明，布格重力異?？臻g變化特征與地殼深部構(gòu)造和地震的分布規(guī)律有著內(nèi)在的聯(lián)系，重力場時空變化特征則與區(qū)域構(gòu)造活動和地震活動有著密切關(guān)系。祁連山地區(qū)重力異常變化的幾何形態(tài)與布格重力空間分布如此密切相關(guān)，較好地反映了岷縣漳縣6.6級地震前區(qū)域構(gòu)造應(yīng)力增強引起活動斷層物質(zhì)變遷和構(gòu)造變形，在地表產(chǎn)生大空間尺度重力場的有序性變化。2012—2013年甘東南出現(xiàn)較大范圍的局部重力異常及劇烈的重力差異變化，這是局部應(yīng)力集中釋放的表現(xiàn)，較好地突出了岷縣漳縣6.6級地震前震中附近的重力變化前兆信息。另外，區(qū)域重力場的變化對測區(qū)西部的河西走廊地區(qū)發(fā)生的2013年甘肅肅南-青海門源交界5.1級地震也有一定程度的反映。綜上所述，1a尺度的重力場動態(tài)圖像較好地反映了區(qū)域應(yīng)力增強引起的大尺度重力場有序性變化及沿祁連-海原大斷裂的構(gòu)造活動，臨震前震中區(qū)域出現(xiàn)局部重力異常及其伴生的重力變化高梯度帶。

分析地震前后區(qū)域重力場差分動態(tài)圖像（圖3）可以發(fā)現(xiàn)，2013年7月22日岷縣漳縣6.6級地震前后測區(qū)的重力變化具有一定的分區(qū)性，以103°E為界，測區(qū)東南的甘東南地區(qū)2012-05—10期間區(qū)域重力場主要正、負相間地波動性變化，并在岷縣漳縣地震的震中區(qū)附近出現(xiàn)+40×10-8ms-2的重力變化異常區(qū)及重力變化梯度帶；2012-10—2013-05期間，震中西北的臨夏附近地區(qū)持續(xù)重力正值變化，震中附近的岷縣地區(qū)重力變化由上期正值轉(zhuǎn)為負值變化，震中附近重力變化較低，重力變化達-100×10-8ms-2，并在震中附近出現(xiàn)NE向重力變化高梯度帶；2013-05—10期間區(qū)域重力場總體變化與上期反向，地震發(fā)生在NE向重力變化高梯度帶上、重力變化零值線附近和等值線的拐彎部位，這與以往的研究基本一致（祝意青等，2004，2009，2013）。重力場差分動態(tài)圖像較好地反映了岷縣漳縣6.6級地震前后震中附近的重力變化，是一個由“重力正值異常區(qū)向重力劇烈負值變化異常區(qū)以及沿NE向出現(xiàn)重力變化梯度帶向震后反向恢復(fù)變化”的演化過程。

已有研究表明，與特定地震構(gòu)造有關(guān)的一次強震的發(fā)生，會引起地下應(yīng)力的重新排列或分布，導(dǎo)致附近或相鄰斷裂（段）應(yīng)變積累的非線性加速，從而促使那里潛在的強震提前發(fā)生（顧功敘，1997；馬瑾等，2007）。2008年汶川8.0級地震發(fā)生在青藏高原東緣的龍門山斷裂帶上，地震破裂沿向NE方向擴展，與青藏高原東北緣深部殼、幔物質(zhì)運移有著深層次的物質(zhì)與能量的交換和動力作用（滕吉文等，2008）。汶川地震后，西秦嶺北緣斷裂帶地區(qū)出現(xiàn)的顯著重力變化，說明該地區(qū)受汶川大地震的影響顯著（祝意青等，2012a）。2013年的蘆山7.0級地震發(fā)生在NE向龍門山斷裂帶南段，與發(fā)生過2008年汶川8.0級地震破裂的龍門山斷裂帶中-北段相鄰（馬瑾等，2013）。2008年汶川8.0級和2013年蘆山7.0級地震的相繼發(fā)生可能對岷縣漳縣地震具有一定的促震作用（祝意青等，2013）。岷縣漳縣6.6地震發(fā)生在青藏高原東北緣，東昆侖斷裂和西秦嶺北緣斷裂是該地區(qū)幾何特征復(fù)雜多樣的構(gòu)造中2條主要的邊界控制斷裂（鄭文俊等，2013）。受西秦嶺北緣斷裂帶向南側(cè)的擴展和青藏高原向NE擴展過程中東昆侖斷裂帶的NE向擠壓作用的共同影響，岷縣漳縣震中附近出現(xiàn)NE向重力變化高梯度帶，地震發(fā)生在NE向重力變化高梯度帶上、重力變化零值線附近和等值線的拐彎部位。

4 結(jié)論與討論

強震大都發(fā)生在活動板塊的邊界或活動斷裂帶上，2013年7月22日甘肅岷縣漳縣6.6級強震發(fā)生在臨潭-宕昌斷裂帶上。構(gòu)造活動斷裂帶由于其差異運動強烈而構(gòu)造變形非連續(xù)性最強，最有利于應(yīng)力的高度積累而孕育強震。已有研究結(jié)果表明（Chen，1979；顧功敘等，1997；Kuo，1999；祝意青等，2009；梁偉鋒等，2013），流動重力測量反映的是區(qū)域重力場的非潮汐變化，復(fù)測結(jié)果除了能較好地反映地殼厚度、密度空間差異變化的信息外，還包含由于構(gòu)造運動／變形引起的深部物質(zhì)遷移等信息，與地震發(fā)生有關(guān)的重力場變化是分布在地殼所有深度上的地下流體變化的響應(yīng)。由于流動重力觀測主要反映的是深部地下流體運移與變遷的信息，因此干擾較少，資料可信度較高。本文主要從重力場動態(tài)變化的角度，分析研究岷縣漳縣6.6級地震前后青藏高原東北緣的重力場動態(tài)變化特征，獲得了以下主要認識：

（1）青藏高原東北緣重力網(wǎng)的流動重力觀測資料顯示，2013年岷縣漳縣6.6級地震前出現(xiàn)了較好的中期前兆性變化圖像，即區(qū)域重力場變化先呈現(xiàn)大尺度空間范圍的有序性及與祁連-海原斷裂構(gòu)造帶走向基本一致的重力變化高梯度帶（圖2a）、后呈現(xiàn)震中附近特征性異常及與地震孕育發(fā)生有關(guān)的局部重力異常區(qū)（圖2b），2008汶川地震和2013蘆山地震之前也出現(xiàn)了這種類似現(xiàn)象（祝意青等，2010，2013）。岷縣漳縣地震發(fā)生在重力變化異常區(qū)（反映能量的積累與釋放）及NE向的重力變化高梯度帶（反映有利于地震破裂的部位）上。

（2）重力場差分動態(tài)圖像較好地反映了岷縣漳縣6.6級地震前后震中附近的重力變化，是一個由“重力正值異常區(qū)→重力劇烈負值變化異常區(qū)及沿NE向出現(xiàn)重力變化梯度帶→震后反向恢復(fù)變化”的演化過程，較好地反映了岷縣漳縣6.6級地震震中附近的重力異常變化。

（3）流動重力資料對2013年岷縣漳縣6.6級強震震中地點的判定，進一步證實區(qū)域重力場觀測對未來強震震中位置的判定具有獨到的優(yōu)勢。強震易發(fā)生在重力變化正、負異常區(qū)過渡的高梯度帶上，并須結(jié)合地震構(gòu)造活動情況。這是由于重力變化高梯度帶是物質(zhì)密度增加與減少的過渡地帶，該處產(chǎn)生的物質(zhì)增減差異運動劇烈，易產(chǎn)生剪應(yīng)力而首先破裂，從而誘發(fā)地震。岷縣漳縣6.6級地震發(fā)生在NE向重力變化高梯度帶上、重力變化零值線附近和等值線的拐彎部位。

（4）2013年7月岷縣漳縣6.6級地震后，重力場總體呈現(xiàn)反向恢復(fù)變化過程中，在測區(qū)西部的門源和祁連附近出現(xiàn)-30×10-8ms-2和-40×10-8ms-2兩個局部重力異常區(qū)。2013年9月甘肅肅南-青海門源交界5.1級地震發(fā)生在門源和天祝2個局部異常區(qū)過渡地帶的重力變化梯度帶上，重力變化對肅南-青海門源交界5.1級地震也有較好的反映。

（5）重力異常的分析表明，祁連山中東段的甘青交界地區(qū)重力場在時間演化上變化強烈，空間上差異變化劇烈，因此，存在中-長期強震的危險背景。

致謝 馬瑾院士一直關(guān)心地震預(yù)測預(yù)報工作。在與馬院士的交往中，深深感受到她的平易近人和對我們研究工作的關(guān)心支持。在馬院士80壽辰之際，謹以此文表達深深的敬意！在本課題的研究過程中及本文的資料分析處理過程中，與中國地震局地震研究所李輝研究員進行過多次有益的討論，審稿專家提出了有寶貴的意見，在此一并表示感謝。

顧功敘，Kuo JT，劉克人，等.1997.中國京津唐張地區(qū)時間上連續(xù)的重力變化與地震的孕育和發(fā)生［J］.科學通報，42（18）：1919—1930.

GU Gong-xu，Kuo J T，LIU Ke-ren，et al．1997.Seismogenesis and occurrence of earthquakes as observed by temporally continuous gravity variations in China［J］.Chinese Science Bulletin，42（18）：1919—1930（in Chinese）.

李輝，付廣裕，孫少安，等.2000.滇西地區(qū)重力場動態(tài)變化計算［J］.地殼形變與地震，20（1）：60—66.

LIHui，F(xiàn)U Guang-yu，SUN Shao-an，et al．2000.Computation on dynamic gravity changes in the western area of Yunnan Province［J］.Crustal Deformation and Earthquake，20（1）：60—66（in Chinese）.

梁偉鋒，劉芳，徐云馬，等.2013.青藏高原東緣重力觀測及對蘆山M 7.0地震的反映［J］.地震工程學報，35（2）：266—271.

LIANGWei-feng，LIU Fang，XU Yun-ma，et al．2013.Gravity observations along the eastern margin of the Tibetan plateau and an app lication to the Lushan MS7.0 earthquake［J］.Earthquake Science，26（3-4）：251—257.

馬瑾，馬勝利，劉力強.1995.地震前異常的階段性及其空間分布特征［J］.地震地質(zhì)，17（4）：363—371.

MA Jin，MA Sheng-li，LIU Li-qiang，et al．1995.The stages of anomalies before an earthquake and characteristics of the spatial distribution［J］.Seismology and Geology，17（4）：363—371（in Chinese）.

馬瑾，馬勝利，劉力強，等.2000.交叉斷層的交替活動與塊體運動的實驗研究［J］.地震地質(zhì)，22（1）：65—73.

MA Jin，MA Sheng-li，LIU Li-qiang，et al．2000.Experimental study on alternate slip of intersecting faults and block movement［J］.Seismology and Geology，22（1）：65—73（in Chinese）.

馬瑾，劉力強，劉培洵，等.2007.斷層失穩(wěn)錯動熱場前兆模式：雁列斷層的實驗研究［J］.地球物理學報，50（4）：l141—1149.

MA Jin，LIU Li-qiang，LIU Pei-xun，et al．2007.Thermal precursory pattern of fault unstable sliding：An experimental study of en echelon faults［J］.Chinese Journal of Geophysics，50（4）：1141—1149（in Chinese）.

馬瑾，劉培洵，劉遠征.2013.地震活動時空演化中看到的龍門山斷裂帶地震孕育的幾個現(xiàn)象［J］.地震地質(zhì)，35（3）：461—471.doi：10.3969／j.issn.0253-4967.2013.03.001.

MA Jin，LIU Pei-xun，LIU Yuan-zheng.2013.Features of seismogenic process of the Longmenshan Fault zone derived from analysis on the temporal-spatial evolution of earthquakes［J］.Seismology and Geology，35（3）：461—471（in Chinese）.

滕吉文，白登海，楊輝，等.2008.2008年汶川MS8.0地震發(fā)生的深層過程和動力學響應(yīng)［J］.地球物理學報，51（5）：1385—1402.

TENG Ji-wen，BAI Deng-hai，YANG Hui，et al．2008.Deep processes and dynamic responses associated with the Wenchuan MS8.0 earthquake of 2008［J］.Chinese JGeophys，51（5）：1385—1402（in Chinese）.

鄭文俊，閔偉，何文貴，等.2013.2013年甘肅岷縣漳縣6.6級地震震害分布特征及發(fā)震構(gòu)造分析［J］.地震地質(zhì)，35（3）：604—615.doi：10.3969／j.issn.0253-4967.2013.03.014.

ZHENGWen-jun，M INWei，HEWen-gui，et al．2013.Distribution of the related disaster and the causative tectonic of the Minxian-Zhangxian MS6.6 earthquake on July 22，2013，Gansu，China［J］.Seismology and Geology，35（3）：604—615（in Chinese）.

祝意青，李輝，朱桂芝，等.2004.青藏塊體東北緣重力場演化與地震活動［J］.地震學報，26（增）：71—78.

ZHU Yi-qing，LIHui，ZHU Gui-zhi，et al．2004.Gravity evolution and earthquake activities of the northeastern edgeof Qinghai-Tibet block［J］.Acta Seismologica Sinica，26（Suppl）：71—78（in Chinese）.

祝意青，胡斌，朱桂芝，等.2005.民樂6.1、岷縣5.2級地震前區(qū)域重力場變化研究［J］.大地測量與地球動力學，25（1）：24—29.

ZHU Yi-qing，HU Bin，ZHU Gui-zhi，et al．2005.Research on the variation of gravity field before Minle MS6.1 and M inxian MS5.2 earthquakes［J］.Journal of Geodesy and Geodynam ics，25（1）：24—29（in Chinese）.

祝意青，徐云馬，呂弋培，等.2009.龍門山斷裂帶重力變化與汶川8.0級地震關(guān)系研究［J］.地球物理學報，52（10）：2538—2546.doi：10.3969／j.issn.0001-5733.2009.10.012.

ZHU Yi-qing，XU Yun-ma，LU Yi-pei，et al．2009.Relations between gravity variation of Longmenshan Fault zone and Wenchuan MS8.0 earthquake［J］.Chinese JGeophys，52（10）：2538—2546（in Chinese）.

祝意青，梁偉鋒，徐云馬，等.2010.汶川MS8.0地震前后的重力場動態(tài)變化［J］.地震學報，32（6）：633—640.

ZHU Yi-qing，LIANG Wei-feng，XU Yun-ma，et al．2010.Dynamic variation of gravity field before and after Wenchuan MS8.0 earthquake［J］.Acta Seismologica Sinica，32（6）：633—640（in Chinese）.

祝意青，劉芳，付廣裕，等.2012a.汶川地震前后青藏高原東北緣重力場動態(tài)變化研究［J］.地震，32（2）：88—94.

ZHU Yi-qing，LIU Fang，F(xiàn)U Guang-yu，et al．2012a.Study on gravity change before and after Wenchuan MS8.0 earthquake in the northeastern edge of Qinghai-Tibet block［J］.Earthquake，32（2）：88—94（in Chinese）.

祝意青，梁偉鋒，陳石，等.2012b.青藏高原東北緣重力變化機理研究［J］.大地測量與地球動力學，32（3）：1—6.

ZHU Yi-qing，LIANG Wei-feng，CHEN Shi，et al．2012b.Study on mechanism of gravity field change in northeastern margin of Qinghai-Tibet plateau［J］.Journal of Geodesy and Geodynam ics，32（3）：1—6（in Chinese）.

祝意青，聞學澤，孫和平，等.2013.2013年四川蘆山MS7.0地震前的重力變化［J］.地球物理學報，56（6）：1887—1894.doi：10.6038／cjg20130611.

ZHU Yi-qing，WEN Xue-ze，SUN He-ping，et al．2013.Gravity changes before the Lushan，Sichuan，MS=7.0 earthquake of 2013［J］.Chinese JGeophys，56（6）：1887—1894（in Chinese）.

Chen Y T，Gu H D，Lu Z X.1979.Variations of gravity before and after the Haicheng earthquake，1975 and the Tangshan earthquake，1976［J］.Physics of the Earth and Planetary Interiors，18（4）：330—338.

Kuo JT，Zheng JH，Song SH.1999.Determ ination of earthquake epicentroids by inversion of gravity variation data in the BTTZ region，China［J］.Tectonophysics，312（2）：267—281.

Zhu Y Q，Zhan F B，Zhou J，et al．2010.Gravity measurements and their variations before the 2008 Wenchuan earthquake［J］.Bulletin of the Seismological Society of America，100（5B）：2815—282.

DYNAM IC VARIATION OF GRAVITY FIELD BEFORE AND AFTER THE M INXIAN-ZHANGXIAN MS6.6 EARTHQUAKE
ON JULY 22，2013，GANSU，CH INA

ZHU Yi-qing1，2）ZHAO Yun-feng1）LITie-ming3）LIANGWei-feng1）XU Yun-ma1）GUO Shu-song1）

1）Second Crust Monitoring and Application Center，China Earthquake Adiministration，Xi′an 710054，China
2）State Key Laboratory of Geodesy and Earth's Dynam ics，Institute of Geodesy and Geophysics，Chinese Academy of Sciences，Wuhan 430077，China
3）Institute of Geology，China Earthquake Administration，Beijing 100029，China

Using the observation data of relative gravity measurements from the northeastern edge of the Tibetan p lateau，we have systematically analyzed the spatial-temporal variation of the regional gravity field and their relation to the occurrence of the MS6.6 Minxian-Zhangxian earthquake on July 22，2013.Our research mainly shows that：1）Spatially the anomaly change of the regional gravity field is closely related to the major fault zones in the region，suggesting that the tectonic activities or deformation which happened along the fault zones（or segments）during 2011 to 2013 would have resulted in the variation in the ground gravity observations.2）Before the MS6.6 Minxian-Zhangxiang earthquake，regional gravity anomaly firstly appeared in the whole study region，and then local gravity anomalies and high gradient zone of gravity change appeared at and near the potential source area. Among the areas of gravity anomalies，the anomaly variations near Linxia and Minxian were up to and over 150×10-8ms-2，indicating possibly the tectonic movement or stress enhancement in the study region and around the potential source area，which could be related to the preparation and occurrence of the earthquake.3）The dynamic change patterns and the differential evolution patterns of gravity variation of one year scale suggest that significant anomalies of gravity variations in the mobile gravity measurements appeared indeed in the last two years of the preparation process of the MS6.6 Minxian-Zhangxian earthquake.The earthquake occurred on the NE-trending high gradient zone near the zero curve and the turning part of contours of gravity change.4）To a certain degree，wemade amediumterm forecast before the M inxian-Zhangxian MS6.6 earthquake，especially a forecast for the locality of the earthquake based on the anomalous gravity variations.

northeastern edge of the Tibetan plateau，Minxian-Zhangxian MS6.6 earthquake，gravity observation，gravity variation，tectonic activity

P315.72+6

A

0253-4967（2014）03-0667-10

祝意青，男，1962年生，1983年畢業(yè)于武漢測繪學院大地測量系，獲學士學位，研究員，主要從事重力測量與地球動力學及地震分析預(yù)報研究，E-mail：zhuyiqing＠163.com。

10.3969／j.issn.0253-4967.2014.03.010

2013-12-31收稿，2014-04-28改回。

中國地震局地震行業(yè)科研專項（201308009）、國家自然科學基金（41274083，41374026）和大地測量與地球動力學國家重點實驗室開放基金（SKLGED201331E）共同資助。①中國地震局第二監(jiān)測中心2011—2013年度地震趨勢研究報告。