甘肃清水温泉井与临夏井水温在汶川地震前后的异常变化特征

摘要：分析了甘肃省清水温泉与临夏水温观测井自2007年以来的水温观测数据，得出这两口井在2008年汶川80级地震前后有明显的异常变化。进一步用一元线性回归法、斜率K值法和从属函数法对每口井的水温异常信息进行了定量分析，更好地提取了地震异常信息。

关键词：水温观测；前兆异常；一元线性回归；斜率K值；从属函数；汶川地震

中图分类号：P315723文献标识码：A文章编号：1000-0666（2013）03-0269-06

0引言

地球是一个巨大的热库，由内部向地表释放热量，在地球表面产生热流。地热一般是指地壳热状态，在地震孕育与发生过程中，一方面与地壳的热状态密切相关，另一方面又会引起地壳热状态的改变，因此可通过对反映地壳热状态的井水温度或井孔围岩温度的观测，捕捉地震前地壳热状态异常变化信息，进行地震监测预报的探索。目前地震部门的地下水温度观测传感器置于井水中，观测水温变化，大部分的台站都采用具有高精度的石英温度计。

自从开展水温观测以来，已经积累了较多水温地震前兆异常实例。由于异常多出现在震前几天至几十天时间段内，为短临异常；且异常幅度大，一般为正常起伏度的几至几十倍，异常显著而易识别，因此在地震短临监测中表现出特有的优势（付子忠，1988；车用太，金鱼子，1997）；刘耀伟等，2008）。为我国的地震监测与预测，尤其在短临预测与实现地震预报提供了重要数据。但对于水（地）温异常特征，目前主要是水（地）温异常与地震的对应关系的定性分析，定量的研究则较少。

2008年5月12日在青藏块体东部边缘发生汶川80级地震，地震前，甘肃省清水温泉和临夏井都记录到非常显著的水温异常。笔者对这两井的整点值数据进行一元线性回归分析，然后对日均值数据进行平滑滤波处理，分别用斜率K值法、从属函数法等方法，分析两口井水温异常信息的相似性与差异性。

1台站基本情况及数据处理111观测井简介（1）清水温泉井

清水温泉位于清水县城北东8 km的汤浴河谷内（3445°N，10613°E），海拔高程1 445 m。地势北东高，南西低，属低山丘陵。其构造位置属于陇西系南段、秦岭北缘大断裂的北侧，六盘山断裂带西南侧（图1）。台基岩性为通渭—清水断裂晚古生代花岗岩，该温泉是深部热水沿岩石裂隙运动，在一定地质构造条件下出露于地表而形成的。清水温泉井深4433 m，水温观测仪器为SZW-1A型数字式温度计，水温探头位于井下185 m处。2007年6月仪器入网，产出正常数据产品。

（2）临夏井

临夏井位于临夏市东北方向约7 km的折桥乡（3563°N，10426°E），海拔1 830 m。地形开阔平坦，周围没有明显干扰。第四纪覆盖层厚10～20 m，为亚粘土卵石，下伏第三系基底花岗岩。井深2001 m，2006年11月开始水温数字化观测，2007年入网，2007年6月产出正常数据。观测仪器为SZW-1A型数字化水温仪，水温探头位置在井下185 m处的第三系花岗岩中。图1研究区域构造背景及井点分布图

Fig1Regional tectonic background and distribution

of well points in study area12地热正常动态

赵刚等（2009）对全国近277个观测井的地热正常动态进行了系统研究，将其划分为降温型、升温型、波动型、稳定型、跳变型、长周期型6种基本类型。据统计，在整个地下流体台网中，降温型、升温型和波动型占70 %左右。在了解其正常动态特征基础上，才能够提取地震前兆异常信息。井水温度的背景值，一般取决于井孔所在地的大地热流（q）或地温梯度（dT/dz）。一般来说，热流密度大的地区地温高，地温高的地方水温也高。根据地温成因，一个地区不同深度的地温状态是不同的。一般由地表到一定深度段上，地温受太阳辐射热的影响而随时间有规律起伏，如随昼夜、季节变化等，该变化带称为变温带。变温带的厚度一般为十几米至几十米。变温带以下有一个面，该面上的深度不受太阳辐射的影响，一般常年不变，称为恒温带。恒温带以下水温，主要受地温梯度的控制，即随着深度（Z）、水温（T）升高。

清水温泉井水温自架台产出数据以来基本呈现趋势上升变化（图2a），年变幅度达到005 ℃，长期动态特征属于漂移上升型，探头位置处于恒温带下，所测为深部水温变化，不受地表温度影响。日变幅在002 ℃左右，月变幅约为003 ℃。

