陜西漢中地區(qū)定點(diǎn)形變臺(tái)站觀測(cè)數(shù)據(jù)質(zhì)量評(píng)價(jià)

王穎昌 岳鵬 朱宣明 高余慶

（1.陜西省地震局，陜西 西安 710068；2.安康地震監(jiān)測(cè)中心站，陜西 安康 725000）

0 引言

漢中地區(qū)由龍門山北東向斷裂帶和略陽—洋縣東西向斷裂圍限，盆地與周圍山地的差異運(yùn)動(dòng)強(qiáng)烈，地震活動(dòng)較為頻繁。最大歷史地震為5.5級(jí)地震（洋縣1624年、1635年，漢中1636年，略陽1908年），最大近代地震為1967年南鄭5級(jí)地震，該地區(qū)現(xiàn)有定點(diǎn)形變臺(tái)站兩個(gè)，分別是漢中形變臺(tái)（南鄭813廠）和寧強(qiáng)形變臺(tái)，共架設(shè)各類傾斜應(yīng)變觀測(cè)儀器5臺(tái)套，對(duì)監(jiān)測(cè)和研究該區(qū)域地表變形、斷層活動(dòng)、地震前兆及深部構(gòu)造演化等具有重要的科學(xué)價(jià)值。本研究通過對(duì)分別架設(shè)在兩個(gè)臺(tái)站的兩套DSQ型水管傾斜儀和兩套SS?YⅡ型洞體應(yīng)變伸縮儀近5年的連續(xù)率、完整率、年零漂、長(zhǎng)周期擬合相對(duì)噪聲水平和資料質(zhì)量評(píng)定內(nèi)精度（M2波潮汐因子均方差mγ和M2波潮汐因子相對(duì)中誤差σα/α）等技術(shù)指標(biāo)進(jìn)行計(jì)算，總體評(píng)價(jià)觀測(cè)資料質(zhì)量，以期為無人值守地震臺(tái)站監(jiān)測(cè)效能評(píng)估和地震預(yù)測(cè)預(yù)報(bào)提供參考。

1 臺(tái)站及測(cè)項(xiàng)概況

1.1 南鄭813廠形變臺(tái)基本情況及地震地質(zhì)構(gòu)造背景

漢中形變臺(tái)（南鄭813廠）地傾斜（應(yīng)變）觀測(cè)采用全數(shù)字化觀測(cè)儀器，為無人值守臺(tái)站。該臺(tái)位于漢中市南鄭縣813廠軍用山洞內(nèi)，海拔高度564 m，山洞進(jìn)深183 m，臺(tái)基巖性為偉晶花崗巖，堅(jiān)硬完整，致密均勻，巖層傾角8°左右，地理坐標(biāo)為東經(jīng)107°43′、北緯33°1′。該山洞地處龍門山北東向斷裂帶的北東端，距娘娘山—秦家壩斷層1 km，距漢江3 km，周圍無大型建筑、采礦爆破點(diǎn)、輸變電站等干擾。最大歷史地震為5.5級(jí)地震（洋縣1624年、1635年，漢中1636年，略陽1908年），最大近代地震為1967年南鄭5級(jí)地震。

1.2 漢中形變臺(tái)（南鄭813廠）測(cè)項(xiàng)分布情況

漢中形變臺(tái)（南鄭813廠）現(xiàn)安裝有DSQ型水管傾斜儀、SS?YⅡ型洞體應(yīng)變伸縮儀、VP型寬頻帶傾斜儀。儀器測(cè)項(xiàng)分布如圖1所示。

1.3 寧強(qiáng)形變臺(tái)基本情況及地震地質(zhì)構(gòu)造背景

寧強(qiáng)形變臺(tái)位于陜西省寧強(qiáng)縣漢源鎮(zhèn)滴水鋪村石梁溝，距離寧強(qiáng)縣城約14 km，建成于2011年，現(xiàn)為無人值守臺(tái)站。其地理坐標(biāo)為東經(jīng)106°24′、北緯32°92′，海拔高度865 m，屬北亞熱帶內(nèi)陸性季風(fēng)氣候。現(xiàn)有觀測(cè)山洞為村民開挖的石炭洞，洞內(nèi)已無石炭，為沉積灰?guī)r，巖性致密完整，覆蓋30 m，青岡木植被。2011年經(jīng)優(yōu)化改造后，適用于數(shù)字化傾斜儀和洞體應(yīng)變儀觀測(cè)。臺(tái)站地處秦嶺山脈中段南坡，地貌為山地，構(gòu)造上處于平武－青川斷裂帶東南盤上。

1.4 寧強(qiáng)形變臺(tái)測(cè)項(xiàng)分布概況

