高熹微，張 騫，卜玉菲，趙衛(wèi)紅*，郭宇鑫，李 玉

（1. 江蘇省地震局高郵地震臺，江蘇 高郵 225600；2. 江蘇省地震局，南京 210014）

0 引言

針對地電臺站建設(shè)和觀測儀器入網(wǎng)，中國地震局和電磁學(xué)科技術(shù)管理組編制了一系列的行業(yè)標(biāo)準(zhǔn)和技術(shù)規(guī)范[1-7]，提出了明確的技術(shù)要求，如地震觀測儀器進(jìn)網(wǎng)技術(shù)要求[6]，并給出了相應(yīng)的測試方法，但未將地電儀器連入經(jīng)過批準(zhǔn)的野外觀測場地進(jìn)行實地測試來驗證儀器實際野外觀測性能。因此按照地電建設(shè)規(guī)范[2,5]要求，建設(shè)了高郵地電比測平臺和平?jīng)龅仉姳葴y平臺，在通過認(rèn)證和考核后，承擔(dān)地電學(xué)科儀器的臺站野外比測工作。從環(huán)境適應(yīng)性、觀測一致性、運(yùn)行穩(wěn)定性等方面對地電儀器進(jìn)行對比測試，同時編寫地電觀測儀器技術(shù)比測規(guī)范，形成地電學(xué)科專業(yè)設(shè)備臺站比測完整的流程和方法，并設(shè)置技術(shù)指標(biāo)值來衡量比測儀器，確保地震監(jiān)測地電學(xué)科專業(yè)設(shè)備的性能穩(wěn)定和得到高質(zhì)量的觀測數(shù)據(jù)。

在參考地電行業(yè)標(biāo)準(zhǔn)的基礎(chǔ)上，編寫了地電觀測儀器比測技術(shù)規(guī)范（第一部分：地電阻率儀；第二部分：地電場儀），建立了臺站地電儀器比測流程和方法等內(nèi)容，并將地電場儀和地電阻率儀的《地震監(jiān)測專業(yè)設(shè)備定型測試大綱》（https://www.cencc.ac.cn/cenc/_300 647/352 611/index.html）中設(shè)立的4個技術(shù)指標(biāo)擴(kuò)充為8個來分析比測儀器產(chǎn)出數(shù)據(jù)的穩(wěn)定性和一致性。針對這8個技術(shù)指標(biāo)設(shè)置值是否合理及是否滿足現(xiàn)今在用地電儀器觀測要求的問題，將ZD9A-Ⅱ和GEF-Ⅱ兩種不同型號的儀器連接同一地電場觀測裝置進(jìn)行對比觀測。利用該試驗數(shù)據(jù)來計算地電場比測規(guī)范技術(shù)指標(biāo)值，平?jīng)雠_ZD8M地電阻率儀對比觀測數(shù)據(jù)計算地電阻率比測規(guī)范技術(shù)指標(biāo)值，通過試驗數(shù)據(jù)和歷史觀測數(shù)據(jù)計算所得結(jié)果，來驗證技術(shù)指標(biāo)設(shè)置值的合理性，為地電觀測儀器比測技術(shù)規(guī)范編寫中技術(shù)指標(biāo)值的制定提供基礎(chǔ)支撐。從而為后續(xù)的地電臺站比測提供依據(jù)，并為規(guī)范地電觀測儀器和產(chǎn)出用于地震監(jiān)測、預(yù)報的高質(zhì)量地電觀測數(shù)據(jù)提供保障。

1 數(shù)據(jù)選取

