[发明专利]无上测晶体大气氡修正参数求取方法

权利要求书

1.一种无上测晶体大气氡修正参数求取方法，其特征在于，包括以下步骤：

a.在测试带区域内3000m及以上高度飞行多个不同高度，采集测试带区域内不同高度的放射性谱数据；

b.在测试带区域内宽广水面上空和陆地上空分别飞行多个不同高度，飞行高度小于300m，采集不同高度的放射性谱数据；

c.在航空放射性测量模型上，使用飞行器的航空放射性测量设备采集放射性谱数据；

d.利用步骤a采集的不同高度放射性谱数据求取低能铀窗飞机本底和宇宙射线修正系数；

e.根据步骤b采集的水面上空的放射性谱数据提取低能铀窗数据，利用步骤d求得的低能铀窗飞机本底和宇宙射线修正系数进行飞机本底与宇宙射线修正；

f.用上一步修正后的低能铀窗数据除以经过飞机本底和宇宙射线修正后的铀窗数据得到大气氡修正参数C1；

g.利用步骤b采集的放射性谱数据求取各飞行高度纯陆地核素谱数据；

h.根据上一步得到的各飞行高度纯陆地核素谱数据求取纯陆地核素谱数据的各道高度衰减系数；

i.利用步骤c采集的放射性谱数据求取钾、铀、钍各核素1m高度单位谱；

j.利用步骤h的结果和步骤i的结果求取各核素不同高度单位谱；

k.利用铀核素不同高度单位谱求取低能铀窗与铀窗数据之比即为不同高度大气氡修正参数C2，利用钍核素不同高度单位谱求取低能铀窗与钍窗数据之比即为不同高度大气氡修正参数C3。

2.根据权利要求1所述的无上测晶体大气氡修正参数求取方法，其特征在于，步骤a中，在空气清净时，在真高3000m以上，飞行5个或大于5个不同的高度，高度间隔为300m，每个高度飞行不少于10分钟。

3.根据权利要求1所述的无上测晶体大气氡修正参数求取方法，其特征在于，在步骤b中，在空气对流条件较好时，在宽广水面及陆地上空，飞行低于300mde八个高度，每个高度不少于8分钟。

4.根据权利要求1所述的无上测晶体大气氡修正参数求取方法，其特征在于，在步骤c中，将装载航空放射性测量设备的飞机推到本底模型上利用铯源(¹³⁷Cs)和钍源(²⁰⁸Tl)测试每条晶体、每箱晶体和整个探测器系统的分辨率和峰位；然后按本底、钾、钍、铀、混合模型的顺序测试一次，再按混合、铀、钍、钾、本底模型的顺序反向测试一次；要求连续完成且放射性测量设备探测器的中心应位于模型正上方，每个模型每次测试时间为5min；最后在本底模型上再次对每条晶体、每箱晶体和整个探测器系统的分辨率和峰位进行检查。

5.根据权利要求1所述的无上测晶体大气氡修正参数求取方法，其特征在于，在步骤g中，利用步骤b采集的放射性谱数据，求取各高度陆地核素谱数据均值以及各高度水面核素谱数据均值，各高度陆地核素谱数据减去对应高度核素谱数据均值均值即为各高度的纯陆地核素谱数据。

6.根据权利要求1所述的无上测晶体大气氡修正参数求取方法，其特征在于，在步骤h中，求取纯陆地核素谱数据各道高度衰减系数μ；μ值用下式求得：

h₁、h₂——标准大气压、标准温度下的飞行高度，m；

R^h1——实际飞行高度h1时纯陆地核素谱数据各道的计数率，计数每秒；

R^h2——实际飞行高度h2时纯陆地核素谱数据各道的计数率，计数每秒；

μ——高度衰减系数，m^-1。

7.根据权利要求1所述的无上测晶体大气氡修正参数求取方法，其特征在于，在步骤i中，利用下式求取钾、铀、钍各核素1m高度单位谱：

A＝MX^-1

X为模型含量矩阵，下标含义：下标第一个字母表示模型，第二个字母表示核素；

M为计数率矩阵，下标含义：下标字母表示模型，数字表示道数；

A为单位谱矩阵，下标含义：下标字母表示核素，数字表示道数。

8.根据权利要求1所述的无上测晶体大气氡修正参数求取方法，其特征在于，在步骤j中，利用下式求取各核素不同高度单位谱：

R^h＝R¹·e^μ(h-1)

h——标准大气压、标准温度下的高度，m；

R^h——高度h时各核素单位谱各道的计数率，计数每秒；

R¹——高度1m时各核素单位谱各道的计数率，计数每秒；

μ——高度衰减系数，m^-1。