临夏井水温自2007年6月产出数据以来基本呈现上升趋势变化（图2d），年变幅度达到003 ℃，长期动态特征属于漂移上升型，探头位置很深，处于恒温带下，所测温度为深部水温变化，不受地表温度影响。日变幅在001 ℃，水温观测无潮汐效应，月变幅约为002 ℃，该井日变幅较大，几乎为年变幅的30%，和清水温泉井水温相似。

1地震研究136卷第3期蒲小武等：甘肃清水温泉井与临夏井水温在汶川地震前后的异常变化特征13异常提取方法

（1）直线回归分析

变量之间的关系可分为两种：一种是函数关系，另一种叫相关关系，又叫统计相关。水温观测所蕴含的信息成分比较复杂，非单一因素引起，是多种成分的混合、叠加，因而无法找到一个确定的函数关系，但数据的变化还是有一定的规律可循。假定水温是一个随时间变化的量，并与时间变化有相关关系，且按照：y=f（x）（1）的函数形式相关。通过实测的观测数就可以计算出几个数值，用直线来大致表示两个变量的关系：y=a+bx，（2）其中y为水温，x为时间，直线斜率b为回归系数，a为常数，a，b的值计算公式为b=∑N1t=1xtyt-11N∑N1t=1xt∑N1t=1yt1∑N1t=1x2t-11N∑N1t=1x2t，（3）a=y-bx（4）在直线回归分析公式中，自变量x为时间序列，因变量y为水温测值。

（2）斜率K值法

斜率K值法原理上就是趋势速率法，在本文中求取的K值就是水温日变化速率。如果水温数据以一个恒定的速率上升，则计算得出的斜率K值为一个相对固定的值，那么K值曲线则是一条直线；若高于或低于这个值，K值曲线偏离直线的时段都属于异常。

清水温泉井和临夏井长期动态类型都属于长期漂移上升型，但对于一个一直随时间呈趋势上升变化的量，单从原始曲线上无法较直观的观察到这种变化，因此计算两个井水温动态变化的斜率，可以提取地震异常。由于清水温泉、临夏井日变幅较大，对两台日均值数据先进行11日滑动处理，再以11日滑动均值计算斜率。通过反复计算，笔者采用步长为30日的斜率K值，以抑制短周期跳变信息，突出中长周期成分信息：Ki=Ti+m-Ti1m，i=1，2，3……（5）式中，K为斜率，T为水温，i为时间序列，m为步长，取值为30。

（3）从属函数法

前兆观测资料的异常变化形态虽然是多种多样的，如突跳、转折、加速等，但其本质都是观测曲线随时间变化的斜率变化的量。从属函数就是表征一条曲线变化的斜率和跳动程度（叶秀薇，2004）。通过对资料的分析比较，用11点滑动平均值作为数据序列计算的效果较好，明显优于直接使用月均值提取异常的方法（杨兴悦等，2006）。其从属函数为μi=1+α1kiri-1（6）式中，ki为观测值变化斜率；ri 为滑动平均值M（t）与时间t的相关系数，它反应了观测序列内在质量的好坏；α为经验常数。

通过多次程序计算，清水温泉井α取值为0000 3，拟合窗长m值取9，计算出的μ值较佳；临夏台α取值为0000 4，拟合窗长m值取9，计算出的μ值较佳。图2清水、临夏两井水温数据经过处理后的曲线图

（a）清水井11日平滑日均值；（b）清水井斜率K值；（c）清水井从属函数；（d）临夏井11日

平滑日均值；（e）临夏井斜率K值；（f）临夏井从属函数

Fig2The curve of the processed water temperature data recorded by Qingshui and Linxia Well

（a）smooth day mean for 11 days in Qingshui Well；（b）slope Kvalue in Qingshui Well；（c）subordinate function in Qingshui Well；

（d）smooth day mean for 11 days in Linxia Well；（e）slope Kvalue in Linxia Well；（f）subordinate function in Linxia Well2清水温泉井水温异常分析121原始曲线分析2007年6月2日至2008年12月31日，清水温泉井水温一直处于快速上升变化状态，在汶川80级地震前后，没有明显的突变、转折等异常形态出现，但水温上升变化很显著，呈前期加速、后期减缓的异常变化（图2a），且与对数函数曲线形态很相似，这种在汶川地震前后水温曲线呈对数变化形态的水温井还有在甘肃东南部的成县井、武都樊坝井、ZK801井等，由于篇幅有限，此文不再详述。2009年9月1日至2011年3月11日水温变化相对以前速度放慢，日本90级地震发生后，水温测值明显回落，持续了约70天后，于2011年5月20日后恢复正常，数据以比较缓慢而相对固定速度上升，曲线形态呈现直线斜形上升变化。