地電場所用數(shù)據(jù)為2019年10月18日至2020年4月18日期間ZD9A-Ⅱ、GEF-Ⅱ地電場儀連接臺上現(xiàn)有的地電場觀測裝置按照地電學(xué)科現(xiàn)行要求進(jìn)行觀測所產(chǎn)出的觀測數(shù)據(jù)。選取2019年12月22—23日的觀測數(shù)據(jù)（圖1）進(jìn)行對比發(fā)現(xiàn)，兩者對于地電阻率供電干擾的表現(xiàn)方式不一致。在對比觀測中將ZD9A-Ⅱ地電場儀作為參考儀器，用以驗證目前主流的ZD9A-Ⅱ地電場儀和新型的GEF-Ⅱ地電場儀產(chǎn)出的觀測數(shù)據(jù)計算結(jié)果是否滿足于技術(shù)規(guī)范中編寫的技術(shù)指標(biāo)值。其中該臺GEF-Ⅱ地電場儀已符合中國地震局定型測試大綱要求，通過定型測試。現(xiàn)有地電場觀測裝置呈L型布設(shè)，分別為NS向、EW向和N45°E向3個測道。

圖1 ZD9A-Ⅱ與GEF-Ⅱ儀器2 019年12月22—23日觀測數(shù)據(jù)

由于同場地地電阻率供電干擾，高郵地電比測平臺運(yùn)行的地電阻率產(chǎn)出的觀測數(shù)據(jù)量少，因而采用平?jīng)龅仉姳葴y平臺地電阻率觀測資料，選取了2019年5月18日至2019年8月15日兩臺ZD8M地電阻率觀測儀按照地電學(xué)科現(xiàn)行要求運(yùn)行產(chǎn)出的NS、EW向兩測道觀測數(shù)據(jù)（圖2）。從圖2中可以看出，兩者觀測數(shù)據(jù)差異很小，故以平?jīng)雠_現(xiàn)運(yùn)行的ZD8M地電阻率儀作為參考儀器。

圖2 地電阻率儀器2019年5月18至2019年8月15日對比觀測數(shù)據(jù)

2 計算方法及結(jié)果

2.1 地電場技術(shù)指標(biāo)值

2.1.1 儀器之間同裝置觀測結(jié)果每天的日相關(guān)系數(shù)r日

求解測道之間的相關(guān)系數(shù)，是地電場觀測數(shù)據(jù)可靠性的一種檢測方法[8]，多臺比測儀器在連接同度相關(guān)分析的方法[7]，給出r日計算公式（1）～（2）和技術(shù)指標(biāo)值：應(yīng)不小于0.9，符合率C應(yīng)大于95%。一觀測裝置的前提下，其同方位、同極距（同裝置）產(chǎn)出的觀測數(shù)據(jù)應(yīng)具有高度一致的變化趨勢，相關(guān)系數(shù)應(yīng)呈現(xiàn)高度相關(guān)性。因而，設(shè)立儀器之間同裝置觀測結(jié)果的日相關(guān)系數(shù)技術(shù)指標(biāo)和符合率，即求取不同儀器之間同方位、同極距的相關(guān)系數(shù)，來衡量多臺同型號、不同型號比測儀器之間所得的觀測數(shù)據(jù)一致性。引用《地震監(jiān)測專業(yè)設(shè)備定型測試大綱》中臺站比測地電場儀器之間的相關(guān)系數(shù)指標(biāo)值，并結(jié)合行業(yè)標(biāo)準(zhǔn)DB/T 34-2009附錄E中小時尺

