CN104121016A - 一种用于铀矿勘探的伽玛能谱测井仪校正方法 - Google Patents

Abstract

本发明属于核探测仪器校正技术领域，具体涉及一种用于铀矿勘探的伽玛能谱测井仪校正方法。本发明的方法包括以下步骤：设置伽玛能谱测井仪能窗参数并选取基准模型；计算并修正死时间；计算灵敏度系数；校验混合校验模型。本发明解决了现有技术难以对伽玛能谱测井仪进行有效校正的技术问题；采用本发明的方法，能够对伽玛能谱测井仪进行有效校正，校正后的伽玛能谱测井仪可以准确测定矿床钻孔地层中的铀、钍含量。

Description

一种用于铀矿勘探的伽玛能谱测井仪校正方法

技术领域

本发明属于核探测仪器校正技术领域，具体涉及一种用于铀矿勘探的伽玛能谱测井仪校正方法。

背景技术

目前，在铀矿勘查领域中，普遍采用伽玛总量测井来计算钻孔地层中的铀含量，在铀钍混合型矿床上，需要采取部分矿段岩心进行铀、钍含量化验分析，进而对伽玛总量测井获得的铀含量进行校正。

为准确测定铀钍混合型矿床上的铀、钍含量，需采用伽玛能谱测井仪。如何有效校正伽玛能谱测井仪，是本领域技术人员亟需解决的关键问题

发明内容

本发明需解决的技术问题为：现有技术难以对伽玛能谱测井仪进行有效校正。

本发明的技术方案如下所述：

一种用于铀矿勘探的伽玛能谱测井仪校正方法，包括以下步骤：首先采用伽玛能谱测井仪在铀基准模型、钍基准模型和混合校验模型上进行伽玛能谱测量，计算铀基准模型和钍基准模型上各自的铀窗和钍窗死时间，然后根据经死时间校正后的铀窗和钍窗的计数率和铀、钍基准模型的铀、钍含量，计算伽玛能谱测井仪的灵敏度系数；之后根据铀基准模型和钍基准模型上各自的铀窗和钍窗死时间，计算混合校验模型的铀窗和钍窗死时间，然后计算混合校验模型上经死时间校正后的铀窗和钍窗计数率；最后采用上述计算出的灵敏度系数和混合校验模型的铀窗、钍窗计数率，计算混合校验模型的铀、

本发明的方法具体包括以下步骤：

步骤1.设置伽玛能谱测井仪能窗参数并选取基准模型；

步骤2.计算并修正死时间；

步骤2.1.采用双源法和线性相关系数法计算铀基准模型上铀窗死时间T_UU和钍窗死时间T_UTh；

步骤2.2.采用与步骤2.1相同的方法计算钍基准模型上最终铀窗死时间T_ThU和最终钍窗死时间T_ThTh；

步骤2.3.计算混合校验模型的铀窗死时间T_MU和钍窗死时间T_MTh；

步骤2.4.计算混合校验模型上经死时间修正后铀窗计数率和钍窗计数率；

步骤3.计算灵敏度系数；

步骤4.校验混合校验模型。

步骤1具体包括以下步骤：

对待校正的伽玛能谱测井仪，采用“宽两窗”，即U窗和Th窗；选取a个铀基准模型、b个钍基准模型和c个混合校验模型；采用核工业放射性勘查计量站的放射性测井基准模型作为基准模型。

步骤1中，a的取值范围为3～6；b的取值范围为2～5；c的取值范围为2～6。

作为步骤1的优选方案，a=4，b=3，c=2；

基准模型参数如下表所示：

。

步骤2具体包括以下步骤：

步骤2.1具体包括以下步骤：

步骤2.1.1.采用双源法计算铀基准模型上初始铀窗死时间τ_UU和初始钍窗死时间τ_UTh

利用下式计算铀基准模型上初始铀窗死时间τ_UU和初始钍窗死时间τ_UTh：

