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

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

山東青州中黃山流動(dòng)水準(zhǔn)測(cè)值與流體變化關(guān)系研究

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

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

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

流動(dòng)水準(zhǔn)；相關(guān)分析；影響機(jī)理

0 引言

地球表層是流體與固體強(qiáng)烈相互作用的場(chǎng)所，斷層是這種相互作用的主要能量與信息的傳遞通道。一般認(rèn)為，斷層帶上的流體對(duì)于斷層活動(dòng)和大地震的發(fā)生可起到潤(rùn)滑的作用；但對(duì)于地殼淺部，則認(rèn)為流體是斷層活動(dòng)中的干擾因素。

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

流動(dòng)水準(zhǔn)測(cè)量與環(huán)境因子、地殼構(gòu)造變動(dòng)等多種因素有關(guān)。文章以實(shí)際觀測(cè)資料為基礎(chǔ)，分析山東青州中黃山流動(dòng)水準(zhǔn)測(cè)值與降水、水庫水位、地下水位等流體變化間的關(guān)系，為進(jìn)一步探討流動(dòng)水準(zhǔn)測(cè)值與地震活動(dòng)的關(guān)系做基礎(chǔ)研究工作。

1 地質(zhì)環(huán)境與觀測(cè)場(chǎng)地概況

1.1 地質(zhì)背景簡(jiǎn)況

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

1.2 觀測(cè)場(chǎng)地和觀測(cè)資料

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

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

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

2 水準(zhǔn)測(cè)值與流體測(cè)項(xiàng)的年變化

1993—2014年中黃山流動(dòng)水準(zhǔn)及輔助觀測(cè)的年變值如第27頁圖2所示?？梢钥闯?，中黃山水準(zhǔn)測(cè)值的年變規(guī)律較明顯，且受當(dāng)?shù)亟邓?、黑虎山水庫水位和蓄水量、地下水位的影響?/p>

2.1 典型年變分析

選取2013年的資料做年變分析。中黃山水準(zhǔn)測(cè)值、降水量、老井水位、黑虎山水庫水位和蓄水量月均值或月累計(jì)的變化曲線如圖3所示，相關(guān)分析結(jié)果如表1和第28頁表2所示?？梢钥闯觯?/p>

圖1 中黃山流動(dòng)水準(zhǔn)測(cè)線布設(shè)示意圖Fig.1 Arrangement of mobile leveling survey lines in Zhonghuangshan

圖2 1993—2014年青州中黃山流動(dòng)水準(zhǔn)及輔助測(cè)項(xiàng)年變圖Fig.2 Annual variation of mobile leveling and auxiliary items in Zhonghuangshan in Qingzhou from 1993 to 2014

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

(2) 水準(zhǔn)測(cè)值與降水量、老井水位、水庫水位和蓄水量之間，降水量等4個(gè)測(cè)項(xiàng)值越大，水準(zhǔn)測(cè)值就越

圖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) 老井水位受上月降水量的影響明顯，符合水文地質(zhì)學(xué)原理，但兩者的相關(guān)性卻不顯著；水庫水位和蓄水量受上月降水量的影響明顯，但之間的相關(guān)關(guān)系也不成立；水庫的蓄水量與水位之間的相關(guān)關(guān)系較密切(見表1)。

表1 2013年1—12月上月降水量與老井水位、水庫水位、水庫蓄水量相關(guān)分析結(jié)果

(4) 水準(zhǔn)測(cè)值與水庫水位和蓄水量的相關(guān)性較好(見第28頁圖4)，與老井水位的相關(guān)性次之，與降水量的相關(guān)性最差(見表2)。表明，降水為引起水準(zhǔn)測(cè)值變化的成因性因素，水庫水位和蓄水量是最直接的影響因素。此結(jié)論與未做數(shù)據(jù)處理之前的不完全相同，以前主要強(qiáng)調(diào)干旱(即降水)對(duì)水準(zhǔn)觀測(cè)的影響。

表2 2013年1—12月水準(zhǔn)測(cè)值與上月降水量、老井水位、水庫水位和蓄水量相關(guān)分析結(jié)果

