体积分方程法模拟复电阻率三维体电磁响应  

Modeling electromagnetic responses of complex resistivity 3-D body using volume integral equation method

作　　者：张辉[1] 李桐林[1] 董瑞霞 

机构地区：[1]吉林大学地球探测科学与技术学院,吉林长春130026 [2]中油集团东方地球物理公司研究院乌鲁木齐分院,新疆乌鲁木齐830016

出　　处：《煤田地质与勘探》2006年第1期70-74,共5页Coal Geology & Exploration

基　　金：国土资源部重点科技项目(20010212)

摘　　要：利用体积分方程法计算了均匀半空间中复电阻率(激电)效应和电磁效应同时存在时的三维体响应。在计算中,对于需要计算三重积分的张量格林函数电荷项一次部分,应用一种差分近似的方法求解,这种方法在保证计算精度的同时更加便于计算机实现;采用二次剖分的算法解决了计算张量格林函数时的奇异值问题;计算含有贝塞尔函数的积分项时,利用一种结合连分式展开的高斯求积代替常规的快速汉克尔变换方法。验证了计算结果并分析了三维体复电阻率模型(Cole-Cole)参数对正演结果的影响,为三维体复电阻率及其参数反演提供了正演依据。Using volume integral equation method to calculate 3-D body concurrent responses of induced polarization and electromagnetic in homogeneous half-space. The calculation of tensor Green＇s functions is discussed. We adopted difference approximate method fur the charge primary part of tensor Green＇s functions. This method not only can ensure the calculation accuracy but is more convenient for computation as well. Secondly, using the algorithm of secondary division solved singular problem in calculation of tensor Green＇ s functions. Using the Gaussian quadrature and continued fraction calculates the integral with Bessel function. The modeling results are compared with foreign works, all results prove that the procedure is correct. And results are analyzed in different parameter of complex resistivity.

关 键 词：体积分方程法 三维 复电阻率 电磁响应 张量格林函数 

分 类 号：P631.325[天文地球—地质矿产勘探]