$\left\{\begin{array}{c}{\mathrm{\τ}}_{\mathrm{UU}}=\frac{C_{U1}\·N_{2U}-C_{U2}\·N_{1U}}{C_{U1}\·N_{2U}\·N_{1\mathrm{TC}}-C_{U2}\·N_{1U}\·N_{2\mathrm{TC}}}\\ {\mathrm{\τ}}_{\mathrm{UTh}}=\frac{C_{U1}\·N_{2\mathrm{Th}}-C_{U2}\·N_{1\mathrm{Th}}}{C_{U1}\·N_{2\mathrm{Th}}\·N_{1\mathrm{TC}}-C_{U2}\·N_{1\mathrm{Th}}\·N_{2\mathrm{TC}}}\end{array}\right.;$

式中，

τ_UU表示铀模型上计算的初始铀窗死时间，单位：s；

τ_UTh表示铀模型上计算的初始钍窗死时间，单位：s；

C_U1表示1#铀基准模型的铀含量，即铀的质量百分比，单位：%；

C_U2表示2#铀基准模型的铀含量，即铀的质量百分比，单位：%；

N_1U表示1#铀基准模型上采集的谱线U窗的计数率，单位：s^-1；

N_1Th表示1#铀基准模型上采集的谱线Th窗的计数率，单位：s^-1；

N_1TC表示1#铀基准模型上采集的谱线总道的计数率，单位：s^-1；

N_2U表示2#铀基准模型上采集的谱线U窗的总计数率，单位：s^-1；

N_2Th表示2#铀基准模型上采集的谱线Th窗的计数率，单位：s^-1；

N_2TC表示2#铀基准模型上采集的谱线总道的计数率，单位：s^-1；

所述1#铀基准模型、2#铀基准模型为步骤1中任意两个不相同的铀基准模型，即有种组合方式，进而得到组初始铀窗死时间和组初始钍窗死时间

步骤2.1.2.采用线性相关系数法计算铀基准模型上最终铀窗死时间T_UU和最终钍窗死时间T_UTh

步骤2.1.2.1.计算铀基准模型上最终铀窗死时间T_UU

设

其中，m₁优选为10的整数倍；

采用下式计算铀基准模型上最终铀窗死时间T_UU：

$\left\{\begin{array}{c}{N}_{\mathrm{UUi}}^c\left(n_1\right)=\frac{N_{\mathrm{UUi}}}{1-N_{\mathrm{UTCi}}\·T_{d1s}\left(n_1\right)},(i=\mathrm{1,2}...a)\\ \stackrel{\‾}{N_{\mathrm{UU}}^c\left(n_1\right)}=\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{N}_{\mathrm{UUi}}^c\left(n_1\right)/a\\ \stackrel{\‾}{C_U}=\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{C}_{\mathrm{Ui}}/a\\ R_{\mathrm{UU}}\left(n_1\right)=\frac{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}(N_{\mathrm{UUi}}^c\left(n_1\right)-\stackrel{\‾}{N_{\mathrm{UU}}^c\left(n_1\right)})(C_{\mathrm{Ui}}-\stackrel{\‾}{C_U})}{\sqrt{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{(N_{\mathrm{UUi}}^c\left(n_1\right)-\stackrel{\‾}{N_{\mathrm{UU}}^c\left(n_1\right)})}^2}\sqrt{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{(C_{\mathrm{Ui}}-\stackrel{\‾}{C_U})}^2}}\\ R_{\mathrm{UU}}\left(l_1\right)=\max \left\{R_{\mathrm{UU}}\left(n_1\right)\right\}\\ T_{\mathrm{UU}}=T_{d1s}\left(l_1\right)=T_{d1s\min }+l_1\·\mathrm{\Δ}{t}_1\end{array}\right.$

式中，

表示第i个铀基准模型经n₁次死时间校正后的铀窗计数率，单位：s-1；

N_UUi表示第i个铀基准模型的铀窗原始计数率，单位：s^-1；

N_UTCi表示第i个铀基准模型的总道原始计数率，单位：s^-1；

C_Ui表示第i个铀基准模型的铀含量，即铀质量百分比，单位：%；

R_UU(n₁)表示所有铀基准模型经n₁次死时间校正后的铀窗计数率与铀含量之间的相关系数；

相关系数R_UU(n₁)的最大值对应的死时间即为铀基准模型上最终铀窗死时间T_UU；

步骤2.1.2.2.计算铀基准模型上钍窗死时间T_UTh

采用与步骤2.1.2.1相同的线性相关系数法方法计算钍窗死时间T_UTh；

步骤2.2.采用与步骤2.1相同的方法计算钍基准模型上最终铀窗死时间T_ThU和最终钍窗死时间T_ThTh；

步骤2.3具体包括以下步骤：

采用线性插值法计算混合校验模型的铀窗死时间T_MU和钍窗死时间T_MTh：

$\left\{\begin{array}{c}{P}_U=\frac{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{p}_{\mathrm{Ui}}}{a}\\ P_{\mathrm{Th}}=\frac{\underset{j=1}{\overset{b}{\mathrm{\Σ}}}{p}_{\mathrm{THj}}}{b}\\ T_{\mathrm{MU}}=\frac{T_{\mathrm{UU}}-T_{\mathrm{ThU}}}{P_U-P_{\mathrm{Th}}}\·(P_M-P_{\mathrm{Th}})+T_{\mathrm{ThU}}\\ T_{\mathrm{MTh}}=\frac{T_{\mathrm{UTh}}-T_{\mathrm{ThTh}}}{P_U-P_{\mathrm{Th}}}\·(P_M-P_{\mathrm{Th}})+T_{\mathrm{ThTh}}\end{array}\right.;$

式中，

p_Ui表示第i个铀基准模型上铀窗与钍窗的原始计数率比值；

p_Thj表示第j个钍基准模型上铀窗与钍窗的原始计数率比值；

P_M表示混合校验模型上铀窗与钍窗的原始计数率比值；

T_MU表示混合校验模型铀窗死时间，单位：s；

T_MTh表示混合校验模型钍窗死时间，单位：s；

步骤2.4具体包括以下步骤：

采用下式计算和

$\left\{\begin{array}{c}{N}_{\mathrm{MU}}^c=\frac{N_{\mathrm{MU}}}{1-N_{\mathrm{MTC}}\·T_{\mathrm{MU}}}\\ N_{\mathrm{MTh}}^c=\frac{N_{\mathrm{MTh}}}{1-N_{\mathrm{MTC}}\·T_{\mathrm{MTh}}}\end{array}\right.;$

式中，

表示经死时间校正后铀窗计数率，单位s^-1；

表示经死时间校正后钍窗计数率，单位s^-1；

N_MTC表示混合校验模型上总道计数率，单位s^-1。

步骤3具体包括以下步骤：

对基准模型依据下式建立方程组，采用最小二乘法解方程组得到四个灵敏度系数值S_Uu、S_Thu、S_Uth和S_Thth：

$\left\{\begin{array}{c}{S}_{\mathrm{Uu}}\·C_{\mathrm{Mu}}+S_{\mathrm{Uth}}\·C_{\mathrm{Mth}}=N_{\mathrm{MU}}^c\\ S_{\mathrm{Thu}}\·C_{\mathrm{Mu}}+S_{\mathrm{Thth}}\·C_{\mathrm{Mth}}=N_{\mathrm{MTh}}^c\end{array}\right.;$

式中，

S_Uu表示单位含量的铀对铀窗产生的伽玛计数率，即铀对铀窗的灵敏度系数，单位：s^-1·1%eU；

S_Thu表示单位含量的铀对钍窗产生的伽玛计数率，即铀对钍窗的灵敏度系数，单位：s^-1·1%eU；

S_Uth表示单位含量的钍对铀窗产生的伽玛计数率，即钍对铀窗的灵敏度系数，单位：s^-1·1%Th；

S_Thth表示单位含量的钍对钍窗产生的伽玛计数率，即钍对钍窗的灵敏度系数，单位：s^-1·1%Th；

C_Mu表示混合校验模型的铀含量，单位：%；

C_Mth表示混合校验模型的钍含量，单位：%；

表示经死时间校正后铀窗计数率，单位s^-1；

表示经死时间校正后钍窗计数率，单位s^-1。

可以采用下表所示模型作为计算灵敏度系数的基准模型：

。

步骤4具体包括以下步骤：

通过下式计算混合校验模型的铀、钍含量：

$\left\{\begin{array}{c}{C}_u^c=\frac{N_U^c\·S_{\mathrm{Thth}}-N_{\mathrm{Th}}^c\·S_{\mathrm{Uth}}}{S_{\mathrm{Uu}}\·S_{\mathrm{Thth}}-S_{\mathrm{Uth}}\·S_{\mathrm{Thu}}}\\ C_{\mathrm{th}}^c=\frac{N_{\mathrm{Th}}^c\·S_{\mathrm{Uu}}-N_U^c\·S_{\mathrm{Thu}}}{S_{\mathrm{Uu}}\·S_{\mathrm{Thth}}-S_{\mathrm{Thu}}\·S_{\mathrm{Uth}}}\end{array}\right.;$

式中，

表示计算获得的混合校验模型的铀含量，单位：%；

表示计算获得的混合校验模型的钍含量，单位：%；

通过下式计算铀、钍含量误差：

$\left\{\begin{array}{c}{\mathrm{\η}}_u=\frac{C_u^c-C_u}{C_u}\×100\%\\ {\mathrm{\η}}_{\mathrm{th}}=\frac{C_{\mathrm{th}}^c-C_{\mathrm{th}}}{C_{\mathrm{th}}}\×100\%\end{array}\right.$

式中，

C_u—混合校验模型的标称铀含量，单位：%；

C_th—混合校验模型的标称钍含量，单位：%；

η_u—铀含量相对误差，单位：%；

η_th—钍含量相对误差，单位：%。

本发明的有益效果为：

本发明的一种用于铀矿勘探的伽玛能谱测井仪校正方法，能够对伽玛能谱测井仪进行有效校正，校正后的伽玛能谱测井仪可以准确测定矿床钻孔地层中的铀、钍含量。

具体实施方式

下面结对本发明的一种用于铀矿勘探的伽玛能谱测井仪校正方法进行详细说明。

本发明的一种用于铀矿勘探的伽玛能谱测井仪校正方法，首先采用伽玛能谱测井仪在铀基准模型、钍基准模型和混合校验模型上进行伽玛能谱测量，计算铀基准模型和钍基准模型上各自的铀窗和钍窗死时间，然后根据经死时间校正后的铀窗和钍窗的计数率和铀、钍基准模型的铀、钍含量，计算伽玛能谱测井仪的灵敏度系数；之后根据铀基准模型和钍基准模型上各自的铀窗和钍窗死时间，计算混合校验模型的铀窗和钍窗死时间，然后计算混合校验模型上经死时间校正后的铀窗和钍窗计数率；最后采用上述计算出的灵敏度系数和、混合校验模型的铀窗、钍窗计数率，计算混合校验模型的铀、钍含量，与其标称含量进行比较，计算出铀、钍含量相对误差即可。一般铀、钍含量的相对误差均小于10%即可视为校正通过。

具体而言，本发明的方法包括以下步骤：

步骤1.设置伽玛能谱测井仪能窗参数并选取基准模型

对待校正的伽玛能谱测井仪，采用“宽两窗”，即U窗和Th窗。本实施例中，U窗范围为1.05～2.00MeV，Th窗范围为2.00～2.90MeV。

选取a个铀基准模型、b个钍基准模型和c个混合校验模型，a的取值范围为3～6；b的取值范围为2～5；c的取值范围为2～6。本实施例中，采用核工业放射性勘查计量站的放射性测井基准模型作为基准模型，a=4，b=3，c=2，其相关参数如表1所示。

表1

步骤2.计算并修正死时间

步骤2.1.计算铀基准模型上铀窗死时间T_UU和钍窗死时间T_UTh

步骤2.1.1.采用双源法计算铀基准模型上初始铀窗死时间τ_UU和初始钍窗死时间τ_UTh

利用下式计算铀基准模型上初始铀窗死时间τ_UU和初始钍窗死时间τ_UTh：

$\left\{\begin{array}{c}{\mathrm{\τ}}_{\mathrm{UU}}=\frac{C_{U1}\·N_{2U}-C_{U2}\·N_{1U}}{C_{U1}\·N_{2U}\·N_{1\mathrm{TC}}-C_{U2}\·N_{1U}\·N_{2\mathrm{TC}}}\\ {\mathrm{\τ}}_{\mathrm{UTh}}=\frac{C_{U1}\·N_{2\mathrm{Th}}-C_{U2}\·N_{1\mathrm{Th}}}{C_{U1}\·N_{2\mathrm{Th}}\·N_{1\mathrm{TC}}-C_{U2}\·N_{1\mathrm{Th}}\·N_{2\mathrm{TC}}}\end{array}\right..$

式中，

τ_UU表示铀模型上计算的初始铀窗死时间，单位：s；

τ_UTh表示铀模型上计算的初始钍窗死时间，单位：s；

C_U1表示1#铀基准模型的铀含量，即铀的质量百分比，单位：%；

C_U2表示2#铀基准模型的铀含量，即铀的质量百分比，单位：%；

N_1U表示1#铀基准模型上采集的谱线U窗的计数率，单位：s^-1；

N_1Th表示1#铀基准模型上采集的谱线Th窗的计数率，单位：s^-1；

N_1TC表示1#铀基准模型上采集的谱线总道的计数率，单位：s^-1；

N_2U表示2#铀基准模型上采集的谱线U窗的总计数率，单位：s^-1；

N_2Th表示2#铀基准模型上采集的谱线Th窗的计数率，单位：s^-1；

N_2TC表示2#铀基准模型上采集的谱线总道的计数率，单位：s^-1。

所述1#铀基准模型、2#铀基准模型为步骤1中任意两个不相同的铀基准模型，即有种组合方式，进而得到组初始铀窗死时间和组初始钍窗死时间本实施例中，共有6种组合方式，进而得到6组初始铀窗死时间和6组初始钍窗死时间

步骤2.1.2.采用线性相关系数法计算铀基准模型上最终铀窗死时间T_UU和最终钍窗死时间T_UTh

步骤2.1.2.1.计算铀基准模型上最终铀窗死时间T_UU

设其中，优选m₁为10的整数倍。

采用下式计算铀基准模型上最终铀窗死时间T_UU：

$\left\{\begin{array}{c}{N}_{\mathrm{UUi}}^c\left(n_1\right)=\frac{N_{\mathrm{UUi}}}{1-N_{\mathrm{UTCi}}\·T_{d1s}\left(n_1\right)},(i=\mathrm{1,2}...a)\\ \stackrel{\‾}{N_{\mathrm{UU}}^c\left(n_1\right)}=\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{N}_{\mathrm{UUi}}^c\left(n_1\right)/a\\ \stackrel{\‾}{C_U}=\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{C}_{\mathrm{Ui}}/a\\ R_{\mathrm{UU}}\left(n_1\right)=\frac{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}(N_{\mathrm{UUi}}^c\left(n_1\right)-\stackrel{\‾}{N_{\mathrm{UU}}^c\left(n_1\right)})(C_{\mathrm{Ui}}-\stackrel{\‾}{C_U})}{\sqrt{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{(N_{\mathrm{UUi}}^c\left(n_1\right)-\stackrel{\‾}{N_{\mathrm{UU}}^c\left(n_1\right)})}^2}\sqrt{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{(C_{\mathrm{Ui}}-\stackrel{\‾}{C_U})}^2}}\\ R_{\mathrm{UU}}\left(l_1\right)=\max \left\{R_{\mathrm{UU}}\left(n_1\right)\right\}\\ T_{\mathrm{UU}}=T_{d1s}\left(l_1\right)=T_{d1s\min }+l_1\·\mathrm{\Δ}{t}_1\end{array}\right.$

式中，

表示第i个铀基准模型经n1次死时间校正后的铀窗计数率，单位：s^-1；

N_UUi表示第i个铀基准模型的铀窗原始计数率，单位：s^-1；

N_UTCi表示第i个铀基准模型的总道原始计数率，单位：s^-1；

C_Ui表示第i个铀基准模型的铀含量，即铀质量百分比，单位：%；

R_UU(n₁)表示所有铀基准模型经n₁次死时间校正后的铀窗计数率与铀含量之间的相关系数。

相关系数R_UU(n₁)的最大值对应的死时间即为铀基准模型上最终铀窗死时间T_UU。

步骤2.1.2.2.计算铀基准模型上钍窗死时间T_UTh

采用与步骤2.1.2.1相同的线性相关系数法方法计算钍窗死时间T_UTh。

设其中，优选m₂为10的整数倍。

采用下式计算铀基准模型上最终钍窗死时间T_UTh：

$\left\{\begin{array}{c}{N}_{\mathrm{UThi}}^c\left(n_2\right)=\frac{N_{\mathrm{UThi}}}{1-N_{\mathrm{UTCi}}\·T_{d2s}\left(n_2\right)},(i=\mathrm{1,2}...a)\\ \stackrel{\‾}{N_{\mathrm{UTh}}^c\left(n_2\right)}=\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{N}_{\mathrm{UThi}}^c\left(n_2\right)/a\\ \stackrel{\‾}{C_U}=\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{C}_{\mathrm{Ui}}/a\\ R_{\mathrm{UTh}}\left(n_2\right)=\frac{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}(N_{\mathrm{UThi}}^c\left(n_2\right)-\stackrel{\‾}{N_{\mathrm{UTh}}^c\left(n_2\right)})(C_{\mathrm{Ui}}-\stackrel{\‾}{C_U})}{\sqrt{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{(N_{\mathrm{UThi}}^c\left(n_2\right)-\stackrel{\‾}{N_{\mathrm{UTh}}^c\left(n_2\right)})}^2}\sqrt{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{(C_{\mathrm{Ui}}-\stackrel{\‾}{C_U})}^2}}\\ R_{\mathrm{UTh}}\left(l_2\right)=\max \left\{R_{\mathrm{UTh}}\left(n_2\right)\right\}\\ T_{\mathrm{UU}}=T_{d2s}\left(l_2\right)=T_{d2s\min }+l_2\·\mathrm{\Δ}{t}_2\end{array}\right.$

式中，

表示第i个铀基准模型经n2次死时间校正后的钍窗计数率，单位：s^-1；

R_UTh(n₂)表示所有铀基准模型经n₂次死时间校正后的钍窗计数率与铀含量之间的相关系数。

相关系数R_UTh(n₂)的最大值对应的死时间即为铀基准模型上最终钍窗死时间T_UTh。

步骤2.2.计算钍基准模型上最终铀窗死时间T_ThU和最终钍窗死时间T_ThTh

采用与步骤2.1相同的双源法计算钍基准模型上初始铀窗死时间τ_ThU和初始钍窗死时间τ_ThTh；采用与步骤2.1相同的线性相关系数法计算法计算钍基准模型上最终铀窗死时间T_ThU和最终钍窗死时间T_ThTh。

步骤2.2.1.计算钍基准模型上初始铀窗死时间τ_ThU和初始钍窗死时间τ_ThTh

利用下式计算钍基准模型上初始铀窗死时间τ_ThU和初始钍窗死时间τ_ThTh：

$\left\{\begin{array}{c}{\mathrm{\τ}}_{\mathrm{ThU}}=\frac{C_{\mathrm{Th}1}\·{N_{2U}}^{\′}-C_{\mathrm{Th}2}\·{N_{1U}}^{\′}}{C_{\mathrm{Th}1}\·{N_{2U}}^{\′}\·{N_{1\mathrm{TC}}}^{\′}-C_{\mathrm{Th}2}\·{N_{1U}}^{\′}\·{N_{2\mathrm{TC}}}^{\′}}\\ {\mathrm{\τ}}_{\mathrm{ThTh}}=\frac{C_{\mathrm{Th}1}\·{N_{2\mathrm{Th}}}^{\′}-C_{\mathrm{Th}2}\·{N_{1\mathrm{Th}}}^{\′}}{C_{\mathrm{Th}1}\·{N_{2\mathrm{Th}}}^{\′}\·{N_{1\mathrm{TC}}}^{\′}-C_{\mathrm{Th}2}\·{N_{1\mathrm{Th}}}^{\′}\·{N_{2\mathrm{Th}}}^{\′}}\end{array}\right..$

τ_ThU—钍模型上计算的初始铀窗死时间，单位：s；

τ_ThTh—钍模型上计算的初始钍窗死时间，单位：s；

C_Th1—1#钍基准模型的钍含量，即钍的质量百分比，单位：%；

C_Th2—2#钍基准模型的钍含量，即钍的质量百分比，单位：%；

N_1U'—1#钍基准模型上U窗的计数率，单位：s^-1；

N_1Th'—1#钍基准模型上采集的谱线Th窗的计数率，单位：s^-1；

N_1TC'—1#钍基准模型上采集的谱线总道的计数率，单位：s^-1；

N_2U'—2#钍基准模型上采集的谱线U窗的总计数率，单位：s^-1；

N_2Th'—2#钍基准模型上采集的谱线Th窗的计数率，单位：s^-1；

N_2TC'—2#钍基准模型上采集的谱线总道的计数率，单位：s^-1。

所述1#钍基准模型、2#钍基准模型为步骤1中任意两个不相同的钍基准模型，即有种组合方式，进而得到组初始铀窗死时间和组初始钍窗死时间本实施例中，共有3种组合方式，进而得到3组初始铀窗死时间τ_ThU1～τ_ThU3和3组初始钍窗死时间τ_ThTh1～τ_ThTh3。

步骤2.2.2.计算钍基准模型上最终铀窗死时间T_ThU和最终钍窗死时间T_ThTh

步骤2.2.2.1.计算钍基准模型上最终铀窗死时间T_ThU

设其中，优选m₃为10的整数倍。

采用下式计算钍基准模型上最终铀窗死时间T_ThU：

$\left\{\begin{array}{c}{N}_{\mathrm{ThUj}}^c\left(n_3\right)=\frac{N_{\mathrm{ThUj}}}{1-N_{\mathrm{ThTCj}}\·T_{d3s}\left(n_3\right)},(j=\mathrm{1,2}...b)\\ \stackrel{\‾}{N_{\mathrm{ThU}}^c\left(n_3\right)}=\underset{j=1}{\overset{b}{\mathrm{\Σ}}}{N}_{\mathrm{ThUj}}^c\left(n_3\right)/b\\ \stackrel{\‾}{C_{\mathrm{Th}}}=\underset{j=1}{\overset{b}{\mathrm{\Σ}}}{C}_{\mathrm{Thj}}/b\\ R_{\mathrm{ThU}}\left(n_3\right)=\frac{\underset{j=1}{\overset{b}{\mathrm{\Σ}}}(N_{\mathrm{ThUj}}^c\left(n_3\right)-\stackrel{\‾}{N_{\mathrm{ThU}}^c\left(n_3\right)})(C_{\mathrm{Thj}}-\stackrel{\‾}{C_{\mathrm{Th}}})}{\sqrt{\underset{j=1}{\overset{b}{\mathrm{\Σ}}}{(N_{\mathrm{ThUj}}^c\left(n_3\right)-\stackrel{\‾}{N_{\mathrm{ThU}}^c\left(n_3\right)})}^2}\sqrt{\underset{j=1}{\overset{b}{\mathrm{\Σ}}}{(C_{\mathrm{Thj}}-\stackrel{\‾}{C_{\mathrm{Th}}})}^2}}\\ R_{\mathrm{ThU}}\left(l_3\right)=\max \left\{R_{\mathrm{ThU}}\left(n_3\right)\right\}\\ T_{\mathrm{ThU}}=T_{d3s}\left(l_3\right)=T_{d3s\min }+l_3\·\mathrm{\Δ}{t}_3\end{array}\right.$

式中，

表示第j个钍基准模型经n3次死时间校正后的铀窗计数率，单位：s^-1；

N_ThUj表示第j个钍基准模型的铀窗原始计数率，单位：s^-1；

N_ThTCj表示第j个钍基准模型的总道原始计数率，单位：s^-1；

C_Thj表示第j个钍基准模型的钍含量，即铀质量百分比，单位：%；

R_ThU(n₃)表示所有钍基准模型经n₃次死时间校正后的铀窗计数率与钍含量之间的相关系数。

相关系数R_ThU(n₃)的最大值对应的死时间即为钍基准模型上最终铀窗死时间T_ThU。

步骤2.2.2.2.计算钍基准模型上最终钍窗死时间TThTh

设其中，优选m₄为10的整数倍。

采用下式计算钍基准模型上最终钍窗死时间T_ThTh：

$\left\{\begin{array}{c}{N}_{\mathrm{ThThj}}^c\left(n_4\right)=\frac{N_{\mathrm{ThThj}}}{1-N_{\mathrm{ThTCj}}\·T_{d4s}\left(n_4\right)},(j=\mathrm{1,2}...b)\\ \stackrel{\‾}{N_{\mathrm{ThTh}}^c\left(n_4\right)}=\underset{j=1}{\overset{b}{\mathrm{\Σ}}}{N}_{\mathrm{ThThj}}^c\left(n_4\right)/b\\ \stackrel{\‾}{C_{\mathrm{Th}}}=\underset{j=1}{\overset{b}{\mathrm{\Σ}}}{C}_{\mathrm{Thj}}/b\\ R_{\mathrm{ThTh}}\left(n_4\right)=\frac{\underset{i=1}{\overset{b}{\mathrm{\Σ}}}(N_{\mathrm{ThThj}}^c\left(n_4\right)-\stackrel{\‾}{N_{\mathrm{ThTh}}^c\left(n_4\right)})(C_{\mathrm{Thj}}-\stackrel{\‾}{C_{\mathrm{Th}}})}{\sqrt{\underset{j=1}{\overset{b}{\mathrm{\Σ}}}{(N_{\mathrm{ThThj}}^c\left(n_4\right)-\stackrel{\‾}{N_{\mathrm{ThTh}}^c\left(n_4\right)})}^2}\sqrt{\underset{j=1}{\overset{b}{\mathrm{\Σ}}}{(C_{\mathrm{Thj}}-\stackrel{\‾}{C_{\mathrm{Th}}})}^2}}\\ R_{\mathrm{Th}{l}_{4}h}\left(l_4\right)=\max \left\{R_{\mathrm{ThTh}}\left(n_4\right)\right\}\\ T_{\mathrm{ThTh}}=T_{d4s}\left(l_4\right)=T_{d4s\min }+l_4\·\mathrm{\Δ}{t}_4\end{array}\right.$

式中，

表示第j个钍基准模型经n4次死时间校正后的钍窗计数率，单位：s^-1；

N_ThThj表示第j个钍基准模型的钍窗原始计数率，单位：s^-1；

R_ThTh(n₄)表示所有钍基准模型经n₄次死时间校正后的钍窗计数率与钍含量之间的相关系数。

相关系数R_ThTh(n₄)的最大值对应的死时间即为钍基准模型上最终钍窗死时间T_ThTh。

步骤2.3.计算混合校验模型的铀窗死时间T_MU和钍窗死时间T_MTh

根据铀基准模型上铀窗与钍窗的计数率比值、钍基准模型上铀窗与钍窗的计数率比值和混合校验模型上铀窗与钍窗的计数率比值，采用线性插值法计算混合校验模型的铀窗死时间T_MU和钍窗死时间T_MTh：

$\left\{\begin{array}{c}{P}_U=\frac{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{p}_{\mathrm{Ui}}}{a}\\ P_{\mathrm{Th}}=\frac{\underset{j=1}{\overset{b}{\mathrm{\Σ}}}{p}_{\mathrm{Thj}}}{b}\\ T_{\mathrm{MU}}=\frac{T_{\mathrm{UU}}-T_{\mathrm{Thu}}}{P_U-P_{\mathrm{Th}}}\·(P_M-P_{\mathrm{Th}})+T_{\mathrm{ThU}}\\ T_{\mathrm{MTh}}=\frac{T_{\mathrm{UTh}}-T_{\mathrm{ThTh}}}{P_U-P_{\mathrm{Th}}}\·(P_M-P_{\mathrm{Th}})+T_{\mathrm{ThTh}}\end{array}\right..$

式中，

p_Ui表示第i个铀基准模型上铀窗与钍窗的原始计数率比值；

p_Thj表示第j个钍基准模型上铀窗与钍窗的原始计数率比值；

P_M表示混合校验模型上铀窗与钍窗的原始计数率比值；

T_MU表示混合校验模型铀窗死时间，单位：s；

T_MTh表示混合校验模型钍窗死时间，单位：s。

步骤2.4.计算混合校验模型上经死时间修正后铀窗计数率和钍窗计数率

采用下式计算和

$\left\{\begin{array}{c}{N}_{\mathrm{MU}}^c=\frac{N_{\mathrm{MU}}}{1-N_{\mathrm{MTC}}\·T_{\mathrm{MU}}}\\ N_{\mathrm{MTh}}^c=\frac{N_{\mathrm{MTh}}}{1-N_{\mathrm{MTC}}\·T_{\mathrm{MTh}}}\end{array}\right..$

式中，

表示经死时间校正后铀窗计数率，单位s^-1；

表示经死时间校正后钍窗计数率，单位s^-1；

N_MTC表示混合校验模型上总道计数率，单位s^-1。

步骤3.计算灵敏度系数

采用下式计算灵敏度系数：

$\left\{\begin{array}{c}{S}_{\mathrm{Uu}}\·C_{\mathrm{Mu}}+S_{\mathrm{Uth}}\·C_{\mathrm{Mth}}=N_{\mathrm{MU}}^u\\ S_{\mathrm{Thu}}\·C_{\mathrm{Mu}}+S_{\mathrm{Thth}}\·C_{\mathrm{Mth}}=N_{\mathrm{MTh}}^c\end{array}\right..$

式中，

S_Uu表示单位含量的铀对铀窗产生的伽玛计数率，即铀对铀窗的灵敏度系数，单位：s^-1·1%eU；

S_Thu表示单位含量的铀对钍窗产生的伽玛计数率，即铀对钍窗的灵敏度系数，单位：s^-1·1%eU；

S_Uth表示单位含量的钍对铀窗产生的伽玛计数率，即钍对铀窗的灵敏度系数，单位：s^-1·1%Th；

S_Thth表示单位含量的钍对钍窗产生的伽玛计数率，即钍对钍窗的灵敏度系数，单位：s^-1·1%Th；

C_Mu表示混合校验模型的铀含量，单位：%；

C_Mth表示混合校验模型的钍含量，单位：%；

表示经死时间校正后铀窗计数率，单位s^-1；

表示经死时间校正后钍窗计数率，单位s^-1。

考虑到热液型铀矿勘查中铀含量的范围，本实施例中采用表2所示模型作为计算灵敏度系数的基准模型：

表2

根据上表所示基准模型，能够建立下示方程组，采用最小二乘法解此方程组，得到四个灵敏度系数值S_Uu、S_Thu、S_Uth和S_Thth。

$\left\{\begin{array}{c}{S}_{\mathrm{Uu}}\·C_u\left({\mathrm{UF}}_{0.5}\right)+S_{\mathrm{Uth}}\·C_{\mathrm{th}}\left({\mathrm{UF}}_{0.5}\right)=N_U^c\left({\mathrm{UF}}_{0.5}\right)\\ S_{\mathrm{Thu}}\·C_u\left({\mathrm{UF}}_{0.5}\right)+S_{\mathrm{Thth}}\·C_{\mathrm{th}}\left({\mathrm{UF}}_{0.5}\right)=N_{\mathrm{Th}}^c\left({\mathrm{UF}}_{0.5}\right)\\ S_{\mathrm{Uu}}\·C_u\left({\mathrm{UF}}_{0.2}\right)+S_{\mathrm{Uth}}\·C_{\mathrm{th}}\left({\mathrm{UF}}_{0.2}\right)=N_U^c\left({\mathrm{UF}}_{0.2}\right)\\ S_{\mathrm{Thu}}\·C_u\left({\mathrm{UF}}_{0.2}\right)+S_{\mathrm{Thth}}\·C_{\mathrm{th}}\left({\mathrm{UF}}_{0.2}\right)=N_{\mathrm{Th}}^c\left({\mathrm{UF}}_{0.2}\right)\\ S_{\mathrm{Uu}}\·C_u\left({\mathrm{ThF}}_{0.7}\right)+S_{\mathrm{Uth}}\·C_{\mathrm{th}}\left({\mathrm{ThF}}_{0.7}\right)=N_U^c\left({\mathrm{ThF}}_{0.7}\right)\\ S_{\mathrm{Thu}}\·C_u\left({\mathrm{ThF}}_{0.7}\right)+S_{\mathrm{Thth}}\·C_{\mathrm{th}}\left(\mathrm{Th}{F}_{0.7}\right)=N_{\mathrm{Th}}^c\left({\mathrm{ThF}}_{0.7}\right)\\ S_{\mathrm{Uu}}\·C_u\left({\mathrm{ThF}}_{0.05}\right)+S_{\mathrm{Uth}}\·C_{\mathrm{th}}\left({\mathrm{ThF}}_{0.05}\right)=N_U^c\left({\mathrm{ThF}}_{0.05}\right)\\ S_{\mathrm{Thu}}\·C_u\left({\mathrm{ThF}}_{0.05}\right)+S_{\mathrm{Thth}}\·C_{\mathrm{th}}\left({\mathrm{ThF}}_{0.05}\right)=N_{\mathrm{Th}}^c\left({\mathrm{ThF}}_{0.05}\right)\end{array}\right.$

式中，

C_u(UF_0.5)表示铀基准模型UF-0.5-Ⅰ的铀含量，单位%；

C_th(UF_0.5)表示铀基准模型UF-0.5-Ⅰ的钍含量，单位%；

C_u(UF_0.2)表示铀基准模型UF-0.2-Ⅰ的铀含量，单位%；

C_th(UF_0.2)表示铀基准模型UF-0.2-Ⅰ的钍含量，单位%；

C_u(ThF_0.7)表示钍基准模型ThF-0.7-Ⅰ的铀含量，单位%；

C_th(ThF_0.7)表示钍基准模型ThF-0.7-Ⅰ的钍含量，单位%；

C_u(ThF_0.05)表示钍基准模型ThF-0.05-Ⅰ的铀含量，单位%；

C_th(ThF_0.05)表示钍基准模型ThF-0.05-Ⅰ的钍含量，单位%；

表示经死时间T_UU校正后铀基准模型UF-0.5-Ⅰ上的铀窗计数率，单位：s^-1；

表示经死时间T_UTh校正后铀基准模型UF-0.5-Ⅰ上的钍窗计数率，单位：s^-1；

表示经死时间T_UU校正后铀基准模型UF-0.2-Ⅰ上的铀窗计数率，单位：s^-1；

表示经死时间T_UTh校正后铀基准模型UF-0.2-Ⅰ上的钍窗计数率，单位：s^-1；

表示经死时间T_ThU校正后钍基准模型ThF-0.7-Ⅰ上的铀窗计数率，单位：s^-1；

表示经死时间T_ThTh校正后钍基准模型ThF-0.7-Ⅰ上的钍窗计数率，单位：s^-1；

表示经死时间T_ThU校正后钍基准模型ThF-0.05-Ⅰ上的铀窗计数率，单位：s^-1；

表示经死时间T_ThTh校正后钍基准模型ThF-0.05-Ⅰ上的钍窗计数率，单位：s^-1。

步骤4.校验混合校验模型

根据在混合校验模型上采集的伽玛谱线和各灵敏度系数，通过下式计算混合校验模型的铀、钍含量：

$\left\{\begin{array}{c}{C}_u^c=\frac{N_U^c\·S_{\mathrm{Thth}}-N_{\mathrm{Th}}^c\·S_{\mathrm{Uth}}}{S_{\mathrm{Uu}}\·S_{\mathrm{Thth}}-S_{\mathrm{Uth}}\·S_{\mathrm{Thu}}}\\ C_{\mathrm{th}}^c=\frac{N_{\mathrm{Th}}^c\·S_{\mathrm{Uu}}-N_U^c\·S_{\mathrm{Thu}}}{S_{\mathrm{Uu}}\·S_{\mathrm{Thth}}-S_{\mathrm{Thu}}\·S_{\mathrm{Uth}}}\end{array}\right..$

式中，

表示计算获得的混合校验模型的铀含量，单位：%；

表示计算获得的混合校验模型的钍含量，单位：%；

通过下式计算铀、钍含量误差：

$\left\{\begin{array}{c}{\mathrm{\η}}_u=\frac{C_u^c-C_u}{C_u}\×100\%\\ {\mathrm{\η}}_{\mathrm{th}}=\frac{C_{\mathrm{th}}^c-C_{\mathrm{th}}}{C_{\mathrm{th}}}\×100\%\end{array}\right.$

式中，

C_u—混合校验模型的标称铀含量，单位：%；

C_th—混合校验模型的标称钍含量，单位：%；

η_u—铀含量相对误差，单位：%；

η_th—钍含量相对误差，单位：%。

如表1所示，核工业放射性勘查计量站共有2个混合校验模型，一般若每个模型的铀含量相对误差和钍含量相对误差均小于10%，即视为校正通过。

Claims (10)

1.一种用于铀矿勘探的伽玛能谱测井仪校正方法，其特征在于：包括以下步骤：首先采用伽玛能谱测井仪在铀基准模型、钍基准模型和混合校验模型上进行伽玛能谱测量，计算铀基准模型和钍基准模型上各自的铀窗和钍窗死时间，然后根据经死时间校正后的铀窗和钍窗的计数率和铀、钍基准模型的铀、钍含量，计算伽玛能谱测井仪的灵敏度系数；之后根据铀基准模型和钍基准模型上各自的铀窗和钍窗死时间，计算混合校验模型的铀窗和钍窗死时间，然后计算混合校验模型上经死时间校正后的铀窗和钍窗计数率；最后采用上述计算出的灵敏度系数和混合校验模型的铀窗、钍窗计数率，计算混合校验模型的铀、钍含量，与其标称含量进行比较，计算出铀、钍含量相对误差。

2.根据权利要求1所述的用于铀矿勘探的伽玛能谱测井仪校正方法，其特征在于：具体包括以下步骤：

步骤1.设置伽玛能谱测井仪能窗参数并选取基准模型；

步骤2.计算并修正死时间；

步骤2.1.采用双源法和线性相关系数法计算铀基准模型上铀窗死时间T_UU和钍窗死时间T_UTh；

步骤2.2.采用与步骤2.1相同的方法计算钍基准模型上最终铀窗死时间T_ThU和最终钍窗死时间T_ThTh；

步骤2.3.计算混合校验模型的铀窗死时间T_MU和钍窗死时间T_MTh；

步骤2.4.计算混合校验模型上经死时间修正后铀窗计数率和钍窗计数率

步骤3.计算灵敏度系数；

步骤4.校验混合校验模型。

3.根据权利要求2所述的用于铀矿勘探的伽玛能谱测井仪校正方法，其特征在于：步骤1具体包括以下步骤：

对待校正的伽玛能谱测井仪，采用“宽两窗”，即U窗和Th窗；选取a个铀基准模型、b个钍基准模型和c个混合校验模型；采用核工业放射性勘查计量站的放射性测井基准模型作为基准模型。

4.根据权利要求3所述的用于铀矿勘探的伽玛能谱测井仪校正方法，其特征在于：步骤1中，a的取值范围为3～6；b的取值范围为2～5；c的取值范围为2～6。

5.根据权利要求3所述的用于铀矿勘探的伽玛能谱测井仪校正方法，其特征在于：步骤1中，a=4，b=3，c=2；

基准模型参数如下表所示：

。

6.根据权利要求3所述的用于铀矿勘探的伽玛能谱测井仪校正方法，其特征在于：步骤2具体包括以下步骤：

步骤2.1具体包括以下步骤：

步骤2.1.1.采用双源法计算铀基准模型上初始铀窗死时间τ_UU和初始钍窗死时间τ_UTh

利用下式计算铀基准模型上初始铀窗死时间τ_UU和初始钍窗死时间τ_UTh：

$\left\{\begin{array}{c}{\mathrm{\τ}}_{\mathrm{UU}}=\frac{C_{U1}\·N_{2U}-C_{U2}\·N_{1U}}{C_{U1}\·N_{2U}\·N_{1\mathrm{TC}}-C_{U2}\·N_{1U}\·N_{2\mathrm{TC}}}\\ {\mathrm{\τ}}_{\mathrm{UTh}}=\frac{C_{U1}\·N_{2\mathrm{Th}}-C_{U2}\·N_{1\mathrm{Th}}}{C_{U1}\·N_{2\mathrm{Th}}\·N_{1\mathrm{TC}}-C_{U2}\·N_{1\mathrm{Th}}\·N_{2\mathrm{TC}}}\end{array}\right.;$

式中，

τ_UU表示铀模型上计算的初始铀窗死时间，单位：s；

τ_UTh表示铀模型上计算的初始钍窗死时间，单位：s；

C_U1表示1#铀基准模型的铀含量，即铀的质量百分比，单位：%；

C_U2表示2#铀基准模型的铀含量，即铀的质量百分比，单位：%；

N_1U表示1#铀基准模型上采集的谱线U窗的计数率，单位：s^-1；

N_1Th表示1#铀基准模型上采集的谱线Th窗的计数率，单位：s^-1；

N_1TC表示1#铀基准模型上采集的谱线总道的计数率，单位：s^-1；

N_2U表示2#铀基准模型上采集的谱线U窗的总计数率，单位：s^-1；

N_2Th表示2#铀基准模型上采集的谱线Th窗的计数率，单位：s^-1；

N_2TC表示2#铀基准模型上采集的谱线总道的计数率，单位：s^-1；

所述1#铀基准模型、2#铀基准模型为步骤1中任意两个不相同的铀基准模型，即有种组合方式，进而得到组初始铀窗死时间和组初始钍窗死时间

步骤2.1.2.采用线性相关系数法计算铀基准模型上最终铀窗死时间T_UU和最终钍窗死时间T_UTh

步骤2.1.2.1.计算铀基准模型上最终铀窗死时间T_UU

设

采用下式计算铀基准模型上最终铀窗死时间T_UU：

$\left\{\begin{array}{c}{N}_{\mathrm{UUi}}^c\left(n_1\right)=\frac{N_{\mathrm{UUi}}}{1-N_{\mathrm{UTCi}}\·T_{d1s}\left(n_1\right)},(i=\mathrm{1,2}...a)\\ \stackrel{\‾}{N_{\mathrm{UU}}^c\left(n_1\right)}=\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{N}_{\mathrm{UUi}}^c\left(n_1\right)/a\\ \stackrel{\‾}{C_U}=\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{C}_{\mathrm{Ui}}/a\\ R_{\mathrm{UU}}\left(n_1\right)=\frac{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}(N_{\mathrm{UUi}}^c\left(n_1\right)-\stackrel{\‾}{N_{\mathrm{UU}}^c\left(n_1\right)})(C_{\mathrm{Ui}}-\stackrel{\‾}{C_U})}{\sqrt{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{(N_{\mathrm{UUi}}^c\left(n_1\right)-\stackrel{\‾}{N_{\mathrm{UU}}^c\left(n_1\right)})}^2}\sqrt{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{(C_{\mathrm{Ui}}-\stackrel{\‾}{C_U})}^2}}\\ R_{\mathrm{UU}}\left(l_1\right)=\max \left\{R_{\mathrm{UU}}\left(n_1\right)\right\}\\ T_{\mathrm{UU}}=T_{d1s}\left(l_1\right)=T_{d1s\min }+l_1\·\mathrm{\Δ}{t}_1\end{array}\right.$

式中，

表示第i个铀基准模型经n₁次死时间校正后的铀窗计数率，单位：s-1；

N_UUi表示第i个铀基准模型的铀窗原始计数率，单位：s^-1；

N_UTCi表示第i个铀基准模型的总道原始计数率，单位：s^-1；

C_Ui表示第i个铀基准模型的铀含量，即铀质量百分比，单位：%；

R_UU(n₁)表示所有铀基准模型经n₁次死时间校正后的铀窗计数率与铀含量之间的相关系数；

相关系数R_UU(n₁)的最大值对应的死时间即为铀基准模型上最终铀窗死时间T_UU；

步骤2.1.2.2.计算铀基准模型上钍窗死时间T_UTh

采用与步骤2.1.2.1相同的线性相关系数法方法计算钍窗死时间T_UTh；

步骤2.2.采用与步骤2.1相同的方法计算钍基准模型上最终铀窗死时间T_ThU和最终钍窗死时间T_ThTh；

步骤2.3具体包括以下步骤：

采用线性插值法计算混合校验模型的铀窗死时间T_MU和钍窗死时间T_MTh：

$\left\{\begin{array}{c}{P}_U=\frac{\underset{i=1}{\overset{a}{\mathrm{\Σ}}}{p}_{\mathrm{Ui}}}{a}\\ P_{\mathrm{Th}}=\frac{\underset{j=1}{\overset{b}{\mathrm{\Σ}}}{p}_{\mathrm{Thj}}}{b}\\ T_{\mathrm{MU}}=\frac{T_{\mathrm{UU}}-T_{\mathrm{ThU}}}{P_U-P_{\mathrm{Th}}}\·(P_M-P_{\mathrm{Th}})+T_{\mathrm{ThU}}\\ T_{\mathrm{MTh}}=\frac{T_{\mathrm{UTh}}-T_{\mathrm{ThTh}}}{P_U-P_{\mathrm{Th}}}\·(P_M-P_{\mathrm{Th}})+T_{\mathrm{ThTh}}\end{array}\right.;$

式中，

p_Ui表示第i个铀基准模型上铀窗与钍窗的原始计数率比值；

p_Thj表示第j个钍基准模型上铀窗与钍窗的原始计数率比值；

P_M表示混合校验模型上铀窗与钍窗的原始计数率比值；

T_MU表示混合校验模型铀窗死时间，单位：s；

T_MTh表示混合校验模型钍窗死时间，单位：s；

步骤2.4具体包括以下步骤：

采用下式计算和

$\left\{\begin{array}{c}{N}_{\mathrm{MU}}^c=\frac{N_{\mathrm{MU}}}{1-N_{\mathrm{MTC}}\·T_{\mathrm{MU}}}\\ N_{\mathrm{MTh}}^c=\frac{N_{\mathrm{MTh}}}{1-N_{\mathrm{MTC}}\·T_{\mathrm{MTh}}}\end{array}\right.;$

式中，

表示经死时间校正后铀窗计数率，单位s^-1；

表示经死时间校正后钍窗计数率，单位s^-1；

N_MTC表示混合校验模型上总道计数率，单位s^-1。

7.根据权利要求6所述的用于铀矿勘探的伽玛能谱测井仪校正方法，其特征在于：步骤2.1.2.1中，m₁为10的整数倍。

8.根据权利要求6所述的用于铀矿勘探的伽玛能谱测井仪校正方法，其特征在于：步骤3具体包括以下步骤：

对基准模型依据下式建立方程组，采用最小二乘法解方程组得到四个灵敏度系数值S_Uu、S_Thu、S_Uth和S_Thth：

$\left\{\begin{array}{c}{S}_{\mathrm{Uu}}\·C_{\mathrm{Mu}}+S_{\mathrm{Uth}}\·C_{\mathrm{Mth}}=N_{\mathrm{MU}}^c\\ S_{\mathrm{Thu}}\·C_{\mathrm{Mu}}+S_{\mathrm{Thth}}\·C_{\mathrm{Mth}}=N_{\mathrm{MTh}}^c\end{array}\right.;$

式中，

S_Uu表示单位含量的铀对铀窗产生的伽玛计数率，即铀对铀窗的灵敏度系数，单位：s^-1·1%eU；

S_Thu表示单位含量的铀对钍窗产生的伽玛计数率，即铀对钍窗的灵敏度系数，单位：s^-1·1%eU；

S_Uth表示单位含量的钍对铀窗产生的伽玛计数率，即钍对铀窗的灵敏度系数，单位：s^-1·1%Th；

S_Thth表示单位含量的钍对钍窗产生的伽玛计数率，即钍对钍窗的灵敏度系数，单位：s^-1·1%Th；

C_Mu表示混合校验模型的铀含量，单位：%；

C_Mth表示混合校验模型的钍含量，单位：%；

表示经死时间校正后铀窗计数率，单位s^-1；

表示经死时间校正后钍窗计数率，单位s^-1。

9.根据权利要求8所述的用于铀矿勘探的伽玛能谱测井仪校正方法，其特征在于：采用下表所示模型作为计算灵敏度系数的基准模型：

。

10.根据权利要求8所述的用于铀矿勘探的伽玛能谱测井仪校正方法，其特征在于：步骤4具体包括以下步骤：

通过下式计算混合校验模型的铀、钍含量：

$\left\{\begin{array}{c}{C}_u^c=\frac{N_U^c\·S_{\mathrm{Thth}}-N_{\mathrm{Th}}^c\·S_{\mathrm{Uth}}}{S_{\mathrm{Uu}}\·S_{\mathrm{Thth}}-S_{\mathrm{Uth}}\·S_{\mathrm{Thu}}}\\ C_{\mathrm{th}}^c=\frac{N_{\mathrm{Th}}^c\·S_{\mathrm{Uu}}-N_U^c\·S_{\mathrm{Thu}}}{S_{\mathrm{Uu}}\·S_{\mathrm{Thth}}-S_{\mathrm{Thu}}\·S_{\mathrm{Uth}}}\end{array}\right.;$

式中，

表示计算获得的混合校验模型的铀含量，单位：%；

表示计算获得的混合校验模型的钍含量，单位：%；

通过下式计算铀、钍含量误差：

$\left\{\begin{array}{c}{\mathrm{\η}}_u=\frac{C_u^c-C_u}{C_u}\×100\%\\ {\mathrm{\η}}_{\mathrm{th}}=\frac{C_{\mathrm{th}}^c-C_{\mathrm{th}}}{C_{\mathrm{th}}}\×100\%\end{array}\right.$

式中，

C_u—混合校验模型的标称铀含量，单位：%；

C_th—混合校验模型的标称钍含量，单位：%；

η_u—铀含量相对误差，单位：%；

η_th—钍含量相对误差，单位：%。