甘肃清水温泉井与临夏井水温在汶川地震前后的异常变化特征
甘肃清水温泉井与临夏井水温在汶川地震前后的异常变化特征
作者：蒲小武　武银　狄国荣　梅东林　陈彦平　王军燕　叶媛媛
来源：《地震研究》2013年第03期
        摘要：分析了甘肃省清水温泉与临夏水温观测井自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异常提取方法
5分钟前兰州地震了
        （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值

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。