22直线回归分析

按年对两井整点值数据进行线性回归分析，计算出回归方程（表1），2007年直线斜率b值最大，为0000 2，即水温数据基本以每小时0000 2℃的速率在上升；2008年降为0000 1℃，在b值变化大的2007～2008年间，发生了汶川80级地震，2009年、2010年b值相同且较小，2011年b值变为很小的负值，且在2011年3月11日的发生了本90级大地震。2012年b值为正，但b值很小，仅为0000 01。表1井水温直线回归方程a值与b值

Tab1The avalue and bvalue in linear regression

equation of well water temperature

23斜率K值法

2007年6月2日至2008年12月31日，K值以稳定的速率由高值逐渐下降（图2b），但相对其后的K值，此时段K值仍处于高值区，说明在此时段水温总体上升较快，但上升的速率是逐渐减小的；2009年9月1日至2011年3月11日K值下降速度放缓；3月11日日本90级地震发生后，K值出现一个明显的大幅度下降—回升的变化形态；2011年5月20日后K值没有出现明显的上升、下降等变化形态，基本稳定在0000 01附近。

24从属函数法

2007年6月7日至2009年1月23日，从属函数μ值呈现明显的高值变化（图2c），汶川地震前出现的高值变化更多、更剧烈，μ最高值出现在2007年6月初，为062，震后有所衰减，但也很明显，异常变化一直持续到2009年1月24日，μ值接近为0。2009年11月26日至2010年4月5日，μ值又出现高值，最大幅度出现在2009年12月30日，μ为058，随后于2010年4月14日发生了青海玉树71级地震，μ值异常变化显著。在2011年3月11日日本90级地震前6天，μ值在2011年3月4日、3月5日出现明显大于背景值的高值变化，约为035；2012年8月12日新疆和田62级地震前8天μ值也出现了034的较高数值变化。

3临夏井水温异常分析131原始曲线分析临夏水温自2007年6月15日至2007年12月20日基本为缓慢下降变化（图2d），21日后开始由降转升，2008年1月9日开始加速上升，直到6月17日，水温数据累积变化达0083℃，约为正常年变幅的6倍，在此时段发生了汶川80级大地震，异常变化形态非常显著，出现了明显的非对称“V”字型变化（何案华等，2012）。2008年6月18日后水温以指数形式回落，7月7日后开始以缓慢而稳定的速度呈上升变化，其后的2008年12月29日、2009年3月15日有两次比较明显的“V”字型下降变化，这可能是汶川地震的震后效应，但持续时间不长，变化幅度不大。临夏水温自2009年3月15日出现最后一个“V”字型变化后，水温以稳定的速率呈直线斜形上升变化。

32直线回归分析

2007年临夏井b值为负值（表1），2008年转为正值，且b值较大，为0002，在此期间发生汶川80级地震。自2009年至今，b值变化不大。

33斜率K值法

临夏水温自2007年6月15日至2008年6月17日，K值基本为加速上升变化，2008年6月18日后急速回落，下降到一个低值后恢复到正常变化状态，在此时段发生了汶川80级大地震，K值异常变化形态非常显著（图2e）。2008年7月7日后，K值基本稳定在0000 2附近，没有趋势上升和下降变化形态，这种变化持续至今。

34从属函数法

临夏井μ值自2007年6月起基本上呈现上升型高值异常变化，汶川80级地震后不久，于2008年6月21日达到最高值0788，然后从高值快速回落到背景值0（图2f），此后至2009年3月22日期间还出现了3次μ值高值变化。青海玉树地震前没有出现μ值高值变化，2011年3日11日日本90级地震前11天，μ值出现了持续4天的高达032的变化。

4讨论与认识141异常特征的共同性及差异性两水温观测井在汶川地震前都有明显的前兆异常，水温测值震前都出现快速上升变化；但异常持续时间不同，临夏井震后1个多月基本恢复正常，而清水温泉井则持续到2008年底。

正常情况下，一个上升漂移型的井水温会以稳定的速率上升，因而斜率K值几乎保持不变，曲线形态呈现为直线形态；但当有干扰因素或有地震前兆异常出现时，K值就会有大的变化。两井在汶川地震前后斜率K值都处于高值状态，但清水温泉是由高值快速降低，而临夏井则是加速上升变化，震后K值变化的时间长度也不同，清水温泉井持续到2008年底，临夏井震后一个多月恢复正常。

