文章快速检索     高级检索
  大地测量与地球动力学  2021, Vol. 41 Issue (8): 871-876  DOI: 10.14075/j.jgg.2021.08.018

引用本文  

苏维刚, 马茹莹, 马震, 等. 钻孔应变的干扰源及其干扰机制的多方法分析[J]. 大地测量与地球动力学, 2021, 41(8): 871-876.
SU Weigang, MA Ruying, MA Zhen, et al. The Multi-Method Analysis of Interference Source and Interference Mechanism of Borehole Strain[J]. Journal of Geodesy and Geodynamics, 2021, 41(8): 871-876.

项目来源

中国地震局“三结合”课题(3JH-202001092);中国地震局震情跟踪定向工作任务(2020010302)。

Foundation support

Combination Project with Monitoring, Prediction and Scientific Research of Earthquake Technology, CEA, No. 3JH-202001092; The Earthquake Tracking Task of CEA, No. 2020010302.

第一作者简介

苏维刚,工程师,主要从事地震前兆资料分析研究,E-mail:yangyuwan_su@163.com

About the first author

SU Weigang, engineer, majors in analysis of earthquake precursor data, E-mail: yangyuwan_su@163.com.

文章历史

收稿日期:2020-10-26
钻孔应变的干扰源及其干扰机制的多方法分析
苏维刚1,2     马茹莹2     马震2     孙玺皓2     
1. 中国科学院青海盐湖研究所,西宁市新宁路18号,810008;
2. 青海省地震局,西宁市兴海路1号,810001
摘要:德令哈台YRY-4钻孔应变观测数据长期受环境干扰影响。本文基于多种方法,通过分析观测资料与气象因素的相关性及环境水体蓄放水的载荷影响,追溯不同水体的地球化学特征,判断钻孔应变所受干扰来源及干扰机制。研究表明,在长时间尺度上,德令哈台钻孔应变存在稳定、显著的年周期性气温、气压干扰和不稳定、细微的地下水渗漏扰动;在短时间尺度上,德令哈台钻孔应变每年4~11月前后受人工湖和水库蓄放水影响,存在载荷效应。
关键词钻孔应变干扰源相关性载荷效应地球化学特征

钻孔应变观测作为一种连续监测地壳形变的重要手段,对小空间尺度上微小而重要的构造应变捕捉有较好的效果,对认识和探索地震孕育过程的动力学过程具有重要意义[1-4]。由于其应变灵敏度较高(分辨率达10-9~10-11),河流、湖泊、仓库和铁路等所有与地面载荷变化有关的干扰源都容易对其观测造成干扰[5-8],同时气象因素的变化、钻孔地质条件和格值的变化也会对钻孔应变观测造成不同程度的影响[5, 9-12]

德令哈台四分量钻孔应变仪架设于2006年,仪器自洽性较好,但受环境干扰影响严重。为分离干扰信息,本文使用多种方法研究台站及周边干扰源,并进行地球化学特征、载荷及气象因素等相关性分析,以期明确德令哈台钻孔应变观测的干扰源及干扰机制。

1 台站及观测系统概况

德令哈地震台位于柴达木盆地北缘,区域大地构造位置处于古亚洲构造域与古特提斯-喜马拉雅构造域的交界处,自新生代以来,盆地演化主要受印度板块与欧亚板块的碰撞作用影响。盆地内部发育一系列新生代高角度逆冲断裂,新构造活动迹象十分明显[13-14],沿活动断裂带发生过多次中强地震。德令哈地区属荒漠半荒漠大陆气候,年降雨量约100 mm。

德令哈台钻孔应变观测井位于德令哈地震台院内,钻孔于2006-09建成,孔径约132 mm,井深约43.5 m,井孔深度0.90 m以上为填土,以下为花岗岩。四分量钻孔应变仪器型号为YRY-4,探头埋深约40 m,具有NS、EW、NE、NW等4个测向,各分量方位角见表 1,数采采样率为1 s。

表 1 各分量方位角 Tab. 1 Azimuth angles of each component
2 研究方法 2.1 相关性分析

