山東青州中黃山流動水準測值與流體變化關系研究

王 錚，閆德橋，趙銀剛，張小東，魏秀才

(1.山東濰坊市地震局，山東 濰坊 261061；2.山東安丘地震臺，山東 安丘 262100)

山東青州中黃山流動水準測值與流體變化關系研究

王 錚1，閆德橋2，趙銀剛1，張小東1，魏秀才1

(1.山東濰坊市地震局，山東 濰坊 261061；2.山東安丘地震臺，山東 安丘 262100)

運用1993—2014年山東青州中黃山流動水準及輔助手段的觀測資料，系統分析中黃山流動水準測值與降水量、水庫蓄水量、水庫水位、地下水位等影響因素的關系。結果表明，水準測值與水庫水位和蓄水量的相關性良好，與老井水位的相關性次之，與降水量的相關性最差。同時討論流體影響水準觀測的力學機理，認為水準測值變化的原因是，在水庫蓄水量載荷改變、地下水孔隙壓力變化和腐蝕的共同作用下，斷層上盤向下或向上滑動的結果。

流動水準；相關分析；影響機理

0 引言

地球表層是流體與固體強烈相互作用的場所，斷層是這種相互作用的主要能量與信息的傳遞通道。一般認為，斷層帶上的流體對于斷層活動和大地震的發生可起到潤滑的作用；但對于地殼淺部，則認為流體是斷層活動中的干擾因素。

長期以來，形變觀測與流體的關系也引起了一些學者的關注，相關研究指出：降水、地下水位與斷層形變、GPS垂直向測值變化間都具有比較密切的相關關系[1-2]。

流動水準測量與環境因子、地殼構造變動等多種因素有關。文章以實際觀測資料為基礎，分析山東青州中黃山流動水準測值與降水、水庫水位、地下水位等流體變化間的關系，為進一步探討流動水準測值與地震活動的關系做基礎研究工作。

1 地質環境與觀測場地概況

1.1 地質背景簡況

雙山—李家莊斷裂為觀測地區內主要的地質構造，是一條走向北西的不同地質單元的分界線，臨朐盆地的東北邊界。該斷裂也是濰坊轄區內一條重要的活動斷裂，與上五井斷裂交匯處曾發生1829年6益都地震。斷層長約28 km，走向北35°西，傾向南西，傾角50°～80°，系正斷層。沿斷裂走向的活動性質和特點可分為南段、中段和北段3個明顯不同的段落，南段和中段可見斷裂面出露。中段斷裂規模大、活動時代新，晚更新世晚期有發生斷錯地表的活動。斷裂擠壓破碎帶主要發育在太古界花崗片麻巖(Art)和寒武系地層中，新活動以左旋張扭正斷活動為主，破裂方式以粘滑為主，兼有少量蠕滑。

1.2 觀測場地和觀測資料

觀測場地位于山東青州市彌河鎮中黃山西，測線跨雙山—李家莊斷裂中段，破碎帶寬8～10 m，該斷裂在上黃山村東水庫溢洪口和中黃山村南橋西的沖溝中部有較好的剖面出露，斷裂兩側巖性為石灰巖。

中黃山流動水準觀測點西南8 km處建有黑虎山水庫，水庫位于彌河支流石溝河上，控制流域面積為190 km2，總庫容為5 646萬立方米(見第27頁圖1)。距水準觀測北端點東南約150 m處有一口深20 m左右的百年老井。

中黃山流動水準觀測由安丘地震臺負責，始于1993年，每月中旬野外觀測一次，同時開展降水、水庫蓄水及老井水位等輔助觀測。水準測線與斷層夾角35°，全長268 m，共設6站，高差為0.89 m(見圖1)。測線兩端點均建有鋼管基巖水準標石、不銹鋼標志，北端點孔深8.05 m，南端點孔深7.0 m，標桿管均下到完整基巖中。

2 水準測值與流體測項的年變化

1993—2014年中黃山流動水準及輔助觀測的年變值如第27頁圖2所示。可以看出，中黃山水準測值的年變規律較明顯，且受當地降水、黑虎山水庫水位和蓄水量、地下水位的影響。

