引言

体积式应变仪(简称体应变仪)因工作原理简单、易获取连续可靠的观测资料, 成为国内外主要的钻孔应变观测仪器(中国地震局监测预报司, 2008)。青岛地震台体应变仪采用中国地震局地壳应力研究所研制的TJ-2型体应变仪(苏恺之等, 2003)。2015年以来, 青岛地震台周边相继出现两处较显著的工程施工干扰。2018年需在地震台院内安装CZB钻孔倾斜仪, 进行钻孔施工。2015—2016年降雨偏少, 导致青岛地震台钻孔水位较往年偏低。上述干扰因素对青岛体应变观测资料产生显著影响(高福旺等, 2004；苏恺之等, 2005；李光科等, 2015), 甚至可能掩盖前兆异常信息, 造成异常信息提取困难。本文研究青岛体应变仪长趋势和周期特性, 并分时段分析降雨、钻孔水位、周围工程施工等因素对体应变资料的影响, 得到有益结论, 可为研究地下水位等因素对钻孔体应变观测的影响及今后异常信息的提取提供参考。

1.   台站及观测系统概况

青岛地震台位于山东省青岛市主城区, 青岛市在大地构造上位于华北狼林地块与扬子京畿地块的对接带——胶南超高压变质带上(王红霞等, 2005)。区内主导断裂包括即墨断裂、沧口断裂、劈石口断裂、王哥庄断裂等NNE向构造(见图 1)。青岛地震台紧邻劈石口断裂和王哥庄断裂(张继红等, 2010), 最后一次断错地表活动发生于中更新世末至晚更新世初, 晚更新世以来没有发生活动(栾光忠等, 2009；潘元生等, 2004)。台站基底为破碎花岗岩, 覆盖层厚3m。《山东省历史地震目录》中关于青岛地区史料记载到的历史地震较少, 历史上未发生过M_L5.0以上地震。据山东地震台网记录, 青岛地区现代地震活动呈现活动微弱、分布零散、低频度、低强度等特点, 属于少震、弱震地区(王红霞等, 2005；潘元生等, 2004)。地震台周边区域2007年以来主要以M_L＜4.0地震活动为主, 地震活动较弱, 除东北方向的乳山震群较集中外, 其他地区分布较零散(见图 1)。

图 1  青岛地震台周围地区断裂分布及地震活动

Figure 1.  Fault distribution and seismic activity around Qingdao seismic station

注:右上方插图表示左侧大图在华北地区的位置, 其中蓝色矩形框为左侧大图所指区域, 黑色三角形表示青岛地震台位置；填色圆圈为2007年以来的地震分布, 圆圈大小表示震级(参考图中下侧图例), 填充颜色表示发震时刻(参考图中右侧图例)；黑色实线为断裂分布(晁洪太等, 1997)；F1为即墨断裂；F2为沧口断裂；F3为劈石口断裂；F4为王哥庄断裂

下载: 全尺寸图片 幻灯片

青岛地震台TJ-2型体应变仪于2006年4月13日安装运行, 可进行辅助钻孔水位观测。仪器所在钻孔孔径130mm, 钻孔深度66m, 仪器安装深度55m。钻孔岩石性质为花岗岩(见表 1), 岩体十分坚硬, 但破碎仅在58—59m、54—55m深度处相对完整。岩石节理面与裂隙的倾角近乎垂直, 受地下水位干扰严重(苏恺之, 2006), 如图 2所示。

表 1  青岛体应变钻孔地层分布

Table 1.  Distribution of strata in Qingdao volumetric strain borehole

层号	层底深度/m	地层描述
1	12.0	强中风化花岗岩, 岩石松软破碎, 矿物成分以石英、长石为主, 岩石结构清楚
2	26.2	硅化石英正长岩, 坚硬破碎, 矿物成分以石英为主
3	34.8	微风化花岗岩, 节理裂隙发育, 较破碎, 矿物成分以石英、长石为主
4	41.7	硅化石英正长岩, 坚硬, 破碎, 矿物成分以石英为主, 接触带有糜棱岩化
5	51.9	微风化花岗岩, 节理裂隙发育, 较破碎, 矿物成分以石英、长石为主
6	56.8	微风化花岗岩, 节理少, 较完整, 矿物成分以石英、长石为主
7	66.2	微风化花岗岩, 硅化石英正长岩, 坚硬, 较破碎, 矿物成分以石英为主

