玄松柏 申重陽 邢樂林 李 輝

1)中國地震局地震研究所(地震大地測(cè)量重點(diǎn)實(shí)驗(yàn)室)，武漢 430071

2)中國地震局地殼應(yīng)力研究所武漢科技創(chuàng)新基地，武漢430071

1 引言

太原盆地北起石嶺關(guān)，南至韓侯嶺，東西以邊山斷裂為界，分別受太谷斷裂和交城斷裂控制，呈NESW 走向，為典型的新生代斷陷盆地。地質(zhì)學(xué)和地震學(xué)證據(jù)表明晚第四紀(jì)以來以0.5 ～1.6 mm/a 的速率擴(kuò)展［1，2］和約3 mm/a 的地表沉降［3］。大量現(xiàn)代水準(zhǔn)測(cè)量結(jié)果驗(yàn)證了太原盆地繼承了第四紀(jì)以來的構(gòu)造活動(dòng):胡惠民等［4］的結(jié)果顯示太原盆地是山西斷陷帶沉降最為顯著的斷陷盆地，在10 mm/a 以上;王秀文等［5］給出了太原盆地地表垂直形變形態(tài)，認(rèn)為隨著城市用水對(duì)地下水的開采需求逐漸增大，地表沉降日趨嚴(yán)重;張崇立［3］認(rèn)為斷陷盆地的沉降是上地幔的局部隆起作用的結(jié)果;韓月萍等［6］通過華北地區(qū)水準(zhǔn)結(jié)果獲得山西斷陷帶盆地沉降呈現(xiàn)北端最小，南端最大的空間漸進(jìn)形態(tài);張四新等［7］的水準(zhǔn)剖面結(jié)果顯示垂直形變幅度基本在20 mm 以內(nèi)。

目前，關(guān)于地表垂直形變的研究以水準(zhǔn)和GPS資料和成果居多，而利用重力資料獲得地表形變的研究則較少［8－10］。本文利用山西地震監(jiān)測(cè)網(wǎng)重力復(fù)測(cè)資料和太原基準(zhǔn)地震臺(tái)(以下簡(jiǎn)稱晉祠)絕對(duì)重力資料，結(jié)合地下水水位變化成果，獲得了太原盆地地表垂直形變趨勢(shì)，以檢測(cè)太原盆地的地表沉降。

2 資料與方法

2.1 重力資料處理

太原盆地的重力資料包括相對(duì)重力聯(lián)測(cè)資料和絕對(duì)重力測(cè)量資料。2009年以前，僅實(shí)施每年3 ～4期的相對(duì)重力測(cè)量，而未進(jìn)行絕對(duì)重力觀測(cè)，因此，我們采用盆地內(nèi)每年實(shí)施4 期的重力點(diǎn)平均值的方式來獲取太原盆地整體重力變化［8］。2009年以后，華北地區(qū)重力網(wǎng)改造，相對(duì)重力聯(lián)測(cè)則改為每年2期，部分點(diǎn)位亦有變更，因而，利用晉祠點(diǎn)的絕對(duì)重力資料分析其重力變化，一定程度上可反映太原盆地的地表沉降。

2.1.1 相對(duì)重力聯(lián)測(cè)資料處理

2009年以前，山西地震監(jiān)測(cè)重力網(wǎng)每年施測(cè)3～4 期，其中太原盆地內(nèi)包括13 個(gè)重力點(diǎn)(圖1)，重力復(fù)測(cè)資料處理采用擬穩(wěn)平差的方式，通過資料較完備的4 期平均的方式以減弱季節(jié)性變化的影響。由于擬穩(wěn)點(diǎn)(太原點(diǎn))在平差結(jié)果反映的重力變化較?。?1］，而地下水位變化則較為明顯［12，13］，因此，重力復(fù)測(cè)資料中不包括地下水位變化引起的重力變化?；诖?，為了獲取太原盆地整體的重力變化趨勢(shì)，首先將太原盆地內(nèi)2001—2006年資料完備的10 個(gè)重力點(diǎn)分別取年均值，作為每個(gè)點(diǎn)的當(dāng)年重力值，然后根據(jù)Ekman et al［8］的思路，將10 個(gè)重力點(diǎn)的年均值再進(jìn)行平均的方法，作為盆地整體重力值(表1)，進(jìn)一步采用最小二乘線性擬合得到太原盆地的重力年變化率。

