尤宇星，陳小云，謝文杰，林文峰，吳聲偉，黃圣棕

（1.福建省地震局，福州 350001；2.泉州市地震局，福建 泉州 362000）

0 引言

地下水位變化伴隨著地球內(nèi)部巖石受力變形、破裂，流體遷移的過程，也伴隨著地震孕育和發(fā)生發(fā)展的過程，揭示了許多有用的信息。因此，對地下水位變化的觀測、識別提取地震前兆信息是十分有意義的［1－8］?！熬盼濉?、“十五”和“十一五”期間建設(shè)的泉州地區(qū)數(shù)字化地下流體前兆觀測網(wǎng)自正式觀測以來，為福建省地震預(yù)測提供了豐富的地下流體前兆觀測資料。由于異常是相對于正常動態(tài)而言的，其異常信息量是在與其正常動態(tài)的對比中獲得的。因此，筆者通過具體分析研究，計(jì)算出泉州地區(qū)主要觀測井的正常動態(tài)特征參數(shù)，并對這些觀測井內(nèi)在質(zhì)量進(jìn)行分析評估，以此作為異常對比識別可信度的主要依據(jù)。

1 泉州地區(qū)數(shù)字化水位觀測基本概況

由于福建省內(nèi)80%的強(qiáng)震活動都是分布在長樂－詔安斷裂帶上，該斷裂帶是一條縱貫福建省東南沿海的活動斷裂帶。而作為地震前兆觀測，觀測點(diǎn)必須能夠控制一個地區(qū)的構(gòu)造活動，因?yàn)榛顒訑嗔焉弦恍┯欣臉?gòu)造部位往往是應(yīng)力集中部位，出現(xiàn)前兆異常的概率明顯高于遠(yuǎn)離斷裂處。所以泉州地區(qū)的各個觀測點(diǎn)均分布在NNE－NE 向的長樂－詔安斷裂帶和NW 向的永安－晉江斷裂帶及其次一級斷裂帶上，具體情況見表1。泉州地區(qū)各觀測點(diǎn)的水位觀測均采用LN－3型數(shù)字化水位儀，該儀器由2部分組成，一部分是水位傳感器，另一部分是數(shù)字記錄器，其優(yōu)點(diǎn)是采樣率高、信息量大、自動化程度高、資料傳遞快速等。

2 泉州地區(qū)地下水位資料質(zhì)量分析

資料的完整性、穩(wěn)定性、內(nèi)在質(zhì)量精度是衡量觀測資料質(zhì)量的重要指標(biāo)。觀測資料只有達(dá)到觀測技術(shù)標(biāo)準(zhǔn)要求時，其應(yīng)用價(jià)值才高，才能更好地判斷出觀測資料的正常動態(tài)背景，以利于識別異常。為此有必要對觀測資料進(jìn)行完整性、穩(wěn)定性、內(nèi)在質(zhì)量及其影響因素的檢驗(yàn)［6］。

2.1 觀測資料的完整性分析

利用下述公式計(jì)算統(tǒng)計(jì)了2006—2011年泉州地區(qū)各觀測點(diǎn)自觀測以來水位整點(diǎn)值觀測資料的完整連續(xù)率：

得出泉州地區(qū)各觀測井點(diǎn)除南安2號井、南安梅山井（2010年后停測）、惠安2號井（2011年后停測）外，其它11口觀測井的完整連續(xù)率均大于97.4%，其中2011年均超過99%，影響完整率的主要因素是儀器故障。

2.2 水位觀測資料的內(nèi)在質(zhì)量對比分析

