鄭　兵，王青華，羅海元，馬伶俐，蒲曉霞，李菲菲，蘇　琴

(1.四川省地震局 測繪工程院，四川 雅安 625000；2.云南省地震局，云南 昆明 650300)

2012年寧蒗—鹽源5.7級地震前的重力變化*

鄭兵1，王青華2，羅海元1，馬伶俐1，蒲曉霞1，李菲菲1，蘇琴1

(1.四川省地震局 測繪工程院，四川 雅安 625000；2.云南省地震局，云南 昆明 650300)

摘要：利用2010年7月～2012年3月四川省和云南省的流動重力觀測數(shù)據(jù)，對2012年6月24日寧蒗—鹽源5.7級地震孕育和發(fā)生過程中的重力場動態(tài)演化情況進行分析。結果顯示：地震孕育階段震中區(qū)域重力場出現(xiàn)大面積的重力正值變化高值異常區(qū)域，并形成密集的正值變化梯度帶；震前震中區(qū)域呈重力正變化低值狀態(tài)，處于重力變化0值線附近且等值線稀疏；震前與發(fā)震有關的斷裂帶的跨斷層重力測段段差值變化明顯。

關鍵詞：流動重力測量；重力場；川滇地區(qū)；寧蒗—鹽源5.7級地震

0引言

流動重力測量反映的是區(qū)域重力場的非潮汐變化信息，地殼內部的物質遷移、地殼構造和地震的形成過程等都可以在流動重力復測結果中反映出來。重力變化與構造應力之間有著密切的聯(lián)系(向文，李輝，1999；蘇有錦，泰嘉政，2001)，區(qū)域構造應力場發(fā)生強弱的變化，重力場圖像也會相應地顯示出時空變化，而重力場動態(tài)變化圖像一直是地震監(jiān)測預報研究的基本信息源(陳運泰等，1980；賈民育等，1995；李瑞浩等，1997；李輝等，2009)。通過對重力場動態(tài)變化的深入分析，有可能捕獲某些地震前兆的有效信息，為地震的預測預報進行有效的探索。

2012年6月24日在云南省麗江市寧蒗彝族自治縣、四川省涼山彝族自治州鹽源縣交界地區(qū)(27.7°N，100.7°E)發(fā)生了5.7級地震，震源深度約11 km，屬主—余震型地震，造成了一定的人員傷亡和財產損失。四川和云南的流動重力資料顯示重力場在該地震前有異常變化，本文利用兩省的流動重力資料研究震前區(qū)域重力場時空動態(tài)演化特征，探討重力場的異常變化與地震活動的規(guī)律，對地震孕育過程中捕捉地震前兆進行監(jiān)測預報具有一定的實際意義。

1測區(qū)概況及資料處理

筆者采用2010年7月～2012年3月的流動重力觀測數(shù)據(jù)，對數(shù)據(jù)資料的處理利用中國地震局實用化攻關推廣的應用軟件LGADJ，統(tǒng)一起算基準點，采用絕對重力控制下的相對重力聯(lián)測的弱基準平差，絕對重力點采用四川郫縣、云南下關和昆明3個絕對重力點(圖1)作為起算基準點(3個絕對重力點相對于2011年的重力變化分別為-6.1×10 nm·s-2、-1.8×10 nm·s-2、5.8×10 nm·s-2，重力變化較小，2012年沒有測量，采用的是2011年的測量結果)，對川滇流動重力聯(lián)測資料進行整體平差，得到統(tǒng)一起算基準的絕對重力變化(李輝等，2000；申重陽等，2011)。在了解測量儀器的后驗中誤差后，依據(jù)后驗中誤差在平差計算時合理調整各臺儀器的先驗方差以得到最佳解算結果。川滇的數(shù)據(jù)是分測區(qū)觀測的，2011年上半年在四川省攀枝花市以渡口、河門口和格平8G共3個公共接點將兩個測區(qū)的數(shù)據(jù)聯(lián)為一個整體。由于云南和四川兩個測網(wǎng)單獨計算得出的公共點重力值差別較大，2011～2012年上半年3期云南和四川公共點點值數(shù)值差值渡口：-402 nm·s-2、623 nm·s-2、1 524 nm·s-2；格平8G：-581 nm·s-2、647 nm·s-2、1 186 nm·s-2；河門口：67 nm·s-2、648 nm·s-2、1 335 nm·s-2，所以從2011年起以公共點將云南、四川兩個測區(qū)數(shù)據(jù)連接進行整體平差，減少誤差，得到較好的處理結果。解算出4期數(shù)據(jù)的點值精度分別是151 nm·s-2、110 nm·s-2、96nm·s-2、73 nm·s-2，數(shù)據(jù)的可信度較高，系統(tǒng)誤差較小，成果質量較好，能夠用于分析研究。