2.1.2 絕對(duì)重力資料處理

2009年以來，中國地震局地震研究所在晉祠開始實(shí)施每年兩期的絕對(duì)重力測(cè)量和相應(yīng)的垂直梯度測(cè)量，其精度一般為(1 ～2)×10－8ms－2，至2012年已積累了4年7 期較完備的絕對(duì)重力資料，采用最小二乘線性擬合得到晉祠點(diǎn)的絕對(duì)重力年變化率。與重力復(fù)測(cè)資料不同的是，絕對(duì)重力資料中包含地下水位變化的重力效應(yīng)。

2.2 地表沉降模型

太原盆地周邊地塹的隆升量級(jí)較小，對(duì)于盆地內(nèi)重力測(cè)點(diǎn)來講，其與斷陷盆地沉降的重力效應(yīng)相比可以忽略，若不考慮地下水水位變化的影響，設(shè)地表觀測(cè)重力變化為Δgobs，地殼平均密度ρc，則地表沉降ht(以向下為正)引起的重力變化Δg表示為:

太原盆地內(nèi)因城市地下水開采，地下水水位變化較為明顯，模型中須考慮地下水水位變化的重力效應(yīng)ΔgWT，則式(1)變?yōu)?

地下水水位變化的重力效應(yīng)ΔgWT為［14］

式中，φe為巖石的孔隙度，δz 為地下水水位變化。

綜合式(1)～(3)，取地殼平均密度ρc=2.67 ×103kg/m3，得到地表沉降率為:

3 太原盆地重力變化與地表沉降

3.1 太原盆地重力變化

獲得的太原盆地的重力年變化率為(4.3 ±0.5)×10－8ms－2/a(圖2(a))，因該重力變化基本不包含地下水位變化的影響，反映了太原盆地重力增加的趨勢(shì)。

圖2 太原盆地及晉祠重力變化時(shí)序Fig.2 Sequence of gravity change in Taiyuan Basin and Jinci Area

2009—2012年近4年來，晉祠點(diǎn)絕對(duì)重力變化曲線如圖2(b)，太原臺(tái)重力變化呈現(xiàn)增大的趨勢(shì)，其速率為(9.7 ±0.5)×10－8ms－2/a。雖然水準(zhǔn)測(cè)量結(jié)果揭示太原市內(nèi)以吳家堡為代表的沉降中心沉降量較大［5］，而其與晉祠點(diǎn)相距約20 km，根據(jù)已有的模擬計(jì)算結(jié)果［15］，其對(duì)太原臺(tái)重力效應(yīng)影響可以忽略，可以認(rèn)為絕對(duì)重力變化主要由地下水水位變化和點(diǎn)位附近的垂直形變引起。

表1 太原盆地內(nèi)測(cè)點(diǎn)重力變化Tab.1 Gravity changes of stations in Taiyuan Basin

3.2 太原盆地地表沉降

