九江地震臺洞體應(yīng)變觀測環(huán)境影響分析

趙愛平 周紅艷 羅建明

1）中國南昌330039江西省地震局

2）中國江西342600會昌地震臺

0 引言

依據(jù)連續(xù)地傾斜、地應(yīng)變等形變觀測結(jié)果，地形變可表示為周期應(yīng)變波、構(gòu)造運動、地震前斷層加速擴展和觀測環(huán)境（儀器）等因素引起的聯(lián)合變形（牛安福，2007），其中，觀測環(huán)境變化導(dǎo)致的地形變異常現(xiàn)象是普遍存在的，弄清環(huán)境因素變化如何影響是應(yīng)用地形變資料進行地震前兆異常分析的前提，很多研究者（劉序儼等，1991；陳敏等，2007；王曰鳳，2007；高昂，2013）提出了定點形變觀測存在的干擾因素、識別及排除干擾的方法。目前，采用SS-Y型伸縮儀進行洞體應(yīng)變觀測是常見方法，洞體應(yīng)變觀測所獲得的線應(yīng)變固體潮觀測資料可進行潮汐主應(yīng)變、主方向和剪應(yīng)變的計算（劉序儼等，1988；蔣駿等，1994），從而獲知觀測場地的平面應(yīng)變狀態(tài)，進而分析區(qū)域應(yīng)力應(yīng)變場的變化和地震前兆異常。由于SS-Y型伸縮儀具有較高的觀測靈敏度，觀測資料常記錄到多種觀測環(huán)境變化帶來的影響，大多被當(dāng)成“干擾響應(yīng)”進行剔除，恰反映了該儀器對場地應(yīng)力應(yīng)變場變化的記錄能力。

九江地震臺（以下簡稱九江臺）的洞體應(yīng)變觀測亦受到諸多觀測環(huán)境變化的影響，如集中降雨、氣溫、氣壓、抽水等，隨著數(shù)字化洞體應(yīng)變資料的積累，深入分析資料的需求越來越迫切，從機理上弄清楚環(huán)境因素對洞體應(yīng)變的影響就顯得尤其重要。本文選用九江臺2008—2013年洞體應(yīng)變及觀測環(huán)境變化資料，分析洞體應(yīng)變異常形態(tài)、數(shù)值關(guān)系及其關(guān)聯(lián)性，并進而討論應(yīng)變異常的形成機制。

1 臺站觀測環(huán)境

九江地震臺位于江西省北部，地處廬山西北側(cè)，在大地構(gòu)造上屬于揚子陸塊下?lián)P子地塊中部，區(qū)域構(gòu)造位于廬山山體西北緣邊緣剪切帶的夏家—威家?guī)X左行走滑斷裂帶，該斷裂帶走向NE，延伸約8 km，寬約1—1.5 km，總體傾向320°—330°，傾角60°—70°，在中—新生代以來長期多次活動，繼承了早第四紀(jì)時發(fā)生的斷塊抬升運動，并控制第四紀(jì)的分布（江西地質(zhì)礦產(chǎn)勘查開發(fā)局，1993）。1911年2月6日該斷裂帶發(fā)生5.0級地震，距臺站僅1.5 km。

九江臺形變觀測山洞進深45 m，覆蓋層厚約25 m，洞內(nèi)年溫差<1.0℃，日溫差<0.02℃，洞體基巖為硅質(zhì)灰?guī)r。觀測山體相對高度約100 m，近EW走向，山體呈EW向緩、NS向陡的形狀，表層植被豐富，山體多見巖層裂隙。形變觀測項附近分布4口井，見圖1和表1。九江臺伸縮儀采用含鈮特種因瓦材料做基線材料，其膨脹系數(shù)低于0.33×10-6/℃，山洞NS向基線長9.04 m，EW向基線長11.92 m，2007年8月開始洞體應(yīng)變數(shù)字化觀測，截至2013年底，應(yīng)變觀測資料連續(xù)、可靠，內(nèi)在精度較高。

圖1 九江臺定點形變觀測及井孔位置Fig.1 The sketch map of the fi xed deformation observation points and the wells at Jiujiang Seismic Station

表1 九江臺各井孔基本參數(shù)Table 1 The basic parameters of the wells at Jiujiang Seismic Station

2 儀器觀測原理

2.1 伸縮儀觀測原理

SS-Y型伸縮儀主要用于洞體應(yīng)變固體潮及地震前兆地體應(yīng)變監(jiān)測，也可用于大型精密工程、大型建筑、大壩等應(yīng)變測量。該

