0 引言

地热用于地震前兆观测在我国已有几十年的历史。许多专家对地热(水温)变化的特征和机理进行过专门的研究。陈沅俊等(1992)通过研究1989年大同-阳高地震前的水温异常认为，断层蠕动作为一种剪切摩擦滑动必然产生一定的热量，故有可能利用水温变化来研究断层在震前的加速蠕动并进行地震预测。姚宝树(2004)则通过对三马坊水温井的观测结果进行研究后认为，地震发生后地震波的波动改变了区域地应力场的作用，使岩体周围受力平衡再次发生改变，岩体围压发生变化，岩体内裂隙、节理、孔隙压等也随之发生变化，岩体裂隙、孔隙里的地下流体也发生变化，地下水受压上升，受张下降，流体的温度随之发生变化。杨竹转(2011)研究认为，同一口井的水位同震变化是水温同震变化的必要条件，水温同震响应总是出现在地震波到达和水位同震变化开始之后；且水温同震变化的幅度受震级、井震距、季节、地温梯度、探头放置位置等因素的影响。张彬等(2013)对2008年汶川M_S8.0地震、2010年玉树M_S7.1地震和2013年芦山M_S7.0、岷县M_S6.6地震引起的中国大陆井水温同震响应现象进行了研究，认为水温同震变化的幅度不仅受动力加载作用强弱(距离)的影响，而且很可能与区域构造活动状态有关，应力状态的变化会直接影响到含水层孔隙度等，因而直接导致地下水动力特征的变化。

“十五”数字地震观测网络项目以来，青海省先后建成了地热观测井孔14个，其中，水温观测井孔10个，另外4个井孔没有深层地下水，为地温观测。2015年4月25日尼泊尔8.1级地震时，大部分井孔的水温观测数据没有同震变化反应，只有都兰台深井地温(简称都兰地温，下同)观测数据出现了明显的同震变化。本文即对都兰地温观测数据的同震变化情况进行了总结，分析了同震变化特征，并探讨了地温同震变化发生的机理。

1 基本情况 1.1 井孔

都兰地温观测井位于柴达木盆地南部，可可西里-巴颜喀拉、柴达木2大地震带之间(图 1)。井孔深105.30m，其中，0.50m以上为粉土，含砂砾，0.50~＜44.80m为坡积角砾，44.80~＜105.30m为坡积角砾岩。井中没有地下水，为干井。温度探头置深为100m，以观测井底地层的温度变化(图 2)。

1.2 仪器

都兰地温观测井安装的仪器为中国地震局地壳应力研究所研制的SZW-1A数字式温度计，安装时间为2007年6月20日。仪器自从安装以来观测数据一直不正常，背景噪声非常大，数据经常出现大幅度的台阶和突跳，动态极不稳定。2013年5月仪器出现故障，经维修后，2013年6月1日起恢复正常。

1.3 长期观测动态变化特征

2013年6月1日仪器恢复正常工作后，都兰地温观测曲线的背景噪声明显变小，动态变化也相对稳定，呈现出缓慢上升态势，并在缓慢上升的过程中叠加了一些小幅度的起伏和毛刺。背景干扰的幅度由约0.006℃下降到0.0001℃，日变幅约为0.0003℃，年变幅约为0.021℃(图 3)。

尽管2013年6月1日以后都兰地温曲线较以往变化平稳，但从日变化形态来看，仍存在一些短时间的小幅度变化，变化形态以向下—回升居多，但也有少量上升—下降的变化(图 4)。

2 同震变化

2013年6月1日都兰地温仪器正常观测以来，尽管存在许多不明原因的短期变化，但仍在震后多次观测到了明显的同震变化。2013年6月1日~2015年5月31日，全球共发生M_S≥7.0地震48次，我国共发生M_S≥5.0地震85次。对这132次地震(其中1次地震在2个目录中均有)发生前后都兰地温观测数据进行分析发现，有8次地震记录到了明显的同震变化过程(图 5、6)。这8次地震的参数和分布情况见表 1。

