北京百善深井鉆孔應變儀觀測資料可靠性分析*

吳立恒 李 濤 陳 征

(中國地震局地殼應力研究所，北京 100085)

北京百善深井鉆孔應變儀觀測資料可靠性分析*

吳立恒 李 濤 陳 征

(中國地震局地殼應力研究所，北京 100085)

通過參考元件穩(wěn)定性，固體潮潮汐因子分析，相互垂直元件計算的面應變相關性三個方面，分析了北京百善深井鉆孔應變觀測系統(tǒng)的可靠性，分析說明北京百善深井鉆孔應變觀測系統(tǒng)具有較高的一致性和可靠性。

深井鉆孔應變儀;可靠性;系統(tǒng)自檢;固體潮調和分析;面應變

1 引言

鉆孔應變儀是揭示中短期(從數秒到數月)地殼連續(xù)變形的理想儀器，尤其在短臨前兆觀測中占有重要位置［1］。因此，鉆孔應變儀的研制與應用得到了國內外同行的高度重視。

對于鉆孔應變儀觀測資料的可靠性分析，不少學者進行過討論。如對同一或相近臺站安裝不同類型的鉆孔應變儀或傾斜儀，對其觀測資料進行對比分析和驗證其資料的可靠性［2］;對分量式鉆孔應變儀采用兩組相互垂直的元件計算的面應變進行自檢［3］;對分量式鉆孔應變儀通過實地絕對標定與實地相對標定，根據標定參數的一致性和穩(wěn)定性來檢驗鉆孔應變儀觀測資料的可靠性等［4，5］。

北京百善觀測點于2008年10月安裝了RZB多分量鉆孔應變儀。此應變儀由兩個4分量探頭(8個不同方向的工作元件)集成而成，同時裝有參考元件。本文將利用該井近3年的觀測數據對其深井鉆孔應變觀測系統(tǒng)進行可靠性分析。

2 北京百善RZB深井應變

北京百善RZB井位于昌平區(qū)百善鎮(zhèn)，應變儀安裝深度為210 m。此應變儀由兩個4分量的RZB鉆孔應變儀構成，其元件安裝方向如表1所示，探頭與巖石的接觸采用微膨脹水泥使其耦合。該儀器于2008年10月安裝成功并開始連續(xù)記錄。

表1 北京百善RZB深井多分量鉆孔應變儀元件方位Tab.1 The component orientation of RZB deep strain meter in Baishan，Beijing

3 參考元件的穩(wěn)定性分析

圖1是2009年2月各工作元件的應變記錄曲線。圖1表明8個元件都能清晰地記錄到固體潮汐信號以及趨勢性應變變化。

5號元件是RZB鉆孔應變儀研制時設計的一個獨有元件，它作為參考元件放置在探頭內，主要考察在深井環(huán)境下工作的傳感器自身的長期穩(wěn)定性和測量系統(tǒng)的長期穩(wěn)定性。對5號元件的分析結果表明:該元件的年變化量為3×10-9量級(表2)，說明該儀器具有較好的長期穩(wěn)定性。

對5號元件記錄數據的統(tǒng)計分析結果為:觀測數據符合N(37.34×10-8，(3.2×10-10)2)正態(tài)分布，在2σ置信區(qū)間為(37.28×10-8，37.40× 10-8)，置信概率為95.44%;系統(tǒng)噪聲水平在10-10應變量級。

表2 北京百善深井多分量觀測系統(tǒng)參考元件年變化量(單位：10-8)Tab.2 Annual variation of the reference component(unit：10-8)

4 固體潮調和分析

利用中國地震前兆臺網運行評價系統(tǒng)對2010年9—12月的觀測數據進行調和分析［6］，得到的結果如表3、4所示。

由表3結果可知，除1號元件的固體潮調和分析結果的相對誤差略大外，其他受力元件在2010年9—12月都遠小于0.05，觀測數據內精度高，完全符合地殼形變觀測技術的要求。由表4結果可知，百善深井RZB鉆孔應變觀測系統(tǒng)各元件計算的面應變潮汐因子雖然存在差異，但自身非常穩(wěn)定，相對誤差也非常小，說明整個測量系統(tǒng)具有較高的觀測精度和可靠性。

圖1 2009年2月各元件應變記錄月曲線Fig.1 Monthly curves of the strain records of all the components in February 2009

表3 觀測系統(tǒng)線應變調和分析結果Tab.3 Venedikov results of harmonic analysis of the linestrain

表4 面應變調和分析結果Tab.4 Venedikov results of tidal harmonic analysis of plain train

5 觀測系統(tǒng)自檢

根據鉆孔加襯模型，當探頭和圍巖耦合處于理想狀態(tài)時，鉆孔θ方向的孔徑相對變化為［7，8］:

如果將兩個相互垂直的孔徑相對變化的測值相加，有

通過式(2)可以對觀測數據進行自檢。百善深井鉆孔應變儀安裝上相互垂直的元件有1號和3號，2號和4號，6號和8號，7號和9號，可以通過四組垂直元件計算的面應變的結果一致性，說明探頭和圍巖的耦合情況，考察記錄資料的可靠。用此方法，對2010年12月北京百善數據進行垂直元件相加求面應變處理，曲線如圖2所示，由此可見;通過4組不同垂直元件所得面應變結果，在變化趨勢上具有較好的一致性。此外，從圖2上我們能很直觀的發(fā)現(xiàn)，其中運用1號、3號和2號、4號元件的數據分別計算的兩組面應變結果相似度很高。因此使用此兩組計算的面應變結果會更可靠，更有分析研究價值。在此基礎上采用排列組合的方式進行各組相關性分析，相關系數如表5所示。從表5可知，北京深井鉆孔觀測系統(tǒng)耦合優(yōu)良，數據可靠。

