一种放射性地球物理测量数据融合方法与流程

本发明属于铀矿勘查技术领域，具体涉及一种放射性地球物理测量数据融合方法。

背景技术：

在铀矿勘查中，放射性地球物理测量主要包括伽马能谱测量、氡及其子体测量、β测量、中子测量、x荧光测量等等，其中伽马能谱测量中铀含量计算是利用氡子体特征能量峰计算获取，所以伽马能谱测量、氡及其子体测量所获得的铀含量和土壤氡浓度数据受氡的影响较大，氡主要指土壤或地层中所含²²²rn，由于氡的季度性变化以及不同年度不同仪器校准等因素产生误差等原因，在进行不同年度或季度的铀含量或土壤氡数据处理时，往往涉及到不同年度或季度的数据的融合时整体的不一致性和不协调性，严重影响了数据成果图件的效果和解释。在现有的测量规范中，没有涉及到该方面的问题，也没有相应数据处理方法，因此，开展相关的数据融合方法研究非常必要。

技术实现要素：

本发明的目的是提供一种放射性地球物理测量数据融合方法，能够有效解决在不同年度或不同季度中放射性地球物理测量数据整体不一致性，实现放射性地球物理测量数据无缝融合。

本发明所采用的技术方案是：

一种放射性地球物理测量数据融合方法，该方法包括以下步骤：

步骤一：收集新、老放射性地球物理原始实测数据；

步骤二：对新、老放射性地球物理原始实测数据按地质单元进行分类；

步骤三：对新、老放射性地球物理原始实测数据进行预处理；

步骤四：第一次误差分析；

步骤五：线性修正；

步骤六：第二次误差分析；

步骤七：非线性修正；

步骤八：第三次误差分析；

步骤九：将修正后的新、老数据合并融合和网格化处理；

步骤十：进行滤波分析处理、绘图。

如上所述步骤一：收集新、老放射性地球物理原始实测数据；包括：收集工作区不同年度或不同季度放射性地球物理原始实测数据，相对较早的测量数据记为老数据，将对较晚测量的数据记为新数据。

如上所述步骤二：对新、老放射性地球物理原始实测数据按地质单元进行分类，包括：将步骤一中经过预处理后的新、老放射性地球物理实测数据整理，分别形成含地理坐标和对应的放射性地球物理数据值的新、老二维空间散点数据文件，并将该数据按不同地质单元进行分类。

如上所述步骤三：对新、老放射性地球物理原始实测数据进行预处理，包括：

针对步骤二中按不同地质单元分类的新、老放射性地球物理空间散点数据，进行预处理，剔除单点异常，然后进行常规的统计噪声的降噪处理，再开展逐步剔除异常处理；

逐步剔除异常处理方法如下：对同一个地质单元新、老放射性地球物理数据，分别逐步删除大于等于“平均值+3倍均方差”的放射性地球物理异常数据，直至所有删除异常后的放射性地球物理数据值小于删除异常后统计的“平均值+3倍均方差”，然后，分别统计删除异常后的新、老放射性地球物理数据平均值和均方差，其它地质单元新、老数据逐步剔除异常处理方法相同。

如上所述步骤四：第一次误差分析，包括：对步骤三中剔除异常的相同地质单元实测的不同年度或季度放射性地球物理数据进行整体统计分析，将统计后的平均值进行整体相对误差识别，整体相对误差在大于-10％且小于10％的范围内，则无需校正；整体相对误差大于等于10％或小于等于-10％时，需对相对实测时间较早的老数据进行校正。

如上所述步骤五：线性修正，包括：根据步骤四分析结果，对整体相对误差大于等于10％或小于等于-10％的土壤氡浓度数据进行线性修正；

所述线性修正融合方法按公式(1)进行修正

xri＝k，ri(1)

公式(1)中，ri为统计地质单元每个测点的实测放射性地球物理老数据值；xri为修正后的每个测点实测放射性地球物理老数据值，k为线性修正参数；线性修正参数k的确定方法如下：对于老数据所对应的同一地层内含有新数据，则将步骤三中统计删除异常后的老数据平均值除以新数据平均值所得的比值，即为线性修正参数k；