儀器測量地殼表面兩點間的應(yīng)變量（呂寵吾等，2001），即

其中：L為原地殼表面兩點間的距離（即基線長）；L′為變化后地殼表面兩點間的距離，ΔL為基線的變化量，ε為線應(yīng)變量，即單位長度的相對變化量。

2.2 體應(yīng)變、面應(yīng)變與線應(yīng)變

一般，在純壓力（或張力）作用下，物體內(nèi)部將形成正應(yīng)力，這種正應(yīng)力，將給物體造成線應(yīng)變。物體長度為L，在外力的作用下相對于L的形變量為ΔL，則比值為線應(yīng)變；若截面積為S，在外力作用下，相對面積的變量為ΔS，則其比值為面應(yīng)變。若物體體積為V，在外力作用產(chǎn)生的體積變化為ΔV，則比值為體應(yīng)變。

若在同一水平面內(nèi)，測得EW、NS向的線應(yīng)變?yōu)棣臠EW與εLNS，則在該水平面上的面應(yīng)變?yōu)椋骸?/p>

地表處的體應(yīng)變計算公式（劉序儼等，1988）為

上述各式中，ε和θ小于0時表示壓縮，反之表示為拉張。

3 資料分析

3.1 洞體溫度與洞體應(yīng)變

洞體溫度是洞體應(yīng)變觀測的一個輔助測項，觀測精度為0.001℃，溫度傳感器置于觀測基線腔體內(nèi)。從圖2可見，洞體應(yīng)變與洞體溫度均有清晰的年變動態(tài)，經(jīng)計算，線應(yīng)變NS、EW向與洞體溫度相關(guān)系數(shù)分別為0.7和0.85；溫度升高應(yīng)變增大，呈拉張狀態(tài)，溫度降低應(yīng)變減小，呈壓縮狀態(tài)。洞體溫度年變化幅度約1.0℃，應(yīng)變變化在10-6量級。頻譜分析結(jié)果表明，洞體溫度和線應(yīng)變EW向卓越周期為365.2天，線應(yīng)變NS向卓越周期為1 826天（5年），其365.2天的周期成分亦很明顯（圖3）。因此，洞體應(yīng)變長期趨勢與洞體溫度變化相關(guān)。幾年來，線應(yīng)變NS、EW向?qū)Χ大w溫度的平均變化率為1.36×10-6/℃和4.50×10-6/℃，均遠大于伸縮儀基線材料因瓦介質(zhì)的線膨脹系數(shù)，可見線應(yīng)變觀測準(zhǔn)確反映了基線長度的變化。

圖2 2009—2013年洞體應(yīng)變與洞體溫度變化Fig.2 The cave strain and cave temperature from2009 to 2013 at Jiujiang Seismic Station

圖3 九江臺洞體溫度與洞體應(yīng)變頻譜分析結(jié)果Fig.3 The cave temperature and cave strain spectrumresults at Jiujiang Seismic Station

3.2 氣壓持續(xù)上升與洞體應(yīng)變

九江臺氣壓持續(xù)上升時，洞體應(yīng)變會發(fā)生畸變，表現(xiàn)為線應(yīng)變NS、EW向同步減小，氣壓停止上升開始下降時，線應(yīng)變逐步回歸以往趨勢。表2統(tǒng)計了2010—2013年氣壓上升與線應(yīng)變同步變化的有關(guān)參數(shù)，氣壓的持續(xù)上升，線應(yīng)變和體應(yīng)變下降，應(yīng)變呈壓縮變化，反之，氣壓下降，應(yīng)變呈拉張變化。氣壓每增加1 hPa時，NS、EW向線應(yīng)變和體應(yīng)變分析平均下降18.94、10.82和19.84個10-10。

3.3 集中降雨與洞體應(yīng)變

九江臺日降雨量在40 mm以上時，洞體應(yīng)變記錄到明顯的異常變化。降雨開始后，經(jīng)過一定的時間延遲，水位上升，線應(yīng)變NS、EW向同步下降，直至降雨量變小，水位和洞體應(yīng)變同步停止變化。表3統(tǒng)計了2010—2013年降雨引起洞體應(yīng)變與水位同步變化情況，每降雨量1 mm，水位平均上升0.8 mm，NS、EW向線應(yīng)變與體應(yīng)變平均下降13.2、3.3和11個10-10，降雨使應(yīng)變呈壓縮變化。

表2 九江臺氣壓持續(xù)升高與應(yīng)變變化Table 2 The strain resulted in the air pressure raising at Jiujiang Seismic Station

3.4 抽水與洞體應(yīng)變

