CN116611026A - 一种航空 γ 能谱数据融合处理方法和系统 - Google Patents
一种航空 γ 能谱数据融合处理方法和系统 Download PDFInfo
- Publication number
- CN116611026A CN116611026A CN202310598009.6A CN202310598009A CN116611026A CN 116611026 A CN116611026 A CN 116611026A CN 202310598009 A CN202310598009 A CN 202310598009A CN 116611026 A CN116611026 A CN 116611026A
- Authority
- CN
- China
- Prior art keywords
- window
- detector
- nai
- thorium
- uranium
- 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.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 371
- 238000000034 method Methods 0.000 title claims abstract description 97
- 238000007499 fusion processing Methods 0.000 title claims abstract description 27
- 238000005259 measurement Methods 0.000 claims abstract description 73
- 238000012937 correction Methods 0.000 claims abstract description 22
- 229910052776 Thorium Inorganic materials 0.000 claims description 280
- 229910052770 Uranium Inorganic materials 0.000 claims description 280
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 claims description 267
- ZSLUVFAKFWKJRC-IGMARMGPSA-N 232Th Chemical compound [232Th] ZSLUVFAKFWKJRC-IGMARMGPSA-N 0.000 claims description 265
- 229910052700 potassium Inorganic materials 0.000 claims description 191
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 claims description 189
- 239000011591 potassium Substances 0.000 claims description 189
- 238000004364 calculation method Methods 0.000 claims description 92
- 229910052705 radium Inorganic materials 0.000 claims description 50
- HCWPIIXVSYCSAN-UHFFFAOYSA-N radium atom Chemical compound [Ra] HCWPIIXVSYCSAN-UHFFFAOYSA-N 0.000 claims description 50
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 claims description 46
- 238000012545 processing Methods 0.000 claims description 39
- 238000001914 filtration Methods 0.000 claims description 31
- 230000009467 reduction Effects 0.000 claims description 20
- 230000035945 sensitivity Effects 0.000 claims description 20
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 18
- 230000004927 fusion Effects 0.000 claims description 17
- 238000012886 linear function Methods 0.000 claims description 16
- 239000013078 crystal Substances 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 7
- ZSLUVFAKFWKJRC-UHFFFAOYSA-N thorium Chemical compound [Th] ZSLUVFAKFWKJRC-UHFFFAOYSA-N 0.000 claims description 7
- 230000001419 dependent effect Effects 0.000 claims description 6
- 238000003672 processing method Methods 0.000 claims description 6
- DLFWIFNRAUYTHF-UHFFFAOYSA-N [Ra].[U] Chemical compound [Ra].[U] DLFWIFNRAUYTHF-UHFFFAOYSA-N 0.000 claims description 4
- 230000001360 synchronised effect Effects 0.000 claims description 4
- 238000003705 background correction Methods 0.000 claims description 3
- 230000010363 phase shift Effects 0.000 claims description 3
- 238000012935 Averaging Methods 0.000 claims description 2
- 229910004755 Cerium(III) bromide Inorganic materials 0.000 claims 2
- MOOUSOJAOQPDEH-UHFFFAOYSA-K cerium(iii) bromide Chemical compound [Br-].[Br-].[Br-].[Ce+3] MOOUSOJAOQPDEH-UHFFFAOYSA-K 0.000 claims 2
- 238000010606 normalization Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 10
- 230000009977 dual effect Effects 0.000 abstract description 4
- 230000001629 suppression Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 41
- 238000003860 storage Methods 0.000 description 14
- 230000008569 process Effects 0.000 description 13
- 230000005251 gamma ray Effects 0.000 description 11
- 238000004590 computer program Methods 0.000 description 10
- 230000002285 radioactive effect Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 238000009826 distribution Methods 0.000 description 6
- 230000003595 spectral effect Effects 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000001730 gamma-ray spectroscopy Methods 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 235000010755 mineral Nutrition 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000000084 gamma-ray spectrum Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000010006 flight Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- ZOKXTWBITQBERF-UHFFFAOYSA-N Molybdenum Chemical compound [Mo] ZOKXTWBITQBERF-UHFFFAOYSA-N 0.000 description 1
- WCUXLLCKKVVCTQ-UHFFFAOYSA-M Potassium chloride Chemical compound [Cl-].[K+] WCUXLLCKKVVCTQ-UHFFFAOYSA-M 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 229910052768 actinide Inorganic materials 0.000 description 1
- 150000001255 actinides Chemical class 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052750 molybdenum Inorganic materials 0.000 description 1
- 239000011733 molybdenum Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003758 nuclear fuel Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 239000002367 phosphate rock Substances 0.000 description 1
- OJMIONKXNSYLSR-UHFFFAOYSA-N phosphorous acid Chemical compound OP(O)O OJMIONKXNSYLSR-UHFFFAOYSA-N 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 229920001983 poloxamer Polymers 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 239000001103 potassium chloride Substances 0.000 description 1
- 235000011164 potassium chloride Nutrition 0.000 description 1
- 239000012857 radioactive material Substances 0.000 description 1
- 229910052761 rare earth metal Inorganic materials 0.000 description 1
- 150000002910 rare earth metals Chemical class 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/251—Fusion techniques of input or preprocessed data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
- G01T1/362—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry with scintillation detectors
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种航空γ能谱数据融合处理方法和系统，包括：γ能谱各能窗计数由全能峰计数和康普顿散射本底计数组成；以高能量分辨率的CeBr3探测器各能窗范围为基准参考，将高探测效率的大体积NaI(Tl)探测器获取的各能窗计数率数据修正为与CeBr3探测器各能窗具有相同的窄能量范围康普顿散射本底，并将两类探测器各能窗计数率数据进行合并，再采用统一的改正系数对合并后能窗数据进行处理，获取最终的天然核素含量值。本发明的优点是：发挥两种探测器高能量分辨率的双重优势，达到高效率全能峰探测及抑制非全能峰贡献的双重目标，提高整套系统航空γ能谱天然核素测量的精度。
Description
技术领域
本发明涉及电力系统调度技术领域，特别涉及一种基于NaI(Tl)探测器和CeBr3探测器的航空γ能谱数据融合处理方法和系统。
背景技术
自然界中存在着放射性物质衰变，为三个天然放射性系列：钍系、铀系和锕系。在衰变过程中，放射出α、β、γ粒子，形成α射线、β射线、γ射线。α、β射线自身能量很低，穿透性弱，不易测量，γ射线具有一定的能量，可以穿透一定厚度的物质，如：土壤、岩层等进入大气层中，因此可以使用航空方法来测量γ射线的强度。
航空γ能谱测量将航空γ能谱仪安装在飞行器上，在测量地区上空按照预先设计的测线和高度测量岩石或地层以及其他介质中放射性元素衰变放出的γ射线，获得多道γ能谱数据，通过数据处理获得所测量的放射性核素含量或辐射强度，是一种快速、经济、有效的核地球物理-地球化学测量方法。
航空γ能谱测量通过测量地面介质的特征γ射线谱来求得地面岩石或土壤中的铀、钍、钾的含量。可用于寻找放射性矿产、石油和天然气资源，与放射性核素伴生或共生的金属矿(铝、钼、金、稀土等)及钾盐、磷灰岩等其他矿产资源，后较多用于进行岩性地质填图和解决水文地质问题。
由于核试验及核工业、核电站、核燃料的后处理以及人工放射性核素的应用等，环境中也存在人工放射性核素，也是放射线污染物主要来源之一。
航空γ能谱测量是将航空γ能谱仪系统安装在飞机或其他飞行器中，通过按设定比例尺的测线进行低空飞行开展γ能谱测量。航空γ能谱测量为测量周期为1秒动态连续测量，并通过GPS定位系统确定测点的坐标信息，最终获得与各测点坐标对应的测量周期为1秒的多道γ能谱数据。
航空γ能谱测量为保证动态测量过程中各测点足够的粒子计数，降低统计涨落提高测量精度，需采用较高探测效率的探测器；航空γ射线测量过程中各种相互作用以及各种因素造成γ能谱的复杂化，γ谱数据分析感兴趣的只有γ射线的全能峰，只有全部损失其能量的γ光子的贡献才是最重要的，要获得最佳的谱数据，就要使全能峰的份额越大越好，其他成分的贡献越小越好。因此提高探测器的全能峰探测效率和降低其它干扰因素的贡献是主要两种有效手段。全能峰探测效率通过增大探测器的体积来实现；抑制非全能峰的贡献主要通过抑制康普顿散射分布来实现。
单个NaI(Tl)晶体可生成较大尺寸，通过较少的阵列数量即可实现大体积NaI(Tl)晶体探测器。现行航空γ能谱测量主要采用大体积阵列式NaI(Tl)晶体探测器，提高测量过程中各能窗的计数率，从而提高测量过程全能峰的探测效率，降低统计涨落对测量结果的影响，实现高信噪比测量。
随着放射性核素航空测量快速识别需求增加，具有较优的能量分辨率的CeBr3探测器也在航空γ能谱测量中得到应用，而单个CeBr3闪烁晶体难实现大体积生成，并且CeBr3晶体价格相对昂贵，阵列式探测器难以实现大体积。但CeBr3晶体探测器具有较优的能量分辨率，各单能γ射线的全能峰半高宽(fwhm)较窄，从而可实现能窗范围较窄的天然核素特征峰窗口；多道γ能谱特征峰窗口中计数由全能峰计数和康普顿散射本底计数组成，康普顿散射在多道γ能谱产生的本底计数是全能峰的主要干扰因素，而康普顿散射本底在多道谱中呈现连续分布，因此具有较优的能量分辨率的CeBr3晶体探测器在较窄天然核素特征峰窗口中全能峰计数比例较高，康普顿散射本底计数对全能峰干扰较小，达到降低康普顿散射分布影响的效果，即CeBr3晶体探测器各能窗中全能峰窗口中信噪比较优。
大体积NaI(Tl)闪烁体探测器较CeBr3闪烁体探测器在全能峰探测效率方面具有优势，而CeBr3闪烁体探测器较大体积NaI(Tl)闪烁体探测器在抑制康普顿散射本底方面具有优势，因此为最大程度发挥NaI(Tl)闪烁体探测器和CeBr3闪烁体探测器在航空γ能谱测量中各自的优势，航空γ能谱测量系统可采用NaI和CeBr3组合探测器，从提高全能峰探测效率和抑制康普顿散射本底两方面来实现测量精度的提高；而NaI和CeBr3组合探测器在天然放射性核素测量过程中，两种探测器是分别获得每个测点的多道能谱数据，在后期处理过程中需将两种探测器的数据进行融合处理才能发挥组合探测器在航空γ能谱测量中的优势，获取更高精度的地面天然放射性核素K、U和Th的含量。而现有技术中并没有相应的NaI(Tl)和CeBr3组合探测器航空γ能谱测量数据融合处理方法。
缩略语和关键术语定义
γ能谱测量
γ能谱测量是利用γ能谱仪测量岩石或地层以及其他介质中放射性元素衰变放出的具有特定能量的γ射线，从而来测定辐射强度或放射性元素含量的一种方法。γ能谱测量方法可用来进行岩矿石样品分析，放射性矿产勘查，岩性分类和地质填图，水资源勘查，工程地质中裂隙和断层确定，与放射性元素共生的非放射性矿产或油气勘查，以及放射性辐射环境评价等。
航空γ能谱测量
航空γ能谱测量(airborne gamma-ray spectrometry survey)将航空γ能谱仪系统安装在飞机或其他飞行器中，是按设定比例尺的测线进行低空飞行，进行γ射线能谱测量。具有效率高、速度快、不受地形等因素限制、适合大面积普查的特点。
全能峰
单能γ射线光子进入闪烁体探测器后损失全部能量产生全能脉冲经幅度分析后在γ能谱仪获得的多道γ能谱谱线中形成的峰形，闪烁体探测器获取的γ射线全能峰反映了统计性质，符合正态分布，可用高斯函数来描述。
特征峰
不同放射性核素产生的特定能量的γ射线，不同能量的特征γ射线在γ能谱仪测量的多道谱线上形成的峰。由于γ射线进入闪烁体探测器产生光电效应、康普顿效应和电子对效应，在多道γ能谱谱线中特征峰包括全能峰和康普顿散射本底两部分，全能峰符合高斯型，康普顿散射本底可采用线性函数来描述，因此可用线性函数和高斯函数来描述特征峰。
能窗
闪烁体探测器测得的多道γ能谱数据解析一般采用能窗法，即在多道γ能谱中划分出与目标探测核素的特征γ射线能量相对应的能窗范围，利用各能窗范围内积分计数与标准源的刻度关系，从而求得核素的含量，能窗选择应用能表征目标核素且能区别于其他核素特征峰。在天然放射性核素探测中，通常选取铀系的214Bi的1.76MeVγ射线特征峰能窗、钍系的208Tl的2.62MeVγ射线特征峰能窗和40K的1.46MeV的γ射线特征峰能窗分别来测量铀(238U)、钍(232Th)和钾(40K)的含量。
发明内容
本发明针对现有技术的缺陷，提供了一种航空γ能谱数据融合处理方法和系统。
为了实现以上发明目的，本发明采取的技术方案如下：
一种航空γ能谱数据融合处理方法，包括以下步骤：
一、数据融合
系统获取的钾窗、铀窗和钍窗的计数率分别为同步获取的NaI(Tl)探测器获取的相应能窗计数率与CeBr3探测器获取的相应能窗计数率之和；
NaI(Tl)探测器各特征峰能窗计数率由能窗中全能峰计数率和康普顿散射本底计数率构成。为将NaI(Tl)探测器与CeBr3探测器的能窗数据融合后统一处理，NaI(Tl)探测器的康普顿散射本底各能窗范围采用CeBr3探测器所确定能窗范围；
各航空模型包括：本底模型、钾模型、铀模型、钍模型；
二、确定天然核素特征峰窗口范围；
2.1特征峰曲线拟合
在用于检定航空γ能谱仪的本底模型、钾模型、铀模型、钍模型上分别进行测量，同步获得NaI(Tl)和CeBr3探测器两种512道γ能谱数据，进行平均处理获得归一到单个测量周期的512道能谱数据，即512道γ能谱计数率。分别对NaI(Tl)和CeBr3两种探测器测得的钾模型上的512道能谱计数率扣除本底模型上的512道能谱计数率获得纯40K的512道γ能谱计数率；分别对NaI(Tl)和CeBr3两种探测器测得的铀模型上的512道能谱计数率扣除本底模型上的512道能谱计数率获得纯铀镭系的512道γ能谱计数率；分别对NaI(Tl)和CeBr3两种探测器测得的钍模型上的512道能谱计数率扣除本底模型上的512道能谱计数率获得纯钍系的512道γ能谱计数率。
用线性函数加高斯函数来描述512道γ能谱中40K的特征峰、214Bi的特征峰和208Tl的特征峰。
以道址作为自变量x，以道址x上计数率作为应变量y，分别对NaI(Tl)探测器和CeBr3探测器获取的纯40K 512道γ能谱计数率中的特征峰、纯铀镭系512道γ能谱计数率中的特征峰、纯钍系512道γ能谱计数率中的特征峰采用线性函数加高斯函数进行最小二乘拟合；
2.2确定特征峰能窗范围
确定CeBr3探测器和NaI探测器的航空γ能谱数据各能窗范围，分别为:钾窗的起始道址，钾窗的结束道址；铀窗的起始道址，铀窗的结束道址；钍窗的起始道址，钍窗的结束道址。
三、确定康普顿散射系数；
航空γ能谱数据处理的能窗分别为钾窗、铀窗和钍窗，康普顿散射剥离系数是纯40K、纯铀镭、纯钍源在两个不同能窗中的计数率之比，分别采用α、β、γ、a、b和g来表示，其中α为纯钍源在铀窗和钍窗上的计数率之比；β为纯钍源在钾窗和钍窗上的计数率之比；γ为纯铀镭源在钾窗和铀窗上的计数率之比；a为纯铀镭源在钍窗和铀窗中的计数率之比，b为纯40K源在钍窗和钾窗上的计数率之比；g为纯40K源在铀窗和钾窗上的计数率之比。康普顿散射系数b和g均为0，计算剩余康普顿散射系数，公式如下：
其中nK,U、nU,U、nTh,U分别为航空模型坪纯铀镭源在钾窗、铀窗和钍窗上产生的计数率，单位为计数每秒(s-1)；nK,Th、nU,Th、nTh,Th分别为航空模型坪纯钍源在钾窗、铀窗和钍窗上产生的计数率，单位为计数每秒(s-1)。
各能窗的计数率为NaI(Tl)探测器和CeBr3探测器能窗计数率之和，为将两类探测器的能窗数据融合，NaI(Tl)探测器各能窗计数率视为能窗中全能峰计数率和康普顿散射本底计数率之和，其中纯铀镭源在NaI(Tl)探测器的钾窗和钍窗上无全能峰计数，只有康普顿散射本底计数；纯钍源在NaI(Tl)探测器的钾窗和铀窗上无全能峰计数，只有康普顿散射本底计数。求解康普顿散射系数所需各能窗计数率计算公式如下：
公式中为航空模型坪纯铀镭源在NaI(Tl)探测器的铀窗中产生的全能峰计数率，单位为计数每秒(s-1)；/>为航空模型坪纯钍源在NaI(Tl)探测器的钍窗中产生的全能峰计数率，单位为计数每秒(s-1)；/>分别为航空模型坪纯铀镭源在NaI(Tl)探测器的钾窗、铀窗和钍窗中产生的康普顿散射本底计数率，单位为计数每秒(s-1)；分别为航空模型坪纯钍源在NaI(Tl)探测器的钾窗、铀窗和钍窗中产生的康普顿散射本底计数率，单位为计数每秒(s-1)；/>分别为航空模型坪纯铀镭源在CeBr3探测器的钾窗、铀窗和钍窗中产生的计数率，单位为计数每秒(s-1)；分别为航空模型坪纯钍源在CeBr3探测器的钾窗、铀窗和钍窗中产生的计数率，单位为计数每秒(s-1)。
四、飞机本底及宇宙射线系数计算；
飞机本底和宇宙射线系数的校准采用五个高高度飞行，每个高度飞行10分钟，获取的能谱数据来计算。
对每个高高度上NaI(Tl)探测器获得512道能谱数据进行累加平均获得每个高高度上归一到单个测量周期的512道γ能谱计数率(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，单位为计数每秒(s-1)。
对每个高高度上CeBr3探测器获得512道能谱数据进行累加平均获得每个高度上归一到单个测量周期的512道γ能谱计数率单位为计数每秒(s-1)。
4.1分别对每个高高度上NaI(Tl)探测器512道能谱计数率降噪处理；
4.2分别计算每个高高度上NaI(Tl)探测器512道能谱康普顿散射本底；
4.3分别计算每个高高度上NaI(Tl)探测器钾窗、铀窗、钍窗计数率；
4.4分别计算每个高高度上CeBr3探测器钾窗、铀窗、钍窗计数率数据；
4.5分别计算每个高高度上整套测量系统钾窗、铀窗和钍窗计数率；
计算整套测量系统在高高度h上钾窗计数率单位为计数每秒(s-1)，h表示高度，h＝3000,3300,3600,3900,4200；
计算整套测量系统在高高度h上铀窗计数率单位为计数每秒(s-1)，h表示高度，h＝3000,3300,3600,3900,4200；
计算整套测量系统在高高度h上钍窗计数率单位为计数每秒(s-1)，h表示高度，h＝3000,3300,3600,3900,4200；
4.6获取飞机本底及宇宙射线系数；
所有高高度进行飞行测量，宇宙射线道中的计数率与各能窗的计数率呈线性关系；
对每个高高度层的宇宙射线道的平均计数率和能窗计数率进行最小二乘法线性拟合，获得线性关系系数，即钾窗飞机本底计数率，单位为计数每秒(s-1)；钾窗宇宙射线影响系数，无纲量；铀窗飞机本底计数率，单位为计数每秒(s-1)；铀窗宇宙射线影响系数，无纲量；钍窗飞机本底计数率，单位为计数每秒(s-1)；钍窗宇宙射线影响系数，无纲量。
五、高度衰减系数及空中灵敏度计算
5.1NaI(Tl)探测器各能窗计数率计算
对动态带各个高度上NaI(Tl)探测器获得的水域部分和陆域部分512道谱数据分别进行累加平均获得每个高度水域上512道γ能谱平均计数率和陆域上512道γ能谱数据计数率/>表示NaI(Tl)探测器在动态带水域部分上空飞行高度为h获得的第i道γ能谱计数率，单位为计数每秒(s-1)，/>表示NaI(Tl)探测器在动态带陆域部分上空飞行高度h获得的第i道γ能谱计数率，单位为计数每秒(s-1)，其中h＝30、60、90、120、150、210、240、300，i＝1、2、3、...、511、512。
和/>作差为NaI(Tl)探测器在高度h上获得动态带纯陆地放射性核素产生的512道γ能谱中第i道的平均计数率/>单位为计数每秒(s-1)：
将每个高度层上的512道γ能谱计数率采用4.1至4.3相同的处理方法，获得NaI(Tl)探测器在每个高度上动态带纯陆地放射性核素产生的γ能谱各能窗计数率：
为NaI(Tl)探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钾窗计数率，单位为计数每秒(s-1)；
为NaI(Tl)探测器在高度h上动态带纯陆地放射性核素产生的γ能谱铀窗计数率，单位为计数每秒(s-1)；
为NaI(Tl)探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钍窗计数率，单位为计数每秒(s-1)。
5.2CeBr3探测器各能窗计数率计算
对动态带各个高度上CeBr3探测器获得的水域部分和陆域部分512道谱数据进行累加平均分别获得每个高度水域上512道γ能谱平均计数率和陆域上512道γ能谱数据计数率/>表示CeBr3探测器在动态带水域部分上空飞行高度为h获得的第i道γ能谱计数率，/>表示CeBr3探测器在动态带陆域部分上空飞行高度h获得的第i道γ能谱计数率，单位为计数每秒(s-1)，其中h＝30、60、90、120、150、210、240、300，i＝1、2、3、...、511、512。
和/>作差为CeBr3探测器飞行高度h获得动态带纯陆地放射性核素产生的512道γ能谱中第i道的平均计数率/>单位为计数每秒(s-1)：
将每个高度层上的512道γ能谱计数率采用4.4相同的处理方法，获得CeBr3探测器在每个高度上动态带纯陆地放射性核素产生的γ能谱各能窗计数率：
为CeBr3探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钾窗计数率，单位为计数每秒(s-1)；
为CeBr3探测器在高度h上动态带纯陆地放射性核素产生的γ能谱铀窗计数率，单位为计数每秒(s-1)；
为CeBr3探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钍窗计数率，单位为计数每秒(s-1)。
5.3整套测量系统各能窗计数率计算
计算整套测量系统在高度h上动态带纯陆地放射性核素产生的γ能谱钾窗计数率单位为计数每秒(s-1)；
计算整套测量系统在高度h上动态带纯陆地放射性核素产生的γ能谱铀窗计数率单位为计数每秒(s-1)；
计算整套测量系统在高度h上动态带纯陆地放射性核素产生的γ能谱钍窗计数率单位为计数每秒(s-1)；
5.4高度衰减系数和空中灵敏度计算
在动态带各高度上纯陆地放射性核素产生的各能窗计数率进行康普顿散射改正；
利用30m、60m、90m、120m、150m、210m、240m、、300八个高度上获取的钾窗计数率、铀窗计数率和钍窗计数率，单位为计数每秒(s-1)，按e指数公式拟合就求得各能量窗的高度衰减系数μK、μU和μTh，单位为每米(m-1)。
5.5空中灵敏度计算
计算120m高度空中系统各能窗灵敏度。
六、测点数据的处理
航空γ能谱测量过程中获得每个测点的NaI(Tl)探测器512道γ能谱计数率数据和CeBr3探测器512道γ能谱计数率数据/>其中i表示道址，i＝1,2,3,...,511,512。
6.1NaI(Tl)探测器各能窗数据计算
每个测点的NaI(Tl)探测器512道γ能谱计数率采用与4.1至4.3中每个高高度层上NaI(Tl)探测器的512道γ能谱计数率/>相同的处理方法，获得NaI(Tl)探测器在每个测点上各能窗计数率：
为NaI(Tl)探测器在每个测点上获得的钾窗计数率，单位为计数每秒(s-1)；
为NaI(Tl)探测器在每个测点上获得的铀窗计数率，单位为计数每秒(s-1)；
为NaI(Tl)探测器在每个测点上获得的钍窗计数率，单位为计数每秒(s-1)。
6.2CeBr3探测器各能窗计数率计算
将每个测点的CeBr3探测器512道γ能谱计数率采用与4.4中每个高高度层上CeBr3探测器的512道γ能谱计数率/>相同的处理方法，获得CeBr3探测器在每个测点上各能窗计数率：
为CeBr3探测器在测点上获得的钾窗计数率，单位为计数每秒(s-1)；
为CeBr3探测器在测点上获得的铀窗计数率，单位为计数每秒(s-1)；
为CeBr3探测器在测点上获得的钍窗计数率，单位为计数每秒(s-1)。
6.3整套测量系统各能窗计数率计算
计算整套测量系统在各测点上获得的钾窗计数率nK，单位为计数每秒(s-1)；
计算整套测量系统在各测点上获得的铀窗计数率nU，单位为计数每秒(s-1)；
计算整套测量系统在各测点上获得的钍窗计数率nTh，单位为计数每秒(s-1)；
6.4飞机本底及宇宙射线本底改正
测点上各能窗计数率数据飞机本底及宇宙射线本底改正计算；
6.5康普顿散射改正
测点上各能窗计数率数据康普顿散射改正计算；
6.6高度衰减改正
计算将测点上各能窗计数率高度改正到120m高度计数率
6.7放射性核素含量计算
计算每个测点地面上钾、铀、钍含量，公式如下：
QK QU QTh分别为测点处钾、铀和钍的含量，单位分别为％、ug/g和ug/g；和分别为钾窗、铀窗和钍窗在空中120m高度的灵敏度，单位分别为1/(s·％)、1/(s·ug/g)和1/(s·ug/g)；/>分别为经飞机本底、宇宙射线本底和康普顿散射改正及高度改正到空中120米高度的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)。
进一步地，步骤2.1中用线性函数加高斯函数来描述512道γ能谱中40K的1.46MeV特征峰、214Bi的1.76MeV特征峰和208Tl的2.62MeV特征峰，高斯函数公式如下：
以道址作为自变量x，以各道上的计数率作为应变量y，分别对NaI(Tl)探测器和CeBr3探测器获取的纯40K 512道γ能谱计数率中的1.46MeV特征峰、纯铀镭系512道γ能谱计数率中的1.76MeV特征峰、纯钍系512道γ能谱计数率中的2.62MeV特征峰采用公式16表示的函数进行最小二乘拟合。获得如下公式17-22。
其中和/>分别为采用公式16函数拟合NaI(Tl)探测器纯40K 512道γ能谱计数率中40K的1.46MeV特征峰的各项系数；
其中和/>分别为采用公式16函数拟合NaI(Tl)探测器纯铀镭系512道γ能谱计数率中214Bi的1.76MeV特征峰的各项系数；
其中和/>分别为采用公式16函数拟合NaI(Tl)探测器纯钍系512道γ能谱计数率中208Tl的2.62MeV特征峰获得的各项系数；
其中和/>分别为采用公式16函数拟合CeBr3探测器纯40K 512道γ能谱计数率中40K的1.46MeV特征峰获得的各项系数；
其中和/>分别为采用公式16函数拟合CeBr3探测器获取纯铀镭系512道γ能谱计数率中214Bi的1.76MeV特征峰获得的各项系数；
其中和/>分别为采用公式16函数拟合CeBr3探测器获取纯钍系512道γ能谱计数率中208Tl的2.62MeV特征峰获得的各项系数。
进一步地，步骤2.2中特征峰能窗范围的确定
CeBr3探测器的航空γ能谱数据各能窗范围为:取整作为钾窗的起始道址/> 取整作为钾窗的结束道址/> 取整作为铀窗的起始道址/> 取整作为铀窗的结束道址chC U eEd；
取整作为钍窗的起始道址/> 取整作为钍窗的结束道址/>
NaI(Tl)探测器的航空γ能谱数据各能窗范围为:取整作为钾窗的起始道址/> 取整作为钾窗的结束道址/> 取整作为铀窗的起始道址/> 取整作为铀窗的结束道址/> 取整作为钍窗的起始道址/> 取整作为钍窗的结束道址/>
进一步地，步骤三中，
3.1计算航空模型坪纯铀镭源在NaI(Tl)探测器的铀窗中产生的全能峰计数公式为：
式中和/>分别为公式18中采用线性函数加高斯函数公式拟合NaI(Tl)探测器纯铀镭系512道γ能谱计数率中214Bi的1.76MeV特征峰获得的系数，/>为铀窗的起始道址和/>为铀窗的结束道址。
3.2计算航空模型坪纯钍源在NaI(Tl)探测器的钍窗中全能峰的计数率公式为：
式中和/>为公式19中采用线性函数加高斯函数公式拟合NaI(Tl)探测器纯钍系512道γ能谱计数率中208Tl的2.62MeV特征峰获得获得的系数；/>NaI(Tl)探测器的钍窗的起始道址。/>为NaI(Tl)探测器的钍窗的结束道址。
3.3计算航空模型坪纯铀镭源在NaI(Tl)探测器的铀窗中产生的康普顿散射本底计数率公式为：
式中为公式18中采用线性函数加高斯函数公式拟合NaI(Tl)探测器纯铀镭系512道γ能谱计数率中214Bi的1.76MeV特征峰获得的系数；/>为CeBr3探测器的铀窗的起始道址，/>为CeBr3探测器的铀窗的结束道址。
3.4计算航空模型坪纯钍源在NaI(Tl)探测器的钍窗中康普顿散射本底计数率公式为：
式中为公式19中采用线性函数加高斯函数公式拟合NaI(Tl)探测器纯钍系512道γ能谱计数率中208Tl的2.62MeV特征峰获得的系数；/>为CeBr3探测器的钍窗的起始道址，/>为CeBr3探测器的钍窗的结束道址。
3.5计算在航空模型坪纯铀镭源上CeBr3探测器铀窗计数率公式为：
式中和/>分别为公式21中采用线性函数加高斯函数公式拟合CeBr3探测器获取纯铀镭系512道γ能谱计数率中214Bi的1.76MeV特征峰获得的系数，/>为CeBr3探测器的铀窗的起始道址，/>为CeBr3探测器的铀窗的结束道址。
3.6计算在航空模型坪纯钍源上CeBr3探测器钍窗计数率公式为：
式中和/>分别为公式22中采用线性函数加高斯函数公式拟合CeBr3探测器获取纯钍系512道γ能谱计数率中208Tl的2.62MeV特征峰获得的系数，/>为CeBr3探测器的钍窗的起始道址，/>为CeBr3探测器的钍窗的结束道址。/>
进一步地，4.1分别对每个高高度上NaI(Tl)探测器512道能谱计数率降噪处理，具体如下：
分别将每个高高度上NaI(Tl)探测器的512道γ能谱计数率视为采样频率为1000Hz的数字信号，采用通带频率为50Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得滤波后的512道γ能谱计数率/>(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，单位为计数每秒(s-1)。
为实现滤波无相位位移，保证特征峰峰位无位移，对滤波后获得的512道γ能谱计数率进行转置，获取512道γ能谱计数率/>(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，单位为计数每秒(s-1)；
将512道γ能谱计数率视为采样频率为1000Hz的数字信号，采用通带频率为50Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得滤波后的512道γ能谱计数率ri h(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，单位为计数每秒(s-1)。
512道γ能谱计数率ri h进行转置，获得各高高度上降噪滤波后的512道γ能谱计数率(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，单位为计数每秒(s-1)；
进一步地，4.2分别计算每个高高度上NaI(Tl)探测器512道能谱康普顿散射本底，具体如下：
将每个高高度层上NaI(Tl)探测器的512道γ能谱计数率采用4.1相似的处理流程，仅在FIR滤波器参数上有差别，选择通带频率为30Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得每个高高度层上的512道γ能谱计数率/>
(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，单位为计数每秒(s-1)；
将每个高高度层上NaI(Tl)探测器的512道γ能谱计数率采用4.1相似的处理流程，仅在FIR滤波器参数上有差别，选择通带频率为10Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得每个高度层上的512道γ能谱计数率/>(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，单位为计数每秒(s-1)；
每个高高度层取的50至210道的能谱计数率数据作为对应道址的康普顿散射本底计数，取/>的211至512道的能谱计数率数据作为对应道址的本底计数，构成每个高度层从50道至512道的γ能谱康普顿散射本底计数率数据/>(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，单位为计数每秒(s-1)：
以道址作为自变量x，以各道上康普顿散射本底计数率作为应变量y，分别对每个高高度层NaI(Tl)晶体探测器的50道至512道的γ能谱散射本底计数率采用幂函数y＝w·xm+t进行最小二乘拟合，获得各高高度层NaI(Tl)探测器50道至512道康普顿散射本底描述函数。/>
进一步地，4.3分别计算每个高高度上NaI(Tl)探测器钾窗、铀窗、钍窗计数率，具体如下：
4.3.1各高高度层上NaI(Tl)探测器钾窗计数率计算方法为公式：
式中为高高度为h的NaI(Tl)探测器钾窗计数率，单位为计数每秒(s-1)；为高高度为h的NaI(Tl)探测器钾窗中全能峰计数率，单位为计数每秒(s-1)；/>为高高度为h的NaI(Tl)探测器钾窗中康普顿散射本底计数率，单位为计数每秒(s-1)；Ni h为获取降噪滤波后高度h的NaI(Tl)探测器第i道γ能谱计数率，单位为计数每秒(s-1)；wh、mh、th为高度为h的NaI(Tl)探测器50道至512道的γ能谱散射本底计数率幂函y＝w·xm+t拟合所得参数；/>为NaI(Tl)探测器钾窗范围起止道址，/>为CeBr3探测器钾窗范围起止道址。
4.3.2各高度层上NaI(Tl)探测器铀窗计数率计算方法为公式：
式中为高高度为h的NaI(Tl)探测器铀窗计数率，单位为计数每秒(s-1)；为高高度为h的NaI(Tl)探测器铀窗中全能峰计数率，单位为计数每秒(s-1)；/>为高高度为h的NaI(Tl)探测器铀窗中康普顿散射本底计数率，单位为计数每秒(s-1)；/>为降噪滤波后高高度为h的NaI(Tl)探测器的第i道γ能谱计数率，单位为计数每秒(s-1)；wh、mh、th为高高度为h的NaI(Tl)探测器50道至512道的γ能谱散射本底计数率幂函y＝w·xm+t拟合所得参数；/>分别为NaI(Tl)探测器铀窗范围起止道址，/>分别为CeBr3探测器铀窗范围起止道址。
4.3.3各高度层上NaI(Tl)探测器钍窗计数率计算方法为公式：
式中为高高度为h的NaI(Tl)探测器钍窗计数率，单位为计数每秒(s-1)；为高高度为h的NaI(Tl)探测器钍窗中全能峰计数率，单位为计数每秒(s-1)；/>为高高度为h的NaI(Tl)探测器钍窗中康普顿散射本底计数率，单位为计数每秒(s-1)；Ni h为降噪滤波后高高度为h的NaI(Tl)探测器的第i道γ能谱计数率，单位为计数每秒(s-1)；wh、mh、th为高高度为h的NaI(Tl)探测器50道至512道的γ能谱散射本底计数率幂函y＝w·xm+t拟合所得参数；/>分别为NaI(Tl)探测器钍窗范围起止道址，/>分别为CeBr3探测器钍窗范围起止道址。
进一步地，所述系统获取的钾窗、铀窗和钍窗的计数率分别为同步获取的NaI(Tl)探测器获取的相应能窗计数率与CeBr3探测器获取的相应能窗计数率之和，计算公式如下：
nK、nU和nTh分别为数据融合后整套测量系统的航空γ能谱钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)；和/>分别为CeBr3探测器航空γ能谱的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)；/>和/>分别为NaI(Tl)探测器航空γ能谱的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)。
本发明还公开了一种航空γ能谱数据融合处理系统，该系统能够用于实施上述的一种航空γ能谱数据融合处理方法，具体的，包括：数据融合模块、天然核素特征峰窗口范围确定模块、康普顿散射系数确定模块、计算模块和测点数据处理模块；
数据融合模块：获取的钾窗、铀窗和钍窗的计数率；
获取NaI(Tl)探测器各特征峰能窗计数率。为将NaI(Tl)探测器与CeBr3探测器的能窗数据融合后统一处理，NaI(Tl)探测器的康普顿散射本底各能窗范围采用CeBr3探测器所确定能窗范围；
天然核素特征峰窗口范围确定模块：用于确定天然核素特征峰窗口范围；
康普顿散射系数确定模块：用于确定康普顿散射系数；
计算模块：用于计算飞机本底、宇宙射线系数、高度衰减系数和空中灵敏度；
测点数据处理模块；用于测点数据的处理。将两类探测器各能窗计数率数据进行合并，再采用统一的改正系数对合并后能窗数据进行处理，获取最终的天然核素含量值。
与现有技术相比，本发明的优点在于：
实现对同步获取的NaI(Tl)探测器航空γ能谱数据和CeBr3探测器航空γ能谱测量数据的融合处理，发挥大体积NaI(Tl)探测器高探测效率和CeBr3探测器高能量分辨率的双重优势，达到高效率全能峰探测及非全能峰贡献(康普顿散射本底)抑制的双重目标，有效提高整套系统航空γ能谱天然核素测量过程的信噪比，实现天然放射性核素高精度航空测量。
附图说明
图1是本发明实施例一种航空γ能谱数据融合处理方法流程图；
图2是本发明实施例NaI(Tl)探测器纯40K 512道γ能谱谱线及1.46MeV特征峰拟合曲线图；
图3是本发明实施例NaI(Tl)探测器纯铀镭系512道γ能谱谱线及1.76MeV特征峰拟合曲线图；
图4是本发明实施例NaI(Tl)探测器纯钍系512道γ能谱谱线及2.62MeV特征峰拟合曲线图；
图5是本发明实施例CeBr3探测器纯40K 512道γ能谱谱线及1.46MeV特征峰拟合曲线图；
图6是本发明实施例CeBr3探测器纯铀镭系512道γ能谱谱线及1.76MeV特征峰拟合曲线图；
图7是本发明实施例CeBr3探测器纯钍系512道γ能谱谱线及2.62MeV特征峰拟合曲线图；
图8是本发明实施例NaI(Tl)探测器512道γ能谱数据降噪滤波处理流程图；
图9是本发明实施例NaI(Tl)探测器512道γ能谱数据滤波处理流程图(通带频率为30Hz)；
图10是本发明实施例NaI(Tl)探测器512道γ能谱数据滤波处理流程图(通带频率为10Hz)；
图11是本发明实施例NaI(Tl)探测器512道γ能谱数据各处理谱图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白，以下根据附图并列举实施例，对本发明做进一步详细说明。
如图1所示，一种航空γ能谱数据融合处理方法流程如下：
一、数据融合
为充分利用NaI(Tl)探测器与CeBr3探测器各自在航空γ能谱测量中的优势，实现将NaI(Tl)探测器与CeBr3探测器同步获取的航空γ能谱数据最终转换为的天然放射性核素含量，需将NaI(Tl)探测器与CeBr3探测器同步测量获得的航空γ能谱数据融合处理，即将NaI(Tl)探测器天然核素各特征峰能窗计数率与CeBr3天然核素各特征峰能窗计数率合并后进行统一处理。
航空γ能谱整套测量系统获取的钾窗、铀窗和钍窗的计数率分别为同步获取的NaI(Tl)探测器获取的相应能窗计数率与CeBr3探测器获取的相应能窗计数率之和，如式1-式3：
nK、nU和nTh分别为数据融合后整套测量系统的航空γ能谱钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)；和/>分别为CeBr3探测器航空γ能谱的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)；/>和/>分别为NaI(Tl)探测器航空γ能谱的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)。
NaI(Tl)探测器各特征峰能窗计数率由能窗中全能峰计数率和康普顿散射本底计数率构成。为将NaI(Tl)探测器与CeBr3探测器的能窗数据融合后统一处理，NaI(Tl)探测器的康普顿散射本底各能窗范围采用CeBr3探测器所确定能窗范围，如式4-式6：
/>
和/>分别为NaI(Tl)探测器航空γ能谱数据的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)；/>和/>分别为NaI(Tl)探测器航空γ能谱数据的钾窗、铀窗和钍窗中全能峰的计数率，单位为计数每秒(s-1)；
和/>分别为CeBr3探测器钾窗、铀窗和钍窗所对应能窗范围内NaI(Tl)探测器γ能谱康普顿散射本底的计数率，单位为计数每秒(s-1)。
二、天然核素特征峰窗口范围的确定
2.1特征峰曲线拟合
在用于检定航空γ能谱仪的本底模型、钾模型、铀模型、钍模型上分别进行5分钟测量，同步获得NaI(Tl)和CeBr3探测器两种512道γ能谱数据，进行平均处理获得归一到单个测量周期(1秒)的512道能谱数据，即512道γ能谱计数率。分别对NaI(Tl)和CeBr3两种探测器测得的钾模型上的512道能谱计数率扣除本底模型上的512道能谱计数率获得纯40K的512道γ能谱计数率(公式7和公式10)；分别对NaI(Tl)和CeBr3两种探测器测得的铀模型上的512道能谱计数率扣除本底模型上的512道能谱计数率获得纯铀镭系的512道γ能谱计数率(公式8和公式11)；分别对NaI(Tl)和CeBr3两种探测器测得的钍模型上的512道能谱计数率扣除本底模型上的512道能谱计数率获得纯钍系的512道γ能谱计数率(公式9和公式12)。
公式7-9中和/>分别为NaI(Tl)探测器获取的纯40K、纯铀镭系、纯钍系的512道γ能谱中第i道计数率(i＝1,2,3,…,511,512)，单位为计数每秒(s-1)；和/>分别为NaI(Tl)探测器在钾模型、铀模型和钍模型上测得的512道能谱数据中第i道计数率(i＝1,2,3,…,511,512)，单位为计数每秒(s-1)；/>为NaI(Tl)探测器在本底模型上测得的512道能谱数据中第i道计数率(i＝1,2,3,…,511,512)，单位为计数每秒(s-1)。
公式10-12中和/>分别为CeBr3晶体探测器获取的纯40K、纯铀镭系、纯钍系的512道γ能谱中第i道计数率(i＝1,2,3,…,511,512)，单位为计数每秒(s-1)；和/>分别为CeBr3探测器在钾模型、铀模型和钍模型上测得的512道能谱数据中第i道计数率(i＝1,2,3,…,511,512)，单位为计数每秒(s-1)；/>为CeBr3探测器在本底模型上测得的512道能谱数据中第i道计数率(i＝1,2,3,…,511,512)，单位为计数每秒(s-1)。
多道γ能谱中特征峰计数包括全能峰计数和康普顿散射本底计数。单能射线γ光子进入闪烁体后损失全部能量产生的全能脉冲呈高斯型幅度分布，闪烁晶体探测器获取的γ射线全能峰反映了统计性质，符合正态分布，采用高斯函数来描述；康普顿散射本底包括高能γ射线康普顿散射贡献及其本身的康普顿散射贡献，采用线性函数来描述；因此用线性函数加高斯函数(公式13)来描述512道γ能谱中40K的1.46MeV特征峰、214Bi的1.76MeV特征峰和208Tl的2.62MeV特征峰。
以道址作为自变量x，以各道上的计数率作为应变量y，分别对NaI(Tl)探测器和CeBr3探测器获取的纯40K 512道γ能谱计数率中的1.46MeV特征峰、纯铀镭系512道γ能谱计数率中的1.76MeV特征峰、纯钍系512道γ能谱计数率中的2.62MeV特征峰采用公式13表示的函数进行最小二乘拟合，获得如下公式14-19，拟合结果如图2至7。
其中和/>分别为采用公式13函数拟合NaI(Tl)探测器纯40K 512道γ能谱计数率中40K的1.46MeV特征峰的各项系数；
其中和/>分别为采用公式13函数拟合NaI(Tl)探测器纯铀镭系512道γ能谱计数率中214Bi的1.76MeV特征峰的各项系数；
其中和/>分别为采用公式13函数拟合NaI(Tl)探测器纯钍系512道γ能谱计数率中208Tl的2.62MeV特征峰获得的各项系数；
其中和/>分别为采用公式13函数拟合CeBr3探测器纯40K 512道γ能谱计数率中40K的1.46MeV特征峰获得的各项系数；
其中和/>分别为采用公式13函数拟合CeBr3探测器获取纯铀镭系512道γ能谱计数率中214Bi的1.76MeV特征峰获得的各项系数；
其中和/>分别为采用公式13函数拟合CeBr3探测器获取纯钍系512道γ能谱计数率中208Tl的2.62MeV特征峰获得的各项系数。/>
2.2特征峰能窗范围的确定
CeBr3探测器的航空γ能谱数据各能窗范围为:取整作为钾窗的起始道址/>(公式20)，/>取整作为钾窗的结束道址/>(公式21)；取整作为铀窗的起始道址/>(公式22)，/>取整作为铀窗的结束道址/>(公式23)；/>取整作为钍窗的起始道址/>(公式24)，/>取整作为钍窗的结束道址/>(公式25)。
NaI探测器的航空γ能谱数据各能窗范围为:取整作为钾窗的起始道址/>(公式26)，/>取整作为钾窗的结束道址/>(公式27)；/>取整作为铀窗的起始道址/>(公式28)，/>取整作为铀窗的结束道址/>(公式29)；/>取整作为钍窗的起始道址/>(公式30)，/>取整作为钍窗的结束道址/>(公式31)。
三、康普顿散射系数的确定
航空γ能谱数据处理的能窗分别为钾窗、铀窗和钍窗，康普顿散射剥离系数是纯40K、纯铀镭、纯钍源在两个不同能窗中的计数率之比，分别采用α、β、γ、a、b和g来表示，其中α为纯钍源在铀窗和钍窗上的计数率之比；β为纯钍源在钾窗和钍窗上的计数率之比；γ为纯铀镭源在钾窗和铀窗上的计数率之比；a为纯铀镭源在钍窗和铀窗中的计数率之比，b为纯40K源在钍窗和钾窗上的计数率之比；g为纯40K源在铀窗和钾窗上的计数率之比。由于纯40K源只辐射1.46MeV一种能量γ射线，因此纯40K源在铀窗和钍窗中无计数，康普顿散射系数b和g均为0，其他康普顿散射系数计算公式如下：
其中nK,U、nU,U、nTh,U分别为航空模型坪纯铀镭源在钾窗、铀窗和钍窗上产生的计数率，单位为计数每秒(s-1)；nK,Th、nU,Th、nTh,Th分别为航空模型坪纯钍源在钾窗、铀窗和钍窗上产生的计数率，单位为计数每秒(s-1)。
各能窗的计数率为NaI探测器和CeBr3探测器能窗计数率之和，为将两类探测器的能窗数据融合，NaI(Tl)探测器各能窗计数率视为能窗中全能峰计数率和康普顿散射本底计数率之和，其中纯铀镭源在NaI(Tl)探测器的钾窗和钍窗上无全能峰计数，只有康普顿散射本底计数；纯钍源在NaI(Tl)探测器的钾窗和铀窗上无全能峰计数，只有康普顿散射本底计数。求解康普顿散射系数所需各能窗计数率计算公式36-41如下：
公式36-41中为航空模型坪纯铀镭源在NaI(Tl)探测器的铀窗中产生的全能峰计数率，单位为计数每秒(s-1)；/>为航空模型坪纯钍源在NaI(Tl)探测器的钍窗中产生的全能峰计数率，单位为计数每秒(s-1)；/>分别为航空模型坪纯铀镭源在NaI(Tl)探测器的钾窗、铀窗和钍窗中产生的康普顿散射本底计数率，单位为计数每秒(s-1)；/>分别为航空模型坪纯钍源在NaI(Tl)探测器的钾窗、铀窗和钍窗中产生的康普顿散射本底计数率，单位为计数每秒(s-1)；/>分别为航空模型坪纯铀镭源在CeBr3探测器的钾窗、铀窗和钍窗中产生的计数率，单位为计数每秒(s-1)； 分别为航空模型坪纯钍源在CeBr3探测器的钾窗、铀窗和钍窗中产生的计数率，单位为计数每秒(s-1)。
1)航空模型坪纯铀镭源在NaI(Tl)探测器的铀窗中产生的全能峰计数率计算方法为公式42：/>
式中和/>为2.1中拟合公式15中系数，/>和/>为2.2中公式28和公式29确定的NaI(Tl)探测器的铀窗的起止道址。
2)航空模型坪纯钍源在NaI(Tl)探测器的钍窗中全能峰的计数率计算方法为公式43：
式中和/>为2.1中拟合公式16中系数，/>和/>为2.2中公式30和公式31确定的NaI(Tl)探测器的钍窗的起止道址。
3)航空模型坪纯铀镭源在NaI(Tl)探测器的铀窗中产生的康普顿散射本底计数率计算方法为公式44：
式中为2.1中拟合公式15中系数，/>和/>为2.2中公式22和公式23确定的CeBr3探测器的铀窗的起止道址。
4)航空模型坪纯钍源在NaI(Tl)探测器的钍窗中康普顿散射本底计数率计算方法为公式45：
式中为2.1中拟合公式16中系数，/>和/>为2.2中公式24和公式25确定的CeBr3探测器的钍窗的起止道址。
5)航空模型坪纯铀镭源在NaI(Tl)探测器的钾窗中康普顿散射本底计数率计算方法为公式46：
式中为2.1中公式8计算所得的NaI(Tl)探测器纯铀镭系的512道γ能谱中第j道计数率，/>和/>为2.2中公式20和公式21确定的CeBr3探测器的钾窗的起止道址。
6)航空模型坪纯铀镭源在NaI(Tl)探测器的钍窗中康普顿散射本底计数率计算方法为公式47：
/>
式中为2.1中公式8计算所得的NaI(Tl)探测器纯铀镭系的512道γ能谱中第j道计数率，/>和/>为2.2中公式24和公式25确定的CeBr3探测器的钍窗的起止道址。
7)航空模型坪纯钍源在NaI(Tl)探测器的钾窗中康普顿散射本底的计数率计算方法为公式48：
式中为2.1中公式9计算所得的NaI(Tl)探测器纯钍系的512道γ能谱中第j道计数率，/>和/>为2.2中公式20和公式21确定的CeBr3探测器的钾窗的起止道址。
8)航空模型坪纯钍源在NaI(Tl)探测器的铀窗中康普顿散射本底的计数率计算方法为公式49：
式中为2.1中公式9计算所得的NaI(Tl)探测器纯钍系的512道γ能谱中第j道计数率，/>和/>为2.2中公式22和公式23确定的CeBr3探测器的铀窗的起止道址。
9)在航空模型坪纯铀镭源上CeBr3探测器铀窗计数率计算方法为公式50：
式中和/>为2.1中拟合公式18中系数，/>和/>为2.2中公式22和公式23确定的CeBr3探测器的铀窗的起止道址。
10)在航空模型坪纯钍源上CeBr3探测器钍窗计数率计算方法为公式51：
式中和/>为2.1中拟合公式19中系数，/>和/>为2.2中公式24和公式25确定的CeBr3探测器的钍窗的起止道址。
11)在航空模型坪纯铀镭源上CeBr3探测器钾窗计数率计算方法为公式52：
式中为2.1中公式11计算所得的CeBr3探测器纯铀镭系的512道γ能谱中第j道计数率，/>和/>为2.2中公式20和公式21确定的CeBr3探测器的钾窗的起止道址。
12)在航空模型坪纯铀镭源上CeBr3探测器钍窗计数率计算方法为公式53：/>
式中为2.1中公式11计算所得的CeBr3探测器纯铀镭系的512道γ能谱中第j道计数率，/>和/>为2.2中公式24和公式25确定的CeBr3探测器的钍窗的起止道址。
13)在航空模型坪纯钍源上CeBr3探测器钾窗计数率计算方法为公式54：
式中为2.1中公式12计算所得的CeBr3探测器纯钍系的512道γ能谱中第j道计数率，/>和/>为2.2中公式20和公式21确定的CeBr3探测器的钾窗的起止道址。
14)在航空模型坪纯钍源上CeBr3探测器铀窗计数率计算方法为公式55：
式中为2.1中公式12计算所得的CeBr3探测器纯钍系的512道γ能谱中第j道计数率，/>和/>为2.2中公式22和公式23确定的CeBr3探测器的铀窗的起止道址。
四、飞机本底及宇宙射线系数计算
飞机本底和宇宙射线系数的校准采用3000m、3300m、3600m、3900m和4200m五个高高度飞行，每个高高度飞行10分钟，获取的能谱数据来计算。
对每个高高度上NaI(Tl)探测器获得512道能谱数据进行累加平均获得每个高度上归一到单个测量周期(1秒)的512道γ能谱计数率(i表示道址，i＝1,2,3,...,511,512；h表示高度，h＝3000,3300,3600,3900,4200)，单位为计数每秒(s-1)。
对每个高高度上CeBr3探测器获得512道能谱数据进行累加平均获得每个高高度上归一到单个测量周期(1秒)的512道γ能谱计数率(i表示道址，i＝1,2,3,...,511,512；h表示高度，h＝3000,3300,3600,3900,4200)，单位为计数每秒(s-1)。
4.1每个高高度上NaI(Tl)探测器512道能谱计数率降噪处理方法(如图8)
分别将每个高高度上NaI(Tl)探测器的512道γ能谱计数率视为采样频率为1000Hz的数字信号，采用通带频率为50Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得滤波后的512道γ能谱计数率/>(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)。
为实现滤波无相位位移，保证特征峰峰位无位移，对滤波后获得的512道γ能谱计数率进行转置，获取512道γ能谱计数率/>(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，单位为计数每秒(s-1)，见公式56。
将512道γ能谱计数率视为采样频率为1000Hz的数字信号，采用通带频率为50Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得滤波后的512道γ能谱计数率ri h，单位为计数每秒(s-1)。
512道γ能谱计数率ri h进行转置，获得各高高度上降噪滤波后的512道γ能谱计数率单位为计数每秒(s-1)，见公式57。
4.2每个高高度上NaI(Tl)探测器512道能谱康普顿散射本底计算
将每个高高度NaI(Tl)探测器的512道γ能谱计数率采用4.1相似的处理流程，仅在FIR滤波器参数上有差别，选择通带频率为30Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得新的512道γ能谱计数率/>(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，单位为计数每秒(s-1)，处理流程见图9。
将每个高高度层上NaI(Tl)探测器的512道γ能谱计数率采用4.1相似的处理流程，仅在FIR滤波器参数上有差别，选择通带频率为10Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得每个高高度层上的512道γ能谱计数率/>(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，处理流程见图10。
每个高高度层取的50至210道的能谱计数率数据作为对应道址的康普顿散射本底计数，取/>的211至512道的能谱计数率数据作为对应道址的本底计数，构成每个高高度层从50道至512道的γ能谱康普顿散射本底计数率数据，公式58和59：
如图11所示，以道址作为自变量x，以各道上康普顿散射本底计数率作为应变量y，分别对每个高高度层NaI(Tl)晶体探测器的50道至512道的γ能谱散射本底计数率采用幂函数y＝w·xm+t进行最小二乘拟合，获得拟合公式60，作为50道至512道康普顿散射本底描述函数。
其中wh、mh、th为各高度数据拟合获得参数，其中h表示高度，h＝3000,3300,3600,3900,4200。
4.3每个高高度上NaI(Tl)探测器钾窗、铀窗、钍窗计数率计算
NaI(Tl)探测器各能窗计数率为能窗中全能峰计数率和康普顿散射本底计数率之和，其中各能窗康普顿散射本底的范围采用CeBr3探测器所确定能窗范围。
1)各高高度层上NaI(Tl)探测器钾窗计数率计算方法为公式61：
式61中为高度为h的NaI(Tl)探测器钾窗计数率；/>为高度为h的NaI(Tl)探测器钾窗中全能峰计数率，单位为计数每秒(s-1)；/>为高度为h的NaI(Tl)探测器钾窗中康普顿散射本底计数率；/>为4.1中公式57获取降噪滤波后高度h的NaI(Tl)探测器第i道γ能谱计数率；wh、mh、th为4.2中高度为h的NaI(Tl)探测器50道至512道的γ能谱散射本底计数率幂函y＝w·xm+t拟合所得参数；/>分别为2.2中确定的NaI(Tl)探测器钾窗范围起止道址，/>分别为2.2中确定的CeBr3探测器钾窗范围起止道址。
(2)各高高度层上NaI(Tl)探测器铀窗计数率计算方法为公式62：
式62中为高度为h的NaI(Tl)探测器铀窗计数率，单位为计数每秒(s-1)；为高度为h的NaI(Tl)探测器铀窗中全能峰计数率，单位为计数每秒(s-1)；/>为高度为h的NaI(Tl)探测器铀窗中康普顿散射本底计数率；/>为4.1中公式57获取降噪滤波后高度为h的NaI(Tl)探测器的第i道γ能谱计数率；wh、mh、th为4.2中高度为h的NaI(Tl)探测器50道至512道的γ能谱散射本底计数率幂函y＝w·xm+t拟合所得参数；/>分别为2.2中确定的NaI(Tl)探测器铀窗范围起止道址，/>分别为2.2中确定的CeBr3探测器铀窗范围起止道址。
3)各高高度层上NaI(Tl)探测器钍窗计数率计算方法为公式63：
式63中为高度为h的NaI(Tl)探测器钍窗计数率；/>为高度为h的NaI(Tl)探测器钍窗中全能峰计数率；/>为高度为h的NaI(Tl)探测器钍窗中康普顿散射本底计数率；/>为4.1中公式57获取降噪滤波后高度为h的NaI(Tl)探测器的第i道γ能谱计数率；wh、mh、th为4.2中高度为h的NaI(Tl)探测器50道至512道的γ能谱散射本底计数率幂函y＝w·xm+t拟合所得参数；/>分别为2.2中确定的NaI(Tl)探测器钍窗范围起止道址，/>分别为2.2中确定的CeBr3探测器钍窗范围起止道址。
4.4每个高高度上CeBr3探测器各能窗计数率数据计算
分别将每个高高度上CeBr3探测器的512道γ能谱计数率数据采用4.1中NaI(Tl)探测器512道能谱计数率降噪滤波计算方法获得CeBr3探测器降噪滤波后的512道γ能谱计数率(h表示高度，h＝3000,3300,3600,3900,4200；i表示道址，i＝1,2,3,...,511,512)，单位为计数每秒(s-1)。
(1)各高高度层上CeBr3探测器钾窗计数率计算方法为公式64：
式64中为高度为h的CeBr3探测器钾窗计数率，单位为计数每秒(s-1)；分别为2.2中确定的CeBr3探测器钾窗范围起止道址。/>
(2)各高高度层上CeBr3探测器铀窗计数率计算方法为公式65：
式65中为高度为h的CeBr3探测器铀窗计数率，单位为计数每秒(s-1)；分别为2.2中确定的CeBr3探测器铀窗范围起止道址。
(3)各高高度层上CeBr3探测器钍窗计数率计算方法为公式66：
式66中为高度为h的CeBr3探测器钍窗计数率，单位为计数每秒(s-1)；分别为2.2中确定的CeBr3探测器钍窗范围起止道址。
4.5每个高高度上整套测量系统各能窗计数率计算
(1)整套测量系统在高高度h上钾窗计数率计算公式为式67
其中为高度为h的NaI(Tl)探测器钾窗计数率，单位为计数每秒(s-1)，/>为高度为h的CeBr3探测器钾窗计数率，单位为计数每秒(s-1)。
(2)整套测量系统在高高度h上铀窗计数率计算公式为式68
其中为高度为h的NaI(Tl)探测器铀窗计数率，单位为计数每秒(s-1)，/>为高度为h的CeBr3探测器铀窗计数率，单位为计数每秒(s-1)。
(3)整套测量系统在高高度h上钍窗计数率计算公式为式69
其中为高度为h的NaI(Tl)探测器钍窗计数率，单位为计数每秒(s-1)，/>为高度为h的CeBr3探测器钍窗计数率，单位为计数每秒(s-1)。
4.6飞机本底及宇宙射线系数
在进行高高度飞行测量时，宇宙射线道中的计数率与各能窗的计数率呈线性关系(公式70-72)：
分别为4.5中整套测量系统在高高度h上钾窗、铀窗、钍窗的平均计数率，单位为计数每秒(s-1)；/>为在高高度h上宇宙射线道的平均计数率，单位为计数每秒(s-1)；aK和bK为宇宙射线道平均计数率与钾窗平均计数率呈线性关系的系数；aU和bU为宇宙射线道平均计数率与铀窗平均计数率呈线性关系的系数；aTh和bTh为宇宙射线道平均计数率与钍窗平均计数率呈线性关系的系数。
对五个高高度层(3000、3300、3600、3900和4200)的宇宙射线道的平均计数率和能窗计数率进行最小二乘法线性拟合，获得上述线性关系系数，即aK为钾窗飞机本底计数率，单位为计数每秒(s-1)；bK为钾窗宇宙射线影响系数，无纲量；aU为铀窗飞机本底计数率，单位为计数每秒(s-1)；bU为铀窗宇宙射线影响系数，无纲量；aTh为钍窗飞机本底计数率，单位为计数每秒(s-1)；bTh为钍窗宇宙射线影响系数，无纲量。
五、高度衰减系数及空中灵敏度计算
通过动态测试带来确定高度衰减系数和空中窗灵敏度。动态测试是在动态测试带上空30m、60m、90m、120m、150m、210m、240m、300m的8个不同高度飞行测试，动态测试包括水域和相邻的陆地两部分。
5.1NaI(Tl)探测器各能窗计数率计算
对动态带各个高度上NaI(Tl)探测器获得的水域部分和陆域部分512道谱数据分别进行累加平均获得每个高度水域上512道γ能谱平均计数率和陆域上512道γ能谱数据计数率/> 表示NaI(Tl)探测器在动态带水域部分上空飞行高度为h获得的第i道γ能谱计数率，单位为计数每秒(s-1)，/>表示NaI(Tl)探测器在动态带陆域部分上空飞行高度h获得的第i道γ能谱计数率，单位为计数每秒(s-1)，其中h＝30、60、90、120、150、210、240、300，i＝1、2、3、...、511、512。
和/>作差为NaI(Tl)探测器在高度h上获得动态带纯陆地放射性核素产生的512道γ能谱中第i道的平均计数率/>单位为计数每秒(s-1)：
将每个高度层上的512道γ能谱计数率采用4.1至4.3相同的处理方法，获得NaI(Tl)探测器在每个高度上动态带纯陆地放射性核素产生的γ能谱各能窗计数率：
为NaI(Tl)探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钾窗计数率，单位为计数每秒(s-1)；
为NaI(Tl)探测器在高度h上动态带纯陆地放射性核素产生的γ能谱铀窗计数率，单位为计数每秒(s-1)；
为NaI(Tl)探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钍窗计数率，单位为计数每秒(s-1)。
5.2CeBr3探测器各能窗计数率计算
对动态带各个高度上CeBr3探测器获得的水域部分和陆域部分512道谱数据进行累加平均分别获得每个高度水域上512道γ能谱平均计数率和陆域上512道γ能谱数据计数率/> 表示CeBr3探测器在动态带水域部分上空飞
行高度为h获得的第i道γ能谱计数率，表示CeBr3探测器在动态带陆域部分上空飞行高度h获得的第i道γ能谱计数率，单位为计数每秒(s-1)，其中h＝30、60、90、120、150、210、240、300，i＝1、2、3、...、511、512。/>
和/>作差为CeBr3探测器飞行高度h获得动态带纯陆地放射性核素产生的512道γ能谱中第i道的平均计数率/>单位为计数每秒(s-1)：
将每个高度层上的512道γ能谱计数率采用4.4相同的处理方法，获得CeBr3探测器在每个高度上动态带纯陆地放射性核素产生的γ能谱各能窗计数率：
为CeBr3探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钾窗计数率，单位为计数每秒(s-1)；
为CeBr3探测器在高度h上动态带纯陆地放射性核素产生的γ能谱铀窗计数率，单位为计数每秒(s-1)；
为CeBr3探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钍窗计数率，单位为计数每秒(s-1)。
5.3整套测量系统各能窗计数率计算
整套测量系统在高度h上动态带纯陆地放射性核素产生的γ能谱钾窗计数率计算公式为式75，其中h＝30、60、90、120、150、210、240、300。
整套测量系统在高度h上动态带纯陆地放射性核素产生的γ能谱铀窗计数率计算公式为式76，其中h＝30、60、90、120、150、210、240、300。
整套测量系统在高度h上动态带纯陆地放射性核素产生的γ能谱钍窗计数率计算公式为式77，其中h＝30、60、90、120、150、210、240、3。
5.4高度衰减系数和空中灵敏度计算
整套测量系统在动态带各高度上纯陆地放射性核素产生的各能窗计数率进行康普顿散射改正公式为式78-80：
其中分别为系统在高度h上动态带纯陆地放射性核素产生的经过康普顿散射改正后的钾窗、铀窗和钍窗计数率，单位为计数每秒(s-1)；α、β、γ、a为步骤三中获得的康普顿散射系数，无纲量。
动态带纯陆地放射性核素产生的经过康普顿散射改正后的钾窗计数率铀窗计数率/>和钍窗计数率/>随高度h(单位为米)的变换呈e指数关系变化：/>
其中分别为空中120m高度纯陆地放射性核素产生的经过康普顿散射改正后的钾窗计数率、铀窗计数率和钍窗计数率，单位为计数每秒(s-1)；μK、μU、μTh分别为钾窗计数率、铀窗计数率和钍窗计数率高度衰减系数，单位为每米(m-1)；
利用30m、60m、90m、120m、150m、210m、240m、300m八个高度上获取的经过康普顿散射改正后的钾窗计数率、铀窗计数率和钍窗计数率，单位为计数每秒(s-1)，分别以经过康普顿散射改正后的各能窗计数率为因变量，以高度h为自变量，按e指数公式拟合就可求得各能量窗的高度衰减系数μK、μU和μTh，单位为每米(m-1)，以及120m高度纯陆地放射性核素产生的经过康普顿散射改正后的钾窗计数率铀窗计数率/>和钍窗计数率/>
5.5空中灵敏度计算
120m高度空中系统各能窗灵敏度计算公式如下：
式中和/>分别为钾窗、铀窗和钍窗在空中120m高度的灵敏度，单位分别为1/(s·％)、1/(s·ug/g)和1/(s·ug/g)；/>和/>分别为系统在高度120m动态带上纯陆地放射性核素产生的经过康普顿散射改正后的钾窗、铀窗和钍窗计数率，单位为计数每秒(s-1)；qK、qU和qTh分别为已知的动态带陆地部分钾(40K)、铀(238U)和钍(232Th)的含量，单位分别为％、ug/g和ug/g。
六、测点数据的处理
航空γ能谱测量过程中获得每个测点的NaI(Tl)探测器512道γ能谱计数率数据和CeBr3探测器512道γ能谱计数率数据/>其中i表示道址，i＝1,2,3,...,511,512。
6.1NaI(Tl)探测器各能窗数据计算
每个测点的NaI(Tl)探测器512道γ能谱计数率采用与4.1至4.3中每个高高度层上NaI(Tl)探测器的512道γ能谱计数率/>相同的处理方法，获得NaI(Tl)探测器在每个测点上各能窗计数率：
为NaI(Tl)探测器在每个测点上获得的钾窗计数率，单位为计数每秒(s-1)；
为NaI(Tl)探测器在每个测点上获得的铀窗计数率，单位为计数每秒(s-1)；/>
为NaI(Tl)探测器在每个测点上获得的钍窗计数率，单位为计数每秒(s-1)。
6.2CeBr3探测器各能窗计数率计算
将每个测点的CeBr3探测器512道γ能谱计数率采用与4.4中每个高度层上CeBr3探测器的512道γ能谱计数率/>相同的处理方法，获得CeBr3探测器在每个测点上各能窗计数率：
为CeBr3探测器在测点上获得的钾窗计数率，单位为计数每秒(s-1)；
为CeBr3探测器在测点上获得的铀窗计数率，单位为计数每秒(s-1)；
为CeBr3探测器在测点上获得的钍窗计数率，单位为计数每秒(s-1)。
6.3整套测量系统各能窗计数率计算
整套测量系统在各测点上获得的钾窗计数率nK计算公式为式87：
整套测量系统在各测点上获得的铀窗计数率nU计算公式为式88：
整套测量系统在各测点上获得的钍窗计数率nTh计算公式为式89：
6.4飞机本底及宇宙射线本底改正
测点上各能窗计数率数据飞机本底及宇宙射线本底改正计算方法如下:
其中和/>分别为测点上经飞机本底及宇宙射线本底改正后钾窗、铀窗和钍窗计数率，单位为计数每秒(s-1)；nK、nU和nTh分别为整套测量系统在各测点上获得的钾窗、铀窗和钍窗计数率，单位为计数每秒(s-1)；aK、aU和aTh分别为4.6中获得的钾窗、铀窗和钍窗的飞机本底计数率，单位为计数每秒(s-1)，bK、bU和bTh分别为4.6中获得的钾窗、铀窗和钍窗的宇宙射线系数，无量纲；NCos为测点上获得的宇宙射线道计数率，单位为计数每秒(s-1)。
6.5康普顿散射改正
测点上各能窗计数率数据康普顿散射改正计算方法如下:
其中分别为测点上经飞机本底、宇宙射线本底和康普顿散射改正后的钾窗、铀窗和钍窗计数率，单位为计数每秒(s-1)；α、β、γ、a、b和g为步骤三中获得的康普顿散射系数，无量纲。
6.6高度衰减改正
将测点上各能窗计数率高度改正到120m高度计数率的计算方法如下:
其中分别为经飞机本底、宇宙射线本底和康普顿散射改正及高度改正到空中120米高度的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)；
分别为在测点上空实际飞行高度h上经飞机本底、宇宙射线本底和康普顿散射改正的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)；
μK、μU、μTh分别为2.2.6.4中获得的钾窗计数率、铀窗计数率和钍窗计数率高度衰减系数，单位为每米(m-1)；h为该测点的飞行高度，单位为米(m)。
6.7放射性核素含量计算
经过上述各项修正后，每个测点地面上钾、铀、钍含量的计算公式如下：
QK QU QTh分别为测点处钾(40K)、铀(238U)和钍(232Th)的含量，单位分别为％、ug/g和ug/g；和/>分别为钾窗、铀窗和钍窗在空中120m高度的灵敏度，单位分别为1/(s·％)、1/(s·ug/g)和1/(s·ug/g)；/>分别为经飞机本底、宇宙射线本底和康普顿散射改正及高度改正到空中120米高度的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)。
本发明再一个实施例中，提供了一种航空γ能谱数据融合处理系统，该系统能够用于实施上述的一种航空γ能谱数据融合处理方法，具体的，包括：数据融合模块、天然核素特征峰窗口范围确定模块、康普顿散射系数确定模块、计算模块和测点数据处理模块；
数据融合模块：获取的钾窗、铀窗和钍窗的计数率；
获取NaI(Tl)探测器各特征峰能窗计数率。为将NaI(Tl)探测器与CeBr3探测器的能窗数据融合后统一处理，NaI(Tl)探测器的康普顿散射本底各能窗范围采用CeBr3探测器所确定能窗范围；
天然核素特征峰窗口范围确定模块：用于确定天然核素特征峰窗口范围；
康普顿散射系数确定模块：用于确定康普顿散射系数；
计算模块：用于计算飞机本底、宇宙射线系数、高度衰减系数和空中灵敏度；
测点数据处理模块；用于测点数据的处理。将两类探测器各能窗计数率数据进行合并，再采用统一的改正系数对合并后能窗数据进行处理，获取最终的天然核素含量值。
本发明再一个实施例中，提供了一种终端设备，该终端设备包括处理器以及存储器，所述存储器用于存储计算机程序，所述计算机程序包括程序指令，所述处理器用于执行所述计算机存储介质存储的程序指令。处理器可能是中央处理单元(Central ProcessingUnit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor、DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable GateArray，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等，其是终端的计算核心以及控制核心，其适于实现一条或一条以上指令，具体适于加载并执行一条或一条以上指令从而实现相应方法流程或相应功能；本发明实施例所述的处理器可以用于上述一种航空γ能谱数据融合处理方法的操作。
本发明再一个实施例中，本发明还提供了一种存储介质，具体为计算机可读存储介质(Memory)，所述计算机可读存储介质是终端设备中的记忆设备，用于存放程序和数据。可以理解的是，此处的计算机可读存储介质既可以包括终端设备中的内置存储介质，当然也可以包括终端设备所支持的扩展存储介质。计算机可读存储介质提供存储空间，该存储空间存储了终端的操作系统。并且，在该存储空间中还存放了适于被处理器加载并执行的一条或一条以上的指令，这些指令可以是一个或一个以上的计算机程序(包括程序代码)。需要说明的是，此处的计算机可读存储介质可以是高速RAM存储器，也可以是非不稳定的存储器(non-volatile memory)，例如至少一个磁盘存储器。
可由处理器加载并执行计算机可读存储介质中存放的一条或一条以上指令，以实现上述实施例中有关一种航空γ能谱数据融合处理方法的相应步骤；计算机可读存储介质中的一条或一条以上指令由处理器加载并执行。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的实施方法，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。
Claims (9)
1.一种航空γ能谱数据融合处理方法，其特征在于，包括以下步骤：
一、数据融合
系统获取的钾窗、铀窗和钍窗的计数率分别为同步获取的NaI(Tl)探测器获取的相应能窗计数率与CeBr3探测器获取的相应能窗计数率之和；
NaI(Tl)探测器各特征峰能窗计数率由能窗中全能峰计数率和康普顿散射本底计数率构成；为将NaI(Tl)探测器与CeBr3探测器的能窗数据融合后统一处理，NaI(Tl)探测器的康普顿散射本底各能窗范围采用CeBr3探测器所确定能窗范围；
各航空模型包括：本底模型、钾模型、铀模型、钍模型；
二、确定天然核素特征峰窗口范围；
2.1特征峰曲线拟合
在用于检定航空γ能谱仪的本底模型、钾模型、铀模型、钍模型上分别进行测量，同步获得NaI(Tl)和CeBr3探测器两种512道γ能谱数据，进行平均处理获得归一到单个测量周期的512道能谱数据，即512道γ能谱计数率；分别对NaI(Tl)和CeBr3两种探测器测得的钾模型上的512道能谱计数率扣除本底模型上的512道能谱计数率获得纯40K的512道γ能谱计数率；分别对NaI(Tl)和CeBr3两种探测器测得的铀模型上的512道能谱计数率扣除本底模型上的512道能谱计数率获得纯铀镭系的512道γ能谱计数率；分别对NaI(Tl)和CeBr3两种探测器测得的钍模型上的512道能谱计数率扣除本底模型上的512道能谱计数率获得纯钍系的512道γ能谱计数率；
用线性函数加高斯函数来描述512道γ能谱中40K的特征峰、214Bi的特征峰和208Tl的特征峰；
以道址作为自变量x，以各道上的计数率作为应变量y，分别对NaI(Tl)探测器和CeBr3探测器获取的纯40K 512道γ能谱计数率中的特征峰、纯铀镭系512道γ能谱计数率中的特征峰、纯钍系512道γ能谱计数率中的特征峰采用线性函数进行最小二乘拟合；
2.2确定特征峰能窗范围
确定CeBr3探测器和NaI(Tl)探测器的航空γ能谱数据各能窗范围，分别为:钾窗的起始道址，钾窗的结束道址；铀窗的起始道址，铀窗的结束道址；钍窗的起始道址，钍窗的结束道址；
三、确定康普顿散射系数；
航空γ能谱数据处理的能窗分别为钾窗、铀窗和钍窗，康普顿散射剥离系数是纯40K、纯铀镭、纯钍源在两个不同能窗中的计数率之比，分别采用α、β、γ、a、b和g来表示，其中α为纯钍源在铀窗和钍窗上的计数率之比；β为纯钍源在钾窗和钍窗上的计数率之比；γ为纯铀镭源在钾窗和铀窗上的计数率之比；a为纯铀镭源在钍窗和铀窗中的计数率之比，b为纯40K源在钍窗和钾窗上的计数率之比；g为纯40K源在铀窗和钾窗上的计数率之比；康普顿散射系数b和g均为0，计算剩余康普顿散射系数，公式如下：
其中nK,U、nU,U、nTh,U分别为航空模型坪纯铀镭源在钾窗、铀窗和钍窗上产生的计数率，单位为计数每秒(s-1)；nK,Th、nU,Th、nTh,Th分别为航空模型坪纯钍源在钾窗、铀窗和钍窗上产生的计数率，单位为计数每秒(s-1)；
各能窗的计数率为NaI(Tl)探测器和CeBr3探测器能窗计数率之和，为将两类探测器的能窗数据融合，NaI(Tl)探测器各能窗计数率视为能窗中全能峰计数率和康普顿散射本底计数率之和，其中纯铀镭源在NaI(Tl)探测器的钾窗和钍窗上无全能峰计数，只有康普顿散射本底计数；纯钍源在NaI(Tl)探测器的钾窗和铀窗上无全能峰计数，只有康普顿散射本底计数；求解康普顿散射系数所需各能窗计数率计算公式36-41如下：
公式中为航空模型坪纯铀镭源在NaI(Tl)探测器的铀窗中产生的全能峰计数率，单位为计数每秒(s-1)；/>为航空模型坪纯钍源在NaI(Tl)探测器的钍窗中产生的全能峰计数率，单位为计数每秒(s-1)；/>分别为航空模型坪纯铀镭源在NaI(Tl)探测器的钾窗、铀窗和钍窗中产生的康普顿散射本底计数率，单位为计数每秒(s-1)；分别为航空模型坪纯钍源在NaI(Tl)探测器的钾窗、铀窗和钍窗中产生的康普顿散射本底计数率，单位为计数每秒(s-1)；/>分别为航空模型坪纯铀镭源在CeBr3探测器的钾窗、铀窗和钍窗中产生的计数率，单位为计数每秒(s-1)；分别为航空模型坪纯钍源在CeBr3探测器的钾窗、铀窗和钍窗中产生的计数率，单位为计数每秒(s-1)；
四、飞机本底及宇宙射线系数计算；
飞机本底和宇宙射线系数的校准采用五个高高度飞行，每个高高度飞行10分钟，获取的能谱数据来计算；
对每个高高度层上NaI(Tl)探测器获得512道能谱数据进行累加平均获得每个高高度上归一到单个测量周期的512道γ能谱计数率i表示道址，i＝1,2,3,...,511,512；h表示高度，单位为计数每秒(s-1)；
对每个高高度层上CeBr3探测器获得512道能谱数据进行累加平均获得每个高高度上归一到单个测量周期的512道γ能谱计数率单位为计数每秒(s-1)；
4.1分别对每个高高度上NaI(Tl)探测器512道能谱计数率降噪处理；
4.2分别计算每个高高度上NaI(Tl)探测器512道能谱康普顿散射本底；
4.3分别计算每个高高度上NaI(Tl)探测器钾窗、铀窗、钍窗计数率；
4.4分别计算每个高高度上CeBr3探测器钾窗、铀窗、钍窗计数率数据；
4.5分别计算每个高高度上整套测量系统钾窗、铀窗和钍窗计数率；
计算整套测量系统在高高度h上钾窗计数率单位为计数每秒(s-1)，h表示高度，h＝3000,3300,3600,3900,4200；
计算整套测量系统在高高度h上铀窗计数率单位为计数每秒(s-1)，h表示高度，h＝3000,3300,3600,3900,4200；
计算整套测量系统在高高度h上钍窗计数率单位为计数每秒(s-1)，h表示高度，h＝3000,3300,3600,3900,4200；
4.6获取飞机本底及宇宙射线系数；
在所有高高度进行飞行测量，宇宙射线道中的计数率与各能窗的计数率呈线性关系；
对每个高高度层的宇宙射线道的平均计数率和能窗计数率进行最小二乘法线性拟合，获得线性关系系数，即钾窗飞机本底计数率，单位为计数每秒(s-1)；钾窗宇宙射线影响系数，无纲量；铀窗飞机本底计数率，单位为计数每秒(s-1)；铀窗宇宙射线影响系数，无纲量；钍窗飞机本底计数率，单位为计数每秒(s-1)；钍窗宇宙射线影响系数，无纲量；
五、高度衰减系数及空中灵敏度计算
5.1NaI(Tl)探测器各能窗计数率计算
对动态带各个高度上NaI(Tl)探测器获得的水域部分和陆域部分512道谱数据分别进行累加平均获得每个高度水域上512道γ能谱平均计数率和陆域上512道γ能谱数据计数率/>表示NaI(Tl)探测器在动态带水域部分上空飞行高度为h获得的第i道γ能谱计数率，，单位为计数每秒(s-1)，/>表示NaI(Tl)探测器在动态带陆域部分上空飞行高度h获得的第i道γ能谱计数率，单位为计数每秒(s-1)，其中h＝30、60、90、120、150、210、240、300，i＝1、2、3、...、511、512；
和/>作差为NaI(Tl)探测器在高度h上获得动态带纯陆地放射性核素产生的512道γ能谱中第i道的平均计数率/>单位为计数每秒(s-1)：
将每个高度层上的512道γ能谱计数率采用4.1至4.3相同的处理方法，获得NaI(Tl)探测器在每个高度上动态带纯陆地放射性核素产生的γ能谱各能窗计数率：
为NaI(Tl)探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钾窗计数率，单位为计数每秒(s-1)；
为NaI(Tl)探测器在高度h上动态带纯陆地放射性核素产生的γ能谱铀窗计数率，单位为计数每秒(s-1)；
为NaI(Tl)探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钍窗计数率，单位为计数每秒(s-1)；
5.2CeBr3探测器各能窗计数率计算
对动态带各个高度上CeBr3探测器获得的水域部分和陆域部分512道谱数据进行累加平均分别获得每个高度水域上512道γ能谱平均计数率和陆域上512道γ能谱数据计数率/>表示CeBr3探测器在动态带水域部分上空飞行高度为h获得的第i道γ能谱计数率，/>表示CeBr3探测器在动态带陆域部分上空飞行高度h获得的第i道γ能谱计数率，单位为计数每秒(s-1)，其中h＝30、60、90、120、150、210、240、300，i＝1、2、3、...、511、512；
和/>作差为CeBr3探测器飞行高度h获得动态带纯陆地放射性核素产生的512道γ能谱中第i道的平均计数率/>单位为计数每秒(s-1)：
将每个高度层上的512道γ能谱计数率采用4.4相同的处理方法，获得CeBr3探测器在每个高度上动态带纯陆地放射性核素产生的γ能谱各能窗计数率：
为CeBr3探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钾窗计数率，单位为计数每秒(s-1)；
为CeBr3探测器在高度h上动态带纯陆地放射性核素产生的γ能谱铀窗计数率，单位为计数每秒(s-1)；
为CeBr3探测器在高度h上动态带纯陆地放射性核素产生的γ能谱钍窗计数率，单位为计数每秒(s-1)；
5.3整套测量系统各能窗计数率计算
计算整套测量系统在高度h上动态带纯陆地放射性核素产生的γ能谱钾窗计数率
计算整套测量系统在高度h上动态带纯陆地放射性核素产生的γ能谱铀窗计数率
计算整套测量系统在高度h上动态带纯陆地放射性核素产生的γ能谱钍窗计数率
5.4高度衰减系数和空中灵敏度计算
在动态带各高度上纯陆地放射性核素产生的各能窗计数率进行康普顿散射改正；
利用30m、60m、90m、120m、150m、210m、240m、、300八个高度上获取的钾窗计数率、铀窗计数率和钍窗计数率，单位为计数每秒(s-1)，按e指数公式拟合就求得各能量窗的高度衰减系数μK、μU和μTh，单位为每米(m-1)；
5.5空中灵敏度计算
计算120m高度空中系统各能窗灵敏度；
六、测点数据的处理
航空γ能谱测量过程中获得每个测点的NaI(Tl)探测器512道γ能谱计数率数据和CeBr3探测器512道γ能谱计数率数据/>其中i表示道址，i＝1,2,3,...,511,512；
6.1NaI(Tl)探测器各能窗数据计算
每个测点的NaI(Tl)探测器512道γ能谱计数率采用与4.1至4.3中每个高度层上NaI(Tl)探测器的512道γ能谱计数率/>相同的处理方法，获得NaI(Tl)探测器在每个测点上各能窗计数率：
为NaI(Tl)探测器在每个测点上获得的钾窗计数率，单位为计数每秒(s-1)；
为NaI(Tl)探测器在每个测点上获得的铀窗计数率，单位为计数每秒(s-1)；
为NaI(Tl)探测器在每个测点上获得的钍窗计数率，单位为计数每秒(s-1)；
6.2CeBr3探测器各能窗计数率计算
将每个测点的CeBr3探测器512道γ能谱计数率采用与4.4中每个高度层上CeBr3探测器的512道γ能谱计数率/>相同的处理方法，获得CeBr3探测器在每个测点上各能窗计数率：
为CeBr3探测器在测点上获得的钾窗计数率，单位为计数每秒(s-1)；/>为CeBr3探测器在测点上获得的铀窗计数率，单位为计数每秒(s-1)；/>为CeBr3探测器在测点上获得的钍窗计数率，单位为计数每秒(s-1)；
6.3整套测量系统各能窗计数率计算
计算各测点上获得的钾窗计数率nK
计算各测点上获得的铀窗计数率nU
计算各测点上获得的钍窗计数率nTh
6.4飞机本底及宇宙射线本底改正
测点上各能窗计数率数据飞机本底及宇宙射线本底改正计算；
6.5康普顿散射改正
测点上各能窗计数率数据康普顿散射改正计算；
6.6高度衰减改正
计算将测点上各能窗计数率高度改正到120m高度计数率
6.7放射性核素含量计算
计算每个测点地面上钾、铀、钍含量，公式如下：
QK QU QTh分别为测点处钾、铀和钍的含量，单位分别为％、ug/g和ug/g；和/>分别为钾窗、铀窗和钍窗在空中120m高度的灵敏度，单位分别为1/(s·％)、1/(s·ug/g)和1/(s·ug/g)；/>分别为经飞机本底、宇宙射线本底和康普顿散射改正及高度改正到空中120米高度的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)。
2.根据权利要求1所述的一种航空γ能谱数据融合处理方法，其特征在于：步骤2.1中用线性函数加高斯函数来描述512道γ能谱中40K的1.46MeV特征峰、214Bi的1.76MeV特征峰和208Tl的2.62MeV特征峰，线性加高斯函数公式如下：
以道址作为自变量x，以各道上的计数率作为应变量y，分别对NaI(Tl)探测器和CeBr3探测器获取的纯40K 512道γ能谱计数率中的1.46MeV特征峰、纯铀镭系512道γ能谱计数率中的1.76MeV特征峰、纯钍系512道γ能谱计数率中的2.62MeV特征峰采用公式13表示的函数进行最小二乘拟合。
3.根据权利要求2所述的一种航空γ能谱数据融合处理方法，其特征在于：步骤2.2中特征峰能窗范围的确定
CeBr3探测器的航空γ能谱数据各能窗范围为:取整作为钾窗的起始道址取整作为钾窗的结束道址/>取整作为铀窗的起始道址/>取整作为铀窗的结束道址/> 取整作为钍窗的起始道址/>取整作为钍窗的结束道址/>
NaI(Tl)探测器的航空γ能谱数据各能窗范围为:取整作为钾窗的起始道址/>取整作为钾窗的结束道址/>取整作为铀窗的起始道址/>取整作为铀窗的结束道址/>取整作为钍窗的起始道址/>取整作为钍窗的结束道址/>
4.根据权利要求3所述的一种航空γ能谱数据融合处理方法，其特征在于：步骤三中：
3.1计算航空模型坪纯铀镭源在NaI(Tl)探测器的铀窗中产生的全能峰计数公式为：
式中和/>分别为采用高斯函数公式拟合NaI(Tl)探测器纯铀镭系512道γ能谱计数率中214Bi的1.76MeV特征峰的系数；，/>为铀窗的起始道址和/>为铀窗的结束道址；
3.2计算航空模型坪纯钍源在NaI(Tl)探测器的钍窗中全能峰的计数率公式为：
式中和/>为采用高斯函数公式拟合NaI(Tl)探测器纯钍系512道γ能谱计数率中208Tl的2.62MeV特征峰获得的系数；/>NaI(Tl)探测器的钍窗的起始道址；/>为NaI(Tl)探测器的钍窗的结束道址；
3.3计算航空模型坪纯铀镭源在NaI(Tl)探测器的铀窗中产生的康普顿散射本底计数率公式为：
式中为采用高斯函数公式拟合NaI(Tl)探测器纯铀镭系512道γ能谱计数率中214Bi的1.76MeV特征峰的系数；/>为CeBr3探测器的铀窗的起始道址，/>为CeBr3探测器的铀窗的结束道址；
3.4计算航空模型坪纯钍源在NaI(Tl)探测器的钍窗中康普顿散射本底计数率公式为：
式中为采用高斯函数公式拟合NaI(Tl)探测器纯钍系512道γ能谱计数率中208Tl的2.62MeV特征峰获得的系数；/>为CeBr3探测器的钍窗的起始道址，/>为CeBr3探测器的钍窗的结束道址；
3.5计算在航空模型坪纯铀镭源上CeBr3探测器铀窗计数率公式为：
式中和/>分别为采用高斯函数公式拟合CeBr3探测器获取纯铀镭系512道γ能谱计数率中214Bi的1.76MeV特征峰获得的系数，/>为CeBr3探测器的铀窗的起始道址，/>为CeBr3探测器的铀窗的结束道址；
3.6计算在航空模型坪纯钍源上CeBr3探测器钍窗计数率公式为：
式中和/>分别为采用高斯函数公式拟合CeBr3探测器获取纯钍系512道γ能谱计数率中208Tl的2.62MeV特征峰获得的系数，/>为CeBr3探测器的钍窗的起始道址和/>为CeBr3探测器的钍窗的结束道址。
5.根据权利要求4所述的一种航空γ能谱数据融合处理方法，其特征在于：4.1分别对每个高高度上NaI(Tl)探测器512道能谱计数率降噪处理，具体如下：
分别将每个高高度上NaI(Tl)探测器的512道γ能谱计数率视为采样频率为1000Hz的数字信号，采用通带频率为50Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得滤波后的512道γ能谱计数率/>h表示高度；i表示道址，i＝1,2,3,...,511,512；
为实现滤波无相位位移，保证特征峰峰位无位移，对滤波后获得的512道γ能谱计数率进行转置，获取512道γ能谱计数率/>h表示高度；i表示道址，单位为计数每秒(s-1)；
将512道γ能谱计数率视为采样频率为1000Hz的数字信号，采用通带频率为50Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得滤波后的512道γ能谱计数率/>单位为计数每秒(s-1)；
512道γ能谱计数率进行转置，获得各高高度上降噪滤波后的512道γ能谱计数率单位为计数每秒(s-1)。
6.根据权利要求5所述的一种航空γ能谱数据融合处理方法，其特征在于：4.2分别计算每个高高度上NaI(Tl)探测器512道能谱康普顿散射本底，具体如下：
将每个高高度NaI(Tl)探测器的512道γ能谱计数率采用采样频率为1000Hz的数字信号，采用通带频率为30Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得新的512道γ能谱计数率/>h表示高度；i表示道址，i＝1,2,3,...,511,512，单位为计数每秒(s-1)；
将每个高高度层上NaI(Tl)探测器的512道γ能谱计数率采用4.1相似的处理流程，仅在FIR滤波器参数上有差别，选择通带频率为10Hz、阻带频率为300Hz的FIR滤波器进行滤波，获得每个高度层上的512道γ能谱计数率/>
每个高高度层取的50至210道的能谱计数率数据作为对应道址的康普顿散射本底计数，取/>的211至512道的能谱计数率数据作为对应道址的本底计数，构成每个高高度层从50道至512道的γ能谱康普顿散射本底计数率数据：
以道址作为自变量x，以各道上康普顿散射本底计数率作为应变量y，分别对每个高高度层NaI(Tl)晶体探测器的50道至512道的γ能谱散射本底计数率采用幂函数y＝w·xm+t进行最小二乘拟合，获得拟合作为50道至512道康普顿散射本底描述函数。
7.根据权利要求6所述的一种航空γ能谱数据融合处理方法，其特征在于：4.3分别计算每个高高度上NaI(Tl)探测器钾窗、铀窗、钍窗计数率，具体如下：
4.3.1各高高度层上NaI(Tl)探测器钾窗计数率计算方法为公式：
式中为高度为h的NaI(Tl)探测器钾窗计数率；/>为高度为h的NaI(Tl)探测器钾窗中全能峰计数率，单位为计数每秒(s-1)；/>为高度为h的NaI(Tl)探测器钾窗中康普顿散射本底计数率；/>为获取降噪滤波后高度h的NaI(Tl)探测器第i道γ能谱计数率；wh、mh、th为高度为h的NaI(Tl)探测器50道至512道的γ能谱散射本底计数率幂函y＝w·xm+t拟合所得参数；/>为NaI(Tl)探测器钾窗范围起止道址，/>为CeBr3探测器钾窗范围起止道址；
4.3.2各高高度层上NaI(Tl)探测器铀窗计数率计算方法为公式：
式中为高度为h的NaI(Tl)探测器铀窗计数率，单位为计数每秒(s-1)；/>为高度为h的NaI(Tl)探测器铀窗中全能峰计数率，单位为计数每秒(s-1)；/>为高度为h的NaI(Tl)探测器铀窗中康普顿散射本底计数率；/>为降噪滤波后高度为h的NaI(Tl)探测器的第i道γ能谱计数率；wh、mh、th为高度为h的NaI(Tl)探测器50道至512道的γ能谱散射本底计数率幂函y＝w·xm+t拟合所得参数；/>分别为NaI(Tl)探测器铀窗范围起止道址，/>分别为CeBr3探测器铀窗范围起止道址；
4.3.3各高高度层上NaI(Tl)探测器钍窗计数率计算方法为公式：
式中为高度为h的NaI(Tl)探测器钍窗计数率；/>为高度为h的NaI(Tl)探测器钍窗中全能峰计数率；/>为高度为h的NaI(Tl)探测器钍窗中康普顿散射本底计数率；/>为降噪滤波后高度为h的NaI(Tl)探测器的第i道γ能谱计数率；wh、mh、th为高度为h的NaI(Tl)探测器50道至512道的γ能谱散射本底计数率幂函y＝w·xm+t拟合所得参数；分别为NaI(Tl)探测器钍窗范围起止道址，/>分别为CeBr3探测器钍窗范围起止道址。
8.根据权利要求7所述的一种航空γ能谱数据融合处理方法，其特征在于：所述步骤一中，系统获取的钾窗、铀窗和钍窗的计数率分别为同步获取的NaI(Tl)探测器获取的相应能窗计数率与CeBr3探测器获取的相应能窗计数率之和，计算公式如下：
nK、nU和nTh分别为数据融合后整套测量系统的航空γ能谱钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)；和/>分别为CeBr3探测器航空γ能谱的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)；/>和/>分别为NaI(Tl)探测器航空γ能谱的钾窗、铀窗和钍窗的计数率，单位为计数每秒(s-1)。
9.一种航空γ能谱数据融合处理系统，其特征在于：该系统能够用于实施权利要求1至8其中一项所述的一种航空γ能谱数据融合处理方法；
航空γ能谱数据融合处理系统包括：数据融合模块、天然核素特征峰窗口范围确定模块、康普顿散射系数确定模块、计算模块和测点数据处理模块；
数据融合模块：获取的钾窗、铀窗和钍窗的计数率；
获取NaI(Tl)探测器各特征峰能窗计数率；为将NaI(Tl)探测器与CeBr3探测器的能窗数据融合后统一处理，NaI(Tl)探测器的康普顿散射本底各能窗范围采用CeBr3探测器所确定能窗范围；
天然核素特征峰窗口范围确定模块：用于确定天然核素特征峰窗口范围；
康普顿散射系数确定模块：用于确定康普顿散射系数；
计算模块：用于计算飞机本底、宇宙射线系数、高度衰减系数和空中灵敏度；
测点数据处理模块；用于测点数据的处理；将两类探测器各能窗计数率数据进行合并，再采用统一的改正系数对合并后能窗数据进行处理，获取最终的天然核素含量值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310598009.6A CN116611026B (zh) | 2023-05-25 | 2023-05-25 | 一种航空γ能谱数据融合处理方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310598009.6A CN116611026B (zh) | 2023-05-25 | 2023-05-25 | 一种航空γ能谱数据融合处理方法和系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116611026A true CN116611026A (zh) | 2023-08-18 |
| CN116611026B CN116611026B (zh) | 2024-01-09 |
Family
ID=87685038
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310598009.6A Active CN116611026B (zh) | 2023-05-25 | 2023-05-25 | 一种航空γ能谱数据融合处理方法和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116611026B (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0704368D0 (en) * | 2007-03-07 | 2007-04-11 | Kokkas Nikolaos | Geodata fusion for automated 3D city modelling |
| CN101078776A (zh) * | 2007-06-29 | 2007-11-28 | 西安奥华电子仪器有限责任公司 | 核测井数据多尺度处理方法 |
| CN108825220A (zh) * | 2018-05-25 | 2018-11-16 | 东华理工大学 | 融合自然γ能谱与中子时间谱的铀矿测井仪及铀定量方程 |
| CN108825219A (zh) * | 2018-05-25 | 2018-11-16 | 东华理工大学 | 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法 |
| CN109630091A (zh) * | 2018-11-02 | 2019-04-16 | 中国石油天然气股份有限公司 | 一种碳氧比测井中优化数值模拟能谱的方法 |
| CN112083471A (zh) * | 2020-09-15 | 2020-12-15 | 成都理工大学 | 反康普顿航空伽马能谱测量系统及其测量方法 |
| FR3104737A1 (fr) * | 2019-12-11 | 2021-06-18 | Orano Mining | Procédé et système d’évaluation d’un paramètre représentatif de la concentration massique en uranium d’un échantillon de matériau uranifère par spectrométrie gamma |
| ES2877772A1 (es) * | 2020-05-14 | 2021-11-17 | Consejo Superior Investigacion | Dispositivo de deteccion, identificacion, cuantificacion y/o localizacion simultanea de fuentes de radiacion gamma y de neutrones |
| CN114236596A (zh) * | 2021-12-30 | 2022-03-25 | 成都理工大学 | 基于双模探测器系统的核废物包装体自适应扫描方法 |
| CN114386497A (zh) * | 2021-12-31 | 2022-04-22 | 核工业北京地质研究院 | 面向铀成矿构造的航空高光谱和伽玛能谱数据融合方法 |
| CN114895346A (zh) * | 2022-05-05 | 2022-08-12 | 中国自然资源航空物探遥感中心 | 适用于全地形的航空γ能谱测量系统及其测量方法 |
-
2023
- 2023-05-25 CN CN202310598009.6A patent/CN116611026B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0704368D0 (en) * | 2007-03-07 | 2007-04-11 | Kokkas Nikolaos | Geodata fusion for automated 3D city modelling |
| CN101078776A (zh) * | 2007-06-29 | 2007-11-28 | 西安奥华电子仪器有限责任公司 | 核测井数据多尺度处理方法 |
| CN108825220A (zh) * | 2018-05-25 | 2018-11-16 | 东华理工大学 | 融合自然γ能谱与中子时间谱的铀矿测井仪及铀定量方程 |
| CN108825219A (zh) * | 2018-05-25 | 2018-11-16 | 东华理工大学 | 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法 |
| CN109630091A (zh) * | 2018-11-02 | 2019-04-16 | 中国石油天然气股份有限公司 | 一种碳氧比测井中优化数值模拟能谱的方法 |
| FR3104737A1 (fr) * | 2019-12-11 | 2021-06-18 | Orano Mining | Procédé et système d’évaluation d’un paramètre représentatif de la concentration massique en uranium d’un échantillon de matériau uranifère par spectrométrie gamma |
| ES2877772A1 (es) * | 2020-05-14 | 2021-11-17 | Consejo Superior Investigacion | Dispositivo de deteccion, identificacion, cuantificacion y/o localizacion simultanea de fuentes de radiacion gamma y de neutrones |
| CN112083471A (zh) * | 2020-09-15 | 2020-12-15 | 成都理工大学 | 反康普顿航空伽马能谱测量系统及其测量方法 |
| CN114236596A (zh) * | 2021-12-30 | 2022-03-25 | 成都理工大学 | 基于双模探测器系统的核废物包装体自适应扫描方法 |
| CN114386497A (zh) * | 2021-12-31 | 2022-04-22 | 核工业北京地质研究院 | 面向铀成矿构造的航空高光谱和伽玛能谱数据融合方法 |
| CN114895346A (zh) * | 2022-05-05 | 2022-08-12 | 中国自然资源航空物探遥感中心 | 适用于全地形的航空γ能谱测量系统及其测量方法 |
Non-Patent Citations (2)
| Title |
|---|
| 毛鑫: "基于集成分类器的多路γ能谱测量软件开发", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
| 花永涛等: "基于数据融合的航空γ能谱特征峰信息提取技术研究", 《核技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116611026B (zh) | 2024-01-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Grasty et al. | A Guide to the Techncial Specifications for Airborne Gamma-ray Surveys | |
| Pitkin et al. | Design parameters for aerial gamma-ray surveys | |
| Minty | Airborne gamma-ray spectrometric background estimation using full spectrum analysis | |
| CN112083471B (zh) | 反康普顿航空伽马能谱测量系统及其测量方法 | |
| Germanenko et al. | High-energy photons connected to atmospheric precipitations | |
| Grasty et al. | The analysis of multichannel airborne gamma-ray spectra | |
| CN108318910A (zh) | 基于7Be峰的航空伽玛能谱仪稳谱方法 | |
| JP2016033459A (ja) | 放射性物質のモニタリング装置及び放射性物質のモニタリング方法 | |
| Srinivas et al. | Assessment of background gamma radiation levels using airborne gamma ray spectrometer data over uranium deposits, Cuddapah Basin, India–A comparative study of dose rates estimated by AGRS and PGRS | |
| CN112034533B (zh) | 无上测晶体大气氡修正参数求取方法 | |
| Letaw et al. | Satellite observation of atmospheric nuclear gamma radiation | |
| CN116611026B (zh) | 一种航空γ能谱数据融合处理方法和系统 | |
| CN104504273B (zh) | 一种提高γ射线峰探测灵敏度的能谱处理方法 | |
| Baranwal et al. | Mapping of caesium fallout from the Chernobyl accident in the Jotunheimen area | |
| Grasty | Direct snow-water equivalent measurement by air-borne gamma-ray spectrometry | |
| Fortin et al. | Airborne gamma-ray spectrometry in 2017: solid ground for new development | |
| CN106125126B (zh) | 采用溴化镧探测器测量环境中的钾40的方法 | |
| Baranwal et al. | Reprocessing of airborne gamma-ray spectrometry data in Norway for mapping of Cs-137 deposition from the Chernobyl accident | |
| Grasty et al. | Developments in the standardization and analysis of airborne gamma ray data | |
| Grasty et al. | Back calibration and reprocessing an airborne gamma-ray survey, Malaysia | |
| JP7378377B2 (ja) | 放射線分析装置およびダストモニタ装置 | |
| Jurza et al. | Use of 214Pb photopeaks for radon removal: utilising current airborne gamma-ray spectrometer technology and data processing | |
| Grasty | Snow-water equivalent measurement using natural gamma emission | |
| Guild-Bingham | Aerial Radiological Survey of Defense-Related Uranium Mines in the Shirley Basin, Crooks Gap, and Gas Hills Regions in Wyoming | |
| Grasty et al. | The standardisation of airborne gamma-ray surveys in Australia |
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 |