ABSTRACT

The borehole inclinometer in Yanshou has been running well since its official observation, and obtained continuous, complete and high-quality data. The operating status of the Yanshou drilling tilt was evaluated. Based on past experience, the main factors of disturbing the borehole inclinometer in Yanshou County seismic station are instrument drift, voltage sudden change, simultaneous-earthquake, atmospheric pressure, rainfall, capacity change of reservoir etc., through analyzing related data during the period of 2005-2015. The results show that the main factors affecting the instrument’s curve shape and observation quality are zero drift of the instrument, and the change of reservoir capacity and short-term rainfall have little effect on borehole tilt observation.

Keywords:Borehole Inclinometer, Interference Analysis, Load Calculation, Capacity Change of Reservoir

延寿地震台钻孔倾斜干扰因素分析

孙强¹，王强¹，李继业²，孙伟¹，贾乃圣元¹，朱家明³

¹哈尔滨市阿城地震台，黑龙江 哈尔滨

²黑龙江省地震局，黑龙江 哈尔滨

³延寿地震台，黑龙江 哈尔滨

收稿日期：2018年7月3日；录用日期：2018年7月20日；发布日期：2018年7月27日

摘 要

延寿地震台钻孔倾斜自正式观测以来运行良好，取得了连续完整且高质量的数据资料。通过对延寿钻孔倾斜的运行状态进行评价，并结合以往的研究结果，归纳总结出可对延寿地震台竖直摆倾斜观测造成干扰的主要因素为仪器零漂干扰、电压突变干扰、同震干扰、气压干扰、降水干扰和水库库容变化干扰等。通过对2005~2015年期间延寿钻孔倾斜及相关资料进行数据分析，结果表明影响仪器曲线形态和观测质量的主要因素为仪器零漂，水库库容变化、短期内降雨对钻孔倾斜观测影响不大。

关键词 :钻孔倾斜，干扰分析，荷载计算，水库库容变化

Copyright © 2018 by authors and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

1. 引言

延寿地震台钻孔倾斜于2004年12月正式开始观测，观测仪器为CZB-1型竖直摆倾斜仪。CZB-1型竖直摆倾斜仪是高精度地倾斜观测仪器，分辨率达10⁻¹⁰ m，灵敏度高，具有较好的应震能力。延寿台钻孔倾斜经过十多年的观测，期间出现多次转折、曲线斜率加速等异常变化特征，且异常出现之后基本都有近震随与之对应。但观测多年以来，同样发现竖直摆倾斜仪受一些干扰因素的影响，如人为引起的干扰，仪器零漂干扰，周围环境和天气干扰等 [1] 。通常来说，地壳形变观测资料主要受降雨、地下水、气压、温度、海潮、水库蓄水等因素的影响，这些因素造成的干扰程度会在不同的环境下有不同的影响程度，在某些观测环境下某种因素占主导地位，在另一种观测环境下可能只是次要因素而已 [2] [3] 。

2. 台址概况

延寿地震台位于哈尔滨市延寿县长寿乡，地理位置位于松嫩平原和东部山区之间的丘陵地带，属于张广才岭余脉。台址周围岩性主要以花岗岩类为主，台基岩性是完整的燕山期花岗岩。该台大地构造位置属兴安岭–内蒙古地槽褶皱、伊春–延寿地槽褶皱北段。观测井位于依兰–舒兰断裂带的东部，依兰–舒兰断裂带贯穿整个东北，由北向南延伸，经过延寿县，属挤压型断裂带。延寿县西部还发育有北西向断裂和北北东向铁力–尚志断裂 [4] 。

3. 仪器工作状态

延寿地震台钻孔倾斜观测仪器为河南省地震局研制的CZB-1型竖直摆倾斜仪(图1)，该观测井井深71.6 m，井孔直径152 mm，传感器距离井口距离为61 m，仪器输出信号线性范围≤5，数字记录的范围为±2 V，测量精度1‰ ms，分辨率10⁻¹⁰ m。

M2波潮汐因子中误差是一项用来评定固体潮观测内在质量优劣的定量指标。中误差可以反映数据观测质量和观测环境的影响，同时观测环境的变化也可以影响中误差，对异常的判别有参考作用。延寿竖直摆潮汐因子、潮汐因子中误差、相对误差和相位滞后变化幅度都比较小，相对稳定。通过对2007年仪器工作状态较好的数据进行计算分析发现：EW、NS两个分项M2波的潮汐因子中误差分别可达0.0021和0.0024，属一类台站观测标准 [5] 。12个月的潮汐因子中误差中，除了一月、五月、六月以外其他月份都低于5% (潮汐因子低于5%为一类台站标准)，延寿竖直摆2007年达到一类台站的标准百分比为79.2%。