相對(duì)重力聯(lián)測(cè)與絕對(duì)重力測(cè)量結(jié)果均顯示太原盆地重力變化呈現(xiàn)增大的趨勢(shì)，若不考慮地下水水位變化的重力效應(yīng)，將重力變化全部歸結(jié)為地表沉降的影響，相對(duì)重力聯(lián)測(cè)資料顯示太原盆地整體以21.8 mm/a 的速率沉降(2001—2006年)，與精密水準(zhǔn)測(cè)量結(jié)果［6］相近。是因?yàn)橄鄬?duì)重力聯(lián)測(cè)資料采用擬穩(wěn)平差的方式，其中未包括地下水水位變化的影響，考慮到2001—2006年太原臺(tái)地下水水位以速率約為2 m/a［12，13］下降的重力效應(yīng)，則2001—2006年太原盆地的整體重力變化率約為8.5 × 10－8ms－2/a，與2009—2012年期間絕對(duì)重力年變化率基本一致，說明該時(shí)段太原盆地地表沉降是地下水開采和繼承性沉降運(yùn)動(dòng)的綜合效應(yīng)，地下水開采的重力效應(yīng)與繼承性沉降的重力效應(yīng)相當(dāng)。

用絕對(duì)重力觀測(cè)資料進(jìn)行地下水水位變化改正，以獲取因垂直形變的重力效應(yīng)。晉祠鄰近區(qū)域含水層為石灰?guī)r、白云巖等碳酸鹽類，其孔隙度一般約為5%。2009年以來，太原臺(tái)地下水水位約以4 m/a 的速率持續(xù)上升［13］，其重力效應(yīng)約為8.4 ×10－8ms－2/a，則2009—2012年由地表沉降引起的重力變化為1.3 ×10－5ms－2，表現(xiàn)為6.6 mm/a 的速率沉降。

4 結(jié)論與討論

1)通過太原盆地內(nèi)10 個(gè)點(diǎn)的重力復(fù)測(cè)資料得到2001—2006年太原盆地重力變化為(4.3 ±0.5)×10－8ms－2/a，地表沉降速率為21.8 mm/a;

2)太原臺(tái)2009年以來絕對(duì)重力年變化率為(9.7 ±0.5)×10－8ms－2/a，進(jìn)一步結(jié)合地下水水位變化成果得到太原盆地地表沉降速率為6.6 mm/a。

2001—2006年地表沉降遠(yuǎn)遠(yuǎn)大于2009年以后的結(jié)果，究其原因，可能是地下水等的大量開采造成地表大幅沉降。2009年以后，隨著汾河清水復(fù)流工程的開展和大量飲用黃河水，地下水水位升高并逐漸恢復(fù)，則現(xiàn)今垂直運(yùn)動(dòng)與地質(zhì)和近期水準(zhǔn)測(cè)量結(jié)果［3，7］基本一致，主要表現(xiàn)為繼承性沉降。

需要指出的是，因山西斷陷帶是印度板塊與歐亞板塊碰撞和太平洋板塊俯沖作用的過渡部位，地震活動(dòng)頻繁，且文中未對(duì)2008年5月12日汶川地震和2011年3月11日日本地震的影響進(jìn)行分析，若利用重力資料精確獲取地表垂直形變特征，則以上兩次強(qiáng)震和山西斷陷帶地震活動(dòng)對(duì)垂直形變產(chǎn)生的間接作用的分析研究十分必要。

1 Zhang Y Q，Mercier J L and Vergély P.Extension in the graben systems around the Ordos(China)，and its contribution to extrusion tectonics of south China with respect to Gobi-Mongolia［J］.Tectonophysics，1998，285(1－2):41－75.

2 Xu X and Ma X.Geodynamics of Shanxi rift system，China［J］.Tectonophysics，1992，208(1－3):325－340.

3 張崇立.山西斷陷帶垂直形變特征及其成因初探［J］.中國地震，1993，9(4):327－333.(Zhang C L.Character-istics and origin of vertical deformation in Shanxi rift system，Shanxi province［J］.Earthquake Research in China，1993，9(4):327－333)

4 胡惠民.山西及其鄰近地區(qū)地殼垂直形變與地震危險(xiǎn)性分析［J］.山西地震，1993，73(2):23－28.(Hu H M.Crustal vertical deformation and analysis for earthquake risk in Shanxi province and its vicinity［J］.Earthquake Research in Shanxi，1993，73(2):23－28)

