[发明专利]重力反演黄土层厚度分布的方法及装置

权利要求书

1.一种重力反演黄土层厚度分布的方法，其特征在于，包括：

采集重力勘探数据；

根据重力勘探数据确定布格重力异常数据和重力区域背景场；

根据布格重力异常数据和重力区域背景场，确定浅层黄土层剩余重力异常数据；

根据浅层黄土层剩余重力异常数据和重力勘探数据中的测点高程数据，进行三维重力密度反演计算，得到重力反演的三维密度数据；

根据重力反演的三维密度数据，确定黄土层厚度分布。

2.如权利要求1所述的方法，其特征在于，根据重力勘探数据确定布格重力异常数据，包括：

计算重力勘探数据的布格重力异常，中间层密度选取第一设定值并进行地形改正计算，获得布格重力异常数据；其中，在进行地形改正计算时不进行变密度改正；所述第一设定值大于等于2.1g/cm³。

3.如权利要求1所述的方法，其特征在于，根据重力勘探数据确定重力区域背景场，包括：

利用向上延拓方法取延拓高度为第二设定值，计算重力区域背景场；所述第二设定值小于等于1km。

4.如权利要求1所述的方法，其特征在于，确定浅层黄土层剩余重力异常数据之前，还包括：对布格重力异常数据进行网格化；

根据布格重力异常数据和重力区域背景场，确定浅层黄土层剩余重力异常数据，包括：

将网格化的布格重力异常数据减去重力区域背景场，得到浅层黄土层剩余重力异常网格文件；

进行三维重力密度反演计算之前，还包括：对测点高程数据进行网格化，得到高程数据网格文件；

根据浅层黄土层剩余重力异常数据和重力勘探数据中的测点高程数据，进行三维重力密度反演计算，包括：

根据浅层黄土层剩余重力异常网格文件和高程数据网格文件，进行三维重力密度反演计算。

5.如权利要求4所述的方法，其特征在于，根据浅层黄土层剩余重力异常网格文件和高程数据网格文件，进行三维重力密度反演计算，包括：

根据浅层黄土层剩余重力异常网格文件和高程数据网格文件，在设定反演密度上下限差、设定三维空间纵向剖分和设定重力反演数据拟合限差的条件下，以已知钻井黄土层底界深度作为井约束，进行三维重力密度反演计算。

6.如权利要求4所述的方法，其特征在于，根据重力反演的三维密度数据，确定黄土层厚度分布，包括：

从重力反演的三维密度数据中提取密度小于或等于第三设定值的网格节点，将提取的网格节点中密度最接近第三设定值的网格节点的纵向深度确定为平面上网格节点的黄土层厚度值；

将顶层网格节点中密度大于第三设定值的网格节点的黄土层厚度值确定为0；

根据平面上各网格节点的黄土层厚度值，确定平面上黄土层厚度分布。

7.如权利要求6所述的方法，其特征在于，根据各网格节点的黄土层厚度值，确定黄土层厚度分布，包括：

对平面上各网格节点的黄土层厚度值进行等值线成图，得到黄土层厚度分布等值线图。

8.一种重力反演黄土层厚度分布的装置，其特征在于，包括：

数据采集模块，用于采集重力勘探数据；

布格重力异常数据确定模块，用于根据重力勘探数据确定布格重力异常数据；

重力区域背景场确定模块，根据重力勘探数据确定重力区域背景场；

剩余重力异常数据确定模块，用于根据布格重力异常数据和重力区域背景场，确定浅层黄土层剩余重力异常数据；

重力反演模块，用于根据浅层黄土层剩余重力异常数据和重力勘探数据中的测点高程数据，进行三维重力密度反演计算，得到重力反演的三维密度数据；

黄土层厚度确定模块，用于根据重力反演的三维密度数据，确定黄土层厚度分布。

9.如权利要求8所述的装置，其特征在于，布格重力异常数据确定模块，具体用于：

计算重力勘探数据的布格重力异常，中间层密度选取第一设定值并进行地形改正计算，获得布格重力异常数据；其中，在进行地形改正计算时不进行变密度改正。