耿杰++張玲

摘要：2015年11月至2016年8月諸城井完成了水溫不同深度對(duì)比觀測(cè)儀器安裝調(diào)試與產(chǎn)出數(shù)據(jù)的試驗(yàn)，在對(duì)產(chǎn)出數(shù)據(jù)進(jìn)行直觀和定性分析的基礎(chǔ)上，對(duì)不同深度、同層位、不同通道的2套地?zé)崴痪C合觀測(cè)儀產(chǎn)出數(shù)據(jù)進(jìn)行了趨勢(shì)、方差和均值的一致性檢驗(yàn)。結(jié)果表明：諸城井不同深度水溫觀測(cè)數(shù)據(jù)與地溫梯度測(cè)試結(jié)果基本一致，顯示水溫產(chǎn)出數(shù)據(jù)是可靠的；經(jīng)一致性檢驗(yàn)，同層位水溫變化趨勢(shì)、方差一致，但均值不一致，表明綜合觀測(cè)儀器產(chǎn)出數(shù)據(jù)間有一定的系統(tǒng)誤差，數(shù)據(jù)變化本底值有一定的差異。

關(guān)鍵詞：諸城井；水溫觀測(cè)；一致性檢驗(yàn)

中圖分類號(hào)：P315723文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1000-0666（2017）02-0226-07

0前言

2015年3月25日，山東省地震局監(jiān)測(cè)預(yù)報(bào)處組織專家對(duì)諸城井流體對(duì)比觀測(cè)實(shí)施方案進(jìn)行了論證，初步確定在諸城井不同深度安裝8個(gè)水溫傳感器和2個(gè)水位傳感器，實(shí)施數(shù)字化水溫對(duì)比觀測(cè)，以期為地震分析預(yù)報(bào)與研究提供有價(jià)值的觀測(cè)數(shù)據(jù)。

按照論證后的實(shí)施方案和國(guó)家相關(guān)規(guī)定，山東省地震局完成了儀器（DRSW-2 地?zé)崴痪C合觀測(cè)儀）招標(biāo)采購，在此基礎(chǔ)上于2015年11月至2016年8月先后完成了設(shè)備安裝調(diào)試、井房及避雷網(wǎng)改造、地溫梯度測(cè)試、儀器設(shè)備安裝、主機(jī)軟件更新升級(jí)、UPS電源及氣象三要素儀器安裝，并在分析不同儀器、不同通道產(chǎn)出數(shù)據(jù)的基礎(chǔ)上，進(jìn)行了儀器主機(jī)工作狀態(tài)及不同深度產(chǎn)出數(shù)據(jù)的試驗(yàn)觀測(cè)。本文對(duì)諸城井安裝完成后儀器所產(chǎn)出的水溫?cái)?shù)據(jù)進(jìn)行了系統(tǒng)分析研究，并對(duì)同層位對(duì)比觀測(cè)的2套水溫儀器產(chǎn)出數(shù)據(jù)，進(jìn)行了趨勢(shì)、方差和均值的一致性檢驗(yàn)。

1井孔及觀測(cè)基本情況

諸城井由勝利油田地質(zhì)鉆井公司第二分公司鉆探完成，于1997 年 12 月 24日開始，至1999 年 2 月 8 日成井，是石油部門為揭露諸城凹陷地層層序、巖性特征，探查諸城凹陷含油氣情況而鉆探的。

諸城井位于諸城市昌城鎮(zhèn)大下泊村東約900 m，構(gòu)造位置處在膠萊坳陷西南端的諸城凹陷辛興鼻狀構(gòu)造帶上，井深5 005 m，上部套管3 138 m，井孔自上而下有新生界第四系平原組，中生界白堊系王氏組、青山組、萊陽組等地層，地層巖性主要是砂巖、礫巖、凝灰?guī)r等。觀測(cè)含水層為中生界白堊系萊陽組。

