王　靜，王利兵，宋　昭，解　真，羅　娜，趙志遠，李細順

(河北省地震局紅山基準地震臺，河北　邢臺　054000)

?

辛集冀21井數(shù)字化水位觀測內(nèi)在質(zhì)量分析

王靜，王利兵，宋昭，解真，羅娜，趙志遠，李細順

(河北省地震局紅山基準地震臺，河北邢臺054000)

對辛集冀21井模擬、數(shù)字化水位觀測精度變化進行對比研究，結(jié)果表明數(shù)字化水位資料觀測精度高、儀器運行穩(wěn)定。文章分析并總結(jié)了影響數(shù)字化水位觀測精度的因素，為今后觀測資料分析及異常判定提供依據(jù)，也為進一步提高水位觀測內(nèi)在質(zhì)量奠定基礎(chǔ)。

水位觀測；觀測精度；內(nèi)在質(zhì)量；影響因素

0　引言

地下水位能夠隨潮汐引力的大小發(fā)生升降變化，表現(xiàn)較明顯的地球固體潮效應(yīng)。辛集冀21井觀測的水位數(shù)據(jù)日變動態(tài)均呈現(xiàn)顯著規(guī)律性，較多地呈雙峰雙谷型，分析井點對固體潮記錄的響應(yīng)能力可以判別該井點的靈敏度。

固體潮具有豐富的潮汐分量，由月球引力引起的O1波和M2波受氣壓影響的程度較小，理論上都是分析水位固體潮較理想的指標。但由于O1波幅度較小，信息往往隱含在噪聲中難以識別，而M2波幅度為水位潮汐信息中最大的一個諧波，因此M2波是分析井水位固體潮以及含水層特征最有效的潮汐分量[1]。通過對不同時間長度采樣數(shù)據(jù)的頻譜分析表明，M2波幅度與潮汐長度變化密切相關(guān)，選取不同時間長度的采樣數(shù)據(jù)對M2波幅度的差異性較大。一個月時間尺度的采樣數(shù)據(jù)包含2個大潮和小潮，且有足夠的數(shù)據(jù)長度可以使計算結(jié)果相對穩(wěn)定[2]。

本文選取2015年辛集冀21井模擬和數(shù)字化地下水位數(shù)據(jù)，采用維尼迪可夫(Venedlkov)調(diào)和分析方法，逐月分析計算水位對固體潮M2波的響應(yīng)能力，以此檢驗井水位數(shù)據(jù)的內(nèi)在質(zhì)量。

1　觀測井概況

辛集冀21井觀測站位于河北辛集市北偏西7 km處，屬邢臺—河間地震帶。井深2 052.4 m，觀測層為震旦亞界白云巖，埋深1 539.40 m至2 052.37 m，厚度512.97 m。站點的地質(zhì)—水文條件清楚，套管止水措施效果良好且不受降雨直接滲入補給影響，但受地下熱水開采影響明顯。目前，進行正式觀測的測項包括模擬水位、數(shù)字化水位和氣象三要素。觀測儀器采用SW-40型模擬水位儀、LN-3A型數(shù)字式水位儀及WYY-1型氣象三要素觀測儀。

辛集冀21井模擬觀測采用SW-40型日記水位計，圖紙記錄方式，人工讀取整點值數(shù)據(jù)。數(shù)字化觀測采用LN-3A型數(shù)字水位儀，使用中國地震前兆臺網(wǎng)數(shù)據(jù)管理系統(tǒng)收取儲存數(shù)據(jù)，用中國地震前兆臺網(wǎng)數(shù)據(jù)處理系統(tǒng)對數(shù)據(jù)進行預處理。該儀器采用高精度壓力式水位傳感器，無機械隨動系統(tǒng)，觀測精度不受水位埋深和井斜影響，有極高的水位跟蹤速度。其工作原理是利用測量液體內(nèi)部壓強的方法來測量水位。水位傳感器的量程范圍為0～10 m，水位分辨率為1 mm[3]。

2　辛集冀21井模擬、數(shù)字化水位觀測資料內(nèi)在質(zhì)量對比分析