將ZD9A-Ⅱ、GEF-Ⅱ所測數(shù)據(jù)中地電阻率供電時段的干擾數(shù)據(jù)舍棄。因為同場地觀測的地電阻率中存在供電干擾，ZD9A-Ⅱ使用設(shè)置的門限來抑制供電干擾[9]，GEF-Ⅱ沒有門限設(shè)置，因而計算兩者的技術(shù)指標(biāo)值時舍棄供電時段的數(shù)據(jù)。從儀器面板觀察，ZD9A-Ⅱ為當(dāng)前時間觀測后，存儲為下一分鐘觀測值；GEF-Ⅱ為當(dāng)前時間觀測后，存儲為當(dāng)前時間的觀測值，因而將ZD9A-Ⅱ儀器數(shù)據(jù)超前1分鐘參與計算。考慮同場地地電阻率供電對地磁、地電產(chǎn)生干擾，高郵地電比測平臺的地電阻率觀測設(shè)置為每天工作2～4次不等。這不固定的供電干擾摻雜在地電觀測數(shù)據(jù)中造成一定影響，因而在2019年10月17日停止觀測，地電場技術(shù)指標(biāo)計算所選取的數(shù)據(jù)從2019年10月18日開始。按照公式（1）～（3）計算日相關(guān)系數(shù)和符合率，得到兩套儀器各個方位長極距與長極距、短極距與短極距之間的每日相關(guān)系數(shù)（圖3）。兩臺儀器共六測道的日相關(guān)系數(shù)結(jié)果在0.8～1.0之間波動（異常點(diǎn)除外），統(tǒng)計六測道計算結(jié)果，兩種儀器之間的同裝置相關(guān)系數(shù)大于0.9的個數(shù)為340，符合率為33.51%。

圖3 ZD9A-Ⅱ與GEF-Ⅱ儀器同方位同極距日相關(guān)系數(shù)

對照圖3可看出，在12月22日兩臺儀器六測道的相關(guān)系數(shù)低至0.4附近。從圖1a～f中可見，兩者觀測數(shù)據(jù)并無數(shù)據(jù)突跳等畸變。根據(jù)兩者面板觀察，兩者存儲的數(shù)據(jù)相位應(yīng)存在1分鐘的差異，但從圖1a1中可以看出，12月22日23時20—59分的觀測數(shù)據(jù)兩者相位一致，兩臺儀器存在時間差；從圖1a2中看出，12月23日兩者的觀測數(shù)據(jù)存在著明顯的相位差異，其中GEF-Ⅱ儀器每天進(jìn)行4次自動校時。因而將12月21—23日的ZD9A-Ⅱ儀器觀測數(shù)據(jù)與GEF-Ⅱ儀器觀測數(shù)據(jù)按正常時間序列計算儀器之間同裝置相關(guān)系數(shù)（表1），12月22日正常時間序列計算結(jié)果優(yōu)于超前1分鐘計算結(jié)果，這說明兩種儀器存在時間差。在實際的觀測過程中，ZD9A-Ⅱ儀器使用時間長達(dá)十年多，雖然定時對兩臺儀器校時，也難以確保兩者的儀器時間完全一致，時間的不同步導(dǎo)致觀測數(shù)據(jù)變化不同步，進(jìn)而使兩者同方位同極距的相關(guān)系數(shù)低。

表1 ZD9A-Ⅱ與GEF-Ⅱ相關(guān)系數(shù)計算統(tǒng)計表

為了進(jìn)一步探究該技術(shù)指標(biāo)值的設(shè)置是否合理，采用3臺GEF-Ⅱ儀器記錄2019年5月8日至2019年6月7日比測期間的觀測數(shù)據(jù)按照公式（1）～（2）計算該技術(shù)指標(biāo)值。考慮地電阻率同場地供電對地電場觀測的干擾，選取3臺GEF-Ⅱ儀器記錄的5月10日中0—1時觀測數(shù)據(jù)（圖4），圖中5—14分、34—42分時間段分別是臺上和比測平臺的地電阻率供電時間。從圖4可以看出，受供電干擾的3臺儀器觀測數(shù)據(jù)變化幅度不一致，且遠(yuǎn)大于正常數(shù)據(jù)的變化幅度，形態(tài)也略有差異。GEF-Ⅱ儀器產(chǎn)出的分?jǐn)?shù)據(jù)是由其秒數(shù)據(jù)取均值得到，因GEF-Ⅱ儀未進(jìn)行GPS授時，3臺儀器時間難以確保秒同步。因此在供電干擾時段內(nèi)，3臺儀器存在秒值數(shù)據(jù)上的差異，進(jìn)而表現(xiàn)分?jǐn)?shù)據(jù)供電時段內(nèi)差異，這種已知的供電干擾影響3臺儀器同方位同極距的相關(guān)計算結(jié)果。因此，為了得到更準(zhǔn)確的結(jié)果，地電阻率供電時段數(shù)據(jù)不參與計算。