若老数据所对应的同一地层内没有新数据，则选取老数据测点总量5％进行重复测量，然后按照步骤三中的方法对重复数据和对应测点的老数据进行逐步删除异常，然后、分别统计重复数据和对应测点的老数据的平均值，将重复数据的平均值除以对应测点的老数据平均值获得比值即为线性修正参数k值，然后，将k带入公式(1)，利用公式(1)对老数据进行修正。

如上所述步骤六：第二次误差分析，包括：对步骤五中经过线性校正后的老数据和新数据进行整体相对误差分析；整体相对误差大于-10％且小于10％，且线性校正后的老数据和新数据合并后，经过网格化处理的栅格数据没有跳跃性变化，则进行步骤九数据融合处理，否则，按步骤七进行非线性修正。

如上所述步骤七：非线性修正，包括：对步骤六中融合处理效果不明显时，则对数据非线性修正；

所述非线性修正步骤如下：

步骤(7.1)、利用离散余弦变换将步骤五中修正后的老数据和新数据，获得在离散余弦域内的新、老频谱数据；

所述步骤(7.1)中的离散余弦正变换所用的公式如下面式(2)所示：

式(2)中，{x(n),n＝1,2,…,n}为修正后的老数据或新数据，n为数据的行数，{y(k):k＝1,2,…n}为修正后的老数据或新数据的离散正余弦变换值，变换系数a(k)的定义如式(3)所示：

步骤(7.2)、在离散余弦域内，按照步骤三的处理方法对新、老频谱数据进行逐步剔除异常后，再分别求取新、老频谱的平均值和均方差，将新、老频谱数据的平均值加上三倍均方差后取整，分别记为e2n、e1n；

步骤(7.3)、将余弦域内老数据频谱大于e1n及新数据频谱大于e2n所有y(k)赋值为零，再将老频谱所有数据乘以步骤(7.2)中计算的新、老频谱的平均值的比值，最终形成新的新、老频谱数据；

步骤(7.4)、将步骤(7.3)中处理后的新的新、老频谱数据进行离散反余弦变换，将其从离散余弦域内转换为二维空间域；

所述步骤(7.4)中进行离散反余弦变换时所使用的公式如下面式(4)所示：

式(4)中各项含义与式(2)、式(3)中相同，获得修正后的新、老散点数据。

如上所述步骤八：第三次误差分析，包括：将步骤(7.4)中修正后的新老散点数据进行整体误差分析，整体对误差大于-10％且小于10％，且非线性校正后的新、老数据合并后经过网格化处理的栅格数据没有跳跃性变化，则进行步骤九数据融合处理，否则，返回步骤三开始重新循环处理。

如上所述步骤十：进行滤波分析处理、绘图，包括：对步骤九中网格化处理后的栅格数据，进行滤波分析处理，然后绘制成等值线平面图。

本发明的有益效果是：

本发明的一种放射性地球物理测量数据融合方法，采用线性和非线性相结合方法修正新、老土壤氡测量数据，能够有效解决不同年度或季度放射性实测的放射性地球物理数据整体的不一致性，实现了新、老放射性地球物理测量数据的无缝融合，为工程制图和成果解释提供了技术支撑。

附图说明

图1：本发明所述的一种放射性地球物理测量数据融合方法的流程图；

图2：为采用本发明处理前土壤氡浓度等值图；

图3：为采用本发明处理后的土壤氡浓度等值图。

具体实施方式

下面结合附图和具体实施例对本发明作进一步详细说明。

如图1至3所示，一种放射性地球物理测量数据融合方法，该方法包括以下步骤：

步骤一：收集新、老放射性地球物理原始实测数据；

步骤二：对新、老放射性地球物理原始实测数据按地质单元进行分类；

步骤三：对新、老放射性地球物理原始实测数据进行预处理；

