基于快速汉克尔变换算法的CSAMT一维正演
Linear dimension forward solution of CSAMT by Hankel transformation method
-
摘要: 水平层状介质上的电偶源产生的电磁场数值解是含有零阶和一阶贝塞尔函数的无穷积分,这类积分实质上是汉克尔变换式。通过3种方案的计算实例,说明快速汉克尔变换滤波系数的选取对正演结果的影响。计算结果显示:加长滤波系数能明显提高计算精度,使正演结果真实地接近实际;相反,滤波系数区间过短会出现场值的畸变,视电阻率会受到较大影响。因此,合理选取滤波系数可以提高计算精度和速度。Abstract: The numerical solution of an electic dipole on a layered medium is materially Hankel-transform.By working three models of numerical solutions,which explain how the electromagnetic fields are effected by the Hankel transformation coefficient,then we know adding filtering factor can increase computation precision apparently and make the result of forward solution reflect actual situation reliably;otherwise,if the filtering range is too short the field data and apparenty resistivity may be impacted largely.Therefore,select correct filtering factor can improve the computation precision and computation speed.