选取钻孔应变和气象三要素2015~2019年的日值数据,其中降水量测项由于WYY-1型仪器数采死机导致数据多次断迹,2017-11-15更换为RTP-Ⅱ型仪器,因此降水数据选取时间为2018~2019年。

通过SPASS软件对不同尺度的观测资料与温度和气压进行相关性分析,以气象要素为变量x,观测值为变量y,分别计算Pearson相关系数rp值,结果见表 2

$ r=\frac{\sum\limits_{i=0}^{n}\left(x_{i}-\bar{x}\right)\left(y_{i}-\bar{y}\right)}{\sqrt{\sum\limits_{i=0}^{n}\left(x_{i}-\bar{x}\right)^{2}+\sum\limits_{i=1}^{n}\left(y_{i}-\bar{y}\right)^{2}}} $
表 2 钻孔应变与气象三要素的相关关系 Tab. 2 The correlation between borehole strain and meteorological elements

式中,r值在-1~+1之间,n为样本量。

2.2 载荷分析

黑石山水库及西海公园人工湖与钻孔应变的观测距离均小于2 km(图 1),根据环境水体干扰源面积与最小“安静”距离的关系[6],水库及人工湖蓄放水可能会对应变观测造成干扰影响。现有钻孔应变的荷载模型主要有地表集中荷载模型和三维有限线段荷载模型[15-16],黑石山水库呈长条形,人工湖呈圆形,水库蓄水对钻孔应变观测的干扰问题可近似使用三维集中载荷模型进行分析:

$ \begin{gathered} e_{x}=-\frac{3 F x^{2} z}{4 G {\rm{ \mathsf{ π} }} R^{5}}+\frac{F z}{4 G {\rm{ \mathsf{ π} }} R^{3}}+ \\ \frac{F x^{2}}{4 {\rm{ \mathsf{ π} }} R^{2}(z+R)^{2}(G+\lambda)}+ \\ \frac{F x^{2}}{4 {\rm{ \mathsf{ π} }} R^{3}(z+R)(G+\lambda)}-\frac{F}{4 {\rm{ \mathsf{ π} }} R(z+R)(G+\lambda)} \end{gathered} $
图 1 德令哈采样位置 Fig. 1 Sampling location at Delhi

式中,F为库容重量,λG为拉梅常数,x为干扰源至钻孔的直线距离,z为钻孔深度。

2.3 地球化学分析

考虑汛期和非汛期水库水量的变化可能会造成差异性影响,分别于2020-04和2020-07采集德令哈水温井(距应变井约5 m)、西海公园人工湖及黑石山水库中的水样进行常量元素的测试分析,取样量为0.1 L。常量元素送至中科院青海盐湖研究所进行测试,测试结果见表 3(单位mg/L)。同时于2020-07采集3端元的水样进行镭同位素测试分析。

表 3 不同水体的常量元素测试结果 Tab. 3 The results of constant element in different water bodies

现场采集50 L的水体样品,利用锰纤维单柱富集法进行镭同位素富集[17],装填15 g锰纤维于富集柱,并在富集柱前面接入1 μm的过滤桶,富集时的水流速控制在450~550 ml/min左右,之后将锰纤维吹干至其含水率在1.5%以内[18],并通过四通道同步延时计数器RaDeCC进行223Ra和224Ra的比活度测试,测试结果见表 4(单位dpm/L)。

表 4 不同水体的Ra活度 Tab. 4 Ra activities of different water bodies
3 钻孔应变观测影响因素分析

德令哈钻孔应变四分量资料连续性较好,从长趋势变化来看,各个分量均存在不同幅度的趋势转折,转折节点与省内6级以上强震的发生存在一定的响应关系(图 2),表明德令哈钻孔应变仪可能位于构造敏感位置。但由于该测项受干扰影响严重,各分量均存在显著的张压性起伏变化,需分析其干扰源及干扰机制,以有效提取异常信息。

图 2 德令哈钻孔应变日均值时序曲线 Fig. 2 The daily average temporal curve of Delhi borehole strain
3.1 气象因素相关性分析

