1. 引言
地下水位(流量)和水温观测一直作为地震前兆观测的重要测项,在我国已有几十年的历史,取得了许多宝贵的观测资料和一些震例资料。数字化下水位和水温观测在我国已有十多年的历史,“九五”和“十五”期间中国地震局对水位和水温观测台网进行了数字化改造,甘肃地区也建立起比较完整的地下流体观测网。
国内外地震学者曾以大量的观测事实论证了水位动态不仅可以直接反映受力状态变化引起的孔隙压力的变化,而且还可以反映出含水层地下水径流的速度与流量等的变化;不仅水位对含水层受力有所反映,对井水微温度也会有灵敏的响应。井–含水层系统所受的动力加载作用方式有多种,地震波作用就是重要的一种,地下流体观测中水位与水温对大震的同震响应是地震波作用于井–含水层系统最直接的体现 [1] [2] [3] 。近几年来,大震引起地震地下流体同震响应的现象受到不少学者的关注,刘耀炜等 [1] [4] 、孙小龙等 [2] 、张昱等 [5] 、刘成龙等 [6] 、周志华等 [7] 分别分析了不同地震水位与水温的同震响应特征。
据前人研究,地下水位的震时和震后效应,从形态上分为两类,一类是急速的阶跃变化,这主要是近震的井水位效应;另一类是水位的振荡变化,即为水震波效应,这种变化一般在震中距大于2000 km的井孔中比较常见 [4] - [9] 。水温同震响应主要是地震波作用引起的含水层介质形变,使其孔隙压力发生变化,这种孔隙压力的变化会导致水流速度及水体与围岩间的热量变化,而各水体间及水体与围岩的热量交换既有加强也有减弱,从而导致观测井内水体温度产生变化 [2] [9] [10] 。另外,井水温度的变化与井–含水层系统参数改变密切相关。
本文分析了2013年7月22日甘肃岷县–漳县6.6级地震前后甘肃地区的水位、流量及水温资料,结果发现在这次地震前两个测点的资料存在明显的震前异常 [8] ,部分测点也出现不同程度的同震响应。
2. 本次地震的基本参数及观测资料概况
2.1. 地震基本参数
本次地震的发震时刻为:2013年07月22日07时45分,震中位置位于甘肃省岷县漳县交界,地理坐标为:北纬34.5˚,东经104.2˚,震源深度20 km;发生在西秦岭断裂带南侧的临潭–宕昌断裂上,极震区烈度为VIII,VIII度区西北自岷县中寨镇,东南至岷县禾驮乡东南、东北自岷县禾驮乡东北、西南至岷县禾驮乡西南,长轴为40 km,短轴为21 km,面积706 km2。
该地震灾区主要涉及甘肃定西市岷县、漳县、临洮县、渭源县、陇西县,甘南州临潭县、卓尼县、迭部县、舟曲县,陇南市宕昌县、礼县,临夏州康乐县和天水市武山县等13个县区。灾区总面积约16,432 km2,灾区人口约2,359,126人、589,782户,其中,乡村失去住所人数为279,684人,69,921户,城镇失去住所人数为34,436人,8609户。
2.2. 观测资料概况
甘肃地震地下流体(水物理)观测至2013年7月22日地震前的资料包括水位、流量、水温3种观测手段,42个测项;其中流量3项,水位16项,水温23项;观测资料最早为1985年1月开始的人工流量观测、之后为“九五”、“十五”数字化改造、直至汶川地震灾后恢复重建,观测最短的为2011年11月至12月开始。水位均为静水位观测,流量观测为流满1升水需要的时间。表1为测点与震中的大概距离,图1为震中、观测点及同震测点位置。
![](Images/Table_Tmp.jpg)
Table 1. Epicenter distance of measuring points
表1. 测点与震中的距离
![](//html.hanspub.org/file/8-2410576x10_hanspub.png)
Figure 1. Distribution of water level (flow), water temperature measuring points and the co-seism points
图1. 水位(流量)、水温观测点及同震测点分布图
3. 6.6级地震震前异常分析讨论
对所有观测资料进行了分析,发现只有清水李沟流量和清水温泉水温出现比较明显的震前异常。
3.1. 清水李沟流量
清水李沟流量观测井位于甘肃省清水县城以西6 km,处于西秦岭北缘断裂的北侧。观测资料始于1985年1月,多年来观测资料整体呈下降趋势,观测以来,附近地区的多次中强震前出现比较明显的异常,而以往震例为在高值下降的基础上发生地震,只是以往异常幅度较小,而2003年甘肃民乐–山丹6.1级地震和2013年甘肃岷县–漳县6.6级地震异常幅度明显较大。
由于2008年汶川地震引起清水流量的同震效应幅度较大,我们截取2009年以来的资料进行分析 [4] ,并对观测资料进行了趋势拟合分析(图2),从图中看出,无论是原始观测曲线还是趋势分析后的曲线,2012年4月出现的高值异常均比较明显,持续时间为3个月,至2013年7月下旬结束(在这期间,进行过现场调查,未发现环境干扰和人为影响,异常可靠),总体上升幅度达到26秒(每天观测为流满
1 升
水需要的时间),以往该测点出现异常持续结束后,到发震最长时间为9个月,这次在异常结束一年后发生地震,属中长期异常,距离较近,为180 km,地震之后的上升为同震引起。
3.2. 清水温泉水温
清水温泉水温观测点位于甘肃省清水县城东北的汤峪河峡谷内,观测井深443.3 m,观测井水位、水温不受降水影响。该测点为“十五”数字化观测,观测资料始于2007年7月,从图中看出(图3),多年来比较稳定,整体呈上升趋势。2011年初观测资料出现过低值过程,但附近地区并未发生5.0级以上地震,岷县6.6级地震前11天观测资料出现下降,下降幅度为0.005℃,显示出震前短临异常(据调查,观测仪器正常,环境无干扰)。
4. 水位(流量)观测同震响应特征及分析
本次地震,所有水位(流量)观测资料中,只有清水李沟流量、礼县流量、平凉柳湖水位3个测项记录到比较明显的同震现象,另外5个测项记录到很微弱的同震现象,其它11个测项没有记录到同震响应。分析发现,这些测点与震中位置、方向、距离无明显规律。
![](//html.hanspub.org/file/8-2410576x11_hanspub.png)
Figure 2. Observation of water flow discharge at Qingshui
图2. 清水流量观测曲线
清水李沟流量2013年以来观测资料比较平稳,每天早上8点钟观测,地震为7时45分,所以,当天观测时就观测到明显的同震效应,上升幅度为15S/L,以往观测资料在正常时段(环境无干扰、无异常)不会出现这么大幅度的变化,持续时间为40天左右,之后略有恢复(图4)。
礼县流量观测井位于甘肃省礼县石桥镇,处于西秦岭北缘断裂的北侧。观测含水层的封闭性较好。以往观测资料日变幅较大,每天观测时间为早上8点钟,2013年观测资料总体较平稳,地震之前观测资料比较稳定,当天观测是在地震后,观测到31S/L的上升变化,持续时间较短,20天左右恢复到上升前的水平(图4)。
平凉柳湖井原为自流井,2007年断流,2008年汶川地震之后灾后重建,架设水位观测,现有观测资料始于2011年11月。该井水位观测资料虽然在趋势上受当地开采的影响,还是能记录到一定的潮汐,在以往大震时水位和水温也记录到同震响应现象。这次地震时记录到比较明显的同震现象,地震之后观测资料快速下降,3小时之后基本恢复原来的上升趋势(图4,其中图4上、中时间坐标为日值,图4下时间坐标为分钟值)。
地下水位(流量)的同震变化,反应了地壳形变和地面震动引起地下介质贮层变形、孔隙疏通、裂缝的清理、产生裂缝等变化。地震引起水位同震变化主要有振荡和阶变两种形态,也有振荡中出现上升或下降变化的,也上升或下降之后随之发生振荡的。振荡型变化是指在地震波作用下水位快速来回波动,地下水位出现类似地震波的高频振荡,地震波经过后水位很快平静下来。阶变型的同震变化是指在地震振动的作用下,地下水位阶变式的上升或下降变化,它反映了地下介质的孔隙、裂隙被疏通或地下水力学特征发生改变,可能是塑性变化的
![](//html.hanspub.org/file/8-2410576x12_hanspub.png)
Figure 3. Observation of hot spring water temperature at Qingshui
图3. 清水温泉水温观测曲线
![](//html.hanspub.org/file/8-2410576x13_hanspub.png)
Figure 4. Observation of flow discharge at Qingshui and Lixian as well as water level at Liuhu and Pingliang
图4. 清水和礼县流量与平凉柳湖水位同震图
结果。这种变化通常需要十几分钟至几个月的时间才能恢复,也有恢复不了,为台阶式永久性的变化 [7] 。
5. 水温同震响应及分析讨论
水温观测资料中,只有平凉柳湖水温和成县水温在这次地震中出现比较明显的同震响应,其余资料基本不明显。这两个测点与震中位置、方向、距离无明显规律
平凉柳湖水温观测时间为2011年11月,观测资料变幅较大,但地震前后资料比较平稳,地震后快速下降,8时11分左右下降至最低,下降幅度为0.025℃,之后缓慢恢复,两小时之内基本恢复正常(图5)。成县水温观测资料始于2007年6月,多年来观测资料比较稳定,整体呈上升趋势。地震之后缓慢下降,下午14时之后快速恢复,20时开始基本恢复至正常水平,下降幅度为0.011℃ (图5)。
水位同震变化已有比较多的研究,而水温的同震响应现象近年来才引起地震学家的关注,其机理的研究目前也处于探讨阶段,而且比较复杂。另外,前人研究的水温同震响应特征 [2] [5] [9] ,多数以水温阶变下降为主,也有波动或上升变化,这也说明水温的同震响应受多种因素的制约。
地震波作用引起的含水层介质形变可促使其孔隙压力发生变化,这种孔隙压力的变化会导致水流速度及水体与围岩间的热量变化。如果这种作用持续时间较长,则观测井内水体温度是一持续变化的过程;如果井–含水层系统在地震波作用下发生了塑性形变,那么井–含水层系统会在短时间内达到新的平衡,则井水温度会很快达到新的稳定状态。当然,井水温度的变化与井–含水层系统参数改变密切相关。地震发生时,地震波作用于含水层系统后,使其介质发生了变形所致,这种应力作用一方面激活了孔隙、裂隙中的充填物(如气体、滞水等),使得在空隙内运移的水体流动状态发生变化,导致各水体间及水体与围岩的热量交换加强或减弱,从而导致观测井内水体温度产生变化 [9] 。
6. 结论与讨论
这次6.6级地震的震中距离多数观测点在400 km以内,从以上的讨论知道,出现明显异常的测项只有两项,这说明地震前兆异常的复杂性,当然,也可能关系到多个方面的问题,比如,观测仪器的稳定性、观测环境的好坏、还有观测井点的条件等等,但是,只要某种观测手段在震前能观测到一些异常,我们可以与其它观测手段和学科的资料综合,也可能对地震预报有一定的借鉴意义。
该地震有10个测项记录到不同程度的同震响应,其中5个测项反映比较明显,另外5个测项记录到比较微
![](//html.hanspub.org/file/8-2410576x14_hanspub.png)
Figure 5. Measured water temperature at Liuhu, Pingliang and Cheng Xian
图5. 平凉柳湖水温与成县水温同震效应图
弱的同震响应,其余测项没有记录到,而这些井点中也有在其它大震时记录到过同震响应的例子,这也说明地下水观测的同震响应受发震地点、震级、震中位置、地震波到观测井的传播路径及发震构造等多种因素的影响。
由于水位动态的物理意义明显,含水层受压时水位上升,受拉张时下降,因此震后阶跃上升的水井水位可能包含有区域应力场的信息,水位的同震阶变可视为大震后区域应力调整的表现,水位阶跃上升集中区可能也是区域压应力相对集中区,对未来该区域发生较显著的地震具有空间上的指示意义。
基金项目
中国地震局地震预测研究所兰州科技创新基地基本科研业务费专项(2017IESLZ07)。