由表 1和图 5、6可见，都兰地温同震响应的地震都集中在中国大陆及周边，距离最远的为2015年5月30日日本小笠原群岛8.0级地震(Δ=4079km)。出现同震响应的最小震级地震为2014年10月2日乌兰5.3级地震(Δ=31km)，最大的为2015年4月25日尼泊尔8.2级地震(Δ=1553km)。这些地震均为相同震中距范围内震级最大或较大的地震。另外，2013年8月1日18：24墨西哥海岸远海7.5级(Δ=13852km)、2013年11月17日17：04斯科舍海7.8级(Δ=16352km)2次地震后地温数据也有变化，但是墨西哥7.5级地震在都兰地震台及周边台站均未记录到振动波形，表明该地震的地震波到达都兰台时已经完全衰减，无法引起同震变化，故同时段的地温变化应该不是该地震的同震变化；斯科舍海7.8级地震后的地温变化比地震的面波晚到将近半个小时，而其他各次地震的同震变化开始时间均在面波(乌兰地震为S波)到达后1~10min之内，故此次变化应该也不是斯科舍海7.8级地震的同震变化。

3 典型同震变化及特征 3.1 典型同震变化 3.1.1 甘肃岷县6.7级地震

2013年7月22日07：52，都兰地温数据出现明显下降，08：12数据下降减缓，至17：04数据下降到最低值，累计下降幅度为0.0008℃。之后，数据变化转平，并在低值持续约5hr后重新开始缓慢上升。测震记录波形显示，都兰台岷县6.7级地震P波到时07：47，S波到时07：48，面波到时07：48。都兰地温值的下降，是在面波到达4min后开始的，其变化应是面波到达后所引起的变化(图 7(a))。

3.1.2 新疆于田7.3级地震

2014年2月12日17：31，都兰地温观测数据出现下降，至17：34累计下降0.0003℃，17：34~17：47数据持平，到17：48数据小幅度回升并转平，18：37起再次缓慢小幅度下降至19：18转平，之后恢复数据缓慢上升的正常动态变化。测震记录波形显示，都兰台于田7.3级地震P波到时17：22，S波到时17：25，面波到时17：27。都兰地温值的下降，是在面波到达4min后开始的，其变化应是面波到达后所引起的变化(图 7(b))。

3.1.3 尼泊尔8.2级地震

2015年4月25日14：20，都兰地温观测值出现快速下降，至14：28下降到最低值，下降幅度达0.0028℃。14：31数值开始回升，到16：13基本恢复到下降前的水平。整个过程持续114min。测震记录波形显示，都兰台尼泊尔8.1级地震P波到时14：14，S波到时14：17，面波到时14：19。都兰地温值的下降，是在面波到达1min后开始的，其变化应是面波到达后所引起的变化(图 7(c))。

3.1.4 日本小笠原群岛8.0级地震

2015年5月30日19：41，都兰地温观测值出现明显下降，至19：52下降幅度达0.0003℃，之后缓慢恢复，20：16恢复到下降前的水平。整个过程持续35min。测震记录波形显示，都兰台日本8.0级地震P波到时19：29，S波到时19：34，面波到时19：37，最大面波在19：43~19：47出现。都兰地温值的下降，是在面波到达后开始的，其变化应是面波到达4min后所引起的变化(图 7(d))。

3.2 同震变化特征