2.1 典型年變分析

選取2013年的資料做年變分析。中黃山水準測值、降水量、老井水位、黑虎山水庫水位和蓄水量月均值或月累計的變化曲線如圖3所示，相關分析結果如表1和第28頁表2所示。可以看出：

圖1 中黃山流動水準測線布設示意圖Fig.1 Arrangement of mobile leveling survey lines in Zhonghuangshan

圖2 1993—2014年青州中黃山流動水準及輔助測項年變圖Fig.2 Annual variation of mobile leveling and auxiliary items in Zhonghuangshan in Qingzhou from 1993 to 2014

(1) 水準測值、降水量、老井水位、水庫水位和蓄水量都有明顯的年變特點。其中，降水等變化為“雨季突升”型，水準測值變化為“雨季突降”型。

(2) 水準測值與降水量、老井水位、水庫水位和蓄水量之間，降水量等4個測項值越大，水準測值就越

圖3 2013年1—12月青州中黃山降水量、老井水位、黑虎山水庫水位和蓄水量變化曲線圖Fig.3 Change of precipitation in Zhonghuangshan, level of old wells, water level and storage capacity of Heihushan Reservoir from January to December in 2013

小；反之則越大。

(3) 老井水位受上月降水量的影響明顯，符合水文地質學原理，但兩者的相關性卻不顯著；水庫水位和蓄水量受上月降水量的影響明顯，但之間的相關關系也不成立；水庫的蓄水量與水位之間的相關關系較密切(見表1)。

表1 2013年1—12月上月降水量與老井水位、水庫水位、水庫蓄水量相關分析結果

(4) 水準測值與水庫水位和蓄水量的相關性較好(見第28頁圖4)，與老井水位的相關性次之，與降水量的相關性最差(見表2)。表明，降水為引起水準測值變化的成因性因素，水庫水位和蓄水量是最直接的影響因素。此結論與未做數據處理之前的不完全相同，以前主要強調干旱(即降水)對水準觀測的影響。

表2 2013年1—12月水準測值與上月降水量、老井水位、水庫水位和蓄水量相關分析結果

圖4 2013年1—12月水準測值與水庫蓄水量相關圖Fig.4 Correlation between leveling values and reservoir storage capacity from January to December in 2013

2.2 逐年變化分析

對1993—2014年水準測值與降水量、水庫蓄水量的年變關系逐年做分析。得出：

(1) 水準測值與上月降水量之間，絕大多數年份的相關性一般或較差。1993—2014年，兩者的相關系數最小為0.007 6，最大為0.551，平均為0.263 7。僅2014年的相關性較好，相關系數為0.634 2(樣本數n-2=10，可信度α=0.05的r值為0.576)。

(2) 水準測值與水庫蓄水量之間存在較好的相關關系(見表3)。在22年中，有15年的相關性良好，占總數的68.2%；7年的相關性較差或反常(見圖5)，占31.8%，說明水庫蓄水量對水準測值的影響更為直接。

表3 1993—2014年水準測值與水庫蓄水量逐年相關分析結果

圖5 2005年水準測值與降水、水庫蓄水量關系曲線圖Fig.5 Relationship of leveling values with rainfall and reservoir storage capacity in 2005

(3) 水庫蓄水量的大小，受降水量和人為用水量雙重因素的影響。“干旱”實為降水量的顯著偏少，而人為用水量又與降水量和經濟發展狀況兩方面因素有關，且關系復雜。

2.3 逐年水準增量最小值與降水量、水庫蓄水增量關系分析

為研究水準值增量與降水量、水庫蓄水增量的關系，采用列表分析法。列出1年中≥100 mm的月降水量，顯著變化前后的水庫蓄水增量，1年中水準測值增量的最小值(見表4)。對比得出：

(1) 在大多數年份，水準測值增量最小值出現在7至9月，少數出現在6月、10月，極個別在11月(2010年)、12月(2014年)。

(2) 水準測值增量的減少，最小為-0.11，最大為-3.74。

(3) 在大多數年份，水準測值與降水量、水庫蓄水量之間為正相關關系。