寧強(qiáng)形變臺(tái)現(xiàn)安裝有DSQ型水管傾斜儀、SS?YⅡ型洞體應(yīng)變伸縮儀兩套儀器。儀器測(cè)項(xiàng)分布圖如圖2所示。

2 觀測(cè)資料質(zhì)量評(píng)價(jià)

2.1 資料選取與評(píng)價(jià)指標(biāo)

本研究選取了漢中形變臺(tái)（南鄭813廠）和寧強(qiáng)形變臺(tái)水管儀2017—2021年DSQ型水管傾斜儀和SS?YⅡ型洞體應(yīng)變伸縮儀的觀測(cè)數(shù)據(jù)產(chǎn)品，包括整點(diǎn)值、日均值和五日均值等，參考地殼形變學(xué)科組擬定的技術(shù)指標(biāo)進(jìn)行對(duì)比研究，選定的洞體傾斜觀測(cè)質(zhì)量評(píng)價(jià)指標(biāo)有連續(xù)率、完整率、年零漂、M2波潮汐因子均方差mγ、相對(duì)噪聲水平，洞體應(yīng)變觀測(cè)質(zhì)量指標(biāo)有連續(xù)率、完整率、年零漂、M2波潮汐因子相對(duì)中誤差σα/α、相對(duì)噪聲水平。

2.2 數(shù)據(jù)連續(xù)率、完整率和年零漂

觀測(cè)數(shù)據(jù)連續(xù)、完整是地震研究對(duì)觀測(cè)數(shù)據(jù)的基本要求，也是檢驗(yàn)觀測(cè)儀器運(yùn)行質(zhì)量的重要條件，按照地殼形變學(xué)科組評(píng)分標(biāo)準(zhǔn)一般大于99.5%視為優(yōu)秀，連續(xù)率和完整率的計(jì)算公式分別為式（1）和式（2）。

按式（1）和（2）分別計(jì)算漢中形變臺(tái)（南鄭813廠）和寧強(qiáng)形變臺(tái)DSQ型水管傾斜儀和SS?YⅡ型洞體應(yīng)變伸縮儀2017—2021年的數(shù)據(jù)連續(xù)率（原始數(shù)據(jù)）和完整率（預(yù)處理數(shù)據(jù)），結(jié)果如表1、表2所示。由表1可知，數(shù)據(jù)連續(xù)率方面，漢中臺(tái)水管儀除2020年因更換整套觀測(cè)系統(tǒng)導(dǎo)致連續(xù)率為98.2%以外，其他年份數(shù)據(jù)連續(xù)率均在99.5%以上，寧強(qiáng)臺(tái)水管儀2017—2021年數(shù)據(jù)連續(xù)率均在99.6%以上，達(dá)到優(yōu)秀標(biāo)準(zhǔn)；數(shù)據(jù)完整率方面，201—2021年漢中臺(tái)水管儀完整率在97.7%～99.7%之間，平均值98.8%，寧強(qiáng)臺(tái)水管儀完整率在97.9%～99.7%之間，平均值99.0%。由表2可知，漢中臺(tái)和寧強(qiáng)臺(tái)伸縮儀2017—2021年數(shù)據(jù)連續(xù)率為97.9%～100%，2019年寧強(qiáng)臺(tái)伸縮儀兩分量數(shù)據(jù)均完整連續(xù)，無缺記、斷記；數(shù)據(jù)完整率方面，漢中臺(tái)在2018年對(duì)伸縮儀觀測(cè)系統(tǒng)進(jìn)行了整體更換，儀器安裝調(diào)試過程產(chǎn)生中產(chǎn)生了較多的錯(cuò)誤數(shù)據(jù)，造成當(dāng)年數(shù)據(jù)完整率偏低，其他年份均在98%以上，寧強(qiáng)臺(tái)伸縮儀2017—2021年數(shù)據(jù)完整率僅為89%～99%，整體來看DSQ水管傾斜儀運(yùn)行質(zhì)量寧強(qiáng)臺(tái)略優(yōu)于漢中臺(tái)，SS?YⅡ型洞體應(yīng)變伸縮儀運(yùn)行質(zhì)量漢中臺(tái)優(yōu)于寧強(qiáng)臺(tái)，經(jīng)查詢觀測(cè)日志可知，影響觀測(cè)數(shù)據(jù)連續(xù)性和完整性的主要原因有：供電系統(tǒng)故障，數(shù)據(jù)采集器故障，傳感器故障和人員進(jìn)洞調(diào)零，維修儀器干擾寧強(qiáng)臺(tái)伸縮儀以及受到山洞內(nèi)鼠患影響。

