[发明专利]一种基于库伦规范的三维电性源数值模拟方法

说明书

本发明公开了一种基于库伦规范的三维电性源数值模拟方法，在进行背景场计算时，提前计算了bessel积分中的层状系数，并在求解二次电磁场时采用傅里叶变换，将三维偏微分方程转化为波数下的一维常微分方程，最后采用紧算子修正求解过程中的电场。该方法基于二次场法，无需进行散度校正，利用这个优势，利用库伦规范将麦克斯韦方程组转换为库伦标量位矢量位的控制方程，在计算背景场时采用一种快速方法、二次场计算时利用傅里叶变换技术，极大地提高了三维电性源电磁场数值模拟的计算效率。本发明应用于数值模拟领域。

技术领域

本发明数值模拟涉及技术领域，具体是一种基于库伦规范的三维电性源数值模拟方法。

背景技术

人工源电磁法能克服大地电磁场源信号弱的缺点，抗干扰能力强，在矿山资源、天然气、地热和水文环境工程勘探等领域中得到了广泛的应用。复杂介质中电性源电磁场的解析解难以求得，所以频率域三维正演主要以数值模拟为主，目前常用的三维正演方法主要是有限差分法(Streich R.2009.3D finite-difference frequency-domain modelingofcontrolled-source electromagnetic data:direct solution and optimization forhigh accuracy.Geophysics,74(5):

F95-F105.)、有限单元法(刘颖，李予国，韩波.2017.可控源电磁场二维自适应矢量有限元正演模拟.地球物理学报，60(12)：4874-4886,)、有限体积法(彭荣华，胡祥云，韩波等.2016.基于拟态有限体积法的频率域可控源三维正演计算.地球物理学报，59(10)：3927-3939.)和积分方程法(汤井田，周峰，任政勇等.2018.复杂地下异常体的可控源电磁法积分方程正演.地球物理学报，61(4)：1549-1562)。近年来，三维正演的效率仍然是研究热点，高效、高精度的电性源三维电磁法正演技术仍有很大的发展空间。

发明内容

针对上述现有技术中的不足，本发明提供一种基于库伦规范的三维电性源数值模拟方法，解决了现有技术中复杂介质中电性源电磁场的解析过程计算量大、存储需求高、计算效率低等问题。

为实现上述目的，本发明提供一种基于库伦规范的三维电性源数值模拟方法，包括如下步骤：

步骤1，对研究区域沿x、y、z方向进行六面体剖分，其中，在x、y方向上为均匀剖分，在z方向上为均匀或非均匀剖分，x、y、z方向的剖分节点数量分别为N_x、N_y、N_z；

步骤2，获取地电模型参数，包括背景层参数、计算频率、异常导电率参数，并基于计算频率、源的位置和给定的背景层参数得到背景电场和背景磁场；

步骤3，将麦克斯韦方程组采用库伦规范转换为关于矢量位和标量位的方程，再利用二次场方法得到二次场矢量位标量位满足的控制方程，将背景电场与背景磁场的乘积作为控制方程中右端项的散射电流；

步骤4，将控制方程进行水平方向二维傅里叶变换，得到空间-波数域矢量位标量位的常微分方程，采用高斯偏移波数，对每个波数下的常微分方程用伽辽金法加权、采用有限元求解，得到每个波数下的空间-波数域二次场矢量位标量位；

步骤5，基于空间-波数域二次场矢量位标量位与电场在空间-波数域之间的关系式，得到空间-波数域二次电场，再利用水平方向二维反傅里叶变换，得到空间域二次电场；

步骤6，将背景电场与空间域二次电场相加得到空间域总电场，并判断空间域总电场是否满足给定的迭代收敛条件：

若满足，则基于空间-波数域二次场矢量位标量位得到空间-波数域二次磁场，水平方向二维反傅里叶变换后得到空间域二次磁场，加上背景磁场，得到空间域总磁场，再输出空间域总电场和空间域总磁场；