(4) 月降水量顯著增加(即≥100 mm)，可使下一月或再下月水庫的蓄水增量顯著增加(不一定轉為正值)，即蓄水量顯著增加滯后于降水量增加一、二個月。蓄水量顯著增加后，可使當月(個別年份為次月)的水準值增量最大幅度下降而達到最小值。

(5) 2001、2005、2006、2009、2014年的變化規律反常，占總年數的22.7%。如，2005年8月水準測值與水庫蓄水量的關系，由負相關轉變為正相關(見圖5)；2014年全年未出現≥100 mm，但11月份卻出現水準測值增量的最小值-0.11(見表4)。

表4 1993—2013年部分水準增量最小值與月降水量、水庫蓄水量對比

3 水準測值增量與流體測項的多年變化

1994—2014年，中黃山水準年均值增量、年降水量、黑虎山水庫水位和蓄水量年均值增量的變化曲線如第30頁圖6所示，之間的相關分析結果如第30頁表5所示。分析得出：

(1) 降水是水庫水位和蓄水量最重要的影響因素(水庫水的補給來源)，計算結果表明，年降水量與水庫水位和蓄水量年均值增量之間的相關性較差。

(2) 1994—2012年，水準年均值增量與水庫水位和蓄水量年均值增量間的相關性較好(呈負相關)。

(3) 1994—2013年、1994—2014年，水準年均值增量與水庫蓄水量年均值增量間不相關。這與2013年水準測值增量太大、2014年水準測值增量出現反常下降等有關。

(4) 1994—2014年，水準年均值增量與年降水量之間的相關性都較差。

以上分析表明，降水與水準測值變化間的相關性并不良好，降水是通過水庫水位和蓄水量的變化來引起水準測值變化的。

圖6 1994～2014年降水量年累計、水庫水位和蓄水量年均值增量、水準年均值增量的變化曲線圖Fig.6 The change of annual total precipitation, average annual increment of reservoir water level and storage capacity, average annual increment of leveling from 1994 to 2014

4 影響機理分析

當大氣中降水落到地面形成地表水和地下水后，才能影響流動水準的測量。從幾種因素的相互關系上看，降水的作用是雙重和連鎖式的，可造成水庫蓄水量和人為用水量的增減，人為用水需求量的改變又直接影響老井水位和水庫的庫容量。

上述分析表明，水庫水位和蓄水量、井水位與水準測值之間有著直接的密切關系。下面討論其相互影響的機理。

(1) 水庫(地表水)對水準測值的影響。

水庫的水體載荷改變(包括水庫兩岸被河水飽和或疏干后巖層載荷的改變)對地層產生加壓或減壓的作用，進而引起斷層上盤滑動的方向和速度的變化，這是最主要的影響。若降水量增大，對黑虎山水庫的補給量亦增大，并可因降水量充足而減少水庫的輸出量，故水庫蓄水量增加。其結果是，作用于斷層上盤(南

表5 水庫水位和蓄水量、水準年增量與年降水量的相關分析結果

盤)的水庫荷載增大，使上盤下降，水準測值顯著減小。反之，若降水量減小，對黑虎山水庫的補充量亦減少，又因干旱使水庫輸出量增大，故水庫蓄水量減小。其結果是，作用于斷層上盤(南盤)的水庫荷載減輕，使上盤上升，水準測值亦顯著上升。

近年來，降水量減小和水庫輸出量增大，斷層上盤上升，水準測值的增加更顯著(見圖2)。

(2) 地下水對水準測值的影響。

當降水補給引起地下水位上升時，含水層巖石的孔隙壓力增加，斷層面上的有效正壓力減小，有利于斷層滑動；反之，當水位下降時，孔隙壓力減小，則可抑制斷層的滑動[3-4]。另外，水位的升降變化，通過化學腐蝕作用，可導致巖石強度的減小或增強，從而有利于斷層滑動或者抑制斷層滑動。

需說明的是，在流動水準測線一帶，雙山—李家莊斷裂兩側巖性為石灰巖，即斷層上下盤巖性無明顯差異，斷層兩盤含水層巖石的孔隙度、儲水性能的差異不大。因而由于水位升降、巖體總重量作用于斷層面上的附加剪切分力的變化不大，故地下水位變化而導致的含水層巖石的荷載作用可不予考慮。