通过表 2可以看出,德令哈钻孔应变主要受年际尺度的气象因素(气温和气压)影响,并存在显著的负相关关系,气温和气压的变化具有稳定的年变周期,因此选用MAPSIS软件通过傅里叶滑动去掉钻孔应变资料中的年周期干扰[19],获得的钻孔应变观测曲线见图 3

图 3 德令哈钻孔应变去年周期后的日均值时序曲线 Fig. 3 The daily average temporal curve after canceling the annual variation period of Delhi borehole strain

图 3可以看出,去除年周期干扰影响之后,德令哈钻孔应变总体比较平稳,但在年尺度内仍然存在相对比较显著的张压性变化,且该变化不具有稳定的周期性。因此,需进一步考虑台站周边长期存在的环境干扰因素,如黑石山水库和西海公园人工湖等。

3.2 水库及人工湖蓄放水的加卸载作用

黑石山水库每年4月汛期开始蓄放水,蓄水量变化较为显著,11月初水库库容量最低且蓄水量较为稳定。人工湖具体蓄放水时间不确定,根据现场异常核实调查,每年4月中旬左右蓄水,7月左右水量达到最大值,9月初水量开始减少。选取2015年德令哈台钻孔应变数据变化幅度较大的3个时段,计算水库和人工湖蓄放水的载荷效应,得到理论最大干扰值,并利用德令哈台钻孔应变经傅里叶滑动去掉年周期干扰之后观测数据计算出该时段的最大主应变,详细参数与计算结果分别见表 56

表 5 载荷模型计算参数 Tab. 5 The computational parameters of the load model

表 6 载荷模型计算结果与实测计算结果 Tab. 6 The calculated results of load model and measured results

通过表 6可见,3个时段中,模型计算值均为实测计算值的3倍左右,说明德令哈钻孔应变在观测中确实受到周边环境水体蓄放水的影响,且在5~9月受人工湖蓄放水和水库蓄放水2个方面的载荷影响,9~12月主要受水库蓄放水的影响。

3.3 德令哈台及周边水体地球化学特征分析 3.3.1 水体常量元素分析

由测试结果可以看出,虽然德令哈井、黑石山水库和海西公园人工湖水样矿化度有所区别,但2个时段的3组水样在水-岩平衡图(图 4)中均位于“未成熟水”区域,并集中分布在Mg端元附近,表明这些水体与周围环境的水-岩反应程度均较弱,循环深度较浅[20]。由此推断,德令哈台井水的补给来源主要为大气降水和浅层的地下水,较少有深部来源的物质。

图 4 德令哈井及周边水样水-岩平衡图 Fig. 4 Water-rock equilibrium map of Delhi well and surrounding water samples

每年4月初水库和人工湖近乎处于干涸状态,7月处于蓄水量高值。通过图 5可以看出,水库和人工湖水体中Ca2+、HCO3-及SiO2的浓度在4月和7月存在显著的差异性变化,但应变井水体离子浓度均保持稳定,无明显变化。综合分析后初步判定,水库和人工湖水量的大幅变化未对井水体造成直接影响,水库和人工湖对应变井产生渗漏的可能性较小。

图 5 不同水体的常量元素浓度特征对比 Fig. 5 Comparison of the concentration characteristics of constant elements in different water bodies
3.3.2 水体镭同位素分析

镭是一种天然放射性的同位素,自然水体中的镭同位素(223Ra、224Ra、226Ra和228Ra)来源于不溶于水却在沉积物和岩石中广泛分布的钍的衰变。其中,223Ra和224Ra的半衰期分别为11.4 d和3.66 d,分别为研究10~50 d和1~10 d时间尺度的水体混合与循环的理想示踪剂。地下水中的镭同位素活度普遍较高,一般比地表水(河水、湖水、海水等)的活度高出1~2个数量级[21-23]

