胡玉良 程冬焱 李惠玲 穆慧敏 郭宇 劉俊芳

1）山西省地震局，山西省太原市舊晉祠路二段69號(hào) 030021

2）太原大陸裂谷動(dòng)力學(xué)國(guó)家野外科學(xué)觀測(cè)研究站，太原 030025

3）山西省地震局定襄地震臺(tái)，山西定襄 035400

0 引言

國(guó)內(nèi)外學(xué)者以大量的觀測(cè)事實(shí)論證了水位動(dòng)態(tài)不僅可以直接反映含水層應(yīng)力應(yīng)變的變化，還可反映井-含水層水動(dòng)力條件的變化（國(guó)家地震局科技監(jiān)測(cè)司，1994），因此水位動(dòng)態(tài)研究是研究地球某些動(dòng)力過(guò)程和監(jiān)測(cè)地震孕育的重要途徑（汪成民等，1988）。

隨著 “九五”、“十五”、“十一五”項(xiàng)目的全面建設(shè)，流體數(shù)字化觀測(cè)儀器在全國(guó)開(kāi)始投入運(yùn)行，為水位動(dòng)態(tài)研究提供了豐富的數(shù)據(jù)（廖麗霞等，2013），以此為基礎(chǔ)的地震觀測(cè)技術(shù)由模擬觀測(cè)向數(shù)字化觀測(cè)發(fā)展（耿杰，2008）?！笆濉逼陂g，山西省地震局對(duì)鎮(zhèn)川、朔州、太原、靜樂(lè)、祁縣、介休、孝義、漫水和東郭流體井進(jìn)行了技術(shù)改造，統(tǒng)一安裝了由中國(guó)地震局地震預(yù)測(cè)研究所研制的LN-3A水位儀?！笆晃濉逼陂g，山西省地震局對(duì)榆社紅崖頭流體井和忻州鴉兒坑流體井進(jìn)行了技術(shù)改造，統(tǒng)一安裝了由中國(guó)地震局地殼應(yīng)力研究所研制的SWY-Ⅱ水位儀。

“九五”、“十五”項(xiàng)目建設(shè)完成后，耿杰等（2001）、黃春玲等（2010）曾做過(guò)數(shù)字化水位觀測(cè)資料與模擬觀測(cè)資料的對(duì)比分析，主要針對(duì)LN-3型數(shù)字水位儀和SW 40-1日記式水位儀或LN-3A型數(shù)字水位儀（以下簡(jiǎn)稱(chēng) LN-3A水位儀）和 SW 40-1日記式水位儀（以下簡(jiǎn)稱(chēng)SW 40-1水位儀）。由于SWY-Ⅱ型數(shù)字化水位儀（以下簡(jiǎn)稱(chēng)SWY-Ⅱ水位儀）為新研制的儀器（何案華等，2012），尚未對(duì)SWY-Ⅱ水位儀和LN-3A水位儀進(jìn)行對(duì)比分析研究，因此，本文通過(guò)在靜樂(lè)井對(duì)新安裝的1套SWY-Ⅱ及原有的LN-3A水位儀開(kāi)展為期3個(gè)多月的并行觀測(cè)，從穩(wěn)定性、數(shù)據(jù)一致性、數(shù)據(jù)精度、系統(tǒng)漂移、水震波、氣壓效應(yīng)、自然環(huán)境干擾這7方面進(jìn)行對(duì)比分析。

1 對(duì)比分析

1.1 儀器運(yùn)行

2015年4月17日將 SWY-Ⅱ水位儀安裝到靜樂(lè)井，從4月18日開(kāi)始投入運(yùn)行，臺(tái)站值班人員按照流體學(xué)科規(guī)范要求做好儀器維護(hù)、數(shù)據(jù)預(yù)處理、觀測(cè)日志填寫(xiě)、儀器校測(cè)等工作，確保SWY-Ⅱ水位儀產(chǎn)出數(shù)據(jù)連續(xù)、穩(wěn)定可靠，同時(shí)對(duì)靜樂(lè)井原有的LN-3A水位儀的觀測(cè)數(shù)據(jù)進(jìn)行跟蹤分析。利用中國(guó)地震前兆臺(tái)網(wǎng)運(yùn)行評(píng)價(jià)系統(tǒng)對(duì)兩套儀器的數(shù)據(jù)完整率進(jìn)行統(tǒng)計(jì)（表1）。