步骤四：第一次误差分析；

步骤五：线性修正；

步骤六：第二次误差分析；

步骤七：非线性修正；

步骤八：第三次误差分析；

步骤九：将修正后的新、老数据合并融合和网格化处理

步骤十：进行滤波分析处理、绘图。

所述步骤一：收集新、老放射性地球物理原始实测数据；包括：收集工作区不同年度或不同季度放射性地球物理原始实测数据，相对较早的测量数据记为老数据，将对较晚测量的数据记为新数据。

所述步骤二：对新、老放射性地球物理原始实测数据按地质单元进行分类，包括：将步骤一中经过预处理后的新、老放射性地球物理实测数据整理，分别形成含地理坐标和对应的放射性地球物理数据值的新、老二维空间散点数据文件，并将该数据按不同地质单元进行分类。

所述步骤三：对新、老放射性地球物理原始实测数据进行预处理，包括：针对步骤二中按不同地质单元分类的新、老放射性地球物理空间散点数据，进行预处理，剔除单点异常，然后进行常规的统计噪声的降噪处理，再开展逐步剔除异常处理；

逐步剔除异常处理方法如下：对同一个地质单元新、老放射性地球物理数据，分别逐步删除大于等于“平均值+3倍均方差”的放射性地球物理异常数据，直至所有删除异常后的放射性地球物理数据值小于删除异常后统计的“平均值+3倍均方差”，然后，分别统计删除异常后的新、老放射性地球物理数据平均值和均方差，其它地质单元新、老数据逐步剔除异常处理方法相同。

所述步骤四：第一次误差分析，包括：对步骤三中剔除异常的相同地质单元实测的不同年度或季度放射性地球物理数据进行整体统计分析，将统计后的平均值进行整体相对误差识别，整体相对误差在大于-10％且小于10％的范围内，则无需校正；整体相对误差大于等于10％或小于等于-10％时，需对相对实测时间较早的老数据进行校正。

所述步骤五：线性修正，包括：根据步骤四分析结果，对整体相对误差大于等于10％或小于等于-10％的土壤氡浓度数据进行线性修正；

所述线性修正融合方法按公式(1)进行修正

xri＝k·ri(1)

公式(1)中，ri为统计地质单元每个测点的实测放射性地球物理老数据值；xri为修正后的每个测点实测放射性地球物理老数据值，k为线性修正参数；线性修正参数k的确定方法如下：对于老数据所对应的同一地层内含有新数据，则将步骤三中统计删除异常后的老数据平均值除以新数据平均值所得的比值，即为线性修正参数k；

若老数据所对应的同一地层内没有新数据，则选取老数据测点总量5％进行重复测量，然后按照步骤三中的方法对重复数据和对应测点的老数据进行逐步删除异常，然后、分别统计重复数据和对应测点的老数据的平均值，将重复数据的平均值除以对应测点的老数据平均值获得比值即为线性修正参数k值，然后，将k带入公式(1)，利用公式(1)对老数据进行修正。

所述步骤六：第二次误差分析，包括：对步骤五中经过线性校正后的老数据和新数据进行整体相对误差分析；整体相对误差大于-10％且小于10％，且线性校正后的老数据和新数据合并后，经过网格化处理的栅格数据没有跳跃性变化，则进行步骤九数据融合处理，否则，按步骤七进行非线性修正。

所述步骤七：非线性修正，包括：对步骤六中融合处理效果不明显时，则对数据非线性修正；

所述非线性修正步骤如下：

步骤(7.1)、利用离散余弦变换将步骤五中修正后的老数据和新数据，获得在离散余弦域内的新、老频谱数据；

所述步骤(7.1)中的离散余弦正变换所用的公式如下面式(2)所示：

式(2)中，{x(n),n＝1,2,…,n}为修正后的老数据或新数据，n为数据的行数，{y(k):k＝1,2,…n}为修正后的老数据或新数据的离散正余弦变换值，变换系数a(k)的定义如式(3)所示：