通过表 4可以看出,井水中223Ra活度最高,为0.022 dpm/L,水库次之,人工湖最低;而224Ra活度井水为1.7 dpm/L,水库和人工湖的224Ra活度均低于检测线。根据文献[24],河水平均Ra活度较低,大气降水Ra活度低于检测线,而水库水源主要为河流和大气降水,因此水库223Ra活度低于河水平均Ra活度,而224Ra半衰期更短,自身Ra的衰变影响更大,导致224Ra活度低于检测线;由于人工湖水源主要为水库引水,水体没有流动性,无法及时更新,Ra活度衰变较为彻底,因此人工湖223Ra和224Ra均低于检测线。

采集德令哈水温观测井井水样品时发现,水样中有明显腐烂气味,一般认为这种水样流动性差,Ra活度衰变会较为彻底,但其Ra活度与地下水平均活度相当,且224Ra活度还高于地下水平均224Ra活度,因此认为有其他高Ra活度的地下水体混入,从而导致井水Ra活度增大。

4 结语

通过气象因素相关性、载荷效应和不同端元水体的地球化学特征分析对德令哈台钻孔应变观测资料的干扰来源和干扰机制进行深入研究。结果发现,德令哈台钻孔应变观测受多种干扰因素影响,在长时间尺度上存在稳定、显著的年周期性气温、气压干扰和不稳定、细微的地下水渗漏扰动;在短时间尺度内受人工湖和水库蓄放水的载荷干扰,时间为每年的4~11月。

随着人类活动的日益加剧,观测环境受干扰影响的程度日趋严重,且干扰因素互相叠加出现,互相影响,给观测资料的分析和异常核实工作造成很大的困难。在观测资料的异常核实和异常提取过程中可应用多种方法多方面探寻可能存在的干扰源和干扰机制,避免单一分析方法给研究结果带来的不确定性和偏差。

致谢: 感谢中国科学院青海盐湖研究所孔凡翠博士在地球化学分析方面提供的帮助!