4. 竖直摆倾斜主要干扰因素分析

钻孔倾斜观测作为一种高精度的形变观测手段，河流、湖泊、仓库(货场、列车编组站)和铁路是几种主要的影响钻孔倾斜观测的地面载荷变化干扰源 [2] 。竖直摆倾斜观测是利用一个竖直悬挂的重力摆来检测地球表面的倾斜变化，钻孔倾斜观测也会不同程度地受到气象、仪器系统噪声等因素的影响。

4.1. 仪器零漂

延寿地震台钻孔倾斜2004年开始观测以来，由于受周围环境、地质环境、井孔介质特性等原因，仪器出现较大幅度的零漂现象。表1列出了2005~2015年各年度的漂移速率变化(通过每年年末的月均值和

图1. 竖直摆倾斜仪

表1. 2005~2015年竖直摆零漂和年幅度变化

年初月均值来求得)。通过计算可以得出2005~2015年期间，NS向年漂零值由2005年6.55降到2015年2.45，年幅度变化由2005年11.7降到2015年5.93，EW向年漂零值由2005年7.87降到2015年3.88，年幅度变化由2005年8.18降到2015年3.54，研究期间段内，两个分向整体上呈逐步下降趋于稳定的曲线形态。

4.2. 气压干扰

气压的变化体现在负荷效应上，气压使得地表负荷改变致使对井孔的垂直附加应力的改变，从而干扰观测结果。经过多年的观测发现，气压对延寿竖直摆的观测影响并不是十分显著，究其原因是延寿竖直摆倾斜安装于井下50~70米深，因此气压的变化对于观测的干扰非常小。在观测期间多次出现气压发生巨大变化的情况，但同时间段或者滞后时间段，竖直摆曲线并无出现相应异常现象，如2009年5月18~19日观测气压变幅可达17.6 hpa，但对应时段的倾斜固体潮未受其影响(如图2，图3)。当然气压对竖直摆的干扰程度也会因台站的岩石介质、地理位置不同会有不同的响应。

4.3. 电压干扰

电压的突然变化或不稳定也是影响观测结果的因素之一。由于台站位于偏远山区，经常会出现突然停电现象或雷电天气，都会造成台站电压的突变或不稳定，从而影响仪器观测质量。但电压造成的干扰一般是短时间或瞬间的，通过对观测数据的分析可以发现，电压造成干扰的主要表现形式为“阶跃”式(如图4，图5)，日常观测中对电压造成的干扰比较容易识别。

4.4. 水库库容变化造成荷载变化影响因素分析

荷载是指地表重物负荷的变化，如物资集散、建筑兴废和雨雪河流的变化等等。在仪器观测井西南方约100米有一个新城水库，水库的库容量达到2500万m³。本文依据邱泽华研究的三维集中荷载分析方法建立模型，从而求得水库库容变化造成的荷载变化这种干扰源对竖直摆的影响程度。这种模型的特点是只考量水库重量变化程度，而不考量水库的面积和分布情况等。我们把库容的变化直接定义成重量

图2. 气压变化量

图3. 气压变化当日竖直摆曲线

图4. 电压不稳定造成的干扰

图5. 雷电造成的干扰

的变化，在这种情形下，可以利用以下的位移解公式 [2] ：

$u=\frac{Fxz}{4\pi G{R}^3}-\frac{Fx}{4\pi \left(G+\lambda \right)R\left(R+z\right)}$ (1)

式中， $x$ 为观测井到水库中心的距离， $y$ = 0， $z$ 为探头深度， $u$ 表示在 $x$ 的水平方向上的位移， $F$ 为集中力， $\lambda$ 和 $G$ 是拉梅常数， $R=\sqrt{x^2+y^2+z^2}$ 。我们只需计算水平方向的倾斜 $w_x=\frac{\partial u}{\partial x}$ 即：

$w_x=-\frac{3Fx{z}^2}{4\pi G{\left(x^2+y^2+z^2\right)}^{5/2}}-\frac{F\left(2G+\lambda \right)x}{4\pi G\left(G+\lambda \right){\left(x^2+y^2+z^2\right)}^{3/2}}$ (2)

对于E = 5.0 × 10¹⁰ Pa，v = 0.25， $G=\frac{E}{2\left(1+v\right)}$ = 2.0 × 10¹⁰ Pa， $\lambda =\frac{vE}{\left(1+v\right)\left(1-2v\right)}$ = 2.0 × 10¹⁰ Pa。式中E是杨氏模量，v是岩石泊松比。探头深度z = 60 m，可绘制出 ${\omega }_X$ 随x变化的曲线图，如图6所示：