下载: 导出CSV 

| 显示表格

图 2  青岛体应变(去除二次多项式趋势项)及钻孔水位观测

Figure 2.  Qingdao volumetric strain (detrend the term of quadratic polynomial) and borehole water level observation

下载: 全尺寸图片 幻灯片

2.   降雨对钻孔水位的影响

降雨是造成钻孔水位变化的常见因素, 由于青岛地震台无降雨观测手段, 所以从青岛市水文局收集降雨观测数据。该观测站与青岛地震台直线距离约250m, 可反映青岛地震台体应变仪所在区域的降水情况。降雨观测属于非连续型数据, 直接观测无法反映降雨对体应变连续的累积影响机制。对降雨资料做以下预处理:

① 将原始雨段观测值转化为每日降雨量序列${P_n}(n = 1, 2, \quad 3, \ldots, m)$, m为序列长度；

② 对每日降雨量序列做累加处理, 即

③ 将$P_n^\prime $去除线性趋势, 获得累积降雨线性拟合残差序列$P_{n}^{r}$。可见, $P_{n}^{r}$考虑了累积降雨的连续影响, 并反映了各时刻降雨量与正常雨量的偏离程度。$P_{n}^{r}$序列与体应变及钻孔水位观测的对比(见图 3)表明, 青岛地震台钻孔水位与降雨情况(通过累积降雨线性拟合残差$P_{n}^{r}$表示)具有明显的同步特征, 其中2015—2016年降雨偏少, 导致青岛地震台钻孔水位较往年低。可见青岛地震台地下水主要由降雨补给, 降雨下渗改变了钻孔水位, 并影响体应变资料的变化。

图 3  降雨对青岛体应变及水位的影响

Figure 3.  Influence of rainfall on Qingdao volumetric strain and water level

下载: 全尺寸图片 幻灯片

3.   工程施工对钻孔应变的影响

岩石孔隙中富含流体时, 会产生孔隙流体压力效应, 岩石孔隙内流体压力称为孔隙压力, 孔隙压力可抵消岩石构造围压作用(曾佐勋等, 2013)。岩石中的有效应力(佘成学等, 2010)可描述为:

其中, σ为岩石中的构造应力；β为孔隙压力作用系数；P_w为孔隙水压力。台站周围区域地震活动(见图 1)表明青岛地震台近期地震活动和构造应力较平稳, 而钻孔水位和钻孔环境的变化对β、P_w产生影响, 因此可认为青岛体应变所反映的应力变化主要由式(2)中的βP_w项产生。

2015年后青岛地震台周边开始两处施工:一处为台站向北150m处的青岛市立医院住院大楼基坑开挖, 另一处为台站向北850m处的青岛市地铁2号线(距青岛地震台最近的地铁线)燕儿岛站施工。青岛市立医院住院大楼工程总建筑面积86000m², 设计地上15层、地下3层, 于2015—2016年进行基坑开挖。青岛市地铁2号线呈东西走向, 根据青岛新闻网记者2015年5月30日的报道(吴帅等, 2015), 地铁2号线浮燕区间采用盾构机施工, 该区间地质条件复杂, 存在富水砂层、破碎带及上软下硬地层等不良地质, 于2016年7月29日完成土建施工。这两处大规模的工程施工涉及基坑开挖, 改变台站周边的岩石应力应变状态, 影响岩石孔隙压力作用系数, 对应变仪观测造成影响。

2018年1月17—18日进行青岛地震台距体应变观测井8m的CZB钻孔倾斜仪钻孔, 先用XY-3-500型岩芯钻机装备金刚石钻头钻进, 但由于台站岩石坚硬, 钻进速度缓慢, 后改用400型潜孔锤钻机负责钻孔钻进和护井管安装。该钻机利用压缩空气作为循环介质, 并作为驱动孔低冲击器的能源, 使之上下往复运动, 并连续不断地对下部钻具施加一定频率的冲击载荷, 进行冲击回转钻进(王达等, 2014)。钻进过程中, 潜孔锤利用对钻头的往复冲击作用破碎岩石, 并将岩屑排出至地表。潜孔锤钻机较岩芯钻机钻进速率快, 但钻孔孔径大, 且对钻孔岩体具有更大的破坏性。CZB钻孔和体应变钻孔均穿过了含水层, 钻孔施工可能对区域岩石孔隙环境造成影响, 引起体应变钻孔水位的永久阶变(见图 4)。当岩石孔隙裂缝增多, 岩石渗透率增大时, 直接扩大孔隙压力作用系数, 同样的水位上升变化引起更大的孔隙水压力, 导致体应变观测出现加速上升, 可解释2018年3月以来水位小幅上升引起体应变显著压性上升的现象。最新的观测数据(见图 5)表明, 2018年1月17日(钻孔后)至2019年1月(青岛体应变仪工作状态开始改变)受近距离钻井的影响, 钻孔应变井周边岩石孔隙环境及水位对体应变的影响作用显著增强(见图 5), 验证了上述推测。