圖4 多臺GEF-Ⅱ同裝置05月10日0時0—10分觀測數(shù)據(jù)

將3臺GEF-Ⅱ觀測數(shù)據(jù)的相關(guān)系數(shù)計算結(jié)果繪制成圖5，圖中的圓形點(diǎn)為儀器1和儀器2同方位、同極距計算結(jié)果，其中儀器之間六測道的日相關(guān)系數(shù)大于0.9的個數(shù)為183，符合率為98.39%；圖中三角形點(diǎn)為儀器1和儀器3同方位、同極距計算結(jié)果，其中儀器之間六測道的日相關(guān)系數(shù)大于0.9的個數(shù)為179，符合率為96.24%；圖中正方形點(diǎn)為儀器2和儀器3同方位、同極距計算結(jié)果，其中儀器之間六測道的日相關(guān)系數(shù)大于0.9的個數(shù)為186，符合率為100%。3臺GEF-Ⅱ儀器之間的同裝置觀測結(jié)果日相關(guān)系數(shù)滿足于設(shè)定的技術(shù)指標(biāo)值。

圖5 多臺GEF-Ⅱ同裝置日相關(guān)系數(shù)

根據(jù)ZD9A-Ⅱ儀器與GEF-Ⅱ儀器觀測數(shù)據(jù)、3臺GEF-Ⅱ儀器觀測數(shù)據(jù)的計算結(jié)果，分析認(rèn)為前者所得計算結(jié)果不滿足于設(shè)置的技術(shù)指標(biāo)值的原因在于兩種儀器之間的時間差。

2.1.2 儀器之間同裝置觀測數(shù)據(jù)差值絕對值的日均值ΔX

連接同一觀測裝置進(jìn)行比測的觀測儀器，由于儀器本身存在的儀器差的原因，不同儀器產(chǎn)出同裝置（同方位、同極距）的觀測數(shù)據(jù)做差值分析，差值應(yīng)在相對穩(wěn)定的范圍內(nèi)變化。因此設(shè)立儀器之間同裝置觀測數(shù)據(jù)差值絕對值的日均值，并給出該技術(shù)指標(biāo)值和符合率，從差值角度來衡量多臺儀器觀測數(shù)據(jù)變化趨勢，確定儀器是否存在漂移現(xiàn)象。參考《地震監(jiān)測專業(yè)設(shè)備定型測試大綱》中臺站比測地電場儀器之間的同測道觀測數(shù)據(jù)日差值技術(shù)指標(biāo)值，給出ΔX技術(shù)指標(biāo)值：應(yīng)不大于0.5 mV/km，符合率C應(yīng)大于95%，ΔX按公式（4）～（5）計算。

式中：xi、yi為不同儀器的觀測數(shù)分鐘值序列；n為參與計算的數(shù)據(jù)個數(shù)，不大于1 440；NΔX、nΔX分別為ΔX的總個數(shù)、ΔX超過0.5 mV/km的總個數(shù)。

將ZD9A-Ⅱ、GEF-Ⅱ觀測數(shù)據(jù)，舍棄地電阻率供電時段的干擾數(shù)據(jù)，按照公式（4）進(jìn)行計算，得到兩套儀器各個方位長極距與長極距、短極距與短極距之間的觀測數(shù)據(jù)差值絕對值的日均值（圖6）。統(tǒng)計6個測道的計算結(jié)果，兩種儀器之間的同裝置觀測數(shù)據(jù)差值絕對值的日均值大于0.5的個數(shù)為0，符合率為100%，滿足于設(shè)定的技術(shù)指標(biāo)值。