从属函数是对测值斜率变化的另一种反映方式，虽然和K值变化有一定的相似性，但只要参数取值合理，异常变化则会更显著。两水温观测井μ值在汶川地震前出现了多次大于05的异常变化，异常变化都非常明显，震后都出现了μ值的高值变化；但清水温泉井震后μ值逐渐减小，而临夏井在震后μ值则呈现增大趋势，一个多月后达到最高值然后快速减小。

42认识

（1）对于长期漂移上升（下降）型且日动态变化很大的水温观测井，采用一定的数学方法进行异常的提取较合理，斜率K值法、从属函数法是两种非常好的异常提取方法。

（2）汶川80级地震前清水温泉井和临夏井都出现比较明显的异常变化。但两井的异常形态具有差异性，清水温泉井地震异常形态呈现为“对数变化形态”，而临夏井则呈现为非对称“V”字型变化形态。震后两口井都有继承性异常变化，但结束时间不同步。

（3）清水温泉井在青海玉树71级地震前出现了测值跳动、K值波动，μ值高值变化等异常变化形态，临夏井却没有出现异常。日本90级地震前两口井都出现幅度较小的异常变化，但震后形态不同，清水温泉井震后测值明显下降，K值降低，而临夏井在震后立即恢复正常。

43讨论

（1）清水温泉井和临夏井在汶川地震前都出现了明显的地震前兆异常。虽然两井都距汶川约500 km，但都在青藏块体的东北缘，因而印度洋板块对欧亚板块的挤压所引起的青藏地块应力场的增加，也同样在这两个地方能反映出来。在强震活动的高潮阶段，构造块体内各地震带（包括未参与强震活动的地震带）都有应力增强的过程（耿鲁明，石耀霖，1993）。如果把青藏地块应力场的增加看做是“因”，则汶川地震的发生，和清水、临夏应力场的增加则是3个不同形态的“果”。两口井在汶川地震后都继承性的延续了一定时段震前异常形态，这可能是汶川地震后累积的应力并未完全释放完毕，后续地应力释放引起了区域应力场的变化。

（2）清水温泉井和临夏井相距约290 km，而青海玉树71地震前只有清水温泉井出现了地震异常。这可能与两井所处构造部位不同有关，清水温泉井位于NNW、NWW向断层的交汇部位附近，处在鄂尔多斯地块、华南地块、青藏地块的交接部位。因为其构造部位的特殊性，按照“多点场应力集中模型”（马宗晋，1980），清水温泉井所处位置可能是应力集中点，对区域应力场的变化反应更灵敏。

（3）日本90级地震前后，清水温泉井原始曲线、从属函数异常都比较明显，震后不久水温数据恢复正常，变化速率也即斜率K值基本保持在一个基值附近。临夏井只有从属函数出现异常，但不显著。日本地震距这两口井3 000 km，这样的远兆异常在地下流体观测中不断被发现（车用太等，1999）。

参考文献：

车用太，鱼金子，刘五洲，等1999张北—尚义地震的地下流体异常场及其成因分析[J].地震学报，21（2）：194-201

车用太，鱼金子1997源兆、场兆、远兆及其在地震预报中的意义[J].地震，17（3）：283-289

付子忠1988地热动态观测与地热前兆[J].地壳构造与地壳应力，1（1）：1-5

耿鲁明，石耀霖1993地震孕育发生的场源关系初步研究[J].中国地震，9（4）：310-319

何案华，赵刚，刘成龙，等2012青海玉树与德令哈地热观测井在汶川与玉树地震前的异常特征[J].地球物理学报，55（4）：1 261-1 268

刘耀炜，孙小龙，王世芹，等2008井孔水温异常与2007年宁洱64级地震关系分析[J].地震研究，31（4）：347-353

马宗晋1980华北地壳的多（应力集中）点场与地震[J].地震地质，2（1）：39-47

杨兴悦，王燕，闫万生，等2006武山22号井水氡中期异常与地震关系的探讨[J].西北地震学报，28（4）：378-380

叶秀薇2004粤闽地区地下流体从属函数异常与地震关系的初步研究[J].防灾减灾工程学报，24（2）：195-201

赵刚，何案华，马文娟，等2011不同动态背景的地热对比观测研究[J].地震学报，33（1）：51-61

赵刚，何案华，王军2009地热正常动态特征的研究[J].地震，29（3）：109-116