图 4  钻孔施工前后青岛体应变及钻孔水位、降雨观测

Figure 4.  Qingdao volumetric strain and borehole water level, rainfall observation before and after drilling

下载: 全尺寸图片 幻灯片

图 5  青岛体应变(去除一次趋势项)及钻孔水位观测

Figure 5.  Qingdao volumetric strain (detrend the one degree term) and borehole water level observation

下载: 全尺寸图片 幻灯片

4.   工程施工对钻孔水位与体应变相关性的影响

青岛体应变存在稳定的长趋势变化特征(见图 2), 分离长趋势变化对资料趋势异常识别至关重要。通常以多项式(3)确定仪器的长趋势项:

已有研究表明, 钻孔水位是体应变常见的趋势性干扰因素(马震等, 2011；卢双苓等, 2012), 可改变观测的短时趋势, 将体应变观测资料的长趋势项正确分离, 可提高体应变与钻孔水位的相关性。因此, 可利用钻孔水位与体应变的相关性检验确定多项式的阶数。考虑长期特征及稳定特性, 尽量通过低阶多项式进行描述, 当低阶多项式不能描述时增加阶数。通过取n=1, 2, 3, 4, 5逐次对观测数据进行多项式拟合, 并分析分离长趋势项后的体应变与钻孔水位的相关性, 通过最优相关系数确定多项式阶数。考虑2015年1月以来青岛体应变辅助钻孔水位出现显著变化, 分2009—2018年和2015—2018年两个时段对体应变长趋势项进行研究。由于长趋势项不考虑固体潮因素, 为减少数据量, 将整点值观测数据转为日值数据进行相关性分析。本文利用线性相关系数表示相关性, 并对其进行T检验, 计算结果(见表 2)表明, 青岛体应变2009年以来长趋势项呈二次多项式分布, 而2015年以来的长趋势项呈一次多项式分布, 可利用式(3)对青岛体应变分离长趋势项, 并研究钻孔水位的影响。

表 2  青岛体应变去除趋势后与钻孔水位的相关系数

Table 2.  Correlation coefficients between the detrended Qingdao volumetric strain and borehole water level

多项式阶数	2009年1月1日至2018年4月25日(样本数3394)		2015年1月1日至2018年4月25日(样本数1210)
	相关系数	T检验显著性/%	相关系数	T检验显著性/%
1	0.12	100.0	0.35	100.0
2	0.43	100.0	0.07	98.5
3	0.27	100.0	0.33	100.0
4	0.20	100.0	0.29	100.0
5	0.20	100.0	0.29	100.0

下载: 导出CSV 

| 显示表格

由图 2、表 2中的相关性分析可知, 青岛体应变受钻孔水位的影响。2015年后, 青岛体应变与钻孔水位的相关性出现显著变化, 开始呈现具有滞后特征的相关性(见表 3)。

表 3  青岛体应变与钻孔水位滞后相关性

Table 3.  Lag correlation between Qingdao volumetric strain and borehole water level

序号	时间段	相关系数	样本数	滞后特征	最优相关系数
1	2009-01-01—2015-01-01	-0.09	2191	不显著	—
2	2015-01-02—2018-01-17	0.57	1112	显著(105d)	0.76
3	2018-01-18—2018-04-25	0.89	97	不显著	—

下载: 导出CSV 

| 显示表格

为研究滞后性, 对于钻孔水位观测数据, 进行一定的滞后改正后分析相关性, 即在一定时间范围内不断调整滞后改正时间, 通过收敛至最优相关系数确定最优滞后时间。如果给定时间范围内未收敛至最优相关系数, 则认为滞后特征不显著。考虑季节性周期因素, 滞后时间范围设定为0—180d。考虑2015年1月以来青岛体应变辅助钻孔水位出现显著变化、2015—2017年青岛降雨量减少及台站周边存在工程施工干扰、2018年1月17日在距体应变钻孔8m处存在显著的钻孔施工干扰等情况, 分3个时间段进行分析计算, 时间段划分及计算结果如表 3所示。

