CN109975884B - 一种放射性地球物理测量数据融合方法 - Google Patents
一种放射性地球物理测量数据融合方法 Download PDFInfo
- Publication number
- CN109975884B CN109975884B CN201711457720.0A CN201711457720A CN109975884B CN 109975884 B CN109975884 B CN 109975884B CN 201711457720 A CN201711457720 A CN 201711457720A CN 109975884 B CN109975884 B CN 109975884B
- Authority
- CN
- China
- Prior art keywords
- data
- old
- new
- radioactive
- geophysical
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000002285 radioactive effect Effects 0.000 title claims abstract description 73
- 238000005259 measurement Methods 0.000 title claims abstract description 26
- 238000007500 overflow downdraw method Methods 0.000 title claims description 9
- 238000012937 correction Methods 0.000 claims abstract description 67
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000004458 analytical method Methods 0.000 claims abstract description 34
- 238000012545 processing Methods 0.000 claims abstract description 23
- 238000001914 filtration Methods 0.000 claims abstract description 11
- 238000007781 pre-processing Methods 0.000 claims abstract description 11
- 230000004927 fusion Effects 0.000 claims abstract description 10
- 238000001228 spectrum Methods 0.000 claims description 38
- 229910052704 radon Inorganic materials 0.000 claims description 31
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims description 31
- 239000002689 soil Substances 0.000 claims description 26
- 230000005856 abnormality Effects 0.000 claims description 19
- 238000012217 deletion Methods 0.000 claims description 13
- 230000037430 deletion Effects 0.000 claims description 13
- 238000003672 processing method Methods 0.000 claims description 10
- 230000002159 abnormal effect Effects 0.000 claims description 9
- 238000002441 X-ray diffraction Methods 0.000 claims description 8
- 230000009191 jumping Effects 0.000 claims description 7
- 238000007499 fusion processing Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 4
- 230000008030 elimination Effects 0.000 claims description 4
- 238000003379 elimination reaction Methods 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 4
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 claims description 3
- 238000004064 recycling Methods 0.000 claims description 3
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 claims description 3
- 229910052770 Uranium Inorganic materials 0.000 abstract description 6
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于铀矿勘查技术领域,具体涉及一种放射性地球物理测量数据融合方法,包括以下步骤:步骤一:收集新、老放射性地球物理原始实测数据;步骤二:对新、老放射性地球物理原始实测数据按地质单元进行分类;步骤三:对新、老放射性地球物理原始实测数据进行预处理;步骤四:第一次误差分析;步骤五:线性修正;步骤六:第二次误差分析;步骤七:非线性修正;步骤八:第三次误差分析;步骤九:将修正后的新、老数据合并融合和网格化处理;步骤十:进行滤波分析处理、绘图。本发明能够有效解决不同年度或季度放射性实测的放射性地球物理数据整体的不一致性,实现了新、老放射性地球物理测量数据的无缝融合。
Description
技术领域
本发明属于铀矿勘查技术领域,具体涉及一种放射性地球物理测量数据融合方法。
背景技术
在铀矿勘查中,放射性地球物理测量主要包括伽马能谱测量、氡及其子体测量、β测量、中子测量、X荧光测量等等,其中伽马能谱测量中铀含量计算是利用氡子体特征能量峰计算获取,所以伽马能谱测量、氡及其子体测量所获得的铀含量和土壤氡浓度数据受氡的影响较大,氡主要指土壤或地层中所含222Rn,由于氡的季度性变化以及不同年度不同仪器校准等因素产生误差等原因,在进行不同年度或季度的铀含量或土壤氡数据处理时,往往涉及到不同年度或季度的数据的融合时整体的不一致性和不协调性,严重影响了数据成果图件的效果和解释。在现有的测量规范中,没有涉及到该方面的问题,也没有相应数据处理方法,因此,开展相关的数据融合方法研究非常必要。
发明内容
本发明的目的是提供一种放射性地球物理测量数据融合方法,能够有效解决在不同年度或不同季度中放射性地球物理测量数据整体不一致性,实现放射性地球物理测量数据无缝融合。
本发明所采用的技术方案是:
一种放射性地球物理测量数据融合方法,该方法包括以下步骤:
步骤一:收集新、老放射性地球物理原始实测数据;
步骤二:对新、老放射性地球物理原始实测数据按地质单元进行分类;
步骤三:对新、老放射性地球物理原始实测数据进行预处理;
步骤四:第一次误差分析;
步骤五:线性修正;
步骤六:第二次误差分析;
步骤七:非线性修正;
步骤八:第三次误差分析;
步骤九:将修正后的新、老数据合并融合和网格化处理;
步骤十:进行滤波分析处理、绘图。
如上所述步骤一:收集新、老放射性地球物理原始实测数据;包括:收集工作区不同年度或不同季度放射性地球物理原始实测数据,相对较早的测量数据记为老数据,将对较晚测量的数据记为新数据。
如上所述步骤二:对新、老放射性地球物理原始实测数据按地质单元进行分类,包括:将步骤一中经过预处理后的新、老放射性地球物理实测数据整理,分别形成含地理坐标和对应的放射性地球物理数据值的新、老二维空间散点数据文件,并将该数据按不同地质单元进行分类。
如上所述步骤三:对新、老放射性地球物理原始实测数据进行预处理,包括:
针对步骤二中按不同地质单元分类的新、老放射性地球物理空间散点数据,进行预处理,剔除单点异常,然后进行常规的统计噪声的降噪处理,再开展逐步剔除异常处理;
逐步剔除异常处理方法如下:对同一个地质单元新、老放射性地球物理数据,分别逐步删除大于等于“平均值+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中网格化处理后的栅格数据,进行滤波分析处理,然后绘制成等值线平面图。
上面结合附图和实施例对本发明作了详细说明,但是本发明并不限于上述实施例,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。本发明中未作详细描述的内容均可以采用现有技术。
Claims (7)
1.一种放射性地球物理测量数据融合方法,其特征在于:该方法包括以下步骤:
步骤一:收集新、老放射性地球物理原始实测数据;
步骤二:对新、老放射性地球物理原始实测数据按地质单元进行分类;
步骤三:对新、老放射性地球物理原始实测数据进行预处理;
步骤四:第一次误差分析;
步骤五:线性修正;包括:根据步骤四分析结果,对整体相对误差大于等于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)中相同,获得修正后的新、老散点数据;
步骤八:第三次误差分析;
步骤九:将修正后的新、老数据合并融合和网格化处理
步骤十:进行滤波分析处理、绘图。
2.根据权利要求1所述的一种放射性地球物理测量数据融合方法,其特征在于:所述步骤一:收集新、老放射性地球物理原始实测数据;包括:收集工作区不同年度或不同季度放射性地球物理原始实测数据,相对较早的测量数据记为老数据,将对较晚测量的数据记为新数据。
3.根据权利要求2所述的一种放射性地球物理测量数据融合方法,其特征在于:所述步骤二:对新、老放射性地球物理原始实测数据按地质单元进行分类,包括:将步骤一中经过预处理后的新、老放射性地球物理实测数据整理,分别形成含地理坐标和对应的放射性地球物理数据值的新、老二维空间散点数据文件,并将该数据按不同地质单元进行分类。
4.根据权利要求3所述的一种放射性地球物理测量数据融合方法,其特征在于:所述步骤三:对新、老放射性地球物理原始实测数据进行预处理,包括:针对步骤二中按不同地质单元分类的新、老放射性地球物理空间散点数据,进行预处理,剔除单点异常,然后进行常规的统计噪声的降噪处理,再开展逐步剔除异常处理;
逐步剔除异常处理方法如下:对同一个地质单元新、老放射性地球物理数据,分别逐步删除大于等于“平均值+3倍均方差”的放射性地球物理异常数据,直至所有删除异常后的放射性地球物理数据值小于删除异常后统计的“平均值+3倍均方差”,然后,分别统计删除异常后的新、老放射性地球物理数据平均值和均方差,其它地质单元新、老数据逐步剔除异常处理方法相同。
5.根据权利要求4所述的一种放射性地球物理测量数据融合方法,其特征在于:所述步骤四:第一次误差分析,包括:对步骤三中剔除异常的相同地质单元实测的不同年度或季度放射性地球物理数据进行整体统计分析,将统计后的平均值进行整体相对误差识别,整体相对误差在大于-10%且小于10%的范围内,则无需校正;整体相对误差大于等于10%或小于等于-10%时,需对相对实测时间较早的老数据进行校正。
6.根据权利要求1所述的一种放射性地球物理测量数据融合方法,其特征在于:所述步骤八:第三次误差分析,包括:将步骤(7.4)中修正后的新老散点数据进行整体误差分析,整体对误差大于-10%且小于10%,且非线性校正后的新、老数据合并后经过网格化处理的栅格数据没有跳跃性变化,则进行步骤九数据融合处理,否则,返回步骤三开始重新循环处理。
7.根据权利要求1所述的一种放射性地球物理测量数据融合方法,其特征在于:所述步骤十:进行滤波分析处理、绘图,包括:对步骤九中网格化处理后的栅格数据,进行滤波分析处理,然后绘制成等值线平面图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711457720.0A CN109975884B (zh) | 2017-12-28 | 2017-12-28 | 一种放射性地球物理测量数据融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711457720.0A CN109975884B (zh) | 2017-12-28 | 2017-12-28 | 一种放射性地球物理测量数据融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109975884A CN109975884A (zh) | 2019-07-05 |
CN109975884B true CN109975884B (zh) | 2020-10-23 |
Family
ID=67074458
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711457720.0A Active CN109975884B (zh) | 2017-12-28 | 2017-12-28 | 一种放射性地球物理测量数据融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109975884B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113050187B (zh) * | 2019-12-26 | 2024-03-01 | 中国石油天然气集团有限公司 | 滤波方法及装置、计算机设备及计算机可读存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4964089A (en) * | 1989-12-22 | 1990-10-16 | Conoco Inc. | Method for derivation of internal velocities from post-migration parts |
EP0477806A2 (en) * | 1990-09-25 | 1992-04-01 | Honeywell Inc. | Method and apparatus for producing an accurate time base |
JPH09245074A (ja) * | 1996-03-11 | 1997-09-19 | Nec Corp | 回路接続情報の比較装置 |
CN1683939A (zh) * | 2004-04-13 | 2005-10-19 | 西门子公司 | 在磁共振断层造影中融合加权成像的运动校正的多点方法 |
US7489928B2 (en) * | 2002-12-31 | 2009-02-10 | Smith Brian K | Adaptive RF link failure handler |
CN103499836A (zh) * | 2013-07-18 | 2014-01-08 | 中国石油集团西北地质研究所 | 空变-多时窗融合高精度剩余静校正方法 |
CN105652316A (zh) * | 2014-12-01 | 2016-06-08 | 北京石大创新石油科技有限公司 | 一种基于裂缝模型的智能优化地震多属性融合方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5022240B2 (ja) * | 2008-01-08 | 2012-09-12 | インターナショナル・ビジネス・マシーンズ・コーポレーション | テープ媒体へのデータの上書き方法、及びその方法を具備した磁気テープ装置 |
-
2017
- 2017-12-28 CN CN201711457720.0A patent/CN109975884B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4964089A (en) * | 1989-12-22 | 1990-10-16 | Conoco Inc. | Method for derivation of internal velocities from post-migration parts |
EP0477806A2 (en) * | 1990-09-25 | 1992-04-01 | Honeywell Inc. | Method and apparatus for producing an accurate time base |
JPH09245074A (ja) * | 1996-03-11 | 1997-09-19 | Nec Corp | 回路接続情報の比較装置 |
US7489928B2 (en) * | 2002-12-31 | 2009-02-10 | Smith Brian K | Adaptive RF link failure handler |
CN1683939A (zh) * | 2004-04-13 | 2005-10-19 | 西门子公司 | 在磁共振断层造影中融合加权成像的运动校正的多点方法 |
CN103499836A (zh) * | 2013-07-18 | 2014-01-08 | 中国石油集团西北地质研究所 | 空变-多时窗融合高精度剩余静校正方法 |
CN105652316A (zh) * | 2014-12-01 | 2016-06-08 | 北京石大创新石油科技有限公司 | 一种基于裂缝模型的智能优化地震多属性融合方法 |
Non-Patent Citations (1)
Title |
---|
多源遥感影像融合技术在地质灾害调查中的应用;薛东剑等;《遥感技术与应用》;20111031;第26卷(第5期);第664-669页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109975884A (zh) | 2019-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Puetz et al. | Statistical analyses of global U-Pb database 2017 | |
Davis et al. | The European modern pollen database (EMPD) project | |
Satkoski et al. | Likeness among detrital zircon populations—An approach to the comparison of age frequency data in time and space | |
Charro et al. | Chemometric interpretation of vertical profiles of radionuclides in soils near a Spanish coal-fired power plant | |
CN111045072B (zh) | 一种适用于CeBr3探测器的γ能谱迭代解谱方法 | |
CN105158796B (zh) | 确定toc含量的方法和装置 | |
Fan et al. | Quantitative stratigraphy of the Wufeng and Lungmachi black shales and graptolite evolution during and after the Late Ordovician mass extinction | |
CA2456459C (en) | Method for enhancing depth and spatial resolution of one and two dimensional residual surfaces derived from scalar potential data | |
Tucaković et al. | 137Cs in mushrooms from Croatia sampled 15–30 years after Chernobyl | |
Link et al. | An automatized XKS-splitting procedure for large data sets: Extension package for SplitRacer and application to the USArray | |
CN103678778A (zh) | 一种放射性物化探信息综合方法 | |
Fagel et al. | Mineralogical signatures of Lake Baikal sediments: Sources of sediment supplies through Late Quaternary | |
CN109975884B (zh) | 一种放射性地球物理测量数据融合方法 | |
Biswas et al. | Curvelet transform to study scale-dependent anisotropic soil spatial variation | |
CN102692651A (zh) | 速度空变的初至波剩余静校正方法 | |
Ma et al. | A new global Rayleigh‐and Love‐wave group velocity dataset for constraining lithosphere properties | |
Zhong et al. | Radiocarbon and OSL dating on cores from the Chaoshan delta in the coastal South China Sea | |
Schmidt et al. | A chronological and palaeoenvironmental re‐evaluation of two loess‐palaeosol records in the northern Harz foreland, Germany, based on innovative modelling tools | |
CN201747364U (zh) | 一种声波时差曲线重构设备 | |
Green | Magnetic profile analysis | |
Conaway et al. | Quantitative uranium determinations from gamma-ray logs by application of digital time series analysis | |
SARMADIAN et al. | Mapping of spatial variability of soil organic carbon based on radial basis functions method | |
Li et al. | Progress in recent airborne gamma ray spectrometry measurement technology | |
CN111830562B (zh) | 一种油气储层渗透率预测方法及装置 | |
CN112926260A (zh) | 一种岩土材料物理力学参数准确估值新方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |