盛艷蕊 張子廣 張素欣 楊歧焱 于春頌

1)河北省地震局，石家莊市槐中路262號(hào) 050021

2)河北省地震局秦皇島中心臺(tái)，河北 秦皇島 066100

0 引言

地下水可靈敏地記錄到地殼微小的體應(yīng)變，反映地震孕育信息，但水位的變化也包含其它干擾因素。干擾排除、前兆異常識(shí)別及異常特征的分析，是進(jìn)行震情判定工作的基礎(chǔ)。井水位異常的主要特征有多年趨勢(shì)上升、下降以及水位破年變或固體潮畸變等(國(guó)家地震局科技監(jiān)測(cè)司，1995)，其中井水位破年變異常是中強(qiáng)震前中短期異常的重要表現(xiàn)形式之一，分析水位破年變異常對(duì)震情判定具有重要意義。

井水位破年變異常是指水位年變形態(tài)遭到破壞(國(guó)家地震局科技監(jiān)測(cè)司，1995;周翠英等，2010)，總結(jié)《中國(guó)震例》中水位破年變主要類型有年變幅改變、年變規(guī)律消失、改變多年基值、年動(dòng)態(tài)改變(反向變化、上升或下降時(shí)間提前或推遲)、上升或下降速率改變等，其中反向變化、速率改變的破年變異常對(duì)應(yīng)地震的信度較高。井水位破年變異常多出現(xiàn)在地震前十幾天～1年左右，屬中短期異常，一般在地震發(fā)生后異常逐步結(jié)束;破年變異常對(duì)應(yīng)的地震震中距多在 300km范圍左右(張肇誠(chéng)等，1988、1990a、1990b、1999、2000、2002;陳棋福等，2002a、2002b、2003)。

分析昌黎井水位年動(dòng)態(tài)變化，發(fā)現(xiàn)2000年以來(lái)井水位破年變異常頻次增多。近年來(lái)由于經(jīng)濟(jì)的發(fā)展，環(huán)境干擾對(duì)觀測(cè)井的影響增多，井水位破年變是環(huán)境干擾引起還是構(gòu)造活動(dòng)引起的地震前兆異常，需要進(jìn)一步分析。本文利用動(dòng)態(tài)年周期方法識(shí)別2000～2013年昌黎井水位破年變異常，基于前人研究成果，主要從降雨、井孔觀測(cè)條件改變、附近井孔抽水、孕震機(jī)理等方面分析引起昌黎井水位破年變異常的原因。

1 昌黎井井孔概況

昌黎井地處燕山山前平原，飲馬河沖積扇上，距離昌黎-寧河斷裂約1.5km(圖1)，位于斷裂的西北盤(上升盤)，為地?zé)岙惓^(qū)溫泉鉆孔，成井深度301.04m，主要出水段為172～180m(圖2)，地下水類型為裂隙、孔隙混合承壓水(曹新來(lái)等，2001)。昌黎井處于低山丘陵水文地質(zhì)區(qū)，北部山區(qū)到井區(qū)的地形坡度較大，是該井裂隙水的補(bǔ)給區(qū)(車用太等，2004a)。昌黎井水位具有明顯的年變特征:每年3月底至4月初水位下降，5月份前后一段時(shí)間下降速率較大，雨季來(lái)臨水位開(kāi)始回升，一般在降雨充分的7～8月水位回升較快，隨后緩慢上升，直至翌年3月份達(dá)年最高水位，雨季水位回升幅度受降雨量大小影響(車用太等，2004b);這種變化的動(dòng)態(tài)說(shuō)明含水層受春灌影響水頭降低，雨季水位又迅速回升則說(shuō)明補(bǔ)給區(qū)不是很遠(yuǎn)。2004年7月洗井后水位仍具有年變特征。

圖1 昌黎井周邊地質(zhì)構(gòu)造簡(jiǎn)圖

2 昌黎井水位破年變異常判別

本文利用動(dòng)態(tài)年周期方法識(shí)別昌黎井水位破年變異常，依據(jù)《中國(guó)震例》中總結(jié)的破年變異常特征進(jìn)行分析。

2.1 動(dòng)態(tài)年周期方法

h(t)為井水位年動(dòng)態(tài)變化為某一相對(duì)時(shí)段井水位平均值，作為標(biāo)準(zhǔn)年動(dòng)態(tài)，用該時(shí)段的標(biāo)準(zhǔn)年動(dòng)態(tài)和實(shí)際觀測(cè)曲線對(duì)比，當(dāng)兩者變化相差較大時(shí)，作為異常。標(biāo)準(zhǔn)年動(dòng)態(tài)曲線根據(jù)地震前2～3年井水位變化較平穩(wěn)時(shí)期的觀測(cè)值計(jì)算得到。