九江臺抽水井距山洞中心約20 m，我們于2013年7月30—31日分別進行2小時抽水實驗。兩次實驗結(jié)果重復(fù)性好，見圖4，有關(guān)數(shù)據(jù)見表4。從圖4可見，抽水實驗時兩日線應(yīng)變NS、EW向變化趨勢一致，抽水井蓄水池的水位降深也基本一致，而流體井水位未見明顯變化。線應(yīng)變NS向隨抽水期趨勢上升，抽水結(jié)束停止上升，線應(yīng)變EW向隨抽水期下降，在停止抽水后繼續(xù)下降約2小時。可以看出，抽水時呈現(xiàn)NS向拉張、EW向壓縮的反向變化，兩者在恢復(fù)到抽水前狀態(tài)時基本同步。兩日的體應(yīng)變變化幅度接近，總體應(yīng)變呈拉張狀態(tài)。

表3 九江臺降雨時應(yīng)變與水位變化統(tǒng)計Table 3 The strain and water level resulted in the rainfalls at Jiujiang Seismic Station

圖4 九江臺抽水與應(yīng)變變化實驗曲線

表4 九江臺抽水實驗時應(yīng)變變化統(tǒng)計Table 4 The strain variations from the pumping water tests at Jiujiang Seismic Station

3.5 承壓含水層卸載時水位與洞體應(yīng)變

2008年9月8日起，距流體井10 m處，新鉆一口水化井，孔徑為130 mm，終孔71m，11月3日停止鉆孔。9月23日打穿第1個含水層，17—22 m，有少量水溢出，水位和線應(yīng)變無明顯變化；10月13日進入第2個含水層，58—63 m，大量水涌出，流量可達15 t/小時，水位開始大幅下降，線應(yīng)變NS、EW向同步減小，至10月25日，水位和線應(yīng)變停止下降，見圖5。根據(jù)水位下降速率的差異，可把這次承壓含水層卸載分為突然釋放、平穩(wěn)釋放和剩余釋放3個階段，水位變化速率之比為9.66/2.95/1，而體應(yīng)變—水位速率之比為0.59/0.72/1，表明水位導(dǎo)致應(yīng)變變化是穩(wěn)定的，水位平均每下降1 mm，可引起體應(yīng)變下降12.34個10-10，見表5。承壓含水層水位下降，導(dǎo)致洞體應(yīng)變降低，即呈壓縮變化。

圖5 承壓含水層卸載引起的水位和應(yīng)變變化Fig.5 The strain and water level variations during the confi ned aquifer unloading at Jiujiang Seismic Station

表5 承壓含水層卸載引起的水位與應(yīng)變變化統(tǒng)計Table 5 The strain and water level variations resulted in unloading of the confi ned aquifer at Jiujiang Seismic Station

4 洞體應(yīng)變異常

4.1 異常量級

由九江臺觀測環(huán)境因素變化對應(yīng)變的影響分析可知，臺站抽水導(dǎo)致的應(yīng)變量在10-8量級，集中降雨和氣壓變化引起的應(yīng)變量在10-9—10-7，承壓含水層卸載和洞體溫度年變化引起的應(yīng)變量量級為10-6，根據(jù)各種環(huán)境變化引起的應(yīng)變量繪制變化范圍，見圖6。由于應(yīng)變固體潮日變幅度在10-8量級，因此，這些觀測環(huán)境因素的變化均能使應(yīng)變固體潮觀測記錄疊加干擾，導(dǎo)致固體潮發(fā)生畸變，在對該資料排除進行異常分析時需干擾。

4.2 靜力學(xué)模型的建立

圖6 觀測環(huán)境變化引起的應(yīng)變量級范圍Fig.6 The strain magnitude ranges resulted inthe observation conditions changing

伸縮儀實質(zhì)上是對山洞里固定基線長度即一維變量進行測量，通過格值轉(zhuǎn)換計算，獲得地表的線應(yīng)變。材料的長度變化，取決于材料所處熱狀態(tài)和應(yīng)力狀態(tài)，如果該材料為含水介質(zhì)，則還與介質(zhì)含水度有關(guān)，洞體應(yīng)變觀測基線，既滿足一般材料的屬性，還應(yīng)考慮介質(zhì)含水度變化的影響。把九江臺洞體應(yīng)變觀測抽象成一個簡單的靜力學(xué)模型，見圖7（a），來分析觀測環(huán)境因素變化是如何導(dǎo)致應(yīng)變異常變化的。山洞作為被測體，上部是山體，由灰?guī)r和覆蓋其上的粘土層或砂礫層組成，表層植被豐富，山體含有裂隙，被測體下部依次是較薄的潛水層及由上、下底板夾持的承壓含水層構(gòu)造，承壓含水層厚達數(shù)米，儲水量豐富。根據(jù)化繁為簡的處理方法，進行某個環(huán)境因素的應(yīng)變異常分析，應(yīng)盡量選擇其他因素不變時，以消除相互影響，從靜力學(xué)角度解釋觀測結(jié)果的成因。

