国外钻孔应变观测的发展现状
-
摘要: 简要论述了钻孔应变观测在人类对地壳运动的研究中扮演的越来越重要的角色. 这种观测已经被美国、日本、秘鲁、澳大利亚、南非、冰岛、意大利等国家用于对板块运动、地震、火山、水库、油田等众多方面的观测研究. 指出这种观测日益成熟,迅速向多分量和综合探头发展. 介绍了根据这种观测结果,地震应变阶、慢地震、地震前兆以及火山喷发预测等方面的研究取得的重要进展.
-
引言
20世纪60年代以来的研究表明,虽然经过世界各国研究者的长期努力,地震预测仍然处在经验性预测阶段. 此阶段的研究特点主要有以下几方面: ① 广泛地积累地震特别是大地震的震例,这要求不断地从大量观测资料中科学地、 全面地分析和提取可能与地震孕育过程相关联的信息; ② 不断地提高对观测信息的物理机制的研究水平,以提高人们分析和识别与地震相关信息的能力; ③ 不断地改进和提高观测技术水平,提高观测的科学性,使获取的观测信息有可靠的科学基础和研究价值.
从这个意义上说,2008年5月12日发生在四川汶川的MS8.0地震,其震例研究具有十分重要的科学价值. 汶川MS8.0地震至今已经过去4年多的时间,在对其进行的科学总结以及各类科学论文有关震前异常的综合分析中(关华平等,2008; 杜学彬,2010; 中国地震局监测预报司,2010)几乎都会提到震前距震中仅35 km的成都(郫县)地电台(以下简称成都台)地电阻率长趋势下降异常,被认为是此次大震前最显著、 最清晰的异常之一. 但这些研究目前还停留在图像识别和对比阶段,缺少深入的分析和论证. 对于这样的8级大地震的震例研究,停留在这个阶段是远远不够的,例如关于汶川震前成都台地电阻率的异常参数,包括起始时间、 形态特征、 变化幅度等重要参数的认定,目前观点并不一致. 特别是该台NE测线地电阻率观测存在明显的年变化,而且地电阻率变化与地下水位的变化有极好的相关性. 但目前所提取的地电阻率的异常变化期间,地下水位的变化态势如何尚无人涉及,从而使地电阻率异常参数判定的科学性会受到一定的质疑. 此外,成都台地电阻率与汶川大地震相关联的异常变化,在长达30年左右的全程观测中的地位,还没有受到应有的重视. 本文将以这些视角,从成都台建台以来的长期观测的全程动态分析开始,深入探讨汶川MS8.0地震震例研究中成都台地电阻率异常的特点及其科学意义.
1. 成都台地电阻率变化全程分析
成都地电台原位于四川灌县,后搬迁至郫县. 本文数据出自中国地震台网中心数据库; 本文所指的全程分析,严格说来是作者根据该数据库的数据进行的分析. 1978年8月是该数据库所能提取到成都台数据的最早期限,因此姑且把1978年8月作为全程分析的开始. 图 1示出成都台测区电测深曲线和观测装置系统(引自中国地震局监测预报司,2010(有删略)),其中布设有NE(供电极距AB=736 m,测量极距MN=226 m)和NW(AB=846 m,MN=270 m)两条测线. 为了进行全程分析,图 2给出了两条测线1979—2008年地电阻率全程观测结果(年均值曲线)以及年变化距平分析结果. 两条测线中,NE测线电阻率观测稳定,信噪比高,地电阻率测值总体变化较大,数十年间年均值最大值与最小值之差达7.43 Ω·m,相对变化幅度为12.4%,但整体变化趋势清晰,有利于异常的趋势分析; 而NW测线地电阻率总体变化较小,数十年年均值最大值与最小值之差为1.25 Ω·m,相对变化幅度仅为2.9%,且电阻率观测长期不够稳定,信噪比低,其原因尚有待深入研究. 此外,NE测线存在比较清晰的年变化,根据多年测值的距平分析结果,其年变化平均幅度约为3%; NW测线年变化幅度小,仅为0.8%,远小于NE测线. 二者相比,可能表明测区下方结构的特殊影响. 因此,本文将主要对NE测线地电阻率的全程变化进行深入分析和讨论.
图 2 成都台地电阻率1979—2008年全程观测结果(a) NE和NW测线年均值; (b) 两测线年变化距平分析比较Figure 2. Results of the whole course of apparent resistivity observation at Chengdu station during 1979—2008(a) Annual mean values in the configurations of NE and NW directions; (b) Comparison between seasonal variations in the configurations of the two directions from statistical analysis图 3给出成都台NE测线1979年1月—2009年9月地电阻率月均值曲线,其中标注了台站300 km范围内所有M>6.5的地震. 从图中可以明显看出地电阻率趋势变化的一些特点. 总体上看,全程地电阻率变化大致可以分为3个阶段: ① 1979—1990年前,趋势上升; ② 1990—1999年前后,平稳变化; ③ 1999年以后到汶川MS8.0地震后数月,在经过了小幅度趋势上升的背景下,从2004年前后开始出现快速大幅度持续下降,下降速度之快、 幅度之大,在该台历史上是前所未有的. 所谓与汶川MS8.0地震有关的成都台地电阻率异常正是发生在这样一个阶段和背景之下.
第一阶段,从建台到1990年以前,这是我国1900年以来第四次地震活动高潮期(1966—1976年)之后全国性地震活动的相对平静期.不过,1981年1月距成都台280 km的四川道孚发生了M6.9地震. 震后,成都台NE测线地电阻率发生了第一次趋势性转折,呈上升态势,且一直持续到1989年9月小金M6.6地震(距成都台仅130 km)为止. 从图像分析来看,成都台的第一阶段变化似乎与该台附近这两次重大震情相关. 此外,正是在这一阶段,成都台开展了一次十分有意义的多极距观测试验. 其分析结果对于今天认识汶川MS8.0地震的震前异常有很好的借鉴,第2节中将有进一步阐述.
第二阶段(1990—1999年),自1989年9月小金M6.6地震发生后,成都台地电阻率趋势发生明显转折,进入了长达10年左右的平稳变化时期. 这一阶段成都台周围200 km甚至300 km范围内地震活动水平较低,没有6级以上地震发生.
第三阶段(1999—2008年),从1999年2月起,经过了将近5年幅度约为2%的小幅度上升期; 从2004年1月开始,地电阻率趋势开始转为下降和加速下降态势. 后期由于下降速率太大,正常年变化的形态开始畸变,到 2007年已经完全无法识别年变化的形态特征,到2008年2月下降的态势似乎终止,前后下降幅度可达6.9%; 之后曲线转为上升态势,到2008年6月,即汶川MS8.0地震后的两个月达到高值,幅度约为3%,与图 2b所示的NE测线多年统计的年变化幅度相当. 笔者分析认为,所谓汶川MS8.0地震前成都台的地电阻率异常正是被包含在这一阶段的过程中. 应当说明,2008年6,7月间,位于NE测线供电电极的区域开始成为灾民的安置场所,大批金属板材活动板房群落就设置在电极附近数米到数十米范围内,对观测产生重大干扰; 大约在2009年9月前后,成都市区到都江堰的轻轨列车开始运行,其干扰使得正常观测更无法进行. 因此,这一阶段的有效资料主要在2008年6月以前.
大地震前后震中附近台站的地电阻率趋势发生明显转折的例子,在以往的大地震震例研究中曾多次提到(Qian et al,1984,1996; 钱家栋,1993),是在地震预测处于经验性的图像识别阶段论证观测参数(包括地电阻率在内)异常变化与地震相关性的重要论据之一. 作为对介质电学属性的测量,地电阻率观测较其它观测(例如地电场、 地磁场观测等)稳定,加之在观测技术上有相当完备的台站、 省级和国家级电学计量标准传递管理规范,通过这些严格的管理规范的实施,可以保障观测仪器性能的及时检定. 此外,在国家台网中心资料库中也未见成都台装置系统变迁变更的记录. 这些都为地电阻率的长期观测对比提供了基本保障(Zhao et al,2011).
本文重点分析图 3所示的成都台NE测线地电阻率第三阶段的变化. 由图 3可以明显看出,地电阻率变化具有两个显著特点,即年变化和大幅度的趋势下降变化.
研究表明,地电阻率年变化,以1年为主要周期,与自然界季节性气象因素(如降雨、 温度、 地下水活动等)以及人类生产和生活受到这些要素的影响而呈现出的季节性行为(如抽水、 灌溉、 农田施肥等),有着很高的统计相关性(钱家栋等,1985; 金安忠,1981). 由这些季节性因素影响而产生的年变化出现在许多地电阻率台站的多年观测记录中,一年四季周而复始,所以也称之为地电阻率的正常变化; 它们的出现一般与孕震过程无关,因此在地震研究中,相对于地震异常而言,这些变化与各类人为干扰一道又统称为“非震变化”. 当然,数学上的相关关系并不等同于物理上的因果关系. 经过深入的研究,在众多因素中,由自然界和人类生产、 生活的综合影响导致的测区浅层潜水的水位变化,才是地电阻率年变化的关键因素和具有实质联系的物理成因(钱家栋,曹爱民,1998).
另一方面,我国还积累了诸多7级以上大地震的震例,在近震中台站上记录到以地电阻率趋势下降现象为主的重要特征. 这些台站一般在震中附近300 km甚至200 km以内,即大震震源尺度的两倍左右的范围内. 由于具有多个震例良好的重复性,这种现象被认为可能与大震的孕震过程相关联. 与地电阻率的正常变化相对照,在地震预报处于经验性阶段的今天,这类变化被统称为“异常变化”(钱家栋等,1985; 钱家栋,1993; 中国地震局监测预报司,2010).
2. 地电阻率法两大理论支柱与成都台地电阻率年变化的分析
由图 3可以看出,汶川MS8.0地震震例中,地电阻率的正常变化和异常变化这两种变化叠加在一起出现在成都台观测记录中. 为了更清晰地提取异常信息,必须充分理解和认识这两类变化的特点. 基于地电阻率长期观测实践、 实验研究和理论研究,形成了以扩展的阿契定律(Archie’s law)和影响系数理论作为地电阻率法两大理论支柱的构想(钱家栋等,1991; 钱家栋,1993). 前者用于讨论介质电阻率的变化与其所处的环境条件变化的关联,解释介质电阻率为什么会变化,属于地电阻率法的微观机理的范畴; 后者用于讨论地电阻率的变化与介质内部电阻率变化的关联,解释测区内部介质电阻率的变化对地表布设的装置系统所测得的地电阻率的影响,或布设在地表的装置系统测得的地电阻率变化对测区下方介质电阻率变化的响应,也就是解决“地上”和“地下”的关系问题,属于地电阻率法的宏观机理的范畴. 它们为本文关于汶川地震成都台异常分析提供了理论依据.
2.1 地电阻率法微观机理: 扩展的阿契定律
根据大量承载岩石电阻率试验研究,钱家栋等(1991)提出将物探电法中阿契定律加以扩展,用以作为解释介质电阻率变化与其环境条件变化关联的重要工具. 物探电法中的阿契定律,是1942年阿契依据实验室实验提出的一个半经验公式,用于表述不同结构之间的岩石(特别是沉积岩)电阻率与其结构参数之间的定量关系(Archie,1942)
式中,α为常数,是与介质类型有关的常数; ρR表示介质整体电阻率; ρ0表示孔隙中孔隙流体的电阻率; φ表示孔隙度,即孔隙体积(VP)在介质总体积(VT)中所占的比例,φ=VP/VT; s表示孔液流体饱和率,即孔隙中流体体积VW所占的比例,s=VW/VP; m,n是两个常数,分别表示与孔隙通道和含水状况相关的结构指数,通常m取值为1.3—1.95,n取值为2(Параснис,1965). 从式(1)可以看出,5个主要参数(ρ0,η,s,m,n)中任何一个发生改变,对岩石整体电阻率都会产生影响. 由此对式(1)两边取对数并求微分,描述阿契定律中5个主要参数(ρ0,η,s,m,n)在环境条件变化时介质电阻率发生的变化. 将其用于讨论介质电阻率变化的微观物理机制,在地电阻率方法中被称为扩展的阿契定律(钱家栋等,1991; 钱家栋,2002; 中国地震局监测预报司,2010; 陆阳泉等,1998),即
对式(1)做适当处理,可得
式中,第一项与孔隙流体电阻率有关,描述介质孔隙中流体电阻率变化对介质整体电阻率变化的影响; 第二项与孔隙流体的含量有关,描述孔隙流体体积变化的影响; 第三项与介质中孔隙体积有关,描述介质中孔隙体积变化的影响; 第四项与结构指数m有关,描述介质孔隙排列的导电通道弯曲性程度改变产生的影响(Qian,1985; 陆阳泉等,1998; 中国地震局监测预报司,2010). 总的来说,扩展的阿契定律为研究介质电阻率变化与其环境条件的关系提供了重要工具. 当所研究的环境条件属于与非震要素关联时,它是研究非震变化的工具; 当所研究的环境条件属于与孕震过程关联时,它是研究异常变化的工具.
研究表明,在4种效应中,孔隙含水量的变化和结构指数变化的效应,影响的量级大. 例如,在第二项中孔隙含水量增大/减小20%,介质的电阻率将减小/增大40%; 在第四项中,设η=0.01,Δm=-0.1,则介质电阻率将减小约46%. 它们是地电阻率法研究中最值得注意的效应(Qian et al,1996; 钱家栋,2002; 中国地震局监测预报司,2010). 前者是年变化即非震效应研究中十分重要的理论工具,因为测区下方潜水水位的涨落,导致表层含水量的增减,是表层介质电阻率变化的主要原因; 后者是研究近震中台站地电阻率异常变化机理的重要工具. 根据Барсуков(1979)的研究,结构指数m的减小,相当于岩石内部的等效导电通道截面增大. 这与大量承载岩石实验后期(在加载应力超过1/2破裂应力水平下)非弹性膨胀阶段新裂隙增长,导致外界水迅速向新裂隙补充,或者新生裂隙内形成表面大面积水膜覆盖,从而记录到岩石电阻率大幅度下降现象的物理过程(Brace,Orange, 1968a,b; Brace,1975)密切相关. 扩展阿契定律中关于结构指数m变化所反映的物理过程,可以为我国7级以上大地震震例中近震中台站地电阻率趋势下降异常的特征,提供比较清晰和有力的解释.
2.2 地电阻率宏观机理: 影响系数理论
由于地学研究的特殊性,几乎所有的地球物理方法,只能通过置于地表或浅层的装置系统的观测来了解地下介质的结构或地下深处发生的物理过程. 因此,研究地表或浅层观测结果与介质内部结构或地下深处发生过程的物理参数的关系,即所谓“地上”和“地下”的关系,是地球科学永恒的主题,也是作为地震科学研究重要分支的地电阻率方法研究所必须关注的内容. 其中,研究地表装置系统测得的地电阻率随时间的变化,与地下介质的结构参数由于各种物理过程引起(真)电阻率变化之间的关联,是地电阻率法中“地上”和“地下”关系的主要研究对象. 地电阻率方法宏观机理的影响系数理论(钱家栋等,1985),正是对两者关系的具体表述. 基于影响系数理论,地表装置系统的地电阻率变化量与测区下方分区均匀介质的(真)电阻率变化量的关系,可以用下式描述:
式中,N表示测点下方分区均匀介质的分区总数; i是介质分区的序号; Δρi/ρi(t)表示在时刻t到t+Δt期间第i个分区介质(真)电阻率的相对变化量; Δρs/ρs(t)表示在时刻t到t+Δt期间地表的装置系统所测到的地电阻率相对变化量; 比例系数Si被称为地电阻率影响系数,它与测区地下介质的电性结构以及地表装置系统的参数有关,表示第i区介质真电阻率的变化对地表设定的装置系统观测到的地电阻率变化的影响程度,反映第i区介质真电阻率的变化对地电阻率变化的贡献.
2.3 成都台地电阻率与地下水位关系
如前所述,在我国地电阻率观测的历史上,对成都台年变化的研究曾经为认识地电阻率年变化的物理机制有过重要的贡献(赵和云,钱家栋,1987; 钱家栋,赵和云,1988; 薛顺章等,1994). 这些研究成果对于汶川MS8.0地震前地电阻率异常的深入分析,有很好的借鉴.
在地下水位变化与地电阻率相关性的研究中,必须区分潜水水位和承压水层水位的作用. 因为承压水层一般来说是饱和的含水层,用于测量水位的水位计实际上是含水层压力的测量器件,测量结果本质上是源头水头高度和含水岩层压力变化的反映,不能表征地电阻率测区下方岩层中含水量的变化,因而承压水层的水位变化并不能作为地电阻率年变化的物理成因(钱家栋,曹爱民,1998).
成都台除地电阻率观测外,还有一项长期与之同步的测区潜水地下水位(水面深度)连续观测和完整记录. 这为开展地电阻率年变化与地下水位的变化关系研究提供了极好的条件. 从微观机理分析的角度,测区下方潜水面水位的上升意味着表层含水量的增加,从而引起表层电阻率的减小; 反之,潜水水位的下降,意味着表层含水量的减小,从而引起表层电阻率的增大. 表层电阻率随着潜水水位季节性的涨落而发生季节性变化,这才是地电阻率年变化的基本原因. 例如图 4a所示嘉峪关地电台的结果,地电阻率与地下水位年变化呈现“镜像”特点(钱家栋等,1985),即两者统计上呈“负相关”. 笔者曾计算过,在图 4a所示的时间段内,嘉峪关台的地电阻率和地下水位之间的相关系数达到-0.94. 然而,研究表明,成都台地电阻率年变化与地下水位年变化的关系,却呈现相反的现象,表现出地电阻率测值与测区水位准同步升降,两者呈现“正相关”,如图 4b所示. 也就是说,随着潜水水位的上升,表层(真)电阻率下降,但地电阻率上升; 反之,潜水水位下降引起表层电阻率上升,而地电阻率却下降. 由于这类现象在我国地电阻率观测中比较少见,加之测区局部介质电阻率减小/增大,而地电阻率反而增大/减小,使得两者变化形态相反的事实从直观上难以理解,从而引起了地电研究人员的高度关注和深入思考. 对成都台这类现象的深入研究,使人们对地电阻率年变化的局域化特点有了全新的认识. 所谓局域化,是指地电阻率年变化除了受季节性因素造成的地下水位变化影响外,测区局部地下介质电性结构和观测系统的耦合是一个制约其影响的重要因素(钱家栋,赵和云,1988),必须应用地电阻率宏观机理的思路才能得到清楚的解释.
基于这一理论思路,钱家栋等(1985)及赵和云和钱家栋(1987)提出了地电阻率观测应当从单一极距向多极距观测系统转变的建议,以便根据地电阻率时间变化的观测结果,通过反演研究,提取介质内部不同区域电阻率变化的信息,作为识别和排除表层干扰与提取深部可能与地震孕育过程相关信息的参考.
2.4 成都台地电阻率年变化与多极距观测试验的启示
1984年,兰州地震研究所薛顺章等与成都台合作,在分析该台电性结构的基础上,在测区沿NE方向架设了5个极距不等的四极对称装置(图 1b),开始了长达7年的地电阻率多极距观测系统的观测实验,受到了国内从事地电研究者的高度重视,开创了用多极距地电阻率观测数据反演介质内部不同部位电阻率变化的新的尝试(赵和云,钱家栋,1987; 薛顺章等,1994; Lu et al,2004). 图 5a由笔者在薛顺章等(1994)给出的成都台多极距观测试验与地电阻率反演结果的基础上重新绘制而成; 图 5b(薛顺章等(1994))示出底层(第四层)真电阻率时间变化与小金地震的关系. 不同文献给出的反演结果有些差异,这与反演方法和反演工作本身的特点有关. 但通过这些反演研究,可以清楚地显示出如下3个特点:
图 5 成都台多极距观测试验与地电阻率反演结果(a) 成都台上部三层真电阻率变化与其地下水位变化的对比; (b) 底层真电阻率时间变化与小金地震Figure 5. Inversion results with apparent resistivity data observed in the multi-separation experiments at Chengdu station(a) Temporal changes of true resistivity at the top 3 layers versus variations of water table at Chengdu station; (b) Temporal changes of true resistivity at the deepest layer before the Xiaojin M6.6 earthquake in 19891)测区下方表层介质真电阻率,由于受地下水动态的影响明显地出现与潜水水位升降相反的年变化,即表层电阻率的减小与水位的上升过程相关; 反之亦然(图 5a). 这与地电阻率微观机理分析的结论并无矛盾,即两者呈现负相关.
2)导致地电阻率与地下水位呈正相关的真实原因是,成都台地电阻率年变化与表层电阻率变化形态相反,即与依据式(4)计算出的表层影响系数S1<0有关. 例如赵和云和钱家栋(1987)给出的结果表明,对于台站原有装置AB=736 m的极距,S1=(Δρs/ρs)/(Δρ1/ρ1)=-0.119. 因此,当Δρ1/ρ1 < 0,即由水位上升引起的表层真电阻率下降的物理过程出现时,Δρs/ρs>0,地电阻率变化量反而增加; 反之亦然. 这可以较好地解释成都台地电阻率年变化与地下水位之间呈正相关的事实. 这一关系之所以奇特,是因为类似成都台的现象很少.
3)从反演结果中(薛顺章等,1994),基本上看不出底层介质真电阻率年变化的态势,但趋势性变化更为明显. 1989年9月在距成都台130 km的M6.6小金地震前,底层真电阻率ρ4呈上升趋势,幅度约为7.4%,而在震后趋势发生明显转折(图 5b). 这种变化态势似乎与此次地震孕育过程相关. 从图 5还可进一步看出,将底层真电阻率的变化趋势量7.44%与薛顺章等(1994)给出的地表装置系统(AB=736 m)地电阻率测值,在消除年变化影响后的趋势变化幅度2.26%相比较,其比例系数约达3.28 倍. 此结果表明,成都台NE测线地电阻率变化主要由两部分影响因素组成: 一部分源自表层,即地下水位的影响; 一部分源自底层,可能与孕震过程的影响相关联.
2.5 成都台地电阻率与地下水位年变化幅度的定量关系
基于1985年6月—1990年12月的距平分析,可以分别得到成都台NE测线地电阻率和测区地下水位年变化平均形态(图 6),并在此基础上做一些定量的讨论. 从图 6 可以看出,两者之间同步升降趋势明显,均呈夏高冬低的态势. 地电阻率高值点位于8—9月间,地下水位的高值点则位于6—8月,相对于地下水位,地电阻率相位略有滞后; 两者的低值点均在1—2月. 从幅度对比来看,在此期间地下水位年变化幅度为1.04 m,而相应的地电阻率年变化幅度约为3.19%.
3. 汶川MS8.0地震前成都台地电阻率异常变化与地下水位的综合分析
3.1 分析方法: 特征点的认定
由于成都台地电阻率与测区地下水位存在密切的关系,因此,对于汶川MS8.0地震前成都台地电阻率的异常分析,应当同时分析该台测区地下水位的变化,即从地电阻率与地下水观测结果综合分析的角度,判定成都台第三阶段地电阻率变化的异常特征(图 7). 图 7a为2003—2008年间地电阻率与地下水位对比曲线图.
为了进行定量分析,本文选择年变化的低值点作为分析异常变化特征的判定基础. 因为低值点在冬季,没有降雨等不稳定因素对地下水位的影响,比选择夏季的高值点的变化更为可靠. 更重要的是,从图 7a可以看出,这一阶段后期地电阻率快速下降的趋势,在2008年2月开始转折,正是冬季地电阻率的低值点. 因此,无论它是否意味着年变化形态的恢复,以冬季低值点的地电阻率测值描述地电阻率的趋势变化是有意义的. 另一方面,地电阻率年变化的高值在8—9月间,2007年的高值点已经被急剧的下降趋势所“淹没”; 而2008年的高值点,又是在5月12日大地震之后,跨越了地震时段因而无从对比. 这使依据高值点作为特征点的选择也面临着不少困难. 各个特征(时间)点上地电阻率和地下水位的数据如表 1所示. 表中t表示特征(时间)点,h和ρS分别表示相应(时间)点地下水位(水面深度)和地电阻率测值,Δh和ΔρS分别表示两个时间点之间的时段内h以及ρS的变化量(时间段的起终点标识在数字后的括弧中),ΔρS/ρS则表示所选时段内ρS的相对变化量. 地电阻率和地下水位的特征点位置稍有不同,但差异不大,不会影响本文讨论.
表 1 成都台第三阶段地电阻率与地下水位变化参数对比Table 1. Parameters of variations of apparent resistivity and water table in the third stage of the whole course of observation at Chengdu station直观上注视图 3和图 7a中地电阻率的变化曲线,可以看出汶川MS8.0地震前成都台地电阻率的趋势发展,经历了从1999年2月前后的缓慢上升阶段和2004年2月起转为下降阶段,并在2006年2月—2008年2月期间出现了大幅度下降的变化过程,前后约为9年. 震后的2008年异常形态基本恢复. 另一方面,从图 7a可以明显看出,在所认证的地电阻率异常出现大幅度快速下降期间,地下水位年变化形态基本稳定,未出现明显的改变. 因此,成都台地电阻率年变化畸变可以作为地震异常的重要旁证.
3.2 成都台地电阻率异常特征的定量分析: 关于孕震过程影响的推断
然而,细致分析图 7a和表 1,可以将2004年2月—2008年2月大致分为前后两个时段(以图 7a中竖直的粗虚线为界). 前期2004年2月—2006年2月,水位下降1.04 m,地电阻率下降1.3%(即0.82 Ω·m); 后期2006年2月—2008年2月,水位上升0.4 m,地电阻率下降5.6%(即3.39 Ω·m). 即前期水位下降趋势与地电阻率下降同步,而后期两者趋势相反(如图 7a中细虚线所示).
因此,基于成都台地电阻率与地下水位的关系呈正相关的特点,对照图 6中关于1985—1990年观测试验的对比结果(即将地下水位1.04 m的变幅与地电阻率3.19%的变化幅度相等效),依线性关系处理,估算出1.04 m和0.4 m 的水位变化对地电阻率将产生的附加影响,分别约为3.19%和1.23%. 也就是说,可以推断前期地电阻率变化中3.2%左右(减小)、 后期1.2%左右(增大),应归因于地下潜水水位变化引起的表层电阻率变化,分别叠加在实际的观测结果中,它们不属于孕震过程的影响,应当加以修正. 扣除潜水水位的影响后,前期阶段地电阻率应有1.9%(--1.3%-(--3.2%))的上升异常呈现; 而后期应有--6.8%(--5.6%-1.2%)的下降异常,作为底层介质电阻率变化的贡献. 修正后的结果如图 7b中上面一条实线所示. 这一结果与从图 7a中直观看到的地电阻率异常过程有所区别: 直观上看,下降的异常态势始于2004年2月,并持续4年(图 7a); 而在定量地分析地下水位影响之后,得到的异常图像是,先小幅上升,后转为大幅度下降,下降的起始点应该定在2006年2月前后,持续约2年左右(图 7b).
如本文2.4节所述,根据20世纪80年代成都台的多极距观测试验,地表装置系统观测的地电阻率变化与测区底层真电阻率变化的关系,约为3.28倍. 可以推断底层介质真电阻率受孕震过程的影响,在异常后期的两年中出现的下降幅度可以达到22.3%(6.8%×3.28)的量级. 当然,由于多极距观测已经结束多年,不能像20世纪80年代中期那样,通过反演求得测区底层介质(真电阻率)变化的具体过程,因此,只能用这种方式进行推断. 测区下方底层介质电阻率减小达到20%这样的量级,可以用扩展的阿契定律中结构指数m的减小来解释(中国地震局监测预报司,2010),即推断这一时段测区底层介质正处于类似承载岩石临近大破裂后期过程中,微裂隙的广泛发育、 含水介质中水向新生裂隙的掺入、 裂隙排列方向的趋同化以及原有导电通道的弯曲性效应逐步改善,导致介质底层(真)电阻率大幅度下降(Brace,Orange, 1968a,b; Brace,1975; Барсуков,1979). Ma(2002)曾采用3维边界单元法,计算两层结构条件下底层介质中存在的低阻块体的尺度逐步增大对地表地电阻率变化的贡献,模拟深层介质中裂隙发育和排列方向趋同化的过程,得出了同样的结论. 成都台位于龙门山断裂带附近,距离汶川MS8.0地震震中仅35 km,相对于汶川地震断层320 km的尺度,可以认为处于该地震的震源区内. 这里给出的测区底层介质电阻率变化的结果,可以作为汶川MS8.0地震孕震过程研究,特别是临近大地震发生的后期阶段研究的重要约束条件之一,有着重要的科学价值.
4. 结论
1)成都(郫县)地电台自建台以来,NE测线地电阻率的全程观测值,按趋势变化的不同态势可以分为3个阶段,其趋势转折的时间与成都台周围300 km范围内仅有的3次强震(M>6.5)发生的时间密切关联,即1981年1月道孚M6.9地震、 1989年9月小金M6.6地震和2008年5月汶川MS8.0大地震. 汶川MS8.0地震是建台以来发生在该台最近、 震级最大的一次地震,也是成都台全程观测中异常变化最为突出的一次地震.
2)成都台与汶川MS8.0地震有关的地电阻率异常,属于全程观测的第三阶段,即1999—2008年,持续时间约为9年. 由于NE测线地电阻率与测区地下水位存在显著的(正)相关现象,综合对比分析显示,在所认证的地电阻率异常出现大幅度快速下降期间,地下水位年变化形态基本稳定,未出现明显的改变,从而可以将成都台NE测线地电阻率年变化畸变和消失作为地震异常的重要旁证. 通过对这一阶段地电阻率趋势变化的影响进行的定量分析,可以认证的异常持续时间大约为2004年2月—2008年2月,并分为前、 后两个时段. 前期2004年2月—2006年2月,呈小幅上升态势,幅度为1.9%; 后期2006年2月—2008年2月,呈现大幅度快速下降态势,幅度达-6.8%,其幅度和速率,建台30多年以来所未见.
3)推断异常后期成都台地电阻率大幅度下降期间,其测区下方底层介质电阻率减小的量级可达22%,与物理试验中承载岩石样本临近大破裂后期的性状(微裂隙的广泛发育、 含水介质中水向新生裂隙掺入、 裂隙排列方向的趋同化导致岩石电阻率大幅度下降)十分相似. 考虑到成都台所处的特定地理位置,位于大地震的震源区内,这些结果可以作为成都台直接记录到孕震后期过程效应的一个佐证,不仅成为我国大地震震例中地电阻率异常研究库中又一重要的科学资料,也可以作为汶川MS8.0地震孕震物理过程研究中的重要约束条件之一,具有重要的研究价值.
4)成都台地电阻率异常的全面分析,其科学价值不仅仅局限于与以往震例的相似性,提高了地电阻率异常变化在地震预测中的实用性,其重要的科学价值还在于,在以往多次大地震震例研究中,对于地电阻率年变化被趋势性下降过程所扭曲而发生畸变的现象,只能停留在对异常进行形态分析的水平. 也就是说,在所有的震例中凡涉及地电阻率年变化的畸变现象的解释,都是建立在当地导致地电阻率年变化的地下水位的季节性周期变化规律本身没有发生任何改变的“假设”前提之下. 本文关于地电阻率和实测地下水位的综合研究,使观测到年变化畸变的成都台的异常分析,有了实际科学论证的条件. 这表明,开展有针对性的综合观测,例如开展地电阻率测区地下水位的同步观测、 台站电性结构的详查和解译,以及在可能的条件下采取多极距观测等,是实现地电观测科学化进程的重要措施之一.
相对于地电阻率数据,本文收集的地下水资料尚不够完备,有关的讨论只能局限于笔者于2008年赴台站落实异常期间收集到的地下水位资料. 其余资料,早期部分记录在纸介质上,当时尚未能及时入库.
在收集资料的过程中得到中国地震台网中心叶青和成都台张保全等同志以及四川省地震局资料室有关同志的大力支持,在此表示衷心感谢.
-
-
期刊类型引用(75)
1. 高熹微,张宇,郭宇鑫,张骞,李玉,殷翔,胡超. 风力发电设备检修引起地电阻率变化的分析研究. 中国地震. 2025(01): 114-125 . 百度学术
2. 娄晓宇,张宇,张兴国,李钟慧. 基于ADS1256的地电阻率仪数据采集系统设计. 电子器件. 2024(04): 891-896 . 百度学术
3. 张宇,王兰炜,张世中,张兴国,胡哲. 地电阻率观测中信噪比与观测精度. 地震地质. 2024(04): 955-971 . 百度学术
4. 郭建芳,佟鑫,郭骄,阴佳诗,巩洪学,王枭一. 昌黎地电台地电阻率与气象因素、地下水位的相关性分析. 地震地磁观测与研究. 2024(05): 107-113 . 百度学术
5. 张宇,柯浩楠,娄晓宇,胡颢迪. 2023年12月18日甘肃积石山6.2级地震前地电阻率变化分析. 地震. 2024(04): 82-96 . 百度学术
6. 马小溪,王兰炜,朱涛,张宇,张兴国. 地电阻率交流观测中的电感性耦合效应特征模拟和分析. 地球物理学报. 2023(01): 289-300 . 百度学术
7. 王兰炜,张宇,张兴国,胡哲. 地电阻率观测中地铁杂散电流特征. 地震学报. 2023(02): 285-301 . 本站查看
8. 王同利,邢海林,崔博闻,王丽红,孙小入. 小极距井下地电阻率观测对地表环境影响的抑制分析——以北京地区台站为例. 地震. 2023(02): 1-13 . 百度学术
9. 张宇,王兰炜,胡哲,张世中,张兴国,娄晓宇. 加权平均算法在地电阻率日均值计算中的应用. 地震. 2023(02): 25-37 . 百度学术
10. 王同利,崔博闻,王丽红,颜晓晔,李菊珍. 2008年汶川M_S8.0地震前后地电阻率时空变化. 地震研究. 2022(01): 66-74 . 百度学术
11. 叶青,王晓,杜学彬,解滔,范晔,周振贵,刘高川. 中国地震井下地电阻率研究进展. 吉林大学学报(地球科学版). 2022(03): 669-683 . 百度学术
12. 文雯,毛先进. 地电阻率小偶极装置观测的特性. 华北地震科学. 2022(03): 50-54 . 百度学术
13. 张远富,赵斐,张飞丽,张旭燕,尹欣欣. 地电阻率观测装置布设与金属管线类干扰的关系研究——以平凉台为例. 地震工程学报. 2022(04): 945-952+961 . 百度学术
14. 王晓,叶青,王军,范晔,樊俊屹,刘高川. 地震地球物理观测台网地电数据集. 地震地磁观测与研究. 2022(03): 134-139 . 百度学术
15. 崔博闻,王同利,李然,胡毅涛,王丽红,李菊珍. 通州台井下地电阻率测深初步分析. 地震工程学报. 2022(05): 1176-1184 . 百度学术
16. 王晓,叶青,余丹,陶志刚,范晔,刘高川. 2022年6月1日芦山M_S 6.1地震前地电阻率异常分析. 地震地磁观测与研究. 2022(04): 35-40 . 百度学术
17. 张宇,王兰炜,胡哲. 地电阻率观测中两种电极接地电阻测量方法及其对比研究. 地震学报. 2022(06): 1049-1060 . 本站查看
18. 王同利,崔博闻,李然,王丽红,胡毅涛,邢海林,李经山. 通州地电阻率对唐山地区中强地震的映震能力分析. 中国地震. 2022(03): 455-463 . 百度学术
19. 于磊,陈俊,何康,刘泽民,黎哲君. 合肥形变台地电阻率NS向数据异常的解释. 地震. 2022(03): 180-191 . 百度学术
20. 朱涛. 地震电阻率实验研究新进展及展望. 地球与行星物理论评. 2021(01): 61-75 . 百度学术
21. 王兰炜,张宇,张兴国,张世中,胡哲. 井下地电阻率观测装置稳定性研究. 电子学报. 2021(06): 1167-1177 . 百度学术
22. 吴鹤帅,罗娜,赵志远,张波,尹小兵. 降雨对大柏舍地电台地电阻率观测的影响. 地震地磁观测与研究. 2021(03): 173-179 . 百度学术
23. 王丽,叶青,周康雅,丁烨. 南京地电阻率长趋势变化与区域应力调整的关系研究. 地震工程学报. 2020(01): 143-149+214 . 百度学术
24. 黄明威,张琪,李梦莹,杨牧萍. 新城子井下地电阻率观测影响系数分析. 防灾减灾学报. 2020(01): 49-56 . 百度学术
25. 王同利,崔博闻,叶青,李菊珍,王丽红,童琼. 九寨沟M_S7.0地震地电阻率变化时空演化分析. 地球物理学报. 2020(06): 2345-2356 . 百度学术
26. 沈红会,祝涛,孙春仙,李鸿宇,张桂霞. 江苏地电阻率预测指标研究. 中国地震. 2020(01): 115-125 . 百度学术
27. 肖武军,解滔,胡玉良. 冬奥会保障项目井下地电阻率观测专用电缆设计与测试. 中国地震. 2020(01): 146-152 . 百度学术
28. 高曙德. 四川九寨沟7.0级地震前震情跟踪概述及震后总结. 地球物理学进展. 2020(04): 1250-1260 . 百度学术
29. 邹广,陈亮,牛中华. 2020年6月26日新疆于田M_S6.4地震地电场异常分析. 内陆地震. 2020(03): 310-316 . 百度学术
30. 解滔,于晨,王亚丽,李美,卢军. 基于断层虚位错模式讨论2008年汶川M_S8.0地震前视电阻率变化. 中国地震. 2020(03): 492-501 . 百度学术
31. 解滔,于晨,卢军. 开展小极距井下地电阻率观测的可行性分析. 中国地震. 2019(01): 14-24 . 百度学术
32. 肖武军,解滔,张尧. 晋冀蒙交界及附近地区小极距井下地电阻率观测装置设计. 中国地震. 2019(01): 134-143 . 百度学术
33. 毛先进,段炜,庄儒新,杨玲英,赵晋民. 井下地电阻率观测中布极参数的确定方法. 地震研究. 2019(01): 96-101 . 百度学术
34. 王燚坤,何康,李军辉,葛计划,陶媛. 合肥台地电阻率NS向趋势异常变化原因分析. 华南地震. 2019(02): 58-64 . 百度学术
35. 于磊,陈俊,刘泽民,戚浩,孙军,王哲. 合肥形变台地电阻率观测映震能力分析. 华南地震. 2019(02): 79-85 . 百度学术
36. 樊晓春,解滔,吴帆,袁慎杰. 井下地电阻率观测影响系数分析——以江宁地震台为例. 中国地震. 2019(02): 347-358 . 百度学术
37. 王燚坤,何康,李军辉,汪小厉,孙鸿博. 合肥中心台地电阻率地铁干扰分析与抑制处理. 华北地震科学. 2019(03): 67-74 . 百度学术
38. 张国苓,茅远哲,乔子云,贾立峰,罗娜,张素欣. 降雨对昌黎台地电阻率的干扰分析. 地震. 2019(04): 110-117 . 百度学术
39. 李姜,乔子云,张国苓,贾立峰,罗娜. 华北地区地电阻率归一化速率异常分析. 地震地磁观测与研究. 2019(06): 62-71 . 百度学术
40. 钱家栋,张学民,王亚璐,李雪浩. 2018年汶川M_S8.0地震前成都台NE向地电阻率趋势异常的数值模拟. 地震. 2018(02): 107-116 . 百度学术
41. 解滔,刘杰,卢军,李美,姚丽,王亚丽,于晨. 2008年汶川M_S8.0地震前定点观测电磁异常回溯性分析. 地球物理学报. 2018(05): 1922-1937 . 百度学术
42. 马钦忠. 与2008年汶川M_S8.0地震有关的成都台地电场异常信号. 地震学报. 2018(03): 351-363 . 本站查看
43. 王燚坤,陈俊,李发,李军辉,何康,裴红云,李罡风. 安徽省地电阻率台网观测质量分析与改进技术. 华北地震科学. 2018(02): 58-65 . 百度学术
44. 刘川琴,金艳,裴红云,李军辉. 安徽省数字化地电阻率观测数据的干扰分析与处理方法. 地震. 2017(02): 167-178 . 百度学术
45. 安张辉,高悦,范莹莹,刘君,王建军. 地电潮汐力谐振共振波研究回顾与展望. 地球物理学进展. 2017(01): 1-7 . 百度学术
46. 刘素珍,李自红,宁亚灵. 大同地电阻率NS测向异常分析. 地震地磁观测与研究. 2017(04): 105-109 . 百度学术
47. 毛先进,杨玲英,赵晋民,段炜,庄儒新. 井下地电阻率观测的探测深度初步研究. 地震学报. 2017(02): 230-238 . 本站查看
48. 杜学彬,孙君嵩,陈军营. 地震预测中的地电阻率数据处理方法. 地震学报. 2017(04): 531-548 . 本站查看
49. 张斌,朱涛,周建国. 岩石电阻率图像及各向异性变化的实验研究. 地震学报. 2017(04): 478-494+633 . 本站查看
50. 沈红会,王丽,王维,李鸿宇,孙春仙,叶碧文. 地电阻率长趋势变化及其预测意义. 地震学报. 2017(04): 495-505+633 . 本站查看
51. 王立凤,赵国泽,陈小斌,汤吉,蔡军涛,詹艳,韩冰,程远志,赵凌强,肖骑彬. 芦山M_S7.0地震余震期间大地电磁视电阻率变化. 地震学报. 2017(01): 64-77+155 . 本站查看
52. 高曙德,罗维斌,苏永刚,梁子斌,李旭升,韩英,孙海洋. 编码源地电阻率观测试验. 地震学报. 2017(04): 506-519+633 . 本站查看
53. 高曙德,郭安宁,王军燕,武善艺,杨晓鹏,张博,赵斐,李旭升. 四川九寨沟Ms7.0地震前地电阻率的异常变化. 地震工程学报. 2017(04): 645-651 . 百度学术
54. 高曙德. 深井地电观测技术在地震监测中的应用探讨. 地球物理学进展. 2016(05): 2078-2088 . 百度学术
55. 安张辉,詹艳,陈小斌,姜峰,高悦. 滑动自相关方法在地电阻率观测资料分析中的应用初探. 地震地质. 2016(04): 1019-1029 . 百度学术
56. 高曙德,赵国泽,汤吉,苏永刚,詹艳,王立凤. 2008年中国大陆6.0级以上地震前后电磁脉冲异常现象. 地震地质. 2016(04): 987-1004 . 百度学术
57. 刘素珍,李自红,张淑亮,殷志刚,李宏伟,王霞. 2014年大同地电阻率NW测向异常分析. 地震地磁观测与研究. 2016(02): 104-110 . 百度学术
58. 高曙德,狄国荣,苏永刚,高原,张磊,梅东林,陈彦平. ELF电磁仪远程控制和数据传输及编译的实现. 地震工程学报. 2016(03): 471-477 . 百度学术
59. 郭建芳,周剑青,佟鑫,王晶晶,张思远. 河北昌黎台多极距观测数据和地下结构分析. 地震工程学报. 2016(S2): 190-196 . 百度学术
60. 张斌,周建国. 岩石电阻率的实验研究. CT理论与应用研究. 2015(02): 239-249 . 百度学术
61. 杜学彬,刘君,崔腾发,范莹莹,安张辉,闫睿,王丽. 两次近距离大震前成都台视电阻率重现性、相似性和各向异性变化. 地球物理学报. 2015(02): 576-588 . 百度学术
62. 解滔,卢军. 地电阻率三维影响系数及其应用. 地震地质. 2015(04): 1125-1135 . 百度学术
63. 王兰炜,张宇,张世中,颜蕊,王子影,张兴国,胡哲. 我国井下地电阻率观测技术现状分析. 地震地磁观测与研究. 2015(02): 95-102 . 百度学术
64. ZHAO GuoZe,BI YaXin,WANG LiFeng,HAN Bing,WANG Xiao,XIAO QiBin,CAI JunTao,ZHAN Yan,CHEN XiaoBin,TANG Ji,WANG JiJun. Advances in alternating electromagnetic field data processing for earthquake monitoring in China. Science China(Earth Sciences). 2015(02): 172-182 . 必应学术
65. 赵国泽,BI YaXin,王立凤,韩冰,汪晓,肖骑彬,蔡军涛,詹艳,陈小斌,汤吉,王继军. 中国地震交变电磁场观测数据处理技术新进展. 中国科学:地球科学. 2015(01): 22-33 . 百度学术
66. 沙宏武. 南京台地电场观测中电极问题的判定探讨. 信息系统工程. 2015(06): 124-126 . 百度学术
67. 沈红会,李鸿宇,袁慎杰,缪阿丽. 江宁台地电场地地铁干扰试验结果分析. 地震. 2014(03): 125-131 . 百度学术
68. 戴苗,刘坚,蒋玲霞,赵凌云,张念. 巴东地震前介质电性变化特征. 大地测量与地球动力学. 2014(03): 48-51 . 百度学术
69. 王兰炜,张世中,张宇,胡哲,顔蕊. 井下地电阻率观测中装置系数的计算——以天水地震台井下观测为例. 工程地球物理学报. 2014(01): 50-59 . 百度学术
70. 张建国,陈化然,张国岺,贾立峰,乔子云. 电磁方法在地震监测预测中的新思路探讨. 工程地球物理学报. 2014(03): 426-430 . 百度学术
71. 徐锡泉,高昌志,王亮. 内蒙古宝昌台地电阻率长期观测数据研究. 地震工程学报. 2014(02): 405-412 . 百度学术
72. Xie Tao,Lu Jun,Ren Yuexia,Zhao Minyuan. Analysis on Apparent Resistivity Variations of Garzê Station before the 2013 Lushan M_S 7. 0 Earthquake. Earthquake Research in China. 2014(03): 388-402 . 必应学术
73. 康云生,安海静,马可兴,谭大成. 天水地电阻率地表与井下多种观测方式的试验分析. 地震工程学报. 2013(01): 190-195 . 百度学术
74. 曹玲玲,王宗礼,李春燕. 岷县漳县6.6级地震前天水花牛水氡异常可靠性研究. 地震工程学报. 2013(04): 800-807 . 百度学术
75. 解滔,卢军,任越霞,赵民渊. 2013年芦山M_S7.0地震前甘孜台地电阻率变化分析. 震灾防御技术. 2013(04): 377-388 . 百度学术
其他类型引用(5)
计量
- 文章访问数: 1101
- HTML全文浏览量: 35
- PDF下载量: 54
- 被引次数: 80