2區(qū)域重力場的演化特征

2.1重力場動態(tài)變化

實踐表明，分析重力場的動態(tài)變化圖像能夠捕獲地震前后重力場變化特征。將處理后的2010年7月～2012年3月的重力資料做出差分動態(tài)變化等值線圖進行分析研究，如圖2所示。雖然寧蒗—鹽源5.7級地震的震中區(qū)域沒有流動重力實測點位通過，但川滇測區(qū)的流動重力點位鄰近該區(qū)域，重力變化等值線圖仍然能夠較為真實的反映出震前的重力場時空動態(tài)演化特征。

(1)2010年7月～2011年2月

該時段川滇測區(qū)雖然沒有公共點，有一定的誤差，但是以絕對重力點為起算基準，算出的結果能夠反映出重力場變化情況。四川和滇東測區(qū)出現(xiàn)了較大面積的重力負值異常區(qū)域，且變化數(shù)值較大，最大值約-1 200 nm·s-2，異常區(qū)域主要分布在安寧河—則木河斷裂帶區(qū)域，等值線較為密集，這與該時段該區(qū)域發(fā)生較多的地震*四川省地震局—震情速遞.http：//www.scdzj.gov.cn/xwzx/zqsb.是相對應的。四川的寶興—石棉區(qū)域和德陽—成都區(qū)域呈重力正值變化，形成了面積較小、變化數(shù)值較小的重力正值閉合等值線區(qū)域；麗江—攀枝花—昆明區(qū)域西南部的滇西測區(qū)呈重力正值變化，變化數(shù)值不大，等值線較為稀疏。寧蒗—鹽源區(qū)域和寧蒗—鹽源5.7級地震震中位于川滇菱形塊體中部南側，是重力負值變化區(qū)域，西邊鄰近重力變化0等值線，等值線相對密集，等值線負值變化數(shù)值向鹽源方向逐步增大，在鹽源縣城附近重力變化值約-600 nm·s-2，屬于重力負值異常區(qū)域(圖2a)。

構成過腔的不同音樂材料，無意間為過腔構筑起了多種不同的結構樣式，使得過腔充滿了無窮的活力，更重要的是，將過腔的結構與字腔本身固有的頭、腹、尾腔的結構結合起來形成的字腔+（或 0）過腔①所謂“（或0）過腔”是說，不一定每個字腔都要帶過腔，過腔的出現(xiàn)也根據(jù)藝術的需要。即使是0過腔，對每一個字腔來說，這個位置卻是不可替代的。就像劇場中的座位，即使沒人坐，但這個位置照樣存在一樣。這種結構，不僅是對傳統(tǒng)曲牌音樂結構的重大突破，更是世所未見的、嶄新的昆曲曲牌唱調的基本結構。

(2)2011年2月～2011年8月

