文安5.1級地震前后不同時間尺度的重力變化分析*

梁偉鋒，祝意青，徐云馬，劉 練，郭樹松，劉 芳

(中國地震局第二監測中心，陜西西安710054)

文安5.1級地震前后不同時間尺度的重力變化分析*

梁偉鋒，祝意青，徐云馬，劉 練，郭樹松，劉 芳

(中國地震局第二監測中心，陜西西安710054)

系統深入地分析了2005～2009年華北地區流動重力觀測資料，提取2006年河北省文安縣MS5.1地震前后0.5 a、1 a等不同時間尺度可靠的重力場動態變化信息，研究其與地震活動的關系。結果顯示:文安地震孕育階段，測區出現較大范圍的區域性重力異常，并產生與地震孕育有關的局部重力異常區及重力變化高梯度帶。重力場動態圖像較清晰地反映了區域重力場的有序性演化與地震活動。

文安5.1級地震;流動重力觀測;重力變化;重力場

0 引言

2006年7月4日11時56分，河北省文安縣發生MS5.1地震，震中位置 (38.9°N，116.3°E)，震源深度19.3 km。地震發生后，震中區周邊的北京、天津、山西以及山東濟南、淄博、濱州、河北石家莊、保定等大部分地區有明顯震感。河北文安縣位于北京正南方向，因此北京南部地區的震感比較強烈，但對北京沒有造成破壞。

文安地震發生在河北平原強震構造帶上。該構造帶北起盧龍、灤縣，向西南經唐山、寧河、武清、河間、深縣、辛集、寧晉、任縣、邯鄲至磁縣，全長540 km。歷史上該構造帶記錄到MS6.0以上強震19次。在強震構造帶的河間至寧河之間，是缺失6級以上地震最長的一段，地震能量釋放最不充分，文安位于缺震段的中南部。文安MS5.1地震發生在首都圈南部的文霸流動重力測網區域內 (圖1)，震中處于測網的中部，河北省地震局保定測量隊使用LCR-G型重力儀測量，地震前后均積累多期流動重力復測資料。王武星等 (2009)分析了2006年7月4日文安5.1級地震前后文霸流動重力測網的觀測結果，討論了垂直形變和淺層地下水位變化引起的重力改變，并對重力資料進行了高程校正。蘇樹朋等 (2009，2010)和張素靈等 (2009)對1998年1月10日張北6.2級、2002年4月22日邢臺5.4級、2006年7月4日文安5.1級3次地震震前的部分重力點值和震前0.5 a尺度的重力場空間變化進行了分析，但所用資料截止于2007年9月 (國家地震局科技監測司，1986)。本文使用《LGADJ》數據處理系統針對文霸流動重力測網在地震前后的多期流動重力觀測資料 (2005.03～2009.03)進行了統一標準的處理，并分析了不同的時間尺度的重力變化，研究文安地震孕育發生過程中重力場的有序性演化與地震活動。

圖1 文霸流動重力測網圖Fig.1 Flow gravity network in Wenba

1 不同時間尺度的重力變化

1.1 0.5 a尺度的重力變化

參考已有的研究成果，季節變化、地下水位變化等因素對重力變化的影響為:地下水位每變化 (上升或下降)1 m，重力變化在 (1.4～1.5)×10－8m·s－2之間，且基本呈線性函數關系 (何志堂等，2009)。根據國土資源部2010年5月發布的《我國主要城市和地區地下水水情通報》①國土資源部.2010.我國主要城市和地區地下水水情通報.，2008～2009年度華北地區的地下水位基本保持穩定，升降幅度在0.5 m以內，對重力變化的影響小于1×10－8m·s－2，所以地下水位變化對重力變化的影響可以忽略。

(1)震前相鄰兩期的的重力變化

震前以0.5 a尺度相鄰兩期的的重力變化表明，在2005.03～2005.09期間 (圖2a，圖中帶數字的線條為重力變化等值線，數值表示重力變化數值，單位為10－8m·s－2，實線為重力正值變化，虛線為重力負值變化，以下所有圖相同)，震中附近重力變化平緩，在震中西部有 －10×10－8m·s－2的重力負值區;但在2005.09～2006.05期間 (圖2b)，在震中附近的易縣到文安縣之間形成1個北東向的重力變化高梯度帶，最小重力差值高達90×10－8m·s－2，靠近震中附近重力急劇增加，在文安和任丘之間形成1個重力變化正值高值區，最大重力變化值高達70×10－8m·s－2。此為地震孕育期，是地震前斷層處于閉鎖狀態，能量積累階段，地下物質變化劇烈，重力變化相應強烈。

