基干年降水干扰排除云南地区地下水位群体异常研究.docVIP

基于年降水干扰排除的云南地区地下水位群体异常研究 　　摘要：利用云南地区10年以上的20口地下水位观测井孔资料，选出其中受降水影响较大的18口井，分别采用峰―谷幅值法和谷―谷幅值法对其观测资料进行年降水干扰排除后，基于这20口井的资料，提取全省和红河断裂东、西部地区的地下水位群体异常。结果表明：云南水位观测井的群体异常与云南地区6级以上地震有较好的相关性，其中谷―谷幅值方法提取的异常、预测效果较峰―谷幅值方法好，可能表明降水影响的延迟过程对这些井更为重要。分区预测的预报效能不如全省，可能说明井水位异常增多是区域构造活动和应力水平增加的结果，因此将用较小区域的异常对应较小区域的地震，虚报、漏报就会增多，从而降低预报效能。在异常增多区域构造活动增强的背景下，有利于6级以上地震发生，因此将群体异常指标用于云南的年度地震活动水平预测，具有一定的效能 关键词：地下水位；降水影响；群体异常；区域构造活动 中图分类号：P315.723 文献标识码：A 文章编号：1000-0666（2016）04-0545-08 0 前言 地下水作为地壳结构中最活跃的组分之一，会因受到地壳应力状态改变而引起含水层中含水量变化的影响，能够比较灵敏地反映出地壳应力应变的信息（万迪坤等，1993）。震例研究结果表明，在许多中强地震前均记录到了地下水位异常变化，尤其在云南地区，地下水位异常作为中强震预测的重要依据，目前已积累了很多丰富的震例经验（万登堡，1992；张体移等，2012；付虹等，2014） 地下水位作为流体学科重要的物理观测量，是一个包含大气降水、气压、固体潮、地应力场等多种影响因素的复合参数（车用太，鱼金子，2006），其中降水是地下水位变化的主要影响因素，它使得多数井孔所观测记录的水位呈现雨季上升、旱季下降且比较容易识别的年变规律，但当年降水量随着旱、涝年份的不均匀而出现明显差别时，就将成为地下水位异常识别的重要干扰因素 云南大部分水位观测井的井深均小于300 m，且受降水干扰极为突出，在水位的资料分析中排除降水对地下水位观测带来的影响，对提高水位资料的使用率是非常必要的。一些学者对云南水位的降水干扰排除已做过一些研究（白宝荣，付虹，2006；毕青等，2012），得到了一些有益的结果，本文在前人研究的基础上，增加了水位年变幅度与降水量的相关性研究以及异常机理的讨论 1 基于降水干扰排除的地下水位异常提取 1.1 资料选取 云南省目前共有水位观测井孔63口，其中有10年以上观测记录的有29口，排除观测质量较差和基本无震例对应的9口井孔，笔者将观测连续稳定、映震效果较好的20口井的井水位观测资料作为研究对象。其中受降水影响显著（即观测数据年变形态清晰）的井孔有18口，受降水影响不明显的井孔2口，具体井孔分布见图1。雨量数据收集来自于云南省气象局，井孔水位观测起始时间不一致，所以群体异常结果采用异常井孔所占同期全部井点比例表示，具体始测时间见表1 1.2 不同的水位群体异常提取方法 大气降水对地下水位的影响具有滞后效应，因此笔者主要考虑年降水量对水位年变幅度的总体影响，分别利用峰―谷幅值法和谷―谷幅值法对受降水影响明显的18口井排除降水对地下水位观测的影响后进行异常识别和提取，对其余2口不受降水影响的水位观测井孔采取原始曲线人工识别异常，主要识别大幅上升或下降为主的异常形态 1.2.1 峰―谷幅值法 用地下水位同一年的峰值与谷值之差（即年变幅度）和年降水量进行相关分析，将其称为峰―谷幅值法。该方法以每年地下水位观测的年变幅度为纵坐标，降水量为横坐标，横纵坐标的取值存在明显的线性关系，即降水与年变化幅度相关，该趋势变化可由一条直线进行拟合。如果当年变幅超出一倍均方差，出现明显偏离正常范围的情况时，则表示该年份水位处于异常状态（图2）。表2给出利用峰―谷幅值方法排除降水干扰以后的地下水位异常结果 对每年排除降水干扰后的异常井孔比例进行统计，根据震例对应情况，以异常比例大于20%作为异常指标，得到自1998年以来，利用峰―谷幅值法排除了降水影响之后，所提取出的水位群体异常与云南省内未来一年6级以上地震的对应情况，如图3所示。根据历史震例统计，1998年以来，云南地区共发生6级以上地震9次，根据年份可划分为7组，从图3可见，在2009年姚安6.0级和2014年盈江6.1级、鲁甸6.5级、景谷6.6级2组地震前出现了较为显著的水位群体异常，其余5组6级以上地震前均未出现相对应的异常现象，对应率为29% 1.2.2 谷―谷幅值法 对地下水位谷值变化逐年差值与年降水量进行相关分析，将其称为谷―谷幅值法。该方法用每年地下水位观测的最低值（雨季来临前），以当年与前一年的谷值差值为纵坐标，前一年的