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

说明书

本发明属于热液岩型铀矿床重力数据处理技术领域，具体涉及一种热液型铀矿床基底界面反演方法。本发明包括以下步骤：获得勘查区测点经纬度坐标、高程数据、重力数据、以及每一类地层或岩体岩石的平均密度参数；计算测点重力异常△g值；对高程数据及重力异常△g值网格化处理后创建反演目标界面常量网格，对目标界面进行反复进行反演和正演，获得最终基底密度界面网格。本发明能够有效、真实地探测基底界面形态和空间位置。

技术领域

本发明属于热液岩型铀矿床重力数据处理技术领域，具体涉及一种热液型铀矿床基底界面反演方法。

背景技术

经长期铀矿找矿勘查及开采工作，目前近地表可利用的铀资源量已大幅度减少，新一轮的铀矿找矿工作目标逐渐转向矿田深部。在热液岩型铀矿勘探领域，不同类型岩石(如火山碎屑岩、凝灰岩、花岗斑岩和变质岩等)密度差异明显，而变质岩基底与上覆火山岩界面逐渐被关注，并可能成为有利于铀成矿的第二空间，因此基底界面的起伏及其深度信息显得尤为重要，对热液型铀矿床整体结构、深部构造和深部成矿环境的研究具有重要的现实意义。

在地球物理重力数据处理领域，重力数据处理和解释中准确推断铀矿床基底界面起伏形态、深度往往是地球物理勘查的难点。由于火山岩盖层与变质岩基底存在较大密度差异，常规的密度界面自动反演技术如帕克反演法可以自动确定重力界面起伏深度，但由于向下延拓因子的存在，导致了反演过程的不稳定性，且在地质构造复杂区域，原始数据存在误差干扰，向下延拓因子可能导致反演精度降低甚至使反演结果发散，从而最终影响基底界面埋深信息的准确性。

因此，亟需研制一种热液型铀矿床密度界面反演方法，以缓解向下延拓因子对反演结果的影响，在不降低反演精度的前提下提高反演算法的稳定性。从而准确对地下地质信息进行有效的提取，同时反映出基底界面起伏形态和深度的变化。

发明内容

本发明需要解决的技术问题为：提供一种热液型铀矿床基底界面反演方法，能够有效、真实地探测基底界面形态和空间位置。

本发明的技术方案如下所述：

一种热液型铀矿床基底界面反演方法，包括以下步骤：

步骤一

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

步骤二

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

步骤三

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

步骤四

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

步骤五

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

步骤六

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

步骤七

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

步骤八

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

步骤九

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