复杂角域地形的视电阻率组合叠加校正方法.docxVIP

复杂角域地形的视电阻率组合叠加校正方法 0 模拟地形影响方法 20世纪50年代，随着《直接法》在矿产资源调查中的大规模应用，地形影响问题开始受到关注。理论与实践表明, 地形不但可以引起假异常, 而且会掩盖地下由矿体或目标物引起的真异常;如果不能正确地认识和消除, 就会导致完全错误的解释结果。目前获得纯地形视电阻率异常的方法有物理模拟、数值模拟和解析计算等方法。 物理模拟方法:Takeshi借用物理模拟方法研究电阻率曲线的地形影响。20世纪60—70年代, 徐世浙、林崇远、何继善等曾采用水槽模拟、土槽模拟、导电纸模拟等物理模拟方法研究地形影响问题。这一时期基本是定性解释阶段。物理模拟成本大, 速度慢且较为繁琐。 数值模拟:Oppliger利用面积分方程方法实现对充电法和磁电阻率法 (MMR) 3-D地形影响的模拟。徐世浙介绍起伏地形上2.5-D电场计算的边界元方法, 此后黄兰珍、汪晓东、田宪谟等为该方法的发展应用做了大量的贡献。Coggon针对直立导电极化体, 用有限元法计算得到平底、斜边的山谷和山脊2-D地形下电阻率及极化率响应曲线, 之后此方法得到广泛应用。有限差分法于20世纪60年代被用于应用地球物理领域, 90年代, 交错式样网格被广泛用于地电磁场的分析中来, 使有限差分法步入全盛时期。数值模拟法计算纯地形异常较为方便, 但计算量大, 在运用中不可能划分太多单元格, 因而导致精度降低;但随着计算机的发展, 此方法被广泛应用。 解析计算法, 如保角变换法、坐标网法、角域地形叠加法、格林理论、模拟电荷法等, 由于数学公式的复杂性和只能计算简单地形, 应用受到一定限制。 20世纪70年代, 何继善等研究了点源条件下角域地形上电位的解析解, 并绘制了电位畸变值量板。实践表明, 组合叠加法在满足观测精度的要求下, 方法简单, 计算速度快。笔者对组合叠加法的误差进行分析, 给出组合叠加法的应用条件, 并将利用组合叠加法而得到的复杂角域上的视电阻率曲线与ANSYS数值解进行对比。最后给出利用组合叠加方法获得纯地形影响后应用比较法进行地形改正的模型算例。 1 角域地形上m点电位解析 所谓角域, 是指由2个半无限大平面所限的二度区域, 在区域中充满电阻率为ρ的均匀且各向同性导电介质。两半平面之张角φ即为角域之顶角。两侧面与水平面之夹角为侧面之倾角, 分别用β和α表示左右两侧面之倾角。规定斜坡在水平轴下面, 坡角取负值;在上面取正值。则 0φ2πφ+α+β=π}。(1)0φ2πφ+α+β=π}。(1) 显然:当φπ时, 角域即为坡长无限之二度山脊;当φπ时, 角域即为坡长无限之二度山谷。角域是二度地形的基本单元, 在电阻率法的地形改正中很有用处。 现有一角域二度体, 角域的张角为φ, 在角域以下有一电阻率为ρ的均匀介质, 地面以上为空气, 其电阻率为无穷大, 在地面的A点有一电流强度为I的点电源供电,A到坐标原点的距离为r0。为了求出电位的分布, 首先建立坐标系, 将x轴取为过供电点A而垂直于地形两斜面的交线,x轴与地形两侧面交线的交点为O,y轴取为过O点而垂直于地形的水平面,z轴为地形两侧面的交线。坐标轴正向如图1。这样角域地形上M点电位解析解为 U(r,θ,z)=ρIπφrr0√[12Q?12(ξ)+∑m=1∞(?1)mcosmπ(φ?θ)φQmπφ?12(ξ)]。(2)U(r,θ,z)=ρΙπφrr0[12Q-12(ξ)+∑m=1∞(-1)mcosmπ(φ-θ)φQmπφ-12(ξ)]。(2) 式中:ξ=r2+r20+z22rr0ξ=r2+r02+z22rr0;Qv(ξ) 为第二类勒让德函数。 式 (2) 本身不完整, 存在以下问题: 1)Qv(1) 不存在, 故在此时解有奇异性, 不能用式 (2) 计算, 使用下式计算: U=ρIπφr0{π4+∑m=1∞(?1)m[φ(m+12)?φ(mπφ+12)]}。(3)U=ρΙπφr0{π4+∑m=1∞(-1)m[φ(m+12)-φ(mπφ+12)]}。(3) 其中,φ(x)=limn→∞[lnn?1x?1x+1?1x+2???1x+n]φ(x)=limn→∞[lnn-1x-1x+1-1x+2-?-1x+n]。 2) 电极在角域顶点时, 电位公式也可表示为 U=πφU0=πφIρ2πR=Iρ2φr2+z2√。(4)U=πφU0=πφΙρ2πR=Ιρ2φr2+z2。(4) 式中:U0为背景场。 当测量电极M在顶点时, 由电极互换原理得 U=Iρ2φr20+z2√。(5)U=Ιρ2φr02+z2。(5) 角域地形点电源电场电位由 (2) 、 (3) 、 (4) 、 (5) 几个电位公式及其附属公式计算。 为了突出地形影响, 用U/U0表示电位畸变值U(τ) 。在z=0时, 讨论电位畸变值是非