(2)地震期間相鄰兩期的重力變化

在地震發生的2006.05～2006.12期間 (圖2c)，重力場的變化格局發生了改變，2005.09～2006.05期間北負南正的重力變化格局變為了東負西正，圖2b所示文安和任丘之間重力正值變化高值區變為負值，地震發生在重力變化零值線附近。地震發生的過程中釋放能量，重力變化變小。

圖2 文安MS5.1地震前后0.5 a尺度的重力變化(a)2005.03～2005.09;(b)2005.09～2006.05;(c)2006.05～2006.12;(d)2006.12～2007.03Fig.2 Gravity variation in half a year scale before and after Wen'an MS5.1 earthquake

(3)震后相鄰兩期的重力變化

地震發生后的2006.12～2007.03期間 (圖2d)，文安地震震中東西兩側重力變化反向，以震中為界，2006.05～2006.12期間“東負西正”的格局變為了“東正西負”，震中處的重力變化仍在零值線附近，此屬于震后調整階段。

圖3 文安MS5.1地震后0.5 a尺度的重力變化(a)2007.03～2007.09;(b)2007.09～2008.03;(c)2008.03～2008.11;(d)2008.11～2009.03Fig.3 Gravity variation in half a year scale after Wen'an MS5.1 earthquake

華北地區的地殼介質處于一致性較好的統一的應力場中:最大壓應力主軸方位為北東東—南西西向，最小壓應力主軸方位為北北西—南南東向，而且都接近于水平向;兩個最大剪切應力方向分別為北北東—南南西和北西西—南東東向。受近東西向水平擠壓力為主的作用，震后相鄰兩期 (2007.03 ～2007.09、2007.09 ～2008.03、2008.03～2008.11、2008.11～2009.03)(圖2) 的重力變化表現為持續調整階段，隨著能量的釋放結束，地下物質的運動趨于平靜，重力變化由緩慢的調整也逐漸變為平靜，慢慢恢復到地震前的狀態:2007.03～2007.09期間 (圖3a)，重力變化負值區西移，在易縣—保定—饒陽一線形成一南北條帶狀的重力變化負值區，重力負值變化最大－30×10－8m·s－2，震中附近及其它地區重力變化平緩;2007.09～2008.03期間 (圖 3b)，2007.03～2007.09期間的南北條帶狀的重力變化負值區消失，安國—饒陽附近由上期的負值重力變化最大轉為正值最大，震中附近有－10×10－8m·s－2的重力變化;2008.03～2008.11(圖 3c)和2008.11－2009.03(圖3d)期間，震中附近重力變化趨于平靜，接近震前的變化水平。

1.2 1 a尺度的重力變化

(1)震前相鄰兩年的的重力變化

從震前相鄰兩年的重力變化可看出，2004.03～2005.03期間 (圖4a)，屬于震前能量積累，震中附近重力變化平穩;而地震發生前的2005.03～2006.05期間 (圖4b)，震中附近的地下物質變化劇烈，能量積累加劇，重力變化也急劇增加，在容城—文安之間形成了一北東向的重力變化密集帶，北負南正，最大最小重力差值高達90×10－8m·s－2，震中在該重力變化密集帶的法線方向靠近高值區。

(2)地震期間相鄰兩年的的重力變化

地震期間相鄰兩年 (2006.05～2007.03)(圖5)的重力變化可看出，文安5.1級地震使得震源區地下的能量予以釋放，其結果在重力變化上的反映是:上期重力變化的高值區變為了重力變化低值區，震中附近的重力變化由震前的60 ×10－8m·s－2變為 －20 ×10－8m·s－2，變化高達80×10－8m·s－2，震中以南區域這種重力變化更加明顯。

(3)震后相鄰兩年的重力變化

震后相鄰兩年的重力變化 (2007.03～2008.03、2008.03～2009.03) (圖6)表現為震后恢復階段，重力負值變化區由南轉向東西兩側，繼而轉為平緩;震中區附近的重力變化則保持持續平靜。

圖4 文安MS5.1地震前1 a尺度的重力變化(a)2004.03～2005.03;(b)2005.03～2006.05Fig.4 Gravity variation in one year scale before Wen'an MS5.1 earthquake

圖5 文安MS5.1地震期間1 a尺度的重力變化Fig.5 Gravity variation in one year scale during Wen'an M5.1 earthquake

圖6 文安MS5.1地震后1 a尺度的重力變化(a)2007.03～2008.03;(b)2008.03～2009.03Fig.6 Gravity variation in one year scale after Wen'an MS5.1 earthquake

1.3 較長時段的重力變化

