巴東地區均衡重力異常與構造活動特征*

楊光亮 申重陽 談洪波 吳桂桔 王嘉沛

1)中國科學院計算地球動力學實驗室，北京 100049

2)湖北省地震局，武漢 430071

3)地殼運動與地球觀測實驗室，武漢 430071

巴東地區均衡重力異常與構造活動特征*

楊光亮1，2，3)申重陽2，3)談洪波2，3)吳桂桔2，3)王嘉沛2，3)

1)中國科學院計算地球動力學實驗室，北京 100049

2)湖北省地震局，武漢 430071

3)地殼運動與地球觀測實驗室，武漢 430071

基于巴東地區布格重力異常信息，采用艾黎-海斯卡寧均衡模型，計算了該區均衡重力異常。結果顯示，該區均衡重力異常變化范圍為－45×10－5ms－2～45×10－5ms－2，大部分區域基本處于重力均衡狀態，局部存在弱正異常與弱負異常，弱正均衡異常呈條帶狀。研究表明，太行-武陵重力梯級帶基本處于均衡狀態，內部存在差異;巴東地區西北部存在一定程度的地殼均衡調整空間;巴東Ms5.1地震區域處于重力均衡狀態，推測其發震與地殼均衡調整無直接關系。

巴東地區;地震;重力;布格異常;均衡異常

巴東地區處于揚子克拉通北緣與四川盆地東端，北部與秦嶺造山帶相連，東北及東南部分別為華北克拉通及華南塊體，區內主要為一系列北東向弧形褶皺，褶皺軸自西向東漸轉為北東東，后以近東西向嵌入秭歸向斜中［1］。在三峽水庫蓄水前，巴東及鄰區地震活動具有低頻度、高強度的特點，蓄水后地震事件以微小震群特征出現［2］。三峽蓄水前后的精密重力測量和模擬研究均表明［3－5］，庫壩附近水荷載的最大重力變化約 200×10－8ms－2，影響范圍離庫岸約5 km隨后迅速衰減。由此可見，三峽水庫蓄水對該區地震活動存在一定影響，但其影響強度及范圍有限。因此，地震活動特征發生改變與其地殼構造背景密切相關，厘清水庫地震深部成因機制是水庫地震監測的重要方向。

近幾十年來，很多學者在該區開展了大量的地球物理觀測與研究［6－9］，獲得了該區地質、速度結構及密度結構模型。本文利用收集到的布格重力異常信息對該區均衡重力異常及其特征進行分析，了解該區地殼均衡特點及構造背景信息。

1 重力與地形數據

本文研究范圍為巴東地區3°×3°的區域(圖1)，圖1中紅色圓點為2013-12-16的Ms5.1地震震中。從圖1可見，區域重力場整體上為負異常，異常值自北西向南東方向逐漸遞增，存在約120×10－5ms－2的重力梯度變化。該區中部為橫亙中國東部的太行-武陵重力梯級帶的一部分，呈南東-北西向弧形展布。該梯級帶將研究區分隔為東南部的江漢幔隆區和西北部的鄂西幔陷區，該幔隆區和幔陷區布格異常分別集中在 －25×10－5ms－2和 －100×10－5ms－2附近(圖1中布格異常統計分布的兩個峰值)。在遠離重力梯級帶的區域，布格異常變化趨于平緩。該梯級帶的出現反映深部存在規模宏大的構造變異帶，顯示本區在深部構造變異帶的控制下，江漢幔隆向鄂西幔陷過渡的構造特征。該區在歷史上經歷了3次較大的構造運動，導致區內西北部重力異常局部劇烈變化，并導致重力梯級帶扭曲，反映了殼層中淺源構造信息，如基底、侵入巖體、巖性構造、斷陷盆地等導致的物性差異及密度不均勻。

該地區的ASTER GDEM數字高程模型顯示(圖2)，其垂直精度達到7 m，水平精度約30 m。該區地形高差接近3 000 m，西北部地形起伏較大，地勢較高，東南部地形平坦，最低處海拔高度接近0。

2 均衡重力異常

采用艾黎-海斯卡寧均衡模型計算均衡重力異常。按離計算點的距離將地形劃分為近、中、遠3個方形區域，其距離分別為10、30和166.7 km。3個區域分別采用不同的公式［10］，計算獲得該區均衡重力異常分布如圖3。

