直流电场拟解析近似理论实现技术研究
  • 【摘要】

    积分方程方法是基于Maxwell方程的积分形式的数值方法,它往往通过求解目标表面或者目标体内的等效源来间接获得场的解.积分方程方法是一种在解决三维电场数值模拟时有很多优点的方法,它可以避免微分方程方法中场值传递过程中的误差积累问题.出于对利用积分方程法快速进行三维模拟的需要,出现了很多的近似方法,如Born近似方法、局部非线性近似方法和拟线性近似方法等.拟解析近似方法是一种解决电磁场散射问题时提出... 展开>>积分方程方法是基于Maxwell方程的积分形式的数值方法,它往往通过求解目标表面或者目标体内的等效源来间接获得场的解.积分方程方法是一种在解决三维电场数值模拟时有很多优点的方法,它可以避免微分方程方法中场值传递过程中的误差积累问题.出于对利用积分方程法快速进行三维模拟的需要,出现了很多的近似方法,如Born近似方法、局部非线性近似方法和拟线性近似方法等.拟解析近似方法是一种解决电磁场散射问题时提出的快速求解积分方程的近似方法,它适用于强散射和大扰动问题,可以避免传统数值方法中遇到的求解大型代数方程组或大型矩阵问题,将其引入直流电场的数值模拟,可以加快三维直流电场模拟的速度.因此,本文研究了利用拟解析近似方法求解三维直流电场的积分方程的实现.首先本文对孙建国应用拟解析近似方法求解直流电场积分方程的理论公式进行了研究.从直流电场的基本理论出发推导出直流电场满足的积分方程,利用总场与背景场和异常场的关系得到异常电场的积分方程.根据拟线性近似方法的基本思想,推导直流电场的拟线性近似理论公式,对拟线性近似方法进行分析,从而引出直流拟解析近似方法.对直流电场的拟解析近似理论进行研究,在理论公式中的分析中,利用电反射张量的不同假设,将拟解析近似细分为两种近似方法,一种是假设电反射张量为标量,给出了标量拟解析近似的理论公式;另一种假设电反射张量为二阶张量,给出了张量拟解析近似的理论公式.根据拟解析近似理论的初始公式,对其进行改造,引申出迭代拟解析近似,得到了拟解析近似序列,通过拟解析近似序列迭代可求出积分方程n阶拟解析近似解,进而使拟解析近似方法达到更高的精度.接下来我们考虑在拟解析近似理论实现过程中遇到的几个亟待解决的问题.第一个是磁并矢格林函数的问题.在本课题中用到的格林函数是直流电场磁并矢格林函数,它与前人所用到的格林函数不同,尚未出现其数学表达式,因此,需要对其进行研究.我们从矢量位方程出发,利用位与场的关系,推导出全空间位并矢格林函数与磁并矢格林函数的关系,利用已知的全空间标量位格林函数进行推导,得到全空间磁并矢格林函数.第二个是数值积分方法的选取.在利用拟解析近似方法的过程中,一个无法避免的问题是数值积分问题,如何选取合适的数值积分方法成为拟解析近似方法实现的一个关键.通过对高斯积分的研究,我们发现高斯积分具有很多的优点,尤其是它易于推导得到多重求积公式,因此我们选取高斯积分公式来完成数值积分.通过研究发现,当高斯节点数为5时,它既能够保持很高的精度,又具有较快的计算速度.第三个是磁并矢格林函数积分奇异性问题.为了解决这个问题,我们对格林函数的奇异性处理方法进行了研究,对Yaghjian处理张量格林函数奇异性的方法进行了分析和总结后,提出了拟源并矢方法处理磁并矢格林函数的积分奇异性.这种方法将积分区域分为奇异点小邻域和无奇异性区域两部分,无奇异性区域可直接进行数值积分,对奇异点小邻域积分可以通过数学方法得到拟源并矢,然后将两部分积分相加,这样就在积分区域上消除了磁并矢格林函数的积分奇异性,使得格林函数奇异性问题得到解决.在解决了拟解析近似方法实现中遇到的难题后,本文对孙建国应用拟解析近似方法求解直流电场积分方程的理论进行了实现.在研究过程中,分别进行了标量和张量拟解析近似方法的实现,对不同形状异常体模型,不同的电导率变化,不同类型的场源进行拟解析近似模拟.在标量拟解析近似方法实现过程中,首先在均匀场源的情况下选取一系列不同的电导率对比的球体异常体进行计算,求得标量拟解析近似解,利用Born近似方法得到Born近似解.均匀场中的球体异常体模型的异常电场是具有解析解的,因而可以得到异常电场的解析解.然后将标量拟解析近似解、Born近似解分别同解析解进行对比,得到各自的误差.当背景电导率与异常电导率对比很小时,Born近似解和标量拟解析近似解都能较好的逼近解析解,但标量拟解析近似解比Born近似解更接近解析解.当背景电导率与异常电导率对比增大时,Born近似解的误差也逐渐增大,不能再适用于该模型的电场模拟,标量拟解析近似解的误差虽然也增大,但仍能满足误差要求,其解可以用于电场数值模拟.随着背景电导率与异常电导率对比的进一步增大,Born近似解已经完全失去了意义,标量拟解析近似解也超过了允许的误差范围,但我们利用标量拟解析近似迭代方法得到迭代解,一般情况下进行三到五次迭代,得到的结果就能满足误差要求,能够很好的与解析解相拟合.随后本文又对异常体为长方体、椭球体的地电模型进行了数值模拟,得到了长方体异常体和椭球体异常体产生的异常电场的拟解析近似解.接下来选取场源为点源场,我们同样选取了一系列的模型进行了模拟,并将标量拟解析近似解与Born近似解和解析解进行对比,可以看出标量拟解析近似方法在进行大扰动的地电模型模拟时有着很大的优势.同样我们也进行了长方体异常体和椭球体异常体在点源场中产生的异常电场数值模拟,得到了不同背景电导率与异常电导率对比情况下的异常电场.进行张量拟解析近似方法实现研究时,我们采用与标量拟解析近似实现研究相同的方式.可以得出下面的结论,张量拟解析近似方法能够很好的适应于大扰动地电模型的异常电场的数值模拟,具有很高的精度,能够用于计算如球体、长方体和椭球体等形状的异常体在均匀场和点源场中的响应.通过对拟解析近似方法的实现研究,可以看出该方法是一种可以用于求解三维直流电场的积分方程的近似方法,它的应用可以推进三维直流电场正反演模拟的进程,实现三维直流电法勘探数据的快速正反演,因此拟解析近似方法具有很好的应用前景. 收起<<

  • 【作者】

    张金会 

  • 【学科专业】

    地球探测与信息技术

  • 【授予学位】

    博士

  • 【授予单位】

    吉林大学

  • 【导师姓名】

    孙建国

  • 【学位年度】

    2010

  • 【语种】

    chi

  • 【关键词】

    直流电场%数值模拟%拟解析近似方法%磁并矢格林函数