-
构造应力变化可引发裂隙张开、断层滑动,地震断裂带在地震孕育与发生过程中,为深部流体提供了运移通道[1],沿断裂带向地表运移的气体地球化学异常可以反映断裂带活动情况[2]及地震活动性[3-4]。苏联自1966年开展了地球化学断层气方面研究工作以来取得了几次较好的震例[5];日本自1973年在东海、伊豆半岛等地布设观测点,获得了1995年神户7.2级地震的震例[6];1998年,Ciotoli等[7]在意大利富希纳盆地进行的断层气探测结果表明,在1915年发生7.0级地震的断裂带上同步显示出CO2浓度异常。1993年王基华等[8]在北京昌平白浮进行了CO2浓度的观测,研究了其动态特征与地温的关系;林元武等[9-10]研制了断层气CO2快速测定方法,并分析了怀来后郝窑观测点7年的CO2资料,基于9个4级以上地震获得了CO2的前兆异常特征,1998年张北6.2级地震前47天怀来后郝窑断层土壤气体CO2浓度值异常升高,最高值是背景值的10倍左右;2011年任佳等[11]分别在赤城地震台、后郝窑水化观测站和阳原地震台进行CO2实验观测,证明了二氧化碳红外气体传感技术在二氧化碳断层气测量方面的可行性。1976年唐山7.8级地震前,天津棉4井土壤气体CO2浓度值急剧升高数倍,震后在唐山断裂带上进行了CO2跨断裂探测,发现跨断裂CO2浓度值明显高于非断裂地段[12]。川西地区的断层土壤气体CO2观测点在区域内4.0级以上地震短期前兆异常,最高异常峰值是背景值的1.7~10倍[13]。综上,已有观测结果显示CO2浓度值在中强地震前会出现大幅度的升高,可能反映区域构造应力变化[14-16]。因此,断层土壤气体CO2浓度观测越来越受到国内地震研究者的重视[17-20]。
目前,国内主要利用气相色谱仪[21]、CO2浓度含量快速测定法[19]和二氧化碳红外气体传感技术进行二氧化碳气体测试分析。气相色谱仪精确度高,需要采集气体样品后送实验室测试,不能实时测量土壤气体浓度分布;CO2浓度含量快速测定法,采样率为1次/日,观测数据为24小时累积数据,受测量方法限制,难以精确测定浓度梯度变化,数据量不能满足监测短临异常变化的需求;二氧化碳红外气体传感技术连续观测可以产出采样率1次/分钟的观测数据[22],可观测到震前、震时、震后的连续动态变化。本文通过在唐山地区建设CO2浓度连续观测台阵,分析了CO2年变动态受到气温、地温、气压、降雨量等干扰因素的影响,并初步总结了该地区CO2基本动态特征和震前异常特征,为今后CO2浓度连续观测台阵建设和数据分析提供参考。
-
唐山地区地质构造复杂,位于Ⅱ级构造单元华北地台东北隅,北部为燕山隆起带,南部为华北坳陷区,Ⅲ级构造单元有蓟县坳陷,山海关隆起,本区基底经燕山运动褶皱成山,形成褶皱和断裂,称为唐山市基底构造骨架[23]。
研究区内主要分布2组断裂带:一组NE向的丰台野鸡坨断裂、宁河昌黎断裂和唐山断裂带,另一组为NW向的滦县乐亭断裂及蓟运河断裂等[24]。该区断裂活动以张性为主,地震活动频繁,是全国中强地震活动的重要区域。历史上发生的较大地震为1976年7月28日唐山7.8级地震。2018年以来唐山地区小震活动水平较高,共发生M1.0以上地震169次,最大地震为2020年7月12日唐山市古冶区M5.1地震(图1)。
-
为研究唐山地区断层土壤气体CO2与断裂带活动及地震情况,2018年6月进行了土壤气体CO2浓度观测背景勘选。隐伏断层的勘选测线均垂直于断层走向布设,其中,靳庄测线方向东西设置12个观测点(图2),CO2平均浓度8 700 ppm,变化范围1 400~18 600 ppm;赵滩测线方向东西设置8个观测点,CO2平均浓度164 ppm,变化范围7 100~32 500 ppm;北泊测线方向东北设置19个观测点,CO2平均浓度17 200 ppm,变化范围900~89 100 ppm;孙家楼测线方向东西设置15个观测点,CO2平均浓度11 500 ppm,变化范围5 300~29 800 ppm;陡河台测线方向东西设置15个观测点,CO2平均浓度7 300 ppm,变化范围2 400~20 300 ppm。滦县台、林西矿、赵各庄和马家沟均位于活动断层基岩破碎带上,其连续观测点均建在前期背景勘选测线中CO2最高值的位置(表1)。
表 1 唐山地区CO2连续观测点概况
台站 经度/(°E) 纬度/(°N) 断裂 测点情况 埋深/m 观测开始时间 靳庄 118.55 39.41 宁河昌黎断裂 地表松散覆盖土壤中 10 2018-08-15 陡河台 118.30 39.73 陡河断裂 断层基岩破碎带中 10 2018-08-01 滦县台 118.67 39.74 滦县乐亭断裂(北段) 断层基岩破碎带中 10 2018-08-17 赵滩 118.93 39.32 滦县乐亭断裂(南段) 断层破碎带的断层泥中 10 2018-08-15 林西矿 118.45 39.72 唐山古冶断裂(北段) 地表松散覆盖土壤中 10 2018-08-16 北泊 118.01 39.31 唐山古冶断裂(南段) 断层破碎带的断层泥中 10 2018-08-01 孙家楼 118.18 39.56 唐山古冶断裂(中段) 地表松散堆积物中 10 2018-08-01 赵各庄 118.41 39.76 巍山长山南坡断裂(北段) 地表松散覆盖土壤中 10 2018-08-15 马家沟 118.26 39.69 巍山长山南坡断裂(南段) 地表松散覆盖土壤中 10 2018-08-08 -
断层土壤气体CO2连续观测仪器型号为RZW-3A型,具有长期连续、稳定观测的特点(表2)。观测系统采用模块化设计,使用中国地震局数字化地震前兆观测仪器通信协议。仪器安装前为满足野外观测需要依据国家计量检定规程JJG635-2011《一氧化碳、二氧化碳红外气体分析器》对仪器进行检测(表3)。校准结果见表4 。
表 2 RZW - 3A型CO2测量仪技术指标
测量量程 0~10 000 ppm或0~50 000 ppm 分辨率 0~10 000 ppm:10 ppm;0~50 000 ppm:100 ppm 零点重复性 ±1% FS (20℃) 长期零漂 <1 %/月 工作温度 −20 ℃~+70 ℃ 采样率 1次/分钟 供电 50~80 W的太阳能电池供电 表 3 RZW - 3A型CO2测量仪检测结果
检测项目 技术要求 0012# 0025# 0018# 引用误差 二级 ±2% FS 1.3 −1.5 −1.9 重复性 二级 ≤1.0% 0.20 0.35 0.41 响应时间 ≤60 s 20 26 28 漂移 二级 ≤ ±2% FS 0 0.2 0.2 0.2 0.2 0.2 干扰误差 二级 ≤ ±2% FS 0 0.2 −0.2 0.2 −0.2 0.2 表 4 RZW - 3A型CO2测量仪校准结果
标准值/
(mol/mol)测量值/
(μmol/mol)测量值测量结果的不确定度
(k=2)0.201% 2 400 Urel=1.2% 0.405% 4 500 Urel=1.2% 0.803% 8 200 Urel=1.2% 本文使用的CO2连续观测仪气体采集方式为扩散式(逸出式),与人工观测的二氧化碳测试管所放置环境类似。即先在土壤中打10 m深孔,然后埋入底部约2 m长度花孔状的pvc管,将二氧化碳传感器直接竖直放置于花孔位置处。二氧化碳传感器由封装装置、采集装置和信号传输电路3部分组成。整个观测系统采用太阳能供电、无线传输、数据保存于云端服务器(图3)。
-
根据波义尔-马里奥气体定律,CO2浓度变化受温度影响。唐山地区断层土壤气体CO2与气温、地温呈正相关关系。地温相对于气温与CO2浓度的变化趋势一致性更高,相关系数分别为0.68和0.60。根据已有前人对气温影响深度的研究,一般超过2 m就可以忽略气温日变化的影响[25],唐山地区连续观测孔深度10 m,故可以忽略气温日变化的影响。CO2浓度与气压相关性系数为-0.43。
以陡河台为例,由降雨量、气压、气温、地温和CO2浓度值曲线图的h1、h2、h3时段中可看出CO2受降雨影响较大(图4)。2018年、2019年和2020年丰水期时段相较于枯水期CO2浓度值曲线波动幅度大幅减小,曲线上升和下降趋势减缓。在2020年丰水期长、全年降雨量大的情况下,CO2浓度值在丰水期时段波动幅度明显减小,上升趋势减缓趋于平稳,枯水期时段曲线波动幅度明显增大,整体变化趋势明显。h2时段当连续多天强降雨后CO2浓度值曲线出现明显下降,分析认为降雨导致土壤中的空隙被雨水填充,而CO2在水中的溶解度较大 [16],沿断层逸出至地表的CO2部分溶于水,导致浓度降低。综上,研究区测点断层土壤气体CO2浓度值受降雨影响较大,与地温正相关,与气压负相关。
为进一步分析短时降雨对CO2浓度测值的影响,绘制了唐山地区2020年8月份丰水期降雨较多时段CO2连续观测台阵测值变化曲线(图5)。全部观测点降雨过后浓度值较降雨前都出现明显下降,且降雨量越大浓度值下降幅度越大,其中马家沟、靳庄、滦县台和陡河台浓度值变化有滞后性。图5a中马家沟、靳庄和滦县台随降雨量增大CO2浓度值阶跃上升;图5b赵滩、陡河台整体平稳,随降雨量增大CO2浓度测值快速下降,降雨结束后测值快速上升;图5c孙家楼、林西矿CO2浓度整体呈上升趋势,受降雨影响缓慢上升;图5d赵各庄、北泊受降雨影响CO2浓度值波动较大。图5中,2020年8月1—31日整月CO2浓度值与降雨量相关系数为北泊:−0.18;孙家楼:−0.03;马家沟:−0.27;赵各庄:−0.37;赵滩:−0.35;林西矿:−0.06;陡河台:−0.18;滦县台:−0.01;靳庄:−0.10。进一步选取2020年8月9—19日的连续降雨时段内相关系数分别为北泊:−0.56;孙家楼:−0.34;马家沟:−0.35;赵各庄:−0.37;赵滩:−0.38;林西矿:−0.33;陡河台:−0.44;滦县台:−0.48;靳庄:−0.64。不同时间段的降雨量与浓度相关性计算结果差异表明,浓度观测值会受降雨量影响。
进一步使用基于Gamma分布函数建立的组合水箱模型来分析降雨量对CO2浓度值变化的滞后延迟效果[26],映证了CO2浓度值与降雨量的负相关性(图6)。
-
断层气体CO2多年资料动态显示,年动态多为夏高冬低型[13],唐山地区断层土壤气体9个连续观测点呈不同程度的夏高冬低年变动态(图7)。其中,赵各庄CO2浓度曲线形态表现为脉冲式波动,总体呈上升趋势;靳庄CO2浓度曲线形态表现为上升趋势;北泊CO2浓度曲线形态表现为缓升缓降,总体呈上升趋势;孙家楼CO2浓度曲线形态表现为缓升缓降,总体呈下降趋势;马家沟CO2浓度曲线形态表现为前期大幅下降,后期缓升缓降,总体呈下降趋势;赵滩CO2浓度曲线形态表现为下降趋势;林西矿CO2浓度曲线形态表现为振荡,总体呈下降趋势;陡河台、滦县台CO2浓度曲线形态表现为均值附近上下振荡,总体趋势平稳。
进一步分析2018、2019及2020年的连续观测数据(表5),发现处于宁河昌黎断裂上的靳庄断层土壤气体年平均值连续两年上升,而位于唐山古冶断裂上的北泊、孙家楼、林西矿,长山巍山断裂的马家沟、赵各庄、唐山断裂的陡河台和滦县乐亭断裂北段的滦县台CO2平均值没有明显连续变化,但滦县乐亭断裂南段的赵滩两年平均值持续下降。浓度均方差结果发现,靳庄、北泊、赵各庄的CO2浓度值年变化范围逐年增大。由均方差与平均值的比值可知,北泊和马家沟2个观测点明显高于其余7个观测点。
表 5 各观测点数据分析
台站 年份 靳庄 陡河台 滦县台 赵滩 林西矿 北泊 孙家楼 赵各庄 马家沟 最大值(月
份)~最小值
(月份)/
ppm2020 79 000(4)~
98 000(11)34 200(3)~
101 300(11)32 060(2)~
79 660(8)39 000(6)~
50 800(9)10 100(5)~
19 700(10)2 700(3)~
39 900(9)16 000(5)~
26 600(10)44 800(5)~
79 300(10)2 000(3)~
11 800(9)2019 67 000(6)~
102 200(11)35 800(2)~
109 400(11)30 460(2)~
94 060(11)41 700(3)~
55 200(8)11 000(4)~
24 800(10)2 800(2)~
29 800(10)17 900(4)~
31 700(10)44 600(3)~
68 600(10)2 700(3)~
16 200(8)2018 72 200(8)~
84 000(9)34 700(12)~
89 200(10)36 060(12)~
73 860(10)46 900(12)~
55 942(9)20 500(12)~
29 900(9)4 700(12)~
30 100(9)27 600(12)~
37 800(9)45 800(11)~
59 500(9)7 900(12)~
67 000(8)均方差/平均值 2020 0.06 0.22 0.24 0.04 0.22 0.88 0.14 0.07 0.40 2019 0.15 0.24 0.26 0.07 0.26 0.68 0.19 0.07 0.41 2018 0.03 0.16 0.15 0.04 0.10 0.49 0.08 0.03 0.61 均方差 2020 5 762 15 624 13 759 2 187 3 259 11 249 3 103 3 496 2 086 2019 12 821 18 981 17 077 3 464 4 629 8 552 4 625 3 539 3 441 2018 2 925 11 549 9 190 2 328 2 711 7 962 3 016 1 689 13 677 平均值/ppm 2020 88 186 70 932 57 205 44 388 14 304 12 769 20 857 47 720 5 178 2019 81 090 76 573 63 796 47 294 17 363 12 544 24 280 48 210 8 317 2018 77 301 71 410 58 194 52 755 26 732 15 940 34 085 47 210 22 336 唐山地区断层土壤气体CO2连续观测台阵浓度空间分布表现为中部最低,西北部高于东南部,唐山古冶断裂整体断层土壤气体CO2浓度值呈缓升缓降趋势;长山巍山断裂北段振荡,南段缓降;滦县乐亭断裂北段缓升,南段缓升缓降;唐山南部的宁河昌黎断裂呈缓升趋势,陡河断裂整体呈震荡形态。空间上年变可能反映了区域的断裂破碎程度差异,即唐山-古冶断裂整体破碎程度高,多年连续变化规律不明显,滦县乐亭断裂南段破碎程度低于北段,宁河昌黎断裂的浓度测值最高,可能反映了该断裂较为发育。
-
在2018—2020年唐山地区断层土壤气体CO2浓度连续观测期间,唐山地区共发生M3以上地震3次,根据观测时段内的区域地震活动计算对应能量释放情况:
$$ {\rm{Log}} E=9.9+1.9M-0.024M ^{2} $$ (1) 选取地震前后各一个月内数据明显变化的台站分析其与地震可能的关系(图8)。
-
2018年8月5日唐山M3.3地震前3天,赵各庄、北泊观测数据出现连续大幅上升,最高异常峰值是背景值的1.6倍,震后3天恢复正常值(图8b)。同时段内靳庄观测数据平稳,无明显变化。
-
2019年12月5日唐山M4.5地震前8天,赵各庄观测数据出现大幅度阶跃式上升异常,整体降趋势最高异常峰值是背景值的1.6倍,震后7天恢复正常值(图8c)。靳庄震前无明显变化,震后出现小幅度震荡上升,整体呈下降趋势。同时段内北泊观测数据整体平稳,无明显波动。
-
2020年7月12日唐山M5.1地震前20天内赵各庄、靳庄、北泊观测点CO2浓度值先后出现不同程度的震荡上升,最高异常峰值是背景值的3.8倍,震后14天恢复正常值(图8d)。震前多个观测点记录到异常变化,CO2变化最大值为背景值的3.8倍,而孙家楼、滦县台、陡河台、林西矿则无显著变化(图7~8)。这可能与观测点位置、各条断裂破碎程度及震中区构造活动有关。
根据震级越大、震中距越近、CO2异常幅度越大[27],则地震一般发生在异常峰值转折下降过程中的一般规律[13],2018年8月至2020年12月期间,唐山地区断层土壤气体CO2浓度连续观测台阵共记录6次显著异常,为平缓波动基础上的阶跃上升、震荡式上升和突跳3种形态,其中有3次异常结束后的一段时间内发生了地震。根据地震预测指标理论方法计算了唐山地区断层土壤气体3次地震的R值,其中,赵各庄观测点的R(0.86)>R0(0.71),通过R值检验。其余观测点R<R0,但是均大于0,具有参考意义[28]。
-
通过分析2018—2020年唐山地区断层土壤气体CO2连续观测台阵观测结果,得到以下认识。
1)通过2年观测实验发现,唐山CO2连续观测台阵适合于偏远地区或干旱少雨的无人值守观测点及需要加密观测的震情监视区。因为安装简单快速、低成本及连续稳定等特点,具有应用推广前景。CO2连续观测仪在观测过程中有偶尔丢失单个数据及数据传输慢的现象,通过排查发现可能是电源、数据采集和信号传输问题所致。因此,连续观测仪的电源、采集和传输模块需进一步优化。
2)唐山地区断层土壤气体CO2浓度值受降雨量和地温影响较大,与气压负相关。CO2浓度总体呈现夏高冬低的年变动态,间接反映了与气温、地温变化的正相关性。
3)唐山CO2连续观测台阵2年观测实验时段内部分观测点观测到CO2异常,可能是唐山及邻近地区M3.3以上地震的前兆异常,异常出现时间多为地震前1个月内,范围是背景值的1.6~3.8倍。
Analysis of Continuous Observation Array Data of Fault Soil Gas CO2 in Tangshan Area
-
摘要: 利用唐山地区9个断层土壤气体CO2浓度观测点构成的连续观测台阵2018—2020年观测数据,分析发现该地区断层土壤气体CO2年变动态主要受地温及降雨量影响;在2018年唐山古冶M3.3地震、2019年唐山丰南M4.5地震和2020年唐山古冶M5.1地震前一个月左右,赵各庄、靳庄、北泊观测点记录到不同程度的CO2浓度异常变化,异常形态主要有阶跃上升、震荡式上升和突跳3种,异常峰值最大为背景值的3.8倍。连续观测数据的异常变化可能包含地震前兆信息。Abstract: Based on the observation data from 2018 to 2020 of nine continuous observation stations constructed in Tangshan area, it is found that the annual variation of fault soil gas CO2 is mainly affected by ground temperature and rainfall. About one month before the 2018 Tangshan Guye M3.3 earthquake, 2019 Tangshan Fengnan M4.5 earthquake and 2020 Tangshan Guye M5.1 earthquake, different degrees of abnormal changes were recorded in Zhaogezhuang, Jinzhuang and Beipo. The abnormal patterns mainly included step rise, oscillation rise and sudden jump. The maximum abnormal peak value is 3.8 times of the background value. Abnormal changes in continuous observation data may contain earthquake precursory information.
-
Key words:
- Tangshan area /
- fault soil gas /
- CO2 concentration /
- interference factors /
- precursor information
-
表 1 唐山地区CO2连续观测点概况
台站 经度/(°E) 纬度/(°N) 断裂 测点情况 埋深/m 观测开始时间 靳庄 118.55 39.41 宁河昌黎断裂 地表松散覆盖土壤中 10 2018-08-15 陡河台 118.30 39.73 陡河断裂 断层基岩破碎带中 10 2018-08-01 滦县台 118.67 39.74 滦县乐亭断裂(北段) 断层基岩破碎带中 10 2018-08-17 赵滩 118.93 39.32 滦县乐亭断裂(南段) 断层破碎带的断层泥中 10 2018-08-15 林西矿 118.45 39.72 唐山古冶断裂(北段) 地表松散覆盖土壤中 10 2018-08-16 北泊 118.01 39.31 唐山古冶断裂(南段) 断层破碎带的断层泥中 10 2018-08-01 孙家楼 118.18 39.56 唐山古冶断裂(中段) 地表松散堆积物中 10 2018-08-01 赵各庄 118.41 39.76 巍山长山南坡断裂(北段) 地表松散覆盖土壤中 10 2018-08-15 马家沟 118.26 39.69 巍山长山南坡断裂(南段) 地表松散覆盖土壤中 10 2018-08-08 表 2 RZW - 3A型CO2测量仪技术指标
测量量程 0~10 000 ppm或0~50 000 ppm 分辨率 0~10 000 ppm:10 ppm;0~50 000 ppm:100 ppm 零点重复性 ±1% FS (20℃) 长期零漂 <1 %/月 工作温度 −20 ℃~+70 ℃ 采样率 1次/分钟 供电 50~80 W的太阳能电池供电 表 3 RZW - 3A型CO2测量仪检测结果
检测项目 技术要求 0012# 0025# 0018# 引用误差 二级 ±2% FS 1.3 −1.5 −1.9 重复性 二级 ≤1.0% 0.20 0.35 0.41 响应时间 ≤60 s 20 26 28 漂移 二级 ≤ ±2% FS 0 0.2 0.2 0.2 0.2 0.2 干扰误差 二级 ≤ ±2% FS 0 0.2 −0.2 0.2 −0.2 0.2 表 4 RZW - 3A型CO2测量仪校准结果
标准值/
(mol/mol)测量值/
(μmol/mol)测量值测量结果的不确定度
(k=2)0.201% 2 400 Urel=1.2% 0.405% 4 500 Urel=1.2% 0.803% 8 200 Urel=1.2% 表 5 各观测点数据分析
台站 年份 靳庄 陡河台 滦县台 赵滩 林西矿 北泊 孙家楼 赵各庄 马家沟 最大值(月
份)~最小值
(月份)/
ppm2020 79 000(4)~
98 000(11)34 200(3)~
101 300(11)32 060(2)~
79 660(8)39 000(6)~
50 800(9)10 100(5)~
19 700(10)2 700(3)~
39 900(9)16 000(5)~
26 600(10)44 800(5)~
79 300(10)2 000(3)~
11 800(9)2019 67 000(6)~
102 200(11)35 800(2)~
109 400(11)30 460(2)~
94 060(11)41 700(3)~
55 200(8)11 000(4)~
24 800(10)2 800(2)~
29 800(10)17 900(4)~
31 700(10)44 600(3)~
68 600(10)2 700(3)~
16 200(8)2018 72 200(8)~
84 000(9)34 700(12)~
89 200(10)36 060(12)~
73 860(10)46 900(12)~
55 942(9)20 500(12)~
29 900(9)4 700(12)~
30 100(9)27 600(12)~
37 800(9)45 800(11)~
59 500(9)7 900(12)~
67 000(8)均方差/平均值 2020 0.06 0.22 0.24 0.04 0.22 0.88 0.14 0.07 0.40 2019 0.15 0.24 0.26 0.07 0.26 0.68 0.19 0.07 0.41 2018 0.03 0.16 0.15 0.04 0.10 0.49 0.08 0.03 0.61 均方差 2020 5 762 15 624 13 759 2 187 3 259 11 249 3 103 3 496 2 086 2019 12 821 18 981 17 077 3 464 4 629 8 552 4 625 3 539 3 441 2018 2 925 11 549 9 190 2 328 2 711 7 962 3 016 1 689 13 677 平均值/ppm 2020 88 186 70 932 57 205 44 388 14 304 12 769 20 857 47 720 5 178 2019 81 090 76 573 63 796 47 294 17 363 12 544 24 280 48 210 8 317 2018 77 301 71 410 58 194 52 755 26 732 15 940 34 085 47 210 22 336 -
[1] Lombardi S, Voltattorni N. Rn, He and CO2 soil gas geochemistry for the study of active and inactive faults[J]. Applied Geochemistry, 2010, 25(8): 1206-1220. doi: 10.1016/j.apgeochem.2010.05.006 [2] King C Y. Radon emanation on San Andreas fault[J]. Nature, 1978, 271(5645): 516-519. doi: 10.1038/271516a0 [3] Seminsky K Z, Demberel S. The first estimations of soil-radon activity near faults in Central Mongolia[J]. Radiation Measurements, 2013, 49: 19-34. doi: 10.1016/j.radmeas.2012.12.013 [4] 王江, 李营, 陈志. 口泉断裂断层气地球化学变化特征及断层活动性[J]. 地震, 2017, 37(1): 39-51. doi: 10.3969/j.issn.1000-3274.2017.01.005 [5] 朱旭, 王子影, 张宇. 红外传感器用于断层气体测量[J]. 大地测量与地球动力学, 2011, 31(5): 156-159. doi: 10.14075/j.jgg.2011.05.033 [6] Tsun U, 常丕兴. 日本神户地震时地下水的前兆化学变化[J]. 国外地质与勘测, 1996(4): 31-32,36. [7] Ciotoli G, Guerra M, Lombardi S, et al. Soil gas survey for tracing seismogenic faults: a case study in the Fucino Basin, central Italy[J]. Journal of Geophysical Research:Solid Earth, 1998, 103(B10): 23781-23794. doi: 10.1029/98JB01553 [8] 国家地震局地质研究所. 首都圈地区1993年断层气监测及1994年震情趋势分析[M]. 北京: 地震出版社, 1993. [9] 林元武, 王基华, 高松升. 断层气CO2测定新方法与张北-尚义6.2级地震预报[J]. 地震, 1998, 18(4): 353-357. [10] 林元武, 刘五洲, 王基华, 等. 张北-尚义地震现场CO2测量与震后趋势判断[J]. 地震地质, 1998, 20(2): 117-121. [11] 任佳, 张晓刚, 王长江, 等. 数字化二氧化碳观测仪的研制和观测[J]. 地震研究, 2013, 36(4): 536-539. doi: 10.3969/j.issn.1000-0666.2013.04.020 [12] 车用太, 鱼金子. 地震地下流体学[M]. 北京: 气象出版社, 2006. [13] 官致君, 闻学泽, 程万正. 断层气CO2快速测定方法在四川地震预测中的应用[J]. 地震研究, 2003, 26(S1): 118-125. [14] 高清武, 范树全. 夏垫断裂断层气在地震前的异常反映[J]. 地震, 1992(6): 73-76. [15] 鱼金子, 车用太, 刘五洲. 地壳中的CO2及其释放与地震短临预测[J]. 国际地震动态, 1998(8): 8-12. [16] Mazzini A, Svensen H, Akhmanov G G, et al. Triggering and dynamic evolution of the LUSI mud volcano, Indonesia[J]. Earth and Planetary Science Letters, 2007, 261(3/4): 375-388. [17] 高小其, 王海涛, 赵建政, 等. 地震危险区短临跟踪的一种有效观测手段-地下岩土气体CO2观测在新疆地区的应用[J]. 华北地震科学, 2004, 22(1): 31-37. doi: 10.3969/j.issn.1003-1375.2004.01.007 [18] 车用太, 鱼金子, 刘五洲. 地壳放气动态监测与张北-尚义MS6.2级地震预报[J]. 地质论评, 1999, 45(1): 59-65. doi: 10.3321/j.issn:0371-5736.1999.01.010 [19] 高小其, 王道, 许秋龙, 等. CO2快速测定方法的应用研究[J]. 内陆地震, 2002, 16(1): 76-83. doi: 10.3969/j.issn.1001-8956.2002.01.012 [20] 谷茂龙, 东得淼. 宁河台逸出气二氧化碳观测数据分析[J]. 甘肃科技, 2019, 35(23): 42-43. [21] 朱方保, 李敏莉, 束克明, 等. SP-2304A型气相色谱仪测定气体中注意的几个问题[J]. 地震地磁观测与研究, 2003, 24(5): 120-123. doi: 10.3969/j.issn.1003-3246.2003.05.025 [22] 张晓刚, 任佳, 冯振玮, 等. 怀4井数字化溶解气CO2观测资料分析[J]. 华北地震科学, 2015, 33(3): 68-72. doi: 10.3969/j.issn.1003-1375.2015.03.013 [23] 雷明堂, 蒋小珍, 李瑜. 岩溶塌陷预测评价系统及其应用——以唐山岩溶塌陷为例[J]. 中国岩溶, 1997, 16(2): 97-104. [24] 刘亢, 张学文, 李津津, 等. 唐山地震震害分布与地质构造关系探讨[J]. 灾害学, 2014, 29(2): 25-28. [25] 金旭, 陈晓冬, 管彦武. 气候变化对浅层地温测量影响的改正[J]. 地球学报, 2004, 25(5): 579-582. doi: 10.3321/j.issn:1006-3021.2004.05.016 [26] 王旭升, 王广才, 董建楠. 断裂带地下水位的降雨动态模型及异常识别[J]. 地震学报, 2010, 32(5): 570-578. doi: 10.3969/j.issn.0253-3782.2010.05.006 [27] 上官志冠. 地热流体溶解CO2总量参数的地震前兆特征[J]. 科学通报, 1995, 40(9): 814-817. doi: 10.3321/j.issn:0023-074X.1995.09.015 [28] 唐杰, 张素欣, 张子广, 等. 怀4井水汞短临预测指标研究[J]. 山西地震, 2021(3): 13-17. doi: 10.3969/j.issn.1000-6265.2021.03.004 -