內蒙古阿拉善左旗MS5.8震源區及周邊重力布格異常的計算與分析

朱文武　秦　昆　王同慶　鄭智江　高艷龍　張　珂

1　中國地震局第一監測中心，天津市耐火路7號，300180 2　中國地質大學土地科學技術學院，北京市學院路29號，100083

?

內蒙古阿拉善左旗MS5.8震源區及周邊重力布格異常的計算與分析

朱文武1秦昆2王同慶1鄭智江1高艷龍1張珂1

1中國地震局第一監測中心，天津市耐火路7號，300180 2中國地質大學土地科學技術學院，北京市學院路29號，100083

利用EGM2008自由空氣重力異常數據及具有1′分辨率的DEM數據，借助FA2Boug思想，對內蒙古阿拉善左旗MS5.8震源區及周邊(105°～110°E,36°～41°N)進行重力布格異常值的計算，獲取了震源區及周邊的重力布格異常信息。結果表明，震源區位于重力布格異常值變化的高梯度帶上；同時，在鮮有震例的鄂爾多斯塊體內部存在著不一致的作用力。

阿拉善左旗MS5.8地震；鄂爾多斯塊體；重力布格異常；FA2Boug；DEM

內蒙古阿拉善盟阿拉善左旗MS5.8地震發生于2015-04-15 15：39，震中位于北緯39.78°、東經106.34°。相關數據處理結果顯示，主震及主要余震基本分布在磴口-本井斷裂兩側，其中主震所在的磴口-本井斷裂西側分布最多、最集中。矩張量反演結果顯示，此次地震為略帶逆沖分量的走滑型地震，地震的矩心深度約為19 km，其機制類型與周邊的構造體系相吻合(http://www.cea-igp.ac.cn/cxdt/272066.shtml)。

Fernandez等[1]在2008年利用FA2Boug模型計算亞特蘭大-地中海過渡帶的重力布格變化異常情況，付廣裕等[2]利用該方法對2015年尼泊爾地震開展相關的研究分析工作，取得了很好的效果。本文利用最新版本(V23)的EGM2008(模型空間分辨率為5′，該模型提供高達2 160階的重力異常場球諧系數，在中國大陸上的重力異常平均精度優于10.5 mGal)自由空氣重力異常數據，再結合具有1′網格分辨率的DEM數據(V18.1)，借助FA2Boug思想進行研究區的重力布格異常計算，在此基礎上分析研究區內的重力場區域特征。

1　研究區歷史震例與地質背景

如圖1所示，巴彥烏拉山斷裂與磴口-本井斷裂之間的區域多分布正走滑型地震[3]。本次地震的發震位置位于磴口-本井斷裂附近，與該斷裂的最近距離約13 km。結合自1970年以來的歷史震例來看，以云霧山斷裂、羅山東麓斷裂、黃河-靈武斷裂、磴口-本井斷裂、鄂爾多斯北緣斷裂、包頭斷裂為界，研究區內大多數歷史地震的發震地點位于上述斷裂帶以西或以北的位置。進一步可以看出，歷史震例的發震地點更多地集中于斷層附近。而在上述斷裂另一側，即在鄂爾多斯塊體內部，歷史上幾乎沒有地震發生。

圖1　研究區歷史震例及地質構造Fig.1　History samples of earthquake and geological structure in the research region

2　重力布格異常值計算原理及FA2Boug數學模型

2.1重力布格異常計算原理

重力布格異常的表達式為[4]：

(1)

式中，g觀為實測重力值，g高為高度校正，g中為中間層校正，g0為理論重力值，g形為地形校正。

儀器在地面各點的實測重力值經高度校正、中間層校正、參考橢球面上相應的重力理論值改正以及地形校正后，便可得到重力布格異常信息。重力布格異常普遍顯示出與地形的相關性，山區的異常值往往為負，海洋地區的異常值則往往為正，而地殼在重力作用下的重力均衡及其補償運動便構成了地殼的各種復雜構造運動，導致地震不斷發生。

2.2FA2Boug數學模型

實質上，EGM2008自由空氣重力異常數據已經對g高、g中、g0進行了改正。但正如式(1)所示，還需要進行地形校正才能獲得重力布格值。

本文采用Fa2Boug算法[5-6]，結合DEM數據，對EGM2008自由空氣重力異常數據進行改正，進而獲得重力布格異常值。該算法以待計算的布格重力異常點位為中心，根據距離不同將其周圍的自由空氣重力異常數據/高程數據劃分為3個區域，即遠區、中區、近區。以待算點為中心，3個區域的劃分方式如下：當距待算點的距離(d)為20

1)遠區：

(2)

2)中區：

(3)

3)近區:

(4)

(5)

(6)

根據遠、中、近區的重力布格異常計算結果，結合式(1)，得出待算點的最終重力異常布格值：

(7)

3　重力布格異常計算結果及分析

根據此次地震的發震地點，本文選擇105°～110°E、36°～41°N范圍作為研究區，在研究區內改正后重力布格異常值約在-245～-76 mGal之間，重力異常高值區主要分布于研究區東經108°以東的全部范圍內，以及震中西北和東南兩側附近地區；而重力異常低值區則位于研究區西南側以及震中以北的范圍內。在幾乎沒有歷史震例的鄂爾多斯塊體內部，其西側重力布格異常值明顯低于東側的重力布格異常值。這說明，穩定的鄂爾多斯塊體內部可能存在著一定程度的擠壓和拉張作用，原因似乎是由于塊體內部西側受到拉張的作用致使殼內物質密度下降，而東側則受到擠壓導致密度提升進而導致重力布格異常值升高。