地下水位能夠能隨潮汐引力的大小而發(fā)生升降變化，體現(xiàn)較明顯的地球固體潮效應(yīng)已經(jīng)是眾所皆知的事實(shí)。具體表現(xiàn)在泉州地區(qū)各測點(diǎn)觀測的水位數(shù)據(jù)日變動態(tài)均呈現(xiàn)顯著規(guī)律性，較多地呈雙峰雙谷型，分析各測點(diǎn)對固體潮記錄的響應(yīng)能力可以判別各測點(diǎn)的靈敏度。固體潮具有豐富的潮汐分量，由月球引力引起的O1波和M2波受氣壓影響的程度較小，理論上O1波和M2波對于分析水位固體潮都是較理想的指標(biāo)。但O1波的幅度較小，信息往往掩埋在噪聲中，而M2波幅度為水位潮汐信息中最大的一個諧波，因此M2波是分析井水位固體潮以及含水層特征最有效的潮汐分量。通過對不同時間長度采樣數(shù)據(jù)的頻譜分析表明：M2波幅度與潮汐長度變化密切相關(guān)，選取不同時間長度的采樣數(shù)據(jù)對M2波幅度的差異性較大，一個月時間尺度的采樣數(shù)據(jù)，包含兩個大潮和小潮，且有足夠的數(shù)據(jù)長度，可以使計(jì)算結(jié)果相對穩(wěn)定［11］。所以，筆者選取2011年泉州地區(qū)地下水位數(shù)據(jù)，采用維尼迪可夫（Venedlkov）調(diào)和分析方法逐月分析計(jì)算水位對固體潮M2波的響應(yīng)能力，以此來檢驗(yàn)泉州地區(qū)上述井點(diǎn)井水位數(shù)據(jù)的內(nèi)在質(zhì)量［9］。

圖1 泉州地區(qū)地下水位日均值動態(tài)變化圖

表2 2011年泉州地區(qū)地下水位M2 潮汐因子計(jì)算結(jié)果

表3 2011年泉州地區(qū)地下水位M2 潮汐因子相對誤差計(jì)算結(jié)果

從表2及表3可以看出：

（1）永春1號井潮汐因子變化在9.38～17.88之間，平均值為13.10，其幅度與泉州地區(qū)其它觀測井相比較大，表明該井的靈敏度很高；相對誤差在0.028～0.046之間，平均值為0.038，說明動水位觀測相對靜水位觀測要靈敏很多。

（2）晉江地辦井、惠安地辦井、惠安凈峰井的潮汐因子基本保持在0.3以下，這3個觀測井的潮汐效應(yīng)不明顯，說明其觀測資料不能較好反映該地區(qū)地下構(gòu)造應(yīng)力變化，其可利用價(jià)值不高。南安地辦井的潮汐因子0.489～0.825，起伏較大，其潮汐因子的相對誤差也是最大的，說明該井的觀測資料可利用價(jià)值也不太理想。

（3）不同井點(diǎn)的潮汐可能受多種因素的影響，如井的封閉性、井孔觀測層巖性及所處觀測位置等。泉州局井、晉江深滬井、晉江安海井、泉港一中井、石獅赤湖井、永春一號井等井點(diǎn)的潮汐因子均大于2.0，其M2波的相對誤差也較小，均有較強(qiáng)的記錄固體潮汐效應(yīng)的能力。安溪1 號井的潮汐因子不大，但相對誤差較小，資料亦比較穩(wěn)定可靠。

柔性復(fù)合材料是由紡織材料經(jīng)過機(jī)織、針織、非織造、編織等加工成基材作為增強(qiáng)材料，再經(jīng)涂層、壓延、復(fù)合等加工制成的具有柔性，以及高強(qiáng)、質(zhì)輕、耐用、透光、耐腐蝕等優(yōu)良機(jī)械和物理化學(xué)性能的復(fù)合材料，可廣泛應(yīng)用于建筑頂棚、交通運(yùn)輸、充氣材料、農(nóng)用材料等領(lǐng)域。

2.3 氣壓對水位動態(tài)變化的影響分析

泉州地區(qū)的井孔都屬于深井，且封閉性較好，受降雨垂直滲入補(bǔ)給的影響較小。然而，“在水位觀測中不受氣壓影響的井，幾乎是不存在的，只是影響程度不同罷了?！币虼嗽诘卣鸱治鲱A(yù)報(bào)當(dāng)中，人們采用了不同的數(shù)理方法，希望將它從水位觀測值中扣除得越干凈越好［12］。筆者研究了泉州地區(qū)帶有氣象三要素的測點(diǎn)水位氣壓效應(yīng)，發(fā)現(xiàn)各測點(diǎn)水位與氣壓之間規(guī)律性不強(qiáng)。但經(jīng)過分析整個地區(qū)井網(wǎng)2011年1—12月水位日均值與氣壓日均值的相關(guān)關(guān)系，發(fā)現(xiàn)泉州井、晉江安海井、晉江深滬井、石獅赤湖井及泉港一中井上述5個（永春一號井其井孔類型為動水位）具有較強(qiáng)記錄固體潮能力的測點(diǎn)，與氣壓基本呈負(fù)相關(guān)關(guān)系，氣壓系數(shù)均較大，但氣壓系數(shù)本身并非常數(shù)，其變化幅度也較大，沒有固定規(guī)律，其中以表4中最有代表性的石獅赤湖井為例。而晉江地辦井、惠安地辦井以及惠安凈峰井等記錄固體潮能力較弱的井點(diǎn)在與氣壓的相關(guān)關(guān)系分析中，基本表現(xiàn)為不相關(guān)。由此可見，記錄固體潮能力愈強(qiáng)的井孔往往受到的氣壓效應(yīng)越大，如何去除水位觀測值中氣壓的影響，需要我們正確地選擇時段以及運(yùn)用合理的數(shù)理方法。