各次同震变化的具体特征见表 2。由表 2可见，都兰地温的同震变化具有以下特征：① 变化形态主要有2种：下降—回升(6次)、下降—转平(2次)。都兰地温同震变化第1阶段的变化形态均为下降，之后的变化形态则有低值回升和低值转平2种。② 尼泊尔8.2级地震的同震变化的幅度最大，达0.0028℃，其他几次地震的同震变化幅度均小于0.001℃，其中，乌兰5.3级地震的同震变化仅在其背景变化基础上下降了0.0001℃，这应与震级较小有关系。③ 除乌兰5.3级地震外，其他地震几乎都是在面波到达后才出现同震变化反应。其原因应是面波振幅远大于体波，比体波携带了更多的能量，因而其更易于引起同震变化。乌兰5.3级地震震中距仅有31km，未记录到面波，尽管震级不大，但其体波也对都兰地温造成了小幅度的影响。④ 各次同震变化与该次地震面波(乌兰地震为S波)到达台站的时间差为1~10min。其中，4个地震的同震变化发生在面波到达后4min，2个地震分别为3、5min，尼泊尔8.2级、巴基斯坦7.9级地震则分别为1、10min。⑤ 下降—回升形态变化的持续时间大多为几十分钟，仅有尼泊尔8.2级地震因震级较大，变化幅度大，持续时间近2hr。下降—转平形态变化的持续时间均较长，为5~7hr。

4 同震变化幅度特征分析 4.1 同震变化幅度与地震波振幅的关系

表 3统计了都兰台测震仪记录到的各次地震的最大振幅，并分别计算了最大振幅与同震变化幅度之间的拟合相关系数。由表 3可见，各分量振幅与同震变化幅度的拟合相关系数均大于0.800，表明同震变化幅度与地震波最大振幅之间有较好的线性关系，其中，水平分量合成振幅的拟合情况最好，拟合系数为0.924，三分量合成和NS分量的相似程度也较好，拟合系数都大于0.88。

上述分析均间接证明，都兰地温一系列的以下降为主要特征的异常变化，确实是上述地震的地震波振动所引起的同震变化。同震变化的幅度，更多地受水平方向合成振动的影响。

4.2 同震变化幅度与震级的关系

由同震变化幅度与震级间的关系可见(表 3)，同震变化的幅度对于同一震中距而言，随着震级的增大而增大；震级相近时，则随着震中距的增大而减小。如2015年4~5月的3次尼泊尔地震，震中距相同，同震变化的幅度随着震级的变化而不同。尼泊尔8.2级、日本小笠原群岛8.0级地震，震级相差虽然不大，但前者比后者震中距要小，其同震变化的幅度也比后者要大。巴基斯坦7.9级、日本小笠原群岛8.0级地震，前者震级和震中距均略小于后者，二者的同震变化幅度基本相同。

由此可见，同震变化幅度与震级间的线性关系并不明显。由于单台计算地震震级的大小不仅取决于单台记录到的地震波的振幅大小，而且取决于震中距的远近，因此，二者的拟合相关系数不大，仅为0.437。

5 变化形态及原因分析

于田7.3级、岷县6.7级2次地震与其他地震的都兰地温同震变化形态有所不同，均为下降—转平形态(图 7)。这表明在地温同震下降后，探头处温度未恢复，而是达到了新的热力学平衡。

从空间分布来看，于田7.3级地震发生在巴颜喀拉地块西边界附近的硝尔库勒盆地南缘，为巴颜喀拉地块向东挤出的构造响应和应变调整所致(吴传勇等，2014)；岷县6.7级地震震中位于巴颜喀拉块体东北缘附近的临潭-宕昌断裂带中东段，青藏高原东北部东昆仑断裂的向北挤压和向东的运动是该地区构造应力集中的主要动力(郑文俊等，2013)。因此，这2次地震都是青藏高原北部块体向北和向东的挤压作用所导致的。都兰井位于青藏高原北部的东昆仑-柴达木地块内，地块南缘与巴颜喀拉块体相邻，与这2次地震震中一样都位于青藏高原北部的重要活动断裂——库玛断裂带北侧附近。因而，这2次地震所产生的应力变化，都对都兰井造成了直接的影响，导致其地温同震变化后无法恢复原状，发生了类似“塑性”的变化。

