一种基于能谱测井交叉谱段的铀矿定量剥离系数求法
阅读说明：本技术 一种基于能谱测井交叉谱段的铀矿定量剥离系数求法 (Uranium ore quantitative stripping coefficient solving method based on energy spectrum logging cross spectrum section ) 是由 王海涛 汤彬 王仁波 刘志锋 黄凡 张丽娇 周书民 张雄杰 张焱 陈锐 刘琦 于 2020-11-24 设计创作，主要内容包括：本发明公开了一种基于能谱测井交叉谱段的铀矿定量剥离系数求法。剥离系数可表示为将单位含量的其它元素在某个特征谱段产生的计数率剥离掉的系数。本发明所述基于能谱测井交叉谱段的铀矿定量剥离系数求法是以自然伽马能谱曲线中钍、铀-镭、钾元素所对应的特征峰为对象,将自然伽马能谱曲线划分为多个相互交叉的特征谱段,求得单位含量的不同元素在某个特征谱段的计数率,两种元素在该特征谱段计数率的比值即为剥离系数。本发明所述基于能谱测井交叉谱段的铀矿定量剥离系数求法,可在确保钍、铀-镭、钾三种自然伽马放射性元素含量解释精度的同时,使自然伽马能谱测井的速度达到与自然伽马总量测井相当的水平。(The invention discloses a uranium ore quantitative stripping coefficient calculation method based on an energy spectrum logging cross spectrum section. The stripping factor can be expressed as a factor that strips off the count rate produced by other elements of the unit content in a certain characteristic spectral band. The method for determining the quantitative stripping coefficient of the uranium ore based on the energy spectrum logging cross spectrum section takes characteristic peaks corresponding to elements of thorium, uranium-radium and potassium in a natural gamma energy spectrum curve as objects, the natural gamma energy spectrum curve is divided into a plurality of mutually crossed characteristic spectrum sections, the counting rates of different elements with unit content in a certain characteristic spectrum section are obtained, and the ratio of the counting rates of the two elements in the characteristic spectrum section is the stripping coefficient. The uranium ore quantitative stripping coefficient solving method based on the energy spectrum logging cross-spectral band can ensure the interpretation precision of the contents of three natural gamma radioactive elements of thorium, uranium-radium and potassium, and simultaneously enables the logging speed of the natural gamma energy spectrum to reach the level equivalent to the logging speed of the natural gamma total amount.)
技术领域
本发明属于核辐射探测领域，应用该方法可在铀矿探采行业中通过快速的自然伽马能谱测井实现放射性元素定量，特别适用于铀钍混合型或铀钍钾混合型矿产。
背景技术
自然伽马测井是一种常见的钻井地球物理方法，也是铀矿勘探的基本方法。它是通过探测天然衰变系(铀系、钍系、锕铀系等)和钾(40K)的伽马射线总量或能谱计数率(它们均正比于衰变率)，以此推算地层岩石中由起始核素表征的铀、钍或钾元素的含量。自然伽马测井是将伽马总量测井仪或者伽马能谱测井仪放置到钻孔中，测量井壁岩矿石的天然伽马照射量率，并根据沿井深的伽马照射量率曲线确定钻孔所穿过的放射性地层位置、厚度，以及放射性元素(铀、钍、钾)含量。目前，在铀矿床、铀钍混合矿床的勘探中，伽马测井已成为储量计算的主要方法，特别是在钻孔中的岩芯采取率不高时，基于定量放射性元素的伽马测井显得尤为重要。中国的《γ测井规范》仅要求采用技术成熟的自然伽马总量测井方法，以此作为地层岩石中铀定量的主要依据。
自然界中的铀系和钍系在放射性平衡状态下，系内核素的原子核数的比例关系是确定的，因此不同能量伽马射线的相对强度也是确定的，可以分别在两个系中选出某种核素的特征核素伽马射线的能量来识别铀和钍。特征核素所发射的伽马射线的能量称为特征能量，在自然伽马能谱测井中，如选用铀系中的214Bi发射的1.76MeV的伽马射线来识别铀，选用钍系中的208Tl发射的2.62MeV的伽马射线来识别钍，用1.46MeV的伽马射线来识别钾。如果把伽马射线按所选定的特征能量分别计数，即为测谱。把这些粒子发射的伽马射线的能量画在坐标系中，横坐标表示为伽马射线的能量，纵坐标表示为相应的该能量的伽马射线的强度，就得到了伽马射线的能量和强度的关系图，这个图称为自然伽马射线的能谱图或能谱曲线图。因而将测得的自然伽马能谱转换成地层的铀、钍、钾的含量，并以连续测井曲线的形式输出，就是自然伽马能谱测井。
与自然伽马总量测井相比，自然伽马能谱测井不仅可以实现总量测井的功能，而且可以获取更多的有用信息，确定地层中铀、钍、钾的含量，以便更详细地划分地层和研究与放射性元素分布有关的多种地质问题。
相比于自然伽马总量测井，自然伽马能谱测井的相对计数率(铀、钍、钾计数率)低，放射性统计涨落误差大，为了提高曲线质量，必须加大伽马射线探测器(晶体)体积、降低测速，这往往又与实际生产要求相冲突。而且由于地层中铀、钍、钾的含量不同，导致铀、钍、钾曲线的质量也不同。自然伽马能谱测井曲线的质量不但取决于测井仪本身的性能及技术水平，而且还受测井环境(井眼和地层)、测速、采样间隔等因素的影响。
目前急需研究一种新的铀矿领域的快速自然伽马能谱测井方法，可在确保放射性元素含量解释精度的同时，使自然伽马能谱的测井速度达到与自然伽马总量测井相当的水平，满足生产应用要求。基于交叉谱段法的自然伽马能谱测井方法有望解决这一问题，本发明所述的一种基于能谱测井交叉谱段的铀矿定量剥离系数求法是实现快速自然伽马能谱测井方法的关键。迄今为止，尚未见到将该项方法直接应用于铀矿自然能谱测井的报道。
发明内容
本发明的目的是为了解决在铀矿探采行业中通过快速的自然伽马能谱测井实现放射性元素定量，而提出的一种基于能谱测井交叉谱段的铀矿定量剥离系数求法。
本发明的目的是通过下述技术方案实现的：
(1)根据不同放射性元素的特征峰所在位置，将自然伽马能谱曲线其划分为多个相互交叉的特征谱段：
1)选取钍系的i个特征峰，划分为对应的i个钍系特征谱段，每个特征谱段的伽马射线能量范围均从400keV起至对应的特征峰为止，且至少有一个特征谱段应包含从400keV至能量为2.62MeV的钍系特征峰，
2)选取铀-镭系的j个特征峰，划分为对应的j个铀-镭系特征谱段，每个特征谱段的伽马射线能量范围均从400keV起至对应的特征峰为止，且至少有一个特征谱段应包含400keV至能量为1.76MeV的铀-镭系特征峰，
3)选取钾系的k个特征峰，划分为对应的k个钾系特征谱段，每个特征谱段的伽马射线能量范围均从400keV起至对应的特征峰为止，且至少有一个特征谱段应包含400keV至能量为1.46MeV的钾系特征峰；
(2)在自然伽马放射性标准模型上求取钍系、铀-镭系、钾系的各个特征谱段的计数率：
1)在自然伽马放射性本底标准模型上求取钍系特征谱段的计数率铀-镭系特征谱段的计数率钾系特征谱段的计数率
2)在标称含量为QTh的自然伽马放射性钍元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QTh的自然伽马放射性钍元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QTh的自然伽马放射性钍元素标准模型上求取钾元素所对应的各个特征谱段的计数率
3)在标称含量为QU的自然伽马放射性铀元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QU的自然伽马放射性铀元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QU的自然伽马放射性铀元素标准模型上求取钾元素所对应的各个特征谱段的计数率
4)在标称含量为QK的自然伽马放射性钾元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QK的自然伽马放射性钾元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QK的自然伽马放射性钾元素标准模型上求取钾元素所对应的各个特征谱段的计数率
(3)求取基于交叉谱段的自然伽马能谱剥离系数：
1)钍元素对自身各个特征谱段的剥离系数均为1、铀-镭元素对自身各个特征谱段的剥离系数均为1、钾元素对自身各个特征谱段的剥离系数均为1；
2)利用自然伽马放射性本底标准模型、标称含量为QTh的自然伽马放射性钍元素标准模型和标称含量为QU的自然伽马放射性铀元素标准模型，求取铀-镭元素对钍元素i个特征谱段的剥离系数：
3)利用自然伽马放射性本底标准模型、标称含量为QTh的自然伽马放射性钍元素标准模型和标称含量为QK的自然伽马放射性钾元素标准模型，求取钾元素对钍元素i个特征谱段的剥离系数：
4)利用自然伽马放射性本底标准模型、标称含量为QU的自然伽马放射性铀元素标准模型和标称含量为QTh的自然伽马放射性钍元素标准模型，求取钍元素对铀-镭元素j个特征谱段的剥离系数：
5)利用自然伽马放射性本底标准模型、标称含量为QU的自然伽马放射性铀-镭元素标准模型和标称含量为QK的自然伽马放射性钾元素标准模型，求取钾元素对铀-镭元素j个特征谱段的剥离系数：
6)利用自然伽马放射性本底标准模型、标称含量为QK的自然伽马放射性钾元素标准模型和标称含量为QTh的自然伽马放射性钍元素标准模型，求取钍元素对钾元素k个特征谱段的剥离系数：
7)利用自然伽马放射性本底标准模型、标称含量为QK的自然伽马放射性钾元素标准模型和标称含量为QU的自然伽马放射性铀-镭元素标准模型，求取铀-镭元素对钾元素k个特征谱段的剥离系数：
8)综合步骤(3)中的步骤1)-7)，求得基于交叉谱段的自然伽马能谱剥离系数：
可用公式表示为：
其中x、y表示任意的自然伽马放射性元素，x/y表示元素x对元素y，m表示各个特征谱段，m＝i+j+k，m∈y。
本发明的优点在于：利用基于交叉谱段法的自然伽马能谱测井剥离系数求取方法，能够在进行某种放射性元素含量解析过程中，剥离其他自然伽马放射性元素的影响，实现钍、铀-镭、钾等放射性元素的准确定量；同时由于采用交叉谱段法，有效利用了自然伽马能谱曲线的测量计数率，相对提高了自然伽马能谱曲线的信噪比，进而可提高自然伽马能谱测量的速度。若应用到铀矿自然伽马测井过程中，可以在确保放射性元素含量解释精度的同时，使自然伽马能谱测井的速度达到与自然伽马总量测井相同的水平，满足生产应用要求。
附图说明
图1为本发明实施例1的剥离系数求取流程图；
图2为本发明实施例1的含有钍、铀-镭、钾放射性元素自然伽马能谱曲线的交叉谱段划分方法示例；
图3为本发明实施例1的针对铀矿测井中钍、铀-镭、钾放射性元素定量过程示意图；
图4为本发明实施例1的在同一个井孔中不同测井速度下的铀-镭含量解释结果对比图；
图5为本发明实施例1的在铀含量为800ppm的模型井中采用自然伽马总量测井和自然伽马能谱测井时的铀-镭含量解释结果对比图。
图中：1-自然伽马能谱测井曲线，2-划分钍、铀-镭、钾的相互交叉的特征谱段，3-选取本底模型和已知含量的钍、铀、钾模型，4-各模型所测自然伽马能谱曲线数据，5-各交叉谱段的计数率，6-剥离系数。
具体实施方式
下面通过附图和具体实施方式对本发明作更为详细的描述。
本发明的基本思想是以自然伽马能谱曲线中钍、铀-镭、钾元素所对应的特征峰为对象，将自然伽马能谱曲线划分为多个相互交叉的特征谱段，求得单位含量的不同元素在某个特征谱段的计数率，两种元素在该特征谱段计数率的比值即为剥离系数。
本发明所述的基于能谱测井交叉谱段的铀矿定量剥离系数求法，步骤如下：
(1)钍系、铀-镭系、和钾元素含有的伽马核素并不多，然而它们放出上百种能量的特征伽马射线。其中，放射几率较大、能量较高的特征伽马射线及其对应的伽马核素如表1所示。现有自然伽马能谱测井仅仅能够分辨出为数不多的十几种特征伽马射线，即放射几率>0.01、能量>0.4MeV且无重叠峰的特征伽马射线。这些特征伽马射线在能谱曲线上称为特征峰，如图2所示。根据不同放射性元素的特征峰所在位置，将自然伽马能谱曲线其划分为多个相互交叉的特征谱段：
1)选取钍系的i个特征峰，划分为对应的i个钍系特征谱段，每个特征谱段的伽马射线能量范围均从400keV起至对应的特征峰为止，且至少有一个特征谱段应包含从400keV至能量为2.62MeV的钍系特征峰，
2)选取铀-镭系的j个特征峰，划分为对应的j个铀-镭系特征谱段，每个特征谱段的伽马射线能量范围均从400keV起至对应的特征峰为止，且至少有一个特征谱段应包含400keV至能量为1.76MeV的铀-镭系特征峰，
3)选取钾系的k个特征峰，划分为对应的k个钾系特征谱段，每个特征谱段的伽马射线能量范围均从400keV起至对应的特征峰为止，且至少有一个特征谱段应包含400keV至能量为1.46MeV的钾系特征峰；
(2)在自然伽马放射性标准模型上求取钍系、铀-镭系、钾系的各个特征谱段的计数率：
1)在自然伽马放射性本底标准模型上求取钍系特征谱段的计数率铀-镭系特征谱段的计数率钾系特征谱段的计数率
2)在标称含量为QTh的自然伽马放射性钍元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QTh的自然伽马放射性钍元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QTh的自然伽马放射性钍元素标准模型上求取钾元素所对应的各个特征谱段的计数率
3)在标称含量为QU的自然伽马放射性铀元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QU的自然伽马放射性铀元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QU的自然伽马放射性铀元素标准模型上求取钾元素所对应的各个特征谱段的计数率
4)在标称含量为QK的自然伽马放射性钾元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QK的自然伽马放射性钾元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QK的自然伽马放射性钾元素标准模型上求取钾元素所对应的各个特征谱段的计数率
(3)求取基于交叉谱段的自然伽马能谱剥离系数：
1)钍元素对自身各个特征谱段的剥离系数均为1、铀-镭元素对自身各个特征谱段的剥离系数均为1、钾元素对自身各个特征谱段的剥离系数均为1；
2)利用自然伽马放射性本底标准模型、标称含量为QTh的自然伽马放射性钍元素标准模型和标称含量为QU的自然伽马放射性铀元素标准模型，求取铀-镭元素对钍元素i个特征谱段的剥离系数：
3)利用自然伽马放射性本底标准模型、标称含量为QTh的自然伽马放射性钍元素标准模型和标称含量为QK的自然伽马放射性钾元素标准模型，求取钾元素对钍元素i个特征谱段的剥离系数：
4)利用自然伽马放射性本底标准模型、标称含量为QU的自然伽马放射性铀元素标准模型和标称含量为QTh的自然伽马放射性钍元素标准模型，求取钍元素对铀-镭元素j个特征谱段的剥离系数：
5)利用自然伽马放射性本底标准模型、标称含量为QU的自然伽马放射性铀-镭元素标准模型和标称含量为QK的自然伽马放射性钾元素标准模型，求取钾元素对铀-镭元素j个特征谱段的剥离系数：
6)利用自然伽马放射性本底标准模型、标称含量为QK的自然伽马放射性钾元素标准模型和标称含量为QTh的自然伽马放射性钍元素标准模型，求取钍元素对钾元素k个特征谱段的剥离系数：
7)利用自然伽马放射性本底标准模型、标称含量为QK的自然伽马放射性钾元素标准模型和标称含量为QU的自然伽马放射性铀-镭元素标准模型，求取铀-镭元素对钾元素k个特征谱段的剥离系数：
8)综合步骤(3)中的步骤1)-7)，求得基于交叉谱段的自然伽马能谱剥离系数：
可用公式表示为：
其中x、y表示任意的自然伽马放射性元素，x/y表示元素x对元素y，m表示各个特征谱段，m＝i+j+k，m∈y。
(4)求取基于交叉谱段的自然伽马能谱换算系数：
1)利用自然伽马放射性本底标准模型和标称含量为QTh的自然伽马放射性钍元素标准模型，求取钍元素对应各个特征谱段的换算系数：
2)利用自然伽马放射性本底标准模型和标称含量为QU的自然伽马放射性铀-镭元素标准模型，求取铀-镭元素对应各个特征谱段的换算系数：
3)利用自然伽马放射性本底标准模型和标称含量为QK的自然伽马放射性钾元素标准模型，求取钾元素对应各个特征谱段的换算系数：
综合步骤(4)中的步骤1)-3)，求得基于交叉谱段的自然伽马能谱换算系数：
可用公式表示为：
其中y表示任意的自然伽马放射性元素，Qy表示放射性标准模型元素y的标称含量，m表示各个特征谱段，m＝i+j+k，m∈y。
(5)根据如下公式求取放射性元素的含量qx：
可用公式表示为：
其中x、y表示任意的自然伽马放射性元素(分别为钍、铀-镭和钾)，x/y表示元素x对元素y，qx表示不同的放射性元素含量，m表示各个特征谱段，m＝i+j+k，m∈y。
若i、j、k都为1，则钍元素、铀-镭元素、钾元素都选用一个特征谱段，由实施例1的图3所示过程完成自然伽马能谱测井并求得钍、铀-镭、钾的放射性元素含量，得到在同一井孔中不同测井速度情况下的铀-镭放射性元素解释含量对比，效果如图4所示。并通过在模型井中分别进行自然伽马总量测井和基于交叉谱段法的自然伽马能谱测井，得到总量测井解释结果和能谱测井解释结果中钍、铀-镭、钾的放射性元素含量对比，效果如图5所示。从图4、图5的对比效果可看出，通过利用基于交叉谱段法的自然伽马能谱测井剥离系数求取方法，当自然伽马能谱测井速度达到6m/min时，仍可实现钍、铀-镭、钾等放射性元素的准确定量，满足生产应用要求。
表1.天然放射性衰变的伽马核素数据表(仅列出放射几率与能量均较高的特征伽马射线)
注：仅仅放射几率>0.001(指单次放射性衰变的放射几率)、能量>0.4MeV的钍系、铀-镭系和钾放出的特征伽马射线及其伽马核素的数据被列于表中。
- 上一篇：一种医用注射器针头装配设备
- 下一篇：移动侦测的防误触发方法、装置、设备及存储介质