由表1可知，兩套水位儀運(yùn)行穩(wěn)定，均未發(fā)生故障，數(shù)據(jù)連續(xù)率均達(dá)99.99%以上。

表1 靜樂(lè)井?dāng)?shù)字水位儀觀測(cè)數(shù)據(jù)完整率對(duì)比

1.2 數(shù)據(jù)一致性分析

首先對(duì)靜樂(lè)井SWY-Ⅱ、LN-3A兩套水位儀曲線(xiàn)形態(tài)進(jìn)行定性分析，提取2015年4月18日～7月31日的分鐘值數(shù)據(jù)，如圖1所示，從曲線(xiàn)形態(tài)上看，兩套儀器的趨勢(shì)是一致的，從圖上也能明顯看出在水震波記錄的幅度上有差異（見(jiàn)本文1.6節(jié)），為了更詳細(xì)地說(shuō)明兩套儀器記錄數(shù)據(jù)的差異是否隨時(shí)間而變化，將兩套儀器產(chǎn)生的分鐘值數(shù)據(jù)相減，得到的差值如圖2所示。兩套儀器除了在地震波記錄的幅度上有差異外，在其余時(shí)間內(nèi)差值最小值是0.091m，最大值是0.115m，平均差值為0.103m，差值變化幅度不超過(guò)0.024m。平均差值不在零附近主要是由于兩套儀器之間存在系統(tǒng)誤差，通過(guò)修正SWY-Ⅱ或LN-3A的參數(shù)值即可解決，差值曲線(xiàn)形態(tài)和變化幅度說(shuō)明兩套儀器的一致性非常好。

圖1 靜樂(lè)井2015年4月18日～7月31日數(shù)字水位儀分鐘值數(shù)據(jù)對(duì)比曲線(xiàn)圖

圖2 靜樂(lè)井2015年4月18日～7月31日數(shù)字水位儀分鐘值數(shù)據(jù)差值曲線(xiàn)圖

標(biāo)準(zhǔn)偏差越小，表明測(cè)量的分散度越小，測(cè)量的精度越高，用標(biāo)準(zhǔn)偏差可以表征測(cè)量精度。在無(wú)系統(tǒng)誤差的情況下，一般把標(biāo)準(zhǔn)偏差作為測(cè)量誤差（何平，1994）。由于兩套儀器均放置在水井中，當(dāng)井水位發(fā)生變化時(shí)兩套儀器的產(chǎn)出數(shù)據(jù)也會(huì)隨之發(fā)生變化。為了計(jì)算兩套儀器的標(biāo)準(zhǔn)偏差，從而確定兩套儀器的精度，需要按照一定的條件對(duì)兩套儀器產(chǎn)出的分鐘值數(shù)據(jù)進(jìn)行篩選，條件為：①選取的時(shí)間段必須短；②確保兩套儀器觀測(cè)數(shù)據(jù)要平穩(wěn)，波動(dòng)要盡可能地小。基于以上兩點(diǎn)，選取了2015年4月23日7：00～8：00的分鐘值，使用 SPSS軟件做單因素方差分析（羅納德，2010；里斯，2011），計(jì)算了兩套儀器的標(biāo)準(zhǔn)偏差。

為了使用 SPSS軟件對(duì)兩套儀器時(shí)均值數(shù)據(jù)做配對(duì)樣本 t檢驗(yàn)（羅納德，2010；里斯，2011），需要消除兩套儀器之間存在的系統(tǒng)誤差。本文對(duì)LN-3A水位儀時(shí)均值進(jìn)行了修正，具體修正過(guò)程如下：SWY-Ⅱ、LN-3A水位儀 2015年 4月 18日零點(diǎn)數(shù)據(jù)分別是 9.693和9.787，兩者相差0.094，故對(duì)LN-3A水位儀從2015年4月18日～7月31日之間所有時(shí)均值數(shù)據(jù)都減去0.094即可。

1.2.1 單因素方差分析

方差分析就是分析不同水平下各個(gè)總體的均值是否有顯著的差異。統(tǒng)計(jì)推斷方法是計(jì)算F統(tǒng)計(jì)量，進(jìn)行F檢驗(yàn)。單因素方差分析指的是自變量只有一個(gè)的方差分析，計(jì)算公式為（邱軼兵，2008）

其中，VA為因素偏差平方和；Ve為誤差平方和；FA為一個(gè)統(tǒng)計(jì)量，服從自由度為（FA，fe）的F分布，其中FA為因素偏差平方和自由度；fe為誤差平方和自由度。給定α＝0.05，如果FA≥F0．05（fA，fe），就說(shuō)明因素A變化的影響大于誤差的影響，即該因素變化顯著。

