-
公开(公告)号：CN113779816A
公开(公告)日：2021-12-10
申请号：CN202111325242.4
申请日：2021-11-10
Applicant: 中南大学
IPC: G06F30/20 , G06F111/10
Abstract: 本申请涉及一种基于微分法的三维直流电阻率法数值模拟方法。所述方法包括：通过将基于微分法的异常电位方程式沿着水平方向做二维傅里叶变换，将三维偏微分方程转换为多个一维常微分方程进行求解，在求解复杂地质条件模型时，可有效减少计算量和存储需求，计算效率高，且具有较好的并行性；另外，在得到各个剖分节点当前的总电场后，通过紧算子对总电场进行修改，得到修改后的总电场，紧算子的引入使得高精度的数值解快速得到。本发明给出了一种基于微分法的三维直流电阻率法数值模拟方法，相对于积分方程法，微分法公式简单，适应性更强。本发明所提数值模拟方法具有高效率、高精度以及适应性强的特点。
-
公开(公告)号：CN113656750A
公开(公告)日：2021-11-16
申请号：CN202111218085.7
申请日：2021-10-20
Applicant: 中南大学
Abstract: 基于空间波数域的强磁介质的磁感应强度计算方法，剖分包含异常体的初始棱柱体模型，得到一系列的节点，对每个节点进行磁化率赋值，磁化率为标量；设定高斯参数，计算离散偏移波数，计算空间域背景场磁场强度，获得磁化强度模型；将空间域异常场磁位和磁化强度满足的方程转为空间波数域一维常微分方程；结合设定的空间波数域异常场磁位的边界条件，将空间波数域异常场磁位满足的边值问题模型等价转化为变分问题模型并求解，得到空间波数域异常场磁位、空间波数域异常场磁场强度后将其转换为空间域异常场磁位及空间域异常场磁场强度，根据输出的空间域异常场磁场强度求得空间域磁感应强度。本发明能更准确的对强磁性介质进行磁场数值模拟。
-
公开(公告)号：CN111967169B
公开(公告)日：2021-01-01
申请号：CN202011127944.7
申请日：2020-10-21
Applicant: 中南大学
IPC: G06F30/20 , G06F17/14 , G06F111/10
Abstract: 本发明提出了一种二度体重力异常积分解数值模拟方法、装置、计算机设备和存储介质。本发明通过重力模型表示、波数域重力异常公式推导、波数域非均匀采样、空间域重力异常场计算等步骤，利用非均匀采样快速傅里叶变换实现了重力异常场高效、高精度数值模拟。在本发明中，将非均匀采样傅里叶变换应用于重力异常模拟，充分融合了快速傅里叶变换的效率优势和高斯傅里叶变换的精度优势，有效兼顾了重力异常数值模拟的计算精度与计算效率，解决了波数域重力异常数值模拟只适合规则测线、数值模拟方法不能同时兼顾计算效率和计算精度的问题。
-
公开(公告)号：CN114004127B
公开(公告)日：2024-08-09
申请号：CN202111303885.9
申请日：2021-11-05
Applicant: 中南大学
IPC: G06F30/23 , G06F17/13 , G06F111/10
Abstract: 二维主轴各向异性强磁场数值模拟方法、装置、设备及介质对于沿y轴方向无限延伸的异常体构建异常体模型并对其剖分，考虑地下地质体磁化率主轴各向异性给其磁化率张量赋值；根据磁化率张量、空间域背景场磁场强度、空间域异常场磁场强度，得到空间域磁化强度；然后将空间域异常场磁位和空间域磁化强度满足的二维偏微分方程转为空间波数混合域一维常微分方程，将空间波数混合域异常场磁位满足的边值问题模型转化为等价的变分问题模型；通过求解变分问题模型，得到空间波数混合域异常场磁位、磁场强度后将其分别转换为空间域异常场磁位、磁场强度，迭代计算直至收敛。本发明考虑到地下地质体磁化率主轴各向异性，与实际地质情况更符合。
-
公开(公告)号：CN114065585B
公开(公告)日：2024-05-10
申请号：CN202111386036.4
申请日：2021-11-22
Applicant: 中南大学
IPC: G06F30/23 , G06F111/10
Abstract: 本发明公开了一种基于库伦规范的三维电性源数值模拟方法，在进行背景场计算时，提前计算了bessel积分中的层状系数，并在求解二次电磁场时采用傅里叶变换，将三维偏微分方程转化为波数下的一维常微分方程，最后采用紧算子修正求解过程中的电场。该方法基于二次场法，无需进行散度校正，利用这个优势，利用库伦规范将麦克斯韦方程组转换为库伦标量位矢量位的控制方程，在计算背景场时采用一种快速方法、二次场计算时利用傅里叶变换技术，极大地提高了三维电性源电磁场数值模拟的计算效率。本发明应用于数值模拟领域。
-
Patent Agency Ranking
公开(公告)号：CN115017782B
公开(公告)日：2022-12-16
申请号：CN202210941065.0
申请日：2022-08-08
Applicant: 中南大学
IPC: G06F30/23 , G06F17/13 , G01R29/08 , G06F111/10
Abstract: 本发明公开了一种考虑介质各向异性的三维天然源电磁场计算方法，采用基于二次场方法的偏微分方程的求解方法，引入洛伦兹规范，使最终得到的控制方程里仅有矢量位为未知量；在求解偏微分方程时，引入水平方向二维傅里叶变换，将三维偏微分方程转化为一维常微分方程，并利用伽辽金法最终合成三个五对角矩阵。本发明中最后合成的对角方程能用追赶法快速求解，计算效率高。使得一个复杂的各向异性介质三维天然源电磁场数值模拟转化为很多个常微分对角方程组的求解问题，简化了计算，又考虑了介质的各向异性。并提供了相应的计算设备，为大规模正演计算和反演成像提供了新的技术支持。
-
公开(公告)号：CN114036745A
公开(公告)日：2022-02-11
申请号：CN202111313171.6
申请日：2021-11-08
Applicant: 中南大学
IPC: G06F30/20 , G06F111/10
Abstract: 本发明公开了一种各向异性介质三维大地电磁正演方法及系统，该系统包括计算区域设定单元、计算区域剖分单元、导电率赋值单元、频率赋值单元、背景场计算单元、二次场计算单元、视电阻率和相位计算单元与输出单元。该系统在二次场计算过程中采用基于紧算子的迭代框架将二次电场逐步逼近真实的二次电场，得到满足误差条件的总电场后，再利用空间‑波数域磁场积分算法计算二次磁场。考虑了介质的各向异性，在正演时，采用了一个压缩紧算子，使得算法对于任意有损介质稳定收敛，能进行快速、准确的各向异性介质三维大地电磁数值模拟，能有效地提高各向异性介质三维大地电磁数据反演的精度和效率。本发明应用于导电率各向异性介质电磁数值模拟领域。
-
公开(公告)号：CN113656750B
公开(公告)日：2021-12-17
申请号：CN202111218085.7
申请日：2021-10-20
Applicant: 中南大学
Abstract: 基于空间波数域的强磁介质的磁感应强度计算方法，剖分包含异常体的初始棱柱体模型，得到一系列的节点，对每个节点进行磁化率赋值，磁化率为标量；设定高斯参数，计算离散偏移波数，计算空间域背景场磁场强度，获得磁化强度模型；将空间域异常场磁位和磁化强度满足的方程转为空间波数域一维常微分方程；结合设定的空间波数域异常场磁位的边界条件，将空间波数域异常场磁位满足的边值问题模型等价转化为变分问题模型并求解，得到空间波数域异常场磁位、空间波数域异常场磁场强度后将其转换为空间域异常场磁位及空间域异常场磁场强度，根据输出的空间域异常场磁场强度求得空间域磁感应强度。本发明能更准确的对强磁性介质进行磁场数值模拟。
-
公开(公告)号：CN113051779B
公开(公告)日：2021-08-10
申请号：CN202110596835.8
申请日：2021-05-31
Applicant: 中南大学
IPC: G06F30/20 , G06F111/10
Abstract: 本发明公开了一种三维直流电阻率法数值模拟方法，通过设置直流电法勘探区域和背景电阻率参数、三维异常体展布范围和电阻率参数、电流密度计算、二维离散傅里叶变换、波数域异常电位计算、二维反傅里叶变换、迭代收敛判断等步骤，实现了三维直流电法勘探高效、高精度的数值模拟。解决了目前直流电法勘探数值模拟中计算大规模模型时数据量大、存储要求高、计算时间慢，导致无法满足大规模直流电法数据精细反演的问题，有助于实现野外复杂地形或地下复杂结构下直流电法实测数据的精细反演和解释。
-
公开(公告)号：CN108197389A
公开(公告)日：2018-06-22
申请号：CN201810006227.5
申请日：2018-01-04
Applicant: 中南大学
IPC: G06F17/50
CPC classification number: G06F17/5009
Abstract: 本发明提供一种二维强磁性体磁场的快速、高精度数值模拟计算方法，其通过二维强磁性体模型表示、高斯参数设计、离散偏移波数计算、波数域加权系数计算、磁化强度计算、一维离散傅里叶反变换、空间域磁异常计算、迭代收敛判断等步骤，实现了二维强磁性体磁场数值模拟在效率和精度上的统一。本发明解决了目前强磁性体磁场数值模拟方法计算精度低，无法满足大规模强磁测数据精细反演成像的问题，有助于开展大规模强磁测数据二维磁化率精细反演成像、人机交互建模和解释的研究。
-
-
-
-
-
-
-
-
-
Search Results
26
records
(took 0.02 seconds)