通過對潮汐效應(yīng)的分析計算可獲得水位動態(tài)及潮汐響應(yīng)函數(shù)的基本特征參數(shù)，即為觀測井的特征參數(shù)。其可以客觀反映觀測井的基本“素質(zhì)”和周邊環(huán)境的影響，可用于定量評價各個井的內(nèi)在質(zhì)量優(yōu)劣。水位數(shù)據(jù)的觀測精度為中誤差與潮汐因子比值，固體潮相對誤差越小，數(shù)字化水位觀測越穩(wěn)定。

水位觀測內(nèi)在質(zhì)量是全國觀測資料的質(zhì)量單項評比與綜合評比中的一項，評比分數(shù)為20分。根據(jù)2015年新修訂的水位觀測資料評分標準，水位數(shù)據(jù)內(nèi)在質(zhì)量占20分：當觀測精度<0.01得15分，每增大 0.01扣0.35分；當觀測精度<0.04，潮汐因子≥1.00(動水位≥0.08)時，得5分，每降低0.1(動水位0.05)扣0.05分；觀測精度>0.4時，潮汐因子不得分。

采用Venedlkov調(diào)和分析方法對辛集冀21井模擬和數(shù)字化水位觀測資料進行分析。表1為模擬水位2012至2015年(2011年無數(shù)據(jù))與數(shù)字水位2011至2015年潮汐因子、中誤差、觀測精度變化年均值對比表，表2為2015年1至12月模擬、數(shù)字化水位觀測資料潮汐因子、中誤差、觀測精度變化月均值對比表。

表1　辛集冀21井模擬、數(shù)字化水位2011年至2015年潮汐因子、觀測精度變化年均值對比表

表2　辛集冀21井模擬、數(shù)字化水位2015年1月至12月潮汐因子、觀測精度變化月均值對比表

由表2可以看出，辛集冀21井模擬水位觀測資料12個月的潮汐因子在2.61～2.98之間，均大于1.00，不扣分；觀測精度9月為0.007 8(<0.01)，不扣分，其余每月均超過0.01。圖1為辛集冀21井模擬、數(shù)字化水位觀測資料2015年1至12月觀測精度曲線對比圖，經(jīng)分析發(fā)現(xiàn)，模擬水位觀測資料的觀測精度低于數(shù)字化的，且后者相對穩(wěn)定。數(shù)字化水位觀測資料潮汐因子在3.07～3.14之間，觀測精度在0.006～0.012之間。由表1、表2可知，辛集冀21井不論是模擬還是數(shù)字化水位觀測資料，固體潮效應(yīng)都很明顯，較好地反映了地球潮汐變化。由于模擬水位觀測過程中，需每24小時更換一張記錄紙，數(shù)據(jù)靠人工讀取，存在換紙造成的資料缺數(shù)和人為誤差及記錄筆污染圖紙和記錄曲線虛的弊端[4]。因此，數(shù)字化觀測數(shù)據(jù)的固體潮效應(yīng)反映都優(yōu)于相應(yīng)的模擬觀測數(shù)據(jù)，數(shù)字化觀測數(shù)據(jù)的觀測精度也優(yōu)于相應(yīng)的模擬觀測數(shù)據(jù)。

圖1　辛集冀21井數(shù)字化、模擬水位觀測資料2015年1月至12月觀測精度曲線圖Fig.1　Observation accuracy of digitized water level and simulated water level in Xinjiji Well 21 from January to December in 2015

3　影響數(shù)字化水位觀測內(nèi)在質(zhì)量的因素

3.1周圍井孔地下水開采造成干擾

辛集冀21井水位觀測主要干擾源為每年冬季采暖期周圍井孔的地熱開采。采暖期開始就會引起地下水趨勢性下降，采暖期結(jié)束后開始緩慢上升。井水位的年起伏完全受控于周圍井孔的地熱開采，影響水位的正常動態(tài)。如圖2所示，水位自2011年至2015年已累計下降約26.5 m。

圖2　辛集冀21井2011年至2015年水位日均值圖Fig.2　Daily mean values of water level in Xinjiji Well21 from 2011 to 2015

3.2儀器零漂

據(jù)張道口-1井LN-3A型水位儀2001年6月17日始測至2003年1月8日進行零點檢查，1年半時間零漂達0.2 m，日零漂幅度0.35 mm[5]。據(jù)寧德一號井數(shù)字化水位儀零漂問題的跟蹤校測與分析結(jié)果看，數(shù)字化水位儀在觀測過程中存在明顯的零漂現(xiàn)象，水位零漂幅度受儀器室溫度、井水溫、傳感器埋深多種因素綜合作用影響[6]。由此可知，儀器的零漂降低了數(shù)據(jù)的可靠性，影響數(shù)據(jù)的觀測精度。目前，大部分臺站缺乏零點檢查這方面工作的實施記錄，建議每半年對LN-3A型水位儀進行水位探頭零點的歸零處理，并做記錄。

