利索能及
我要发布
收藏
专利号: 2024103319883
申请人: 山东科技大学
专利类型:发明专利
专利状态:已下证
更新日期:2025-10-27
缴费截止日期: 暂无
联系人

摘要:

权利要求书:

1.基于深度三视图的重力异常智能反演方法,其特征在于,包括以下步骤:步骤1,构建重力异常样本集;

步骤2,基于U‑Net网络构建重力异常反演模型,利用训练集中的重力异常样本数据训练重力异常反演模型后,再用测试集中的重力异常样本数据测试重力异常反演模型,得到测试合格的重力异常反演模型;

步骤3,根据已勘探区域的重力异常数据,利用测试合格的重力异常反演模型计算得到已勘探区域中地下异常体的三视图深度数据后,利用模拟退火算法还原地下异常体的三维结构,结合重力异常反演模型反演计算的三视图深度数据优化模拟退火算法;

步骤4,将待还原区域的重力异常数据,利用重力异常反演模型计算得到待还原区域中地下异常体的三视图深度数据后,再利用优化后的模拟退火算法还原待还原区域中地下异常体的三维结构,得到待还原区域中地下异常体的三维立体图;

所述步骤1中,包括以下步骤:

步骤1.1,利用计算机程序设计软件在三维坐标系中随机生成多个地下异常体,所述地下异常体为长方体模型,分别设置各地下异常体的剩余密度,获取各地下异常体的质心坐标;

步骤1.2,针对各地下异常体分别进行正演计算,确定各地下异常体的重力异常数据,包括重力异常值和重力梯度异常值;

步骤1.3,分别针对各地下异常体,对地下异常体进行网格化处理,将地下异常体划分为多个立方体,获取地下异常体的三视图,将地下异常体主视图、俯视图和左视图中的立方体数量作为三视图深度数据并生成标签,再将同一地下异常体的标签与重力异常数据相匹配形成重力异常样本数据,得到多个重力异常样本数据;

步骤1.4,基于重力异常样本数据构建重力异常样本集,所述重力异常样本集包括训练集、测试集和验证集,按照预设的数量比将重力异常样本数据随机分配至训练集和测试集。

2.根据权利要求1所述的基于深度三视图的重力异常智能反演方法,其特征在于,所述步骤1.2中,利用重力异常及重力梯度异常正演计算公式分别计算各地下异常体的重力异常数据;

所述重力异常及重力梯度异常正演计算公式为:

      (1)

    (2)

    (3)

     (4)

     (5)

     (6)

     (7)

式中,为三维坐标系中 方向上离散网格点的序号,为三维坐标系中 方向上离散网格点的序号,为三维坐标系中 方向上离散网格点的序号; 为重力异常值; 为万有引力常量, 为长方体模型单元密度;、 、均为观测点的坐标; 为 方向上观测点与地下异常体中第 个离散网格点的距离差值, 为 方向上观测点与地下异常体中第个离散网格点的距离差值, 为 方向上观测点与地下异常体中第 个离散网格点的距离差值, 为观测点到地下异常体中离散网格点的距离, ;

为重力势能在 方向上的二阶偏导数, 为重力势能在 方向上的二阶偏导数,为重力势能在 方向上的二阶偏导数, 为重力势能在 方向上的二阶偏导数,为重力势能在 方向上的二阶偏导数, 为重力势能在 方向上的二阶偏导数。

3.根据权利要求1所述的基于深度三视图的重力异常智能反演方法,其特征在于,所述步骤2中,包括以下步骤:步骤2.1,基于U‑Net网络构建重力异常反演模型;

步骤2.2,设置重力异常反演模型的损失函数、学习率和损失函数预设值;

步骤2.3,从训练集中随机选取重力异常样本数据输入至重力异常反演模型中,利用重力异常反演模型根据重力异常样本数据中的重力异常值和重力梯度异常值计算得到三视图深度数据,并与重力异常样本数据标签中的三视图深度数据进行对比,计算损失函数值并与损失函数预设值进行对比;

步骤2.4,若训练过程中损失函数值小于损失函数预设值,则结束对重力异常反演模型的训练,得到训练后的重力异常反演模型并进入步骤2.5中,否则,则返回步骤2.3中继续训练重力异常反演模型;

步骤2.5,从测试集中随机选取重力异常样本数据输入至训练后的重力异常反演模型中,利用训练后的重力异常反演模型根据重力异常样本数据中的重力异常值和重力梯度异常值计算三视图深度数据,并与重力异常样本数据标签中的三视图深度数据进行对比,计算损失函数值并与损失函数预设值进行对比;

