一种时域有限差分的三维感应-极化双场数值模拟方法技术

本发明专利技术涉及一种时域有限差分的三维感应‑极化双场数值模拟方法，目的在于快速计算三维模型的感应‑极化双场的电磁响应。主要包括采用逆拉普拉斯变换获得迪拜模型电导率的时域表达式，构建电导率参数的e指数辅助方程，通过梯形积分法获得欧姆定律时域离散递推表达式，由四维数值运算降低为三维运算，再将其代入无源Maxwell旋度方程中，基于三维时域有限差分方法推导电场和磁场的迭代方程，进而完成三维模型的感应‑极化双场电磁响应数值计算。本发明专利技术目的在于解决欧姆定律的时域卷积运算耗时长、内存占用大等问题，最终实现三维模型的感应‑极化双场电磁响应快速数值计算。

【技术实现步骤摘要】

本专利技术涉及一种地球物理勘探领域中的时域电磁场数值模拟方法，尤其适用快速计算三维模型的感应-极化双场电磁响应。
技术介绍
瞬变电磁法(Time-domain Electromagnetic，简称TEM)是基于电磁感应涡流效应为主的地下近地表探测方法。当低频时忽略位移电流和极化电荷产生的电磁场，仅利用地下介质的涡流效应观测二次瞬变场。时域激发极化法(Time-domain InducedPolarization，简称IP)是基于地下介质极化效应为主的地下浅层探测方法，主要采用电性源施加直流或交流电磁场激励下，忽略传导电流产生的电磁场以及高频介电常数引起的极化响应，仅利用时域电流断开后低频激励下地下介质产生的极化场，观测极化电荷产生的二次电位变化过程。在时域激发极化测量中总是去除感应场的耦合。但是，无论是瞬变电磁还是时域激发极化方法，都是以电磁场麦克斯韦方程为理论依据并进行波场简化近似，实际探测中感应场和极化场同时存在，互相伴生。无论瞬变电磁还是时域激发极化方法，将两种电流场割裂开分别进行研究，是不符合实际电磁场的传播规律的，可能导致实测数据与理论计算结果不符，甚至导致解释结果出现错误。为此，只有准确计算地下结构的感应和极化双场电磁响应，理论模拟结果与实际测量数据互相吻合。中国专利CN104408021A公开了一种电偶源三维时域有限差分正演成像方法，对海洋空气、海水和海底大地三分空间均建立麦克斯韦方程组和本构方程，采用时域有限差分法得到海水和海底大地中任意时刻电磁场的分布。中国专利CN105277980A公开了一种高精度空间和时间任意倍数可变网格 有限差分正演方法，通过建立地下介质的正演速度模型，对正演速度模型中的声波波场进行二维网格离散化，同时对最佳匹配层边界条件进行网格离散化，进而通过声波波动方程进行时域有限差分正演模拟。加拿大专利CA2388271公开了一种基于时域有限差分方法的电磁场计算方法，并提出了一种基于二维导体的网格划分方法和数据储存方式，但只涉及感应场的数值计算。加拿大，D Marchant(2015，University of British Columbia)采用有限元方法，在时间域直接求解并模拟计算了三维空间下存在激发极化效应的电磁响应。以上所述方法公布了基于时域有限差分(FDTD)方法的电磁响应计算方法，国内外专利还未涉及时域有限差分方法计算三维模型的感应和极化双场的电磁响应方面的，为此，本专利技术基于FDTD方法直接计算三维模型的感应和极化双场电磁响应，并解决欧姆定律的时域卷积计算耗时长、内存占用大等问题，最终实现三维模型的感应-极化双场电磁响应快速数值计算。
技术实现思路
本专利技术所要解决的关键问题是提供一种时域有限差分的三维感应-极化双场数值模拟方法，先获得迪拜模型的电导率时域形式，通过构建电导率参数e指数辅助方程，并采用梯形积分法获得欧姆定律的时域离散递推形式，最后基于FDTD方法，实现了时域内快速计算三维模型的感应-极化双场电磁响应。本专利技术是这样实现的，一种时域有限差分的三维感应-极化双场数值模拟方法包括：1)、基于迪拜模型(频率相关系数c＝1时的Cole-Cole复电阻率模型)，先获得迪拜模型电导率的频域形式，再通过逆拉普拉斯变换，得到电导率的时域表达式σ(t)；2)、将电导率的时域表达式σ(t)代入欧姆定律中，获得欧姆定律时域卷积积分表达式，结合电导率σ(t)近似于αe-βt形式的特点，构建电导率参数的e指数辅助方程r(t)，代入欧姆定律时域卷积积分表达式，将其转换成线性积分形式；3)、将计算时间进行密集剖分，利用梯形积分法，获得欧姆定律时域离散递推表达式；4)、将欧姆定律时域离散递推表达式代入Maxwell方程中，基于三维时间域有限差分方法，在时间和空间上进行中心差分离散，推导电场E(t)、磁场H(t)的迭代方程；5)、采用非均匀三维Yee氏网格对计算区域进行剖分，设置电性参数，时间步长，并对各个网格进行电导率参数赋值；6)、计算初始场，对三维模型进行电场E(t)、磁场H(t)的迭代，加载狄利克雷边界条件，完成电磁场数值计算；7)、计算结束后，提取磁场或电场的各分量响应进行成图。进一步地，步骤2中，根据欧姆定律的频域微分形式，先通过频时变换得到欧姆定律时域卷积积分形式；再将电导率时域表达式代入，获得欧姆定律时域卷积积分表达式，如式(1)所示；根据电导率σ(t)近似于αe-βt的形式的特点，构建电导率参数的e指数辅助方程r(t)，如式(2)所示，将r(t)的表达式代入欧姆定律的卷积积分表达式中，得到欧姆定律的线性积分形式，如式(3)所示：式中J(t)是电流密度，E(t)为电场值，t为时间；电导率绝对值表达式为τ1为时间常数、η为极化率，σ∞为频率趋近无穷大时的电导率值。进一步地，步骤3中，将计算时间进行密集剖分，利用梯形积分法，得到积分求和ξ(tn)，提取ξ(tn)最后一项构建递推公式，如式(4)所示，进而实现降维计算；最后将式(4)代入式(3)，并对时间进行离散，最终得到n时刻欧姆定律时域离散递推表达式，如式(5)所示：式中n表示第n个时刻，r(tn)、E(tn)、ξ(tn)为第n个时刻所对应的数值，Δt(n)为第n个时间步长；接着在时刻对式(4)进行离散，再代入(5)式中，进行整理，得到时刻欧姆定律时域离散递推表达式，如式所示：进一步地，步骤4中，利用步骤3中的推导结果，将式(6)中的项，采用平均值近似法进行计算，有代入无源Maxwell方程，采用时域有限差分方法，构建出电场E(tn+1)与E(tn)的迭代方程；对于中的ξ(tn)项，先按照步骤3中式(4)的递推公式进行计算，获得n时刻的ξ(tn)值，再将其代入电场迭代公式，最终获得n+1时刻电场的迭代方程，如式(8)所示：式中有(i,j,k)为Yee式网格的坐标点，为n+1时刻x轴方向电场分量，Ex(n)(m)为n时刻x轴方向电场分量，为时刻对应坐标点上的y轴方向的磁场分量，为时刻对应坐标点上的z轴方向的磁场分量，Δyj、Δzk为Yee氏网格的y、z方向的步长。进一步地，步骤5中，根据空间、时间的先后顺序进行离散，先对三维Yee网格的σ∞(x,y,z)赋值，将每个网格σ∞(x,y,z)代入到电导率时域表达式中，得到σ(x,y,z,t)；再结合电流密度的计算时刻特点，设定时间步长Δt(n)，将时间密集剖分成N+1个时刻，计算时间半步长再对时间进行时刻剖分，对应为t(0)、t(n)，共2N+1个时刻，在2N+1个时刻分别对r(t)、进行赋值后，存储r(t)的2N+1个值。对于的数据，先舍弃n时刻的值，再保留时刻的值，这样只需存储N个数据。进一步地，步骤6中，在三维模型的基础上，建立地面中心回线发射与接收线圈的系统，计算出初始场后，代入到控制方程中，计算x、y、z方向的电场与x、y方向的磁场，为了确保计算结果的唯一性，采用散度方程计算z方向的磁场，应用狄利克雷边界条件，完成三维模型的感应-极化双场电磁响应 数值计算。本专利技术与现有技术相比,有益效果在于，将时域电磁探测中的感应和极化双场同时进行理论计算，通过构建辅助方程，推导欧姆定律的时域离散递推表达式，对计算场量进行降维处理，避免了卷积和的求取，大大减少了内存的使用本文档来自技高网...