川滇測區(qū)的重力變化情況與上一個時段相比，發(fā)生了一個較大的逆轉。四川和滇東測區(qū)在攀枝花—鹽源區(qū)域之間有一小塊區(qū)域繼續(xù)保持數(shù)值較小的重力變化負值，其它區(qū)域幾乎全部是重力變化正值區(qū)域，出現(xiàn)了多處重力正值變化異常：在石棉—康定區(qū)域是重力正值閉合等值線，出現(xiàn)數(shù)值較小的重力變化正值異常區(qū)域，等值線分布稀疏；汶川地震余震區(qū)出現(xiàn)重力正值變化異常區(qū)域，等值線分布較為密集，在都江堰—北川—青川方向形成重力等值線變化梯度帶，青川沒有實測的流動重力點通過，但從等值線在北川地區(qū)變密的趨勢來看，青川地區(qū)的重力變化應該是較大的，該梯度帶及其附近發(fā)生了多個地震，震級較大的是2011年11月1日青川文縣5.4級地震和2011年12月26日彭州4.8級地震；滇東測區(qū)重力正值變化較小，沒有出現(xiàn)重力正值變化異常區(qū)域。滇西測區(qū)則由上一個時段的重力正值變化轉為重力負值變化，變化數(shù)值不大，等值線分布稀疏。寧蒗—鹽源區(qū)域的重力變化值由上一個時段的負值變化轉為正值變化，變化幅度較大，在鹽源縣城西南區(qū)域正值變化數(shù)值最大，達800 nm·s-2，出現(xiàn)了等值線較為密集的高值異常梯度帶，該梯度帶延伸到了西昌、普格、寧南和雷波地區(qū)，其中寧南—雷波梯度帶鄰近2012年9月7日彝良5.7、5.6級地震震中區(qū)域，可能是該地震的活動對此梯度帶的形成施加了影響。寧蒗—鹽源—西昌—普格—寧南梯度帶緊靠重力變化0值線，0值線在鹽源地區(qū)呈圓弧形的轉折，該梯度帶及其附近地區(qū)寧蒗—鹽源5.7、5.6級地震的震中在此時段等值線數(shù)值相對較小，密度相對稀疏(由于震中位置沒有實測點位通過)，但是處于高值異常梯度帶上，說明這個時段內寧蒗—鹽源區(qū)域的構造斷裂帶活動較大，使該地區(qū)的地殼結構狀態(tài)和密度發(fā)生變化，直接表現(xiàn)在該地區(qū)的重力變化在短期之內較大，同時也說明該地區(qū)的重力正值變化始于2011年下半年(圖2b)。

(3)2011年8月～2012年3月

這個時段四川測區(qū)的康定—瀘定小塊區(qū)域繼續(xù)重力正值變化并形成閉合等值線區(qū)域，康定地區(qū)的等值線分布較為密集；理縣—汶川—都江堰區(qū)域小塊區(qū)域是閉合的，是數(shù)值較小的重力正值變化等值線區(qū)域；川滇測區(qū)交界的中甸—木里—普格—米易—東川—武定—永仁—云龍所包圍的區(qū)域呈數(shù)值較小的重力正值變化且等值線較為稀疏，曲靖以南、路南以東區(qū)域也是低值重力正值變化區(qū)域。整個川滇測區(qū)重力正值變化區(qū)域面積較小且分散，主要分布在川滇交界地區(qū)和滇東地區(qū)，重力負值變化區(qū)域占主要面積，圖形變化情況與圖2b大致相反，類似于圖2a，只是重力負值變化數(shù)值較小。寧蒗—鹽源地區(qū)繼續(xù)上一個時段的重力正值變化，變化的數(shù)值變小，約200 nm·s-2，處于重力正負變化值交界區(qū)域附近，分布著稀疏的重力正值變化等值線，寧蒗—鹽源5.7級地震的震中區(qū)域鄰近重力變化0值線，其重力變化數(shù)值約100 nm·s-2，這說明在臨震前約3個月寧蒗—鹽源地區(qū)的重力變化雖然持續(xù)增加，但是增加量已經變小(圖2c)。

(4)2011年2月～2012年3月