2007年4月，諸城井開展數(shù)字化地下流體水溫、水位觀測(cè)，采用SZW-1A型數(shù)字式溫度計(jì)進(jìn)行地?zé)嵊^測(cè)。改造前水溫探頭投放深度為600 m，采用LN-3A數(shù)字水位儀進(jìn)行數(shù)字化水位觀測(cè)后，探頭投放深度為30094 m。

自投入觀測(cè)以來，諸城井水溫、水位的長(zhǎng)期動(dòng)態(tài)特征基本清晰，無年動(dòng)態(tài)變化，變化有一定規(guī)律，不受當(dāng)?shù)亟涤旰偷叵滤_采、灌溉影響。諸城井水溫觀測(cè)資料顯示水位變化呈現(xiàn)穩(wěn)定性較好的趨勢(shì)上升變化，在數(shù)次大震時(shí)未發(fā)現(xiàn)同震效應(yīng)。無水位固體潮記錄，記錄水震波能力較差。水位觀測(cè)同震響應(yīng)不顯著，僅對(duì)部分強(qiáng)震有響應(yīng)。諸城井含水層簡(jiǎn)單，干擾因素相對(duì)較少，數(shù)據(jù)變化信度較高，該井地下流體觀測(cè)資料對(duì)震情判定有一定的作用，在震情會(huì)商中得到了有效應(yīng)用。[HJ]

2水溫對(duì)比觀測(cè)的實(shí)施

[KG（0.2mm]按照《地震及前兆數(shù)字觀測(cè)技術(shù)規(guī)范——地下流體觀測(cè)（試行）》（中國(guó)地震局，2001）及《地震地下流體觀測(cè)方法井水位觀測(cè)》（DB/T 48—2012）的相關(guān)規(guī)定，2016年3月5日，山東省地震局對(duì)諸城井進(jìn)行了地溫梯度測(cè)試，起始測(cè)點(diǎn)為井口以下75 m，測(cè)點(diǎn)間距為50 m，每50 m進(jìn)行 30 min的溫度測(cè)量，共測(cè)量了19個(gè)不同層位的水溫值，最深測(cè)試深度為975 m。地溫梯度測(cè)試結(jié)果顯示：諸城井全井段處在正梯度，全井段溫度梯度處在（25～4） ℃/hm區(qū)間；400 m以上井段，溫度梯度振蕩幅度較大，400 m以下，溫度梯度振蕩幅度明顯減?。蝗螠囟忍荻纫?0067 3 ℃/hm的速率下降，即隨著深度的增加，溫度梯度減小。地溫梯度的測(cè)試結(jié)果與該井“十五”數(shù)字化建設(shè)改造地溫梯度測(cè)試結(jié)果一致。

在對(duì)諸城井地溫梯度測(cè)試結(jié)果進(jìn)行分析的基礎(chǔ)上（車用太等，2002；孫亞強(qiáng)等，2001），結(jié)合“諸城井流體對(duì)比觀測(cè)實(shí)施方案”論證討論結(jié)果以及現(xiàn)場(chǎng)工作實(shí)際情況，2016年3月6日至4月19日，山東省地震局完成了2套DRSW-2 地?zé)崴痪C合觀測(cè)儀和RTP-2型氣象三要素儀器設(shè)備安裝與調(diào)試（表1），同時(shí)完成了山特3KVA不間斷電源的安裝。為了保證儀器供電電源的連續(xù)性，2套綜合主機(jī)各接入同一原控。安裝完成后的2套水位、8套水溫及氣象三要素的全部?jī)x器均已接入山東省地震局前兆未參評(píng)儀器服務(wù)器，數(shù)據(jù)開始采集入庫，選取綜合主機(jī)2水位（設(shè)備號(hào)1528）和600 m深度水溫（設(shè)備號(hào)986）作為參評(píng)儀器，接入安丘節(jié)點(diǎn)前兆服務(wù)器。

3深層水溫對(duì)比研究