式中，m為相對(duì)時(shí)段內(nèi)所含的年數(shù)，nj為第j年的觀測(cè)次數(shù)，hjt是第j年度t時(shí)刻的水位觀測(cè)值(劉喜蘭等，2003)。

圖2 昌黎井井孔柱狀圖

1983～1999年昌黎井水位多年水位均值在0.9m左右波動(dòng)，選取水位年動(dòng)態(tài)清晰、變化較為穩(wěn)定的1997～1998年的水位平均值作為標(biāo)準(zhǔn)年動(dòng)態(tài)，利用標(biāo)準(zhǔn)年動(dòng)態(tài)和實(shí)際觀測(cè)值對(duì)比，發(fā)現(xiàn)2001、2002年水位在0.5m以下波動(dòng)，未達(dá)到往年的水位高值，水位年變幅改變。2004年洗井之后水頭升高，2007年以后水位基本穩(wěn)定，在1.2m左右波動(dòng)變化，但年動(dòng)態(tài)未改變，與標(biāo)準(zhǔn)年動(dòng)態(tài)(為符合實(shí)際水位動(dòng)態(tài)變化，有些年份標(biāo)準(zhǔn)年動(dòng)態(tài)為上述時(shí)段計(jì)算的水位值+0.4m)相比，2007、2008年水位趨勢(shì)上升，上升速率改變;2009年雨季水位未回升;2010年2月中旬～5月反向上升;2012年4月反向上升;2013年10月底反向下降(圖3)。

3 昌黎井水位破年變異常分析

分析昌黎井水位破年變成因，主要從降雨、井孔觀測(cè)條件改變、附近井孔抽水以及孕震機(jī)理等方面進(jìn)行。

3.1 水位破年變與降雨關(guān)系

昌黎井含水層上部有較厚的粘土隔水層，不直接接受降雨補(bǔ)給，該井北部為燕山山脈南麓，地勢(shì)陡峭，出露大面積的花崗巖和片麻巖，主要通過(guò)地下徑流對(duì)井水位補(bǔ)給。多年觀測(cè)資料表明，正常情況下，昌黎井水位在雨季(或稍有推遲)有上升變化，因此，分析水位的破年變，可從降雨滲入量的變化入手。

圖3 昌黎井水位觀測(cè)曲線與標(biāo)準(zhǔn)年動(dòng)態(tài)水位曲線的對(duì)比

1983～2000年昌黎地區(qū)年降雨量為332～848mm(表1)，平均年降雨量589mm，雨季大幅降雨多集中在7～8月，也有年份提前到5～6月，或推遲到9～10月。1999～2002年連續(xù)4年年降雨量小于多年平均年降雨量，2001、2002年雨季最大月降雨量小于200mm(圖4)，對(duì)含水層補(bǔ)給較少，通過(guò)地下徑流對(duì)井水位補(bǔ)給也會(huì)偏少，徑流補(bǔ)給對(duì)井水位有滯后效應(yīng)(一般情況下，雨季之后水位翌年3月才能達(dá)到最高值)，對(duì)比1989、1992、1993年降雨量低于多年平均降雨量時(shí)，出現(xiàn)雨季之后水位回升較少、年變幅改變的破年變異常。因此，可以認(rèn)為2001～2002年昌黎井水位年變幅改變異常與1999年以來(lái)的降雨量顯著偏少有關(guān)。

由此可知，第二條控制性結(jié)構(gòu)面以上坡體的位移較大，以下區(qū)域基本沒(méi)有位移產(chǎn)生，第二條控制性結(jié)構(gòu)面對(duì)邊坡開(kāi)挖后的變形起到了控制作用；第一條結(jié)構(gòu)面以上坡體的位移最大，這說(shuō)明第一條結(jié)構(gòu)面以外巖體的變形最大。此外，在外界誘發(fā)因素條件下，第1條結(jié)構(gòu)面以外巖體將首先沿該結(jié)構(gòu)面發(fā)生失穩(wěn)滑動(dòng)，此外，當(dāng)外界因素繼續(xù)增加，第2條結(jié)構(gòu)面以外巖體為潛在滑體，第2條結(jié)構(gòu)面為潛在滑帶，位于此結(jié)構(gòu)面外側(cè)巖體最容易沿此結(jié)構(gòu)面發(fā)生滑動(dòng)破壞。

2007、2008年昌黎井水位趨勢(shì)上升，且上升速率較快。2000～2005年雨季降雨量均值343.5mm，2006～2008年雨季降雨量分別為468.5、454.6和384.3mm。對(duì)比發(fā)現(xiàn)，2006～2008年雨季降雨量較往年雨季降雨量均值偏大，計(jì)算2006～2008年雨季降雨量與水位高值相關(guān)系數(shù)為0.74406，說(shuō)明降雨對(duì)水位的影響較大。降雨對(duì)井水位的影響有即時(shí)效應(yīng)和滯后效應(yīng)，下面用褶積濾波方法分析降雨量對(duì)昌黎井水位的影響。