為進一步研究重力布格異常信息與此次地震的關系，如圖2中綠色線條所示，沿著與磴口-本井斷裂垂直方向繪制一條剖面(圖3)，該剖面恰好經過本次地震的發震地點，其中紅色五星代表此次地震的發震位置。結合圖2、圖3可以看出，震中位置與研究區磴口-本井斷裂距離最近，而重力布格異常的變化在該斷裂帶的兩側呈現出高梯度的變化特征。這在一定程度上說明，重力布格異常特征可以反映與斷裂構造相關的地殼內剩余密度異常體的分布情況。

圖2　研究區重力布格異常計算結果Fig.2　Calculation of Bouguer anomaly in research region

圖3　重力異常計算結果剖面圖Fig.3　Profile of Bouguer anomly results

4　結　語

本次地震發震地點位于重力布格異常值變化的高梯度帶上，因此建議在研究區內應重點關注斷裂帶附近和斷裂帶以西、以北重力布格異常梯度變化較大的區域。結合歷史震例可以看出，研究區內地震的發生地點多位于斷裂帶上，或位于斷裂帶以西、以北的位置。相反，在鄂爾多斯塊體內部，即在研究區內西至銀川、東至榆林、南至桐川、北至包頭，東西長近250 km、南北長近600 km

的范圍內鮮有地震發生，且最大震級也僅為4.2級。然而從重力布格異常的計算結果可以看出，即使在相對穩定的鄂爾多斯塊體內部，重力布格異常信息仍呈現相反的兩種態勢，這說明塊體內部也存在不一致的作用力。

[1]Fernandez M, Fullea J, Zeyen H. FA2Boug-A FORTRAN 90 Code to Compute Bouguer Gravity Anomalies from Gridded Free Air Anomalies: Application to the Atlantic-Mediterranean Transition Zone[J]. Computers & Geosciences, 2008, 34(12):1 665-1 681

[2]付廣裕, 高尚華, 張國慶, 等. 2015年尼泊爾MS8.1地震的地殼重力均衡背景與地表形變響應特征[J]. 地球物理學報,2015,58(6):1 900-1 908(Fu Guangyu, Gao Shanghua, Zhang Guoqing, et al. Gravitational Isostasy Background and Surface Deformation Response Characteristics of the 2015 Nepal MS8.1 Earthquake[J]. Chinese Journal of Geophysics, 2015,58(6):1 900-1 908)

[3]韓曉明，劉芳，胡博，等. 河套地震帶的震源機制類型時空分布特征[J].吉林大學學報：地球科學版,2015(2):592-601(Han Xiaoming, Liu Fang, Hu Bo,et al. Space-Time Distribution Characteristics of the Focal Mechanism Type in Hetao Seismic Belt[J].Journal of Jilin Unviersity：Earth Science Edition,2015(2):592-601)

[4]周友華. 布格異常及均衡異常物理意義的探討[J].西北地震學報,1995,17(2):84-86(Zhou Youhua. Discussion of Significance of Bouguer Anomly and Residual Anomly[J]. Northwestern Seismological Journal, 1995,17(2):84-86)

[5]Kane M F. A Comprehensive System of Terrain Corrections Using a Digital Computer[J]. Geophysics,1962, 27 (4):445-462

[6]Lopez H R B.FORTRAN Program for Automatic Terrain Correction on Gravity Measurements[J]. Computers & Geosciences, 1990, 16 (2):237-244

About the first author:ZHU Wenwu,postgraduate, majors in process and analysis of deformation data, E-mail:fmccea@163.com.

Calculation and Analysis of Bouguer Anomaly in Alxa Left Banner MS5.8 Epicenter and Sourrounding Regions of Inner-Mongolia

ZHUWenwu1QINKun2WANGTongqing1ZHENGZhijiang1GAOYanlong1ZHANGKe1

1First Crust Monitoring and Application Center, CEA, 7 Naihuo Road,Tianjin 300180, China 2School of Land Science and Technology, China University of Geosciences, 29 Xueyuan Road,Beijing 100083, China

The gravity Bouguer anomaly values in Alxa Left Banner MS5.8 epicenter and surrounding regions of Inner-Mongolia (105°-110°E,36°-41°N) are obtained using the FA2Boug algorithm with EGM2008 free-air gravity anomaly and DEM data. The results show that the epicenter locates at the high gradient brand of the gravity Bouguer anomaly, and there exists different forces inside the Erdos block.

Alxa Left Banner MS5.8 earthquake; Erdos block; gravity Bouguer anomaly; FA2Boug; DEM

Combination Project with Monitoring, Prediction and Scientific Research of Earthquake Technology, CEA, No.163304; Science and Technology Basic Work Program, Ministry of Science and Technology of the People’s Republic of China,No.2015FY210400.

2015-09-28

朱文武，碩士生，研究方向為重力、水準等形變數據處理與分析，E-mail:fmccea@163.com。

10.14075/j.jgg.2016.10.002

1671-5942(2016)010-0851-03

P223

A

項目來源：中國地震局三結合課題(163304)；國家科技部科技基礎性工作專項(2015FY210400)。