圖4 2013年1—12月水準(zhǔn)測(cè)值與水庫蓄水量相關(guān)圖Fig.4 Correlation between leveling values and reservoir storage capacity from January to December in 2013

2.2 逐年變化分析

對(duì)1993—2014年水準(zhǔn)測(cè)值與降水量、水庫蓄水量的年變關(guān)系逐年做分析。得出：

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

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

表3 1993—2014年水準(zhǔn)測(cè)值與水庫蓄水量逐年相關(guān)分析結(jié)果

圖5 2005年水準(zhǔn)測(cè)值與降水、水庫蓄水量關(guān)系曲線圖Fig.5 Relationship of leveling values with rainfall and reservoir storage capacity in 2005

(3) 水庫蓄水量的大小，受降水量和人為用水量雙重因素的影響?！案珊怠睂?shí)為降水量的顯著偏少，而人為用水量又與降水量和經(jīng)濟(jì)發(fā)展?fàn)顩r兩方面因素有關(guān)，且關(guān)系復(fù)雜。

2.3 逐年水準(zhǔn)增量最小值與降水量、水庫蓄水增量關(guān)系分析

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

(1) 在大多數(shù)年份，水準(zhǔn)測(cè)值增量最小值出現(xiàn)在7至9月，少數(shù)出現(xiàn)在6月、10月，極個(gè)別在11月(2010年)、12月(2014年)。

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

(3) 在大多數(shù)年份，水準(zhǔn)測(cè)值與降水量、水庫蓄水量之間為正相關(guān)關(guān)系。

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

(5) 2001、2005、2006、2009、2014年的變化規(guī)律反常，占總年數(shù)的22.7%。如，2005年8月水準(zhǔn)測(cè)值與水庫蓄水量的關(guān)系，由負(fù)相關(guān)轉(zhuǎn)變?yōu)檎嚓P(guān)(見圖5)；2014年全年未出現(xiàn)≥100 mm，但11月份卻出現(xiàn)水準(zhǔn)測(cè)值增量的最小值-0.11(見表4)。

表4 1993—2013年部分水準(zhǔn)增量最小值與月降水量、水庫蓄水量對(duì)比

3 水準(zhǔn)測(cè)值增量與流體測(cè)項(xiàng)的多年變化

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

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

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

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

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

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

圖6 1994～2014年降水量年累計(jì)、水庫水位和蓄水量年均值增量、水準(zhǔn)年均值增量的變化曲線圖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 影響機(jī)理分析

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

上述分析表明，水庫水位和蓄水量、井水位與水準(zhǔn)測(cè)值之間有著直接的密切關(guān)系。下面討論其相互影響的機(jī)理。

(1) 水庫(地表水)對(duì)水準(zhǔn)測(cè)值的影響。

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

表5 水庫水位和蓄水量、水準(zhǔn)年增量與年降水量的相關(guān)分析結(jié)果

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

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

(2) 地下水對(duì)水準(zhǔn)測(cè)值的影響。

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

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

總之，中黃山流動(dòng)水準(zhǔn)觀測(cè)值變化的原因是，水庫蓄水量載荷改變與地下水孔隙壓力變化和地下水腐蝕等共同作用的結(jié)果。

5 結(jié)論與討論

(1) 中黃山水準(zhǔn)測(cè)值與水庫水位和蓄水量的相關(guān)性最好，與老井水位的次之，與降水量的最差。表明，對(duì)于中黃山流動(dòng)水準(zhǔn)觀測(cè)，降水為引起變化的成因性因素，水庫水位和蓄水量是最直接的影響因素。這一點(diǎn)與以前的認(rèn)識(shí)不盡相同，過去只籠統(tǒng)強(qiáng)調(diào)干旱的影響。

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

(3) 水庫蓄水量的大小，受降水量和人為用水量雙重因素的影響，而人為用水量則與降水量和經(jīng)濟(jì)發(fā)展?fàn)顩r兩方面因素有關(guān)，之間的關(guān)系較復(fù)雜。

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

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

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

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

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

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

[4] 王吉易，宋貫一，曹志成，等.地下水誘發(fā)的淺層前兆異常及其機(jī)理與有關(guān)的地震預(yù)報(bào)問題(1)[J].華北地震科學(xué)，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

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

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

P315.63

A