【技术保护点】
一种时域有限差分的三维感应‑极化双场数值模拟方法其特征在于，包括如下步骤：1)、基于迪拜模型(频率相关系数c＝1时的Cole‑Cole复电阻率模型)，先获得迪拜模型电导率的频域形式，再通过逆拉普拉斯变换，得到电导率的时域表达式σ(t)；2)、将电导率的时域表达式σ(t)代入欧姆定律中，获得欧姆定律时域卷积积分表达式，结合电导率σ(t)近似于αe‑βt形式的特点，构建电导率参数的e指数辅助方程r(t)，代入欧姆定律时域卷积积分表达式，将其转换成线性积分形式；3)、将计算时间进行密集剖分，利用梯形积分法，获得欧姆定律时域离散递推表达式；4)、将欧姆定律时域离散递推表达式代入Maxwell方程中，基于三维时间域有限差分方法，在时间和空间上进行中心差分离散，推导电场E(t)、磁场H(t)的迭代方程；5)、采用非均匀三维Yee氏网格对计算区域进行剖分，设置电性参数，时间步长，并对各个网格进行电导率参数赋值；6)、计算初始场，对三维模型进行电场E(t)、磁场H(t)的迭代，加载狄利克雷边界条件，完成电磁场数值计算；7)、计算结束后，提取磁场或电场的各分量响应进行成图。