表1 1999～2002年昌黎地區(qū)年降雨量表

圖4 昌黎井水位曲線與月降雨量曲線的對(duì)比

把井孔水位、含水層看成一個(gè)系統(tǒng)，具有系統(tǒng)函數(shù)R(t)，降雨量CR(t)為輸入信號(hào)，降雨對(duì)井水位的影響值HR(t)為輸出信號(hào)，三者之間的褶積關(guān)系為:HR(t)=R(t)*CR(t)，井水位的趨勢(shì)變化h(t)，則昌黎井水位

降雨對(duì)井水位影響的即時(shí)效應(yīng)，利用一次項(xiàng)擬合，降雨對(duì)井水位的滯后效應(yīng)，利用多項(xiàng)式擬合，則降雨對(duì)井水位影響的數(shù)學(xué)模型為

假定井水位的趨勢(shì)變化h(t)為線性變化，利用線性項(xiàng)B0+B1t擬合，考慮降雨對(duì)井水位的即時(shí)及滯后效應(yīng)(張學(xué)民等，1996)，則

昌黎地區(qū)降雨多集中在每年的7～8月，井水位在7月底8月初左右開(kāi)始回升(即時(shí)效應(yīng))，水位在翌年2月底3月初左右達(dá)到最高值(滯后效應(yīng))，因此認(rèn)為，降雨影響即時(shí)效應(yīng)為1個(gè)月，滯后效應(yīng)為5個(gè)月，降雨對(duì)井水位的影響為6個(gè)月(宋曉磊等，2006)，則t1－t0=1，t2－t1=5。根據(jù)式(4)可知降雨對(duì)昌黎井水位的影響有 8 個(gè)待定系數(shù)，即B0、B1、A0、A1、A2、A3、A4、A5。將昌黎井水位觀測(cè)值及降雨量作為時(shí)間的函數(shù)代入(4)式，通過(guò)多元回歸得到8個(gè)系數(shù)，最后得出降雨對(duì)井水位影響的定量值。

經(jīng)褶積濾波消除降雨對(duì)昌黎井水位的影響之后，井水位變化較平穩(wěn)(圖5)，因此，認(rèn)為2007～2008年井水位年變異常是由于2006～2008年的降雨量變化引起的。

圖5 昌黎井原始水位與降雨量校正對(duì)比曲線

3.2 水位破年變與井孔條件改變的關(guān)系

2002年昌黎井的水位處于觀測(cè)以來(lái)的最低值，存在斷流危險(xiǎn)，于2002年10月對(duì)井口裝置進(jìn)行改造，在出水口加裝了三通，三通的內(nèi)徑是原排水管內(nèi)徑的三分之一，泄流量減少，井水位勢(shì)必會(huì)上升。2003年6月出水口流出黑色、有臭味的粘稠狀物質(zhì)，7月15日～8月l日斷流，雨季過(guò)后又緩慢回升至自流，2004年初，出水口底部由于多年腐蝕，出現(xiàn)嚴(yán)重漏水，于2004年7月3～24日采用空壓機(jī)負(fù)壓洗井(岳秀俠等，2005)，洗掉殘留在井壁及井底的沉淀物，使井壁周邊的含水層導(dǎo)水能力變大，井-含水層系統(tǒng)水量交換增加，因此，水位上升(圖6)。由此可見(jiàn)，昌黎井水位2002年10月～2004年的水位年動(dòng)態(tài)消失的破年變異常是由于井孔條件改變引起的(張素欣等，2003)。

圖6 1983～2007年昌黎井水位日均值曲線

3.3 水位破年變與附近井孔抽水關(guān)系

圖7 昌黎井水位、年降雨量對(duì)比曲線

應(yīng)用地下水均衡理論(薛禹群等，1979)，對(duì)一個(gè)含水層(組)來(lái)說(shuō)，任一時(shí)間段內(nèi)，含水層儲(chǔ)水量的變化為地下水補(bǔ)給量與排泄量之差，以水位的變化表現(xiàn)出來(lái)，含水層水位的變化ΔH以觀測(cè)井孔水位的變化反映出來(lái)(曹新來(lái)等，2001)

式中，ΔH為含水層水位變化幅度，Q補(bǔ)為地下水補(bǔ)給量，Q排為抽水量，μ為儲(chǔ)水系數(shù)，A為含水層面積。