5 王秀文，郭躍宏，范雪芳.山西地區(qū)近期垂直形變場(chǎng)分析［J］.地殼形變與地震，2001，(2):64－69.(Wang X W，Guo Y H and Fan X F.Analysis of resent vertical deformation in Shanxi area［J］.Crustal Deformation and Earthquake，2001，(2):64－69)

6 韓月萍，等.華北北部地區(qū)現(xiàn)今地殼垂直形變特征與地震危險(xiǎn)性分析［J］.大地測(cè)量與地球動(dòng)力學(xué)，2010，(2):25－28.(Han Yueping，et al.Characteristics of present-day crust vertical deformation and analysis for earthquake risk in northern area of North China［J］.Journal of Geodesy and Geodynamics，2010，(2):25－28)

7 張四新，張希.汾渭斷陷帶現(xiàn)今垂直形變與近期地震活動(dòng)性［J］.地震，2012，32(4):123－130.(Zhang S X and Zhang X.Current vertical deformation in the Fenwei faultdepression zone in Shanxi province and future seismicity［J］.Earthquake，2012，32(4):123－130)

8 Ekman M，et al.Gravity change and land uplift in Fennoscandia1966－1984［J］.Bulletin Géodésique，1987，61(1):60－64.

9 Sun W K，et al.Gravity and GPS measurement reveal mass loss beneath the Tibetan Plateau:geodetic evidence of increasing crustal thickness［J］.Geophys Res Lett.，2009，36(2):1－5.

10 邢樂林，等.用拉薩點(diǎn)大地測(cè)量資料檢測(cè)青藏高原地殼的增厚［J］.測(cè)繪學(xué)報(bào)，2011，30(1):41－44.(Xing L L，et al.Present-day crustal thickness increasing in Qinghai-Tibetan plateau by using geodetic data at Lasa station［J］.Acta Geodaetica Cartographica Sinica，2011，30(1):41－44)

11 玄松柏，等.山西斷陷盆地帶及其鄰區(qū)1999—2008年地殼物質(zhì)密度變化［J］.大地測(cè)量與地球動(dòng)力學(xué)，2013，(5):7－10.(Xuan Songbei，et al.Density change of crustal material in Shanxi Rift and its adjacent area during 1999 to 2008［J］.Journal of Geodesy and Geodynamics，2013，(5):7－10)

12 孟彩菊.太原地震臺(tái)短水準(zhǔn)受水位影響分析及其異常識(shí)別［J］.地震地磁觀測(cè)與研究，2010，31(1):77－83.(Meng C J.Effect of water level on short leveling in Taiyuan Seismic Station and identification of anomaly［J］.Seismological and Geomagnetic Observation and Research，2010，31(1):77－83)

13 沈曉松，等.原地震臺(tái)水位異常分析［J］.地震地磁觀測(cè)與研究，2012，33(5/6):178－182.(Shen X S，et al.Analysis of water level anomaly in Taiyuan Seismic Station［J］.Seismological and Geomagnetic Observation and Research，2012，33(5/6):178－182)

14 M Battaglia，Segall P and Roberts C.The mechanics of unrest at Long Valley caldera，California.2.Constraining the nature of the source using geodetic and micro-gravity data［J］.J Volcanol Geotherm Res.，2003，127:219－245.DOI:10.1016/S0377－0273(03)00171－9.

15 劉少明，高榮勝，賈民育.香山絕對(duì)點(diǎn)周圍的局部地殼變形及其對(duì)重力測(cè)量結(jié)果的影響［J］.大地測(cè)量與地球動(dòng)力學(xué)，2002，(2):97－100.(Liu Shaoming，Gao Rongsheng and Jia Mingyu.The local crustal deformation around Xiangshan Absolute Station and its effect on the results of gravity surveying［J］.Journal of Geodesy and Geodynamics，2002，(2):97－100)