圖6 ZD9A-Ⅱ與GEF-Ⅱ同方位同極距差值絕對值的日均值

2.1.3 同一臺比測儀器同一方位不同裝置觀測結(jié)果差值的月均值D月

在理論上，局部區(qū)域內(nèi)地電場是基本穩(wěn)定的，同一方位不同測道觀測數(shù)據(jù)的變化幅度應(yīng)該是一致的[10]，這體現(xiàn)觀測數(shù)據(jù)的內(nèi)在質(zhì)量及觀測儀器的穩(wěn)定性。同時查看參考儀器的運(yùn)行穩(wěn)定程度，設(shè)置差值月均值技術(shù)指標(biāo)值。參照現(xiàn)行地電場觀測資料質(zhì)量評比辦法，對同臺儀器產(chǎn)出的觀測數(shù)據(jù)計算其同方位的長、短極距的差值月均值，D月按公式（6）～（7）計算，D月技術(shù)指標(biāo)值應(yīng)不大于2 mV/km。

式中：Xi,Yj分別為同一觀測方向2個不同裝置的地電場分鐘觀測值；分別為當(dāng)日0時0分至3時59分的觀測值的平均值；n為一天中參與差值計算的地電場分鐘值個數(shù)；N為1月中參與月平均值計算的天數(shù)。

利用ZD9A-Ⅱ與GEF-Ⅱ儀器6個月觀測數(shù)據(jù)按照公式（4）～（5）分別計算，得到其同儀器、同方位、不同裝置的差值月均值（表2）。參考儀器ZD9A-Ⅱ觀測數(shù)據(jù)所有方位的差值月均值均小于2 mV/km，這說明觀測場地和參考儀器運(yùn)行穩(wěn)定，GEF-Ⅱ儀器6個月觀測數(shù)據(jù)所有方位的差值月均值也均小于2 mV/km，兩者的計算結(jié)果均滿足該技術(shù)指標(biāo)設(shè)定值。

表2 同儀器同方位不同裝置的差值月均值統(tǒng)計表

2.1.4 同一臺比測儀器同方位不同極距觀測結(jié)果的日相關(guān)系數(shù)r1日

假定同一測量場地的地下介質(zhì)電性結(jié)構(gòu)是均勻的，地電場布設(shè)方式為L型，同一方位分為長、短極距兩道觀測，則在同方位不同極距正常觀測得到的地電場遠(yuǎn)場變化形態(tài)應(yīng)具有一致性[11-13]。因而，設(shè)置同方位不同極距觀測結(jié)果的日相關(guān)系數(shù)技術(shù)指標(biāo)和符合率，即計算同臺儀器同方位長、短極距的日相關(guān)系數(shù)，衡量儀器產(chǎn)出的觀測數(shù)據(jù)內(nèi)在質(zhì)量，同時也驗證參考儀器和觀測場地的質(zhì)量。參照儀器之間同裝置觀測結(jié)果的日相關(guān)系數(shù)技術(shù)指標(biāo)值，r1日技 術(shù)指標(biāo)值為：應(yīng)不小于0.9，符合率C應(yīng)大于95%，按照公式（1）～（3）計算，其中yi、xi為同儀器同方位長、短極距測道的觀測數(shù)據(jù)。

利用ZD9A-Ⅱ與GEF-Ⅱ儀器6個月的觀測數(shù)據(jù)，得到其同儀器、同方位、不同極距的日相關(guān)系數(shù)和符合率，統(tǒng)計計算結(jié)果見圖7。ZD9A-Ⅱ所有方位長、短極距之間的日相關(guān)系數(shù)大于0.9的個數(shù)為533，符合率96.56%；GEF-Ⅱ所有方位長、短極距之間的日相關(guān)系數(shù)大于0.9的個數(shù)為551，符合率99.82%。兩者求得結(jié)果均滿足設(shè)置的技術(shù)指標(biāo)值。