圖7 九江臺洞體應(yīng)變觀測靜力學(xué)模型示意Fig.7 The static mechanical model sketch map of the cave strain observation at Jiujiang Seismic Station

4.3 應(yīng)變異常變化與溫度效應(yīng)

洞體溫度受季節(jié)變化控制，太陽光輻射和地球運動使地殼表面和大氣的溫度發(fā)生變化，是形變觀測資料具有年變動態(tài)的直接原因（陳德福，1993）。把山洞觀測基線看成一種彈性材料，九江臺觀測到的洞體溫度年變化范圍為（16.7±0.5）℃，溫度傳感器安置于測量基線腔體內(nèi)，基線所含介質(zhì)溫度即所測量的洞體溫度，基線介質(zhì)材料不會發(fā)生相變，應(yīng)遵循熱脹冷縮原理。相關(guān)分析和譜分析結(jié)果表明，線應(yīng)變與洞體溫度變化密切相關(guān)，經(jīng)計算，2009—2013年NS、EW向基線熱膨脹系數(shù)分別為1.36×10-6/℃和4.50×10-6/℃，而實際灰?guī)r在常溫下熱膨脹系數(shù)約（3—5）×10-6/℃，二者具有可比性，可見，基線介質(zhì)的熱效應(yīng)可以使基線介質(zhì)產(chǎn)生這種幅度變化的應(yīng)變。因此，九江臺洞體應(yīng)變具有年變動態(tài)的直接原因是洞體溫度效應(yīng)，其形態(tài)規(guī)律性強，可用線性回歸辦法來消除這種環(huán)境因素帶來的影響。

4.4 體應(yīng)變異常變化與靜力效應(yīng)

由應(yīng)力應(yīng)變關(guān)系可知，在介質(zhì)彈性范圍內(nèi)，體應(yīng)變與體應(yīng)力有以下關(guān)系式（考慮到應(yīng)變受壓為負(fù)）

其中，θ為體應(yīng)變，Θ為體應(yīng)力，μ為泊松比，E為彈性模量，且有σx、σy、σz分別為3個方向的主應(yīng)力。給體應(yīng)力一個增量ΔΘ，會產(chǎn)生Δθ，這時有

可得

當(dāng)九江臺氣壓經(jīng)數(shù)天持續(xù)升高時，使體應(yīng)變降低，表現(xiàn)為壓縮，見圖7（b）。可以視為在山體3個方向上增加同樣的主應(yīng)力，主應(yīng)力增量可以用氣壓增量表示。

一般認(rèn)為，降雨時通過巖體節(jié)理或裂隙下滲含水層和降雨積水荷載效應(yīng)影響區(qū)域應(yīng)力場（趙小茂等，2009）。九江臺山體為灰?guī)r區(qū)，存在較多的裂隙或孔隙，為上部流體的滲入創(chuàng)造了條件，雨水一部分通過地表徑流散失，一部分直接入滲山體巖層，被測體受雨水荷載附加垂向應(yīng)力，應(yīng)變表現(xiàn)為受壓，見圖7（c）。考慮山體承受降水面積為50×50 m2，觀測基線NS、EW向分別長9.04、11.92 m，山體降水20%為有效水體，其余通過地表徑流散失，且僅考慮垂向作用力。

九江臺抽水井是泉水井，抽水試驗期間兩個小時排水量約60 m3，一部分出水來自潛水面以下，一部分應(yīng)來自山體巖層孔隙，被測體上部必然有一部分重量損失。由于存在抽水漏斗，使補給區(qū)到井眼處水力梯度依次增大，線應(yīng)變EW向在抽水停止2小時后才達到最大變形量，表明EW基線方向是主要的補給通道，續(xù)補過程顯然是基線介質(zhì)不斷失水的過程，失水后的收縮效應(yīng)使線應(yīng)變EW從抽水開始后一直下降。抽水過程中體應(yīng)變上升，表明被測體綜合表現(xiàn)為受拉張影響，見圖7（d）。考慮總排出水量的50%來自于山體，被測體以20×20 m2計算體應(yīng)力面積，并僅考慮垂向作用力。