由于川滇測區(qū)重力正值變化是從2011年2月開始的，所以從1 a的尺度來分析震前的重力變化情況。從圖2d可以看出，從2011年初到2012年初，川滇測區(qū)的重力主要以正值變化為主，在馬爾康—寶興—康定北部—道孚區(qū)域、北川—都江堰—成都東側和滇西區(qū)域是重力負值變化區(qū)域，且數(shù)值不大。在重力正值變化區(qū)域出現(xiàn)異常的是康定地區(qū)和西昌—鹽源區(qū)域：康定地區(qū)重力正值變化等值線密集，形成閉合區(qū)域，靠近0值線；西昌—鹽源區(qū)域重力變化正值等值線密集，形成閉合區(qū)域，等值線密集方向向寧蒗方向延伸，寧蒗—鹽源5.7級地震的震中區(qū)域位于西昌—鹽源區(qū)域附近，鄰近重力等值線0值線，等值線分布稀疏，重力變化數(shù)值較小，約100 nm·s-2，但這是由其附近川滇重力實測點重力變化輻射過去的，從重力變化等值線的分布趨勢來看，如果該處有實測點位，等值線的分布也會比較密集。這說明在2012年上半年川滇測區(qū)雖然以重力負值變化為主，但是長時間尺度來看還是重力變化增加；寧蒗—鹽源地區(qū)在震前重力正值變化增加速度變慢，形成了重力變化等值線梯度帶，聚集了較多的能量。

2.2重力段差時序變化

重力段差時序變化能較好的反映跨斷裂地段的重力動態(tài)變化。雖然寧蒗—鹽源5.7級地震震中區(qū)域沒有重力實測點位通過，但在震中西北約9 km有一條小活動斷裂延伸到滇西測區(qū)，太安—白漢場新測段跨過該斷裂。震中附近的斷裂帶還有：東北方向約30 km的博科—木里斷裂，東南方向約33 km的小金河斷裂及一些小活動斷裂，西南約45 km和58 km的洪門口斷裂和大具斷裂，由于這些斷裂帶沒有布設流動重力點，不能直接研究震前各斷裂帶的活動變化，但是相鄰斷裂帶的活動是互相影響的，所以在川滇測區(qū)各選取與震中附近斷裂帶相鄰的活動斷裂帶兩側的重力點位，通過研究這些點位的段差變化來研究震中附近斷裂帶的重力動態(tài)變化。鹽源弧形斷裂西南段與博科—木里斷裂和小金河斷裂接壤，其跨斷層測段是鹽源—黃草；洱源—鶴慶斷裂與小金河斷裂接壤，其跨斷層測段是龍山—麗江新。選取上述點位測段和與這些點位相鄰而又不跨斷裂帶的流動重力點位測段(黃草—格薩拉、麗江新—太安、飛崖—龍山、麗江臺—麗江新)的段差求出以起始變化值為0的相對變化量，以下盤相對上盤的變化為基準作出重力段差變化圖(圖3)進行分析研究。

從圖3可以看出，2010-07～2011-08時段跨斷裂段差變化曲線的變化方向不一致，其變化值沒有明顯變化，但在震前2011-08～2012-03時段，段差值變化的方向卻一致增大，太安—白漢場新測段的變化值有明顯變化，約660 nm·s-2，其它段差變化值則較小。從圖3b可以看出，在2010-07～2011-02時段，非跨斷層變化段差變化值方向雖不同，但是變化數(shù)值都較?。?011-02～2011-08段差變化值朝減小方向發(fā)展，飛崖—龍山段差變化值相對較大，但不超過300 nm·s-2；在2011-08～2012-03時段，各段差變化繼續(xù)向減小方向發(fā)展，麗江新—太安測段變化較為顯著，約500 nm·s-2，而其余測段的段差值變化量相對較小。由于地震發(fā)生之前，橫跨發(fā)震構造的測段，重力段差值可能發(fā)生明顯的變化，與發(fā)震構造無關的測段，重力變化值不大(吳兵，馬麗，1982)。寧蒗—鹽源5.7級地震震中西北小活動斷裂帶的太安—白漢場新測段在震前重力段差值變化明顯，其相鄰測段麗江新—太安測段段差雖然沒有跨過斷層，但受太安點影響，段差變化也較為顯著，而同處斷裂帶一側的非跨斷層測段麗江臺—麗江新測段段差值變化數(shù)值則變化相對較小，不超過300 nm·s-2，這說明該斷裂是與發(fā)震有關的構造，該斷裂帶和震中附近斷裂帶的活動對地震的孕育、形成和發(fā)生起到了很大的作用；其余鹽源弧形斷裂和洱源—鶴慶斷裂由于跨斷層和非跨斷層段差變化都不是特別顯著，所以其構造運動與該地震無關，反而該地震的孕震活動通過小金河等斷裂帶的活動對其構造運動產生了影響，引起了它們的重力場的變化，使這些跨斷裂測段段差值在震前一致增加，非跨斷裂測段段差值則在震前一致減小。