總之，中黃山流動水準觀測值變化的原因是，水庫蓄水量載荷改變與地下水孔隙壓力變化和地下水腐蝕等共同作用的結果。

5 結論與討論

(1) 中黃山水準測值與水庫水位和蓄水量的相關性最好，與老井水位的次之，與降水量的最差。表明，對于中黃山流動水準觀測，降水為引起變化的成因性因素，水庫水位和蓄水量是最直接的影響因素。這一點與以前的認識不盡相同，過去只籠統強調干旱的影響。

(2) 月降水量增加(即≥100 mm)后一、二個月，水庫蓄水量亦顯著增加；蓄水量顯著增加之后，可使當月(個別年份為次月)的水準值增量值出現一年中的最大幅度減小。

(3) 水庫蓄水量的大小，受降水量和人為用水量雙重因素的影響，而人為用水量則與降水量和經濟發展狀況兩方面因素有關，之間的關系較復雜。

(4) 水準測值受影響的機理是，降水引起老井水位、水庫水位和蓄水量(還受人為用水的影響)的變化，水庫的蓄水量載荷改變，對地層產生加壓或減壓，致使斷層上盤向下或向上滑動，進而引起水準測值變化；地下水位的升降變化，通過改變孔隙壓力的大小和腐蝕作用的強弱來促進或抑制斷層滑動，也影響水準測值的改變。

(5) 存在問題是，少數年份水準測值與水庫蓄水量年變化的相關性、極少數年份水準年均值增量與水庫水位和蓄水量年均值增量之間的相關性一般或較差。這與降水量、水庫水位和蓄水量、人為用水量彼此之間的關系復雜有關，但目前仍未尋找到具體的原因，有待今后開展此方面的研究。

(6) 流體引起斷層滑動的深度，是否存在淺層和深部都有，不得而知，需加強研究。

[1] 黃輔瓊，陳 顆，白長清，等.八寶山斷層的變形行為與降雨及地下水的關系[J].地震學報,2005，27(6)：637-646.

[2] 顧申宜，王慧琳，張 帆.海南瓊中臺GPS垂直分量變化與降水、地下水位關系研究[J].大地測量與地球動力學，2014,43(6)：37-40.

[3] 郭增建，秦保燕.震源物理[M].北京:地震出版社，1979：20-21，159.

[4] 王吉易，宋貫一，曹志成，等.地下水誘發的淺層前兆異常及其機理與有關的地震預報問題(1)[J].華北地震科學，2002，20(2)：27-41.

Relationship between Measured Values of Mobile Leveling andFluid Change in Zhonghuangshan in Qingzhou, Shandong

WANG Zheng1, YAN De-qiao2, ZHAO Yin-gang2, ZHANG Xiao-dong1, WEI Xiu-cai1

(1.Earthquake Administration of Weifang City, Weifang, Shandong 261041, China;2.Anqiu Seismological Station, Anqiu, Shandong, 262100, China)

The relationship of measured values of the moblie leveling in Zhonghuangshan and influence factors of precipitation, storage capacity of the reservoir, water level of the reservoir and groundwater level is analyzed by using Zhonghuangshan mobile leveling observation data and auxiliary observation data from 1993 to 2014 in Qingzhou, Shandong. The results show that the measured values of leveling has the best correlation with reservoir water level and water storage capacity, and better correlation with the water level of old wells, the poor correlation with precipitation. The mechanical mechanism of fluid influence to leveling observation is discussed. It is considered that the reason for the change of measured values is the sliding up or down of the hanging wall of the fault under the influence of load change of reservoir storage capacity, change of pore pressure and corrosion.

Mobile leveling; Correlation analysis; Effect mechanism

1000-6265(2017)01-0026-06

2016-05-20

山東省地震局合同制科研項目(15Y33)。

王 錚(1977— )，女，山東省濰坊市人。1999年畢業于中國地質大學(北京)，工程師。

P315.63

A