[发明专利]一种热液型铀矿床基底界面反演方法

权利要求书

1.一种热液型铀矿床基底界面反演方法，其特征在于：包括以下步骤：

步骤一

在勘查区开展RTK测量，获得测点经纬度坐标和高程数据；

步骤二

在勘查区开展高精度重力测量，获得重力数据；

步骤三

采集勘查区不同地层、岩体岩石标本，对其进行密度测量，统计得到每一类地层或岩体岩石的平均密度参数；

步骤四

对重力数据进行高程校正、布格重力异常改正和纬度改正，得到测点重力异常△g值；

步骤五

对步骤一得到的测点高程数据进行网格化处理，对步骤四所得到的重力异常△g值进行数据平滑处理，消除高频干扰，得到光滑的重力异常曲线并对所有光滑处理后的△g值进行网格化处理；

步骤六

创建反演初始模型，并把模型分为两层，对两层模型进行密度赋值，创建反演目标界面常量网格、反演观测面选择步骤五得出的高程网格；

步骤七

利用帕克反演公式对目标界面进行反演，在反演过程中加入余弦滤波器，反演出密度界面网格；

步骤八

对步骤七得到的密度界面网格利用帕克正演公式进行正演计算，修改模型参数、余弦滤波器参数；

步骤九

重复步骤七、步骤八，直至拟合误差达到要求，获得最终基底密度界面网格。

2.根据权利要求1所述的热液型铀矿床基底界面反演方法，其特征在于：

步骤四中，所述布格重力异常改正，采用式(1)进行计算：

δ_gB＝(0.3086－0.419{σ}{h}_m) 式(1)

式中，

δ_gB为布格重力异常值；

σ为中间层密度取值，本实施例中，σ取2.67×10^-5g/cm³；

h为测量高程。

3.根据权利要求2所述的热液型铀矿床基底界面反演方法，其特征在于：

步骤七中，所述帕克反演方法采用式(2)进行计算：

式中，

为密度界面起伏的傅里叶变换；

△g为重力异常；

J为上下界面密度差值；

S为径向波数；

为每次迭代计算出的高程值累加值。

4.根据权利要求3所述的热液型铀矿床基底界面反演方法，其特征在于：

步骤九中，所述帕克正演方法采用式(3)进行计算：

式中，

为密度界面起伏的傅里叶变换。

5.根据权利要求1所述的热液型铀矿床基底界面反演方法，其特征在于：步骤二中，重力测量总精度高于0.05mgal，勘查区选择地形平坦的区域。

6.根据权利要求1所述的热液型铀矿床基底界面反演方法，其特征在于：步骤三中，把待测标本切割成5×5×5cm的规则正方体，利用天秤对其进行质量测量，求取出每块标本的密度，并对同类地层或岩体岩石密度的标本进行统计，计算出平均值作为平均密度参数。

7.根据权利要求1所述的热液型铀矿床基底界面反演方法，其特征在于：步骤五中，高程数据网格化处理后得到的高程网格文件、与重力异常△g值网格化处理得到的布格重力异常网格文件，二者规格相同；步骤六中，目标界面常量网格的规格与步骤五得到的高程网格文件、布格重力异常网格文件规格一致。

8.根据权利要求7所述的热液型铀矿床基底界面反演方法，其特征在于：步骤五得到的高程网格文件、布格重力异常网格文件，步骤六得到的目标界面常量网格，网格均为402×403。

9.根据权利要求1所述的热液型铀矿床基底界面反演方法，其特征在于：步骤六中，据步骤三中不同地层或岩体岩石平均密度参数对两层模型进行密度赋值；创建反演目标界面常量网格时网格常量为-1000。

10.根据权利要求1所述的热液型铀矿床基底界面反演方法，其特征在于：步骤七中，所述余弦滤波器选择低通滤波，其上限截止波长为0，下限截止波长为0.02。