地磁方法在地震预报中的应用.doc

 地 震 地 磁 观 测 与 研 究 SE ISM OL O G ICA L A N D GEOM A GN E T IC OB SERV A T ION A N D R E SEA R CH 第 20 卷 第 6 期 1999 年 12 月 V o l 20 N o 6 D ec 1999 地磁方法在地震预报中的应用 林云芳1) 曾小苹1) 续春荣2) 明1) 琪1) 赵 李 1) 中国北京 100081 中国地震局地球物理研究所 2) 中国北京 100085 中国地震局地壳应力研究所 摘要 1989 年以来, 我们陆续提出转换函数法, 空间相关和加权差分法, 以及加卸载响应比法 等 3 种地磁方法, 将其用于 200 多个震例分析, 其中有 106 个地震作了正式的年度预报, 在 97 个次危险区预报中, 有 15 个地震的时间、 地点和震级三要素的预报较为成功。 关键词 地磁方法; 地震预报; 预报评估 引言 地磁现象可反映上自日地空间, 下至大气层和固体地球中发生的与电磁有联系的各种 物理过程。 对主磁场及其变化的研究, 是了解地球深部地幔和地核内发生的物理过程的重 要途径。 利用变化磁场在地球内部的电磁感应现象, 是研究地球内部电性结构及其变化的 一种好方法。 地震磁现象的观测研究, 迄今已有近百年的历史。本世纪 50 年代末, 精度高, 稳定性 好的质子旋进式磁力仪 (又简称“核旋仪”) 的问世和广泛使用, 改变了对地震磁现象的观测 状况, 获得了较为准确的地震磁现象的观测数据和较为可信的震例结果。 近年来, 震磁关系的研究取得了一些进展。1987 年以来, 我们通过对中国地区地球磁 场局部异常现象的研究, 一边探索震磁关系, 一边进行实验性的地震预报。1989 年开始, 我 们陆续选用转换函数法, 空间相关和加权差分法, 以及加卸载响应比法等 3 种地磁方法, 分 析了 1966 年以来的 200 多个震例。 对其中 106 个地震进行了年度预报, 在 97 个次危险区 趋势预报中, 有 15 个地震的时空强三要素与实际发生的情况比较接近。 1 地磁转换函数法 1955 年 R ik itak e Yo k u yam a (1955) , 1959 年 P a rk in so n 证实了地磁短周期变化的垂 直分量 ?Z , 水平分量 ?H 和磁偏角 ?D 具有以下稳定线性关系 ?Z = A ?H + B ?D (1) 其中 A 和 B 就是转换函数。 中国地震局地球物理研究所论著编号 99A C 2060 利用天然磁场作为入射信号, 其特点是频带宽, 分布广, 地面各测点记录到的天然磁场信号的分析得到的转换函数 A 近似平面波, 有穿透力。 对 和 B , 是地点 ( Κ, Υ) , 深度 (h ) , 频率 (Ξ) 或周期 (T ) 和时间 ( t) 的函数, 即 A = A (Κ, Υ, h , Ξ, t) B = B (Κ, Υ, h , Ξ, t) (2) 由电磁波理论的趋肤效应可知, 某一测点处 (Κ0 , Υ0 ) 某一频率 Ξj 的磁场 F (Ξj ) 的穿透深 度为 h j = (2?Ξj ΛΡ) 1?2 (3) 此处 Λ 和 Ρ 分别为介质的磁导率和电导率。不同频率的电磁波对地球的穿透深度不同, 因此, 可以利用各种不同频率的变化场的电磁感应获得地下不同深度的电导率的信息。 分 析和研究 A 和 B 及其参量随时间 t 的变化, 则可了解和监视地震活动区地下电性结构和环 境的变化, 以期在地震发生前捕捉到这种变化, 进行地震危险区的预报。 1987 年以来, 通过分析 100 多个震例后, 我们认为复转换函数的模|A | 和|B | , 解矩阵 方程 (1) 的总方差 Ρz , 以及 P a rk in so n 矢量的磁方位角 Αp 及其变化等几个参量的变化异常与 测点附近发生的地震有关。 这些参量的定义为 A = A r + iA i iB i B = B r + (4) 下角标 r 和 i 分别表示实部和虚部。 |A | = (A 2 + A 2 ) 1?2 r i |B | = (B 2 + B 2 ) 1?2 (5) (6) r i Αp = a rc tan (B r ?A r ) n 1?2 ∑ m = 1 Ρz j = (7) n ?Z ∑ (?Z jm ?n ) (8) m = 1 其中 m 为第 j 周期的事件个数。 计算了全国