1. 引言
氡观测是国际上公认的地震监测手段之一,也是我国地震地下流体观测台网中重要测项,在地震趋势分析与短临震情研判中发挥着重要作用 [1] - [6] ,主要通过测氡仪对深层地下水(或温泉)中溶解气、逸出气以及断裂带土壤气体中的氡气进行连续观测,捕捉地震 [7] [8] 。氡被岩石、土壤颗粒表面吸附与解附,或溶解于地下水并随地下水迁移,在地震孕育与发生过程中,随着岩石应力与热力状态的变化及深部物质运移,可能会导致岩石氡射气系数改变、氡溶解度变化。但不可否定,气氡观测也存在许多影响观测资料质量的因素 [8] [9] ,如井水流量、脱气系统、观测室温度、湿度、气压等环境因素,水井周边环境的干扰因素等 [10] 。因此用气氡资料进行地震的分析预报,必须对异常开展充分的核实,来判断异常是地震前兆还是干扰 [11] 。本文以九江2井氡值异常成因探讨为例,对气氡异常核实进行探讨,为气氡的异常核实过程及异常分析提供实例,为地下流体的异常核实拓宽思路和方法提供具体的实例。
2. 概述
2.1. 异常概述
2021年4月5日九江2井气氡持续升高,日变化幅度在10 Bq/L左右,6日内累积上升84 Bq/L,氡值上升幅度达43%,如图1所示。初步怀疑可能为脱气装置和气路问题导致,进行气路检查并处理气路中的凝积水等,氡值仍然上升。
2.2. 中心站观测环境概况
九江中心站位于赣北部的庐山西北侧,地理位置为东经116.01˚,北纬29.65˚,海拔110米,构造上属于扬子陆块下扬子地块中部,区域构造位于庐山西北缘边缘剪切带上的夏家–威家岭左行走滑断裂带,该断裂带呈北东走向,延伸约8千米,宽约1~1.5千米,总体倾向320˚~330˚,倾角60˚~70˚,在中–新生代以来长期多次活动,继承了早第四纪时发生的断块抬升运动,并控制第四纪的分布,1911年2月6日,在该断裂带上曾发生九江5.0级地震,该地震距中心站1.5千米 [12] [13] 。九江中心站形变山洞进深45米,覆盖层厚约35米,山洞年温差<1.0℃,日温差<0.02℃,洞体基岩为硅质灰岩。山体相对高度约70米,近EW走向,呈EW缓、NS陡形状,表层植被丰富,山体岩层多见裂隙 [12] [13] 。
2.3. 九江地震观测井参数
九江1井(流体物理量观测)和九江2井(流体化学量观测)均为构造承压自流井孔,总出水量约100 m3/d,主要观测含水层埋深约58~63 m,从历年观测证明,两井系连通管结构,为主要含水层相通,九江中心站井孔参数如表1所示 [14] [15] 。
![](//html.hanspub.org/file/3-1771490x8_hanspub.png?20230221083902425)
Figure 1. Data curve of gas radon anomaly in Well Jiujiang 2
图1. 九江2井气氡异常数据曲线图
![](Images/Table_Tmp.jpg)
Table 1. Parameters of seepage property for sandstone
表1. 九江台井孔参数表
3. 异常调查分析
3.1. 异常初步判断
九江中心站对同井水氡加密取样测试,水氡氡值变幅10%以内,属正常范围内。4月12日22时在3号测点气氡仪后串接一台新BG2015R测氡仪,进行串联对比观测实验,及时收集两台测氡仪氡值数据,发现同步持续上升,说明气氡仪工作状态正常。
气氡异常现场核实:1) 检查九江2号井及井口装置变化情况;2) 检查中心站BG2015R型测氡仪、FD-125氡钍分析仪、SWY-II型水位仪、RTP-1气象三要素仪和新安装BG2015R型测氡仪等仪器工作状态,确认仪器均正常工作;3) 检查测氡仪、气汞和气汞G水路及气路,特别检查气路是否通畅等,经检测均属正常;4) 分析上述仪器近期观测资料,均表现为日变动态正常,水位固体潮形态良好。另外,近半年来水位仪、水氡仪、气氡仪校测资料,表明仪器工作状态良好,观测正常;5) 检查与此次异常相关的工作日志与观测日志;6) 检查中心站外围,未发现有明显变化。
3.2. 中心站观测环境调查
对九江中心站的九江1井及井房、抽水泵房、白龙荡水坑及山洞查看,均未发现其它异常现象。对中心站周边环境进行调查,没发现较大的施工工程开工,未发现有打井、抽水等明显的干扰源。本次气氡异常涉及地下水化学组分的含量变化,其他显性的水位埋深、地倾斜、应变均没有明显的变化,故对中心站周边环境调查只限于一般性的了解,中心站200米外的5727井房井孔没有发现异常,周边村民生活井亦无异常。
3.3. 加密监测实验
除在测点气氡后串接一台新BG2015R测氡仪,进行串联对比实验,再进行以下实验排查工作:
1) 在气氡观测室内安装温湿度计进行温度、湿度等环境监控。
2) 考虑观测系统问题,对同井出水口并联安装一套气氡观测系统,进行同井气氡并联对比实验,进一步验证气氡异常的可靠性。4月19日20时在九江2井房主井管外水氡水样取水点处(接主井管的出水管)安装一个三通,一路水用于水氡采样,一路用于并联新观测系统气氡观测用水,选择脱气装置为自然吸气脱气装置,型号与测点测氡仪脱气装置一致,脱出气体连接一台P2000F测氡仪,开展同井并联对比观测实验。
3) 分别在2井主井管出水口、井房外泄水池、氡比测平台出水口采集三个水样,封装邮寄中国地震局地震预测研究所,进行地下水氢氧同位素等水质测量实验,每3天采集一次水样,共采集4次。
4) 由于地震氡观测仪器检测平台之前一直处于实验状态,在3月下旬停止实验,并把水路管道和恒流箱水排空,4月出气氡出现异常。考虑可能是氡检测平台系统影响,于5月11日开展检测平台抽水实验把管路和恒流箱体处于满管状态,验证气氡异常是否与水路管道挤压有关。
4. 异常分析
4.1. 观测资料对比分析
对中心站所有地球物理测项观测数据对比分析,随着春天雨季到来,3月份雨量以三波逐渐增大的方式到来,水位亦同步逐级升高,最高上升19 cm,同期水温逐渐升高至目前维持在相对高值阶段,而同井观测的水氡稳定变化没有出现异常,水位有一定的同步异常变化;2井气汞和气象三要素,气汞值维持比较低水平,变幅不大,观测室温度和室外环境温度缓慢上升。
气氡与定点形变相关观测资料对比,从气氡异常与倾斜、洞体应变对比看,形变各测项均未有明显的同步异常变化,而在3月17~19日降雨量过程中,形变各测项均有较好的响应。总体而言,除降雨、水位资料可能存在同步异常变化外,其它的地下流体观测、形变观测均未见明显异常 [12] [13] 。
4.2. 观测资料关联性分析
4月5日~12日测点气氡上升,上升幅度达43%,九江1号井水位4月5日~7日上升约3 cm,后至12日平稳变化;降雨量在3月30日至4月4日连续小雨,分别降雨5 mm、6.7 mm、12.6 mm、10.3 mm、2.7、2.4 mm总降雨量合计39.7 mm,4月7日降雨16.5 mm,4月11、12日分别降雨9.3 mm、6 mm;一般水位升高说明井孔流量增大,而本次异常变化与2015年6月24日气氡异常不同之处在于,水位4月7日~12日不再升高而气氡值仍维持固定速率持续上升,因此,难以把此次异常与水位上升确定为直接相关。针对气氡异常与水位存在一定的同步性问题,回顾降雨–水位–气氡三者间多年来统计关系,见图2,具体对应表2。
![](//html.hanspub.org/file/3-1771490x9_hanspub.png?20230221083902425)
Figure 2. Corresponding curves of water level, gas radon anomaly and rainfall in Jiujiang Central Station
图2. 九江中心站水位、气氡异常和降雨量对应曲线图
![](Images/Table_Tmp.jpg)
Table 2. Statistical table of water level, gas radon anomaly and rainfall in Jiujiang Central Station
表2. 九江中心站水位、气氡异常和降雨量统计表
从表2看出2021年04月01日~04月14日(G)降雨量75 mm、水位变化0.089 m,氡值变化率为43.68%,与2015/6/21~6/22(B)有点相似,从九江中心站气氡值变化率、水位和降雨对应图来看,相关性较好。但与A、C、D、E和F时间段对比,降雨量和水位变化均大于G时间段,氡值变化率较低,即使受降雨和水位影响也不会引起这么大的变化幅度和这么长的持续时间,因此难以说明本次氡值异常与水位、降雨有直接影响。
4.3. 化学量对比实验分析
按测氡仪比测相关技术要求,4月12日晚10时,在3测点测氡仪出气口串接另一台贝谷公司生产BG2015R测氡仪,采样率为每10分钟采样一次,二者的观测数据和趋势吻合度较高,变化速率与变化形态基本一致,如图3(a)所示。表明在此观测条件下,3测点测氡仪观测到气氡观测数据是真实的,仪器工作状态正常。按测氡仪比测相关技术要求,同井并联气氡观测数据与3测点测氡仪观测对比如图3(b)所示,从图可以看出并联气氡数据稳定,变化较小,而在线观测气氡先升高后降低,说明可能在线观测系统或环境有干扰。
采集中心站三个位置4次水样邮寄中国地震局地震预测研究所进行水化组分测试,根据测试结果分别做Piper图和Schoeller图,Piper图表明水中阴阳离子组成的常用分析图,根据样点在Piper图上位置,能够指示水的化学组成类型,明确不同水体之间的混合作用 [16] ,如图3(c)所示常量离子的含量较为接近,说明补给源一致。Schoeller图是由一组垂直平行线组成,横坐标上每条线代表一种化学成分,纵坐标为离子含量,用折线将各种离子含量的点连接在一起 [17] 。Schoeller图可直观反映研究区内不同水体之间化学成分变化及分布的相关关系,以及同一水体不同时间的变化特征 [16] ,如图3(d)能够清晰地说明主要离子成分的含量基本接近,是同一补给来源 [17] ,与Piper图的结果一致。图3(c)和图3(d)可以看出除cl-离子有微小变化,其他基本没有变化,进一步说明在线观测气氡可能受观测系统或环境干扰。
![](//html.hanspub.org/file/3-1771490x11_hanspub.png?20230221083902425)
![](//html.hanspub.org/file/3-1771490x10_hanspub.png?20230221083902425)
![](//html.hanspub.org/file/3-1771490x13_hanspub.png?20230221083902425)
![](//html.hanspub.org/file/3-1771490x12_hanspub.png?20230221083902425)
Figure 3. Experimental curve of chemical quantity comparison
图3. 化学量对比实验曲线图
4.4. 环境实验分析
新温湿度环境监控与在线气氡观测在同一观测室内如图4(b)。两者数据如图4(a)所示,可以看出4月28至5月8日气氡出现间断性突跳,观测室内温度也存在同样现象,但气氡稍有滞后现象,可看出气氡跟温度存在一定关联,在开空调使室内温度稳定后,突跳现象逐渐消除。
5月11日开展检测平台抽水实验把水路管路和恒流箱处于水满管状态如图4(d)所示,之后气氡数据开始下降并处于稳定,数据恢复到异常前的数值,如图4(d)所示,说明气氡异常可能水路与恒流箱的排空,而后再恢复流水所造成的。因为排空后再注水用于观测,在注水过程后水路的管道内有一些空间是地下水不可能立即充满,留有一些空间,在地下水流动过程中地下水中的溶解氡会慢慢地释放出来变成逸出气氡,这些逸出气氡就充填到没有被地下水充满的空间中累积,然后就慢慢地参与到地下水中流到脱气装置,而后进入传感器。随后这些空间被地下水慢慢地充满,逸出气氡混合到地下水的量慢慢地达到最高,而后也就慢慢地降低,恢复到正常的状态,水路挤压导致压力缓慢变化导致。
![](//html.hanspub.org/file/3-1771490x14_hanspub.png?20230221083902425)
Figure 4. Corresponding curve of environmental experiment
图4. 环境实验对应曲线图
4.5. 气氡异常的震例分析
九江中心站气氡资料做回顾性震例分析得到以中心站100公里以内有两次ML4.0级以上地震,2011年9月10日江西九江ML4.9级地震震前表现为较好的异常对应。2015年6月17日江西九江ML4.2级地震未记录到异常;100~200公里以内有安徽安庆ML4.8、桐城ML4.1、霍山ML4.4、江西浮梁ML4.2,其中桐城ML4.1级地震似有一定异常显示,其它地震均未见明显异常。
本次气氡异常与2011年9月10日前的异常现象有一定的相似性,后者是在2011年8月10日起逐级升高,至19日升高至高值平台后持续高值约12天后,氡值下降,之后一周内发震,本次气氡从5日起逐级升高目前仍在升高,从同井气氡和水氡加密实验得出仪器真实值保持稳定,排除是地震前兆。
5. 讨论和结论
经过对中心站观测仪器检查、对比观测实验、中心站外围2公里范围走访调查等现场工作,就进行数据资料分析,认为2021年4月以来九江气氡异常不是地震前兆异常,引起气氡异常主要成因为以下方面:
1) 由于2井3号测点气氡观测水路和氡检测平台水路并排进行观测室,把氡平台水路排空,氡值出现上升,恢复氡平台水路和恒流箱处于水满管状态时,气氡数值下降,并恢复到异常前的氡值,水路与恒流箱的排空,后再恢复流水所造成的,水路挤压导致脱气压力缓慢变化导致。
2) 九江2井气氡变化与降雨、水位变化存在一定影响和关联,可能性较小。
3) 九江2井3号测点气氡观测系统在观测室内,气氡与观测室内温度变化存在一定关联,在开空调使室内温度稳定后,突跳现象逐渐消除。
通过九江2井气氡的异常核实,认识到异常核实除了要调查周边环境干扰、观测仪器是否正常工作、了解以往震例及以往出现类似异常的原因外,还应设计恰当的对比观测实验和设计环境实验来验证推论,分析气氡变化及特征,必要时还应充分应用异常测项本身的物理、化学性质及水井的基础资料进行综合判定,为气氡异常核实跟踪提供实验支撑,从而提高气氡观测效能,并介绍异常核实的步骤和经验,供他人讨论和借鉴。
基金项目
地震科技星火计划项目(XH21017Y)、中国地震局三结合课题(3JH-20221402)共同资助。
NOTES
*通讯作者。