利用SPSS軟件對(duì) SWY-Ⅱ水位儀和 LN-3A水位儀2015年4月23日7：00～8：00之間的分鐘值數(shù)據(jù)進(jìn)行單因素方差計(jì)算，計(jì)算結(jié)果如表2。

表2 描述統(tǒng)計(jì)

由表2可知，LN-3A水位儀標(biāo)準(zhǔn)偏差是0.000279，SWY-Ⅱ水位儀標(biāo)準(zhǔn)偏差是0.00022，說(shuō)明兩套儀器測(cè)量精度基本無(wú)差異，并且均優(yōu)于0.001m。

1.2.2 配對(duì)t檢驗(yàn)基本步驟

首先，提出兩配對(duì)樣本t檢驗(yàn)的假設(shè)為 H0，兩樣本均值 μ1、μ2無(wú)顯著差異，即 H0＝μ1-μ2＝0。

利用SPSS軟件對(duì)SWY-Ⅱ水位儀和LN-3A水位儀2015年4月18日～7月31日的時(shí)均值數(shù)據(jù)進(jìn)行配對(duì)t檢驗(yàn)，計(jì)算結(jié)果如表3。

表3 配對(duì)t檢驗(yàn)結(jié)果

從表3可以得到SWY-Ⅱ水位儀和LN-3A水位儀數(shù)據(jù)之差的均值、標(biāo)準(zhǔn)差、均值標(biāo)準(zhǔn)誤差、95%的置信區(qū)間以及t檢驗(yàn)值、自由度和雙側(cè)概率值。通過(guò)分析可知，雙側(cè)概率值小于0.05，說(shuō)明兩套儀器均值存在差異；兩套儀器均值相差0.009m，標(biāo)準(zhǔn)偏差為0.005m，說(shuō)明兩套儀器均值雖有差異，但均在觀測(cè)允許的誤差范圍內(nèi)，因此，從觀測(cè)的角度來(lái)說(shuō)，可認(rèn)為兩套儀器均值是一致的。

1.3 數(shù)據(jù)內(nèi)在質(zhì)量

井水位潮汐是固體地球潮汐的一種次生效應(yīng)（汪成民，1988），由于觀測(cè)井一般深度只有數(shù)百米或數(shù)千米，處于地球表層，因此地球表層潮汐應(yīng)力分布會(huì)影響靜水位反周期性變化（國(guó)家地震局科技監(jiān)測(cè)司，1995）。呂芳等（2012）認(rèn)為靜樂(lè)井的井-含水層系統(tǒng)對(duì)地殼應(yīng)力-應(yīng)變有較強(qiáng)的響應(yīng)能力，本文通過(guò)對(duì)靜樂(lè)井兩套數(shù)字水位儀分鐘值數(shù)據(jù)和理論固體潮進(jìn)行曲線(xiàn)對(duì)比及計(jì)算M2波潮汐因子及其觀測(cè)精度，也印證了這一點(diǎn)。

從圖3可以看出，兩套水位儀均能記錄到完整的日波與半日波，相位同步清晰，幅度較大，與理論固體潮曲線(xiàn)相位一致。

圖3 2015年5月1日SWY-Ⅱ、LN-3A水位及理論固體潮時(shí)均值數(shù)據(jù)曲線(xiàn)

使用兩套水位儀分鐘值數(shù)據(jù)，利用中國(guó)地震前兆臺(tái)網(wǎng)運(yùn)行評(píng)價(jià)系統(tǒng) 2011版軟件計(jì)算兩套水位儀M2波潮汐因子及其觀測(cè)精度（劉春國(guó)等，2015），并進(jìn)行對(duì)比分析，其結(jié)果如表4所示。

從表4可以看出，對(duì)于M 2波潮汐因子和觀測(cè)精度，SWY-Ⅱ水位儀均優(yōu)于LN-3A水位儀。

表4 水位儀M 2波潮汐因子及其觀測(cè)精度對(duì)比

1.4 系統(tǒng)漂移

