Key technology and application of DRY-1B capacitive component borehole strain gauge
-
摘要:
文章简述了DRY-1B型电容分量式钻孔应变仪（简称“钻孔应变仪”）的理论基础，攻关了微位移感知、降噪、控温、性能测试、标定等关键技术，并通过了性能测试，标定结果表明：该应变仪达到了高分辨率（≥5×10−11ε）、宽频带（10～100 Hz可选）、大动态范围（≥1×10−3ε）、24位AD记录、低功耗（<3 W）等技术指标，其性能优于同期美国PBO和日本同类钻孔应变仪，是一款国际领先的地壳运动长期观测仪器，基本能够满足地壳长期应变缓慢积累的蠕变运动和短期应变快速变化的地震火山活动等观测需求。2008年以来，通过20余个地应力台站的应用，该钻孔应变仪记录到大量的地壳形变、断裂活动、同震应变波、应变阶跃、矿压活动等应变信息，并以北长山地应力台站应变监测数据自洽性检验和土耳其地震映震能力分析为例发现：北长山地应力台站电容传感器1#+3#和2#+4#元件应变曲线总体平稳，相关系数R2达到0.95；1#－3#和2#－4#元件的差应变年变化速率为10−8量级，反映出长岛地区以剪切应力为主，且处于地震活动相对高的应力环境；利用该应变仪观测到2023年2月6日土耳其M 7.8级和M 7.5级两次地震明显的同震应变响应，尤其是获取了M 7.8级主震面波周期为50~60 s，呈现出面波异常，理论上可分辨出100 km范围M 0.74级地震产生的应变波，达到了应用示范效果。该钻孔应变仪在地球动力学研究、内动力地质灾害监测等领域具有较好的推广价值和应用前景。
Abstract:This paper briefly describes the theoretical basis of the DRY-1B capacitive component drilling strain gauge (drilling strain gauge). It discusses critical technologies such as micro-displacement sensing, noise reduction, temperature control, performance testing, and calibration. The strain gauge achieved high resolution (≥5×10−11ε), wideband (10–100 Hz optional), extensive dynamic range (≥1×10−3ε), 24-bit AD recording, low power consumption (< 3W), and other technical indicators. Its performance is better than the United States PBO and Japan borehole strain gauge of the same period, and it is an international leading long-term observation instrument for crustal movement, which can basically meet the observation requirements of creep movement with slow accumulation of long-term strain and seismic and volcanic activity with a rapid change of short-term strain. Since 2008, through the application of more than 20 geostress stations, the borehole strain gauge has recorded a large amount of strain information, such as crustal deformation, fault activity, co-seismic strain wave, strain step, and ore pressure activity. Based on the results of the self-consistency test of strain monitoring data of the geostress station in the Beichangshan Mountain and the analysis of seismic mapping capacity of the Turkey earthquake, it is found that the strain curves of the 1#+3# and 2#+4# capacitance sensors at the Beichangshan Mountain stress station are generally stable, and the correlation coefficient R2 is 0.95. The annual variation rate of the differential strain of 1#–3# and 2#–4# elements is 10−8 magnitude, which reflects that the shear stress is dominant in the Long Island area and the stress environment of seismic activity is relatively high. The strain gauge was used to observe the apparent co-seismic strain response of Turkey’s M 7.8 and M 7.5 earthquakes on February 6, 2023. In particular, it obtained the M 7.8 main seismic surface wave period of 50–60 s, presenting an out-facing wave anomaly. Theoretically, the strain wave generated by the M 0.74 earthquakes within 100 km can be distinguished, and the application demonstration effect has been achieved. The borehole strain gauge has good popularization value and application prospects in geodynamics research and internal dynamic geological disaster monitoring.
-
0. 引言
近年来，随着煤矿开采深度不断加深，动力灾害的发生日趋频繁。其中，冲击地压灾害的破坏力巨大，造成的后果严重影响了煤矿安全开采。从工作面的尺度出发，煤层在开采过程中常伴随上覆坚硬岩层的突然破裂，引发岩体强烈的冲击破坏，表征为一定范围内的破坏现象（纪洪广等，2003）。近年来，专家学者着眼于工作面尺度之外更大的范围去研究煤矿相关动力灾害的发生机制及预警，将研究范围扩展至区域应力场的范畴（李铁等，2003，刘金海等，2014；张月征等，2016a）。甘肃华亭煤矿于2014年发生的一次冲击地压灾害，其释放的能量达到了108 J量级，造成 1500 m巷道底鼓1.5 m，与工作面尺度内的顶板结构失稳相比，该灾害表现出破坏力强、波及范围广等特点，其根本原因在于周边区域内煤岩体为冲击地压灾害的发生提供了必要的力能条件（张月征等，2016b）。因此，在工作面开采期间，应从更广的区域范围入手，寻求不同尺度的动力异常响应特征，以期探索深部煤矿开采过程中动力灾害的发生及预警机制。
天然地震和矿震都是地层内部的动力学现象，可将天然地震超前预报的观测手段应用到矿震或冲击地压灾害的监测中，其中地壳形变观测是地震地面监测体系中重要的前兆观测项目（邱泽华和石耀霖，2004），监测台站观测到地应力拉伸脉冲异常与地壳地震活动显著相关（邱泽华，2017）。采用S变换和高通滤波从钻孔应变数据中提取短周期信号，能够发现汶川地震前钻孔应变观测出现异常（刘琦和张晶，2011）。借助钻孔应变观测仪器实现对地壳变形的精确记录，可及时捕捉地壳内部如固体潮汐量级的微小变化（苏恺之，2003；邱泽华，2010；张聪聪，2015）。钻孔应变仪因其频带宽、精度高，适用于监测煤矿深部开采动力灾害异常响应特征。为了揭示煤层开采过程中区域应力场的演变规律，需要考虑开采扰动与区域应力场协同作用对矿震活动性的相互影响，往往需要捕捉长短周期信号（云龙等，2021）。相比监测煤矿动力灾害的其他设备，钻孔应变仪因其监测频带宽的特性（池顺良，2013），所采集的数据含有高频短周期的噪声和低频长周期的应变固体潮频段以及能够反映矿震活动规律的特定频段，为探究工作面采动环境中矿震孕育与诱发机制，需滤去噪声和应变固体潮频段。
针对时域内难以有效区分煤矿动力灾害异常响应信息的问题，文章以鲍店煤矿七采区东翼不规则区域（73上11工作面）开采期间观测的钻孔应变数据为基础，采用变分模态分解（VMD）方法，获得不同频率的固有模态函数（IMF）分量。通过构建带通滤波通道，滤去高斯噪声信号和固体潮分量，重构反映矿震活动的IMF分量，利用希尔伯特（Hilbert）变换将时域内看似杂乱无章却能反映矿压活动的应变信息转换成频域信息，求解该分量的瞬时能量和以天为单位的能量累计变化，分析信号异常处矿震的发生与区域应力场之间的相互作用效应，实现对矿震或冲击地压灾害前兆信息的识别和解析，以期对相似地质条件下煤矿动力灾害的预测预报提供一定指导意义。
1. 工程概况及资料选取
鲍店煤矿73上11工作面位于井田东北部边界，东部与东滩煤矿四采区（43上09采空区）、十四采区（14316采空区和14302采空区）相接，北部与兴隆庄煤矿四采区相接，相邻矿井的3个采区回采工作已完成，西部为鲍店煤矿73上02采空区和73上04采空区，南部为十采区、七采区（73上09采空区）南翼采区，相邻采区已基本采完（图1）。因此，工作面之间相互关系呈现孤岛工作面开采。其中，七采区东翼南侧已回采，待回采的73上11工作面处于不规则三角形区域（图1）。73上11工作面主采3#煤层，煤层采深超过600 m，且距煤层上方100 m处有厚硬砂岩层，其破断后释放大量弹性应变能，该厚硬岩层的失稳对工作面安全开采影响显著。根据地质勘探结果，已探明鲍店煤矿七采区厚硬砂岩最小厚度为41.75 m，最大厚度超过373.55 m。由采区东北部边界向东西翼分界线逐渐增加，73上11工作面上覆厚硬砂岩层厚度大于300 m，部分区域超过350 m，工作面周边采空区覆岩运动以及应力场调整势必增加73上11工作面应力分布的复杂性。综上所述，73上11工作面开采期间主要受孤岛工作面导致的应力集中和上覆厚硬岩层稳定性的影响，故73上11工作面在掘进和回采过程中动力异常响应较为突出，易引发矿震或冲击地压相关动力灾害。
图 1 鲍店煤矿73上11工作面邻近关系示意图（海拔 –546.8 m）Figure 1. Neighboring relationship diagram of the working face position (The altitude is −546.8 m)为了揭示矿震发生前的异常响应特征，选取能量为105J的矿震事件作为研究对象，截取5月份和10月份每次矿震发生前20天和震后10天作为一个完整的时间序列，基于VMD识别和提取钻孔应变数据。
2. 基于VMD识别和提取钻孔应变数据
工作面开采期间，钻孔应变仪因其具有高分辨率、高精度等特性，采集的数据会受到所处环境中干扰因素的影响，如高斯噪声的即时干扰和应变固体潮的周期性干扰等（文勇等，2012）。消除噪声常见的滤波方法主要有小波变换、傅里叶变换或短时变换（穆峰等，2015），上述方法在处理平稳信号时所表现的效果是满足滤波需求的，但是在处理非平稳信号时误差较大，因非平稳信号的特性是其信号自身频率随时间变化而变化，傅里叶变换在处理相关问题时只能给出总体效果，并不能够反应出某一时刻信号的局部信息。在此基础上，Huang et al.（1998）提出经验模态分解（EMD）来处理非平稳、非线性信号，其特点是将原始信号按照频率从高到低的顺序分解为若干个固有模态函数分量，分解后的各模态函数具备其自身固有的属性；然后，再根据希尔伯特变换将信号从时域角度转换至频域角度，这在面对时域内非线性、非平稳信号分析效果差的环境下，提供了有力的技术支撑。但随着EMD分解方法在信号处理方面的广泛应用，其缺点也暴露出来，如其端点效应和模态混叠严重，不具备严格的数学模型。基于此问题，在EMD基础上又发展出很多改进算法，如集合经验模态分解（EEMD；Huang and Wu，2008；陈可等，2010）、补充集合经验模态分解（CEEMD；Yeh et al.，2010）和改进集合经验模态分解（MEEMD；郑近德等，2013）等。EEMD是通过向原始信号中加入高斯白噪声来抑制经验模态分解过程中存在的模态混叠现象，CEEMD是通过向目标信号添加成对的符号相反的白噪声，来确保信号分解具有真实的物理意义，MEEMD结合CEEMD与排列熵（PE）算法在抑制模态混叠方面取得理想的结果（邢昀和荣剑，2018）。上述相关改进算法虽然在一定程度上解决了模态混叠和端点效应问题（康佳星等，2016），但本质上仍是时域范畴内的递归求解。随后，Dragomiretskiy and Zosso（2014）提出的一种自适应信号处理方法（VMD），通过迭代不断更新各模态函数及中心频率，得到若干具有一定带宽的模态函数，获取变分模态的最优解。VMD相较于EMD是基于频域的非递归求解算法，弥补了EMD的缺点。
井下工作环境、地质环境复杂多变（谭元隆等，2021），钻孔应变仪采集的数据记录了不同来源、不同特性的混合时频信号，仅仅从应变数据波形特征的时变角度难以精确分析局部细节信息。但诸如噪声、应变固体潮等某个单频信号蕴含着固有的优势频率，采用VMD的方法，能够将具有不同来源、不同频率的波动信号分离出来，进而针对不同的研究问题，选择敏感频率段的数据作为研究对象。
文章采用四分量钻孔应变仪分析采动环境中应变的相对变化规律，四分量钻孔应变仪通过记录元件电容值实现对应变的观测。为界定观测数据质量是否满足研究问题的要求，需要依靠质量评价体系评定数据的可靠性，而钻孔应变数据质量评价体系则是自洽方程（式）。四分量钻孔应变仪采集的数据经过应变换算后得到面应变Sa（邱泽华等，2009），依靠自洽方程既可以评判数据自身的可靠性，也可以用来识别四分量钻孔应变仪的工作状态。若数据段内整体变化趋势一致，则满足自洽方程：
式中：${S_{\text{1}}}$—元件1的观测数值；${S_{\text{2}}}$—元件2的观测数值；${S_{\text{3}}}$—元件3的观测数值；${S_{\text{4}}}$—元件4的观测数值。
动力灾害的孕育及其诱发与区域应力场的力学行为密切相关，这里的动力灾害主要涉及冲击地压、矿震以及工作面来压期间发生的顶板异常动力显现事件。针对能级跃迁较大的煤矿动力灾害，钻孔应变曲线走势表现为显著的“阶跃”（张月征等，2016b），即数据曲线在动力灾害发生时刻产生台阶式上升或下降。此外，还存在一些能级跃迁较小的动压事件。5月4日—5月31日四分量钻孔应变仪分量1和分量3根据自洽方程拟合后的观测数据（图2a），9月22日—10月13日分量2和分量4根据自洽方程拟合后的观测数据（图2b），均满足自洽方程，图中仅绘制拟合后的观测曲线。观察两个不同时间序列上钻孔应变数据的幅值变化，可知某一时刻发生动力异常事件的典型状态失稳现象在应变数据波动特征上并不能够有效直观地显现，应力传递路径受钻孔应变仪安装位置、震源处煤岩介质力学属性等影响，动力异常事件发生前后微小的响应信号大都淹没在固体潮背景下，故有必要根据时频变换识别或提取区域应力场中细微变化与动力灾害的前兆响应特征。
图 2 时频解析前钻孔应变仪采集原始数据曲线Figure 2. Raw data curves by the borehole strainmeter before time-frequency analysis钻孔应变仪采集的是多频段混合后的数据（欧阳祖熙等，2009），经VMD分解后得到若干个固有模态函数（IMF）。定义有意义的瞬时频率其必要条件是函数关于零均值局部对称且过零点的数目与极点的数目一致。据此，Huang et al.（1998）定义了IMF应满足以下2个条件：①在整个数据区间内，极值点的数目与过零点的数目相等或至多相差1个；②在任意一点处，由局部极大值点定义的包络以及由局部极小值点定义的包络的均值为零。此外，为了消除噪声和固体潮等因素对要研究数据段的影响，在VMD分解的基础上，采用希尔伯特（Hilbert）变换依次分析IMF分量，将时域信息转换至频域信息范畴。按照上述理论对5月1日—5月31日的数据段进行分解，获得各频段IMF曲线（图3）。
图 3 5月1日—5月31日钻孔应变数据VMD分解曲线图Figure 3. VMD decomposition curves (IMF1—IMF13 ) of borehole strain data from May 1st to May 31st如表1所示，IMF1—IMF5周期皆小于1 h，推测可能是设备自身或周围环境中产生的高斯噪声（Barman et al., 2016），故剔除IMF1—IMF5分量，以排除噪声对原始数据的干扰。 IMF9和IMF10的频率分别为2.24×10−5 Hz与1.16×10−5 Hz，转换至周期后分别为12.4 h和23.9 h，时间维度上与固体潮严格匹配，确定为固体潮的半日波和日波分量。此外，与IMF9和IMF10邻近的IMF 11和IMF12为固体潮谐波分量。
表 1 钻孔应变数据VMD分量频率周期表Table 1. Frequency periodic table of VMD component of borehole strain dataIMF 1 2 3 4 5 6 频率/Hz 3.98×10−3 2.62×10−3 1.19×10−3 6.64×10−4 3.17×10−4 2.41×10−4 周期/h 0.07 0.10 0.23 0.41 0.89 1.15 IMF 7 8 9 10 11 12 频率/Hz 1.53×10−4 6.37×10−5 2.24×10−5 1.16×10−5 1.51×10−6 7.49×10−7 周期/h 1.8 4.3 12.4 23.9 185 370 为探究噪声频段对所研频段数据造成的干扰，综上所述，将混杂的各频段数据解析后重构，对比分析降噪前与降噪后的频段曲线（图4）。降噪信号段与原始信号背景趋势走向一致，对钻孔应变数据的局部异常突跳点捕捉效果良好，降噪信号能够从局部微小信号上反映钻孔应变的波动特征。
利用VMD对含噪信号按照频率从高到低依次分解，得到若干个固有模态分量。通过构建带通滤波器，滤去高频噪声项和低频固体潮背景趋势项后，将其余IMF分量重构。
通过上述方法，分析用于监测矿山动力灾害和区域应力场变化的钻孔应变数据。依据变分模态分解原理，既能排除原始信号中的高频噪声，又能剔除低频固体潮背景数据。因此，该方法充分保留了特定信号段非线性、非平稳的固有属性，经频谱分析得到IMF1—IMF12频率与周期的优势分布，对比每个IMF的固有属性和信号时频差异性，从而实现钻孔应变数据的降噪处理。
3. 基于希尔伯特（Hilbert）变换的区域应力异常提取方法
冲击地压和矿震都属于煤矿常见的动力灾害，两者既有相似之处，也有其内在的差异。通常由于积聚的弹性应变能瞬间释放致使周围煤岩体的破坏，从而引起冲击地压和矿震。前者多发生在巷道或工作面以及开采扰动近源区，并伴随有煤岩体抛出的动力现象，而后者是开采扰动与区域应力场共同作用诱发的地震，多发生在工作面开采扰动远源区，其在表征动力破坏程度上弱于冲击地压，且不发生煤岩体的抛出现象。此外，从二者的发生频次上分析，冲击地压灾害数要远小于矿震，但当满足一定的临界条件时，矿震会诱发冲击地压灾害的发生。因此，区域应力场监测是实现矿震预测、预警的基础和前提，且有利于揭示冲击地压的致灾机理（姜耀东等，2014）。
为了得到区域应力场的实时变化，将经过VMD分解后得到的精细信号段进一步分析和处理。常用的方法有傅里叶变换和希尔伯特变换。其中，傅里叶变换的优势表现为对平稳或线性信号的处理，若使用该方法处理非平稳、非线性的信号，则需进行特定的变形后方可使用，如：小波变换和短时傅里叶变换。但经过小波变换和短时傅里叶变换解析的信号会出现负频率，负频率在分析相关问题时不具有任何物理意义，而希尔伯特变换在解析过程中虚数部分较实数部分滞后$\text{π} /2$，借用此特性可消除负频率的影响。因此，希尔伯特变换在解析非线性、非平稳信号时，有其自身特有的优势。
根据矿震事件发生的时间，以分钟为单位提取钻孔应变数据，经希尔伯特变换后得到矿震分量的瞬时幅值、瞬时能量和每日累计能量（图5，图6）。由5月1日—5月31日和9月22日—10月13日的钻孔应变瞬时幅值时程变化曲线可知（图5a，图6a），在时间序列内统计矿震发生频次共计6次，分别为5月21日、5月26日和5月29日及9月24日、10月4日和10月11日，从瞬时幅值和累计能量两个角度均反应了相同的变化规律。其中，日累计能量的变化趋势说明区域应力场内能量在不断的积聚、释放和耗散（图5b，图6b）。
图 5 5月1日—5月31日重构分量动力异常响应幅值和能量时程变化Figure 5. Dynamic anomaly response amplitude and energy time–history change of reconstructed components from May 4th to May 31st
图 6 9月22日—10月13日重构分量动力异常响应幅值和能量时程变化Figure 6. Dynamic anomaly response amplitude and energy time–history change of reconstructed components from September 22nd to October 13rd3.1 从幅值曲线波动形态上分析
矿震与地震都是地下岩层运动后产生形变进而导致弹性应变能的积累并释放（蒋骏等，2000；彭华等，2008，2011；王凯英，2018），两者在发生机制上具有相似的特征。Mescherikov（1968）、Yoichiro（1974）通过对大量地震资料的归纳总结发现，在地震的孕育和诱发过程中存在α、β和γ三个阶段，运用该理论在后续的地震跟踪预测中取得了显著效果。为解释地震形变三阶段的物理实质，冯德益（1983）认为α阶段以岩层的微小破裂为主，持续扩展乃至相互贯通，导致最终发生非震滑动或岩层的蠕变；β阶段则是因断层活化而引起的突变异常，并最终进入γ阶段，表征为大量积蓄弹性应变能的释放。如图7所示，地震形变三阶段的曲线走势与岩石破坏的应力应变曲线相似，运用该理论分析煤矿矿震的活动规律，以期探究矿震发生的前兆特征。以5月21日矿震事件为例，应变幅值“有升有降”走势契合地震发生的前期α、中期β和后期γ三个形变阶段。α阶段：应变稳定变化，该阶段虽然也会造成煤岩破裂，但区域应力场总体稳定，能量积聚速率与能量释放速率趋于一致，不足以引发动力灾害。β阶段：应变急剧变化，表现形式为“震荡”+“突跳”，这是相关区域岩层破裂产生非稳态应力调整现象的表征，能量积聚与能量释放处于临界状态，如5月21日矿震前期先发生“震荡”后发生“突跳”，而5月26日矿震β阶段应力异常表征为“震荡”、“震荡”+“突跳”2种应变异常表现形式受矿震震源处和“力链”范围内煤岩介质属性、应力传递路径、开采扰动影响范围，以及震源、工作面和钻孔应变安装设备三者相对位置的影响，区分2种异常形式有利于了解区域应力场的赋存状态，并为未来的煤矿开采工作提供技术指导。γ阶段：大量积聚弹性应变能瞬间释放，造成煤岩体破裂并引发矿震。
图 7 地震形变三阶段物理模型（宋治平等，2004）Figure 7. Three-stage physical model of seismic deformation (Song et al., 2004)3.2 从形变异常开始时间分析
钻孔应变数据经VMD识别并提取能够反映矿震的IMF分量。该IMF分量经希尔伯特变换后所得时间幅值变化关系，亦是矿震震源、工作面开采位置和钻孔应变仪安装位置三部分组成的系统所形成的一定空间范围内区域应力场变化的表征。不同于天然地震发生前2年应变出现异常的规律（宋先月等，2003)，解析出矿震形变异常早于其发生2～7天（表2，表3），这是天然地震与煤矿矿震的孕震环境不同所造成的。工作面开采扰动促使区域应力场不断的调整，对比观察开采扰动作用下的区域应力场幅值随时间的变化，形变的三个阶段中稳定变化阶段（α）的应变格值为0.5×10−10ε，形变快速变化阶段（β）的应变格值为稳定变化阶段的3倍以上，呈现为 “突跳” 的变化特征。动力灾害发生前，区域应力场内应力（应变）显著升高，积蓄大量的弹性应变能，当相关区域的可释放弹性应变能超过岩石表面应变能时，即可沿破裂面瞬间释放。
表 2 5月份形变异常前兆模式特征表Table 2. Characteristics of deformation anomaly precursory mode in May发生时间 形变异常
开始时间应力调整
持续时间表现形式 5月21日 5月11日 10 d “震荡”+“突跳” 5月26日 5月24日 2 d “震荡” 5月29日 5月27日 2 d “突跳”+“震荡” 表 3 10月份形变异常前兆模式特征表Table 3. Characteristics of deformation anomaly precursory mode in October发生时间 形变异常
开始时间应力调整
持续时间表现形式 9月24日 9月22日 2 d “突跳” 10月4日 9月27日 7 d “突跳” 10月11日 10月6日 5 d “突跳”+“震荡” 基于鲍店煤矿典型动力灾害−矿震事件，通过分析钻孔应变数据对区域应力场异常的响应，认为应变幅值由形变稳定变化阶段向快速变化阶段的幅值阶跃为形变异常开始时间。随后，区域应力场在开采扰动作用下不断调整直至进入震前失稳阶段，能量释放呈现出动力灾害形式。因此，对应变数据进行VMD和Hilbert变换后，可以识别并提取能够反映煤矿动力异常响应的IMF分量，并在地震形变三阶段的理论基础上，统计分析其幅值和累积能量的时程变化规律，进而尝试建立煤矿动力灾害的预测预报模式。
4. 存在的不足与展望
（1）基于鲍店煤矿所在区域地应力监测数据，采用Hilbert变换方法进行了解析，获得了煤矿动力灾害发生前明显的前兆异常响应特征，但不能代表不同构造环境下的应变前兆特征，未来还需进一步丰富不同构造环境下的应变监测数据采集与分析工作，以期归纳总结出煤矿动力灾害的预测预报模式。
（2）在对分解出的若干个IMF分量进行异常信号提取和识别工作中，由于煤矿井下环境复杂多变，设备采集到的原始数据受到严重干扰，以上研究中所采用的方法仅是针对具体问题的定向研究，在一定程度上删除了干扰信号的影响，但是使用的算法有其自身的局限性，存在部分数据段无法消除干扰的问题，后期还需进一步优化算法，以期提高数据分析的准确性。
5. 结论
采用VMD和Hilbert变换的方法，分析了鲍店煤矿开采期区域应力场变化的钻孔应变监测数据，结论如下：
（1）将VMD方法应用到钻孔应变数据处理中，能够将混频信号分解为若干个具有特定频段的固有模态分量，为后续动力异常响应特征的识别和提取提供无模态混叠的数据源。
（2）通过构建带通滤波器，消除高斯噪声和应变固体潮干扰项，重构信号获得能够反映区域应力场对矿震、冲击地压等动力灾害响应的IMF分量，并采用Hilbert变换方法处理特定频段IMF分量。
（3）依据地震形变三阶段理论，定义形变异常开始时间为区域应力场进入异常状态时刻，并经过形变快速变化阶段的“震荡”和“突跳”后，释放出积聚的弹性应变能，满足岩石强度失效和变形破坏准则时，即能引发动力灾害。
（4）基于钻孔应变仪对区域应力场的精确观测，将形变异常开始时间作为动力灾害事故发生前的响应时刻，为未来煤矿开采工程的安全性和稳定性评价提供了一种切实可行的前兆信息判别模式。
-
图 1 钻孔应变仪探头结构与电容传感器原理图
σ1—最大水平主应力，MPa；θ—测量轴线与地理北的夹角，（°）；φ—最大水平主应力与地理北的夹角，（°）；d1、d2—极板间距，μm
Figure 1. Diagram showing the working principle of capacitance displacement sensor
σ1–maximum horizontal principal stress, MPa; θ–angle between measurement axis and geographic north, (°); φ–angle between maximum horizontal principal stress and geographic north, (°)
图 2 电容位移传感器结构示意图
σ1—最大水平主应力，MPa；σ2—最小水平主应力，MPa；u-钻孔孔径位移量（即极板的位移量），μm；a、b、c—极板编号；d1—极板a和极板b之间的距离，μm；d2—极板b和极板c之间的距离，μm
Figure 2. Schematic showing the structure of capacitance displacement sensor
σ1–maximum horizontal principal stress, MPa; σ2–minimum horizontal principal stress, MPa; u–displacement of the pole plate, μm; a, b and c–pole plate numbers; d1–distance between pole plate a and pole plate b, μm; d2–distance between pole plate b and pole plate c, μm
图 3 双层驱动屏蔽环的电场整形作用
a—无屏蔽环−电极边缘场凌乱；b—有屏蔽环−电极边缘场稳定
Figure 3. The electric field shaping effect of the double-layer drive shielding ring
(a) Messy edge field of the no shielding ring-electrode; (b) Stable edge field of the shielded ring-electrode
图 5 分量式钻孔应变传感器线性标定曲线
Figure 5. Linear calibration curve of component-type borehole strain sensor
图 6 电容传感器微调与标定装置
a—电机驱动楔缝定位器（磁致伸缩标定器）；b—菱形螺杆定位器（压电陶瓷标定器）
Figure 6. Micro position-adjustment and calibration devices of capacitor sensor
(a) Motor driven wedge slider positioner (magnetostrictive calibrator); (b) Rhomboid screw positioner (piezoelectric ceramic calibrator
图 8 钻孔应变监测系统网络拓扑图
a—井下部分；b—井上部分
Figure 8. Network topology diagram of the borehole strain gauge monitoring system
(a) Under the shaft; (b) On the ground
图 9 长岛北长山地应力监测台站DRY-1B钻孔应变仪各分向元件安装方位图
Figure 9. Element orientation diagram of DRY-1B borehole strain gauge at Beichangshan in-situ stress station, changdao, Shandong
图 10 北长山台站四个分量应变趋势图（向上为压）
Figure 10. Curves of 4-component strain observed at the Beichangshan station (upward pressure)
图 11 2016年4月北长山台站1#元件（N15°）月应变曲线
Figure 11. Monthly strain curve of element 1# (N15°) at the Beichangshan station in April 2016
图 12 北长山台站应变趋势图（向上为压）
Figure 12. Strain observed curves of the Changdao station ( upward pressure)
图 13 2023.2.6土耳其两次地震同震应变曲线
Figure 13. Coseismic strain curves of the two earthquakes occured in Turkey on February 6, 2023
图 14 北长山台站记录的土耳其M 7.8级地震应变波Sn波到时（2月6日9时45分42秒）
Figure 14. Seismic strain Sn wave of the Turkey M 7.8 earthquake recorded by the Beichangshan station (9:45:42 on February 6 )
图 15 北长山台站记录的土耳其M 7.5级地震应变波Sn波到时（2月6日18时52分30秒）
Figure 15. Seismic strain Sn wave of the Turkey M 7.5 earthquake recorded by the Beichangshan station (18:52:30 on February 6)
表 1 系统技术指标
Table 1. Technical specifications of the borehole strain monitoring system
技术指标 技术参数 技术指标 技术参数 供电电压 12～48 V DC 数据传输模式/波特率 RS485/9600BPS 井下功耗 总功耗<3 W 通道/角度 4分量/45°分布 AD位数 内置数据采集，24位 观测分辨率 ≥5×10−11ε 采样速率 所有分量10～100 Hz 观测动态范围 ≥1×10−3ε 注：电子罗盘、温度、孔隙压问答式查询返回数据 
下载: 导出CSV
-
[1] BENIOFF H, PRESS F, SMITH S, 1961. Excitation of the free oscillations of the earth by earthquakes[J]. Journal of Geophysical Research, 66(2): 605-619. doi: 10.1029/JZ066i002p00605 [2] CHEN N, WANG C H, CHEN P Z, et al. , 2021. Reunderstanding of ground stress field characteristics of right bank of Baihetan using borehole caving data [J]. Journal of Geomechanics, 27(3): 430-440. (in Chinese with English abstract) [3] CHEN Y Q, ZHANG S L, LI H W, 2021. Analysis on characteristics of regional strain field based on four-component borehole strain[J]. Journal of Geodesy and Geodynamics, 41(10): 1095-1100. (in Chinese with English abstract) [4] CHI S L, 1977. Electric capacity stress meter of automatic continuous recording[M]//Selected papers of the national conference on stress measurement (part 2). 369-373. (in Chinese) [5] CHI S L, 1993. The results of experimental observations of a shallow-hole-installed YRY-2 borehole strainmeter at 8 stations in North China[J]. Acta Seismologica Sinica, 15(2): 224-231. (in Chinese with English abstract) [6] CHI S L, 2007. Deep-hole broad-band strain-seismograph and high-frequency seimology: the hope to successful earthquake prediction[J]. Progress in Geophysics, 22(4): 1164-1170. (in Chinese with English abstract) [7] CHI S L, ZHAGN J, CHI Y, 2014. Failure of self-consistent strain data before Wenchuan, Ludian and Kangding earthquakes and its relation with earthquake nucleation[J]. Recent Developments in World Seismology(12): 3-13, 47. (in Chinese with English abstract) [8] EVERTSON D W, 1977. Borehole strainmeters for seismology[R]. Austin: Applied Research Lab, University of Texas. [9] FAN J Y, LIU G C, 2021. Comparative analysis of coseismic response characteristics from three borehole strain gauges recorded for the Maduo MS7.4 earthquake[J]. Earhquake Research in China, 37(3): 564-573. (in Chinese with English abstract) [10] FAN Z N, WAN Y G, 2020. Study on co-seismic signal of the M7.8 earthquake in Sumatra, Indonesia, recorded by strainmeter[J]. Journal of Geodesy and Geodynamics, 40(8): 849-853, 876. (in Chinese with English abstract) [11] FENG C J, Li B, Li H, et al. 2022. Estimation of geostress field and discussion on structural stability in Nangabawa area [J]. Journal of Geomechanics, 28(6): 919-937. (in Chinese with English abstract) [12] FURUYA I, FUKUDOME A, 1986. Characteristics of borehole volume strainmeter and its application to seismology[J]. Journal of Physics of the Earth, 34(3): 257-296. doi: 10.4294/jpe1952.34.257 [13] FURUYA I, YAMASATO H, SEINO M, 1991. Rayleigh waves recorded by the volumetric strainmeters at the time of a rather small volcanic eruption[J]. Papers in Meteorology and Geophysics, 42(3): 93-103. doi: 10.2467/mripapers.42.93 [14] GLADWIN M T, 1984. High-precision multicomponent borehole deformation monitoring[J]. Review of Scientific Instruments, 55(12): 2011-2016. doi: 10.1063/1.1137704 [15] HE D Q, LI Y F, SHE R C, et al. , 2014. Near field anomalies of Yiliang M5.7, 5.6 earthquakes on Sep. 7, 2102[J]. Journal of Seismological Research, 37(1): 1-8. (in Chinese with English abstract) [16] Li B, ZHANG W, WEN R, 2022. Geostress measurement results and engineering geological significance analysis of hydraulic fracturing method for extra-long highway tunnel in southern Shaanxi Province [J]. Journal of Geomechanics, 28(2): 191-202. (in Chinese with English abstract) [17] LI F Z, REN T X, CHI S L, et al. , 2021. The dynamic Coulomb stress changes caused by remote earthquakes based on the borehole strainmeter data[J]. Chinese Journal of Geophysics Chinese Journal of Geophysics, 64(6): 1949-1974. (in Chinese with English abstract) [18] LUO M J, GU M L, LI A Y, et al. , 1989. Borehole strain-stress calibration in the same place by use of tidal generation force[J]. Crustal Deformation and Earthquake, 9(4): 51-62. (in Chinese with English abstract) [19] PENG H, WU Z H, MA X M, 2006. Unmanned in-situ stress monitoring stations along the Qinghai-Tibet railway[J]. Journal of Geomechanics, 12(1): 96-104. (in Chinese with English abstract) [20] PENG H, MA X M, JIANG J J, 2008. Analysis of the volume strain data from the Shandan in-situ stress monitoring station[J]. Journal of Geomechanics, 14(2): 97-108. (in Chinese with English abstract) [21] PENG H, MA X M, JIANG J J. et al. , 2011a. Strain characteristics of WFSD stress stations and its co-seismic effects analysis: a case study upon Mw9.0 earthquake in Japan[J]. Journal of Geomechanics, 17(1): 1-14. (in Chinese with English abstract) [22] PENG H, MA X M, JIANG J J, 2011b. In-situ stress measurement by differential strain analysis method in WFSD-1[J]. Journal of Geomechanics, 17(3): 249-261. (in Chinese with English abstract) [23] QIU Z H, SHI Y L, OUYANG Z X, 2005a. Absolute calibration of 4-component borehole strainmeters[J]. Earthquake, 25(3): 27-34. (in Chinese with English abstract) [24] QIU Z H, SHI Y L, OUYANG Z X, 2005b. Relative in-situ calibration of 4-component borehole strain observation [J]. Journal of Geodesy and Geodynamics, 25(1): 118-122. (in Chinese with English abstract) [25] QIU Z H, 2010. A review of component borehole observation of stress-strain in China[J]. Journal of Geodesy and Geodynamics, 30(5): 42-47. (in Chinese with English abstract) [26] QIU Z H, 2017. Theory and application of borehole strain observation[M]. Beijing: Earthquake Publishing House. (in Chinese) [27] QIU Z H, TANG L, ZHAO S X, et al. , 2020. Fundamental principle to determine seismic source moment tensor using strain seismographs[J]. Chinese Journal of Geophysics, 63(2): 551-561. (in Chinese with English abstract) [28] QUAN J J, LAI J S, CHEN S H, et al. , 2021. Analysis of observation quality and earthquake reflecting ability of borehole body strain meter at Xiaotao seismic station[J]. South China Journal of Seismology, 41(1): 26-33. (in Chinese with English abstract) [29] SACKS I S, SUYEHIRO S, EVERTSON D W, 1971. Sacks-evertson strainmeter, its installation in Japan and some preliminary results concerning strain steps[J]. Proceedings of the Japan Academy, 47(9): 707-712. doi: 10.2183/pjab1945.47.707 [30] Science and Technology Monitoring Department of the National Seismological Bureau, 1995. Seismic deformation observation technology[M]. Beijing: Earthquake Press: 229-250. (in Chinese) [31] SU K Z, 2003. Relative measurement method of crustal stress[M]//SU K Z, LI H L, ZHANG J, et al. New progress in borehole strain observation. Beijing: Earthquake Publishing House: 252-260. (in Chinese) [32] SU K Z, 2018. The four-component borehole strainmeter calculation formula and strain anomalies of two earthquakes[J]. Seismological and Geomagnetic Observation and Research, 39(5): 101-111. (in Chinese with English abstract) [33] SUN Y M, ZHA N, REN X, et al. , 2021. A review of application research on the borehole body strain gauge in earthquake precursory observation[J]. Journal of Disaster Prevention and Reduction, 37(4): 69-74. (in Chinese with English abstract) [34] WU L H, CHEN Z, LI T, et al. , 2010. Application of capacitive tilt sensor in crustal deformation measurement[J]. Transducer and Microsystem Technologies, 29(10): 129-131. (in Chinese with English abstract) [35] WU L H, XIONG Y Z, LI T, et al. , 2018. In-situ linear calibration technique of the RZB-type borehole strain meter: including a laboratory experiment[J]. Journal of Geodesy and Geodynamics, 38(4): 437-440. (in Chinese with English abstract) [36] WU L H, LI H, CHEN Z, et al. , 2022. Research on the stability of the RZB-type borehole strain sensor with function of the range expansion and in-situ calibration[J]. Journal of Geodesy and Geodynamics, 42(3): 318-320. (in Chinese with English abstract) [37] ZHENG Y D, ZHANG J J, ZHANG B, 2022. Two pillars theory of structural geology in the new century: Maximum effective moment criterion and deformation configuration decomposition [J]. Journal of Geomechanics, 28(3): 319-337. (in Chinese with English abstract) [38] ZHU S Y, CAO J W, FENG C J, et al. 2023. Study on fault instability risk induced by recharge water in Gaoyang low uplift deep hot reservoir, Hebei Province [J]. Journal of Geomechanics, 29(2): 220-235. (in Chinese with English abstract) [39] 陈念, 王成虎, 陈平志, 等, 2021. 利用钻孔崩落数据再认识白鹤滩右岸地应力场特征[J]. 地质力学学报, 27(3): 430-440. doi: 10.12090/j.issn.1006-6616.2021.27.03.039 [40] 陈永前, 张淑亮, 李宏伟, 2021. 基于四分量钻孔应变资料的区域应变场特征分析[J]. 大地测量与地球动力学, 41(10): 1095-1100. [41] 池顺良, 1993. 一种浅孔安装的YRY-2型钻孔应变仪在中国华北地区8个台站试验观测的结果[J]. 地震学报, 15(2): 224-231. [42] 池顺良, 2007. 深井宽频钻孔应变地震仪与高频地震学: 地震预测观测技术的发展方向, 实现地震预报的希望[J]. 地球物理学进展, 22(4): 1164-1170. doi: 10.3969/j.issn.1004-2903.2007.04.023 [43] 池顺良, 张晶, 池毅, 2014. 汶川、鲁甸、康定地震前应变数据由自洽到失洽的转变与地震成核[J]. 国际地震动态(12): 3-13, 47. [44] 樊俊屹, 刘高川, 2021. 3种钻孔应变仪记录的玛多Ms7.4地震同震响应特征对比分析[J]. 中国地震, 37(3): 564-573. doi: 10.3969/j.issn.1001-4683.2021.03.002 [45] 范智旎, 万永革, 2020. 应变仪记录的印尼苏门答腊海域7.8级地震的同震信号研究[J]. 大地测量与地球动力学, 40(8): 849-853, 876. [46] 丰成君李滨, 李惠, 等, 2022. 南迦巴瓦地区地应力场估算与构造稳定性探讨[J]. 地质力学学报, 28(6): 919-937. [47] 国家地震局科技监测司, 1995. 地震地形变观测技术[M]. 北京: 地震出版社: 229-250. [48] 何德强, 李亚芳, 佘如昌, 等, 2014. 2012年9月7日彝良5.7、5.6级地震的近场异常[J]. 地震研究, 37(1): 1-8. [49] 李彬, 张文, 文冉, 2022. 陕南特长公路隧道水压致裂法地应力测量结果及工程地质意义分析[J]. 地质力学学报, 28(2): 191-202. [50] 李富珍, 任天翔, 池顺良, 等, 2021. 基于钻孔应变观测资料分析远震造成的动态库仑应力变化[J]. 地球物理学报, 64(6): 1949-1974. [51] 骆鸣津, 顾梦林, 李安印, 等, 1989. 用引潮力进行钻孔的应变-应力的原地标定[J]. 地壳形变与地震, 9(4): 51-62. [52] 潘立宙, 1981. 测量地应力的钻孔变形法[M]//中国地质科学院地质力学研究所, 国家地震局地震地质大队. 地应力测量的原理和应用. 北京: 地质出版社: 1-47. [53] 彭华, 吴珍汉, 马秀敏, 2006. 青藏铁路无人值守地应力综合监测站[J]. 地质力学学报, 12(1): 96-104. [54] 彭华, 马秀敏, 姜景捷, 2008. 山丹地应力监测站体应变仪的地震效应[J]. 地质力学学报, 14(2): 97-108. doi: 10.3969/j.issn.1006-6616.2008.02.001 [55] 彭华, 马秀敏, 姜景捷, 等, 2011a. WFSD地应力台应变特征及其同震效应分析: 以日本Mw9.0级特大地震为例[J]. 地质力学学报, 17(1): 1-14. [56] 彭华, 马秀敏, 姜景捷, 2011b. 差应变法地应力测量: 以汶川地震断裂带科学钻探WFSD-1钻孔为例[J]. 地质力学学报, 17(3): 249-261. [57] 邱泽华, 石耀霖, 欧阳祖熙, 2005a. 四分量钻孔应变观测的实地绝对标定[J]. 地震, 25(3): 27-34. [58] 邱泽华, 石耀霖, 欧阳祖熙, 2005b. 四分量钻孔应变观测的实地相对标定[J]. 大地测量与地球动力学, 25(1): 118-122. [59] 邱泽华, 2010. 中国分量钻孔地应力-应变观测发展重要事件回顾[J]. 大地测量与地球动力学, 30(5): 42-47. [60] 邱泽华, 2017. 钻孔应变观测理论和应用[M]. 北京: 地震出版社. [61] 邱泽华, 唐磊, 赵树贤, 等, 2020. 用应变地震观测求解震源矩张量的基本原理[J]. 地球物理学报, 63(2): 551-561. [62] 全建军, 赖见深, 陈珊桦, 等, 2021. 小陶地震台钻孔体应变仪观测质量与映震能力分析[J]. 华南地震, 41(1): 26-33. [63] 苏恺之, 2003. 地应力相对测量方法[M]//苏恺之, 李海亮, 张钧, 等. 钻孔地应变观测新进展. 北京: 地震出版社: 252-260. [64] 苏恺之, 2018. 四分量钻孔应变仪计算式梳理及震例应用[J]. 地震地磁观测与研究, 39(5): 101-111. [65] 孙艺玫, 查楠, 任雪, 等, 2021. 钻孔体应变仪在地震前兆观测中的理论和应用研究综述[J]. 防灾减灾学报, 37(4): 69-74. doi: 10.13693/j.cnki.cn21-1573.2021.04.011 [66] 吴立恒, 陈征, 李涛, 等, 2010. 电容式倾斜传感器在地壳形变测量中的应用[J]. 传感器与微系统, 29(10): 129-131. [67] 吴立恒, 熊玉珍, 李涛, 等, 2018. 一种RZB分量式钻孔应变仪原地线性标定技术及试验[J]. 大地测量与地球动力学, 38(4): 437-440. [68] 吴立恒, 李宏, 陈征, 等, 2022. 具有量程扩展和原位标定功能的RZB钻孔应变仪传感器稳定性研究[J]. 大地测量与地球动力学, 42(3): 318-320. [69] 郑亚东, 张进江, 张波, 2022. 新世纪构造地质学两大支柱理论: 最大有效力矩准则与变位形分解[J]. 地质力学学报, 28(3): 319-337. doi: 10.12090/j.issn.1006-6616.2021118 [70] 朱思雨, 曹佳文, 丰成君, 等 , 2023. 河北高阳低凸起深部热储层回灌注水诱发断层失稳危险性探讨[J]. 地质力学学报, 29(2): 220-235. -
下载:
下载:
计量
- 文章访问数: 1097
- HTML全文浏览量: 135
- PDF下载量: 65
- 被引次数: 0