圖7 ZD9A-Ⅱ、GEF-Ⅱ同方位不同極距的日相關(guān)系數(shù)

2.2 地電阻率技術(shù)指標(biāo)值

2.2.1 小時觀測值3倍相對均方差的日均值K(σn1)日

地電阻率觀測數(shù)據(jù)為小時觀測值，每次觀測多個數(shù)值，取其平均作為該小時的觀測值，同時產(chǎn)出該次測量的均方差值。目前ZD8M地電阻率儀中產(chǎn)出的均方差已乘以3，即小時觀測值3倍均方差。在此基礎(chǔ)上求取其相對均方差，取其平均值，得到小時觀測值3倍相對均方差的日均值用來分析儀器的穩(wěn)定性，同時利用參考儀器的結(jié)果來確保觀測場地的穩(wěn)定。將小時觀測值3倍相對均方差的日均值作為技術(shù)指標(biāo)，指標(biāo)值設(shè)置為應(yīng)不大于0.3%，以此來分析參考儀器和比測儀器每次觀測內(nèi)所得到的數(shù)據(jù)的離散程度。用參考儀器的結(jié)果確保觀測場地穩(wěn)定的基礎(chǔ)下，查看比測儀器的運(yùn)行穩(wěn)定程度。引用《地震監(jiān)測專業(yè)設(shè)備定型測試大綱》中臺站比測地電阻率儀器的小時觀測值3倍相對均方差的日均值的技術(shù)指標(biāo)值，給出K(σn1)日技術(shù)指標(biāo)值：應(yīng)不大于0.3%。K(σn1)日按公式（8）～（9）計算。

式中：n為 每天觀測次數(shù)，一般為24；K(σn1)i為第i小時觀測值的3倍相對均方差。

將參考儀器和比測儀器小時觀測值3倍相對均方差的日均值計算結(jié)果繪制成圖8。從圖中可看出，參考儀器和比測儀器的NS向計算結(jié)果和EW向計算結(jié)果在小范圍內(nèi)波動，處于0.07%以下。表明兩者產(chǎn)出觀測數(shù)據(jù)穩(wěn)定，同時兩者的計算結(jié)果均小于0.3%，滿足設(shè)置的技術(shù)指標(biāo)值。

圖8 小時觀測值3倍相對均方差的日均值

2.2.2 觀測值日均值的3倍相對實驗標(biāo)準(zhǔn)偏差K日

比測臺站的環(huán)境需穩(wěn)定、優(yōu)良，在此基礎(chǔ)上建設(shè)地電臺站比測觀測裝置，再經(jīng)過相關(guān)認(rèn)證和考核，形成比測平臺，因而比測平臺的觀測環(huán)境是處于相對穩(wěn)定的狀態(tài)。在相對穩(wěn)定狀態(tài)下的觀測場地內(nèi)，地電阻率儀產(chǎn)出的觀測數(shù)據(jù)也應(yīng)在一定范圍內(nèi)出現(xiàn)波動。因此，設(shè)立觀測值日均值的3倍相對實驗標(biāo)準(zhǔn)偏差，以考察比測儀器觀測的穩(wěn)定性，同時利用參考儀器的該技術(shù)指標(biāo)值來體現(xiàn)觀測場地的環(huán)境。觀測值日均值的3倍相對實驗標(biāo)準(zhǔn)偏差K日按 公式（10）計算，K日應(yīng)不大于0.3%。

式中：n為每天觀測次數(shù)，一般為24；ρsi為第i小時ρs觀測值；是ρs日均值。

