嘉峪关气氡浓度异常与202...多M_S7.4地震关系分析_曹玲玲.pdf

第 卷第期 年月地震工程学报 ，收稿日期：基金项目：甘肃省地震局地震科技野外站基金（，）；地震科技星火计划项目（）第一作者简介：曹玲玲（），女，甘肃会宁人，副研究员，主要从事地下流体预测地震研究。：。通信作者：苏鹤军（），副研究员，主要从事构造地球化学及同位素地球化学研究。：。曹玲玲，苏鹤军 嘉峪关气氡浓度异常与 年玛多地震关系分析 地震工程学报，（）：，（）：嘉峪关气氡浓度异常与 年玛多 地震关系分析曹玲玲，苏鹤军，（甘肃兰州地球物理国家野外科学观测研究站，甘肃 兰州 ；甘肃省地震局，甘肃 兰州 ；中国地震局兰州岩土地震研究所，甘肃 兰州 ）摘要：嘉峪关气氡浓度在多年上升的背景下，于 年出现了转折异常变化，但祁连山地震带内一直没有发生与异常幅度和持续时间相匹配的地震，直到 年月 日，在距离嘉峪关气氡测点 的玛多发生地震。为了判断嘉峪关气氡浓度异常与玛多地震的关系，结合震例，从异常信度、震前异常特征、震后异常变化及地质构造背景等方面对其进行深入分析。结果发现：嘉峪关气氡浓度的异常特征与地震所处的构造有关，发生在祁连山地震带内的地震，气氡浓度异常表现为年或半年尺度的年畸变，发生在该地震带以外远距离的地震，气氡浓度异常表现为多年的趋势性变化；玛多地震发生后，嘉峪关气氡浓度下降速率减缓之后转折，呈恢复状态；嘉峪关气氡测点与玛多地震都位于青藏高原东北部，具有相同的动力背景且在构造上具有关联性。综合分析认为嘉峪关气氡多年趋势异常与玛多地震有关，研究对建立可靠的异常指标体系、提高地震预测水平具有重要意义。关键词：嘉峪关气氡浓度；趋势异常；玛多地震；相关性中图分类号：文献标志码：文章编号：（）：，（.，；.，；.，）：，：；引言地震的有效预测是减轻地震灾害的前提。目前预测地震的主要依据是以震例为基础建立的预测指标体系，但由于地震前兆异常的低重现性和震例统计的概率性特性，当测项出现异常时，依据现有的预测指标还无法对未来地震的三要素做出更为准确的预报。因此，当一次地震，特别是发生概率更低的大震发生后，进行前兆异常与地震关系的深入分析，对建立更加可靠的地震预测指标体系、提高地震预测水平具有重要意义。年月开始，甘肃省嘉峪关气氡浓度出现了加速上升变化，年月达到峰值后转折下降，截至 年月，测值未恢复到上升前的水平（图）。年、年和 年，工作人员对此变化 都 进行 过 现场 核实，判定结果为：气 氡浓度 年月 年月的上升存在周围环境改变的影响，年月以来的趋势下降为地震前兆异常。如此大幅度且持续时间长达几年的异常在该测点是首次出现，由于没有震例，无法对地震三要素做出较为准确的预测。年月 日，距离嘉峪 关 气 氡 测 点 的 青 海 玛 多 县 发 生 了地震。根据 中国震例，级以上地震地球物理异常的分布范围超过 ，据此判断嘉峪关气氡浓度的这次异常可以作为玛多地震的前兆异图嘉峪关气氡浓度月均值曲线图 地震工程学报 年常。但地球物理场异常与地震之间的关系不只受距离约束，还与构造、异常形态和区域背景等有关，仅根据震中距来判断异常与地震的关系，依据并不充分，也不可靠。基于此，本文从异常的可靠性、不同区域地震前异常的差异性、震后异常变化，以及构造等方面，对嘉峪关气氡浓度 年以来的异常与玛多地震之间的关系进行深入分析，给出判断两者关系的可靠依据，以期为异常指标体系的完善和地震的有效预测提供依据。测点概况嘉峪关气氡测点位于嘉峪关断层上。该断裂北起榆树沟山，向南延伸至文殊沟口，全长约 ，总体呈 方向展布，倾向，倾角 ，晚第四纪以来新的构造活动以挤压逆冲为主兼右旋走滑。断层切割奥陶系、白垩系、新近系和上更新统，在剖面上表现为盘上升、盘下降。全新世有古地震活动，第四纪晚期以来的平均滑动速率为 （图）。图嘉峪关断层地质构造图 为了监测嘉峪关断层的活动，年在该断层某处（，）建立了监测点，进行断层土壤气氡浓度的观测。如图所示，在监测点开挖一深、直径 的垂直测孔，测孔下部 为集气坑，上部 为一倒置的长颈玻璃漏。漏斗颈内布设个直径为 的导气管：个导气管延伸至孔底，进行气氡浓度主样观测；个下至漏斗颈口，进行气氡浓度比对观测。导气管通过转接铜管与乳胶管相接并从测孔引出，乳胶管口与大气连通。气氡浓度观测至今，都为人工测量，每天测量一次。图嘉峪关气氡测点实验装置示意图 气氡浓度主要影响因素分析玛多地震前，嘉峪关气氡浓度未出现短期异常，因此本文只分析年尺度和长趋势的干扰因素。已有研究表明，气氡浓度受到气温和气压的影响，无震年份会出现年变变化。嘉峪关气氡浓度在正常背景下具有明显的夏高冬低的年变形态（图），由此可以判断，气温和气压为气氡浓度的固定干扰因素。但是气温和气压之间具有相关性，究竟哪项对气氡浓度的影响更为显著？由于相关系数的大小可以表征变量之间的相关程度，因此可以通过计算气氡浓度与气温和气压的相关系数来判断。由于中强以上地震发生的年份，气氡浓度变化受构造应力的影响会偏离正常动态，因此 年以来祁连山地震带以上地震所发生的年份不参与计算（图、表）。具体步骤：先计算出气氡浓度、气氡测点气温及气压的月均值，再利用式（）分别计算出气氡浓度与气压、气温的相关系数：（?）（?）（?）（?）（）式中：为气氡浓度月均值，?为气氡浓度自然年月均值的平均值；为气温（气压）月均值，?为气温（气压）自然年月均值的平均值。气氡浓度与气压、气温的相关系数计算结果如第 卷 第期曹玲玲，等：嘉峪关气氡浓度异常与 年玛多地震关系分析图 年以来祁连山地震带以上地震空间分布 表 年以来祁连山地震带以上地震基本参数 震中震级（）发震日期北纬（）东经（）震中距 甘肃景泰 甘肃肃南 甘肃玉门 青海德令哈 甘肃民乐 青海德令哈 甘肃肃南 青海海西 青海海西 甘肃金塔 青海门源 青海门源 甘肃肃南 图所示。从图中可以清晰地看出，气氡浓度与气温呈正相关，与气压呈负相关；同期气氡浓度与气温的相关系数都大于其与气压的相关系数，并且前者的值都在以上，而后者仅有个值超过。据此可以判断，在年尺度上，气温是气氡浓度变化的主要影响因素。图 嘉峪关气氡浓度与气温和气压的相关系数曲线 多次现场核实结果表明，周围环境的变化也会使气氡浓度值出现趋势性的异常变化，对此文献已进行过详细分析，这里不再赘述。气氡浓度异常特征分析破年变异常与地震的对应关系畸变型破年变异常与地震对应关系相关系数计算结果（图）显示，嘉峪关气氡浓度主要受气温的影响，在祁连山地震带及其周边无 以上地震的年份，气温与气氡浓度具有很好的正相关关系，因此气氡浓度的正常年变形态为夏高冬低，峰值出现在每年的月，谷值出现在每年 月次年月。当气温与气氡浓度的相关系数降低或变为负值时，气氡浓度年变会出现畸变，因此可以利用相关系数判断气氡浓度是否出现了畸变型破年变异常。由于测点缺少 年的气象资料，考虑地震工程学报 年到测点与地震震中的构造关系，选取 年以来祁连山地震带发生的 以上地震为震例，震中距不做限制，分别分析震前年和震前个月气氡浓度与气温相关系数的变化，依此确定破年变异常的阈值。所选地震的空间分布见图，具体参数列于表。为了定量给出年变畸变异常的阈值，采用以下判定方法：个月尺度的异常，采用个月窗长，个月步长；年尺度的异常，采用 个月窗长，个月步长；再利用式（），向后滑动计算出气氡浓度与气温月均值的相关系数，以报准率（漏报率虚报率）比值的最大值作为异常阈值，结果为。因此本文将 作为畸变型破年变异常的判定阈值，小于此值则为异常。图为所选地震前年和前个月气氡浓度与气温的相关系数柱状图，横坐标为具体地震及其震中距。从图中可以看出，震前年，小于的所有地震前，气氡浓度与气温的相关系数均未达到异常阈值，气氡浓度未出现异常；次以上地震中，有次地震前气氡浓度与气温的相关系数达到异 常 阈值，气 氡 浓 度 出现了异常，对应 率高达，只有震中距超过 的景泰地震前气氡浓度未出现异常。震前个月，小于的所有地震中，只有震中距最大的门源地震前气氡浓度和气温的相关系数达到异常阈值，气氡浓度出现了异常；大于的次地震中，有次地震前两者的相关系数达到异常阈值，气氡浓度出现了异常，未出现异常的两次地震为距离最远的景泰地震和次远的门源地震。图祁连山地震带以上地震前嘉峪关气氡浓度与气温相关系数 综合以上分析，对于祁连山地震带以下地震，无论震中距大小，嘉峪关气氡浓度通常不会出现年尺度或个月尺度的畸变型破年变异常；对于以上地震，气氡浓度是否出现此类异常，与震中距大小有关。已有的震例分析结果显示，对应的震中距范围与 中国震例 的规定相当。年变幅增大的破年变异常与地震的对应关系破年变异常不仅包含年变畸变异常，还包括年变幅增大的异常。此类异常可以利用年标准曲线法提取，即首先计算出表征台站固有属性的年标准值，以倍均方差线作为异常的阈值，超过阈值则达到异常。具体计算方法参见文献。为了与畸变型破年变异常的分析结果进行对比，分析区间依然选择 年，计算结果见图，图中粗线条为观测值曲线，细线条为年标准曲线的一倍方差线，也就是阈值线。从图可以看出，年以来，有个年份嘉峪关气氡浓度的年变幅增大并超过阈值线，且都在 年之后。本文研究的 次地震中，所有以下地震前都未出现 此 类 异 常，以 上 地 震 中 只 有 门 源地震前出现了异常。据此可以判断，嘉峪关气氡浓度单个年变幅增大的异常不能作为地震预测指标。长趋势转折异常与地震的对应关系图显示嘉峪关气氡观测以来，持续多年的长趋势上升转折异常共出现次，第次为 年，第次为 年。第次异常结束后发生了景泰地震和昆仑山西口地震。第次异常期间，除玛多地震外，祁连山地震带及周边没有发生震级与异常相匹配的地震。第 卷 第期曹玲玲，等：嘉峪关气氡浓度异常与 年玛多地震关系分析图嘉峪关气氡浓度测值与阈值曲线图 讨论与分析气氡浓度趋势异常可靠性分析 年嘉峪关气氡浓度开始出现大幅上升，至今进行过次现场核实，其中次在异常上升期间，次在异常下降期间。核实结果为：上升期间存在环境变化的影响，影响程度无法量化；转折下降过程没有干扰因素，也不存在测点更换和观测仪器故障等情况。上述相关性分析结果显示，气温对气氡浓度的影响显著，二者呈正相关。从图可以看出，年至今气温的变化并不大，仅 年的峰值略低，这与 年开始的气氡浓度加速上升不一致；年之后，气温变化非常平稳，但气氡浓度出现了明显的转折下降。综合现场核实结果判断，年以来气氡浓度的转折是区域应力变化的响应。图嘉峪关气氡浓度与气温对比曲线 嘉峪关气氡趋势转折异常特征与玛多地震关系分析不同区域地震前气氡浓度异常特征差异性分析上述分析结果显示，嘉峪关气氡浓度异常主要有三类：第一类为畸变型破年变异常；第二类为年变完整但年变幅显著增大的异常；第三类为持续多年的趋势异常。其中，第二类异常不能作为地震预测的依据。第一类异常是祁连山地震带地震前的典型异常，对应震级为以上。从已有震例来看，对应以下地震的最大震中距不超过 ，异常开始距地震的发生时间通常在一年以内，主要表现为气温和气氡浓度相关系数的减小，而嘉峪关气氡浓度 年以来未出现畸变型破年变异常，因此其目前的异常与祁连山地震带地震前的异常不一致。第三类异常观测以来共出现次（图）。第次异常：嘉峪关气氡浓度连续上升年后转折。转折后 个月，即 年月日在距离测点 处发生了景泰地震；转折后的 个月，即 年 月 日在距离测点 处发生了昆仑山口西地震。前面的分析表明，祁连山地震带地震前嘉峪关气氡浓度为畸变型破年变异常，且震中距基本符合 中国震例 给定的异常统计地震工程学报 年范围，而景泰地震震中位于祁连山地震带中东段，且震中距达 ，因此这次趋势变化作为昆仑山口西地震对应的异常更为合理。第次异常，也就是目前存在的异常。在趋势上升阶段存在环境干扰的影响，由于无法量化，无法判断除环境干扰外是否还包含区域构造应力变化引起的异常。如果包含，这次异常过程与昆仑山口西地震前的异常过程非常相似，因此从异常形态来看，两次地震前的变化高度相似（图）。异常持续时间上，昆仑山西口地震对应异常开始距离发震共 个月，这次异常开始距离玛多地震共 个月，符合震级越大异常开始距离地震发生时间越长的特征。构造块体上，两次地震震中与嘉峪关气氡测点都不在一个二级活动地块，且相距较远。发震断层上，玛多地震发生在东昆仑断裂带以南约 的江错断裂，该断裂向西延伸，与昆仑山口西地震的发震断层昆仑山口断裂相接，因此两次地震的发震断层是相关的，且都为走滑型。由于大震发生的概率低，在该区域并没有更多的震例可供分析，但现有的震例分析结果显示，嘉峪关气氡在相关断层的大震前出现了相近的异常，即异常的重复性特征，这表明这次异常与玛多地震有关。玛多地震后嘉峪关气氡浓度的异常变化震后异常是否改变是判断异常与地震关系的主要依据，根据图无法直观判断玛多地震后嘉峪关气氡异常是否改变。为此计算了差值，即将相邻两年气氡的同期月均值相减（如将 年月的测值减去 年月的测值），根据差值的变化判断气氡在玛多地震前后是否发生了变化。由于气氡转折后相对平稳的下降出现在 年后，因此起始时间选为 年，用前一年的测值减去后一年的测值，结 果见图。可以看出，年的差值在玛图 嘉峪关气氡浓度月均值相邻年份同期差值曲线 多地震前下降幅度较大，地震后下降幅度逐步减小，且 年 月份的测值已经高于 年同期的测值；其他两条曲线并没有这种特征，这表明震后气氡异常趋势发生了改变。据此可以判断嘉峪关气氡的趋势异常与玛多地震有关。大震前不同震中距流体异常特征车用太等发现强震前有三种不同的流体异常，即来自震源的源兆、与地震有区域构造关系的场兆，以及远距离的远兆，华北强震前会出现几百甚至近千公里的流体远兆异常。王世芹等 对云南地区地下水动态的地震源兆与场兆特征分析后发现级地震最远异常井的平均距离为 ，且场兆异常以中期异常居多。高小其等 对新疆级以上地震前地下流体源兆、场兆和远震特征分析后发现强震场兆和远兆都以中期异常为主，异常形态上，场兆中期异常一般是以趋势性大幅升高变化为主，且以水化学测项居多。由此可见，流体远距离的异常在大陆强震前普遍存在，且通常没有阶段性特征，而玛多地震震级达，距离嘉峪关气氡测点 ，异常阶段性上属于中期异常，与其他学者的研究结果一致。因此，嘉峪关气氡浓度的趋势异常作为玛多地震的前兆异常在异常形态、阶段性及震中距等方面是有震例依据的。嘉峪关气氡测点与玛多地震构造关系分析印度板块向北俯冲挤压欧亚板块的长期应力积累和中下地壳低速物质的共同作用是玛多地震的动力来源。嘉峪关气氡测点位于祁连山西段的嘉峪关文殊山断裂，是青藏高原东北缘向北扩展的前缘，也是对高原变形响应最敏感的地形之一。该断裂也分配了印度板块对欧亚板块的北向推挤作用力，因此它们具有相同的动力来源。块体相互关系上，玛多地震位于巴颜喀拉地块，嘉峪关气氡测点位于柴达木祁连地块北边界。巴颜喀拉块体北边界为东昆仑断裂带，西部边界为阿尔金断裂带西南段尾部张性断裂 ；柴达木祁连地块为阿尔金断裂、海源祁连山断裂和东昆仑断裂围陷的一个相对独立的活动地壳块体。因此，嘉峪关气氡测点和玛多地震虽然位于不同的二级地块，但具有共同的块体边界。不同地块之间的构造变形在大区域框架下具有协调性，在这种区域性构造应力作用下，地块相互作用并发生差异运动，继而出现变形和断裂失稳，引发地震，同时引起不同块体的氡浓度出现异常变化。构造上，玛多地震的发震断层向西延第 卷 第期曹玲玲，等：嘉峪关气氡浓度异常与 年玛多地震关系分析伸，与 年地震的发震断层昆仑山口断裂相接。阿尔金断裂的西南段斜切昆仑山，北东段斜切祁连山西段构造带 ，包括嘉峪关断裂在内的祁连山西段的逆冲和走滑断裂吸收了阿尔金断裂的左旋位移。嘉峪关断裂右行走滑活动和嘉峪关隆起形成的区域构造应力场，是整个青藏块体向北推挤和阿尔金断裂东段扩展及其引发的局部应力造成的。因此在构造上，发震断裂应力的积累也会引起气氡测点所在断裂的应力变化，继而使气氡浓度出现异常变化。结论本文从嘉峪关气氡浓度的主要干扰因素、不同异常形态对应地震的区域性特征、震后异常变化，以及构造背景等方面对嘉峪关气氡浓度异常与玛多地震之间的关系进行分析，得出以下主要结论：（）嘉峪关气氡浓度映震震级通常在以上，祁连山构造带内以上地震前的异常为年或个月左右的畸变型破年变异常，异常阈值为气氡浓度和气温的相关系数，因此嘉峪关气氡目前的异常不符合祁连山地震带震前异常的变化特征。（）对比分析了与玛多地震发震构造相关的昆仑山口西地震前嘉峪关气氡浓度的变化，发现两者不仅在异常形态上非常相似，持续时间上也符合震级越大异常持续时间越长的特征。同时分析了玛多地震后嘉峪关气氡浓度的变化，结果显示地震后气氡浓度的下降速率明显减小，表明震后异常发生了改变。（）玛多地震发震断层和嘉峪关气氡测点所在断层不论在区域动力、板块构造，还是在构造关系上，都存在关联。综合分析认为，可以将嘉峪关气氡浓度 年以来的趋势异常作为玛多地震的前兆异常，但还需要在以后的研究中，随着有效震例的增加，进行进一步的验证。参考文献（）苏鹤军，王宗礼，曹玲玲，等 断裂带土壤气测量方法在断层活动性研究中的应用：以嘉峪关断层为例 中国地质，（）：，：，（）：蒋海坤，付虹，杨马陵 中国震例 北京：地震出版社，：，何文贵，袁道阳，王爱国，等嘉峪关断层中段的新活动特征 中国地震，（）：，（）：杜建国，宇文欣，李圣强，等八宝山断裂带逸出氡的地球化学特征及其映震效能 地震，（）：，（）：王博，黄辅琼，简春林 嘉峪关断层带土壤气氡的影响因素及映震效能分析中国地震，（）：，（）：国家地震局地质研究所，国家地震局兰州地震研究所 祁连山河西走廊活动断裂系北京：地震出版社，：，王岩鸿 标准年动态曲线形态分析方法在地下流体数据及前兆异常信息处理中的应用研究兰州：中国地震局兰州地震研究所，：，潘家伟，白明坤，李超，等 年月 日青海玛多地震地表破裂带及发震构造地质学报，（）：，（），（）：车用太，鱼金子地下流体的源兆、场兆、远兆及其在地震预报中的意义地震，（）：，（）：王世芹，刘丽芳，付虹 云南地区地下水动态的地震源兆与场兆特征分析 地震研究，（）：，地震工程学报 年（）：高小其，陈华静，高国英，等新疆地下流体源兆、场兆、大震远兆特征的统计与分析地震，（）：，（）：徐志 国，梁 姗 姗，张 广 伟，等 年月 日 青 海 玛 多地震发震构造分析 地球物理学报，（）：，（）：郑文俊河西走廊及其邻区活动构造图像及构造变形模式 北京：中国地震局地质研究所，：，邓起东，程绍平，马冀，等 青藏高原地震活动特征及当前地震活动形势 地球物理学报，（）：，（）：徐锡伟，郭婷婷，刘少卓，等活动断层避让相关问题的讨论地震地质，（）：，（）：张培震，邓起东，张国民，等中国大陆的强震活动与活动地块中国科学辑，（增刊）：，（），（）：袁道阳，冯建刚，郑文俊，等 青藏地块区大地震迁移规律与未来主体活动区探讨地震地质，（）：，（）：，（）：，：罗浩，徐锡伟，刘小利，等 阿尔金断裂东段的构造转换模式 地质学报，（）：，（）：陈柏林，刘建生，张永双，等嘉峪关断裂晚第四纪新构造活动特征地质力学学报，（）：，（）：（本文编辑：赵乘程）第 卷 第期曹玲玲，等：嘉峪关气氡浓度异常与 年玛多地震关系分析