對(duì)比諸城井水溫?cái)?shù)據(jù)與同深度地溫梯度測(cè)試結(jié)果（圖1），發(fā)現(xiàn)二者測(cè)值相近。不同深度水溫短期日動(dòng)態(tài)均顯示為短周期隨機(jī)波動(dòng)變化（汪成民等，1988；王吉易等，1997），2套水溫觀測(cè)儀器在600 m深度綜合主機(jī)2產(chǎn)出的數(shù)據(jù)與原“十五”投測(cè)觀測(cè)儀器產(chǎn)出數(shù)據(jù)的趨勢(shì)變化形態(tài)和上升速率一致；綜合主機(jī)1在800 m深度產(chǎn)出的觀測(cè)數(shù)據(jù)，同時(shí)顯示為趨勢(shì)下降變化，且趨勢(shì)變化形態(tài)和下降速率基本一致（圖1 a，b）。

[JP2]為了進(jìn)一步考察2套主機(jī)及其通道的工作狀態(tài)，于2016年8月30日將綜合主機(jī)2第二通道的600 m探頭接入主機(jī)1第三通道，將綜主機(jī)1第三通道的800 m探頭接入主機(jī)2第二通道，并完成相關(guān)參數(shù)配置，實(shí)施了主機(jī)1、主機(jī)2不同通道不同層位探頭產(chǎn)出數(shù)據(jù)（圖1 c，d）的進(jìn)一步分析對(duì)比檢驗(yàn)。[JP]

[BT（2]31同層位水溫觀測(cè)產(chǎn)出數(shù)據(jù)的一致性檢驗(yàn)[BT）]

在對(duì)諸城井水位和不同深度多層位水溫產(chǎn)出數(shù)據(jù)進(jìn)行直觀對(duì)比、定性描述的基礎(chǔ)上（蔣駿等，2000；賈化周等，1995），選取2016年7月和9月諸城井深層水溫600 m和800 m兩個(gè)不同層位的水溫產(chǎn)出數(shù)據(jù)，選擇同時(shí)刻的觀測(cè)值（整點(diǎn)值）進(jìn)行對(duì)比計(jì)算，以檢驗(yàn)其同時(shí)刻及較長(zhǎng)時(shí)段觀測(cè)資料的一致性。

對(duì)諸城井600 m和800 m深度同層位并行觀測(cè)的2套綜合主機(jī)產(chǎn)出數(shù)據(jù)進(jìn)行了變化趨勢(shì)、方差及均值一致性的統(tǒng)計(jì)檢驗(yàn)（李慶海，陶本藻，1982；謝中華，2010）。如完全通過檢驗(yàn)，則表明兩者統(tǒng)計(jì)性狀完全一致，資料可進(jìn)行無縫連接使用；若趨勢(shì)及方差一致但均值不一致，則表明2套儀器產(chǎn)出數(shù)據(jù)的本底值不同；若趨勢(shì)一致但均值及方差不一致，則表明2套儀器產(chǎn)出數(shù)據(jù)相關(guān)性較好，即數(shù)據(jù)變化形態(tài)一致性較好，但數(shù)據(jù)本底值、變化幅度等有所不同；若三者皆不一致，則須對(duì)儀器產(chǎn)出數(shù)據(jù)的可靠性作進(jìn)一步分析和檢測(cè)。

32同層位水溫一致性檢驗(yàn)方法

（1）變化趨勢(shì)一致性檢驗(yàn)——相關(guān)系數(shù)檢驗(yàn)法

令xi及yi（i=1，2，…，n，n為統(tǒng)計(jì)天數(shù)）分別表示同深度2套儀器相同時(shí)刻的觀測(cè)值，則相關(guān)系數(shù)r為：[KH*2/3D]

[JP4]r[KG-*4]=[SX（]n∑[DD（]ni=1[DD）]xiyi-∑[DD（]ni=1[DD）]xi∑[DD（]ni=1[DD）]yi[KG-*3][KF（][JB<2*[]n[KG-*2]∑[DD（]ni=1[DD）]y2i[KG-*3]-[KG-*3][JB<2*（]∑[DD（]ni=1[DD）]yi[JB>2*）]2[KG-*2][JB>2*]][KG-*3]·[KG-*2][JB<2*[]n[KG-*2]∑[DD（]ni=1[DD）]x2i[KG-*3]-[KG-*3][JB<2*（]∑[DD（]ni=1[DD）]xi[JB>2*）]2[KG-*2][JB>2*]][KF）][SX）] [JY]（1）[JP][KH*4/3]

