[发明专利]一种放射性地球物理测量数据融合方法

权利要求书

1.一种放射性地球物理测量数据融合方法，其特征在于：该方法包括以下步骤：

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

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

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

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

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

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

XR_i＝k·R_i (1)

公式(1)中，R_i为统计地质单元每个测点的实测放射性地球物理老数据值；XR_i为修正后的每个测点实测放射性地球物理老数据值，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)、在离散余弦域内，按照步骤三的处理方法对新、老频谱数据进行逐步剔除异常后，再分别求取新、老频谱的平均值和均方差，将新、老频谱数据的平均值加上三倍均方差后取整，分别记为E_2n、E_1n；

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

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

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

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

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

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

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

2.根据权利要求1所述的一种放射性地球物理测量数据融合方法，其特征在于：所述步骤一：收集新、老放射性地球物理原始实测数据；包括：收集工作区不同年度或不同季度放射性地球物理原始实测数据，相对较早的测量数据记为老数据，将对较晚测量的数据记为新数据。

3.根据权利要求2所述的一种放射性地球物理测量数据融合方法，其特征在于：所述步骤二：对新、老放射性地球物理原始实测数据按地质单元进行分类，包括：将步骤一中经过预处理后的新、老放射性地球物理实测数据整理，分别形成含地理坐标和对应的放射性地球物理数据值的新、老二维空间散点数据文件，并将该数据按不同地质单元进行分类。