稳定电流场拟解析近似精度分析论文.pdfVIP

稳定电流场拟解析近似精度分析 孙建国 (吉林大学地球探测科学与技术学院地球物理系，长春，130026) 本文的目的是分析利用拟解析近似解决稳定电流场数值模拟词题时的误差来源和量级，并在此基础 上加深对拟解析近似的应用条件和适用范围的理解，为进一步的研究和发展基于拟解析近似的正反演方 案提供依据。 与其它方法不同，拟解析近似法通过引入电位反射函数和电场反射张量的概念将对电位、电荷、电 场和磁场的求解问题转化为对相应的反射函数和反射张量的求解问题。因此，对拟解析法的近似精度研 究要通过对反射函数和反射张量的计算精度研究来实现。 根据出发方程的不同，可以得到不同的反射函数和反射张量。对于电位，有下列三种反射函数：(I) 体电位反射函数：(2)面电位反射函数；(3)面电荷反射函数。根据定义，体电位反射函数只在异常地 质体的电阻率分布不均匀时才有定义，而面电位反射函数和面电荷反射函数和只在地质体的电阻率是常 数时有定义。对于电场，有下列两类电场反射张量：(I)基于磁并矢格林函数的反射张量；(2)基于电 并矢格林函数的反射张量。因此，在分析拟解析近似法的计算精度时要针对不同的反射函数和反射张量 分别讨论。 对于一个给定的电场反射张量或反射函数，通过构造其积分表达式与相应的拟解析近似之间的差， 可以得到一个描述拟解析近似误差的积分。不同的电场反射张量(函数)有着不同的误差积分。利用 Taylor展开的方式，可以将诸误差积分化成一系列与待求反射张量及其导数有关的级数，该级数的截断 误差代表着拟解析近似的误差。 对于拟解析级数，对其近似误差的分析可以利用积分方程的叠代理论来实现。事实上，如果取初值 为零，则可以将拟解析级数化为常规的积分方程叠代解。 ， 根据电磁理论，基于磁并矢格林函数的电场反射张量九(r)满足下列积分方程： 聃E舻掣驯一掣E舯掣[、业盟掣dr’ crtr) or(r) o-jr)“ Lr， ∥Oo 与此相对应的误差项为： o-0(r)fcro(r)G(r，r)·[E。(r)一E。(r)】J，· 丽)如——一。—历万广——Ⅲ 因此，拟解析近似的误差依赖于矢量差函数万E。(r)=E。(r)一E。(r)在异常体内的变化。如果在异 常体内万E。(r)很小，则拟解析近似的精度就会很高。否则，拟解析方法给出结果的近似程度就会很i氐。 这些结论同样适用于基于电格林函数的电场反射张量和基于标量格林函数的各种电位反射函数。 为了对拟解析法的近似精度有一个定量的了解，我们对均匀球体表面的电位和电荷反射函数进行 了定量分析。结果表明，与电荷和电位相对应的电位反射函数拟解析近似的精度是不同的。因此，应该 根据具体情况选取合适的电位反射函数。 参考文献 l field and Kaufman，A methods，Part I l A，Geophysicaltheory A，Gravitational，electric，and fields，AcademicPress，1992· magnetic theelectrical tensorin d．c．electricfield 12lSun，J．，On reflectivity 1 l一42 1，2005 Prospecting，V01．53，4 ·639·