定性地看，r值越趨近于1，則二者的變化趨勢(shì)越一致，取定顯著性水平α=001或α=005或α=01，自由度f=n-2，進(jìn)而可確定相應(yīng)于一定顯著性水平α及自由度f的最小相關(guān)系數(shù)r0，若r>r0則檢驗(yàn)通過，否則不通過。

（2）樣本方差的一致性檢驗(yàn)——F檢驗(yàn)

xi、yi分別為同深度2套儀器相同時(shí)刻的觀測(cè)值，令其長(zhǎng)度分別為nx及ny，統(tǒng)計(jì)量F定義為：[KH*3/4]

F=[SX（]S2xS2y[SX）]-F（nx-1，ny-1）[JY]（2）[KH*3/4D]

其中：S2x、S2y分別為2組數(shù)據(jù)的樣本方差。建立統(tǒng)計(jì)量F時(shí)已認(rèn)定分子大于等于分母，即始終有F≥1。取定顯著性水平α=001或α=005或α=01，可得相應(yīng)于顯著水平為α、第一自由度為（nx-1）、[JP3]第二自由度為（ny-1）所對(duì)應(yīng)的Fα（nx-1，ny-1）；[JP]若F>Fα（nx-1，ny-1），則否定一致性假設(shè)；若F

若前述方差一致性檢驗(yàn)得以通過，則可看作等精度觀測(cè)，否則屬不等精度觀測(cè)。

（3）均值的一致性檢驗(yàn)——t檢驗(yàn)

均值一致性檢驗(yàn)統(tǒng)計(jì)量t定義為：[KH*1]

t=[SX（]（[AKx-D]-[AKy-D]）[KF（][SX（]S2xnx[SX）]-[SX（]S2yny[SX）][KF）][SX）] [JY]（3）[KH*1D]

取定顯著性水平α=001或α=005或α=01，自由度f=nx+ny-2，得t分布臨界值tα（nx+ny-2），使得P{[JB<1|]t[JB>1|]≥tα（nx+ny-2）}=α。進(jìn)而可進(jìn)行檢驗(yàn)：若[JB<1|]t[JB>1|]>tα（nx+ny-2），則檢驗(yàn)不能通過，二者均值不一致；若[JB<1|]t[JB>1|]

33同層位水溫一致性檢驗(yàn)結(jié)果

對(duì)2016年7月和9月諸城井2套綜合觀測(cè)主機(jī)產(chǎn)出的600 m和800 m深度水溫同時(shí)段數(shù)據(jù)進(jìn)行統(tǒng)計(jì)檢驗(yàn)，結(jié)果見表2，以上對(duì)比檢驗(yàn)取顯著性水平α=005（耿杰等，2001；張煒等，1988）。計(jì)算結(jié)果表明，同深度不同主機(jī)產(chǎn)出的水溫?cái)?shù)據(jù)有較一致的變化趨勢(shì)和方差，但2個(gè)層位不同時(shí)段t檢驗(yàn)均未通過，顯示2個(gè)層位的2套不同主機(jī)及不同通道的產(chǎn)出數(shù)據(jù)均值不一致，即不同主機(jī)或觀測(cè)系統(tǒng)間存在一定的系統(tǒng)誤差。

[BT1+*2/3]4結(jié)論

綜合上述對(duì)諸城深層水溫對(duì)比觀測(cè)產(chǎn)出數(shù)據(jù)的研究與分析結(jié)果，得出以下分析結(jié)論：

（1）諸城井不同深度水溫觀測(cè)產(chǎn)出數(shù)據(jù)與地溫梯度測(cè)試結(jié)果基本一致，顯示不同層位水溫產(chǎn)出數(shù)據(jù)是可靠的，產(chǎn)出資料可用于地震分析預(yù)報(bào)。