SWY-Ⅱ和LN-3A水位儀傳感器均采用擴(kuò)散硅半導(dǎo)體壓力器件。一方面半導(dǎo)體的溫度特性使硅壓力傳感器的零點(diǎn)和靈敏度隨溫度而發(fā)生漂移（孫鳳玲等，2007；張艷鋒等，2008），導(dǎo)致測(cè)量精度會(huì)下降（何平等，2008）；另一方面，由于傳感器埋深（即傳感器探頭至水面的距離）會(huì)隨著時(shí)間發(fā)生變化，也會(huì)導(dǎo)致測(cè)量精度發(fā)生變化（楊鼎鴻等，2013），因此有必要計(jì)算兩套數(shù)字化水位儀在工作環(huán)境下的系統(tǒng)漂移。

由于每天均對(duì)靜樂(lè)井SW 40-1水位儀進(jìn)行校測(cè)，因此可以認(rèn)為SW 40-1水位儀測(cè)量值是準(zhǔn)確的，故以SW 40-1水位儀2015年4月18日～6月30日的日均值為基準(zhǔn)，與兩套數(shù)字化水位儀同時(shí)段的日均值分別相減，得到兩組差值，以此來(lái)確定每套儀器隨時(shí)間變化漂移情況。差值對(duì)比曲線(xiàn)如圖4所示，SW 40-1水位儀日均值曲線(xiàn)如圖5所示。

圖4 靜樂(lè)井兩套數(shù)字水位儀系統(tǒng)誤差對(duì)比分析曲線(xiàn)

由圖4可知，SWY-Ⅱ水位儀正向最大漂移量為0.019m，LN-3A水位儀為0.011m；SWY-Ⅱ水位儀反向最大漂移量為-0.011m，LN-3A水位儀為-0.013m；SWY-Ⅱ水位儀漂移量平均值為0.0014m，LN-3A為-0.0048m。從曲線(xiàn)形態(tài)上來(lái)看，2條曲線(xiàn)較為一致，只是幅度略有差異；對(duì)2條曲線(xiàn)做一階差分后可以看出，SWY-Ⅱ水位儀的漂移量稍大于 LN-3A，最大相差0.009m，但漂移量都在mm級(jí)，這可能是由于兩套水位儀壓力傳感器的特性差異造成的?？傮w上來(lái)說(shuō)，兩套水位儀測(cè)量值均會(huì)隨時(shí)間發(fā)生不同程度漂移，并且漂移量隨時(shí)間發(fā)生非線(xiàn)性變化，因此兩套數(shù)字化水位儀均需按照規(guī)范要求定期進(jìn)行校測(cè)。

對(duì)比圖4和5可知，LN-3A和SWY-Ⅱ兩套水位儀系統(tǒng)漂移量最大時(shí)，靜樂(lè)井實(shí)際水位觀測(cè)值最小，由此說(shuō)明水位的上下波動(dòng)對(duì)兩套水位探頭的壓力響應(yīng)系數(shù)有一定的影響。

圖5 靜樂(lè)井SW 40-1水位曲線(xiàn)

1.5 氣壓效應(yīng)

氣壓對(duì)地下水影響具有普遍性（汪成民等，1988；耿杰等，2001），氣壓每變化1hPa所引起的水位變化稱(chēng)為井水位氣壓效率，單位為mm／hPa，利用SPSS軟件對(duì)靜樂(lè)井兩套水位儀日均值與氣壓數(shù)據(jù)進(jìn)行皮爾遜相關(guān)系數(shù)計(jì)算后發(fā)現(xiàn)，SWY-Ⅱ水位儀與氣壓相關(guān)系數(shù)為0.307，LN-3A水位儀與氣壓相關(guān)系數(shù)為0.293，相關(guān)程度均較弱，兩套水位儀測(cè)值與氣壓的相關(guān)系數(shù)較為接近。

1.6 水震波效應(yīng)

2015年4月～7月，靜樂(lè)井SWY-Ⅱ和 LN-3A水位儀記錄到6級(jí)以上地震4次（表5）。

表5 靜樂(lè)井?dāng)?shù)字水位儀記錄地震情況對(duì)照表

地下水位水震波受觀測(cè)層的巖性、埋深、井孔-含水層間的導(dǎo)水能力、井孔水柱高度、地震地點(diǎn)、震級(jí)（耿杰等，2001、2008）及觀測(cè)儀器（廖麗霞等，2009；廖麗霞等，2013）等諸多因素的影響。分析結(jié)果表明，LN-3A和SWY-Ⅱ水位儀在分鐘值采樣下均能清晰記錄到水震波，但震蕩時(shí)間和震蕩幅度均有明顯差別，具體原因是由于SWY-Ⅱ水位儀內(nèi)部采用秒采樣，然后將60s數(shù)據(jù)進(jìn)行平均后獲得分鐘值數(shù)據(jù)，這樣相當(dāng)于SWY-Ⅱ水位儀進(jìn)行了平均濾波，信號(hào)平滑度較好，但降低了靈敏度；而LN-3A水位儀是在每個(gè)分鐘某一個(gè)時(shí)刻進(jìn)行數(shù)據(jù)采集獲得分鐘值數(shù)據(jù)。