表4 石獅赤湖井2011年水位與氣壓日均值相關(guān)關(guān)系計(jì)算結(jié)果

2.4 周圍環(huán)境對水位觀測的影響

泉州地區(qū)上述水位觀測較為完整連續(xù)的11個測點(diǎn)中，根據(jù)歷年各測點(diǎn)的年終資料會商，泉州局井、惠安地辦井及安溪1號井這3個測點(diǎn)的水位觀測值長期受到周邊環(huán)境干擾較大。泉州局井、惠安地辦井均處在各自的新興開發(fā)區(qū)內(nèi)，自開始觀測以來附近多處工地多次斷斷續(xù)續(xù)地施工，而安溪1號井東面的城東大水閘不定期的畜水與放水對測井的水位資料動態(tài)變化造成一定的影響。

3 泉州地區(qū)地下水位資料動態(tài)特征分析

3.1 多年趨勢動態(tài)對比分析

對觀測井不少于一年正常觀測數(shù)據(jù)的計(jì)算，可以獲得水位動態(tài)及潮汐響應(yīng)函數(shù)的基本特征，這些參數(shù)可稱之為觀測井的特征參數(shù)。這些特征參數(shù)可以客觀地反映觀測井的基本素質(zhì)和周邊環(huán)境的影響，因而也可用于定量評價(jià)各個井的內(nèi)在質(zhì)量的優(yōu)劣［6］。然而由于各觀測點(diǎn)的觀測值高低不一，故筆者對2006—2011年泉州地區(qū)地下水位數(shù)字資料相對年變幅做了分析統(tǒng)計(jì)，在每個井點(diǎn)全年365個日均值中找出正負(fù)最大值后，將兩者的絕對值相減，即得該井點(diǎn)水位變化的年變幅，再以年變幅除以年均值得到相對年變幅。并對這幾年的相對年變幅做了均方差的分析計(jì)算（表5），得到泉州地區(qū)上述測井?dāng)?shù)字化水位多年及年趨勢變化的幾個特征：

（1）從表5得出：若相對年變幅較大，說明該井的年動態(tài)變化起伏較大，反之，則較小。如2006年年變規(guī)律起伏最大的惠安地辦井，相對年變幅達(dá)到2.0238 倍，2010年年變規(guī)律起伏最小的泉港一中井，其相對年變幅度為0.106倍。

（2）通過對2006—2011年泉州地區(qū)上述測點(diǎn)地下水位數(shù)字資料相對年變幅的均方差計(jì)算可以得出其多年動態(tài)趨勢變化規(guī)律。多年相對年變幅的均方差越小，說明多年趨勢變化越穩(wěn)定，根據(jù)表5 和圖1，筆者對2006—2011年泉州地區(qū)地下水位數(shù)字資料各觀測點(diǎn)各自多年動態(tài)變化穩(wěn)定性進(jìn)行分析認(rèn)為：永春1號井、石獅赤湖井、泉港一中井穩(wěn)定性最好，屬于平穩(wěn)變化型；泉州局井、晉江地辦井、晉江深滬井、晉江安海井、惠安凈峰井的多年動態(tài)起伏相對較大，基本屬于下降—上升—下降型；安溪1號井則屬于下降—平穩(wěn)—漸升型。

（3）泉州局井、惠安地辦井及安溪1號井受干擾較多，其相對年變幅的均方差亦較大。

表5 泉州地區(qū)水位觀測資料多年穩(wěn)定性分析

3.2 年變化動態(tài)特征對比分析

