2023, Vol. 44
表1(Table 1)
1976年7月28日唐山发生7.8级地震，震后约3.5 h、15.0 h，分别在蓟运河断裂和滦县—乐亭断裂发生了宁河6.2级地震，滦县—乐亭断裂发生了滦县7.1级地震（张素欣等，2017）, 唐山7.8级地震及余震造成24万多人死亡，43万多人受伤，给人民的生命、财产带来了巨大的损失。唐山7.8级主震后余震一直持续至今，这使得该区域与主震前相比至今保持着较高的地震活动性（仲秋等，2012）。区域动力仍在持续，在区域应力场的作用下，唐山中强地震活动频繁。2000年以来，唐山地震老震区共发生ML≥4.0地震23次，其中，ML 4.0—4.9地震21次；ML 5.0—5.9地震2次，现今唐山地区地震活动具有继承性，主要受唐山地震老震区的构造环境和断裂模式控制（刘燕翔等，2021）。采用有效的手段来追踪断裂的活动状态及余震的活动性，对唐山地区地震危险性评价有重要意义。
1 研究背景气体赋存在地球各圈层，断裂带是地球深部气体在壳幔活动过程中向上逸出的良好通道（杜建国等，1994），监测断裂带附近气体逸出情况，是探测地震活动的重要手段之一，可为地震预测提供一定依据。近年来，随着断层土壤气观测技术的不断完善与日趋成熟，利用断层土壤气确定断层特性及地震活动等方面的研究较多。李营等（2009）利用延怀盆地422个测点的土壤气测量结果，得出土壤气中气体组分的地球化学背景特征与区域地震地质特征有关，其在一定程度上反映了研究区地震活动的背景。盛艳蕊等（2015）通过在怀来地区新保安—沙城断裂测量了断层土壤气H2、He、CO2、Rn、Hg浓度及CO2、Rn、Hg的通量，证实该断裂带有较高的脱气强度，与区域地震活动性间有一定的关系。郑海刚等（2016）在郯庐断裂带安徽段布设了4条跨断层土壤气测量剖面，测量土壤气Rn、Hg、CO2的浓度，结果表明气体在断裂带附近较富集，对断层位置有一定指示作用。王喜龙等（2017）通过首都圈地区活动断裂土壤气测量，得出首都圈地区土壤气Rn、Hg、CO2的浓度和通量均呈东高西低的变化趋势，与首都圈地区由西至东应力水平增加、地壳逐渐减薄、沉积层厚度增大、地震活动逐渐增强等趋势相对应。蒋雨函等（2022）采用断层土壤气测量方法，揭示了研究区内库松木契克山前断裂、那拉提断裂、霍尔果斯—吐谷鲁断裂具有较强的构造活动性。王江等(2022, 2023)对河北平原的隐伏断裂带使用土壤气浓度观测数据计算了浓度强度，并结合垂直断裂剖面土壤气浓度、通量分析了隐伏断裂分段活动性，发现平原隐伏区土壤气空间分布主要受断裂展布控制。
地下气体Rn、Hg、CO2等组分能够客观地、灵敏地反映地壳应力状态和地震活动（王喜龙等，2017）。其中，断层土壤气CO2是地球内部生成的众多流体组分中最有可能大量迁移至地表并在地表某点集中释放的气体之一，监测其浓度可作为一种有效的手段来追踪断裂活动状态、反映地震活动性（王江等，2017；杨江等，2019；盛艳蕊等，2020）。断层气CO2浓度在一些MS≥4.0地震之前有明显的异常反应，尤其是在MS≥5.0地震前出现的异常较显著，这些异常不仅表现出信噪比高，而且多出现在震前半个月到1个半月的时段内，是典型的短临异常，对破坏性地震短临预报阶段发震时间的判定具有重要意义（林元武等，1998）。1998年在环渤海地区沿活动断层布设的断层CO2浓度观测点，在这些测点CO2浓度变化较好地对应了1998年6月8日北黄海4.4级地震、1998年8月15日渤海4.2级地震（耿杰等，2000）。2017年8月9日精河MS 6.6地震后，朱成英等（2022）通过断层土壤气测量，发现土壤气CO2浓度在空间上有明显的分段特征，在主震、余震震中集中区附近达到最高值，进一步梳理后发现，土壤气CO2浓度空间分布与震后地震烈度空间分布趋势基本一致。本文研究区唐山地震老震区是由NE向宁河断裂、丰台—野鸡坨断裂、唐山断裂及NW向滦县断裂—乐亭断裂、蓟运河断裂围限的菱形块体，且中、新生代活动强烈（张肇诚，1990）。1976年唐山7.8级地震后，中小地震活动频繁，唐山地震老震区一直也是华北地区的重点监视区。2018年在唐山地震危险区内历史发震断层及多条重点监视断层布设了10个土壤气CO2浓度连续观测台阵（表 1、图 1）。通过分析土壤气CO2浓度变化特征、与地温间的关系及对本区域的映震能力，可判断唐山断裂带的活动性和地震活动，进而对该区域地震危险性进行评价。
|
表 1 唐山断裂带土壤气CO2台阵观测环境 Table 1 Tangshan Soil gas CO2 platform array observation ring |
|
图 1 唐山地区土壤气CO2浓度观测台阵布设点 Fig.1 Soil gas CO2 layout point in Tangshan area |
在唐山地震危险区内历史发震断层及多条重点监视断层布设由测点观测和台网接收部分组成的无人值守全自动CO2浓度观测台阵（图 2），台站观测部分由GPRS无线数据通讯系统、CO2观测系统（包括CO2探头、温度探头）、太阳能供电系统等组成，可实现CO2观测、数据管理、数据通讯、系统供电、供电管理等功能；台网接收部分由GPRS无线数据通讯系统、交换机（或路由器）、数据接收计算机、数据处理计算机等组成，实现数据接收、数据库管理、观测台站管理、观测日志自动形成。台阵包括三山院、滦县台、赵各庄矿、马家沟矿、陡河台、孙家楼、滦南靳庄、赵滩、林西矿、北泊村等10个CO2浓度观测台站，由于数据库技术问题和因农耕遭到破坏等原因，CO2浓度观测台阵中有的数据并不完整，从中国地震前兆台网数据处理系统选取2018年8月至2019年10月滦县台、赵各庄矿、马家沟矿、陡河台、滦南靳庄、赵滩、林西矿、北泊村等8个CO2浓度观测台站产出的分钟值数据，分析土壤气CO2浓度数据的年变规律、与地温间的相关性、映震效果等。
|
图 2 无人职守全自动CO2浓度观测台阵构成模式示意图 Fig.2 Schematic diagram of the composition pattern of the unattended and fully automatic CO2 observatory array |
图 3为唐山断裂带土壤气CO2浓度日均值。由图 3可见，CO2浓度存在正常的日变、年变，年变规律较明显，具有短期正常动态特征，即夏高冬低的特点，高低值之间的变化为渐变关系，而不是大幅度的突变，按照形态可分为平稳型和锯齿型2种类型。
|
图 3 唐山土壤气CO2浓度日均值 （a）陡河台；（b）北泊村；（c）乐亭赵滩；（d）滦南靳庄村；（e）滦县地震台；（f）马家沟矿；（g）赵各庄矿；（h）林西矿 Fig.3 Tangshan soil gas CO2 platform array data |
测定土壤气CO2浓度时干扰因素分为人为干扰和自然干扰，人为干扰可以通过严格的操作、规范的方法降至较低水平，而自然干扰中地温因素尤为重要，所以分析地温与土壤气CO2浓度间的相关性，可为土壤气CO2浓度资料的利用和地震前异常信息的提取提供一种可行的方法。
（1）北泊村地温具有明显的年变特征，通过对比土壤气CO2浓度与地温间的关系分析地温对土壤气CO2浓度的影响（表 2、图 4）。由表 2、图 4可见，1—4月，土壤气CO2浓度与地温间相关性不大；5—12月，二者相关性较大，CO2浓度年变化特征清晰，具有夏高冬低的年变特征，年变最大幅度约11 036×10-6/d，二者变化时间较同步，没有明显的时间滞后。经计算北泊村土壤气CO2浓度与地温间总体相关性为0.903。
|
|
表 2 北泊村观测层地温与断层土壤气CO2浓度间的关系 Table 2 Comparison of soil gas CO2 relationship between the North Moor observation layer and fault soil |
|
图 4 北泊村土壤气CO2浓度和地温 Fig.4 Soil gas carbon dioxide and ground temperature in Beibo |
（2）林西矿地温总体变化幅度不大，约11.8—16.0℃，土壤气CO2浓度年变特征明显，总体呈夏高冬低的特点，与地温呈正的相关，相关系数为0.813（图 5）。
|
图 5 林西矿土壤气CO2浓度和地温 Fig.5 Soil gas carbon dioxide and ground temperature in Linxi Mine |
（3）滦南靳庄地温年变特征明显，土壤气CO2浓度随温度升高而升高，随温度降低而下降，总体呈夏高冬低的趋势，但趋势并不明显，与地温间的相关系数为0.769（图 6）。
|
图 6 滦南靳庄村土壤气CO2浓度和地温 Fig.6 Soil gas carbon dioxide and ground temperature in Jinzhuang Village, South Luan |
（4）马家沟矿地温年变特征明显，土壤气CO2浓度随地温变化的特征并不明显，总体上具有夏高冬低的特点，与地温间的相关系数为0.546（图 7）。
|
图 7 马家沟矿土壤气CO2浓度和地温 Fig.7 Soil gas carbon dioxide and ground temperature in Majiagou mine |
（5）乐亭赵摊、赵各庄矿、滦县台、陡河台的土壤气CO2浓度变化曲线为锯齿型，同时，具有夏高冬低的特点，但是与地温间的相关系数并不大，分别是0.315、0.282、0.103、0.139（图 8—11）。
|
图 8 乐亭赵滩土壤气CO2浓度和地温 Fig.8 Soil gas carbon dioxide and ground temperature in Zhaofang, Leting |
|
图 9 赵各庄矿土壤气CO2浓度和地温 Fig.9 Soil gas carbon dioxide and ground temperature in Zhaogezhuang Mine |
|
图 10 滦县地震台土壤气CO2浓度和地温 Fig.10 Soil gas carbon dioxide and ground temperature in Tai, Luanxian County |
|
图 11 陡河台土壤气CO2浓度和地温 Fig.11 Soil gas carbon dioxide and ground temperature in Douhe Platform |
2018年8月至2019年2月，唐山地震老震区（39°—41.5°N，117°—119°E）发生ML 3.0以上地震4次（表 3），按照历史上划分的唐山地震老震区3个区域，即东、中、西区来分析唐山断裂带CO2浓度观测台阵的映震情况。
|
|
表 3 2018年8月至2019年2月唐山地震老震区ML≥3.0地震 Table 3 The ML 3.0 earthquake in the old earthquake area of Tangshan from August 2018 to February 2019 |
2018年8月唐山地震老震区中区发生2次ML 3.0以上地震，分别为河北古冶5日ML 3.3地震、10日ML 3.0地震，2次地震间隔时间仅为5天，且震中都位于唐山—古冶断裂附近。地震前林西矿土壤气CO2浓度日均值持续升高，直至8月14日出现下降趋势；陡河台土壤气CO2浓度日均值呈波动上升趋势，在小幅度上升和下降的转折点发生地震（图 12）。说明林西矿、陡河台土壤气CO2浓度对唐山地震老震区中区尤其是唐山—古冶断裂附近地震映震效果较好。
|
图 12 林西矿、陡河台土壤气CO2浓度日均值 Fig.12 Daily mean of soil gas in Linxi Mine and Douhetai |
2018年11月29日ML 3.2地震发生在唐山丰南，震中位于蓟运河断裂、唐山断裂、宁河断裂交汇处东部，即唐山地震老震区西区。地震前16天马家沟矿土壤气CO2浓度日均值出现小幅上升，上升持续7天、数值趋于平稳6天后发生了ML 3.2地震；滦南靳庄村土壤气CO2浓度日均值在地震前19天出现下降趋势，然后呈锯齿状的上升和下降，在小幅度上升和下降的转折点发生地震，地震发生后马家沟矿、滦南靳庄村土壤气CO2浓度日均值恢复平稳状态（图 13）。其他台站数据无明显的波动，因此马家沟矿、滦南靳庄村台站对唐山地震老震区西区地震映震效果较好。
|
图 13 马家沟矿、滦南靳庄村土壤气CO2浓度日均值 Fig.13 Daily average of soil gas carbon dioxide in Majiagou Mine and Jinzhuang Village, Luannan |
2019年2月14日ML 3.0地震发生在唐山地震老震区东区，震中位于滦县—乐亭断裂附近。在地震发生前10天，乐亭赵滩、滦县台土壤气CO2浓度日均值均出现不同程度的大幅度震荡式上升和下降，在数值上升的最高点发生地震（图 14），其他台阵数据无明显波动。乐亭赵滩、滦县台站对唐山地震老震区东区地震映震效果较好。
|
图 14 乐亭赵滩、滦县地震台土壤气CO2浓度日均值 Fig.14 Daily mean value of soil gas and carbon dioxide in Zhaojiatan and Luanxian Tai |
通过分析2018年8月至2019年2月唐山断裂带土壤气CO2浓度日均值，得到以下结论。
（1）唐山断裂带土壤气CO2浓度年变规律明显，主要受观测地温的影响，并伴随观测地温的变化而变化，即在气温较高的夏季，大都呈升高的趋势，且有明显的年变规律。
（2）通过计算土壤气CO2浓度与地温间的相关系数，发现北泊村、林西矿、滦南靳庄、马家沟矿土壤气CO2浓度与地温间的相关系数较大，土壤气CO2浓度曲线属于平稳型；乐亭赵摊、赵各庄矿、滦县台、陡河台的土壤气CO2浓度与地温间的相关系数不大，土壤气CO2浓度曲线属于锯齿型。
（3）在活动断裂及周边地震的震中区，土壤气CO2浓度日均值存在异常变化，变化的幅度与台站位置、对应区域的构造活动有关，与土壤气CO2浓度和地温间的相关系数大小、CO2浓度曲线是否为平稳型或锯齿型关系不密切。总体上，土壤气CO2浓度对唐山地区ML 3.0以上地震有一定的映震效果。
杜建国, 王先彬, 谢鸿森. 深部物质运动的气体地球化学特征[J]. 地球科学进展, 1994, 9(3): 48-52. |
耿杰, 薛革, 魏焕, 等. 环渤海地区断层气CO2监测网观测结果的初步分析[J]. 华南地震, 2000, 20(1): 28-32. |
蒋雨函, 高小其, 杨朋涛, 等. 新疆北天山地区断裂带断层土壤气的地球化学特征[J]. 地震地质, 2022, 44(6): 1597-1614. |
李营, 杜建国, 王富宽, 等. 延怀盆地土壤气体地球化学特征[J]. 地震学报, 2009, 31(1): 82-91. |
林元武, 王基华, 髙松升. 断层气CO2测定新方法与张北—尚义6.2级地震预报[J]. 地震, 1998, 18(4): 353-357. |
刘燕翔, 王妍, 李金, 等. 2019年唐山ML 4.9地震与1976年唐山7.8级地震的相关性分析[J]. 地震地磁观测与研究, 2021, 42(1): 42-47. |
盛艳蕊, 张子广, 周晓成, 等. 新保安—沙城断裂带土壤气地球化学特征分析[J]. 地震, 2015, 35(4): 90-98. |
盛艳蕊, 张子广, 丁志华, 等. 唐山断裂带土壤气地球化学特征分析[J]. 震灾防御技术, 2020, 15(2): 452-462. |
王江, 李营, 陈志. 口泉断裂断层气地球化学变化特征及断层活动性[J]. 地震, 2017, 37(1): 39-51. |
王江, 陈志, 张帆, 等. 基于土壤气体地球化学的雄安新区活动断裂空间展布及活动性探讨[J]. 地震研究, 2022, 45(2): 264-274. |
王江, 陈志, 张帆, 等. 雄安新区主要断裂带土壤气体的Rn与CO2脱气特征[J]. 地震地质, 2023, 45(3): 735-752. |
王喜龙, 李营, 杜建国, 等. 首都圈地区土壤气Rn, Hg, CO2地球化学特征及其成因[J]. 地震学报, 2017, 39(1): 85-101. |
杨江, 李营, 陈志, 等. 唐山断裂带南西段和北东段土壤气Rn和CO2浓度特征[J]. 地震, 2019, 39(3): 61-70. |
张素欣, 王晓山, 陈婷, 等. 唐山老震区40年地震时空演化特征分析[J]. 华北地震科学, 2017, 35(1): 32-37. |
张肇诚. 中国震例(1976—1980)[M]. 北京: 地震出版社, 1990: 60.
|
郑海刚, 方震, 周晓成, 等. 郯庐断裂带安徽段土壤气体的地球化学特征[J]. 中国地震, 2016, 32(4): 642-652. |
仲秋, 史保平. 1976年唐山MS 7.8地震余震序列持续时间及对地震危险性分析的意义[J]. 地震学报, 2012, 34(4): 494-508. |
朱成英, 闫玮, 麻荣, 等. 2017年8月9日精河MS 6.6地震宏观烈度及其余震分布的断层气体地球化学表征[J]. 地震地质, 2022, 44(5): 1225-1239. |