（2）諸城井同層位水溫對(duì)比觀測(cè)資料的一致性檢驗(yàn)結(jié)果顯示，同層位水溫變化趨勢(shì)、方差一致，但均值不一致，表明2套數(shù)字化綜合觀測(cè)儀器產(chǎn)出數(shù)據(jù)間有一定的系統(tǒng)誤差，其變化本底值、變化幅度有一定的差異。

（3）開展諸城井同型號(hào)儀器多層位水溫對(duì)比觀測(cè)為分析水溫異常變化提供了可靠的依據(jù)，同井不同深度的水溫觀測(cè)為分析不同層位之間水溫影響如水溫的對(duì)流、彌散、傳導(dǎo)的影響方式和影響程度提供了可靠的第一手觀測(cè)數(shù)據(jù)。

（4）對(duì)比諸城井2套觀測(cè)主機(jī)產(chǎn)出水溫?cái)?shù)據(jù)，發(fā)現(xiàn)不同主機(jī)采集通道不同，對(duì)產(chǎn)出數(shù)據(jù)有一定的影響，有待協(xié)商儀器生產(chǎn)廠家進(jìn)行進(jìn)一步檢查和調(diào)試。[KH*1D]

參考文獻(xiàn)：

車用太，孔令昌，陳華靜，等2002地下流體數(shù)字化觀測(cè)技術(shù)[M].北京：地震出版社

耿杰，任桂香，王忠民，等2001山東數(shù)字化水位觀測(cè)資料與模擬觀測(cè)資料的對(duì)比分析[J].地震地質(zhì)，23（3）：455-462

賈化周，張煒，董守玉，等1995地震地下水手冊(cè)[M].北京：地震出版社

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

李慶海，陶本藻1982概率統(tǒng)計(jì)原理和在測(cè)量中的應(yīng)用[M].北京：測(cè)繪出版社，109-268

孫亞強(qiáng)，宋彥云，王方建，等2001地震前兆臺(tái)網(wǎng)技術(shù)—山東數(shù)字試驗(yàn)臺(tái)網(wǎng)[M].北京：地震出版社

汪成民，車用太，萬迪坤，等1988地下水微動(dòng)態(tài)研究[M].北京：地震出版社

王吉易，董守玉，陳建民，等1997地下流體地震預(yù)報(bào)方法[M].北京：地震出版社

謝中華2010MATLAB統(tǒng)計(jì)分析與應(yīng)用[M].北京：航天航空大學(xué)出版社，144-156

張煒，王吉易，鄂秀滿，等1988水文地球化學(xué)預(yù)報(bào)地震原理與方法[M].北京：教育科學(xué)出版社，215-249

中國(guó)地震局2001地震及前兆數(shù)字觀測(cè)技術(shù)規(guī)范——地下流體觀測(cè)（試行）[M].北京：地震出版社

DB/T 48—2012，地震地下流體觀測(cè)方法井水位觀測(cè)[S].[ZK）][HJ][FL）][LM]

[STHZ][WT4HZ][JZ]Research on Comparative Observation of Water[JZ]Temperature in Zhucheng Well in Shandong

[WT5B1][STBZ][JZ]GENG Jie，ZHANG Ling

[WT5"B1X][JZ]（Earthquake Administration of Shandong Province，Jinan 250014，Shandong，China）

[WT5HZ][JZ]Abstract[WTB1]

From November 2015 to August 2016，we have completed the installation and debugging of contrast observation equipment as well as the test of output data at water temperature of different depth in Zhucheng well Based on the intuitive and qualitative analysis of the output data，we have tested the consistency of the trend，variance and mean of output data from two sets of different channels at water temperature of different depth of the same layer in the well The observation data of water temperature in different depth of the well is basically consistent with the result of geothermal gradient，which shows that the output data of water temperature is reliable The results of consistency test with the observation data of water temperature at contrast layers in the well confirmed the same horizon temperature change trend and variance，but the mean is not same，which shows that the output data from comprehensive observation instruments showing the system error and that the background value of data change is certain different

[WTHZ]Keywords：[WTB1]Zhucheng well；water temperature observation；consistency check