公开(公告)号：CN114065528B

公开(公告)日：2024-08-20

申请号：CN202111373803.8

申请日：2021-11-19

Applicant: 中南大学

Inventor： 王旭龙 ,  柳建新 ,  范平阳 ,  郭荣文

IPC: G06F30/20 ,  G06F17/15 ,  G06F111/10 ,  G06F119/14

Abstract: 本申请涉及一种二度体重力异常快速数值模拟方法、装置和计算机设备。所述方法包括：通过根据待探测的目标体的分布，构建二度体模型，同时将目标体进行网格剖分，得到目标体密度分布模型，其中，目标体密度分布模型的网格间隔与二度体模型一致；在此基础上，给定二度体重力异常公式，离散后计算核函数单元积分系数，最后通过一维快速离散卷积得到某一层模型的重力异常，通过累加可得到整个模型的重力异常。本发明对目标区域和待探测的目标体分别构建了网格模型，在计算重力异常时只需对目标体网格进行积分，极大地降低了内存需求，减少了积分系数计算，在保证原有精度不变的前提下提高了计算效率。

公开(公告)号：CN114065511B

公开(公告)日：2024-08-13

申请号：CN202111350961.1

申请日：2021-11-15

Applicant: 中南大学

Inventor： 李健 ,  柳建新 ,  郭荣文 ,  王永斐 ,  王旭龙 ,  袁洋

IPC: G06F30/20 ,  G06F111/10

Abstract: 起伏地形下大地电磁二维正演数值模拟方法、装置、设备及介质，包括：构建内部包含勘探目标的二维矩形模型；对二维矩形模型进行网格剖分并进行电导率赋值，得到刻画任意电导率分布的二维异常体模型；根据频率参数和二维异常体模型，构建关于观测点上待求解的场的积分表达式；根据积分表达式，加载一维边界条件以及二维大地电磁场源，得到场的控制方程；将二维异常体模型沿Z方向划分为多层子区域，在各层子区域分别沿水平方向构建不同稀疏度的插值算子，所有插值算子合成后得到一个总的稀疏采样算子，基于场的控制方程求解得到修正后的场，进而计算TE模式和TM模式下的阻抗、视电阻率和相位。该方法能够满足大规模电磁数据精细、快速正演成像的需求。

公开(公告)号：CN113642189B

公开(公告)日：2023-09-19

申请号：CN202110978719.2

申请日：2021-08-25

Applicant: 中南大学

Inventor： 王旭龙 ,  柳建新 ,  李健 ,  郭友军 ,  王永斐 ,  范平阳

IPC: G06F30/20 ,  G06F17/14 ,  G06F17/16 ,  G06F111/10 ,  G06F119/02

Abstract: 本申请涉及一种基于积分解的重力梯度张量快速数值模拟方法、装置、计算机设备和存储介质。所述方法包括：构建包含复杂地质体的二度体模型将地下模拟区域均匀剖分成多个规则单元，并根据二度体的截面形状和密度分布设定每个单元的密度；给出二度体重力梯度张量满足的卷积公式，根据上述构建的模型，对卷积公式进行离散，并计算核函数单元积分系数，采用二维离散傅里叶变换实现核函数单元积分系数矩阵和密度的快速相乘；最后，运用二维傅里叶反变换计算得到整个二维剖面的重力梯度。该方法引入二维离散傅里叶变换将重力梯度张量实现了核函数系数矩阵与密度的快速乘积，大大提高了正演的计算效率和精度。

公开(公告)号：CN113466954B

公开(公告)日：2022-05-06

申请号：CN202110814507.0

申请日：2021-07-19

Applicant: 中南大学 ,  中国人民解放军93114部队

Inventor： 陈杭 ,  任利兵 ,  王旭龙 ,  刘利 ,  柳建新 ,  魏登武 ,  雷文太 ,  杨刚强 ,  郭荣文

IPC: G01V3/38

Abstract: 本发明提供了一种自反馈类正则化反演方法，步骤S1：将平面磁异常数据进行化极并成图；步骤S2：根据化极后的磁异常平面图进行圈定，并对圈定数据进行插值；步骤S3：寻找极大值的平面位置，将此平面位置作为反演的初始位置；步骤S4：开始反演迭代过程，首先给初始模型赋值，然后对模型参数进行迭代更新，同时计算拟合差；迭代完毕后，得到输出反演的结果；步骤S5：得到反演模型参数以及吻合度；吻合度满足要求，则输出结果；吻合度不满足，则返回步骤S2，重新圈定数据，直到输出吻合度满足要求为止。本发明的新方法有效避免观测数据中噪声的影响，使反演过程更加稳健，更容易跳出局部极值，得出的反演结果更加接近实际位置。

公开(公告)号：CN114065511A

公开(公告)日：2022-02-18

申请号：CN202111350961.1

申请日：2021-11-15

Applicant: 中南大学

Inventor： 李健 ,  柳建新 ,  郭荣文 ,  王永斐 ,  王旭龙 ,  袁洋

IPC: G06F30/20 ,  G06F111/10

Abstract: 起伏地形下大地电磁二维正演数值模拟方法、装置、设备及介质，包括：构建内部包含勘探目标的二维矩形模型；对二维矩形模型进行网格剖分并进行电导率赋值，得到刻画任意电导率分布的二维异常体模型；根据频率参数和二维异常体模型，构建关于观测点上待求解的场的积分表达式；根据积分表达式，加载一维边界条件以及二维大地电磁场源，得到场的控制方程；将二维异常体模型沿Z方向划分为多层子区域，在各层子区域分别沿水平方向构建不同稀疏度的插值算子，所有插值算子合成后得到一个总的稀疏采样算子，基于场的控制方程求解得到修正后的场，进而计算TE模式和TM模式下的阻抗、视电阻率和相位。该方法能够满足大规模电磁数据精细、快速正演成像的需求。

