CN113250686A - 一种井下伽马能谱直接测铀方法及系统 - Google Patents
一种井下伽马能谱直接测铀方法及系统 Download PDFInfo
- Publication number
- CN113250686A CN113250686A CN202110578002.9A CN202110578002A CN113250686A CN 113250686 A CN113250686 A CN 113250686A CN 202110578002 A CN202110578002 A CN 202110578002A CN 113250686 A CN113250686 A CN 113250686A
- Authority
- CN
- China
- Prior art keywords
- uranium
- spectrum
- stripping
- radium
- characteristic peak
- 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
- 229910052770 Uranium Inorganic materials 0.000 title claims abstract description 244
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 title claims abstract description 244
- 238000001228 spectrum Methods 0.000 title claims abstract description 213
- 238000000034 method Methods 0.000 title claims abstract description 46
- 229910052705 radium Inorganic materials 0.000 claims abstract description 122
- HCWPIIXVSYCSAN-UHFFFAOYSA-N radium atom Chemical compound [Ra] HCWPIIXVSYCSAN-UHFFFAOYSA-N 0.000 claims abstract description 122
- 238000001514 detection method Methods 0.000 claims abstract description 23
- 229910004755 Cerium(III) bromide Inorganic materials 0.000 claims abstract description 13
- 238000012545 processing Methods 0.000 claims description 12
- 230000005477 standard model Effects 0.000 claims description 12
- 230000003595 spectral effect Effects 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 238000005755 formation reaction Methods 0.000 claims description 6
- 238000001730 gamma-ray spectroscopy Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 description 28
- 238000012937 correction Methods 0.000 description 9
- 230000002285 radioactive effect Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 229910052704 radon Inorganic materials 0.000 description 5
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 5
- ZSLUVFAKFWKJRC-IGMARMGPSA-N 232Th Chemical compound [232Th] ZSLUVFAKFWKJRC-IGMARMGPSA-N 0.000 description 4
- 229910014323 Lanthanum(III) bromide Inorganic materials 0.000 description 4
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 description 4
- 229910052776 Thorium Inorganic materials 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000005251 gamma ray Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 239000011591 potassium Substances 0.000 description 4
- 229910052700 potassium Inorganic materials 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- MOOUSOJAOQPDEH-UHFFFAOYSA-K cerium(iii) bromide Chemical compound [Br-].[Br-].[Br-].[Ce+3] MOOUSOJAOQPDEH-UHFFFAOYSA-K 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- XKUYOJZZLGFZTC-UHFFFAOYSA-K lanthanum(iii) bromide Chemical compound Br[La](Br)Br XKUYOJZZLGFZTC-UHFFFAOYSA-K 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 239000000306 component Substances 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- 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
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明涉及一种井下伽马能谱直接测铀方法及系统。该方法包括：根据CeBr3探测器井下伽马能谱探测方法，构建室内纯铀模型与室内铀镭混合模型；根据所述室内纯铀模型和所述室内铀镭混合模型，确定铀镭剥谱比例系数；根据所述铀镭剥谱比例系数采用迭代剥谱方法，确定能量特征峰净面积；获取不同地层铀含量时的效率刻度；根据所述能量特征峰净面积和所述效率刻度，确定铀含量。本发明能够提高铀资源勘查效率、降低勘查成本。
Description
技术领域
本发明涉及铀资源勘查领域，特别是涉及一种井下伽马能谱直接测铀方法及系统。
背景技术
我国铀资源勘查中铀矿定量化评估手段主要采用γ总量或γ能谱测井技术，γ总量测井技术受钍、钾的影响较大，必须结合室内化学分析铀与钍、钾关系，以及铀镭、镭氡修正后方可获取铀含量；而现有γ能谱测井技术虽然解决了钍、钾干扰的影响，具有明显优势，但由于测量的γ射线依然主要来自镭组的214Bi、214Pb等，必须进行铀镭平衡修正、镭氡平衡修正等，尤其在铀镭平衡关系易受破坏的砂岩型铀矿床，必须进行取心分析获取铀镭平衡系数、射气系数等参数加以修正，均属于间接测铀方法，而且岩心与钻孔周围的地质体可能存在差异，修正会引入一定的修正误差，且实际应用中大量取心分析工作，增加了铀资源勘查成本、加长了勘查周期，随着我国找矿和开发深度的增大，钻探施工时矿层的取心难度和成本也会越来越高。
针对上述铀资源勘查间接测铀方法存在的不足，前人研究了基于NaI(Tl)探测器的地面铀矿体直接测铀方法及基于LaBr3(Ce)探测器、利用234mPa特征射线的直接铀定量技术。但该相关研究存在一定的弊端：NaI(Tl)探测器的能量分辨能力不足，一般在137Cs放射源662keV能量特征峰的分辨率约为8.0％，难以获得高分辨率精细伽马能谱数据，难以将可用于直接铀定量的1001KeV能量特征峰与邻近的969KeV、1120keV能量特征峰区分开来；LaBr3(Ce)探测器能量分辨率在662keV能量处接近3.5％，但其探测器自身含有138La及227Ac的放射性影响，带来探测器自身在约1.4MeV能量处存在放射性本底，这对于低能1001keV能量特征峰的探测存在一定的影响，且探测器自身放射性本底会对铀资源勘查中地层低铀含量的准确测量带来较大影响。
另外，中子测铀技术是铀资源勘查中直接有效测量地层铀含量的技术，但该项技术存在两个弊端：一是高通量中子发生器关键技术我国尚未攻克，目前只能依靠进口核心部件，且该部件存在使用寿命，一般为几百小时，如果采用该技术进行大面积铀资源勘查测井应用，成本极高；二是该项在实际应用中测井速度较低，一般为0.5～1.0m/min，测井效率明显低于放射性测井技术，较低工作效率仍会带来高成本。
发明内容
本发明的目的是提供一种井下伽马能谱直接测铀方法及系统，能够提高铀资源勘查效率、降低勘查成本。
为实现上述目的，本发明提供了如下方案：
一种井下伽马能谱直接测铀方法包括：
根据CeBr3探测器井下伽马能谱探测方法，构建室内纯铀模型与室内铀镭混合模型；
根据所述室内纯铀模型和所述室内铀镭混合模型，确定铀镭剥谱比例系数；
根据所述铀镭剥谱比例系数采用迭代剥谱方法，确定能量特征峰净面积；
获取不同地层铀含量时的效率刻度；
根据所述能量特征峰净面积和所述效率刻度，确定铀含量。
可选地，所述根据所述室内纯铀模型和所述室内铀镭混合模型，确定铀镭剥谱比例系数，具体包括：
从所述室内铀镭混合模型的实测谱数据中剥离出纯镭谱数据；
确定所述室内纯铀模型的实测谱数据中纯铀谱数据；
根据所述纯铀谱数据和所述纯镭谱数据进行1001keV能量特征峰处铀镭比例关系计算，得到1001keV能量特征峰处铀镭剥谱比例系数。
可选地，所述根据所述铀镭剥谱比例系数采用迭代剥谱方法，确定能量特征峰净面积，具体包括：
在剥谱处理中利用2615keV剥离232Th在1001keV能量特征峰的面积，得到第一次剥谱剩余面积；
利用1460keV剥离40K在1001keV能量特征峰的面积，得到第二次剥谱剩余面积；
利用1120keV剥离226Ra在1001keV能量特征峰的面积，得到第三次剥谱剩余面积；
根据所述第三次剥谱剩余面积和所述铀镭剥谱比例系数，确定剥离后的1001keV能量特征峰处的剥谱剩余面积，所述剥谱剩余面积为能量特征峰净面积。
可选地，所述获取不同地层铀含量时的效率刻度，具体包括：
利用标准模型UF-0.01、UF-0.03、UF-0.1、UF-0.2、UF-0.5及UF-1.0确定不同地层铀含量时的效率刻度。
可选地，所述根据所述能量特征峰净面积和所述效率刻度，确定铀含量，具体包括：
将剥谱处理后的1001keV能量特征峰作为计算铀含量的目标峰，利用剥谱处理后的1001keV能量特征峰净面积进行不同地层铀含量时的效率刻度并计算地层铀含量。
一种井下伽马能谱直接测铀系统包括：
室内纯铀模型/室内铀镭混合模型构建模块，用于根据CeBr3探测器井下伽马能谱探测方法，构建室内纯铀模型与室内铀镭混合模型；
铀镭剥谱比例系数确定模块，用于根据所述室内纯铀模型和所述室内铀镭混合模型，确定铀镭剥谱比例系数；
能量特征峰净面积确定模块，用于根据所述铀镭剥谱比例系数采用迭代剥谱方法，确定能量特征峰净面积；
效率刻度获取模块，用于获取不同地层铀含量时的效率刻度；
铀含量确定模块，用于根据所述能量特征峰净面积和所述效率刻度，确定铀含量。
可选地，所述铀镭剥谱比例系数确定模块具体包括：
纯镭谱数据确定单元，用于从所述室内铀镭混合模型的实测谱数据中剥离出纯镭谱数据；
纯铀谱数据确定单元，用于确定所述室内纯铀模型的实测谱数据中纯铀谱数据；
铀镭剥谱比例系数确定单元，用于根据所述纯铀谱数据和所述纯镭谱数据进行1001keV能量特征峰处铀镭比例关系计算，得到1001keV能量特征峰处铀镭剥谱比例系数。
可选地，所述能量特征峰净面积确定模块具体包括：
第一次剥谱单元，用于在剥谱处理中利用2615keV剥离232Th在1001keV能量特征峰的面积，得到第一次剥谱剩余面积；
第二次剥谱单元，用于利用1460keV剥离40K在1001keV能量特征峰的面积，得到第二次剥谱剩余面积；
第三次剥谱单元，用于利用1120keV剥离226Ra在1001keV能量特征峰的面积，得到第三次剥谱剩余面积；
能量特征峰净面积确定单元，用于根据所述第三次剥谱剩余面积和所述铀镭剥谱比例系数，确定剥离后的1001keV能量特征峰处的剥谱剩余面积，所述剥谱剩余面积为能量特征峰净面积。
可选地，所述效率刻度获取模块具体包括：
效率刻度获取单元，用于利用标准模型UF-0.01、UF-0.03、UF-0.1、UF-0.2、UF-0.5及UF-1.0确定不同地层铀含量时的效率刻度。
可选地，所述铀含量确定模块具体包括：
铀含量确定单元，用于将剥谱处理后的1001keV能量特征峰作为计算铀含量的目标峰，利用剥谱处理后的1001keV能量特征峰净面积进行不同地层铀含量时的效率刻度并计算地层铀含量。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明无需铀镭平衡修正、镭氡平衡修正，现场可直接获取地层铀含量，提高了生产效率、降低成本，为铀资源勘查提供了一种高效、低成本的铀矿定量探测技术。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明井下伽马能谱直接测铀方法流程图；
图2为本发明井下伽马能谱直接测铀系统结构图；
图3为本发明构建新型CeBr3探测器的井下伽马能谱技术实测的铀镭混合γ谱图；
图4为本发明获取的铀镭剥谱数据图；
图5为本发明针对1001keV能量特征峰的迭代剥谱流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
针对我国铀资源勘查中铀矿定量化评估手段存在的不足，本发明提出一种井下伽马能谱直接测量地层铀含量的方法，通过构建基于新型高分辨率溴化铈(CeBr3)探测器的伽马能谱探测技术，以1001keV能量特征峰为目标峰，利用室内标准模型获取镭(226Ra)谱中1001keV能量特征峰与1120keV能量特征峰比例关系，进而通过铀镭钍钾混合谱数据剥谱获取1001keV能量特征峰的净面积，进而通过标准模型效率刻度来实现地层铀含量的测量。由于1001keV能量特征伽马射线来自234mPa，在漫长地质时期可认为234mPa与238U始终处于放射性平衡状态，所以提出的基于新型CeBr3探测器、通过1001keV能量特征峰测量地层238U含量的本发明，无需铀镭平衡修正、镭氡平衡修正便可实现井下伽马能谱直接测量地层铀含量。图1为本发明井下伽马能谱直接测铀方法流程图，如图1所示，本发明的井下伽马能谱直接测铀方法包括如下步骤：
步骤101：根据CeBr3探测器井下伽马能谱探测方法，构建室内纯铀模型与室内铀镭混合模型。
新型CeBr3探测器具备室温条件下的高能量分辨能力，且无自身放射性本底干扰，是本发明开展弱分支比1001keV能量特征伽马射线探测的硬件基础。基于该新型CeBr3探测器的井下伽马能谱探测技术相对现有基于NaI(Tl)、LaBr3(Ce)等探测器的探测技术优势明显，是本发明区别于现有技术的关键之一。
步骤102：根据所述室内纯铀模型和所述室内铀镭混合模型，确定铀镭剥谱比例系数，具体包括：
步骤1021：从所述室内铀镭混合模型的实测谱数据中剥离出纯镭谱数据；
步骤1022：确定所述室内纯铀模型的实测谱数据中纯铀谱数据；
步骤1023：根据所述纯铀谱数据和所述纯镭谱数据进行1001keV能量特征峰处铀镭比例关系计算，得到1001keV能量特征峰处铀镭剥谱比例系数。
利用室内纯铀模型与室内铀镭混合模型获取了镭(226Ra)谱中1001keV能峰与1120keV能峰比例关系。该比例关系是实现1001keV能量特征峰剥谱处理的关键，也是区别于现有NaI(Tl)、LaBr3(Ce)等探测器的地方，其可将1001KeV能峰与邻近969KeV、1120keV能峰区分开来，方可具备建立剥谱比例系数的条件。另外，该比例系数也是步骤103中迭代剥谱获取1001keV能峰净面积的基础，通过该比例关系可以将铀镭谱数据中的镭剥离出去。
步骤103：根据所述铀镭剥谱比例系数采用迭代剥谱方法，确定能量特征峰净面积，具体包括：
步骤1031：在剥谱处理中利用2615keV剥离232Th在1001keV能量特征峰的面积，得到第一次剥谱剩余面积；
步骤1032：利用1460keV剥离40K在1001keV能量特征峰的面积，得到第二次剥谱剩余面积；
步骤1033：利用1120keV剥离226Ra在1001keV能量特征峰的面积，得到第三次剥谱剩余面积；
步骤1034：根据所述第三次剥谱剩余面积和所述铀镭剥谱比例系数，确定剥离后的1001keV能量特征峰处的剥谱剩余面积，所述剥谱剩余面积为能量特征峰净面积。
步骤104：获取不同地层铀含量时的效率刻度，具体包括：
利用标准模型UF-0.01、UF-0.03、UF-0.1、UF-0.2、UF-0.5及UF-1.0确定不同地层铀含量时的效率刻度。
标准模型是已有的，进行不同地层铀含量的效率刻度在新型CeBr3探测器系统中是需要去做的工作。
步骤105：根据所述能量特征峰净面积和所述效率刻度，确定铀含量，具体包括：
将剥谱处理后的1001keV能量特征峰作为计算铀含量的目标峰，利用剥谱处理后的1001keV能量特征峰净面积进行不同地层铀含量时的效率刻度并计算地层铀含量。
“处理后的1001keV能量特征峰”这个是指测量地层铀含量采用的目标峰或者特征能量射线为1001keV，实际上还有其他可用峰，比如766keV、93keV等；“能量特征峰净面积”是表示一个效率刻度计算时的变量了，限定效率刻度的利用对象。
图2为本发明井下伽马能谱直接测铀系统结构图。如图2所示，一种井下伽马能谱直接测铀系统包括：
室内纯铀模型/室内铀镭混合模型构建模块201，用于根据CeBr3探测器井下伽马能谱探测方法，构建室内纯铀模型与室内铀镭混合模型；
铀镭剥谱比例系数确定模块202，用于根据所述室内纯铀模型和所述室内铀镭混合模型，确定铀镭剥谱比例系数；
能量特征峰净面积确定模块203，用于根据所述铀镭剥谱比例系数采用迭代剥谱方法，确定能量特征峰净面积；
效率刻度获取模块204，用于获取不同地层铀含量时的效率刻度；
铀含量确定模块205，用于根据所述能量特征峰净面积和所述效率刻度，确定铀含量。
所述铀镭剥谱比例系数确定模块202具体包括：
纯镭谱数据确定单元，用于从所述室内铀镭混合模型的实测谱数据中剥离出纯镭谱数据；
纯铀谱数据确定单元，用于确定所述室内纯铀模型的实测谱数据中纯铀谱数据；
铀镭剥谱比例系数确定单元，用于根据所述纯铀谱数据和所述纯镭谱数据进行1001keV能量特征峰处铀镭比例关系计算，得到1001keV能量特征峰处铀镭剥谱比例系数。
所述能量特征峰净面积确定模块203具体包括：
第一次剥谱单元，用于在剥谱处理中利用2615keV剥离232Th在1001keV能量特征峰的面积，得到第一次剥谱剩余面积；
第二次剥谱单元，用于利用1460keV剥离40K在1001keV能量特征峰的面积，得到第二次剥谱剩余面积；
第三次剥谱单元，用于利用1120keV剥离226Ra在1001keV能量特征峰的面积，得到第三次剥谱剩余面积；
能量特征峰净面积确定单元，用于根据所述第三次剥谱剩余面积和所述铀镭剥谱比例系数，确定剥离后的1001keV能量特征峰处的剥谱剩余面积，所述剥谱剩余面积为能量特征峰净面积。
所述效率刻度获取模块204具体包括：
效率刻度获取单元，用于利用标准模型UF-0.01、UF-0.03、UF-0.1、UF-0.2、UF-0.5及UF-1.0确定不同地层铀含量时的效率刻度。
所述铀含量确定模块205具体包括：
铀含量确定单元，用于将剥谱处理后的1001keV能量特征峰作为计算铀含量的目标峰，利用剥谱处理后的1001keV能量特征峰净面积进行不同地层铀含量时的效率刻度并计算地层铀含量。
实施例1：
本实施例提供一种井下伽马能谱直接测铀方法，具体包括：
步骤一：构建新型CeBr3探测器的井下伽马能谱探测技术。新型CeBr3探测器具备室温条件下的高能量分辨能力，且无自身放射性本底干扰，是本发明开展弱分支比1001keV能量特征伽马射线探测的硬件基础。构建的新型CeBr3探测器的井下伽马能谱探测技术，为满足井下钻孔直径限制采用1英寸CeBr3晶体探测器，并采用怀特射极电压放大电路进行信号前置放大处理，其成形时间为1微秒，信号采集采用65MHz高速ADC，为满足井下铀矿体高强度辐照下的γ射线探测需求，多道技术为基于FPGA的快慢双通道成形技术，快成形时间约为0.6微秒、慢成形时间为2微秒，可满足铀矿体铀含量eU1.0％以内的采集需求。另外，基于怀特射极电压放大电路与FPGA信号成形处理技术，实现的新型CeBr3探测器的井下伽马能谱探头在137Cs放射源661.6keV的能量分辨率约为4.0％，其实测γ谱如图3所示，在1001keV具有明显响应谱峰。
步骤二：获取1001keV能峰处铀镭剥谱比例系数。首先利用步骤一构建的新型CeBr3探测器的井下伽马能谱探测硬件及技术，实测室内纯铀模型与铀镭混合模型；然后，利用铀镭混合模型与纯铀模型含量关系，从铀镭混合模型实测谱数据中剥离出纯镭谱，利用纯铀谱与纯镭谱进行1001keV能峰处铀镭比例关系计算，该比例系数即为1001keV能峰处铀镭剥谱比例系数。该比例系数也是步骤三中迭代剥谱获取1001keV能峰净面积的关键，通过该比例关系可以将铀镭谱数据中的镭剥离出去。利用室内模型获取该比例系数，并用于铀镭混合谱中镭的剥离，是本发明的创新之一。
步骤三：针对1001keV能量特征峰采用迭代剥谱的方法获取能峰净面积。利用步骤一的硬件系统及技术、步骤二获取的剥谱比例系数，针对1001keV能峰开展迭代剥谱处理，图4为本发明获取的铀镭剥谱数据图，图5为本发明针对1001keV能量特征峰的迭代剥谱流程图。
第1步：利用步骤一构建的新型CeBr3井下γ能谱技术获取原始谱数据；
第2步：利用2615keV能峰计算232Th，并利用2615keV能峰与1460keV、1120keV及1001keV的比例关系，进行第1次迭代剥谱处理；
第3步：利用1460keV能峰计算40K，并利用1460keV能峰与1120keV及1001keV的比例关系，进行第2次迭代剥谱处理；
第4步：利用1120keV能峰计算226Ra，并利用步骤二中获取的铀镭剥谱比例系数，开展第3次迭代剥谱处理；
第5步：利用前4步迭代剥谱后1001keV能峰谱数据计算1001keV能峰净面积；
步骤四：利用1001keV能量特征峰净面积进行效率刻度与计算。将剥谱处理后的1001keV能量特征峰作为计算铀含量的目标峰，在国防科技工业1313二级计量站标准模型井进行效率刻度。利用标准模型UF-0.01、UF-0.03、UF-0.1、UF-0.2、UF-0.5及UF-1.0开展不同地层铀含量时的效率刻度，获取刻度系数，进而实现在实际应用时利用该系数与步骤三方法获取的1001keV能峰净面积进行铀含量的计算，进而实现本发明的井下伽马能谱直接测铀方法。
本发明与现有NaI(Tl)、LaBr3(Ce)探测器的γ能谱探测技术相比，是一种采用新型CeBr3探测器的井下γ能谱直接测铀技术，创新性通过室内标准模型的建立了1001keV能量特征峰的剥谱系数及迭代剥谱技术，最终实现利用1001keV能量特征峰来测量地层铀含量的目的。本发明无需铀镭平衡修正、镭氡平衡修正，现场可直接获取地层铀含量，提高了生产效率、降低成本，为铀资源勘查提供了一种高效、低成本的铀矿定量探测技术。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种井下伽马能谱直接测铀方法，其特征在于，包括：
根据CeBr3探测器井下伽马能谱探测方法，构建室内纯铀模型与室内铀镭混合模型；
根据所述室内纯铀模型和所述室内铀镭混合模型，确定铀镭剥谱比例系数；
根据所述铀镭剥谱比例系数采用迭代剥谱方法，确定能量特征峰净面积；
获取不同地层铀含量时的效率刻度；
根据所述能量特征峰净面积和所述效率刻度，确定铀含量。
2.根据权利要求1所述的井下伽马能谱直接测铀方法，其特征在于，所述根据所述室内纯铀模型和所述室内铀镭混合模型，确定铀镭剥谱比例系数，具体包括：
从所述室内铀镭混合模型的实测谱数据中剥离出纯镭谱数据；
确定所述室内纯铀模型的实测谱数据中纯铀谱数据；
根据所述纯铀谱数据和所述纯镭谱数据进行1001keV能量特征峰处铀镭比例关系计算，得到1001keV能量特征峰处铀镭剥谱比例系数。
3.根据权利要求1所述的井下伽马能谱直接测铀方法，其特征在于，所述根据所述铀镭剥谱比例系数采用迭代剥谱方法，确定能量特征峰净面积，具体包括：
在剥谱处理中利用2615keV剥离232Th在1001keV能量特征峰的面积，得到第一次剥谱剩余面积；
利用1460keV剥离40K在1001keV能量特征峰的面积，得到第二次剥谱剩余面积；
利用1120keV剥离226Ra在1001keV能量特征峰的面积，得到第三次剥谱剩余面积；
根据所述第三次剥谱剩余面积和所述铀镭剥谱比例系数，确定剥离后的1001keV能量特征峰处的剥谱剩余面积，所述剥谱剩余面积为能量特征峰净面积。
4.根据权利要求1所述的井下伽马能谱直接测铀方法，其特征在于，所述获取不同地层铀含量时的效率刻度，具体包括：
利用标准模型UF-0.01、UF-0.03、UF-0.1、UF-0.2、UF-0.5及UF-1.0确定不同地层铀含量时的效率刻度。
5.根据权利要求1所述的井下伽马能谱直接测铀方法，其特征在于，所述根据所述能量特征峰净面积和所述效率刻度，确定铀含量，具体包括：
将剥谱处理后的1001keV能量特征峰作为计算铀含量的目标峰，利用剥谱处理后的1001keV能量特征峰净面积进行不同地层铀含量时的效率刻度并计算地层铀含量。
6.一种井下伽马能谱直接测铀系统，其特征在于，包括：
室内纯铀模型/室内铀镭混合模型构建模块，用于根据CeBr3探测器井下伽马能谱探测方法，构建室内纯铀模型与室内铀镭混合模型；
铀镭剥谱比例系数确定模块，用于根据所述室内纯铀模型和所述室内铀镭混合模型，确定铀镭剥谱比例系数；
能量特征峰净面积确定模块，用于根据所述铀镭剥谱比例系数采用迭代剥谱方法，确定能量特征峰净面积；
效率刻度获取模块，用于获取不同地层铀含量时的效率刻度；
铀含量确定模块，用于根据所述能量特征峰净面积和所述效率刻度，确定铀含量。
7.根据权利要求6所述的井下伽马能谱直接测铀系统，其特征在于，所述铀镭剥谱比例系数确定模块具体包括：
纯镭谱数据确定单元，用于从所述室内铀镭混合模型的实测谱数据中剥离出纯镭谱数据；
纯铀谱数据确定单元，用于确定所述室内纯铀模型的实测谱数据中纯铀谱数据；
铀镭剥谱比例系数确定单元，用于根据所述纯铀谱数据和所述纯镭谱数据进行1001keV能量特征峰处铀镭比例关系计算，得到1001keV能量特征峰处铀镭剥谱比例系数。
8.根据权利要求6所述的井下伽马能谱直接测铀系统，其特征在于，所述能量特征峰净面积确定模块具体包括：
第一次剥谱单元，用于在剥谱处理中利用2615keV剥离232Th在1001keV能量特征峰的面积，得到第一次剥谱剩余面积；
第二次剥谱单元，用于利用1460keV剥离40K在1001keV能量特征峰的面积，得到第二次剥谱剩余面积；
第三次剥谱单元，用于利用1120keV剥离226Ra在1001keV能量特征峰的面积，得到第三次剥谱剩余面积；
能量特征峰净面积确定单元，用于根据所述第三次剥谱剩余面积和所述铀镭剥谱比例系数，确定剥离后的1001keV能量特征峰处的剥谱剩余面积，所述剥谱剩余面积为能量特征峰净面积。
9.根据权利要求6所述的井下伽马能谱直接测铀系统，其特征在于，所述效率刻度获取模块具体包括：
效率刻度获取单元，用于利用标准模型UF-0.01、UF-0.03、UF-0.1、UF-0.2、UF-0.5及UF-1.0确定不同地层铀含量时的效率刻度。
10.根据权利要求6所述的井下伽马能谱直接测铀系统，其特征在于，所述铀含量确定模块具体包括：
铀含量确定单元，用于将剥谱处理后的1001keV能量特征峰作为计算铀含量的目标峰，利用剥谱处理后的1001keV能量特征峰净面积进行不同地层铀含量时的效率刻度并计算地层铀含量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110578002.9A CN113250686B (zh) | 2021-05-26 | 2021-05-26 | 一种井下伽马能谱直接测铀方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110578002.9A CN113250686B (zh) | 2021-05-26 | 2021-05-26 | 一种井下伽马能谱直接测铀方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113250686A true CN113250686A (zh) | 2021-08-13 |
| CN113250686B CN113250686B (zh) | 2023-07-11 |
Family
ID=77184759
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110578002.9A Active CN113250686B (zh) | 2021-05-26 | 2021-05-26 | 一种井下伽马能谱直接测铀方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113250686B (zh) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0143529A1 (en) * | 1983-10-05 | 1985-06-05 | Halliburton Company | Multi-function natural gamma ray logging |
| US20050199794A1 (en) * | 2004-03-15 | 2005-09-15 | Medhat Mickael | Spectral gamma ray logging-while-drilling system |
| CN102094639A (zh) * | 2009-12-14 | 2011-06-15 | 西安威尔罗根能源科技有限公司 | 重晶石泥浆井自然伽马能谱测井校正的方法 |
| EP2561182A2 (en) * | 2010-04-19 | 2013-02-27 | Baker Hughes Incorporated | Formation evaluation using a bit-based active radiation source and a gamma ray detector |
| CN107605471A (zh) * | 2017-07-21 | 2018-01-19 | 中国石油天然气股份有限公司 | 一种湖相碳酸盐岩岩性测井定量识别方法 |
| CN107656305A (zh) * | 2017-08-15 | 2018-02-02 | 中国原子能科学研究院 | 一种测定铀富集度的方法 |
| CN108457644A (zh) * | 2017-12-26 | 2018-08-28 | 中国石油天然气股份有限公司 | 一种元素俘获能谱测井的伽马能谱解谱方法及装置 |
| CN111045072A (zh) * | 2019-12-27 | 2020-04-21 | 核工业北京地质研究院 | 一种适用于CeBr3探测器的γ能谱迭代解谱算法 |
| CN112523741A (zh) * | 2020-11-24 | 2021-03-19 | 东华理工大学 | 一种基于能谱测井交叉谱段的铀矿定量刻度系数求法 |
| CN112630849A (zh) * | 2020-11-24 | 2021-04-09 | 东华理工大学 | 一种基于能谱测井特征谱峰的铀矿定量剥离系数求法 |
-
2021
- 2021-05-26 CN CN202110578002.9A patent/CN113250686B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0143529A1 (en) * | 1983-10-05 | 1985-06-05 | Halliburton Company | Multi-function natural gamma ray logging |
| US20050199794A1 (en) * | 2004-03-15 | 2005-09-15 | Medhat Mickael | Spectral gamma ray logging-while-drilling system |
| CN102094639A (zh) * | 2009-12-14 | 2011-06-15 | 西安威尔罗根能源科技有限公司 | 重晶石泥浆井自然伽马能谱测井校正的方法 |
| EP2561182A2 (en) * | 2010-04-19 | 2013-02-27 | Baker Hughes Incorporated | Formation evaluation using a bit-based active radiation source and a gamma ray detector |
| CN107605471A (zh) * | 2017-07-21 | 2018-01-19 | 中国石油天然气股份有限公司 | 一种湖相碳酸盐岩岩性测井定量识别方法 |
| CN107656305A (zh) * | 2017-08-15 | 2018-02-02 | 中国原子能科学研究院 | 一种测定铀富集度的方法 |
| CN108457644A (zh) * | 2017-12-26 | 2018-08-28 | 中国石油天然气股份有限公司 | 一种元素俘获能谱测井的伽马能谱解谱方法及装置 |
| CN111045072A (zh) * | 2019-12-27 | 2020-04-21 | 核工业北京地质研究院 | 一种适用于CeBr3探测器的γ能谱迭代解谱算法 |
| CN112523741A (zh) * | 2020-11-24 | 2021-03-19 | 东华理工大学 | 一种基于能谱测井交叉谱段的铀矿定量刻度系数求法 |
| CN112630849A (zh) * | 2020-11-24 | 2021-04-09 | 东华理工大学 | 一种基于能谱测井特征谱峰的铀矿定量剥离系数求法 |
Non-Patent Citations (1)
| Title |
|---|
| 马嫱;肖伟剑;陈勇飞;林碧芬;蔡继杰;: "Gamma能谱无源效率刻度方法在稀土样品放射性检测中的应用", 稀土, no. 01 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113250686B (zh) | 2023-07-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN1158541C (zh) | 测定地层密度的方法及仪器 | |
| Caciolli et al. | A new FSA approach for in situ γ ray spectroscopy | |
| CN108457640B (zh) | 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 | |
| CN106250619B (zh) | 一种确定地层矿物含量的方法和装置 | |
| CN107229080B (zh) | 一种元素测井中子俘获伽马能谱的获取方法 | |
| EP2831646A1 (en) | Neutron porosity based on one or more gamma ray detectors and a pulsed neutron source | |
| CN105093343A (zh) | 一种地层元素俘获伽马能谱的测量方法及装置 | |
| CN103678778B (zh) | 一种放射性物化探信息综合方法 | |
| CN107144868A (zh) | 一种现场多参数快速调查评价煤田放射性辐射环境的方法 | |
| CN111045072A (zh) | 一种适用于CeBr3探测器的γ能谱迭代解谱算法 | |
| CN104074501B (zh) | 一种裂变中子测井修正方法 | |
| Tchorz-Trzeciakiewicz et al. | Seasonal variations of terrestrial gamma dose, natural radionuclides and human health | |
| CN103345002A (zh) | 一种测井资料获取砂岩卤水层钾离子含量的方法 | |
| CN104111482A (zh) | 一种双探测器x射线荧光测井探管及方法 | |
| CN113250686A (zh) | 一种井下伽马能谱直接测铀方法及系统 | |
| CN112630849B (zh) | 一种基于能谱测井特征谱峰的铀矿定量剥离系数求法 | |
| CN112523741B (zh) | 一种基于能谱测井交叉谱段的铀矿定量刻度系数求法 | |
| Petrovic et al. | Scintillators for PGNAA in mineral exploration | |
| Verdoya et al. | Natural gamma-ray spectrometry as a tool for radiation dose and radon hazard modelling | |
| CN112523742A (zh) | 一种基于能谱测井特征谱段的铀矿定量刻度系数求法 | |
| CN106018231B (zh) | 一种原位探测水蚀土壤中137Cs渗透深度的方法和系统 | |
| CN112649888B (zh) | 一种基于能谱测井特征谱峰的铀矿定量刻度系数求法 | |
| AU2022339061A1 (en) | Method for assessing the uranium content in a borehole using gamma spectrometry and associated device | |
| CN204496010U (zh) | 一种地质野外露头伽马仪 | |
| CN115614022B (zh) | 铀裂变瞬发中子测井的井眼水层原位修正方法 |
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 |