2011年观测井周边的新城水库的最大库容变化量为902万m³，最大的变化范围出现在6月9日~6月16日期间，库容变化幅度达到778万m³，我们取这段时间的数据来分析库容的变化对竖直摆造成的影响程度。观测井距离水库中心约100 m，在库容变化F = 778 × 10⁴ t时，利用公式(2)可以计算出倾斜的变化量 ${\omega }_X$ = 0.45 × 10⁻⁹。倾斜干扰限定指标日变化为1.45 × 10⁻⁸ (或3 × 10⁻³) [2] ，通过把计算结果和限定标准进行对比可知，2011年最大库容变化的期间对竖直摆的影响程度低于限定标准，即竖直摆倾斜未受到水库库容变化的干扰。

6月9日~6月16日期间，新城水库受降水影响库容发生巨大变化，变化量达到778万m³，但是竖直摆倾斜固体潮实际变化范围为0.0476~0.328之间(如表2)，低于理论限定标准。以上可以看出新城水库

图6. 倾斜- ${\omega }_X$ 随距离x的变化曲线

表2. 2011年6月9日~6月16日水库库容及倾斜变化范围

库容的增加或减少与倾斜趋势性异常无直接关系。

4.5. 降雨变化影响因素分析

延寿县属于温带大陆性季风气候，冬季这里受来自高纬内陆偏北风的影响，盛行极地大陆气团，气候寒冷干燥。夏季受极地海洋气团或变性热带海洋气团影响，盛行东风和东南风，暖热多雨，雨热同季。本地区年降水量可达500~1000 mm，约有2/3集中于夏季。由于冬季无降雨，所以选择每年的5~10月来进行计算降雨对观测的影响，取其整点值进行相关和对比分析，计算长度选取2009~2011年的每年5~10月，尽量避免随机因素和年周期的影响。

大气降水分析以小时值为基本单位，将中间无间断的连续降雨视为一次降雨过程；采用一次降雨过程量≥2 mm，能补给地下水的大气降水时间段为研究样本。假设在该短时间域内，大气降水因素对倾斜观测的影响是线性的，通过计算得出EW向相关系数0.00136，NS向相关系数−0.15603，可以看出大气降水因素对倾斜观测的影响并不大 [6] 。

大气降水分析以小时值为基本分析单位，将中间无间断的连续降雨视为一次降雨过程；采用一次降雨过程量≥2 mm，能补给地下水的大气降水时间段为研究样本。假设在该短时间域内，大气降水因素对倾斜观测的影响是线性的，取EW向分量观测值资料来计算与大气降水的线性相关系数。当相关系数lrl > 0.8时，视为高度相关；当0.5 ≤ lrl < 0.8时，视为中度相关；当0.3 ≤ lrl < 0.5时视为低度相关。从图7可以看出，在2010年5~10月期间，并没有相关系数lrl> 0.5(中高度相关)的时候，可以表示降雨与竖直摆间并无明显的相关性。

图7. 2010年5~10月竖直摆倾斜EW分量观测值与大气降水的相关系数

5. 结论与讨论

通过对延寿竖直摆的主要干扰因素进行分析，可以得出以下初步结论：

1) 电压变化、同震现象、天气因素、水库蓄水等干扰因素是客观存在的，也分别对观测仪器造成一定的影响，但影响极为有限，均不是主要干扰源。

2) 通过对降雨量和倾斜观测值的线性关系分析发现，短期内降雨对竖直摆倾斜观测的影响不大。

3) 通过三维集中荷载模型对新城水库库容变化对延寿竖直摆影响的定量分析发现：在水库库容量变化最大的期间(变化量达到778万m³)，倾斜的变化量 ${\omega }_X$ = 0.45×10⁻⁹，略低于限定标准1.45 × 10⁻⁸，认定库容变化对仪器是有影响的，但影响程度很小。

4) 影响仪器曲线形态和观测质量的主要因素为仪器零漂，但研究期间，零漂值及年幅度变化量整体上呈下降和趋于稳定的情形发展。

基金项目

中国地震局地震监测、预测、科研三结合课题：CEA-JC/3JH-160802。

文章引用

孙强,王强,李继业,孙伟,贾乃圣元,朱家明. 延寿地震台钻孔倾斜干扰因素分析
Analysis of Interference Factors of Borehole Inclinometer in Yanshou Seismic Station[J]. 测绘科学技术, 2018, 06(03): 223-230. https://doi.org/10.12677/GST.2018.63025

参考文献