一种根据布格重力异常图像提取隐伏断层的方法，该方法首先根据向上延拓的方法对布格重力异常剩余图像进行降噪预处理；再构建U型卷积网络对布格重力异常剩余图像进行训练；其次创建交叉熵损失函数，然后标注布格重力异常剩余图像的前后背景；再对布格重力异常剩余图像中单样本进行数据增强；接着根据快速并行细化算法对布格重力异常剩余图像中表现隐伏断层空间展布的特征体进行骨架提取；再根据最小二乘方法对布格重力异常剩余图像中的骨架进行线性拟合。该方法能够对降噪处理后的布格重力异常剩余图像采用U型卷积网络进行训练，标注布格重力异常剩余图像的前后背景进而提取抽象线性走向，生成线性构造，来减轻识别人员的工作量。
1.一种根据布格重力异常图像提取隐伏断层的方法，其特征在于至少包括如下步骤：(1)根据向上延拓的方法对布格重力异常剩余图像进行降噪预处理；
(2)构建U型卷积网络，并且根据U型卷积网络对布格重力异常剩余图像进行训练；
(3)创建交叉熵损失函数，并且利用交叉熵损失函数表示U型卷积网络实际输出概率与目标输出概率之间的相似度Loss，其计算公式为：式中， 为模型预测样本为正例的概率，y为样本标签，正例样本取值为1，负例样本取值为0；
(4)标注布格重力异常剩余图像的前后背景，将布格重力异常剩余图像前景标注为0，背景标注为1；
(5)对布格重力异常剩余图像中单样本进行数据增强；
(6)根据快速并行细化算法对布格重力异常剩余图像中表现隐伏断层空间展布的特征体进行骨架提取，即通过对符合特定条件的像素进行腐蚀从而只展示骨架；
所述步骤(6)中像素腐蚀步骤为：若像素点满足特定条件，则将该像素点的像素变为0；
ii)该像素点的八连通区域内部像素的像素值从0到1变换次数仅为1；
iii)该像素点的八连通区域内部像素的像素值为1的个数大于2小于6；
iv)该像素点的四连通区域内部上、右、下三个像素中有一个像素的像素值为0；
v)该像素点的四连通区域内部左、下、右三个像素中有一个像素的像素值为0；
(7)根据最小二乘方法对布格重力异常剩余图像中的骨架进行线性拟合。
2.根据权利要求1所述的根据布格重力异常图像提取隐伏断层的方法，其特征在于：所述步骤(1)中向上延拓即为计算高度大于布格重力异常剩余图像的原本高度。
3.根据权利要求1所述的根据布格重力异常图像提取隐伏断层的方法，其特征在于：所述步骤(2)中U型卷积网络包括编码器和解码器，编码器包括修正线性单元Xij，i＝1,2,
3，……，n，j＝1,2,3，……，k、3×3的卷积层，位于最上层的修正线性单元X11通过3×3的卷积层与修正线性单元X12相连，修正线性单元X12通过3×3的卷积层与修正线性单元X13相连，以此类推，直至最上层最后一个修正线性单元X1k被连接；每层的最后一个修正线性单元Xik通过2×2最大池化与下一层第一个修正线性单元Xi1相连，解码器包括修正线性单元yi(k+j)，i＝1,2,3，……，n，j＝1,2,3，……，k、3×3的卷积层，同理，同一层的修正线性单元yi(k+1)通过3×3的卷积层与修正线性单元yi(k+2)相连，且每层的最后一个修正线性单元yi(k+k)通过2×2向上卷积与上一层第一个修正线性单元yi(k+1)相连，每层的最后一个修正线性单元Xik与该层第一个修正线性单元yi(k+1)相连。
4.根据权利要求1所述的根据布格重力异常图像提取隐伏断层的方法，其特征在于，所述步骤(5)中单样本通过如下方式进行数据增强：(5.1)扭曲图像；
一种根据布格重力异常图像提取隐伏断层的方法
技术领域
[0001] 本发明涉及一种根据布格重力异常图像提取隐伏断层的方法，属于重力异常图像空间数据挖掘领域。
背景技术
[0002] 我国南方的存在各类超大型热液矿床和海底喷流沉积矿床，目前通过重力勘探的方法测量与围岩有密度差异的地质体在其周围引起的重力异常﹐以确定这些地质体存在的空间位置﹑大小和形状﹐从而对工作地区的地质构造和矿产分布情况作出判断，同时布格重力异常是重力勘探的基础资料，目前从布格重力异常图像中的隐伏断层都靠目视识别，因此要求识别人员经过专业训练、具有专业知识和素养，且解译质量受解译人员对解译区域的熟悉程度和解译经验等多种因素影响，主观性极大，同时识别人员的目视识别效率低，无法满足大区域、大面积的断层识别需求。
发明内容
[0003] 为解决上述现有技术中存在的技术问题，本发明提供了一种根据布格重力异常图像提取隐伏断层的方法，该方法能够对降噪处理后的布格重力异常剩余图像采用U型卷积网络进行训练，标注布格重力异常剩余图像的前后背景进而提取抽象线性走向，生成线性构造自动解译结果，大大提高识别效率，有利于矿产勘探与成矿预测。
[0004] 实现本发明目的的技术方案为，一种根据布格重力异常图像提取隐伏断层的方法，至少包括如下步骤：
[0005] (1)根据向上延拓的方法对布格重力异常剩余图像进行降噪预处理；
[0006] (2)构建U型卷积网络，并且根据U型卷积网络对布格重力异常剩余图像进行训练；
[0007] (3)创建交叉熵损失函数，并且利用交叉熵损失函数表示U型卷积网络实际输出概率与目标输出概率之间的相似度Loss，其计算公式为：
[0008]
[0009] 式中， 为模型预测样本为正例的概率，y为样本标签，正例样本取值为1，负例样本取值为0；
[0010] (4)标注布格重力异常剩余图像的前后背景，将布格重力异常剩余图像前景标注为0，背景标注为1；
[0011] (5)对布格重力异常剩余图像中单样本进行数据增强；
[0012] (6)根据快速并行细化算法对布格重力异常剩余图像中表现隐伏断层空间展布的特征体进行骨架提取，即通过对符合特定条件的像素进行腐蚀从而只展示骨架；
[0013] (7)根据最小二乘方法对布格重力异常剩余图像中的骨架进行线性拟合。
[0014] 对上述技术方案的进一步改进为：所述步骤(1)中向上延拓即为计算高度大于布格重力异常剩余图像的原本高度。
[0015] 且所述步骤(2)中U型卷积网络包括编码器和解码器，编码器包括修正线性单元Xij，i＝1,2,3，……，n，j＝1,2,3，……，k、3×3的卷积层，位于最上层的修正线性单元X11通过3×3的卷积层与修正线性单元X12相连，修正线性单元X12通过3×3的卷积层与修正线性单元X13相连，以此类推，直至最上层最后一个修正线性单元X1k被连接；每层的最后一个修正线性单元Xik通过2×2最大池化与下一层第一个修正线性单元Xi1相连，解码器包括修正线性单元yi(k+j)，i＝1,2,3，……，n，j＝1,2,3，……，k、3×3的卷积层，同理，同一层的修正线性单元yi(k+1)通过3×3的卷积层与修正线性单元yi(k+2)相连，且每层的最后一个修正线性单元yi(k+n)通过2×2向上卷积与上一层第一个修正线性单元yi(k+1)相连，每层的最后一个修正线性单元Xik与该层第一个修正线性单元yi(k+1)相连。
[0016] 且所述步骤(5)中单样本通过如下方式进行数据增强：
[0017] (5.1)扭曲图像；
[0018] (5.2)以50％的概率旋转；
[0019] (5.3)以80％的概率水平翻转；
[0020] (5.4)以30％的概率垂直翻转；
[0021] (5.5)剪裁50％的区域；
[0022] (5.6)随机生成指定分辨率大小的图像；
[0023] (5.7)从管道里生成指定数量的图像。
[0024] 且所述步骤(6)中像素腐蚀步骤为：若像素点满足特定条件，则将该像素点的像素变为0；若像素点不满足特定条件，则不对该像素点进行处理；
[0025] 特定条件为包含以下任一条件的条件；
[0026] i)该像素点值为1；
[0027] ii)该像素点的八连通区域内部像素的像素值从0到1变换次数仅为1；
[0028] iii)该像素点的八连通区域内部像素的像素值为1的个数大于2小于6；
[0029] iv)该像素点的四连通区域内部上、右、下三个像素中有一个像素的像素值为0；
[0030] v)该像素点的四连通区域内部左、下、右三个像素中有一个像素的像素值为0。
[0031] 由上述技术方案可知：(1)本方法通过向上延拓的方法可以将布格重力异常剩余图像中的规模较小的局部异常进行压制，以突出规模较大的区域异常，以提取对应时期的构造信息。
[0032] (2)重力数据因为专业性和保密性，相较于其他数据较难获取，造成能训练的数据样本较少，本方法采用的U型卷积网络能够融合不同层次的特征，以此增加训练效果。
[0033] (3)本方法使用交叉熵损失函数，用于衡量网络模型输出值与真实值之间的差距，给模型优化指明方向。
[0034] (4)本方法对表现隐伏断层空间展布的特征体进行骨架提取，来生成线性构造，大大提高识别效率，有利于矿产勘探与成矿预测。
附图说明
[0035] 图1为本发明提供的一种根据布格重力异常图像提取隐伏断层的方法的流程图；
[0036] 图2为本发明的U型卷积网络的结构示意图；
[0037] 图3为本发明的骨架提取效果图；
具体实施方式
[0038] 下面结合附图1至3和实施例对本发明进行详细具体说明，本发明的内容不局限于以下实施例。
[0039] 参考图1，一种根据布格重力异常图像提取隐伏断层的方法，包括如下步骤：根据向上延拓的方法对布格重力异常剩余图像进行降噪预处理；
[0040] 重力异常场源自地下不同深度和不同类型源体的位场叠加，包含局部异常和区域异常两种类型，重力异常数据处理时，需要进行重力异常分离，为提取对应时期的构造信息，需要将浅层场源进行压制，将该时期对应高度的场源分离出来，而采用向上延拓的方法主要是突出规模较大的区域异常，而压制规模较小的局部异常，向上延拓的基本方法通常利用波松积分推算相应高度处的重力异常，为了减少计算量，在本实施例中将布格重力异常剩余图像的原本高度上延10km的场值结果作为剩余异常。
[0041] 参考图2，构建U型卷积网络，并且根据U型卷积网络对布格重力异常剩余图像进行训练；U型卷积网络包括编码器和解码器，编码器包括修正线性单元Xij，i＝1,2,3，……，n，j＝1,2,3，……，k、3×3的卷积层，位于最上层的修正线性单元X11通过3×3的卷积层与修正线性单元X12相连，修正线性单元X12通过3×3的卷积层与修正线性单元X13相连，以此类推，直至最上层最后一个修正线性单元X1k被连接；每层的最后一个修正线性单元Xik通过2×2最大池化与下一层第一个修正线性单元Xi1相连，解码器包括修正线性单元yi(k+j)，i＝1,2,3，……，n，j＝1,2,3，……，k、3×3的卷积层，同理，同一层的修正线性单元yi(k+1)通过3×3的卷积层与修正线性单元yi(k+2)相连，且每层的最后一个修正线性单元yi(k+n)通过2×2向上卷积与上一层第一个修正线性单元yi(k+1)相连，每层的最后一个修正线性单元Xik与该层第一个修正线性单元yi(k+1)相连，重力数据因为专业性和保密性，相较于其他数据较难获取，造成能训练的数据样本较少，本方法采用的U型卷积网络能够融合不同层次的特征，以此增加训练效果。
[0042] 创建交叉熵损失函数，并且利用交叉熵损失函数表示U型卷积网络实际输出概率与目标输出概率之间的相似度Loss，，其计算公式为：
[0043]
[0044] 式中，为模型预测样本为正例的概率，y为样本标签，正例样本取值为1，负例样本取值为0。
[0045] 通过使用交叉熵损失函数，用于衡量网络模型输出值与真实值之间的差距，给模型优化指明方向。
[0046] 标注布格重力异常剩余图像的前后背景，将布格重力异常剩余图像前景标注为0，背景标注为1；
[0047] 对布格重力异常剩余图像中单样本进行数据增强；单样本通过如下方式进行数据增强：
[0048] 扭曲图像；
[0049] 以50％的概率旋转；
[0050] 以80％的概率水平翻转；
[0051] 以30％的概率垂直翻转；
[0052] 剪裁50％的区域；
[0053] 随机生成指定分辨率大小的图像；
[0054] 从管道里生成指定数量的图像。
[0055] 单样本数据增强是采用预设的数据变换规则，在现有数据的基础上进行扩增的方法，增强包括几何操作类、颜色变换类。几何操作类包括翻转、旋转、剪裁、变形、缩放等各种操作。翻转和旋转适用于对方向不敏感的任务，比如图像分类；剪裁相对翻转和旋转，会改变图像的大小，缩放变形则会造成图像失真。几何变换类没有改变图像本身的内容，只是对图像某部分像素进行了重分布。颜色变换类则改变了图像内容，常见操作如噪音、模糊、颜色变换、擦除、填充等。数据增强的本质就是增强模型的泛化能力。
[0056] 参考图3，根据快速并行细化算法对布格重力异常剩余图像中表现隐伏断层空间展布的特征体进行骨架提取，即通过对符合特定条件的像素进行腐蚀从而只展示骨架，所述步骤中像素腐蚀步骤为：若像素点满足特定条件，则将该像素点的像素变为0；若像素点不满足特定条件，则不对该像素点进行处理；
[0057] 特定条件为包含以下任一条件的条件；
[0058] i)该像素点值为1；
[0059] ii)该像素点的八连通区域内部像素的像素值从0到1变换次数仅为1；八连通区域指的是从区域内每一象素出发，可通过八个方向，即上、下、左、右、左上、右上、左下、右下这八个方向的移动的组合，在不越出区域的前提下，到达区域内的任意象素。
[0060] iii)该像素点的八连通区域内部像素的像素值为1的个数大于2小于6；
[0061] iv)该像素点的四连通区域内部上、右、下三个像素中有一个像素的像素值为0，四连通区域指的是从区域内每一象素出发，可通过四个方向，即上、下、左、右这四个方向的移动的组合，在不越出区域的前提下，到达区域内的任意象素；
[0062] v)该像素点的四连通区域内部左、下、右三个像素中有一个像素的像素值为0。
[0063] 根据快速并行细化算法对布格重力异常剩余图像中表现隐伏断层空间展布的特征体进行骨架提取，即通过对符合特定条件的像素进行腐蚀从而只展示骨架，[0064] 根据最小二乘方法对布格重力异常剩余图像中的骨架进行线性拟合，最小二乘法(又称最小平方法)是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。
