时间域航空电磁响应数值模拟的精度分析

分析说明了G-S变换和汉克尔变换的计算精度,当n=12或16时,G-S变换适用范围广且能保证精度;而汉克尔变换计算的相对精度达到10-6数量级;并且基于时间域航空电磁法偶极—偶极装置一维正演算法,在选取不同参数时,对航空电磁响应数值模拟精度进行了分析。结果表明G-S变换n值的选取对计算精度影响相对较大,而汉克尔变换中滤波参数的选取对计算精度影响较小。     

汉克尔变换精度分析及在电磁法探测中的应用

汉克尔变换的数值计算方法是计算电磁响应的重要工具.这里重点对比了五种已知汉克尔变换数值滤波方法在不同特性积分函数下的计算精度,以及在电磁法探测计算中的应用.结果表明,五种滤波方法均能较好地应用于电磁探测正演计算,计算精度属Anderson 801滤波较高.     

基于GPU实现汉克尔变换并行计算

地球物理勘探技术日新月异,地球物理勘探数据的处理和解释对高性能计算机的要求越来越高.相比于地震勘探,重力、磁法、电法勘探中的并行计算研究还都处于起步阶段.基于GPU的并行计算能够提供强大的计算能力和存储器带宽,同时具有良好的可编程性、较低的成本和较短的开发周期.这里实现了瞬变电磁法一维正演计算中汉克尔变换基于GPU的并行计算,比较了汉克尔变换串行算法和并行算法的计算耗时,基于GPU技术的并行计算相比串行计算,获得了很高的加速比.     

基于快速汉克尔变换算法的CSAMT一维正演

水平层状介质上的电偶源产生的电磁场数值解是含有零阶和一阶贝塞尔函数的无穷积分,这类积分实质上是汉克尔变换式。通过3种方案的计算实例,说明快速汉克尔变换滤波系数的选取对正演结果的影响。计算结果显示:加长滤波系数能明显提高计算精度,使正演结果真实地接近实际;相反,滤波系数区间过短会出现场值的畸变,视电阻率会受到较大影响。因此,合理选取滤波系数可以提高计算精度和速度。      

Guptasarma线性滤波算法在CSAMT一维正演中的应用

CSAMT法一维正演问题,其核心运算是对汉克尔变换的求解,目前主要使用数值滤波方法进行求解,求解过程中使用给定长度的滤波系数进行累加运算.滤波系数的选取,影响着正演计算的精度和速度.这里使用Guptasarma线性滤波系数,并给出具体实现方法.经过正演模型验证,发现这套方案可以在较短的滤波系数长度基础上,与正演模型拟合...     

汉克尔自适应算法在电磁计算中的应用

对于地球物理电磁计算中常出现的贝塞尔函数积分,对比几种高精度的快速汉克尔长滤波系数的精度,选择出一种对于一般情况都表现较好的汉克尔滤波系数,在其基础上,进行一个自适应的截断,来计算相关核函数的积分问题。这样既能保证计算的精度,又能节省计算时间。通过计算实例表明,该算法在保证精度的前提下,相比汉克尔滤波直接计算方法,能够节省大量的时间,尤其对于积分核函数相对复杂的情况,节省时间量尤其显著。在地球物理数据计算量越来越繁琐的今天,该算法对节省计算机资源,具有可行性和实用性。     

CSAMT一维正演可视化设计与激发极化效应

针对可控源音频大地电磁法(CSAMT)开发了一套CSAMT一维数据正演可视化软件.其界面友好,操作简单;使用MATLAB进行核心算法的编写,正演计算过程引入线性滤波系数对汉克尔积分进行求解,并对不同线性滤波系数方案的精度进行了比较;基于Java语言进行操作界面设计,提供多参数设置界面,自动化程度高;通过软件自带的绘图功...     

频率测深水平磁场的正演计算及应用

本文讨论了频率测深电偶源发射时水平磁场的正演计算问题,所用方法为快速汉克尔变换法。该方法的基本思路是:引入适当的变量替换,将汉克尔积分变为褶积形式,通过抽样再将其离散化为数字滤波形式,这样就可应用快速汉克尔变换滤波方法进行计算。本文讨论了方法的基本原理、计算电偶源发射时水平磁场滤波系数的选取、精度问题,并且讨论了利用电磁场之间的比值定义的视电阻率响应的正演计算问题,给出了一些计算实例。     

频率测深正演问题滤波算法及人机联作反演

频率测深一维问题正演常用的计算方法是数值积分方法,但是这种方法由于利用求和代替积分使得计算速度很慢,精度也较低。本文提出了频率测深一维问题正演计算的一种有效方法——快速汉克尔变换数字滤波法,较系统地讨论了快速汉克尔变换滤波系数的解析计算及其在频率测深正演计算中的应用。最后应用上述正演计算方法进行了频率测深的人机联作反演的研究。     

用快速汉克尔变换计算各向异性层状介质中点电流源的电位响应

针对任意各向异性地层,利用快速汉克尔（Hankel）变换导出了点电流源在地下介质中的电位响应关系,计算了各向异性地层的直流电视电阻率,重点分析了电阻率测深方法对地下各向异性介质的探测能力。这里采用二维傅立叶变换,首先导出了均匀各向异性层的电位一般解函数,利用各层界面电位和电流密度垂向分量的连续性及地下电流源的耦合条件,推导出不同层之间电位和电流密度垂向分量的传递函数,建立了空间电位和电流密度的递归计算关系。针对源点附近电位函数的奇异性,采用柱坐标变换,将二维傅立叶变换转化为汉克尔（Hankel）变换,并用快速汉克尔（Hankel）变换计算所有的位函数值。采用均匀各向异性半空间的解析解,验证了文中算法的正确性,给出了多层各向异性地层模型的视电阻率响应曲线,分析了直流电法探测裂缝性地层,估计裂缝分布性状的可能性。