电阻率层析技术在地震监测中的应用.docxVIP

电阻率层析技术在地震监测中的应用 1 地下电阻率层析技术 直接地源测量一直是我国地震监测最基本的方法之一。每个站的常规观测统一使用一个独特的四极单极装置。到目前为止，几乎所有的伪造性异常都可以通过观测电阻曲线的直观分析来证实。地震的特点是，一些车站对6级以上的大地震反应较慢。这项工作的困难不仅是地震之前的观测异常，而且还有许多干扰。因此，简单的计算和经验判断应该是数据处理的基本方法。虽然观察是用土壤矩阵进行的，但在地震预报时，人们不知道地下电性结构是什么，也不知道物理参数是否发生了具体变化。信息的不足严重制约了工作的深度。 电阻率层析技术在勘探领域的成功为解决地电监测的困难提供了新的思路. 该技术的精髓在于通过阵列电极的扫描性测量和数据叠加,提高了深部探测的信噪比,再通过图像重建把低分辨率的视电阻率测量值转变成高分辨率的介质真电阻率的二维图像. 由于真电阻率的异常幅度(或前兆的信息强度)一般会比视电阻率的大,这对动态检测十分有利. Ushijima曾用垂直电法剖面技术得到流体运动的电位动态图像;自1992年起,Di Maio等在意大利的火山区用真电阻率层析图像进行了动态监测;随后,Lapenna等利用自然电位对意大利的地震区作了监测,在3—7km的深度上对异常电荷的发生概率重建出了层析图像,效果均很好. 本文给出的是我国首次试验结果. 工作开始于1996年底,在天津市宝坻台和河北省昌黎台各自架设了长度近1km的EW和NS两条复测剖面,研究期间曾于1998年4月14日发生ML5.0和ML4.4的两次唐山地震,震前宝坻和昌黎台的电阻率层析图像均出现了异常,震后消失,同时还在其他台站发现了地下水和小震活动的异常,这些变化可能是地震的前兆反应. 从2台同时开展的常规视电阻率资料中,却难以发现这些变化. 2 电流强度 非均匀各向同性介质中稳定直流电的电位U满足Poisson直流传导方程 [CX]Δ[CX]2U=?ρ([CX]Δ[CX]?j+[CX]Δ[CX]U?[CX]Δ[CX](1/ρ)),(1)[CX]Δ[CX]2U=-ρ([CX]Δ[CX]?j+[CX]Δ[CX]U?[CX]Δ[CX](1/ρ)),(1) (1)式的解,即空间中任意一点的电位U,可由Poisson积分给出 U=14π?V?ρ([CX]Δ[CX]?j+[CX]Δ[CX]U?[CX]Δ[CX](1/ρ))RdV,(2)U=14π?V-ρ([CX]Δ[CX]?j+[CX]Δ[CX]U?[CX]Δ[CX](1/ρ))RdV,(2) 这里,j为电流密度,ρ为介质的真电阻率,R为积分元到观测点的距离,V为积分域. (2)式右端被积函数的第1项代表外电流源的贡献,在空间的有源区ΔΔ·j=I0,I0为电流强度,在无源区为零;右端第2项反映介质的真电阻率非均匀分布的贡献,包括电阻率分界面上积累电荷对电位的影响;在均匀介质时ρ为常数,该项为零. (1)式在均匀介质的无源区变为同重力位一样的Laplace方程ΔΔ2U=0. 地电阻率法所给出的量ρa为 ρa=kUM?UNI,(3)ρa=kUΜ-UΝΙ,(3) ρa称之为视电阻率,以别于介质的真电阻率ρ,I为二供电电极(A,B)的供电电流,UM和UN分别为二测量电极处(M,N)的电位,k为与电极坐标参数有关的装置系数. 应指出,由于电流I为已知量,故地电方法的实测量仅仅是电位. 它有两个重要而又常被忽视的性质:(1)实测量U是真电阻率ρ及其梯度的积分之值,视电阻率ρa仅是这一积分差值的简单线性导出量;(2)实测量U是位函数,同重力位一样具有可叠加性. 上述两个基本特点决定了电法的实测量(U或ρa)必然都是数值稳定而对真电阻率的空间分辨很低的物理量. 从理论上讲,即使对于对称四极等一些良好的电法剖面装置,也都无法避免剖面两侧和两端以外任一干扰的影响,因为剖面外围的任一电源(ΔΔ·j)和任何真电阻率的非均匀分布(ΔΔ(1/ρ))都是位函数的两个基本贡献项,都会把干扰位叠加到实测量上. 在地震监测中,电法同地震波和电磁波的观测有着本质的不同:虽然地震波走时t是介质慢度μ的积分量t=∫μdl,电磁波强度E是介质吸收系数β的积分量E=E0exp(-∫βdl),但是后二者都不是位函数,且积分只涉及到发射源-接收点间的射线路径;而(2)式的三维积分却需遍及全空间. 在构造应力的作用下,即使真电阻率ρ在地震发生前出现了很大的变化,但只要在广大积分域(V)内同时分布有正和负的变化,那么积分值U仍然可以基本不变,从而导致视电阻率ρa变化极小甚至没有反应;相反,近距离的浅部干扰(距离R值很小)还会强烈地影响位函数U值. 简言之,根据1个测点的位函数去推断被积函数的多种变化,是得不到惟一解的. 问题的解决不能仅靠提高仪器精度和改进计算方法,关键在于观测系统.