成层半空间深埋电极产生的电位分布.pdfVIP

第 25卷 第 3期(246—255) 中 国地 震 Vo1．25 No．3 2009年9月 EARTHQUAKERESEARCH INCHINA Sep．2009 聂永安 、姚兰予，2009，成层半空间深埋电极产生 的电位分布 ，中国地震 ，25(3)，246—255。 成层半空间深埋电极产生的电位分布 聂永安 姚兰予 天津市地震局 ，天津市河西区友谊路 19号 300201 摘要 首先给出了点电流源在均匀半空间中所产生电位的解析表达式；其次研究了源 电极 位于表层和底层时成层半空间介质中的电位分布 ；最后研究了源电极位于任一层时成层半空间 介质中的电位分布。 关键词 ： 深埋 电极 成层半空间 电位 解析表达式 [文章编号]1001．4683(2009)03．0246—10 [中图分类号]P631 [文献标识码]A 0 前言 40多年的地震前兆观测表 明，大震前一般在震源 区及其附近都会出现视电阻率 的变 化 ，被认为是一种 比较可靠的中短期地震前兆 (张国民等 ，2001)。然而 ，由于地电阻率观测 通常在固定台站的地表进行 ，降雨、气温、地下金属管线等因素都会影响近地表层的电阻率， 进而影响视电阻率的观测结果 ，给前兆 的判定带来很大困难 。另外 ，为了探测到覆盖层下岩 石的电阻率变化 ，通常采用大极距 的地表供 电方式。由于供 电极距较大 ，需要对环境进行保 护的布极区的范围也较大，使得对布极区的地震观测环境保护与当地经济和社会发展对观 测环境的破坏之间的矛盾越来越突出。 解决这些问题的途径之一是将现行观测系统进行深埋 ，进行所谓的全空间观测 (田山 等 ，2009)。我们一般可将实际的地球简化成层状半空间模型，这样 ，研究源电极位于成层半 空间中任一层时的电位分布问题便成为进行深埋电极地电阻率观测研究必须首先解决的问 题 。 许多研究者已研究过点电流源在层状半空间介质中产生的电位分布。由于大多数 电法 勘探是在地表进行的，因此 ，源电极位于地表时层状介质 中的电位分布一直是人们研究的重 点。例如，Raghuwanshi等(1986)研究了地表层电导率为常数，其它奇数序号层 电导率随深度 按指数函数变化 ，而偶数序号层 电导率随深度按幂 函数变化这样一种 2n+1层的介质模 型，当源电极位于地表时的电位分布问题 。Slaoui等 (2001)将多层土介质视作一系列双层介 质的串联 ，利用 Schlumberger装置的对称性及双层介质表面电位分布 的级数解 ，给出了地表 视电阻率的快速计算方法 。随着井 中电法勘探工作的开展 ，Daniels(1978)开始研究多层介 [收稿日期]2009．O1．21；[修定 日期]2009．07．21。 [项目类别]天津市科技支撑计划重点项 目(07ZCGYSF03100)“天津及周边地区地震中短期预测技术研究”资助。 [作者简介 】聂永安，男，1961年生，研究员，曾主要从事工程地震和震害预测研究，现主要从事地震监测预报的研究 与管理工作。E-mail：hie@mail．eq一．ac．cn 3期 聂永安等：成层半空间深埋 电极产生的电位分布 247 质 中深埋电极产生的电位分布问题，但给出的公式不够明确简捷 ，数值计算精度也不高。葛 为中等 (1990、1991)在 Daniels(1978)工作的基础上引入 了电阻率转换 函数 ()，给 出了确定 直立多层介质地面点源电场的递推公式，但这组递推公式应用于水平层状介质地面点源问 题时其结果是不正确的(葛为中，1994)。另外，在某些情况下，利用电阻率转换函数会使计 算无法进行。因此，有必要对这一方法进行改进和完善。 本文给出的汁算在任一层中深埋电极时成层半空间介质中的电位分布，其方法简洁、明 确、易于编程 ，计算稳定、效率高。当各层电阻率相同时，很容易证明其结果就是点电流源在 均匀半空间中所产生电位的解析表达式 。 1 均匀半空间中点电流源产生的电位分布 对于图 1所示的均匀半无 限空间，介质 的电阻率 为 p，供 电极 A