引言

近年来，针对地震活动、断裂构造进行的地球化学探测逐渐兴起并发展（郑海刚等，2016；王喜龙等，2017），氢气作为断裂带主要气体成分之一，被认为与地震构造活动密切相关（杜乐天，2005；高小其等，2017；向阳等，2018）。众所周知，自然界氢气来源主要有以下几点：（1）断层错动：沿地壳下层塑性岩石（蛇纹岩）或硅酸岩，断层新鲜裂隙、断层错动过程中产生高浓度氢气（陈辉，1996；夏林圻等，2008）；（2）应力活动：主要来源为地壳深部岩石孔隙、裂隙中封存的氢气（孙小龙等，2017）；（3）生物化学作用也可产生氢气，但仅有微量，且存在于浅层地表（夏林圻等，2008；车用太等，2002）。地震发生前氢气浓度异常已被众多震例证实，且异常多出现在震前短临阶段（林元武等，1994；张彬等，2018）。国内外已有众多震例证明氢气在震前、震后均发生明显变化，地下氢气是灵敏反映地下应力应变的地球化学成分（高清武，1992；张培仁等，1993；车用太等，2002），断裂带氢气浓度与地震活动释放的能量呈正相关（Ryuichi，1984）。近年来，很多学者利用氢气研究地震、断层，取得了良好的效果（李源等，2018；康健等，2019）。氢气浓度观测是值得广泛推广的地震监测方法之一，也被认为是具有广阔前景的地震监测预报方法（车用太等，2002；刘耀炜等，2006；范雪芳等，2016）。

依兰-伊通断裂是东北地区规模最大的断裂，近年来在断裂通河段发现了1次古地震事件，探槽揭示距今1700年前发生过7级以上地震事件（闵伟等，2011），这一发现彻底颠覆了人们对该断裂“全新世以来不活动”的认识。本文利用地球化学探测方法研究依兰-伊通断裂（图 1），了解依兰-伊通断裂断层氢气时空变化特征，捕捉短期地震构造活动变化信息，建立氢气背景值及氢气与地震活动的耦合关系。

图 1  依兰-伊通断裂北段构造图

Figure 1.  Structural map of the northern section of Yilan-Yitong fault

下载: 全尺寸图片 幻灯片

1.   构造背景

依兰-伊通断裂总体呈NNE-NE向展布，中国境内总长约1100km，是郯庐断裂带渤海以北区段的主要组成部分（徐嘉炜等，1992），也是划分中蒙、中朝和燕山3个活动地块的重要边界断裂（张培震等，2003）。该断裂基本由近于平行的东、西2条分支断裂（东支F4-2和西支F4-1）组成，构成宽10—20km的地堑，部分区段呈北东向沿松花江分布。

依兰-伊通断裂带在黑龙江境内一般分为北段（方正-萝北段）、中段（尚志-延寿-方正段）和南段（伊通-舒兰-尚志段）。依兰-伊通断裂北段地震活动明显强于中段和南段，沿着北段存在整个黑龙江省内最密集的小震活动，显示了较强烈的当前地壳活动特征。就北段而言，地震活动分布不均，萝北地区地震活动强于北段其他区域，该地区小震活跃，历史上发生过多次中强地震，如1963年6月22日M5.8地震。研究发现西支断裂较东支断裂活动性明显，西支通河段发现存在高2m左右的断层陡坎等地貌线性影像，延伸长约70km，逆走滑性质，并在祥顺乡南楼村探槽剖面发现全新世活动，并推测约在1700年前发生过强烈古地震事件（闵伟等，2011；疏鹏等，2014）。

2.   测线布设与测量方法

选择在依兰-伊通断裂北段西支断裂上进行研究，由方正至萝北选出8条剖面垂直断裂进行观测（图 1），由南至北依次为方正西测线、祥顺乡测线、松花江农场测线、香兰农场测线、汤原测线、鹤岗南测线、宝泉岭测线和延军农场4队测线。观测剖面垂直断裂，每条测线等距布设调整间距15—20调整间距个测点，点距一般为调整间距10—50m，具体间距以实际情况为准，可根据地形地貌和断层类型调整。如果未发现断层具体位置，则根据收集到的资料分析断层可能的位置，用于确定最佳观测位置。

观测仪器为便携式氢气分析仪（ATG-300H型），标定后空气检测值一般为0.5×10^-6。测量时需在土壤中打观测孔，测孔直径约3cm，孔深约80cm。首先将取样器（麻花钻）吸气端插入孔中并封住气孔，然后用软管将取样器排气端与仪器相连，并启动仪器开始观测，读取最高值。为减少非构造因素的影响，选取土质条件相同的场地进行观测。

3.   观测结果

3.1   数据