步骤2.6,若测试过程中损失函数值小于损失函数预设值,则结束测试过程,否则,则调整重力异常反演模型的学习率,返回至步骤2.3中。

4.根据权利要求3所述的基于深度三视图的重力异常智能反演方法,其特征在于,基于U‑Net网络构建重力异常反演模型,所述重力异常反演模型包括编码部分和解码部分,编码部分和解码部分之间跳跃连接,编码部分包括输入层和编码器,解码部分包括解码器和输出层,输入层、编码器、解码器和输出层依次连接;

所述输入层用于根据输入的重力异常样本数据获取特征图;

所述编码器内设置有四个依次连接的反卷积块,各反卷积块的结构相同,包括依次连接的转置卷积层、批归一化层和激活函数层,其中,转置卷积层的卷积核尺寸设置为3×3,转置卷积运算的步长设置为2、填充设置为1;

所述解码器内设置有两个依次连接的卷积块,各卷积块的结构相同,包括依次连接的标准卷积层、批归一化层和激活函数层,其中,标准卷积层通过双线性插值上采样调整输入特征图的尺寸,标准卷积层的卷积核尺寸设置为3×3,标准卷积运算的步长设置为1、填充设置为1;

所述输出层内设置有两个依次连接的卷积块;

所述激活函数层中的激活函数设置为ReLU函数;

所述跳跃连接用于将转置卷积层的输出结果与标准卷积层上采样所获取的输入特征图相合并。

5.根据权利要求4所述的基于深度三视图的重力异常智能反演方法,其特征在于,所述损失函数为均方误差损失函数,如公式(8)所示:                (8)

式中, 为重力异常反演模型的损失函数值;为输入重力异常样本数据的序号,为输入重力异常样本数据的总数; 为第 个重力异常样本数据标签中的三视图深度数据, 为重力异常反演模型预测的第 个重力异常样本数据的三视图深度数据。

6.根据权利要求1所述的基于深度三视图的重力异常智能反演方法,其特征在于,所述步骤3中,包括以下步骤:步骤3.1,获取已勘探区域的重力异常数据,对已勘探区域进行网格化处理划分为多个立方体,设置待还原区域三视图深度数据的初始解;

步骤3.2,设置模拟退火算法的能量函数、冷却函数、初始温度、冷却温度、温度下降梯度和迭代总次数;

步骤3.3,将初始解作为当前解、初始温度作为当前温度,利用模拟退火算法进行迭代计算,进入步骤3.4中;

步骤3.4,在当前温度条件下,利用能量函数计算当前解的能量值后,随机生成得到待还原区域三视图深度数据的更新解,并利用能量函数计算更新解的能量值;

步骤3.5,对比更新解的能量值和当前解的能量值;

若更新解的能量值低于当前解的能量值,则利用更新解替代当前解,进入步骤3.6中,否则,则保持当前解不变,进入步骤3.6中;

步骤3.6,基于温度冷却计划按照预设的温度下降梯度降低当前温度值,并将调整后的温度值作为当前温度,利用冷却函数计算当前温度条件下的冷却函数值,进入步骤3.7中;

步骤3.7,判断当前温度是否低于冷却温度以及当前迭代计算次数是否达到预设的迭代总次数;当满足当前温度低于冷却温度或当前迭代计算次数已达到预设迭代总次数,则得到待还原区域三视图深度数据的最终解,停止迭代计算并进入步骤3.8中,否则,则返回步骤3.4中;

步骤3.8,输出待还原区域三视图深度数据的最终解,得到待还原区域中地下异常体的三维立体图。

7.根据权利要求6所述的基于深度三视图的重力异常智能反演方法,其特征在于,所述能量函数用于计算模拟退火算法所还原三视图深度数据与重力异常反演模型所计算三视图深度数据之间的均方误差,如公式(9)所示:        (9)

式中, 为能量值;为三视图深度数据的序号, 为三视图深度数据的总数;

为重力异常反演模型所计算的三视图深度数据, 为模拟退火算法所还原的三视图深度数据; 为重力异常反演模型所计算三视图深度数据的平均值, 为重力异常反演模型所计算三视图深度数据的方差; 为模拟退火算法所还原三视图深度数据的平均值,为模拟退火算法所还原三视图深度数据的方差, 和 均为计算系数。

8.根据权利要求7所述的基于深度三视图的重力异常智能反演方法,其特征在于,所述冷却函数如公式(10)所示:                      (10)

式中, 为冷却函数值,用于表示地下异常体三维结构还原过程对于更新解的接受程度; 为能量函数值,为当前温度值,为自然对数的底数。