將兩者觀測數(shù)據(jù)按照公式（10）計算，結(jié)果繪制成圖9。從圖中可見，參考儀器和比測儀器的NS向觀測值日均值3倍相對實驗標(biāo)準(zhǔn)偏差處于0.05%以下，EW向觀測值日均值3倍相對實驗標(biāo)準(zhǔn)偏差處于0.04%以下。兩者儀器所得結(jié)果皆小于0.3%，滿足于該技術(shù)指標(biāo)的設(shè)置值。

圖9 觀測值日均值的3倍相對實驗標(biāo)準(zhǔn)偏差

2.2.3 月觀測精度kσn

在儀器長時間觀測期間內(nèi)，為確定日與日之間的觀測數(shù)據(jù)波動是否能在穩(wěn)定的范圍內(nèi)，因此以日為單位，求取參考儀器的月觀測精度，作為觀測場地的穩(wěn)定性依據(jù)；再求取比測儀器各個測道的月觀測精度，衡量比測儀器的穩(wěn)定性。kσn按公式（11）計算，kσn應(yīng)不大于0.3%。

式中：n為 每月觀測的天數(shù)；ρsi為 每天ρs的日均值；的ρs月均值。

兩者觀測數(shù)據(jù)的月觀測精度計算結(jié)果見表3，可看出參考儀器和比測儀器的NS向和EW向月觀測精度都小于0.1%，滿足于該技術(shù)指標(biāo)值設(shè)置值0.3%。

表3 月觀測精度統(tǒng)計表

兩臺儀器連接相同的觀測裝置，根據(jù)地電阻率觀測原理，產(chǎn)出的數(shù)據(jù)應(yīng)一致。但由于兩臺儀器的儀器差存在，其兩者的觀測數(shù)據(jù)應(yīng)在一定范圍內(nèi)波動，因此設(shè)立儀器之間觀測差值絕對值的平均值，其應(yīng)不大于（0.2%測值+0.04 Ω·m），符合率C應(yīng)大于95%，以此來考察比測儀器是否存在漂移。兩臺儀器中一臺為平?jīng)雠_現(xiàn)行儀器，其產(chǎn)出數(shù)據(jù)在全國地電資料評比中多次獲得佳績，因此其產(chǎn)出的數(shù)據(jù)可靠性較高，故將其作為參考儀器與另一臺儀器對比分析。引用《地震監(jiān)測專業(yè)設(shè)備定型測試大綱》中臺站比測地電阻率儀器中被測儀器和參考儀器觀測數(shù)據(jù)日均值差值的技術(shù)指標(biāo)值，各個測道的和符合率公式（12）～（14）計算，應(yīng)不大于（0.2%測值+0.04 Ω·m），符合率C應(yīng)大于95%。式中：n為 一天內(nèi)比測儀器觀測次數(shù)；ρsij、ρsik為第j、k臺儀器第i次的觀測值；Δρsijk為為第i次比測儀器j和k觀測值的差值；分別為的總個數(shù)、未達(dá)到要求（0.2%測值+0.04 Ω·m）的總個數(shù)。

將參考儀器的同時段觀測值的平均值作為測值得到具體的技術(shù)指標(biāo)值，然后得到參考儀器與比測儀器每天的觀測數(shù)據(jù)差值絕對值的平均值，兩者進(jìn)行比較（圖10）。從圖中可見，兩者NS向和EW向的差值絕對值的平均值遠(yuǎn)遠(yuǎn)小于其技術(shù)指標(biāo)設(shè)置值。

圖10 儀器之間觀測值差值絕對值的平均值

3 結(jié)果

按照地電觀測儀器比測技術(shù)規(guī)范中的計算方法，對地電場、地電阻率觀測數(shù)據(jù)進(jìn)行計算，并與其中設(shè)置的技術(shù)指標(biāo)值比較，得到以下結(jié)論：

