应用从属函数方法研究安徽水氡的映震特征.doc

应用从属函数方法研究安徽水氡的映震特征 汪世仙作者简介：汪世仙（1982—），男，安徽肥西人，助理工程师，主要从事地震流体观测技术及科研预报工作 作者简介：汪世仙（1982—），男，安徽肥西人，助理工程师，主要从事地震流体观测技术及科研预报工作 基金项目：安徽省地震局科研合同制课题项目资助、中国地震局三结合课题项目资助 金艳2） 朱方保1） 谢庆1） 束克明1）姜薇薇1） 1）中国合肥231511 庐江地震台 2）中国合肥230031安徽省地震局 摘要：利用从属函数法分析安徽省地下流体的水氡观测资料，对其异常部分进行震例分析。并通过对安徽庐江台的水氡观测资料的震例分析，提出郯庐断裂带中南段的中短期异常特征，并进行预报尝试。 关键词：从属函数 水氡 映震特征 前兆 0 引言 从属函数法提取地震前兆异常信息是较为成熟的方法，该方法作为前兆综合研究方法之一在地震预报中取得一定的成效。随后不断有学者对从属函数方法判定中短期地震前兆异常进行尝试和改进。本文采取的水氡观测资料主要是庐江台、蚌埠台、临泉台、香泉台，各台资料均是在安徽省水氡观测网点具有资料连续率好，资料可靠性好等特点。尤其是庐江台的地下流体观测资料具有明显的同震效应、震后效应以及地震前兆异常，通过对庐江台的水氡资料的分析，对郯庐断裂带中南段的未来一年的震情趋势进行讨论。 1 井孔概况 安徽省水氡观测主要有庐江、蚌埠、临泉、香泉等台站。庐江汤池1号井处在郯庐断裂带南段西侧，附近有磨子潭-晓天断裂、南港-龙门冲断裂和北北东向的下扬子破碎带；蚌埠观测井处于郯庐断裂带西侧，近侧通过的断裂有北西向的涡河断裂和北北东的固镇-怀远断裂；临泉观测井处于临泉-刘府断裂带近侧；香泉观测井位于郯庐断裂带以东，下扬子破碎带西侧。 2 从属函数异常判定方法 以庐江台为例，水氡测值在数次地震前兆变化中具有异常时间较短的特点，因此根据实际情况对从属函数的计算方法进行了改良，使用5日均值10点平滑滑动滤波，同时根据震前数据阶梯上升的特点，设置其异常指标为斜率正值。 从属函数值计算：求水氡变化曲线上升段的从属函数值。拟合窗口长度=9。将值置于右端点。其计算式为： （1） 式中，为观测值变化斜率；为滑动平均值与时间的相关系数；为经验常数。在计算庐江汤池1号井水氡从属函数时，取经验常数=0.1。异常判定：若≥0.5, 则视为异常。（张双凤，2005） 从属函数检验各台观测数据 根据各台数据变化的特征，通过从属函数检验，设计适合各台观测资料的从属函数参数，在分析时进行了去阶变处理。临泉台出现从属函数异常有4次（图1），确认为映震变化的为2008年汶川地震、2011年日本地震，，其余两次异常经过查询地震目录无有感地震在异常区域内发生。在台站附近发生的1999年12月30日利辛Ms4.1级地震、2006年4月5日临泉Ms3.0级地震均无前兆异常变化。 蚌埠台数据2001年以前变化基本为平稳型，2001年后有年变化，在处理数据时，去除年变，再经过从属函数处理，较大的异常变化为2000年（图2），经检验为观测条件改变。 图1 临泉台异常分析 图2 蚌埠台异常分析 Fig1. abnormal analysis of Linquan station Fig2. abnormal analysis of bengbu station 图3 1993年-2004年香泉水氡 图4 2010年-2011年香泉水氡 Fig3.radon of Xiangquan from 1993 to 2004 Fig4.radon of Xiangquan from 2010 to 2011 其余异常区域幅度较小，对应时间段内无明显地震对应。2008年汶川地震、2011年日本地震在蚌埠台也未出现较为明显的数据变化。2006年7月26日的定远Ms4.2级地震发生时蚌埠台水氡停测7个月，无法进行异常判定。 由于香泉台多次更换观测井，2004-2009年观测数据连续性较差，因此在处理数据时采取分段处理，1993-2004年为1个时间段（图3）。检索地震目录，发现香泉台在1994-2004年期间出现的数据异常较好的对应台湾Ms6.6级以上地震7次，从属函数值超过0.5的地震3次，另有一次函数高值，映震效应有待进一步研究。 2010-2011年为一个时间段（图4）。在安庆地震前的2011年1月12日水氡数据开始加速上升，1月18日达到震前最高值40.7Bq/L，上升幅度为8.5%，