從圖1可以看出泉州地區(qū)各測點(diǎn)間存在著較為同步的年變規(guī)律，在2012年4月之后，屬多年動態(tài)平穩(wěn)型的永春1號井、石獅赤湖井、泉港一中井，水位趨勢仍保持較為平穩(wěn)的動態(tài)變化，其余測點(diǎn)（除泉州局井外）的水位走勢均呈下降趨勢。對于具有較好記錄固體潮能力的泉州井、晉江安海井、晉江深滬井、永春一號井、石獅赤湖井及泉港一中井，利用泉州地區(qū)自開始觀測以來上述測點(diǎn)水位月均值最高值進(jìn)行統(tǒng)計(jì)得到：泉州地區(qū)上述6個測點(diǎn)水位的月均最高值95%發(fā)生在6—10月（圖2）。結(jié)合圖1可以得出：一般從3月份起各測點(diǎn)水位開始呈上升趨勢，6—10月內(nèi)達(dá)到峰值，10月后開始下降，一般在次年2月降到最低點(diǎn)。而圖1與圖2表明：泉州局井的水位在2011年12月至2012年4月存年畸變異常。對上述6個測點(diǎn)具有一致動態(tài)曲線特征對比分析發(fā)現(xiàn)：泉州局井水位在2011年04月至2012年9月連續(xù)上升，但其余5個測點(diǎn)在7月后有3個下降，2個較為平穩(wěn)（圖1），表明該年變異常并非區(qū)域應(yīng)力場發(fā)生變化所致。事實(shí)上，2012年初，泉州局曾組織相關(guān)人員經(jīng)過對周邊環(huán)境的實(shí)地考察，最后判定2011年4月起至2012年3月附近幾處工地的施工可能是對該井水位造成干擾的原因①關(guān)于泉州市地震局水位異常情況的報(bào)告。

圖2 泉州地區(qū)各測點(diǎn)歷年水位最高月均值統(tǒng)計(jì)

3.3 月變化動態(tài)特征對比分析

對于各觀測點(diǎn)的觀測值高低不一的情況，可以用各測點(diǎn)相對變化幅度的相對均方差對泉州地區(qū)水位觀測資料每月及每日動態(tài)變化特征進(jìn)行分析。相對均方差是給定時間段內(nèi)數(shù)列的均方差與數(shù)列的均值之比［8］：

式中：S為均方差，為均值，選定“日”為單位時間段，那么日相對均方差就是每天水位整點(diǎn)值相對變幅的均方差與日均值之比，用EIS2000軟件［9］計(jì)算，求取各井點(diǎn)每天的相對均方差，再以式（1）計(jì)算得出各測點(diǎn)每天的相對均方差與月均值之比，用該統(tǒng)計(jì)結(jié)果來研究泉州地區(qū)水位觀測資料的曲線變化動態(tài)特征。如果觀測水位的日相對均方差小，說明其日變幅度規(guī)律較好，月動態(tài)規(guī)律性好，其月相對均方差變化較小。

由于泉州局井處在一個新興開發(fā)區(qū)內(nèi)，自該井開始觀測以來附近多處工地多次斷斷續(xù)續(xù)地施工，對該井的水位資料動態(tài)變化造成一定的影響，并從表6看出泉州局井月相對均方差遠(yuǎn)大于泉州地區(qū)其它觀測點(diǎn)，表明測點(diǎn)附近存在干擾的事實(shí)印證了其月動態(tài)規(guī)律不好的現(xiàn)象。

表6 2011泉州地區(qū)水位觀測資料穩(wěn)定性分析

4 結(jié)論

