地震弱短临前兆信息提取方法研究.pdfVIP

笫19卷 第 2期 地 震 研 究 Vo1．19．No．2 l996年 4月 0URNAL OF SE[SMOLOGICAL RESEARCH A9J． 1996 一 ． 地震弱短临前兆信息的提取方法研究 ≥I (It肃省平琼地 医地震局 ，744000) (甘肃省平凉中心地震台，744000) 摘妻：地震短临预报是地震预报的重大难题，它有赖于好的前兆手段和数据处理方法， 本文提 出 “弱前兆信 鼠” 的观点，井研究了提取弱短临前兆信息的方法和定量指标。从1986 年 以来，对平凉地区专业台及地方台包括应力、形变 电阻率、地下水动态和水氡六种手段、 九个台站的多年观测资料选用不同方法进行处理，最后筛选 出非线性 回归方法，井对过去 已 发生的地jl蔓作丁 回顾枪验 ，定出了合理 的异常判定尺度 ，较好地预报 了距本Ⅸ较近 的几次地 震 ，使经验预报、定性预报向概率预报、定量预报转变，并且实现了计算机程序 化 、 实用 主嗣词竺里兰； 亘皿；—tI—gr ， 伯7f／良 狠 怯 1． 问题的提 出及主要思路 在前兆观测资料用于地震预报的数据处理方法中， 目前大多采用线性处理方法，如 实用化攻关系统中的程序。这种方法对于变化形态基本上是直线的资料，其处理效果是 好的，然而对于非线性形态，诸如弓形、波浪形变化的资料却效果很差，以标准差△ (△ 为线性 回归的标准差，垌以区别非线性 回归的标准差 a)的2倍作为异常判剐标准有时 也嫌小，甚至用到3A，以弓形年变资料为例，若以数月以上的资料做处理时，超出2 △的异常点一般 出现于年初或年终的峰 (谷)值区，而这些峰 (谷)值往往并非异常。 为此，寻找一种方法，从观测资料中发现前兆信息，作为一种较为实用的日常前兆资料 分析处理手段，则是我们所希望的。 在方法的选择上，经过多样本处理结果的比较，选用了非线性回归模型 ，事 实 证 明，这种方法对于前兆资料数据有较高的拟台精度，从而有利于从 中提取弱短临前兆信 息。所谓弱前兆信息，是指异常信息的绝对值较小，这是相对于线性处理方法中的异常 判别标准而言的，在非线性处理方法中，其剩余标准差a量值远小于同一资料作线性处 理时的△，对于弓形资料，一般 aO．5A，以全国资料评比三连冠的平凉柳湖泉水氡资 料为例，1989年~1993年的△和 a的处理结果见表 1。 收稿 日期 1994，07，】8· 万存绪，男，48 ，高级工程 师。 2托 万存 结等。地震弱黼临前兆信 息的提取方法研究 4F 由表 1可看出，其 a在0．28A～0．45A间变化，各年 2口均小于A。另外，我们所 谓的前兆信息，并不试 图对三要素都加 描述，与其对三要素都加 粗糙乃至形式的描 述，不如仅对时间逮一因素作出较为细致J}匀描述更具有实 啊价值，而震级及震中位置的 确定可依赖于另外的预报手段，我们认为，前兆手段最主要的功能是 为地震预报提供时 间指标，而前兆手段的优势也在于此。 2． 数学模型及处理方法 在观测实践中，各种观测手段所测得的观测值物理量，一般间接地反映了地壳直力 场及其它与地壳运动有关的场的变化。而欲建立 观测值为 目标函数的数学模型，其中 最关键的反映地壳变化的变量无从得到，例如，在多元线性回归中我们只能建立观测值 与可观测的各相关干扰因素如温度、气压等的模型，这样对排除干扰也有所裨益，但就 数学模型本身而言，事实上承认观测值仅 由干扰因素决定，而与地壳变化无关。 由于入 地无 门，我们无从观测到能够引起观测值变化的最本质的变量，在建立数学模型时我们 必须另辟蹊径。