[发明专利]一种基于超收敛插值逼近的三维电磁慢扩散数值模拟方法

说明书

本发明涉及一种基于超收敛插值逼近的三维电磁慢扩散数值模拟方法，通过将复电导率模型引入频域Maxwell方程组后，电磁场扩散方程中含有复频变量的负分数次幂项，先进行频‑时变换得到含有Caputo分数阶微分项的时间域控制方程；再采用Alikhanov超收敛插值逼近方法，对电场控制方程中Caputo分数阶导数进行超收敛逼近，获得分数阶微分项的非均匀步长离散近似表达式，从而完成时间域分数阶微分项的稳定、高精度直接求解；最后基于有限差分算法对控制方程进行离散，推导出电场和磁场迭代方程，最终实现了三维时域电磁慢扩散的高精度数值模拟。本发明目的在于克服分数阶微分求解的弱奇异不稳定性及误差较大问题，实现三维时域电磁慢扩散的高精度数值模拟。

技术领域

本发明涉及一种基于超收敛插值逼近的三维电磁慢扩散数值模拟方法，适用于时域电磁慢扩散数值模拟，尤其是对磁源电磁慢扩散的数值模拟。

背景技术

磁源瞬变电磁探测方法，选用线圈作为发射源，加载阶跃电流激发地下良导体产生涡流，通过接收机采集二次感应场获得地下丰富的电性信息，由于其具有探测深度大、经济、简便的优点，已经被广泛应用于地质资源探测中。超导量子干涉仪(SQUID，Superconducting Quantum Interference Device)和原子磁力仪的逐渐成熟与实际应用，显著提高了勘探的精度，获得了实际地质中更加准确的电磁衰减信号，不符合理论-5/2幂律的瞬变电磁慢扩散现象逐渐被观测，随着对地质结构特征的不断研究，传统的单一尺度、单一参数的地质电导率模型已经不能满足高精度地质探测的要求。建立更符合地质构造的物理模型是提高电磁探测分辨率的关键技术之一。

三维时域有限差分(Three-dimensional Finite-difference Time-domain)方法是电磁场数值计算的重要方法之一，它的基本思想是场量对时间和空间的一阶偏导数用中心差分近似,通过在时域的递推模拟波的传播过程,从而得出场分布，与一维和二维差分方法相比，三维差分更适用于观察域广、结构复杂的地质探测分析中，目前已被广泛应用于电磁场数值模拟计算中。

Chester J.Weiss等在美国布拉索斯郡利用TEM47观测到时域电磁慢扩散响应，并排除了慢扩散是由分层结构或横向异性导致的可能性，得出慢扩散响应是由粗糙介质引起的一种分数阶扩散。Mark E.Everett等人提出介质具有非均匀性，即粗糙度，并将粗糙参数引入到频率域电导率的表达式中。Detwiler等研究了裂隙中溶质的迁移问题，发现裂隙介质中溶质迁移是裂隙内部的泰勒扩散和裂隙表面的常次扩散的双尺度耦合扩散；在地质问题中，大尺度粗糙介质裂隙介质中电磁波的传播也具有相同的双耦合性质，即常次扩散和次扩散的双尺度耦合扩散。国内主要对电导率随空间随机变化的结构进行研究，很少对电导率随时间变化的慢扩散开展研究。

中国专利CN107766666A公开了一种基于分数阶差分法的三维时域电磁反常扩散模拟方法,通过频时转换获得电磁场的时域分数阶微分-积分表达式；采用Riemann-Liouville分数阶积分和有限差分方法,对扩散方程的积分和微分项进行时域离散,构建电、磁场的时域迭代公式；加载初始条件和边界条件,实现了三维时域电磁反常扩散的数值模拟。实现了三维时域电磁反常扩散的数值模拟。

中国专利CN 106776478A公开了一种反常扩散中的基于分步计算的离散分数阶差分方法，将时间空间反常扩散方程的离散格式中的两个gamma函数的比值看做一个参数,利用gamma函数的递推关系,结合数值模拟中的分布计算的思路和整体考虑的思想减小了gamma函数计算的限制，扩展模拟的点数，提高模拟的效率。

中国专利CN108897052A公开了一种基于分数阶线性近似的三维时域电磁慢扩散模拟方法，将广义电导率引入Maxwell方程组,将分数阶微分项进行线性近似，并进行频时转换得到近似的时间域的整数阶微分，完成从分数阶到整数阶的转换，达到模拟三维时域电磁慢扩散现象的目的。