参考文献
[1]
苏恺之. 地形变连续观测技术基本问题的思考[J]. 大地测量与地球动力学, 2005, 25(1): 131-135 (Su Kaizhi. Discussion on Some Problems of Continual Observation Technique of Crust Deformation[J]. Journal of Geodesy and Geodynamics, 2005, 25(1): 131-135) (0)
[2]
Sturkell E, Einarsson P, Sigmundsson F, et al. Volcano Geodesy and Magma Dynamics in Iceland[J]. Journal of Volcanology and Geothermal Research, 2006, 150(1-3): 14-34 DOI:10.1016/j.jvolgeores.2005.07.010 (0)
[3]
倪友忠, 叶青, 朱国豪, 等. 利用四元件钻孔应变观测资料研究佘山地区主应变及主方向变化特征[J]. 大地测量与地球动力学, 2019, 39(10): 1 091-1 095 (Ni Youzhong, Ye Qing, Zhu Guohao, et al. On Variation Characteristics of Main Strain and Main Direction in Sheshan Area Using Four Element Borehole Strain Observations[J]. Journal of Geodesy and Geodynamics, 2019, 39(10): 1 091-1 095) (0)
[4]
邱泽华, 唐磊, 赵树贤, 等. 用应变地震观测求解震源矩张量的基本原理[J]. 地球物理学报, 2020, 63(2): 551-561 (Qiu Zehua, Tang Lei, Zhao Shuxian, et al. Fundamental Principle to Determine Seismic Source Moment Tensor Using Strain Seismographs[J]. Chinese Journal of Geophysics, 2020, 63(2): 551-561) (0)
[5]
李杰, 刘敏, 邹钟毅, 等. 数字化钻孔体应变干扰机理及异常分析[J]. 地震研究, 2003, 26(3): 230-238 (Li Jie, Liu Min, Zou Zhongyi, et al. Analysis of Disturbance Mechanism and Abnormity of Digital Observation Data of Borehole Body Strain Meters[J]. Journal of Seismological Research, 2003, 26(3): 230-238 DOI:10.3969/j.issn.1000-0666.2003.03.004) (0)
[6]
邱泽华. 钻孔应变观测点到地面载荷干扰源最小"安静"距离的理论分析[J]. 岩石力学与工程学报, 2004, 23(23): 4 063-4 067 (Qiu Zehua. Theoretic Studies on the Shortest 'Quiet' Distance between Borehole Strain Observatory Site and Surficial Mechanical Disturbance Sources[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(23): 4 063-4 067) (0)
[7]
马栋, 陈建国, 郑爽, 等. 荷载对钻孔应变观测影响的实例分析[J]. 地震研究, 2014, 37(1): 79-85 (Ma Dong, Chen Jianguo, Zheng Shuang, et al. Example Analysis of Influence of Load on the Borehole Strain Observation[J]. Journal of Seismological Research, 2014, 37(1): 79-85) (0)
[8]
贾媛, 甘卫军, 李杰, 等. 抽水对泰安地震台钻孔分量应变的影响及其机制[J]. 地震地质, 2019, 41(6): 1 429-1 443 (Jia Yuan, Gan Weijun, Li Jie, et al. Influence of Pumping on Four-Component Borehole Strain at Tai'an Seismic Station and Its Interference Mechanism[J]. Seismology and Geology, 2019, 41(6): 1 429-1 443) (0)
[9]
闫伟, 陈国琴, 杨跃文, 等. 以丽江地震台为例研究潮汐观测仪器格值变化检测及格值归化方法[J]. 地震研究, 2015, 38(4): 613-616 (Yan Wei, Chen Guoqin, Yang Yuewen, et al. Method of Scale Change Detection and Scale Calibration Illustrated by the Case of Tide Observation Instrument at Lijiang Seismic Station[J]. Journal of Seismological Research, 2015, 38(4): 613-616 DOI:10.3969/j.issn.1000-0666.2015.04.013) (0)
[10]
赵玉红, 冯丽丽, 李霞, 等. 青海地区大地电场异常变化特征分析[J]. 国际地震动态, 2019(8): 57 (Zhao Yuhong, Feng Lili, Li Xia, et al. Analysis on the Variation Characteristics of Geoelectric Field Anomalies in Qinghai Region[J]. Recent Developments in World Seismology, 2019(8): 57 DOI:10.3969/j.issn.0253-4975.2019.08.047) (0)
[11]
陈建生, 董海洲, 余波, 等. 利用线热源法研究堤防集中渗漏通道[J]. 地球物理学进展, 2003, 18(3): 400-403 (Chen Jiansheng, Dong Haizhou, Yu Bo, et al. Using Line-Heat-Resource Method to Study Concentration Leakage Passage of Dike[J]. Progress in Geophysics, 2003, 18(3): 400-403 DOI:10.3969/j.issn.1004-2903.2003.03.008) (0)
[12]
马震, 马茹莹, 苏维刚. 德令哈地震台分量钻孔应变干扰分析[J]. 地震地磁观测与研究, 2017, 38(4): 147-152 (Ma Zhen, Ma Ruying, Su Weigang. Analysis of Interference Characteristics of Borehole Strain at Delingha Seismic Station[J]. Seismological and Geomagnetic Observation and Research, 2017, 38(4): 147-152 DOI:10.3969/j.issn.1003-3246.2017.04.025) (0)
[13]
Meng Q R, Fang X. Cenozoic Tectonic Development of the Qaidam Basin in the Northeastern Tibetan Plateau[J]. Special Paper of the Geological Society of America, 2008, 444: 1-24 (0)
[14]
袁道阳. 青藏高原东北缘晚新生代以来的构造变形特征与时空演化[D]. 北京: 中国地震局地质研究所, 2003 (Yuan Daoyang. Tectonic Deformation Features and Space-Time Evolution in Northeastern Margin of the Qinghai-Tibetan Plateau since the Late Cenozoic Time[D]. Beijing: Institute of Geology, CEA, 2003) (0)
[15]
胡卫建, 张俊山, 谢智, 等. 荷载对钻孔应变测值影响的实验及力学解析[J]. 地震, 2002, 22(3): 95-104 (Hu Weijian, Zhang Junshan, Xie Zhi, et al. Experimental Result and Theoretic Analysis Concerning the Borehole Strain Responding to the Load on Earth's Surface[J]. Earthquake, 2002, 22(3): 95-104) (0)
[16]
骆鸣津, 池顺良, 马鸿钧. 地层绝对应力测量与钻孔应变测量[M]. 北京: 地震出版社, 2008 (Luo Mingjin, Chi Shunliang, Ma Hongjun. Measurement of Absolute Stress and Borehole Strain of Formation[M]. Beijing: Seismological Press, 2008) (0)
[17]
Moore W S. Sampling 228Ra in the Deep Ocean[J]. Deep Sea Research and Oceanographic Abstracts, 1976, 23(7): 647-651 DOI:10.1016/0011-7471(76)90007-3 (0)
[18]
Sun Y, Torgersen T. The Effects of Water Content and Mn-Fiber Surface Conditions on 224Ra Measurement by 220Rn Emanation[J]. Marine Chemistry, 1998, 62(3-4): 299-306 DOI:10.1016/S0304-4203(98)00019-X (0)
[19]
蒋骏, 李胜乐, 张雁滨, 等. 地震前兆信息处理与软件系统: EIS 2000[M]. 北京: 地震出版社, 2000 (Jiang Jun, Li Shengle, Zhang Yanbin, et al. Earthquake Precursor Information Disposal and Software System: EIS 2000[M]. Beijing: Seismological Press, 2000) (0)
[20]
刘磊, 马茹莹, 张朋涛, 等. 青海省地震地球化学背景场特征研究[J]. 国际地震动态, 2019(8): 50 (Liu Lei, Ma Ruying, Zhang Pengtao, et al. Study on the Characteristics of Earthquake Hydrogeochemical Background Field in Qinghai Province[J]. Recent Developments in World Seismology, 2019(8): 50 DOI:10.3969/j.issn.0253-4975.2019.08.041) (0)
[21]
Bollinger M S, Moore W S. Radium Fluxes from a Salt Marsh[J]. Nature, 1984, 309(5 967): 444-446 (0)
[22]
Levy D M, Moore W S. 224Ra in Continental Shelf Waters[J]. Earth and Planetary Science Letters, 1985, 73(2-4): 226-230 DOI:10.1016/0012-821X(85)90071-8 (0)
[23]
Burnett W C, Bokuniewicz H, Huettel M, et al. Groundwater and Pore Water Inputs to the Coastal Zone[J]. Biogeochemistry, 2003, 66(1-2): 3-33 DOI:10.1023/B:BIOG.0000006066.21240.53 (0)
[24]
苏维刚, 沙占江, 孔凡翠, 等. 青海湖泉湾镭同位素的空间分布与沿岸地下水排放评价[J]. 湖泊科学, 2019, 31(6): 1 749-1 757 (Su Weigang, Sha Zhanjiang, Kong Fancui, et al. Evaluation of the Groundwater Discharge and the Spatial Distribution of Radium Isotope in Quanwan Bay of Qinghai Lake[J]. Journal of Lake Sciences, 2019, 31(6): 1 749-1 757) (0)
The Multi-Method Analysis of Interference Source and Interference Mechanism of Borehole Strain
SU Weigang1,2     MA Ruying2     MA Zhen2     SUN Xihao2     
1. Qinghai Institute of Salt Lakes, CAS, 18 Xinning Road, Xining 810008, China;
2. Qinghai Earthquake Agency, 1 Xinghai Road, Xining 810001, China
Abstract: The observation data of Delhi YRY-4 have been affected by environmental interference for a long time. Based on various methods, this study judges and analyzes the interference source and mechanism of borehole strain by analyzing the relativity between observation data and meteorological factors, the load influence of environmental water storage and release, and the geochemical characteristics of different water bodies. The results show that there are stable and significant annual periodic air temperature and pressure disturbances, unstable and subtle groundwater leakage disturbances on a long time scale. In a short time scale, there exists load interference of water storage and release in artificial lakes and reservoirs, which last from April to November every year.
Key words: borehole strain; interference source; relativity; effect of load; geochemical characteristics