利用积分方程法的大地电磁三维正演.pdfVIP

2 8 2 V o.l 2 8 N o. 2 2006 6 NORTHW ESTERN SE ISMOLOG ICAL JOU RNAL June, 2006 徐凯军, 李桐林, 张 辉, 李建平 ( 吉林大学地球探测科学与技术学院, 吉林 长春 130026 : 利用积分方程法实现了均匀导电半空 三维大地电磁响应的数值模拟求取张量格林函数 积分时, 采用二次剖分算法解决计算中奇异值问题, 对于含有贝塞尔函数的积分项, 利用结合连分 式展开的高斯求积代替常规的快速汉克尔变换方法, 确保了张量格林函数的正确计算并提高了计 算精度最后通过数值模拟结果的对比及模型试算验证了算法的正确性, 所实现的三维大地电磁 数值模拟算法为理论研究三维地电构造的大地电磁响应的分布规律提供了有效的工具, 也为研究 三维反演算法奠定了基础 : ; ; ; : P 631. 3+ 1 : A : 1000 - 0844 ( 2006 02- 0 104 - 04 Three D im entionalM agnetotelluric Forward M odeling U sing Integral Equation XU Ka-i jun, L I T ong- lin, ZHANG H ui, L I Jian-p ing ( College of GeoExp loration Science and Technology, Jilin University, Changchun 130026, China) Abstract: U sing integral equat ion m ethod, the num erical algo rithm of 3D m agneto tellur ic responses in hom ogeneous ha lf- space is rea lized. W h ile calculat ing tensor Green s functions, the algorithm of seconda- ry d iv ision is adopted to so lve the singu lar problem. T he Gau ssian quadrature and continued fract ion is a- dopted to calcu late the integralw ith Bessel function. T h is m ethod not only can ensure the right calcu la- t ion of tensor G reen s funct ions bu t a lso can get high precision. The correctness o f algor ithm has been proved by m odel com paring and ana lyse. Th is effic ient a lgorithm w ill he lp to study the d istribut ing law s o f 3D m agenetotelluric responses and to setup basis for reaserch of 3D inversion. K ey w ord s: Integral equation ; M agnetote lluric ; 3D orw ard m odeling ; Ten sor Green s unction 0 , , , ,