年零漂用于衡量觀測(cè)儀器及基墩穩(wěn)定程度或地殼繼承性新構(gòu)造運(yùn)動(dòng)，在一定程度上反映了儀器的穩(wěn)定性及所測(cè)方向的微小形變，可采用日均值法和潮汐值法進(jìn)行計(jì)算［1］。本研究采用日均值法計(jì)算漢中形變臺(tái)（南鄭813廠）和寧強(qiáng)形變臺(tái)DSQ型水管傾斜儀和SS?YⅡ型洞體應(yīng)變伸縮儀的年零漂，即將1年內(nèi)12月31日與1月1日日均值的差值作為該年度零漂值。計(jì)算結(jié)果見表1和表2。結(jié)果顯示，寧強(qiáng)臺(tái)水管儀NS向年零漂較小，EW向次之，漢中臺(tái)水管儀兩分量年零漂相對(duì)較大，其中EW向略小于NS向；寧強(qiáng)臺(tái)伸縮儀EW向年零漂波動(dòng)較大，NS向和漢中臺(tái)伸縮儀兩個(gè)分量年零漂較小且相對(duì)穩(wěn)定，漢中臺(tái)2018年和2020年分別對(duì)伸縮儀和水管儀觀測(cè)系統(tǒng)進(jìn)行更換后，年零漂幅度未見明顯變化。

表1 DSQ型水管傾斜儀連續(xù)率、完整率和年零漂統(tǒng)計(jì)

表2 SS?YⅡ型洞體應(yīng)變伸縮儀連續(xù)率、完整率和年零漂統(tǒng)計(jì)

2.3 相對(duì)噪聲水平

長(zhǎng)周期擬合相對(duì)噪聲水平M1是衡量觀測(cè)資料長(zhǎng)期穩(wěn)定性的一項(xiàng)重要指標(biāo)，是用契比雪夫多項(xiàng)式與一個(gè)年度觀測(cè)資料的5日均值作30階擬合得到的［2］，計(jì)算公式為式（3）、式（4）。

式中：Tn(x)稱為契比雪夫多項(xiàng)式；x為五日均值；C0為常數(shù)；Cn為擬合系數(shù)；n為系列號(hào)；m為提取的五日均值總個(gè)數(shù)［3］。

根據(jù)式（3）、式（4）分別計(jì)算漢中形變臺(tái)（南鄭813廠）和寧強(qiáng)形變臺(tái)DSQ型水管傾斜儀和SS?YⅡ型洞體應(yīng)變伸縮儀的相對(duì)噪聲水平M1，計(jì)算結(jié)果見表3和表4，可以看出2017—2021年寧強(qiáng)臺(tái)水管儀兩分量相對(duì)噪聲M1在0.001 1～0.003 7之間，長(zhǎng)期保持較低水平且較為穩(wěn)定，漢中臺(tái)水管儀兩分量相對(duì)噪聲M1在0.000 9～0.063 2之間，相對(duì)噪聲水平偏大；2017—2021年漢中臺(tái)伸縮儀兩分量相對(duì)噪聲M1平均值分別為：NS向0.003 8、EW向0.007 9，均小于寧強(qiáng)臺(tái)，由相對(duì)噪聲水平越低，則儀器穩(wěn)定性越好可知，寧強(qiáng)臺(tái)水管儀穩(wěn)定性優(yōu)于漢中臺(tái)，漢中臺(tái)伸縮儀穩(wěn)定性優(yōu)于寧強(qiáng)臺(tái)，同時(shí)觀測(cè)系統(tǒng)的更新亦未見對(duì)相對(duì)噪聲水平有所提升。

表3 DSQ型水管傾斜儀相對(duì)噪聲M1計(jì)算結(jié)果（單位："）

表4 SS?YⅡ型洞體應(yīng)變伸縮儀相對(duì)噪聲M1計(jì)算結(jié)果（單位：1×10?6）

2.4 資料質(zhì)量評(píng)定內(nèi)精度

DSQ型水管儀觀和SS?YⅡ型洞體應(yīng)變伸縮儀觀測(cè)資料精度分別以調(diào)和分析M2波潮汐因子均方差mγ和調(diào)和分析M2波潮汐因子相對(duì)中誤差σα/α作為資料質(zhì)量評(píng)定內(nèi)精度計(jì)算指標(biāo)。采用Vene?dikov調(diào)和分析方法［4］，利用EIS2000地震前兆信息處理與軟件系統(tǒng)［5］，分別計(jì)算漢中形變臺(tái)（南鄭813廠）和寧強(qiáng)形變臺(tái)DSQ型水管傾斜儀的M2波潮汐因子均方差mγ，如表5所示；分別計(jì)算漢中形變臺(tái)（南鄭813廠）和寧強(qiáng)形變臺(tái)SS?YⅡ型洞體應(yīng)變伸縮儀M2波潮汐因子相對(duì)中誤差σα/α，如表6所示。