該區均衡重力異常變化范圍為(－45～45)×10－5ms－2，呈弱正異常及弱負異常分布，主要集中在±20 ×10－5ms－2以內(圖 3(b))。均衡理論認為，均衡重力異常為0時，地殼處于完全均衡狀態，均衡重力異常為正值表示實際的地殼厚度比均衡時的地殼厚度小，反之則厚度大。通常認為，±20×10－5ms－2時基本處于均衡狀態。據此，本研究區大部分區域基本接近均衡狀態。該區局部均衡正異常以條帶狀出現，主要集中在三峽水庫兩側的神農架林區及五峰、長陽和石門一帶，且主要分布在太行-武陵重力梯級帶的西部，弱正均衡異常包圍了三峽庫首區;三峽庫區、江漢平原、竹山及澧縣的江漢幔隆區則顯示弱負異常。

圖1 巴東地區布格重力異常(a)及統計分布(b)Fig.1 Bouguer gravity anomaly in Badong region(a)and statistical chart(b)

圖2 巴東地區地形圖(ASTER GDEM 2009)Fig.2 Topography of Badong region(ASTER GDEM 2009)

3 討論

為了進一步分析該區重力異常特征，對該區均衡重力異常及布格重力異常圖，分別作沿該區兩個對角線及兩條正交穿過2013-12-16巴東Ms5.1地震震中的縱切剖面。其中剖面A-A1、B-B1分別作布格重力(圖5(a))和均衡重力剖面(圖5(b))。剖面A-A1是太行-武陵重力梯級帶的一部分，顯示的均衡異常變化范圍為(－25～20)×10－5ms－2，存在一定正負異常的起伏變化，但基本處于重力均衡狀態;剖面均衡異常呈現中間高、兩端低，呈弧形。因此，太行-武陵重力梯級帶基本處于均衡狀態，局部存在差異。其布格重力異常變化幅度為－(100～30)×10－5ms－2，沿剖面依次遞增，可見太行-武陵重力梯級帶沿西南向逐漸消減。

剖面B-B1與太行-武陵重力梯級帶基本正交，其均衡重力異常在東北端呈現較為劇烈的變化，在160 km范圍內依次出現弱負異常(－50×10－5ms－2)和弱正異常(40×10－5ms－2)，顯示該區存在一定的地殼結構調整的可能;其布格重力異常變化幅度為120×10－5ms－2左右，在房縣斷裂附近也存在一定幅度的劇烈變化，顯示斷裂帶兩側物質的不均勻性。

2013 -12 -16 巴東Ms5.1地震基本處于地殼均衡至弱負均衡間，附近也無已知斷裂，地質結構較為穩定。圖5中C-C1剖面及D-D1剖面過該地震震中，分別在0.4°范圍內以2 km間距作150 km正交采樣，圖中疊加顯示了各剖面的均衡重力變化及其平均值。震中南北向剖面穿過了均衡異常正負過渡區，平均變化范圍為(－20 ～40)×10－5ms－2，相鄰各切面變化幾乎一致(圖5(c))。震中東西向剖面與三峽水庫走向一致，平均變化范圍為±20×10－5ms－2，但相鄰切面存在差異，顯示局部物質分布不均勻。總體上，該區均衡重力變化幅度不大，基本處于均衡狀態。由此可見，巴東Ms5.1地震無均衡調整導致構造運動的直接證據。

4 結論

圖3 巴東地區均衡重力異常及統計分布Fig.3 Isostatci gravity anomaly in Badong region and statistical chart

圖4 巴東地區重力異常縱切剖面位置Fig.4 Gravity anomaly profiles in Badong region

圖5 巴東地區重力異常縱切剖面Fig.5 gravity anomalies of profile in Badong region

對巴東地區布格重力進行均衡改正，獲得均衡重力異常。均衡改正相當于將有密度差的山根面引起的重力異常從實測布格重力異常中消除。該區均衡重力異常變化范圍為(－45～45)×10－5ms－2，大部分區域基本處于均衡狀態，局部存在弱正異常與弱負異常，局部區域存在地殼結構調整空間，弱正均衡異常呈現條帶狀。研究表明:

1)太行-武陵重力梯級帶基本處于均衡狀態，局部存在差異。布格重力異常變化幅度為－(100～30)×10－5ms－2，沿剖面依次遞增，太行 － 武陵重力梯級帶存在一定的內部差異。

2)巴東地區西北部的房縣一段存在弱正弱負均衡異常過渡帶，存在一定的地殼均衡調整空間。

3)巴東Ms5.1地震區域處于重力均衡狀態，推測其發震與地殼均衡調整無直接關系。這一結論與文獻［11］的研究結果存在差異，其原因可能是本文使用的模型比較簡單所致，實際的巖石圈模型更為復雜，殼幔的橫向、垂向密度均有變化;靜力均衡作用是建立在低密度殼上伏于高密度流體底層上的，地殼下的上地幔很難認為具有流體性或可塑性。