经多期观测，得到断裂8条测线氢气浓度背景值。背景值采用各测点浓度值的几何平均值，这种方法能弱化极端数值对整体数据的影响，从而更真实地反映数据的整体水平，一般可体现研究区域客观氢气浓度背景。计算公式如下：

式中，G为观测点所有样本几何平均值；${x_n}$为观测点第n个样本；n为观测点样本个数。

观测时间段为2017年6月至2018年9月，每条测线观测2—5期。其中，方正西测线、祥顺乡测线、汤原测线完成2期测量；香兰农场测线、宝泉岭测线、延军农场4队测线完成3期测量；鹤岗南测线完成4期测量；松花江农场测线完成5期测量。整体数据如表 1所示。

表 1  测线断层氢气统计结果

Table 1.  Hydrogen statistical results of line faults

测线序号	测线名称	观测日期/年-月-日	测点个数	背景值/ppm
1	方正西	2017-09-12	15	62
		2018-08-19	15	51
2	祥顺乡	2017-09-12	15	26
		2018-08-18	15	260
3	松花江农场	2017-09-13	15	267
		2017-11-02	15	106
		2018-08-18	15	602
		2018-09-13	15	245
		2018-09-19	15	162
4	香兰农场	2018-01-19	30	47
		2018-04-02	30	178
		2018-08-17	15	359
5	汤原	2017-07-18	15	98
		2018-08-15	15	168
6	鹤岗南	2017-07-19	15	744
		2017-11-03	15	33
		2018-08-12	15	214
		2018-09-14	15	158
7	宝泉岭	2017-09-14	15	360
		2018-08-13	15	569
		2018-09-16	15	280
8	延军农场4队	2017-06-12	15	92
		2017-07-20	15	168
		2018-08-13	15	434

下载: 导出CSV 

| 显示表格

3.2   数据分析

将各条测线氢气浓度背景值按月统计，得到2017—2018年部分时段的背景值，如图 2所示。由图 2可知，香兰测线和延军农场4队测线氢气背景值由春季到夏季逐渐升高；松花江农场测线、鹤岗南测线和宝泉岭农场测线氢气背景值由夏季到秋季逐渐降低，9月氢气背景值较8月大幅下降。

图 2  测线背景值分布图

Figure 2.  Distribution map of background values of lines

下载: 全尺寸图片 幻灯片

由图 2可知，春季至夏季为氢气浓度上升段，夏季至秋季为氢气浓度下降段，每年各条测线最高值基本出现在8月前后。杜建国等（1998）发现土壤氡气浓度值存在冬夏高、春秋低的年变规律，并认为断层气在土壤中表现出的夏季高主要受土壤湿度的影响（王博等，2010；孙小龙等，2017），冬季高主要受冻土层的影响（康健等，2019）。土壤中的氡气与氢气源于地壳深部，运移路径与方式相似。图 2中数据虽不完整，但仍可显示氢气浓度具有夏季高、春秋低的特点。

将2018年观测数据与2017年同期数据进行比较，由于条件限制，使每年的观测时间难以保持一致，相对而言，2017、2018年7—9月（夏季）观测最集中，因此以7—9月为基础，最终得到夏季对比数据，如图 3所示。祥顺乡测线、鹤岗南测线、延军农场4队测线氢气浓度同期变化较显著，这可能与构造活动有关。其余5条测线氢气浓度同期变化不大，相关系数达0.87，可见背景值变化较稳定，如松花江农场测线2017年9月13日测得的背景值为267ppm，于2018年9月13日再次观测，测得背景值为245ppm，表明该地区2017、2018年氢气浓度非常稳定。同期变化较大的测线中，鹤岗南测线背景值2018年较2017年明显降低，而祥顺乡测线、延军农场4队测线背景值2018年较2017年明显升高。将依兰-伊通断裂北段观测到的每条剖面进行整体比较（图 3），可知鹤岗南测线2017年背景值最高，为742ppm，祥顺乡测线最低，为26ppm，整条断裂年平均值为279ppm；延军农场4队测线2018年背景值最高，为434ppm，汤原测线最低，为168ppm，整条断裂年平均值为268ppm。综上所述，2017、2018年平均值近乎持平，但从变化幅度上看，2018年较2017年更趋于稳定状态。

图 3  依兰-伊通断裂北段各条测线断层氢气统计结果

Figure 3.  Hydrogen statistical results of line faults of the northern section of Yilan-Yitong fault

下载: 全尺寸图片 幻灯片