因新鉆井承壓含水層逐步發(fā)生應(yīng)力卸載，上層山體和地層總的重力由承壓含水層的浮力和界面上向上的托力合力抵消，從而達到靜力平衡。隨著承壓含水層卸載，靜水壓力減小，向上托力降低，上層不斷產(chǎn)生新的附加垂向作用力，被測體受附加垂直向應(yīng)力的作用而受壓，應(yīng)變大幅降低，見圖7（e）。用承壓含水層靜水壓變化來計算體應(yīng)力，與氣壓效應(yīng)類似，體應(yīng)力應(yīng)為3個方向相同的附加應(yīng)力。

對于九江臺灰?guī)r地區(qū)，取μ=0.29，E=3.71×1010Pa。采用式（6），按上述模型參數(shù)分別對氣壓、降雨、抽水和承壓含水層卸載引起的體應(yīng)變進行計算，結(jié)果見表6，表中的綜合結(jié)果為各種觀測環(huán)境因素變化引起應(yīng)變變化的實測與理論計算均值及其方差，越接近1表明與理論計算越接近。

表6 九江臺觀測環(huán)境變化引起的體應(yīng)變理論與實測結(jié)果對比Table 6 The volume strain academic and practical results resulted in the observation conditions changing at Jiujiang Seismic Station

由上述分析可看出，應(yīng)用彈性介質(zhì)靜力學(xué)模型可以較好地解釋洞體應(yīng)變狀態(tài)問題（壓縮或拉張）。在定量計算方面，氣壓和降雨導(dǎo)致的體應(yīng)變較符合觀測結(jié)果，而抽水和承壓含水層卸載引起的體應(yīng)變有一定差距，可能與采用的模型有關(guān)，可進一步采用含流體的介質(zhì)模型深入研究。

5 結(jié)論

（1）獲得九江臺洞體應(yīng)變受抽水、集中降雨、氣壓上升、承壓含水層卸載和洞體溫度年變化因素影響的應(yīng)變異常數(shù)量級的認(rèn)識，這些因素能直接影響應(yīng)變固體潮觀測。

（2）九江臺洞體溫度年變化主要通過熱脹冷縮效應(yīng)影響洞體應(yīng)變；氣壓上升、集中降雨和承壓含水層卸載主要通過被測體上部附加應(yīng)力加載引起應(yīng)變的下降，抽水過程則是上部附加應(yīng)力卸載引起應(yīng)變的上升。

（3）九江臺抽水、降雨和氣壓變化引起的應(yīng)變量較小，在短期內(nèi)得以恢復(fù)；而承壓含水層卸載引起的應(yīng)變量巨大，區(qū)域應(yīng)變場重新進行了調(diào)整。

衷心感謝周克昌研究員和何案華副研究員的精心指導(dǎo)和幫助。

陳德福.地殼形變動力學(xué)觀測與研究［M］.北京：海洋出版社，1993：42-54.

陳敏，等.重慶臺傾斜潮汐形變干擾異常與震兆異常的識別[J]. 大地測量與地球動力學(xué)，2007，Spec，27：80-83.

高昂.山洞環(huán)境對地震觀測的影響[J]. 地震地磁觀測與研究，2013，34（5/6）：200-203.

蔣駿，等.潮汐線應(yīng)變組合觀測的物理力學(xué)實質(zhì)及其調(diào)和分析[J]. 地球物理學(xué)報，1994，37（Supp I）：204-212.

江西地質(zhì)礦產(chǎn)勘查開發(fā)局.中華人民共和國地質(zhì)圖說明書，海會幅，1993，H-50-89-A.

劉序儼，等.地表的面應(yīng)變和體應(yīng)變固體潮理論值計算及其調(diào)和分析[J].地殼形變與地震，1988，4：354-358.

劉序儼，等.排除形變觀測數(shù)據(jù)中降水干擾的數(shù)學(xué)物理方法的研究[J]. 地殼形變與地震，1991，11（1）：36-40.

呂寵吾，等. SS-Y型短基線伸縮儀及其標(biāo)定裝置[J]. 地殼形變與地震，2001，21（3）：82-88.

牛安福.地殼形變觀測與地震前兆一般性問題的討論[J].國際地震動態(tài)，2007，6（342）：1-6.

王曰鳳，等.懷來地震臺線應(yīng)變觀測中的干擾及消除方法[J]. 地震地磁觀測與研究，2007，28（2）：80-83.

趙小茂，等. “集中降雨”與淺層承壓井水位及區(qū)域地震活動關(guān)系的初步研究[J].西北地震學(xué)報，2009，31（3）：259-264.