其他几次引发同震变化的地震，如乌兰5.3级地震震中距虽近但震级较小，同震的影响也很小；其他地震虽然震级都大于7级，但震中距都较大，发震构造也与都兰井没有直接的联系，其同震变化均为下降，并在之后的较短时间内恢复原状(图 7)，这种变化类似“弹性”的变化。

对比以上2种变化形态可见，当地震发震构造与都兰井之间有相对直接的构造联系且震级较大时，就会导致都兰地温发生“塑性”变化；而当地震的发震构造与都兰井没有相对直接的构造联系或震级较小时，就会导致都兰地温发生“弹性”变化。

6 结论和讨论 6.1 结论

(1) 在井孔周边一定范围内且震级较大的地震发生后，都兰深井地温观测数据出现同震变化，变化形态有下降-回升和下降-转平2种。

(2) 都兰深井地温同震变化均在都兰台所记录到的最大振幅地震波(近距离地震为S波，远距离地震为面波)到达1min后出现。

(3) 都兰地温同震变化幅度与最大地震波幅度有关，相关性大于80%，但与震级的线性相关性不明显。

6.2 地温同震变化机理探讨

水温同震变化是由于地震发生后地震波的波动改变了区域地应力场的作用，使岩体围压发生变化，导致岩体内裂隙、节理、孔隙压也随之发生变化，从而使得岩体裂隙、孔隙里水的流动和汇集方式发生变化所致。但是对于都兰井这种无水的地温观测井，其变化的机理尚需要进一步的探讨。

对于一般的弹性系统的热力学状态，可用应力σ、应变ε和温度T等3个参量描述，写成微分形式(谢锐生，1980)为dσ=Edε+βdT，其中，E为杨氏模量，β为热应力系数。此状态方程经常用来研究热胀冷缩引起的应力问题，但涉及应力引起温度变化的研究则很少。对于等温情况(dT=0)，上式变为胡克定律dσ=Edε，也即为弹性力学侧重研究的内容。在绝热情况下，固体弹性变形的状态方程不再是胡克定律，受热膨胀的材料也会受压升温(对于一维情况：ΔT=aTΔσ，其中，T为初始温度，a为与材料性质有关的热物理常数)。实际上，对于缓慢过程，可近似为等温过程。对于快速变化(如地震)，更接近绝热过程，温度与应力变化相关联, 除了可以通过变形测量获得应力外，还可以通过温度测量来获得应力变化状态(陈顺云等，2013)。尤其是对于弹性模量大的介质，即使应力变化很大，变形也可能较小，应力的变化则可通过温度的变化来体现。

由此可见，都兰深井地温观测数据的同震变化，可能正是由于地震波所携带的能量使得井孔介质的应力状态发生改变，而这种应力状态的改变导致井孔周边的热力场分布发生了明显的变化。

都兰井的井孔基质为坡积角砾岩，质地较坚硬，基岩疏松。这样的介质并不属于上文中提到的弹性模量大的物质。这种介质是否能够更好地感应应力的变化，尚未见相关报道。但从目前的情况来看，青海省其他地热观测井孔都位于基岩之中，却没有明显的同震反应。而都兰井这一介质相对疏松的井孔，却观测到了清晰的同震变化。这与其构造位置有关，还是与其井孔介质有关，需要在以后的研究中逐步深入探讨。

6.3 其他需要进一步研究的问题

都兰地温的动态变化背景存在许多不明原因的短时间上升或下降变化，其中，多数下降变化与同震变化的形态和幅度都非常类似(图 4)。尽管通过上述分析后认为，都兰地温在震后发生的这种变化与地震之间确实存在明显的联系，但仅从变化形态和幅度上进行分析，大部分无法明确地直接判断是否为同震变化，必须结合地震目录及该台地震记录波形才能进行准确判断。因此，对于都兰地温这种不明原因的短期变化，建议开展相关研究，尽可能给出其原因并予以排除。