边界元法在地电预报地震中的研究与应用-地震工程学报.pdfVIP

第# 卷! 第 期 西! 北! 地! 震! 学! 报 BCD$ #! 3C$ ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! #%%/ 年’ 月 3456789:6953 :9;:4=4;?,= @4A53,= EFGH ，#%%/ !!! ! 综! 述 ! !!! 边界元法在地电预报地震中的研究与应用 # 马钦忠 ，钱家栋 （$ 上海市地震局，上海! #%%%# ；# $ 中国地震局分析预报中心，北京! %%%’ ） 在以往地电预报地震研究及电法勘探研究中，二维、三维点电源层状横向非均匀介质条件 摘要： 下的数值模拟方法存在着不少困难问题。利用边界单元法可使得该类二、三维问题中难处理的 问题变得较为简单易行。本文综合分析叙述了边界元法在这类地电断面结构下的理论研究及其 在地电预报地震中的应用情况。探讨了地电学在地震预报研究中的理论及实际应用前景。 关键词：地电法；边界元法；非均匀介质；地震预报 中图分类号：(’)$ *) + #! ! 文献标识码：,! ! 文章编号：%%% - %.// （#%%/ ）% - %%. - % %! 引言 地电阻率法研究的是地震孕育过程中地球介质电学性质变化，探索这些变化与地震发生的时间、地 点、强度之间的关系，找出其时空演化特征与孕震过程的联系。地电阻率法要求通过地表视电阻率的测量 来了解地下介质内不同部位电阻率随时间变化的情况。在实际工作中，由于地下介质常常是非均匀的，加 之孕震过程中地球介质电学性质变化信号与地下电性结构非均匀性变化等其它非震变化信号往往叠加在 一起，使得识别震前地电阻率异常前兆信号的难度大大增加了。所以研究非均匀介质结构对视电阻率的 影响问题是非常重要的。一般而言讨论非均匀介质是较为困难的，理论上能够处理的只有一些特殊的分 区均匀介质，如水平层状介质［］ ［# ］ ［’ ］ 、三层倾斜层介质 以及半无限空间介质中包含另一种球状介质 等比 较规则和简单的情形。因为对于复杂的地电断面问题，例如不规则地形、非均匀介质等问题，很难给出精 确的解析解，所以数值计算方法便成为非常行之有效的手段。目前处理地电断面问题的数值模拟方法主 要有有限元法、有限差分法、边界元法等。由于有限元法的特点是不论对什么问题均在各方向离散，并一 律采用分片低阶多项式插值来逼近各类问题的解函数，这样固然具有通用性强的一面，但又不可避免地带 来自由度多、工作量大的不利一面。因而在这些数值模拟方法中，对地电阻率法中无限区域边值问题而 言，边界单元法已变为良好的数值计算技术。边界单元法应用的发展是由于其所具有的某些优点，如边值 ［/ ］ 问题维数的减少，模拟区域可扩展到无限域并具有较好的数值精度 。我们利用边界单元法这种数值计 算方法对层状介质横向电性不均匀的地电断面模型进行了研究［) 0 . ］。结果表明，对于点电源层状横向非 均匀介质的地电断面而言，边界元法是一种非常有效的数值计算手段。 边界元法在地电阻率法中的研究 #% 世纪% 年代以来发展起来的边界元法，特点是只在求解域的边界上剖分单元，而不是象有限元法 和有限差分法那样在求解区域内剖分网格。边界元法具有以下优点： （ ）能使问题的维数降低一维，并且单元节点少，输入数据简单，占用内存少，计算速度快； （# ）不仅适合有界区域问题，而且更适合研究无界区域问题，并给有场源和集中载荷作用等问题的处 理带来了极大的方便； 收稿日期：#%%’1%#1* 基金项目：中国地震局地震科学联合基金项目（%#%.% ） 作者简介：马钦忠（2)2 - ），男（汉族），河南洛阳人，副研究员，现主要从事地震预报工作。 第$ 期, , ,