（1）泉州數(shù)字化水位觀測目前共有17個測點(diǎn)，其中3個測點(diǎn)為2011年以后新增，觀測時間較長的有14個測點(diǎn)，其中11個測點(diǎn)資料較為連續(xù)完整可靠，利用2006—2011年這11個觀測點(diǎn)的日值圖并結(jié)合相對均方差方法對多年、年、月及日動態(tài)變化進(jìn)行定量分析，可以找出其多年、年、月及日動態(tài)特征。其中泉州井由于受到2011年附近工地施工干擾影響，其月相對均方差偏大，該井該年的月動態(tài)規(guī)律性不好；安溪1號井的多年動態(tài)特征起伏較大，規(guī)律性不強(qiáng)；其余的觀測井?dāng)?shù)字化水位資料多年、年、月及日動態(tài)特征規(guī)律性較強(qiáng)。對于年變規(guī)律清晰、相對動態(tài)變化較為同步且記錄固體潮汐效應(yīng)能力較強(qiáng)的觀測井，比如泉州局井、晉江深滬井、晉江安海井、石獅赤湖井、泉港一中井、永春1號井等，由于這些測點(diǎn)同處于NNE－NE向的長樂－詔安斷裂帶和NW向的永安－晉江斷裂帶上，之間存在的同步變化的動態(tài)變化特征共同反映了該構(gòu)造區(qū)域活動引起的地下流體變化的信息，形成了一個非常直觀的地震前兆監(jiān)控網(wǎng)。所以需要有效識別這些測點(diǎn)水位的正常動態(tài)，才有可能在地震中短期或短期前兆異常監(jiān)測中發(fā)揮一定的效能。

（2）永春井井孔類型屬于動水位類型，其觀測井井孔條件較好、年變規(guī)律清晰、動態(tài)變化平穩(wěn)，是福建省少有的自流溫泉井，不僅能反映含水層孔隙壓力的變化，且由于熱水井或水中含有較多氣體，“溫度效應(yīng)”和“氣體效應(yīng)”又對水位變化產(chǎn)生次生放大作用，所以較泉州地區(qū)其它靜水位觀測要靈敏得多。

（3）泉州地區(qū)的水位觀測井類型（除永春井）均為封閉的承壓含水井，地下水位的觀測值動態(tài)變化中包含了地質(zhì)過程、固體潮、大氣壓力及周圍地下開采與施工等因素。而對含水應(yīng)力層反應(yīng)愈靈敏的井孔（如泉州井、晉江安海井、晉江深滬井、石獅赤湖井及泉港一中井），對固體潮、大氣壓力及周圍環(huán)境變化的記錄愈清晰，這其中除了地質(zhì)過程，其它均可看作干擾因素，應(yīng)采取較為合理的數(shù)理方法去有效地排除地下水位資料中的各種干擾因素，才能提取可靠的地震信息。

［1］ 張昱，李春森，趙小茂，等.甘肅及鄰區(qū)數(shù)字化水位觀測資料的應(yīng)用與地震監(jiān)測效能評估［J］.地震，2008，28（03）：103－109.

［2］ 陳小云，曹軼，張清秀，等.福建省數(shù)字化水位動態(tài)特征研究［J］.華北地震科學(xué)，2009，27（03）：51－55.

［3］ 張翊寧，陳瑩.福建省水氡觀測資料變化規(guī)律與穩(wěn)定性評價(jià)［J］.華南地震，2005，25（04）：66－72.

［4］ 張清秀，鄭小菁，林木金.永春井水位對遠(yuǎn)大地震記震能力及同震效應(yīng)特征分析［J］.華南地震，2011，31（03）：33－40.

［5］ 張素欣，張子廣，劉俊明，等.數(shù)字化水位觀測資料的應(yīng)用研究［J］.地震，2002，22（04）：89－93.

［6］ 唐九安，沈旭章，高安泰，等.承壓水位觀測井觀測數(shù)據(jù)評價(jià)方法探討［J］.大地測量與地球動力學(xué)，2012，32（04）：37－40.

［7］ 張昱，鐘美嬌，楊曉鵬，等.平?jīng)鯟11井地下流體數(shù)字化觀測資料的分析［J］.地震，2008，28（04）：128－135.

［8］ 萬永芳，李健梅.廣東井水位和水溫對汶川8.0級地震的同震響應(yīng)研究［J］.華南地震，2010，30（S）：53－60.

［9］ 蔣駿，李勝樂，張雁濱，等.地震前兆信息處理與軟件系統(tǒng)［M］.北京.地震出版社，2000.

［10］ 晏銳，張立，簡春林.云南曲靖井水位潮汐動態(tài)特征分析［J］.地震學(xué)報(bào)，2012，34（3）：363－373.

［11］ 陳瑩，黃輔瓊，朱石軍，等.采樣數(shù)據(jù)長度對潮汐因子計(jì)算精度的影響［J］.地震地磁觀測與研究，2012，33（5／6）：146－151.

［12］ 董守玉，賈化周，秦清娟，等.水位氣壓系數(shù)的探討［J］.地震研究，1987.10（01）：63－69.