1 閆舉生，譚建民，李明.鄂西巴東縣地質災害發育特征及分布規律初探［J］.華南地質與礦產，2012，28(1):59－63.(Yan Jusheng，Tan Jianmin，Li Ming.Occurrence and distribution characteristics of geological hazards in Badong County，western Hubei Province［J］.Geology and Mineral Resources of South China，2012，28(1):59 －63)

2 李峰，等.三峽水庫巴東地區的地震活動分析［J］.大地測量與地球動力學，2008(4):63 －67.(Li F，et al.Analysis of seismicity in Badong area of Three Gorges reservoir［J］.J Geod Geodyn，2008(4):63 －67)

3 楊光亮，等.三峽水庫蓄水重力效應數值模擬［J］.大地測量與地球動力學，2005(1):19 －23.(Yang Guangliang，et al.Numerical simulation of gravity effect of water impoundment in Three Gorges reservoir［J］.J Geod Geodyn，2005(1):19－23)

4 孫少安，等.三峽工程蓄水前后的精密重力測量［J］.大地測量與地球動力學，2004(2):30 －33.(Sun Shaoan，et al.High precise gravity survey before and after impoundment of Three Gorges project［J］.J Geod Geodyn，2004(2):30 －33)

5 申重陽，等.長江三峽庫首區近期重力場動態變化［J］.大地測量與地球動力學，2004(2):6－13.(Shen Chongyang，et al.Dynamic variations of gravity field in head area of Three Gorges reservoir in recent years［J］.J Geod Geodyn，2004(2):6－13)

6 甘家思，劉鎖旺.黃陵地塊內部北西向霧渡河斷裂的再研究［J］.地殼形變與地震，1996(1):72 －78.(Gan Jiasi，Liu Suowang.In vestigation advance of the Wuduhe fault in the Huangling block［J］.Crustal Deformation and Earthquake，1996(1):72 －78)

7 李安然，等.鄂西三峽區域重力場特征及其地震構造意義［J］.地震地質，1987，9(3):71 － 78.(Li Anran，et al.Characteristics of gravity field in the Three-Gorge region，western Hubei Province and its seismotectonic significance［J］.Seismology and Geology，1987，9(3):71 －78)

8 廖武林，等，三峽地區P波速度層析成像研究［J］.大地測量與地球動力學，2007(3):80 －84.(Liao Wulin，et al.Tomographic imagery of P wave velocity structure in Three Gorges region［J］.J Geod Geodyn，2007(3):80 －84)

9 張毅，等.三峽地區中上地殼密度結構［J］.地球科學，2012，37:213 － 222.(Zhang Yi，et al.Density structure of upper and middle crust in Three Goreges reservoir area［J］.Earth Science，2012，37:213 －222)

10 程振炎，等.均衡重力異常計算方法［J］.物探與化探，1985，9(5):340 － 346.(Cheng Zhenyan，et al.Computing method for isostatic gravity anomalies［J］.Geophysical and Geochemical Exploration，1985，9(5):340 －346)

11 Wei Jin.Analytical study on abnormal change intime_variable gravity at Yichang seismostation before the M5.1 Badong earthquake ［J］.Geodesy and Ceodynamics，2014(1):55－63.

ISOSTATIC GRAVITY ANOMALIES OF BADONG REGION AND TECTONIC ACTIVITY

Yang Guangliang1，2，3)，Shen Chongyang2，3)，Tan Hongbo2，3)，Wu Guiju2，3)and Wang Jiapei2，3)

1)Key Laboratory of Computional Geodynamics，CAS，Beijing 100049
2)Institute of Seismology，CEA，Wuhan 430071
3)Crustal Movement Laboratory，Wuhan430071

Isostatic gravity anomaly in Badong area was calculated using the Airy model.The results show that the isostatic gravity anomaly ranged from －45 ×10－5ms－2to 45 ×10－5ms－2，the most of the region was in the gravity isostatic state.Weak positive anomaly and weak negative anomaly appeaed in local areas，and weak positive isostatic anomaly presents in strip.The study shows that Taihang，Wuling gravity zone basically was in isostatic state with some differences in the internal.The northwestern in Badong area might occur isostatic adjustment.Badong Earthquake region was in the gravity isostatic state，that means the Badong Earthquake does not relate to the isostatic adjustment.

Badong area;earthquake;gravity;Bouguer anomaly;isostatic anomaly

P315.7

A

1671-5942(2014)03-0035-04

2014-03-26

國家自然科學基金青年基金項目(41204014);中國地震局地震預測研究所基本科研業務費專項(2011IES010103，2013IES010103，IS200916004)。

楊光亮，男，1980年生，博士，主要從事重力場變化與地球動力學研究。E－mail:ygl@eqhb.gov.cn。