3討論與結論

重力變化與構造應力之間有著密切的聯(lián)系(向文，李輝，1999)，重力場變化能反映出區(qū)域應力場的微動態(tài)活動。雖然2012年6月24日寧蒗—鹽源5.7級地震震中位置沒有實測的流動重力點位，但是由于斷層活動是相互影響的，本文所采用的重力變化資料就是該地震對鄰近區(qū)域的流動重力產生的影響而反映出地震的孕育過程。從重力場動態(tài)變化圖像和段差時序變化可以看出，寧蒗—鹽源5.7級地震是震中附近斷裂帶活動的結果，而斷裂帶的構造活動又影響了鄰近的鹽源弧形斷裂、寧蒗斷裂、洱源—鶴慶斷裂和玉龍雪山東麓斷裂等的活動，具體表現(xiàn)就是這些斷裂帶附近的流動重力點位的重力場變化。距離震中較近的四川鹽源測區(qū)的鹽源弧形斷裂帶等活動受到的影響較大，使該地區(qū)地殼物質密度發(fā)生變化，從而反映出了震中重力的變化。2011年下半年在寧蒗—鹽源地區(qū)及其附近較大面積的區(qū)域形成了重力正值變化高值異常區(qū)域，重力變化等值線形成了密集的高梯度帶，構造活動聚集了大量的能量，滿足了地震孕育發(fā)生的應力積累條件(申重陽等，2011)，2012年上半年重力場正值變化區(qū)域大面積減少，重力正值高值異常區(qū)域消失，該地區(qū)的重力場正值變化雖然仍在增加，但是增加速度減小了許多，重力變化等值線稀疏，重力變化值減小到約200 nm·s-2，寧蒗—鹽源5.7級地震的震中區(qū)域重力變化數(shù)值約100 nm·s-2，鄰近重力變化0值線的轉折部位，這樣地震破裂部位的條件也具備了(申重陽等，2010)，如果在震中區(qū)域有實測重力點位通過，圖形效果會更明顯；寧蒗—鹽源5.7級地震前，鄰近震中區(qū)域流動重力測段段差值變化相對明顯，這也是地震孕育的跡象(吳兵，馬麗，1982)。重力正值變化顯著的梯度帶及其附近地區(qū)，往往在未來1～2年內有中強震發(fā)生(祝意青等，1999，2003；朱治國等，2012；梁偉鋒等，2012)，但是寧蒗—鹽源5.7級地震從能量積累到發(fā)震，經歷的時間約1 a，時間相對較短，該地震發(fā)生在非主流活動小斷裂附近，斷裂活動力度較大，使重力變化相對較大，在較短的時間積累了發(fā)震的能量并找到了能量爆發(fā)的突破口(震中)，這是一個由非主流活動小斷裂的構造活動引發(fā)的中強震。

從上面的討論分析可以看出，寧蒗—鹽源5.7級地震前的重力變化具有如下特點：

(1)寧蒗—鹽源5.7級地震是震區(qū)斷裂帶構造活動的結果，震前與發(fā)震有關的斷裂帶跨斷層重力測段段差值變化明顯。

(2)雖然寧蒗—鹽源5.7級地震的震中區(qū)域沒有實測流動重力點位，但是仍然較好地反映了地震從孕育到發(fā)生的過程中重力場變化情況：地震孕育階段震中區(qū)域重力場出現(xiàn)大面積的重力正值變化高值異常區(qū)域，形成密集的正值變化梯度帶；震前震中區(qū)域呈重力正變化低值狀態(tài)，處于重力變化0值線附近，且等值線稀疏。

(3)從重力場開始正值變化到發(fā)震經歷的時間約1 a。

