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

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

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

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

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

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

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

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

0引言

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

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

1測區概況及資料處理

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

2區域重力場的演化特征

2.1重力場動態變化

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

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

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

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

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

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

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

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

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

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

2.2重力段差時序變化

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

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

3討論與結論

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

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

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

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

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

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

參考文獻：

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

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

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

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

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

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

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

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

蘇有錦，秦嘉政.2001.川滇地區強震活動與區域新構造運動的關系[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級地震前后的重力場動態圖象特征研究[J].地殼形變與地震，19(1)：71-77.

祝意青，胡斌，朱桂芝，等.2003.河西地區重力場及其動態演化特征[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