【技术特征摘要】
1.一种时域有限差分的三维感应-极化双场数值模拟方法其特征在于，包括如下步骤：1)、基于迪拜模型(频率相关系数c＝1时的Cole-Cole复电阻率模型)，先获得迪拜模型电导率的频域形式，再通过逆拉普拉斯变换，得到电导率的时域表达式σ(t)；2)、将电导率的时域表达式σ(t)代入欧姆定律中，获得欧姆定律时域卷积积分表达式，结合电导率σ(t)近似于αe-βt形式的特点，构建电导率参数的e指数辅助方程r(t)，代入欧姆定律时域卷积积分表达式，将其转换成线性积分形式；3)、将计算时间进行密集剖分，利用梯形积分法，获得欧姆定律时域离散递推表达式；4)、将欧姆定律时域离散递推表达式代入Maxwell方程中，基于三维时间域有限差分方法，在时间和空间上进行中心差分离散，推导电场E(t)、磁场H(t)的迭代方程；5)、采用非均匀三维Yee氏网格对计算区域进行剖分，设置电性参数，时间步长，并对各个网格进行电导率参数赋值；6)、计算初始场，对三维模型进行电场E(t)、磁场H(t)的迭代，加载狄利克雷边界条件，完成电磁场数值计算；7)、计算结束后，提取磁场或电场的各分量响应进行成图。2.按照权利要求1所述的一种时域有限差分的三维感应-极化双场数值模拟方法，其特征在于：步骤2中，根据欧姆定律的频域微分形式，先通过频时变换得到欧姆定律时域卷积积分形式；再将电导率时域表达式代入，获得欧姆定律时域卷积积分表达式，如式(1)所示；根据电导率σ(t)近似于αe-βt的形式的特点，构建电导率参数的e指数辅助方程r(t)，如式(2)所示，将r(t)的表达式代入欧姆定律的卷积积分表达式中，得到欧姆定律的线性积分形式，如式(3)所示：J→(t)=σ∞·E→(t)-∫0tσ^(t-τ)E→(τ)dτ---(1)]]>r(t)=etτ1(1-η)---(2)]]>J→(t)=σ∞·E→(t)-σ^(t)·∫0tr(τ)E→(τ)dτ---(3)]]>式中J(t)是电流密度，E(t)为电场值，t为时间；电导率绝对值表达式为τ1为时间常数、η为极化率，σ∞为频率趋近无穷大时的电导率值。3.按照权利要求1所述的一种时域有限差分的三维感应-极化双场数值模拟方法，其特征在于：步骤3中，将计算时间进行密集剖分，利用梯形积分法，得到积分求和ξ(tn)，提取ξ(tn)最后一项构建递推公式，如式(4)所示，进而实现降...

【专利技术属性】
技术研发人员：嵇艳鞠，吴燕琪，关珊珊，黄廷哲，吴琼，王远，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林;22