[发明专利]一种面向铀矿床定位的瞬时测氡数据异常提取方法

说明书

技术领域

本发明涉及一种面向铀矿床定位的瞬时测氡数据异常提取方法，特别是涉及一种通过分析瞬时测氡数据中由于覆盖层厚度对实测数据带来的不确定影响程度的面向铀矿床定位的瞬时测氡数据异常提取方法。

背景技术

氡是铀的直接衰变物，也是自然界唯一的放射性的气体，由镭衰变产生，镭又是由放射性元素铀衰变而来，铀起了一个氡的永久源作用，因此，在铀矿地质中，土壤和地表瞬时测氡是一种重要的铀矿床攻深找盲方法。大量矿床实例证明，深达几百米的铀矿石可以在地表形成明显的氡异常。

在数据处理和异常提取过程中，受到地质覆盖层岩性及其厚度和测量仪器的误差影响，一般对瞬时测氡数据只进行滤波或者小波变换等噪声滤除处理，数值较大的测量点即被视为异常点，毗邻的多个异常点形成异常区，通过绘制等值线进行分析应用。

不同地区的铀矿埋深是不同的，地质覆盖层的岩性、产状也不尽相同，地貌景观也有所差别，现有数据处理方法对此考虑不全面，只是简单的通过数值大小来判断异常区的位置。

因此亟需提供一种新型的面向铀矿床定位的瞬时测氡数据异常提取方法。

发明内容

本发明要解决的技术问题是提供一种在砂岩型铀矿床定位中获得更高的定位精度的面向铀矿床定位的瞬时测氡数据异常提取方法。

为解决上述技术问题，本发明一种面向铀矿床定位的瞬时测氡数据异常提取方法，依次包括以下步骤：

步骤1、确定待分析的瞬时测氡数据；

每个瞬时测氡布点采集一个测量值，所有测量值组成二位离散数据A；

步骤2、设定均值滤波器；

面向二位离散数据A建立一个均值滤波器B，所述均值滤波器的窗口尺寸为奇数；

步骤3、对测量数据A采用均值滤波器B进行滤波处理，结果记录为C；

步骤4、对C每一点的数值逐点求平方，结果记录为D；

步骤5、对测量数据A的测量数值逐点求平方，并采均值滤波器B进行滤波处理，结果记录为E；

步骤6、对E和D逐点求差，并对差值逐点开平方，结果记录为F；

步骤7、用F数值逐点除以对应的C点数值，获得新的异常提取结果G。

均值滤波器的窗口尺寸为5。

受到地表覆盖层的影响，氡和氡的母体和子体元素在由地下向地表扩散的过程中，其扩散路径受到周围地质体、岩性、产状、厚度、含水性和温度等因素的影响会不断发生变化。与当前仅直接通过测量数值大小确定氡异常的方法相比，本发明考虑到这一影响因素，并利用这种不确定性来提取氡的异常分布，采取本发明具体的步骤及参数，可在砂岩型铀矿床定位中获得更高的定位精度。

具体实施方式

步骤1、确定待分析的瞬时测氡数据；

通常应用于地质矿床定位的瞬时测氡布点构成一个二维矩阵，布点的位置为了便于分析通常沿行和列方向等距，非等距情况下可以通过数据重采样获得，每个布点记录一个测量值，在实际处理中也可以被认为是一幅离散图像并记录为A；

步骤2、设定均值滤波器；

面向二位离散数据A建立一个均值滤波器B，所述均值滤波器的窗口尺寸为奇数；窗口尺寸的大小可根据实测数据和区域地质情况确定,当地质覆盖层比较均匀时，滤波器的窗口尺寸可以适当取大一些；反之，则滤波器的窗口尺寸可以适当取小一些；优选取值5；

步骤3、对测量数据A采用均值滤波器B进行滤波处理，结果记录为C；

步骤4、对C每一点的数值逐点求平方，结果记录为D；

步骤5、对测量数据A的测量数值逐点求平方，并采均值滤波器B进行滤波处理，结果记录为E；

步骤6、对E和D逐点求差，并对差值逐点开平方，结果记录为F；

步骤7、用F数值逐点除以对应的C点数值，获得新的异常提取结果G。