本发明属于铀矿勘查中放射性地球物理勘查领域，具体涉及一种车载伽玛全谱测量系统精细标定方法。
背景技术：
车载伽玛能谱测量是一种比航空伽玛能谱经济比地面伽玛能谱效率高的放射性地球物理勘查方法，适合于地势较平缓的草原、戈壁地区的勘查工作。
车载伽玛能谱测量系统标定是该方法获取核素含量的前提，车载伽玛能谱测量系统探测器离地面距离大约1m，理论上需要建立半无限空间的标准模型，实际上，这是无法实现的。设计模型需要满足饱和度达到95％或以上，按这个要求，设计模型半径要达到15m或以上，而建造模型对材料核素含量及密度均匀性等要求高，建造半径15m或以上的模型难度大，且耗资高，前人根据对称性和场源互换原理设计并建立了边长为7m的正六边形体模型，以模型中心及以六个边为对称轴的六个中心点对称点的伽玛强度代替正六边形主体模型及以正六边形的六个边为边的外围六个模型在中心形成的伽玛强度，其等效半径达到16.85m，满足饱和度不小于95％要求。过去，车载伽玛能谱系统标定基本只在设置铀、钍、钾、本底和混合的正六边形柱体模型上的中心点标定，这样模型的等效半径为6.37m，即模型饱和度不满足要求，也没考虑汽车放置位置的影响的校正，这样产生误差较大，影响测量准确性，对于精细全谱测量不利，为获得精细伽玛全谱数据，在标定时需要考虑模型饱和度及车体位置影响，进行必要的修正，达到精细标定的目的，目前还没有相应的车载伽玛全谱测量系统的精细标定方法。
技术实现要素：
本发明要解决的技术问题是提供一种车载伽玛全谱测量系统精细标定方法，提高伽玛全谱测量准确性。
为解决上述技术问题，本发明一种车载伽玛全谱精细测量系统标定方法，该方法包括以下步骤：
步骤一、选定车载伽玛全谱测量系统；
步骤二、设定车载伽玛能谱测量系统的测量周期和参数及核素能量窗口，选取标准的正六边形体钾模型、铀模型、钍模型、本底模型、混合模型；
步骤三、测量，仪器探头晶体中心对准模型的中心位置o，汽车垂直于模型指定的正六边形的一边，并将车载头正对该边，启动仪器测量t分钟停止测量，保存该测量数据，然后将车行驶至以该指定边为对称轴模型中心点的对称点a1位置停车，车身方向不变，仍与指定正六边形边垂直，测量t分钟停止测量，保存该测量数据，然后车身以该对称点a1为中心，分别转60°、120°、180°、240°、300°停车，依次分别测量t分钟停止测量，仪器分别保存该测量数据；
步骤四、按步骤三分别在铀、钍、钾、本底、混合模型上进行测量，仪器分别自动保存相应的测量数据；
步骤五、计算铀、钍、钾模型a1位置伽玛强度不等对测量产生影响的修正系数，将每个模型a1位置车身分别为0°和180°方向读数平均值与本底模型a1位置车身分别为0°和180°方向读数的平均值比，得修正系数v；
步骤六、计算角响应修正系数，计算铀、钍、钾模型a1位置车身分别为0°、60°、120°、180°、240°和300°方向读数平均值与该模型a1位置车身分别为0°和180°方向读数平均值比，得修正系数w；
步骤七、计算每个模型修正后的净计数率,数值上等于模型中心净计数率加上外围6个模型修正后的净计数率；
步骤八、计算得到灵敏度aij。
所述步骤七中，所述所净计数率计算入下式：
n(o,o)＝n′(o)+∑n′(a)
n′(o)＝n(o)-nb(o)
n′(a)＝w·[n(a)-v·nb(a)]
n(o,o)为模型修正后的净计数率，n(o)、n′(o)为各模型中心点o实测计数率及其净计数率，n(a)、n′(a)为模型外围6个正六边形中心的实测值及其修正值，n(o)，nb(a)为本底模型在o，a点实测值。
所述步骤八中，灵敏度aij计算如下式：
aij表示第j种核素含量对第i道址或窗计数率的灵敏度；cj表示第j种核素的含量；ni表示第i道址或窗口内修正后的净计数率。
本发明的有益技术效果在于：本发明提供的一种用于铀矿勘查的车载伽玛全谱测量系统精细标定方法，通过场地伽玛射线强度差异、车身产生的角响应影响、模型饱和度影响的修正，极大地提高了测量系统参数的准确性，为开展精细伽玛全谱测量和数据处理提供了前提。
具体实施方式
下面结合实施例对本发明作进一步详细说明。
一种车载伽玛全谱精细测量系统标定方法，该方法包括以下步骤：
步骤一、选定车载伽玛全谱测量系统，要求仪器道数大于等于256道，具有天然峰自动稳谱、导航、显示、自动储存功能；
步骤二、设定车载伽玛能谱测量系统的测量周期和参数及核素能量窗口，选取标准的正六边形体钾模型、铀模型、钍模型、本底模型、混合模型；
步骤三、测量，用细绳系住一个小铁锥形体重锤底部中心，将细绳另一端固定在仪器探头晶体中心正对的汽车底盘位置，将汽车开进正六角形模型上，重锤尖端对准模型的中心位置o，汽车垂直于模型指定的正六边形的一边，并将车载头正对该边，启动仪器测量(10～15)分钟停止测量，仪器自动保存该测量数据，然后将车行驶至以该指定边为对称轴模型中心点的对称点a1位置停车，车身方向不变，仍与指定正六边形边垂直，测量(10～15)分钟停止测量，仪器自动保存该测量数据，然后车身以该对称点a1为中心，逆时针分别转60°、120°、180°、240°、300°停车，依次分别测量(10～15)分钟停止测量，仪器分别自动保存该测量数据；
步骤四、按步骤三分别在铀、钍、钾、本底、混合模型上进行测量，仪器分别自动保存相应的测量数据；
步骤五、计算铀、钍、钾模型a1位置伽玛强度不等对测量产生影响的修正系数，将每个模型a1位置车身分别为0°和180°方向读数平均值与本底模型a1位置车身分别为0°和180°方向读数的平均值比，得修正系数v；
步骤六、计算角响应修正系数，计算铀、钍、钾模型a1位置车身分别为0°、60°、120°、180°、240°和300°方向读数平均值与该模型a1位置车身分别为0°和180°方向读数平均值比，得修正系数w；
步骤七、计算每个模型修正后的净计数率,数值上等于模型中心净计数率加上外围6个模型修正后的净计数率，如公式(1)；
n(o,o)＝n′(o)+∑n′(a)(1)
n′(o)＝n(o)-nb(o)(2)
n′(a)＝w·[n(a)-v·nb(a)](3)
公式(1)～公式(3)中，n(o,o)为模型修正后的净计数率，n(o)、n′(o)为各模型中心点o实测计数率及其净计数率，n(a)、n′(a)为模型外围6个正六边形中心的实测值及其修正值，n(o)，nb(a)为本底模型在o，a点实测值。
步骤八、利用公式(4)计算灵敏度aij；
公式(4)中，cj表示第j种核素的含量；ni表示第i道址或窗口内修正后的净计数率。
aij表示第j种核素含量对第i道址或窗计数率的灵敏度
将钾模型、铀模型、钍模型上钾、铀、钍修正后的净计数率及模型钾、铀、钍含量代入公式(4)解方程组即可求得灵敏度aij。
实际测量中，根据车载伽玛能谱测量系统实测钾、铀、钍计数率及灵敏度aij，即可获得钾、铀、钍含量。
本发明的一种车载伽玛全谱测量系统的精细标定方法，具体应用如下
步骤(1)：模型及测量系统选择，采用河北石家庄核工业计量站航空伽玛能谱标准正六角模型及gr-660型256道车载伽玛能谱测量系统；
步骤(2)：标定前将系统标定文件中的钾、铀、钍和总道计数率四个参数的本底(cps)设置为零，钾、铀、钍的灵敏度系数设置为1，矩阵对角线为1，剥离比α、β、γ、a、b、g等其它参数设置为零，并设置钾窗口为：钾窗：(111～126)ch，铀：(135～155)ch，钍：(201～232)ch；
步骤(3)：用细绳系住一个小铁锥形体重锤底部中心，将细绳另一端固定在仪器探头晶体中心正对的汽车底盘位置，将汽车开进钾模型上，重锤尖端对准模型的中心位置o，汽车垂直于模型指定的正六边形的一边，并将车头正对该边，启动仪器测量10分钟停止测量，仪器计数率分别记为nkk(o)、nku(o)、nkt(o)，然后将车行驶至以该指定边为对称轴模型中心点的对称点a1位置停车，车身方向不变，仍与指定正六边形边垂直，测量10分钟停止测量，仪器自动保存该测量数据，仪器计数率分别记为nkk(a1)、nku(a1)、nkt(a1)，然后车身以该对称点a1为中心，逆时针分别转60°、120°、180°、240°、300°停车，依次分别测量10分钟停止测量，仪器分别自动保存该测量数据,仪器计数率分别记为nkk(a2)、nku(a2)、nkt(a2)，nkk(a3)、nku(a3)、nkt(a3)，nkk(a4)、nku(a4)、nkt(a4)，nkk(a5)、nku(a5)、nkt(a5)，nkk(a6)、nku(a6)、nkt(a6)；
步骤(4)：按照步骤(3)分别在铀、钍、混合、本底模型上进行测量，铀模型中心点钾、铀、钍测量值分别记为nuk(o)、nuu(o)、nut(o)，钍模型中心点钾、铀、钍测量值分别记为ntk(o)、ntu(o)、ntt(o)，混合模型中心点钾、铀、钍测量值分别记为nmk(o)、nmu(o)、nmt(o)，本底模型中心点钾、铀、钍测量值分别记为nbk(o)、nbu(o)、nbt(o)，模型中心点的对称点a1逆时针分别转0°、60°、120°、180°、240°、300°位置铀、钍测量值分别记为nuk(a1)、nuu(a1)、nut(a1)，nuk(a2)、nuu(a2)、nut(a2)，nuk(a3)、nuu(a3)、nut(a3)，nuk(a4)、nuu(a4)、nut(a4)，nuk(a5)、nuu(a5)、nut(a5)，nuk(a6)、nuu(a6)、nut(a6)；ntk(a1)、ntu(a1)、ntt(a1)，ntk(a1)、ntu(a2)、ntt(a2)，ntk(a3)、ntu(a3)、ntt(a3)，ntk(a4)、ntu(a4)、ntt(a4)，ntk(a5)、ntu(a5)、ntt(a5)，ntk(a6)、ntu(a6)、ntt(a6)；nmk(a1)、nmu(a1)、nmt(a1)，nmk(a2)、nmu(a2)、nmt(a2)，nmk(a3)、nmu(a3)、nmt(a3)，nmk(a4)、nmu(a4)、nmt(a4)，nmk(a5)、nmu(a5)、nmt(a5)，nmk(a6)、nmu(a6)、nmt(a6)；nbk(a1)、nbu(a1)、nbt(a1)，nbk(a2)、nbu(a2)、nbt(a2)，nbk(a3)、nbu(a3)、nbt(a3)，nbk(a4)、nbu(a4)、nbt(a4)，nbk(a5)、nbu(a5)、nbt(a5)，nbk(a6)、nbu(a6)、nbt(a6)。
步骤(5)：计算钾模型、铀模型、钍模型、本底模型a1位置车身分别为0°和180°方向读数的平均值，再该读数平均值除以本底模型a1位置车身分别为0°和180°方向读数的平均值，得因模型和本底模型外a1点场地钾、铀、钍含量差异对测量产生的影响的修正系数,即钾模型钾、铀、钍修正系数vkk、vku、vkt,铀模型钾、铀、钍修正系数vuk、vuu、vut,钍模型钾、铀、钍修正系数vtk、vtu、vtt；
所述步骤(5)中，钾模型a1场地钾、铀、钍含量差异对测量产生影响修正系数分别为vkk＝[nkk(a1)+nkk(a3)]/[nbk(a1)+nbk(a3)]，vku＝[nku(a1)+nku(a3)]/[nbu(a1)+nbu(a3)]，vkt＝[nkt(a1)+nkt(a3)]/[nbt(a1)+nbt(a3)]，其它类似。
步骤(6)：角度响应修正系数，计算钾、铀、钍模型a1位置车身方向分别为0°、60°、120°、180°、240°和300°时各读数的平均值，再将该平均值除以该模型a1位置车身方向分别为0°和180°时读数平均值，得角度响应修正系数，即钾模型钾、铀、钍角度响应修正系数wkk、wku、wkt，铀模型钾、铀、钍角度响应修正系数wkk、wku、wkt，钍模型钾、铀、钍角度响应修正系数wkk、wku、wkt；
所述步骤(5)中，钾模型a1场地钾、铀、钍角度响应修正系数分别为vkk＝[nkk(a1)+nkk(a2)+nkk(a3)+nkk(a4)+nkk(a5)+nkk(a6)]/[nkk(a1)+nkk(a3)]，vku＝[nku(a1)+nku(a2)+nku(a3)+nku(a4)+nku(a5)+nku(a6)]/[nku(a1)+nku(a3)]，vkt＝[nkt(a1)+nkt(a2)+nkt(a3)+nkt(a4)+nkt(a5)+nkt(a6)]/[nkt(a1)+nkt(a3)]，其它类似。
步骤(7)：计算钾、铀、钍模型修正后的钾、铀、钍净计数率，模型修正后的钾、铀、钍净计数率为模型o点实测值减去本底模型o点的实测值加上a1点车身方向分别为0°、60°、120°、180°、240°和300°实测值修正后的六个值和。
所述步骤(7)中，对钾模型，模型o点钾净计数率表示为
n′kk(o)＝nkk(o)-nbk(o)(5)
模型o点铀净计数率表示为
n′ku(o)＝nku(o)-nku(bg,o)(6)
模型o点钍净计数率表示为
n′kt(o)＝nkt(o)-nbt(o)(7)
a1点车身方向0°实测钾计数率修正后值表示为
n′kk(a1)＝wkk[nkk(a1)-vkknbk(a1)](8)
a1点车身方向0°实测铀计数率修正后值表示为
n′ku(a1)＝wku[nku(a1)-vkunbu(a1)](9)
a1点车身方向0°实测钍计数率修正后值表示为
n′kt(a1)＝wkt[nkt(a1)-vktnbt(a1)](10)
类似，可得到a1点车身方向60°、120°、180°、240°和300°的实测钾、铀、钍计数率修正后值n′kk(ak)、n′ku(ak)、n′kt(ak)，(k＝2,3,4,5,6)。
从而得到钾模型钾修正后的净计数率为
钾模型铀修正后的净计数率为
钾模型钍修正后的净计数率为
对铀模型，模型o点钾净计数率表示为
n′uk(o)＝nuk(o)-nbk(o)(14)
模型o点铀净计数率表示为
n′uu(o)＝nuu(o)-nbu(o)(15)
模型o点钍净计数率表示为
n′tt(o)＝ntt(o)-nbt(o)(16)
a1点车身方向0°实测钾计数率修正后值表示为
n′uk(a1)＝wuk[nuk(a1)-vuknbk(a1)](17)
a1点车身方向0°实测铀计数率修正后值表示为
n′uu(a1)＝wuu[nuu(a1)-vuunbu(a1)](18)
a1点车身方向0°实测钍计数率修正后值表示为
n′ut(a1)＝wut[nut(a1)-vutnbt(a1)](19)
类似，可得到a1点车身方向60°、120°、180°、240°和300°的实测钾、铀、钍计数率修正后值n′uk(ak)、n′uu(ak)、n′ut(ak)，(k＝2,3,4,5,6)。
从而得到铀模型钾修正后的净计数率为
铀模型铀修正后的净计数率为
铀模型钍修正后的净计数率为
对钍模型，模型o点钾净计数率表示为
n′tk(o)＝ntk(o)-nbk(o)(23)
模型o点铀净计数率表示为
n′tu(o)＝ntu(o)-nbu(o)(24)
模型o点钍净计数率表示为
n′tt(o)＝ntt(o)-nbt(o)(25)
a1点车身方向0°实测钾计数率修正后值表示为
n′kk(a1)＝wkk[nkk(a1)-vkknbk(a1)](26)
a1点车身方向0°实测铀计数率修正后值表示为
n′tu(a1)＝wtu[ntu(a1)-vtunbu(a1)](27)
a1点车身方向0°实测钍计数率修正后值表示为
n′tt(a1)＝wtt[ntt(a1)-vttnbt(a1)](28)
类似，可得到a1点车身方向60°、120°、180°、240°和300°的实测钾、铀、钍计数率修正后值n′tk(ak)、n′tu(ak)、n′tt(ak)，(k＝2,3,4,5,6)。
从而得到钍模型钾修正后的净计数率为
钍模型铀修正后的净计数率为
钍模型钍修正后的净计数率为
步骤(8)：利用公式(4)计算灵敏度aij。
将钾模型、铀模型、钍模型上测得的钾、铀、钍修正后的净计数率及模型钾、铀、钍含量代入公式(4)解方程组得方程组(32)～方程组(34)：
方程组(32)～方程组(34)中，akk、aku、akt分别表示单位钾、铀、钍含量对钾窗计数率的灵敏度，auk、auu、aut分别、单位钾、铀、钍含量对铀窗计数率的灵敏度，atk、atu、att分别表示单位钾、铀、钍含量对钍窗计数率的灵敏度，ckk、cku、ckt分别表示钾模型钾、铀、钍含量，cuk、cuu、cut分别表示钾模型钾、铀、钍含量，ctk、ctu、ctt分别表示钾模型钾、铀、钍含量
联立方程组(32)～方程组(34)即可就得六个灵敏度akk、aku、akt、akk、akk、akk、akk、akk、akk
实际测量时，通过六个灵敏度akk、aku、akt、akk、akk、akk、akk、akk、akk，实测净计数率，利用公式(4)即可获得钾、铀、钍含量。