為了進一步分析文安及周邊地區較長趨勢的重力變化，了解重力變化的累計量，我們繪制了2005～2009年較長時期的重力變化等值線圖 (圖7)。

圖7 文安MS5.1地震后相對于2005年的重力變化Fig.7 Gravity variation after Wen'an MS5.1 earthquak relative to that in 2005

由圖7可以看出，4年尺度的重力變化總體趨勢表現為:測區西北的易縣地區有重力負值變化，測區東南部的天津、靜海、文安、任丘等有大范圍重力正值變化，并在河北任丘及天津靜海地區形成 (50～60) ×10－8m·s－2的重力變化異常區，重力變化等值線主要為北東向，與唐縣—涿州斷裂帶走向基本一致。據周海濤等 (2010)的研究，天津、靜海、文安、任丘等地區由于大量開采地下水，形成地面下降，引起地層壓密效應，造成地面重力增加。

2 結論

從文安MS5.1地震前后不同時間尺度的重力變化來看，在地震孕育階段，由于能累積累，地下物質活動劇烈，震前1 a內重力值出現正值急劇變化，在震中附近形成重力變化高梯度帶;地震后由于能量的釋放，震前震中附近的重力變化高值區變成重力變化低值區，隨著震源區附近的應力應變的釋放調整，重力變化趨于平靜。

文安MS5.1地震發生在重力變化高梯度帶的零值線附近，青藏塊體東北緣幾次地震也都發生在重力變化零值線附近，甘肅、滇西地區重力場異常變化形態與地震的關系也有這種規律，這是由于在正負重力異常變化梯度帶的零線附近，是物質密度增加與減少的過渡地帶，該處產生的物質增減差異運動劇烈，能量易于聚集，從而誘發地震 (祝意青等，1999，2001;徐云馬等，2008;梁偉鋒等，2004)。

國家地震局科技監測司.1986.文安—霸縣地殼隆起的研究［M］.北京:地震出版社.

何志堂，閆會杰，雷鵬，等.2009.地下水變化對大地原點重力影響的估算［J］.測繪標準化，27(1):39－41.

梁偉鋒，祝意青，徐云馬，等.2004.民樂MS6.1地震前的重力變化研究［J］.地震研究，27(4):334－338.

蘇樹朋，史彥華，王順昌，等.2010.文安5.1級地震前的一些前兆異常［J］.華北地震科學，26(3):10－12.

蘇樹朋，張躍剛.2009.河北省文霸網流動重力觀測對網區及周邊地震的反應［J］.西北地震學報，31(1):91－96.

王武星，石耀霖，張晶，等.2009.河北文安5.1級地震前后重力場動態變化［J］.地震，29(2):40－46.

徐云馬，祝意青，程宏賓.2008.1998－2004年滇西地區重力場演化與MS≥6.0地震［J］.大地測量與地球動力學，28(2):51－55.

張素靈，趙京軼，王建芳，等.2009.2006年文安5.1級地震的烈度異常區初探［J］.中國地震，25(1):24－30.

周海濤，郭良遷，薄萬舉，等.2010.天津地區地殼垂向形變(GPS)與文安地震［J］.大地測量與地球動力學，30(5):32－37.

祝意青，陳兵，張希，等.2001.景泰5.9級地震前后的重力變化研究［J］.中國地震，17(4):356 －363.

祝意青，胡斌，張永志.1999.永登5.8級地震前后的重力場動態圖像特征研究［J］.地殼形變與地震，19(1):71－77.

Gravity Variation in Different Time Scales before and after the Wen'an MS5.1 Earthquake

LIANG Wei-feng，ZHU Yi-qing，XU Yun-ma，LIU Lian，GUO Shu-song，LIU fang
(Second Crust Monitoring and Application Center，CEA，Xi'an 710054，Shannxi，China)

Analyzing the mobile gravity observational data in North China from 2005 and 2009 systematically and deeply，we derived dynamic variation of the regional gravity field in different time scales before and after Wen'an MS5.1 earthquake，and studied the relationship between gravity anomaly variation and seismic activity.The result showed that there existed regional gravity anomaly in a wide range in the survey area during the preparation of Wen'an MS5.1 earthquake，and there were local gravity anomaly area and high gradient zone of gravity variation referential to the earthquake preparation.The dynamic image of gravity field clearly reflected the evolution orderly and seismic activity of regional gravity field.

Wen'an M5.1 earthquake;mobile gravity observation;gravity variation;gravity field

P315.726

A

1000－0666(2012)03－0318－05

2011－10－24.

地震行業科研重點專項“多種地球物理場動態圖像與強震預測試驗”(201008007)資助.