從昌黎井水位年動(dòng)態(tài)變化來(lái)看，春灌會(huì)引起水位降低，雨季降雨使水位回升，說(shuō)明降雨是該井區(qū)水的主要補(bǔ)給方式，春灌是主要的排泄方式。2009年昌黎地區(qū)年降雨量為580.4mm，接近多年均值降雨量，即2009年降雨對(duì)水的補(bǔ)給量Q補(bǔ)較往年變化不大，雨季之后不是灌溉季節(jié)，在沒(méi)有劇烈地質(zhì)條件變化時(shí)期，含水層面積A和含水層儲(chǔ)水系數(shù)μ基本不變，抽水Q排增加，勢(shì)必引起昌黎井水位的下降。

2013年10月底昌黎井水位下降，且下降速率較快(圖7)，11月底水位整點(diǎn)值曲線呈鋸齒狀周期性波動(dòng)變化(圖8中黑色圈標(biāo)出)，這種變化多出現(xiàn)在白天，周六、周日的觀測(cè)曲線較為平滑。昌黎井位于昌黎縣酒廠內(nèi)，經(jīng)調(diào)查核實(shí)，距離觀測(cè)井大約30m處有一口抽水井，酒廠周一到周五工作時(shí)間抽水，周末休息停止抽水，與昌黎井水位鋸齒狀變化時(shí)間一致，由于種種原因，抽水量無(wú)法獲得，認(rèn)為臨井抽水引起了昌黎井水位下降。圖7顯示昌黎井水位下降時(shí)段仍具有明顯的固體潮效應(yīng)，說(shuō)明水位下降，未掩蓋其固體潮的變化，潮差變化約9cm，表明該井仍能靈敏地記錄到含水層微小的變化。所以昌黎井水位2013年10月底的下降異常雖因抽水影響，但也不能排除包含地震前兆信息的可能。

圖8 2013年12月1～10日昌黎井整點(diǎn)值水位變化曲線

3.4 水位破年變異常映震能力分析

2010年2月下旬，昌黎井水位在趨勢(shì)下降的背景下，轉(zhuǎn)折上升(圖7)，附近溫泉熱水井開(kāi)發(fā)仍在持續(xù)，而且雨季未到來(lái)且是春灌開(kāi)采時(shí)段，水位出現(xiàn)上升的異常變化，認(rèn)為可能與2010年3月6日灤縣4.2級(jí)地震、4月9日豐南4.1級(jí)地震有關(guān)，5月初水位恢復(fù)正常動(dòng)態(tài)。

2011年3月11日的日本9.0級(jí)地震發(fā)生時(shí)，昌黎井水位有同震響應(yīng)，階升幅度約0.063m，地震之后水位緩慢下降，恢復(fù)正常動(dòng)態(tài)。

2012年3月底至4月底，昌黎井水位在正常下降的背景下，發(fā)生上升的破年變異常(圖9)，該時(shí)段未出現(xiàn)大幅降雨，異常結(jié)束后約1個(gè)月，發(fā)生2012年5月28日唐山4.8級(jí)地震，因此，認(rèn)為這種異常變化可能與此次地震有關(guān)。

圖9 地震前后昌黎井水位的觀測(cè)曲線

上述幾次地震前，昌黎井水位破年變異常特征見(jiàn)表2。

4 結(jié)論與討論

(1)利用動(dòng)態(tài)年周期方法分析2000～2013年昌黎井水位破年變異常，結(jié)果顯示，2001、2002、2007、2008年的破年變異常是因?yàn)榻涤暧绊懀?002年10月～2004年破年變異常是因井孔條件改變引起的，2009年異常為臨近井孔抽水引起，2010、2012年的破年變異常為地震前兆異常，2013年破年變異常有臨井抽水影響，但也不排除包含地震前兆異常信息的可能。

表2 地震前昌黎井水位破年變異常特征表

(2)昌黎井水位自觀測(cè)以來(lái)資料連續(xù)性較好，有一定的映震能力，2010、2012年昌黎井水位在正常年動(dòng)態(tài)變化的基礎(chǔ)上出現(xiàn)短期反向破年變的異常變化，對(duì)應(yīng)了唐山地區(qū)的幾次中等地震。但由于近幾年附近井孔抽水的影響，昌黎井水位日變形態(tài)出現(xiàn)周期性波動(dòng)，而且對(duì)年變動(dòng)態(tài)有一定的影響，在以后的資料分析中應(yīng)考慮這些環(huán)境干擾因素對(duì)井水位的影響。

(3)降雨量的變化會(huì)引起昌黎井水位的年變幅、上升速率改變等破年變特征，因此，對(duì)于類似于昌黎井這種具有良好降雨補(bǔ)給條件的井孔，出現(xiàn)年變規(guī)律改變或者出現(xiàn)短期破年變動(dòng)態(tài)時(shí)，首先要分析降雨對(duì)其觀測(cè)數(shù)據(jù)的影響，然后再考慮其它因素。