3.3探頭位置影響觀測精度

LN-3A型水位儀工作原理是利用測量液體內(nèi)部壓強的方法來測量水位。液體內(nèi)部的壓強與其深度的函數(shù)關(guān)系式為P=ρgh，這里密度ρ和重力g都是常數(shù)，壓強P與深度h成正比關(guān)系。當把壓力傳感器固定在水下某一深度，水位變化時壓力傳感器的輸出信號產(chǎn)生與水位同步的變化[7]。分析認為由于儀器自身的原因?qū)е聣簭娕c深度之間的線性關(guān)系不穩(wěn)定。

LN-3A型水位儀在水中的深度與水位觀測精度成反比關(guān)系：在水中越深，觀測精度越差；在水中越淺，觀測精度越好，但目前仍無法給出兩者的定量關(guān)系。因此，放置探頭時應(yīng)盡量靠近水面。

4　結(jié)語

(1) 辛集冀21井水位觀測儀器工作狀態(tài)穩(wěn)定且數(shù)據(jù)連續(xù)可靠，具有一定的研究價值。

(2) 同一井孔模擬水位觀測與數(shù)字化水位觀測相比較，數(shù)字化觀測可以提高水位觀測的內(nèi)在質(zhì)量，且具有較高的穩(wěn)定性和觀測精度。

(3) 辛集冀21井水位潮汐效應(yīng)明顯，水位潮汐因子和觀測精度分別為3.103 3和0.008 2。

(4) 針對LN-3A型水位儀的零點漂移問題，在每半年對水位儀現(xiàn)場檢查時應(yīng)進行零點歸零處理。

[1]尤宇星，陳小云，謝文杰，等.泉州地區(qū)地下水位觀測資料初步分析[J].華北地震,2014(1):65-71.

[2]陳瑩，黃輔瓊，朱石軍，等.采樣數(shù)據(jù)長度對潮汐因子計算精度的影響[J].地震地磁觀測與研究，2012，33(5/6):146-151.

[3]中國地震局地震預測研究所.LN-3A型數(shù)字水位儀使用說明書[N].中國地震局地震預測研究所，2006:7.

[4]解曉靜.瓊海加積井水位模擬觀測與數(shù)字化觀測資料對比分析[J].地震地磁觀測與研究，2010，31(3):83-87.

[5]王建國，劉春國，陳華靜，等.影響張道口—1井數(shù)字化水位內(nèi)在質(zhì)量的因素分析[J].內(nèi)陸地震，2007(2):155-159.

[6]蘭陵，李新勇，夏愛國.新04井數(shù)字化水位觀測與模擬觀測資料的對比分析[J].內(nèi)陸地震，2004，18(4):382-384.

[7]楊鼎鴻，程慶斌，鄧聰，等.數(shù)字化水位觀測儀零漂分析[J].地震地磁觀測與研究，2013，34(3/4):144-149.

Intrinsic Quality of Observation of Digitized Water Level in Xinjiji Well 21

WANG Jing, WANG Li-bing, SONG Zhao, XIE Zhen, LUO Na, ZHAO Zhi-yuan, LI Xi-shun

(Hongshan Reference Seismological Station of Earthquake Administration of Hebei Province, Xingtai,Hebei 054000,China)

The accuracy of simulated water level observation and digital water level observation in Xinjiji well 21 is compared. The results show that the observation accuracy of digital water level data is high, and the operation of the instrument is stable. Several factors affecting the accuracy of digital water level observation are analyzed and summarized. This summary provides the basis for the analysis of observation data and the abnormity judgment in the future, and for further improving the intrinsic quality of water level observation.

Water level observation; Observation accuracy; Intrinsic quality; Influence factor

1000-6265(2016)03-0018-03

2016-03-21

河北省地震局監(jiān)測、預測、科研三結(jié)合課題(DZ201504024092)。

王靜(1983—)，女，河北省邢臺人。2009年畢業(yè)于河北工程大學，工程師。

P332.3

A