-
公开(公告)号:CN114912310B
公开(公告)日:2024-04-12
申请号:CN202210375579.4
申请日:2022-04-11
Applicant: 中南大学
IPC: G06F30/23 , G06F111/10
Abstract: 本发明提供一种基于正则化修正方程的三维大地电磁数值模拟方法,该方法包括:根据目标体圈定研究区域,构建研究模型以及天然电磁场的控制方程,采用矢量有限元离散;用模型的电导率乘以电场强度,连续求散度和梯度,带入天然电磁场的控制方程,获得新的电磁场控制方程;修正模型内部棱边的电场值,达到电流密度的散度为零的条件,并将加权余量方程用节点有限元法离散;采用迭代求解法获得各个测点的电场分量,并求解磁场分量,以此计算视电阻率和相位。本发明的优点在于,在异常体电导率不同的情况下,通过正则化修正可以避免对单元体的积分,并能有效解决单元接触表面电流不连续的问题,简化了积分难度,满足电磁数据精细、反演成像快速的需求。
-
公开(公告)号:CN114912310A
公开(公告)日:2022-08-16
申请号:CN202210375579.4
申请日:2022-04-11
Applicant: 中南大学
IPC: G06F30/23 , G06F111/10
Abstract: 本发明提供一种基于正则化修正方程的三维大地电磁数值模拟方法,该方法包括:根据目标体圈定研究区域,构建研究模型以及天然电磁场的控制方程,采用矢量有限元离散;用模型的电导率乘以电场强度,连续求散度和梯度,带入天然电磁场的控制方程,获得新的电磁场控制方程;修正模型内部棱边的电场值,达到电流密度的散度为零的条件,并将加权余量方程用节点有限元法离散;采用迭代求解法获得各个测点的电场分量,并求解磁场分量,以此计算视电阻率和相位。本发明的优点在于,在异常体电导率不同的情况下,通过正则化修正可以避免对单元体的积分,并能有效解决单元接触表面电流不连续的问题,简化了积分难度,满足电磁数据精细、反演成像快速的需求。
-
公开(公告)号:CN113779818B
公开(公告)日:2022-02-08
申请号:CN202111347154.4
申请日:2021-11-15
Applicant: 中南大学
IPC: G06F30/20 , G06F111/10
Abstract: 三维地质体其电磁场数值模拟方法、装置、设备及介质,包括:构建内部包含勘探目标的三维长方体模型;对三维长方体模型进行网格剖分,剖分成若干个小长方体单元,给每个小长方体单元的电阻率赋值,得到三维异常体模型;根据频率参数和三维异常体模型,构建对应的电场控制方程,并计算电场控制方程的系数矩阵和右端项;将三维异常体模型沿Z方向划分为多层子区域,在各层子区域分别构建对应的插值算子,并施加对应的插值算子给各层子区域;将各层子区域对应的插值算子合成后得到一个总的稀疏采样算子,并计算新的系数矩阵和右端项,由新的系数矩阵和右端项计算得到电磁场。该方法能够满足大规模电磁数据精细、快速正演成像的需求。
-
公开(公告)号:CN113779818A
公开(公告)日:2021-12-10
申请号:CN202111347154.4
申请日:2021-11-15
Applicant: 中南大学
IPC: G06F30/20 , G06F111/10
Abstract: 三维地质体其电磁场数值模拟方法、装置、设备及介质,包括:构建内部包含勘探目标的三维长方体模型;对三维长方体模型进行网格剖分,剖分成若干个小长方体单元,给每个小长方体单元的电阻率赋值,得到三维异常体模型;根据频率参数和三维异常体模型,构建对应的电场控制方程,并计算电场控制方程的系数矩阵和右端项;将三维异常体模型沿Z方向划分为多层子区域,在各层子区域分别构建对应的插值算子,并施加对应的插值算子给各层子区域;将各层子区域对应的插值算子合成后得到一个总的稀疏采样算子,并计算新的系数矩阵和右端项,由新的系数矩阵和右端项计算得到电磁场。该方法能够满足大规模电磁数据精细、快速正演成像的需求。
-
-
-