2018年8月在依兰-伊通断裂北段观测了完整的一期，所有测线均在同一时段内完成观测。图 4中红色数据为2018年8月各条测线氢气浓度分布特征，可知在本文研究的断裂区域内，松花江农场测线、香兰农场测线、宝泉岭测线和延军农场4队测线氢气浓度较高；蓝色数据为2017、2018年氢气浓度平均背景值，可知其变化趋势与氢气浓度变化趋势相似，说明2017、2018年氢气浓度变化较稳定，也可知夏季浓度值普遍高于年平均值。

图 4  2017—2018年依兰-伊通断裂北段氢气浓度特征划分图

Figure 4.  Characteristic partition of hydrogen concentration of the northern section of Yilan-Yitong fault from 2017 to 2018

下载: 全尺寸图片 幻灯片

同一时刻的横向比较更能看出断裂氢气浓度分布特征，因此将整条断裂数据进行横向比较，根据2018年8月观测结果，按浓度值将断裂简单划分为方正-依兰段（平均氢气浓度＜250ppm）、依兰-汤原段（平均氢气浓度＞450ppm）、汤原-鹤岗段（平均氢气浓度＜250ppm）、鹤岗-萝北段（平均氢气浓度＞450ppm），如图 4所示，可见断层氢气在断裂不同地段释放量不同。

4.   机理分析

对断裂上测线进行环比分析，可知氢气浓度呈夏季高、春秋低的特点。土壤氢气浓度具有年变特征，造成这种情况的原因是温度、湿度变化（刘学领等，2011；康健等，2019）。研究表明微生物活动可产生氢气，而温度、湿度影响微生物活跃程度，温度低于10℃时大量微生物停止活动或死亡，温度为20—30℃时微生物最活跃（樊耀亭等，2002；Des Marais，2000；Kasting等，2006；梅冥相，2014；车用太等，2015）。湿度增加会减小土壤孔隙度，孔隙度变小会造成土壤中气体浓度增加（孙小龙等，2017）。因此，夏季温度、湿度最适合微生物生长，而黑龙江地区春秋季温度较低，且干燥少雨，微生物基本停止活动。因此土壤氢气浓度夏季高、春秋低的年变特征主要与微生物活动有关。

对断裂上测线进行同期分析，可知测线2017、2018年平均值近乎持平。特别是松花江农场测线2017、2018年观测日期一致，观测数值也近乎相等，说明该地区断层氢气在2017、2018年输出稳定。个别地区前后2年同期相差较大，分析原因如下：（1）2017、2018年降雨量不同，土壤湿度不同，造成土壤内孔隙比不同，影响氢气浓度变化，且湿度不同还影响微生物的活跃程度；（2）地下地震构造活动存在变化。

对断裂同一时段进行横向比较，可知每个区域氢气浓度变化很大，按浓度特征可简单划分为4个区域。将分区结果结合区内地震活动特征进行分析（图 5），可知地震活动有明显的分段性。图 5（a）显示氢气浓度高低与历史地震活动强弱的对应关系并不明显，图 5（b）显示氢气浓度高低与当前地震活动存在对应关系，说明氢气浓度高低与当前地震活动存在耦合关系（车用太等，2002；李晨桦等，2016）。特别是依兰-汤原段香兰农场测线附近，2018年1月连续发生3.2级和3.8级地震，与地震对应的分区氢气浓度偏高，说明断层气浓度反映当前地震构造活跃情况，与当前地震活动有直接关系，与历史地震甚至地震多发区无关。

图 5  依兰-伊通断裂北段地震分布图

Figure 5.  Earthquake distribution map of the northern section of Yilan-Yitong fault

下载: 全尺寸图片 幻灯片

5.   结论与讨论

通过对依兰-伊通断裂北段氢气流动进行观测，得出以下结论：

（1）从时间纵向比较，依兰-伊通断裂北段2018年背景值与2017年较接近，受地震构造活动等因素影响局部升高，个别有下降；氢气浓度存在年变特征，主要受温度、湿度和微生物活动的影响，一年之中夏季氢气浓度最高，春秋季氢气浓度最低。

（2）从空间横向比较，依兰-伊通断裂北段呈现出分段性特征，氢气浓度在该段断裂上呈现出北高南低的特征，这与当前地震活动特征相符，高值地区对应地震活动较强的地区，低值地区对应地震活动较弱的地区，断层气浓度与当前地震构造活动有很好的耦合关系。

（3）本文没有连续观测数据和显著震例，无法捕捉断层氢气浓度在震前、震后的变化特征，但不难看出构造地球化学方法对当前地震构造活动映震效果明显。如果增加观测手段，多种断层气同步观测，相互印证，将极大提高观测结果的可信度。对重点地区进行连续观测可提高观测精度并及时掌握断层气动态变化，有利于监视地震构造活动。相信随着科技的发展，构造地球化学观测将在分析预报、短临跟踪、断层研究等领域发挥重要作用。