2009年1月1日至2015年1月1日青岛体应变与钻孔水位相关性相对较弱(相关系数为-0.09), 在给定的时间范围内调整滞后改正时间未显著收敛至最优相关系数(见图 6(a)), 未表现出显著的滞后特性, 相关性较低, 使得两者的散点分布非常零散(见图 6(b))。2015年1月2日至2018年1月17日青岛体应变与钻孔水位呈现显著的线性相关性(相关系数为0.57), 但具有滞后特征(见图 7(a)), 在给定的时间范围内调整滞后改正105d时显著收敛至最优相关系数, 此时, 两者的散点分布呈线性分布(见图 7(b)), 可用线性函数y=0.803x+0.130拟合, 相关系数提高至0.76。2018年1月18日至2018年4月25日两者依然显著相关(相关系数为0.89), 且滞后特征消失(见图 8(a)), 两者几乎同步变化(见图 8(b)), 可用线性函数y=1.092x+0.086拟合。上述结果表明青岛体应变与钻孔水位的相关性在2015年前后出现变化, 2015年以来青岛地震台钻孔水位对体应变的影响作用趋于增强。

图 6  2009-01-01—2015-01-01青岛体应变与钻孔水位相关性分析

Figure 6.  Correlation analysis of Qingdao volumetric strain and borehole water level from January 1, 2009 to January 1, 2015

下载: 全尺寸图片 幻灯片

图 7  2015-01-02—2018-01-17青岛体应变与钻孔水位相关性分析

Figure 7.  Correlation analysis of Qingdao volumetric strain and borehole water level from January 2, 2015 to January 17, 2018

下载: 全尺寸图片 幻灯片

图 8  2018-01-18—2018-04-25青岛体应变与钻孔水位相关性分析

Figure 8.  Correlation analysis of Qingdao volumetric strain and borehole water level from January 18, 2018 to April 25, 2018

下载: 全尺寸图片 幻灯片

利用傅里叶变换分析表 3中3个时段青岛体应变及钻孔水位周期特性(见图 9), 周期特性基于整点值数据进行计算。2009年1月1日至2015年1月1日, 青岛体应变除具有固体潮引起的0.5d、1d周期特性外, 还表现出显著的年、半年周期特性。该时段内钻孔水位仅有年周期特性较明显, 两者400d周期以内的振幅谱曲线相关性较低。2015年1月2日至2018年1月17日, 青岛体应变半年周期特性有所衰弱, 半年至一年周期特性增强, 两者400d周期以内的振幅谱曲线相关性显著提高。2018年1月18日以来, 由于观测时间较短, 未能展示完整的周期特性。青岛体应变及钻孔水位的周期特性分析结果也表明两者的相关性在2015年前后出现显著变化。

图 9  青岛体应变及钻孔水位的周期特性

Figure 9.  Periodic characteristics of Qingdao volumetric strain and borehole water level

下载: 全尺寸图片 幻灯片

5.   结论

通过分析降雨、工程施工、钻孔水位等因素对青岛体应变的影响, 得出以下结论:

(1) 青岛地震台钻孔水位反映的区域地下水主要由降雨补给, 降雨下渗改变地下水状态, 进而影响青岛体应变的资料变化。

(2) 青岛体应变与钻孔水位的相关性在2015年前后出现变化, 2015年1月2日至2018年1月17日表现出滞后相关性。两者的周期特性在2015年前后也出现变化, 可能与2015—2016年受降雨偏少引起的低钻孔水位及周边工程施工干扰有关。

(3) 青岛体应变与钻孔水位在CZB钻孔倾斜仪钻孔施工前后几乎同步变化, 很可能是由于青岛地震台CZB钻孔倾斜仪钻孔施工对区域岩石孔隙环境造成了影响。当岩石孔隙裂缝增多, 岩石渗透率增大, 直接扩大孔隙压力作用系数, 同样的水位上升变化将引起更大的孔隙水压力, 导致体应变观测出现加速上升, 也可解释2018年3月以来水位小幅上升引起体应变显著压性上升的现象。