由于SWY-Ⅱ水位儀保存了秒采樣的水位數(shù)據(jù)，這里截取了SWY-Ⅱ水位儀2015年4月25日尼泊爾地震前后記錄的秒采樣數(shù)據(jù)曲線(xiàn)，如圖6所示，通過(guò)圖6可以看到清晰的水震波，震蕩時(shí)間為50m in，和LN-3A水位儀相差無(wú)幾；震蕩幅度達(dá)到1.757m，比LN-3A水位儀高0.895m。基于此，建議將SWY-Ⅱ水位儀秒采樣數(shù)據(jù)接入到前兆數(shù)據(jù)管理系統(tǒng)，這樣獲取的信息更多，有利于進(jìn)一步研究。

圖6 靜樂(lè)井SWY-Ⅱ水位儀2015年4月25日秒數(shù)據(jù)曲線(xiàn)

1.7 自然環(huán)境干擾

利用前兆臺(tái)網(wǎng)（臺(tái)站）觀測(cè)數(shù)據(jù)跟蹤分析軟件查詢(xún)可知，2015年4月18日～7月31日，LN-3A水位儀記錄到自然環(huán)境干擾事件5次，其中降水影響4次，雷電影響1次。以5月12日～17日記錄的事件為例進(jìn)行說(shuō)明，通過(guò)與SWY-Ⅱ水位儀同時(shí)段觀測(cè)數(shù)據(jù)對(duì)比可知（圖7（b）和（c）），2條曲線(xiàn)形態(tài)完全一致，從而說(shuō)明了觀測(cè)數(shù)據(jù)變化與觀測(cè)系統(tǒng)無(wú)關(guān)，為了進(jìn)一步說(shuō)明是何種因素導(dǎo)致的水位上升，增加了忻州鴉兒坑雨量曲線(xiàn)，并且與2014年數(shù)據(jù)進(jìn)行同期對(duì)比，如圖7所示。

圖7 靜樂(lè)井水位與降雨量同期對(duì)比曲線(xiàn)

從圖7可以看出，2014年5月9日～11日與2015年5月9日～10日同期均有降雨，導(dǎo)致靜樂(lè)井水位上升。

2 結(jié)論與討論

兩套數(shù)字化水位儀在并行運(yùn)行期間工作穩(wěn)定，數(shù)據(jù)完整；通過(guò)曲線(xiàn)形態(tài)對(duì)比，兩套儀器一致性非常好；采用相關(guān)統(tǒng)計(jì)學(xué)方法得知，兩套數(shù)字化水位儀觀測(cè)數(shù)據(jù)精度均優(yōu)于0.001m，產(chǎn)出均值一致性較好；利用相關(guān)軟件對(duì)兩套數(shù)字化水位儀時(shí)均值數(shù)據(jù)計(jì)算了M 2波潮汐因子及其觀測(cè)精度，結(jié)果表明，SWY-Ⅱ水位儀觀測(cè)精度優(yōu)于LN-3A水位儀；基于靜樂(lè)井模擬水位數(shù)據(jù)，對(duì)兩套數(shù)字化水位儀系統(tǒng)漂移進(jìn)行了計(jì)算，結(jié)果顯示，兩套數(shù)字化水位儀均出現(xiàn)不同程度系統(tǒng)漂移，因此在數(shù)字化水位儀運(yùn)行過(guò)程中，需對(duì)水位傳感器定期進(jìn)行標(biāo)定校測(cè)，確保水位數(shù)據(jù)真實(shí)可靠；兩套水位儀水位數(shù)據(jù)與氣壓相關(guān)性較為一致；分鐘值采樣下，水震波記錄震蕩時(shí)間及水位振幅有明顯差別；自然環(huán)境干擾方面，兩套水位儀干擾程度較為一致。

通過(guò)對(duì)比分析發(fā)現(xiàn)，兩套數(shù)字化水位儀在系統(tǒng)漂移及水震波記錄方面存在著一些差異，SWY-Ⅱ水位儀秒采樣數(shù)據(jù)記錄地震波信息較為豐富，這為下一步深入研究提供了方向。