圖2 2010年12月數據曲線圖Fig.2 Data curves in December 2010

一般的鉆孔應變儀只有4個受力元件，將1號+3號和2號+4號進行自檢時，只能得到一個相關性系數，無法判別是哪組元件工作狀態(tài)更優(yōu)、測值更真實可靠。但北京百善深井鉆孔應變觀測系統(tǒng)有兩個探頭，可以通過兩兩組合進行相關性分析，以驗證哪組元件的測值更真實可靠，驗證結果如表5所示。從表5可以發(fā)現(xiàn)，8個元件的自檢相關性系數都很高，說明該系統(tǒng)2010年12月的觀測數據質量高。此外，我們認為自檢相關性系數的結果可進一步分析各元件工作狀態(tài)的優(yōu)劣，對于以后數據的分析利用提供權重的參考信息，有利于科學的運用這些數據。從表5中我們發(fā)現(xiàn)7號、9號元件的自檢相關性系數略低于其他組，那么我們認為，在使用7號，9號元件的測量數據進行分析時，可適當降低其信任度。而1號、2號、3號、4號、6號、8號工作狀態(tài)優(yōu)異，測數真實可靠，更適合進一步的數據分析運用。

表5 元件自檢相關性系數Tab.5 Component self-correlation coefficient

6 結論

1)百善深井鉆孔應變儀參考元件的年漂移為3 ×10-9，且年漂移量穩(wěn)定。通過噪聲分析說明此系統(tǒng)的噪聲水平在10-10應變量級。

2)系統(tǒng)線應變、面應變固體潮調和分析說明百善深井鉆孔應變觀測系統(tǒng)具有較高的一致性和可靠性。

3)系統(tǒng)自檢的相關系數在(0.863，0.957)區(qū)間，自檢效果良好，說明百善深井鉆孔應變觀測系統(tǒng)觀測可靠。

1 邱澤華，石耀霖國外鉆孔應變觀測的發(fā)展現(xiàn)狀［J］.地震學報，2004，26(增):481-488.(Qiu Zehua and Shi Yaolin.Developments of borehole strain observation outside China［J］.Acta Seismlogica Sinica，2004，26(Supp.):481 -488)

2 張國紅，等.鉆孔應變資料的可靠性分析［J］.大地測量與地球動力學，2010，(增):62-66.(Zhang Guohong，et al. Reliability analysis of borehole strainmetry observations［J］.Journal of Geodesy and Geodynamics，2010，(Supp.):62-66)

3 李濤，等.RZB型鉆孔應變儀在青藏高原東緣地應力監(jiān)測中的應用［J］.北京大學學報(自然科學版)，2011，47 (4):677-683.(Li Tao，et al.The application of the RZB borehole strainmeter in crustal stress observation at the eastern margin of the Tibet Plateau［J］.Acta Scientiarum Naturalium Unicersitatis Pekinensis，2011，47(4):677-683)

4 邱澤華，石耀霖，歐陽祖熙.四分量鉆孔應變觀測的實地絕對標定［J］.地震，2005，25(3):27-34.(Qiu Zehua，Shi Yaolin and Ouyang Zuxi.Absolute calibration of 4-component borehole trainmeters in situ［J］.Earthquake，2005，25(3): 27-34)

5 邱澤華，石耀霖，歐陽祖熙.四分量鉆孔應變觀測的實地相對標定［J］.大地測量與地球動力學，2005，(1):118-122.(Qiu Zehua，Shi Yaolin and Ouyang Zuxi.Relative insitu calibration of 4-component borehole strain obversation［J］.Journal of Geodesy and Geodynamics，2005，(1):118-122)

6 蔣駿，等.地震前兆信息處理與軟件系統(tǒng)［M］.北京:地震出版社，2000.(Jiang Jun，et al.Earthquake precursor information disposal and software system［M］.Beijing:Seismological Press，2000)

7 潘立宙.測量地應力的鉆孔變形法［A］.中國地質科學院地質力學研究所和國家地震局地震地質大隊編.地應力測量的原理和應用［C］.北京:地質出版社，1981.(Pan Lizou.Borehole deformation method of ground stress measurement［A］.Principles and Application of Ground Stress Measurement［C］.Beijing:Geological Press，1981)

8 潘立宙.與地應力測量有關的幾個公式的推導和討論［A］.中國地質科學院地質力學研究所和國家地震局地震地質大隊編.地應力測量的原理和應用［C］.北京:地質出版社，1981.(Pan Lizou.Derivation of several formulas with their application to in-situ stress measurements［A］.Principles and Application of Ground Stress Measurement［C］.Beijing:Geological Press，1981)

RELIABILITY ANALYSIS OF DEEP BOREHOLE STRAINMETER DATA IN BEIJING BAISHAN

Wu Liheng，Li Tao and Chen Zheng
(The Institute of Crustal Dynamics，CEA Beijing 100085)

Deep borehole strainmeter was installed in Beijing Baishan in October 2008，it has accumulated a lot of observations which provides the basis for the reliability analysis.We analyzed system reliability from stability of the reference component，harmonic tidal factor，and correlation analysis of the plane strain calculated by perpendicular component.In the conclusion，Beijing Baishan deep borehole strainmeter has high consistency and reliability，it can provide reliable and credible observations for geophysics research and earthquake prediction.

deep borehole strainmeter;reliability;system self-test;tidal harmonic analysis;plane strain

1671-5942(2012)03-0041-04

2012-05-03

中國地震局地殼應力研究所基本科研專項基金(ZDJ2010-31-01，ZDJ2011-03)

吳立恒，女，1980年生，碩士，助理研究員，現(xiàn)主要從事地震前兆儀器傳感器的研制.E-mail:smiling-lily8013@163.com

TH762.2

A