1）從儀器之間同裝置的數(shù)據(jù)形態(tài)和差值、儀器的長短極距的數(shù)據(jù)形態(tài)和差值考慮，設(shè)立地電場儀4個技術(shù)指標(biāo)和技術(shù)指標(biāo)值，并給出計算方法。利用ZD9A-Ⅱ、GEF-Ⅱ地電場儀和多臺GEF-Ⅱ地電場儀的觀測數(shù)據(jù)，按照計算方法得到4個地電場儀技術(shù)指標(biāo)值，與設(shè)定的技術(shù)指標(biāo)值對比，其滿足于地電臺站比測技術(shù)規(guī)范中的設(shè)置值，分析認(rèn)為地電臺站比測技術(shù)規(guī)范中地電場的技術(shù)指標(biāo)值設(shè)置合理，能夠用于地電場儀臺站比測，能夠滿足于現(xiàn)今地電場儀器的觀測水平。

2）從儀器之間同裝置的數(shù)據(jù)差值、儀器的觀測精度考慮，設(shè)立地電阻率儀4個技術(shù)指標(biāo)和技術(shù)指標(biāo)值，并給出具體的計算方法。利用平?jīng)霰葴y平臺兩臺ZD8M地電阻率儀的觀測數(shù)據(jù)，按照計算方法得到4個地電阻率儀技術(shù)指標(biāo)值，與設(shè)定的技術(shù)指標(biāo)值對比，其滿足于地電臺站比測技術(shù)規(guī)范中的設(shè)置值，分析認(rèn)為地電臺站比測技術(shù)規(guī)范中地電阻率儀的技術(shù)指標(biāo)設(shè)置合理，能夠用于地電阻率儀臺站比測，能夠滿足于現(xiàn)今地電阻率儀器的觀測水平。

4 分析與討論

1）在計算儀器之間同裝置觀測結(jié)果每天的日相關(guān)系數(shù)時，ZD9A-Ⅱ、GEF-Ⅱ地電場儀所得到的結(jié)果不滿足于該技術(shù)指標(biāo)設(shè)置值，查看其觀測數(shù)據(jù)兩者相位存在滯后，分析認(rèn)為兩臺地電場儀時間不同步是導(dǎo)致計算結(jié)果不佳的主要原因，也可能是ZD9A-Ⅱ、GEF-Ⅱ地電場儀的工作原理、儀器結(jié)構(gòu)的差異造成，如頻率范圍不同，這有待后續(xù)進(jìn)一步研究。在臺站比測中多臺地電場儀之間同裝置的日相關(guān)系數(shù)是基于觀測數(shù)據(jù)形態(tài)的高度一致性，需要多臺地電場儀觀測時間同步，只有在確保時間一致的前提下，儀器之間所得數(shù)據(jù)相位才會保持一致。在以后的臺站比測中，參考儀器應(yīng)是經(jīng)實驗室的定期檢測，或是參考儀器和比測儀器都增加GPS對時模塊，確保時間的準(zhǔn)確性和一致性。

2）通過對現(xiàn)有所得的資料分析，地電臺站比測技術(shù)規(guī)范中地電場儀和地電阻率儀的技術(shù)指標(biāo)值設(shè)置合理，能為地電觀測儀器技術(shù)比測規(guī)范技術(shù)指標(biāo)的編寫和地電觀測儀器臺站比測提供依據(jù)與支撐，從而為規(guī)范地電觀測儀器臺站比測和產(chǎn)出用于地震監(jiān)測、預(yù)報的高質(zhì)量地電觀測數(shù)據(jù)提供保障。由于試驗條件所限，用于分析的觀測數(shù)據(jù)樣本不多，所提出的技術(shù)指標(biāo)值可能存在一定的局限性。在后續(xù)研究中，還需多研究更多的技術(shù)指標(biāo)和其對應(yīng)值，以此來更全面地衡量比測儀器。

致謝感謝應(yīng)急管理部國家自然災(zāi)害防治研究院胡哲高級工程師的對本文的指導(dǎo)，感謝甘肅省平?jīng)雠_提供本文中所用的地電阻率觀測數(shù)據(jù)。