表5 DSQ型水管傾斜儀M2波潮汐因子均方差計(jì)算結(jié)果

表6 SS?YⅡ型洞體應(yīng)變伸縮儀M2波潮汐因子相對(duì)中誤差σα/α計(jì)算結(jié)果

觀測(cè)資料精度與場(chǎng)地環(huán)境、儀器最初架設(shè)情況、安裝點(diǎn)的地理位置、地質(zhì)構(gòu)造、觀測(cè)系統(tǒng)、人為干擾、自然環(huán)境干擾等均密切相關(guān)［6］。從計(jì)算結(jié)果可以看出，總體上看漢中臺(tái)水管儀mγ為0.001 5～0.006 5，寧強(qiáng)臺(tái)mγ為0.001 5～0.004 7，均滿足地殼形變學(xué)科組mγ≤0.02的技術(shù)指標(biāo)，以2021年為例，和陜西省內(nèi)其他同類形變臺(tái)站相比，寧強(qiáng)臺(tái)水管儀兩分量mγ平均值0.002 5，僅高于乾陵臺(tái)（0.002 0），位列全省第二，漢中臺(tái)mγ（0.003 1）與隴縣臺(tái)（0.003 1）持平，優(yōu)于華陰臺(tái)（0.003 7）、西安臺(tái)（0.004 6）、寶雞臺(tái)（0.006 5）和平利臺(tái)（0.011 8）；漢中臺(tái)伸縮儀σα/α為0.005 3～0.012 6，寧強(qiáng)臺(tái)σα/α為0.009 4～0.046 0，同樣符合地殼形變學(xué)科組σα/α≤0.05的技術(shù)指標(biāo)要求。觀測(cè)系統(tǒng)更新對(duì)資料質(zhì)量評(píng)定內(nèi)精度提升同樣不明顯。

3 結(jié)論

綜上，通過對(duì)漢中形變臺(tái)（南鄭813廠）和寧強(qiáng)形變臺(tái)2017—2021年DSQ型水管傾斜儀和SS?YⅡ型洞體應(yīng)變伸縮儀的觀測(cè)數(shù)據(jù)進(jìn)行對(duì)比分析得到以下結(jié)論。

①在儀器運(yùn)行質(zhì)量方面，寧強(qiáng)臺(tái)和漢中臺(tái)的DSQ水管傾斜儀總體運(yùn)行平穩(wěn)，數(shù)據(jù)連續(xù)性好，受觀測(cè)系統(tǒng)故障和人為事件干擾影響，寧強(qiáng)臺(tái)和漢中臺(tái)SS?YⅡ型洞體應(yīng)變伸縮儀觀測(cè)數(shù)據(jù)連續(xù)率和完整率偏低。

②在儀器穩(wěn)定性方面，通過比較年零漂和相對(duì)噪聲水平M1兩個(gè)指標(biāo)，可以看出寧強(qiáng)臺(tái)DSQ水管傾斜儀穩(wěn)定性優(yōu)于漢中臺(tái)，漢中臺(tái)SS?YⅡ型洞體應(yīng)變伸縮儀穩(wěn)定性優(yōu)于寧強(qiáng)臺(tái)。

③觀測(cè)精度方面，漢中臺(tái)和寧強(qiáng)臺(tái)的DSQ水管傾斜儀和SS?YⅡ型洞體應(yīng)變伸縮儀均符合地殼形變學(xué)科組技術(shù)指標(biāo)要求，寧強(qiáng)臺(tái)水管儀觀測(cè)精度在全省同類臺(tái)站中處于優(yōu)秀水平。觀測(cè)系統(tǒng)的更新對(duì)年零漂、相對(duì)噪聲水平M1和觀測(cè)精度的提升作用有限，結(jié)合前人研究成果，觀測(cè)場(chǎng)地地質(zhì)構(gòu)造，洞室環(huán)境和基巖的選擇直接決定儀器內(nèi)在精度。

④整體來看，傾斜觀測(cè)數(shù)據(jù)質(zhì)量寧強(qiáng)臺(tái)優(yōu)于漢中臺(tái)，應(yīng)變觀測(cè)數(shù)據(jù)質(zhì)量，漢中臺(tái)優(yōu)于寧強(qiáng)臺(tái)，為更加科學(xué)全面客觀地評(píng)價(jià)形變臺(tái)站觀測(cè)數(shù)據(jù)質(zhì)量，還應(yīng)引入同類儀器數(shù)據(jù)相關(guān)性、同震響應(yīng)情況等指標(biāo)進(jìn)行分析。