步骤(7.2)、在离散余弦域内，按照步骤三的处理方法对新、老频谱数据进行逐步剔除异常后，再分别求取新、老频谱的平均值和均方差，将新、老频谱数据的平均值加上三倍均方差后取整，分别记为e2n、e1n；

步骤(7.3)、将余弦域内老数据频谱大于e1n及新数据频谱大于e2n所有y(k)赋值为零，再将老频谱所有数据乘以步骤(7.2)中计算的新、老频谱的平均值的比值，最终形成新的新、老频谱数据；

步骤(7.4)、将步骤(7.3)中处理后的新的新、老频谱数据进行离散反余弦变换，将其从离散余弦域内转换为二维空间域；

所述步骤(7.4)中进行离散反余弦变换时所使用的公式如下面式(4)所示：

式(4)中各项含义与式(2)、式(3)中相同，获得修正后的新、老散点数据。

所述步骤八：第三次误差分析，包括：将步骤(7.4)中修正后的新老散点数据进行整体误差分析，整体对误差大于-10％且小于10％，且非线性校正后的新、老数据合并后经过网格化处理的栅格数据没有跳跃性变化，则进行步骤九数据融合处理，否则，返回步骤三开始重新循环处理。

所述步骤十：进行滤波分析处理、绘图，包括：对步骤九中网格化处理后的栅格数据，进行滤波分析处理，然后绘制成等值线平面图。

以下结合具体实施例对本发明作进一步描述：

步骤1、以相山火山盆地2014年和2016年实测的土壤氡测量数据为放射性地球物理原始实测数据，相对较早的2015年土壤氡测量数据记为老数据，将相对较晚的2016年测量的数据记为新数据；

步骤2、将步骤1中经过预处理后的新、老土壤氡实测数据整理，分别形成含地理坐标和对应的土壤氡数据值的新老空间散点数据文件，并将该数据按不同地质单元进行分类，因所选数据区域为同一地层，所以不用分类；

步骤3、针对步骤2中的新、老土壤氡空间散点数据，进行预处理，剔除单点异常，然后进行常规的统计噪声的降噪处理，再开展逐步剔除异常处理；

所述步骤3中，逐步剔除异常处理方法如下：对同一个地质单元新老土壤氡数据，分别逐步删除大于等于“平均值+3倍均方差”土壤氡异常数据，直至所有删除异常后的土壤氡浓度值小于删除异常后统计的“平均值+3倍均方差”，然后，统计删除异常后的土壤氡浓度平均值和均方差；

步骤4、对步骤3中剔除异常的土壤氡浓度进行整体统计分析，实际计算新、老土壤氡浓度整体相对误差为45％，明显大于10％，需对相对实测时间较早的2015年的老数据进行修正；

步骤5、根据步骤4分析结果，对2015年实测土壤氡浓度数据进行线性修正。

所述步骤5中，线性修正融合方法按公式(5)进行修正

xri＝k·ri(5)

公式(5)中，ri为统计地质单元2015年实测的每个测点的土壤氡浓度值；xri为修正后的每个测点的土壤氡值，k为线性修正参数。线性修正参数k的确定方法如下：将步骤(3)中统计删除异常后的2016年的土壤氡浓度平均值除以2015年的土壤氡浓度平均值所得的比值，即为线性修正参数k。计算获取线性修正参数k＝0.87，然后将修正系数k代入公式(5)中，利用公式(5)完成2015年老数据的修正。

步骤6、对步骤5中经过线性修正后的2015年度和2016年度的数据进行整体相对误差分析。统计计算整体对误差为8.9％，数据没有跳跃性变化，无需非线性修正，则按步骤7数据融合处理；

步骤7、将修正后的新老数据合并融合和网格化处理；

步骤8、对步骤7中网格化处理后的栅格数据，进行滤波分析处理，然后绘制成等值线平面图。

上面结合附图和实施例对本发明作了详细说明，但是本发明并不限于上述实施例，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。本发明中未作详细描述的内容均可以采用现有技术。