通過對寧蒗—鹽源5.7級地震的重力場變化的分析和討論，可以看出重力場動態(tài)變化對研究地震活動與地質構造活動的關系具有重要的意義。川滇測區(qū)2012年上半年的流動重力觀測是在6月5日結束的，所以該期數(shù)據(jù)可視為寧蒗—鹽源5.7級地震震前數(shù)據(jù)，對其重力場變化的研究與總結為對中強震乃至更高震級的地震預測預報能夠提供的參考震例和經驗。

參考文獻：

陳運泰，顧浩鼎，盧造勛.1980.1975年海城地震與1976年唐山地震前后的重力變化[J].地震學報，2(1)：21-31.

賈民育，邢燦飛，孫少安.1995.滇西重力變化二維圖像及其與5級(MS)以上地震的關系[J].地殼形變與地震，15(3)：9-19.

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

李輝，申重陽，孫少安，等.2009.中國大陸近期重力場動態(tài)變化圖像[J].大地測量與地球動力學，29(3)：1-10.

李瑞浩，黃建梁，李輝，等.1997.唐山地震前后區(qū)域重力場變化機制[J].地震學報，19(4)：399-407.

梁偉鋒，祝意青，徐云馬，等.2012.文安5.1級地震前后不同時間尺度的重力變化分析[J].地震研究，35(3)：318-322.

申重陽，李輝，孫少安，等.2010.2008年于田MS7.3地震前重力場動態(tài)變化特征分析[J].大地測量與地球動力學，30(4)：1-7.

申重陽、談洪波、郝洪濤，等.2011.2009年姚安MS6.0地震重力場前兆變化機理[J].大地測量與地球動力學，31(2)：17-23.

蘇有錦，秦嘉政.2001.川滇地區(qū)強震活動與區(qū)域新構造運動的關系[J].中國地震，17(1)：24-34.

吳兵，馬麗.1982.斷層活動與重力變化關系的探討[J].地震，(2)：28-30.

向文，李輝.1999.重力場構造與構造應力場內在關系的理論研究[J].地殼形變與地震，19(1)：32-36.

朱治國，劉代芹，李杰，等.2012.2011年新疆伽師MS5.8地震重力場前兆異常特征[J].地震研究，35(2)：201-206.

祝意青，胡斌，張永志.1999.永登5.8級地震前后的重力場動態(tài)圖象特征研究[J].地殼形變與地震，19(1)：71-77.

祝意青，胡斌，朱桂芝，等.2003.河西地區(qū)重力場及其動態(tài)演化特征[J].大地測量與地球動力學，23(4)：44-48.

Using the flow gravity observation data in Sichuan and Yunnan provinces from Jul.，2010 to Mar.，2012，we analyzed the gravity field dynamic evolution in the process of preparation and occurrence of Ninglang-YanyuanMS5.7 earthquake on Jun.24，2012.The results showed that in the process of earthquake preparation，a large area of high value anomaly of gravity positive changes appeared and formed to the intensive gradient belt of gravity positive in the area of epicenter.Before the earthquake，the epicenter region appeared the low value anomaly of gravity positive changes，where was near to the 0 contour line of the gravity field and the sparse contour line，and the variation in segment differences of the gravity across the fault zone which induced the earthquake was obvious.

Key words：flow gravity measuring；gravity field；Sichuan-Yunnan area；Ninglang-YanyuanMS5.7 earthquake

*收稿日期：2014-08-14. 基金項目：四川省重點“項目高精度智能測量機器人跨斷層位移測量實驗”(LY1406)、“基于GPS、重力測量對龍門山斷裂帶蘆山剖面跟蹤分析”(LY1506)和地震科技星火計劃攻關項目“川滇地殼形變特征及在強震中短期預測中的應用研究(XH15037SX)共同資助.

中圖分類號：P315.7

文獻標識碼：A

文章編號：1000-0666(2015)02-0370-07

Gravity Variation before Ninglang-YanyuanMS5.7 Earthquake in 2012

ZHENG Bing1，WANG Qing-hua2，LUO Hai-yuan1，MA Ling-li1，PU Xiao-xia1，LI Fei-fei1，SU Qin1

(1.Surveying Engineering Institute，Earthquake Administration of Sichuan Province，Yaan 625000，Sichuan，China) (2.Earthquake Administration of Yunnan Province，Kunming 620300，Yunnan，China)

Abstract