从水位观测数据中排除非构造因素影响的状态空间分析方法.pdfVIP

第 27卷第 2期 华 北 地 震 科 学 Vo1．27，No．2 2009年 6月 NOR1’’H CHINAEARI、HQUAKESCIENCES JHn．，2009 文章编号 ：1003—1375(2009)02—0017～05 从水位观测数据中排除非构造因素影响的 状态空间分析方法 邓 亮 ，刘春平2，万 飞 ，廖 欣 (1．湖南师范大学资源与环境科学学院，湖南 长沙 410081； 2．防灾科技学院，河北 三河 101601) 摘要 ：为 了更有效地利用地下水位观测数据开展地震预测 ，需要排除其它因素对地下水位 的干扰 。降雨 、气压和潮 汐是影响地下水位的3个重要 的非构造应力因素，对此建立状态空间模型以提取与地震有关的水文异常。这个模 型把观测到的地下水位时间序列分解为气压响应、潮汐响应、降雨响应 、噪声因子和水位残差 5个因子 ，运用状态 空间分析方法 ，从地下水位 的实际观测数据 中分解 出气压 、潮汐和降雨等主要 的非构造 因素对地下水位 的影响，从 而得到构造应力影响下的地下水位数据 ，以便于提取与地震活动有关的异常信息。 关键词 ：地震 ；地下水位 ；状态 空间模 型 ；异常 中图分类号 ：P315．72 文献标识码 ：A 变动 ，已由许多研究证实 t2J。然而 ，影响水位变化 0 引言 的因素很多 ，有些是人为的，有些因素却是来 自自然 地震发生前常有一些异 常现象发生 ，如地下水 界的非构造 因素。来 自自然界的非构造因素影响则 位 、地电、地温、地磁等，其数据容易依靠仪器来观测 须寻求适当方法滤除。过去有多位学者 曾由不同的 记录 ，可通过分析这些数据来发现地震前兆 ，以达到 角度及方法探讨这些非构造 因素所造成 的影响，如 预测地震的目的。地下水位在地震预报的探索中引 最小二乘回归分析法_3J、频变转移函数法 J、BAY— 起了广泛的重视 ，因为地下水位在一定程度上反映 TAP—G方法 (BayesianTidalAnalysisProgram InA 了岩层的应力应变，裂 隙的发育与 串通 ，断层 的闭 GroupingMethod)[53和非 贝叶斯 回归分析法 6【』等。 锁 、蠕动和错动等。地下水存在普遍性和流动性，是 上述方法 中最小二乘 回归分析法简单实用 ，但不能 可以将深部信息直接带到地面的载体，但地下水 的 完全描述气压响应 ；频变转移函数法和 BAYTAP— 水位常因人为抽水或注水 、降雨等影响造成不规则 G方法可估计气压和潮汐 的响应系数 ，但不能估计 的变化 。对封闭型水源的深井而言，虽然没有人为 降雨的响应；非贝叶斯回归分析法能把降雨响应从 抽水或注水 的影响，但仍有大气压力 、地球潮汐及降 地下水位 中分离出来 ，但水位残差 (含构造应力的响 雨等非构造因素 的影响效应 ，这些影响效应 常混淆 应)包含在观测噪声里。考虑到在线弹性含水层 中 来 自地壳形变的水位异常变化 ，以致减弱地下水反 地下水位与各类应力之间为线性关 系，且构造应力 映地震发生前地壳应力改变的监测能力。本文即针 的非周期性和非构造应力 的周期性和周期差异等 ， 对此现象，建立状态空问模型 以除去这些非构造 因 本文将参考前人之研究 ，建立状态空间模型，应用时 子影响地下水位变化的效应 ，提供进一步探讨除去 间序列分析方法把观测到的地下水位分解为3个非 这些影响效应后的水位变化是否与地震有关的基本