公开(公告)号：CN113656976A

公开(公告)日：2021-11-16

申请号：CN202110978629.3

申请日：2021-08-25

Applicant: 中南大学

Inventor： 王旭龙 ,  柳建新 ,  范平阳 ,  郭友军 ,  王永斐

IPC: G06F30/20 ,  G06F17/14 ,  G06F17/16 ,  G06F111/10

Abstract: 本申请涉及一种二维磁梯度张量快速数值模拟方法、装置、计算机设备和存储介质。所述方法包括：通过模型构建、磁化强度计算、磁梯度张量计算公式离散、核函数单元积分系数计算、核函数单元积分系数矩阵和磁化强度的快速相乘、二维网格单元的磁梯度张量计算等步骤实现了二维磁梯度张量快速数值模拟模拟。本发明解决了现有二维磁梯度张量正演方法计算精度和计算效率低、无法实现复杂地质体的高效、精细化磁法勘探的问题。

公开(公告)号：CN113420487A

公开(公告)日：2021-09-21

申请号：CN202110978292.6

申请日：2021-08-25

Applicant: 中南大学

Inventor： 王旭龙 ,  柳卓 ,  柳建新 ,  郭荣文 ,  李健

IPC: G06F30/23 ,  G06F17/15 ,  G06F111/10 ,  G06F119/14

Abstract: 本申请涉及一种三维重力梯度张量数值模拟方法、装置、计算机设备和存储介质。所述方法包括：通过获取待计算区域在三维坐标系中的位置信息，将待计算区域沿x、y、z方向分别进行等间隔剖分，形成多个长方体的网格单元，根据网格剖分和预先获取的核函数三维积分表达式得到核函数单元积分的解析表达式，根据核函数单元积分的解析表达式得到核函数单元系数矩阵，其中核函数单元系数矩阵为Toeplitz矩阵，通过三维离散快速傅里叶变换进行核函数单元系数矩阵和单元体密度值三维离散卷积的快速相乘算法，得到待计算区域的重力梯度张量异常值。本发明通过快速相乘方法，减小了系数矩阵存储和计算问题，大幅提高了正演的精度和效率。

公开(公告)号：CN111967169B

公开(公告)日：2021-01-01

申请号：CN202011127944.7

申请日：2020-10-21

Applicant: 中南大学

Inventor： 王旭龙 ,  赵东东 ,  张钱江 ,  强健科 ,  李昆 ,  陈轻蕊

IPC: G06F30/20 ,  G06F17/14 ,  G06F111/10

Abstract: 本发明提出了一种二度体重力异常积分解数值模拟方法、装置、计算机设备和存储介质。本发明通过重力模型表示、波数域重力异常公式推导、波数域非均匀采样、空间域重力异常场计算等步骤，利用非均匀采样快速傅里叶变换实现了重力异常场高效、高精度数值模拟。在本发明中，将非均匀采样傅里叶变换应用于重力异常模拟，充分融合了快速傅里叶变换的效率优势和高斯傅里叶变换的精度优势，有效兼顾了重力异常数值模拟的计算精度与计算效率，解决了波数域重力异常数值模拟只适合规则测线、数值模拟方法不能同时兼顾计算效率和计算精度的问题。

公开(公告)号：CN113656976B

公开(公告)日：2023-09-01

申请号：CN202110978629.3

申请日：2021-08-25

Applicant: 中南大学

Inventor： 王旭龙 ,  柳建新 ,  范平阳 ,  郭友军 ,  王永斐

IPC: G06F30/20 ,  G06F17/14 ,  G06F17/16 ,  G06F111/10

Abstract: 本申请涉及一种二维磁梯度张量快速数值模拟方法、装置、计算机设备和存储介质。所述方法包括：通过模型构建、磁化强度计算、磁梯度张量计算公式离散、核函数单元积分系数计算、核函数单元积分系数矩阵和磁化强度的快速相乘、二维网格单元的磁梯度张量计算等步骤实现了二维磁梯度张量快速数值模拟模拟。本发明解决了现有二维磁梯度张量正演方法计算精度和计算效率低、无法实现复杂地质体的高效、精细化磁法勘探的问题。

公开(公告)号：CN112287534B

公开(公告)日：2022-05-13

申请号：CN202011128695.3

申请日：2020-10-21

Applicant: 中南大学

Inventor： 赵东东 ,  王旭龙 ,  张钱江 ,  强建科 ,  李昆 ,  范平阳

IPC: G06F30/20 ,  G06F17/14

Abstract: 本申请涉及一种基于NUFFT的二维磁异常快速正演模拟方法、装置、计算机设备和存储介质。本申请通过复杂磁性体模型表示、磁化强度计算、波数域磁位公式推导、一维形函数积分、非均匀波数选取、波数域磁异常场计算、空间域磁异常场计算等主要步骤，利用非均匀快速傅里叶变换和形函数积分实现了二维磁异常场高效、高精度正演模拟。在本申请中，充分利用非均匀快速傅里叶变换的灵活性和基于二次插值形函数积分的高精度优势，解决了波数域磁异常正演模拟局限于规则均匀网格，数值模拟不能同时兼顾计算精度和效率的问题。