CN112523742B - 一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法 - Google Patents
一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法 Download PDFInfo
- Publication number
- CN112523742B CN112523742B CN202011330400.0A CN202011330400A CN112523742B CN 112523742 B CN112523742 B CN 112523742B CN 202011330400 A CN202011330400 A CN 202011330400A CN 112523742 B CN112523742 B CN 112523742B
- Authority
- CN
- China
- Prior art keywords
- natural gamma
- characteristic
- spectrum
- uranium
- characteristic spectrum
- 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.)
- Active
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 190
- 230000002285 radioactive effect Effects 0.000 title claims abstract description 97
- 238000000034 method Methods 0.000 title claims abstract description 38
- ZSLUVFAKFWKJRC-IGMARMGPSA-N 232Th Chemical compound [232Th] ZSLUVFAKFWKJRC-IGMARMGPSA-N 0.000 claims abstract description 91
- 229910052776 Thorium Inorganic materials 0.000 claims abstract description 91
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 claims abstract description 77
- 239000011591 potassium Substances 0.000 claims abstract description 77
- 229910052700 potassium Inorganic materials 0.000 claims abstract description 77
- DLFWIFNRAUYTHF-UHFFFAOYSA-N [Ra].[U] Chemical compound [Ra].[U] DLFWIFNRAUYTHF-UHFFFAOYSA-N 0.000 claims abstract description 65
- 229910052770 Uranium Inorganic materials 0.000 claims abstract description 32
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 claims abstract description 32
- 230000005477 standard model Effects 0.000 claims description 75
- 238000006243 chemical reaction Methods 0.000 claims description 13
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 6
- 239000011707 mineral Substances 0.000 abstract description 6
- 229920006395 saturated elastomer Polymers 0.000 abstract description 4
- VBWSWBQVYDBVGA-NAHFVJFTSA-N uranium-234;uranium-235;uranium-238 Chemical compound [234U].[235U].[238U] VBWSWBQVYDBVGA-NAHFVJFTSA-N 0.000 description 10
- 230000005251 gamma ray Effects 0.000 description 9
- 238000005259 measurement Methods 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000001730 gamma-ray spectroscopy Methods 0.000 description 3
- 238000011002 quantification Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- GFRMDONOCHESDE-UHFFFAOYSA-N [Th].[U] Chemical compound [Th].[U] GFRMDONOCHESDE-UHFFFAOYSA-N 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000000084 gamma-ray spectrum Methods 0.000 description 2
- 230000005258 radioactive decay Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 150000001218 Thorium Chemical class 0.000 description 1
- 150000001224 Uranium Chemical class 0.000 description 1
- YPHLCJMCHHXHTR-UHFFFAOYSA-N [K].[Th].[U] Chemical compound [K].[Th].[U] YPHLCJMCHHXHTR-UHFFFAOYSA-N 0.000 description 1
- 229910052767 actinium Inorganic materials 0.000 description 1
- QQINRWTZWGJFDB-UHFFFAOYSA-N actinium atom Chemical compound [Ac] QQINRWTZWGJFDB-UHFFFAOYSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- ZLMJMSJWJFRBEC-OUBTZVSYSA-N potassium-40 Chemical compound [40K] ZLMJMSJWJFRBEC-OUBTZVSYSA-N 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
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
- E21B47/00—Survey of boreholes 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
- 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
-
- 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)
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法。该方法是以自然伽马能谱曲线中钍、铀‑镭、钾元素所对应的特征峰为对象，将自然伽马能谱曲线划分为多个特征谱段，所述刻度系数表示在饱和矿层中刻度出来的由测井仪响应的一个常数，表征饱和矿层中单位含量的某种放射性元素所放出的全部伽马射线在各个特征谱段响应的计数率。本发明所述基于能谱测井特征谱段的铀矿定量刻度系数求法，可在确保钍、铀‑镭、钾三种自然伽马放射性元素含量解析精度的同时，实现相对快速地自然伽马能谱测井。
Description
技术领域
本发明属于核辐射探测领域，应用该方法可在铀矿探采行业中通过快速的自然伽马能谱测井实现放射性元素定量，特别适用于铀钍混合型或铀钍钾混合型矿产。
背景技术
自然伽马测井是一种常见的钻井地球物理方法，也是铀矿勘探的基本方法。它是通过探测天然衰变系(铀系、钍系、锕铀系等)和钾(40K)的伽马射线总量或能谱计数率(它们均正比于衰变率)，以此推算地层岩石中由起始核素表征的铀、钍或钾元素的含量。自然伽马测井是将伽马总量测井仪或者伽马能谱测井仪放置到钻孔中，测量井壁岩矿石的天然伽马照射量率，并根据沿井深的伽马照射量率曲线确定钻孔所穿过的放射性地层位置、厚度，以及放射性元素(铀、钍、钾)含量。目前，在铀矿床、铀钍混合矿床的勘探中，伽马测井已成为储量计算的主要方法，特别是在钻孔中的岩芯采取率不高时，基于定量放射性元素的伽马测井显得尤为重要。中国的《γ测井规范》仅要求采用技术成熟的自然伽马总量测井方法，以此作为地层岩石中铀定量的主要依据。
自然界中的铀系和钍系在放射性平衡状态下，系内核素的原子核数的比例关系是确定的，因此不同能量伽马射线的相对强度也是确定的，可以分别在两个系中选出某种核素的特征核素伽马射线的能量来识别铀和钍。特征核素所发射的伽马射线的能量称为特征能量，在自然伽马能谱测井中，如选用铀系中的214Bi发射的1.76MeV的伽马射线来识别铀，选用钍系中的208Tl发射的2.62MeV的伽马射线来识别钍，用1.46MeV的伽马射线来识别钾。如果把伽马射线按所选定的特征能量分别计数，即为测谱。把这些粒子发射的伽马射线的能量画在坐标系中，横坐标表示为伽马射线的能量，纵坐标表示为相应的该能量的伽马射线的强度，就得到了伽马射线的能量和强度的关系图，这个图称为自然伽马射线的能谱图或能谱曲线图。因而将测得的自然伽马能谱转换成地层的铀、钍、钾的含量，并以连续测井曲线的形式输出，就是自然伽马能谱测井。
与自然伽马总量测井相比，自然伽马能谱测井不仅可以实现总量测井的功能，而且可以获取更多的有用信息，确定地层中铀、钍、钾的含量，以便更详细地划分地层和研究与放射性元素分布有关的多种地质问题。
相比于自然伽马总量测井，自然伽马能谱测井的相对计数率(铀、钍、钾计数率)低，放射性统计涨落误差大，为了提高曲线质量，必须加大伽马射线探测器(晶体)体积、降低测速，这往往又与实际生产要求相冲突。而且由于地层中铀、钍、钾的含量不同，导致铀、钍、钾曲线的质量也不同。自然伽马能谱测井曲线的质量不但取决于测井仪本身的性能及技术水平，而且还受测井环境(井眼和地层)、测速、采样间隔等因素的影响。
目前急需研究一种新的铀矿领域的快速自然伽马能谱测井方法，可以在确保测量精度的同时，实现相对快速地自然伽马能谱测井，满足生产应用要求。基于特征谱段法的自然伽马能谱测井方法有望解决这一问题，本发明所述的一种基于能谱测井特征谱段是实现快速自然伽马能谱测井方法的关键。迄今为止，尚未见到将该项方法直接应用于铀矿自然能谱测井的报道。
发明内容
本发明的目的是为了解决在铀矿探采行业中通过快速的自然伽马能谱测井实现放射性元素定量，而提出的一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法。
本发明的目的是通过下述技术方案实现的：
(1)根据不同放射性元素的特征峰所在位置，将自然伽马能谱曲线其划分为m个特征谱段，分别为：选取钍系的i个特征峰将其扩展成i个对应的钍系特征谱段，选取铀-镭系的j个特征峰将其扩展成j个对应的钍系特征谱段，选取钾的k个特征峰将其扩展成k个对应的钍系特征谱段，其中m＝i+j+k。具体方法为：
1)钍系特征谱段从伽马射线能量为400keV起，至少至包括能量为2.62MeV的钍系特征峰为止，选取i个钍系特征峰，以特征峰为基准或将其能量范围适当展宽形成i个钍系特征谱段，
2)铀-镭系特征谱段从伽马射线能量为400keV起，至少至包括能量为1.76MeV的铀-镭系特征峰为止，但不包括能量为2.62MeV的钍系特征峰，选取j个铀-镭系特征峰，以特征峰为基准或将其能量范围适当展宽形成j个铀-镭系特征谱段，
3)钾系特征谱段从伽马射线能量为400keV起，至少至包括能量为1.46MeV的钾系特征峰为止，但不包括能量为1.76MeV的铀-镭系特征峰，选取k个钾系特征峰，以特征峰为基准或将其能量范围适当展宽形成k个钾系特征谱段；
(2)在自然伽马放射性本底标准模型上求取m个特征谱段的计数率分别在钍元素、铀-镭元素和钾元素的自然伽马放射性标准模型上求取对应的m个特征谱段的计数率其中x表示自然伽马放射性元素钍、铀-镭、钾中的一种。具体方法为：
1)在自然伽马放射性本底标准模型上求取钍系特征谱段的计数率铀-镭系特征谱段的计数率/>钾特征谱段的计数率/>
2)在标称含量为QTh的自然伽马放射性钍元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QTh的自然伽马放射性钍元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QTh的自然伽马放射性钍元素标准模型上求取钾元素所对应的各个特征谱段的计数率
3)在标称含量为QU的自然伽马放射性铀元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QU的自然伽马放射性铀元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QU的自然伽马放射性铀元素标准模型上求取钾元素所对应的各个特征谱段的计数率
4)在标称含量为QK的自然伽马放射性钾元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QK的自然伽马放射性钾元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QK的自然伽马放射性钾元素标准模型上求取钾元素所对应的各个特征谱段的计数率
(3)求取基于能谱测井特征谱段的铀矿定量能谱剥离系数，元素x对元素y在特征谱段m上的剥离系数可表示为：
其中x、y表示任意的自然伽马放射性元素，x/y表示元素x对元素y，Qx表示放射性标准模型元素x的标称含量，m表示各个特征谱段，m＝i+j+k，m∈y。
具体过程为：
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)，求得基于特征谱段的自然伽马能谱剥离系数：
可用公式表示为：
(4)求取基于能谱测井特征谱段的铀矿定量换算系数，元素y在特征谱段m上的换算系数可表示为：
其中y表示任意的自然伽马放射性元素，Qy表示放射性标准模型元素y的标称含量，m表示各个特征谱段，m＝i+j+k，m∈y；
具体步骤为：
1)利用自然伽马放射性本底标准模型、标称含量为QTh的自然伽马放射性钍元素标准模型，求取钍元素对应各个特征谱段的换算系数：
2)利用自然伽马放射性本底标准模型、标称含量为QU的自然伽马放射性铀-镭元素标准模型，求取铀-镭元素对应各个特征谱段的换算系数：
3)利用自然伽马放射性本底标准模型、标称含量为QK的自然伽马放射性钾元素标准模型，求取钾元素对应各个特征谱段的换算系数：
可用矩阵方式可表示为：
(5)求取基于能谱测井特征谱段的铀矿定量刻度系数：
其中n、x、y表示任意的自然伽马放射性元素，n＝x+y，m表示各个特征谱段，m＝i+j+k，m∈y。
可用矩阵方式可表示为：
(6)根据如下公式求取自然伽马放射性元素的含量qx：
可用公式表示为：
其中x、y表示任意的自然伽马放射性元素(分别为钍、铀-镭和钾)，x/y表示元素x对元素y，m表示各个特征谱段，m＝i+j+k，m∈y。
本发明的优点在于：利用基于能谱测井特征谱段求取自然伽马放射性元素含量的方法，能够在进行某种放射性元素含量解析过程中，剥离其他自然伽马放射性元素的影响，实现钍、铀-镭、钾等放射性元素的准确定量；同时由于采用特征谱段法，有效利用了自然伽马能谱曲线的测量计数率，相对提高了自然伽马能谱曲线的信噪比，进而可提高自然伽马能谱测量的速度。若应用到铀矿自然伽马能谱测井过程中，可在确保放射性元素含量解释精度的同时，实现相对快速地自然伽马能谱测井，满足生产应用要求。
附图说明
图1为本发明实施例1的自然伽马放射性元素含量求取流程图；
图2为本发明实施例1的含有钍、铀-镭、钾放射性元素自然伽马能谱曲线的特征谱段划分方法示例1；
图3为本发明实施例1的含有钍、铀-镭、钾放射性元素自然伽马能谱曲线的特征谱段划分方法示例2；
图4为本发明实施例1的针对铀矿测井中钍、铀-镭、钾放射性元素定量过程示意图；
图5为本发明实施例1的在同一个井孔中不同测井速度下的铀-镭含量解释结果对比图。
图6为本发明实施例1的在铀含量为800ppm的模型井中采用自然伽马总量测井和自然伽马能谱测井时的铀-镭含量解释结果对比图。
图中：1-自然伽马能谱测井曲线，2-根据特征峰划分特征谱段，3-选取本底模型和已知含量的钍、铀、钾模型，4-各模型所测自然伽马能谱曲线数据，5-各特征谱段的计数率，6-剥离系数，7-换算系数，8-刻度系数，9-自然伽马放射性元素含量。
具体实施方式
下面通过附图和具体实施方式对本发明作更为详细的描述。
本发明的基本思想是以自然伽马能谱曲线中钍、铀-镭、钾元素所对应的特征峰为对象，将自然伽马能谱曲线划分为多个特征谱段，刻度系数表示在饱和矿层中刻度出来的由测井仪响应的一个常数，表征饱和矿层中单位含量的某种放射性元素所放出的全部伽马射线在各个特征谱段响应的计数率。
本发明所述的基于能谱测井特征谱段求取自然伽马放射性元素含量的方法，步骤如下：
(1)根据不同放射性元素的特征峰所在位置，将自然伽马能谱曲线其划分为m个特征谱段，分别为：选取钍系的i个特征峰将其扩展成i个对应的钍系特征谱段，选取铀-镭系的j个特征峰将其扩展成j个对应的钍系特征谱段，选取钾的k个特征峰将其扩展成k个对应的钍系特征谱段，其中m＝i+j+k。具体方法为：
1)钍系特征谱段从伽马射线能量为400keV起，至少至包括能量为2.62MeV的钍系特征峰为止，选取i个钍系特征峰，以特征峰为基准或将其能量范围适当展宽形成i个钍系特征谱段，
2)铀-镭系特征谱段从伽马射线能量为400keV起，至少至包括能量为1.76MeV的铀-镭系特征峰为止，但不包括能量为2.62MeV的钍系特征峰，选取j个铀-镭系特征峰，以特征峰为基准或将其能量范围适当展宽形成j个铀-镭系特征谱段，
3)钾系特征谱段从伽马射线能量为400keV起，至少至包括能量为1.46MeV的钾系特征峰为止，但不包括能量为1.76MeV的铀-镭系特征峰，选取k个钾系特征峰，以特征峰为基准或将其能量范围适当展宽形成k个钾系特征谱段；
(2)在自然伽马放射性本底标准模型上求取m个特征谱段的计数率分别在钍元素、铀-镭元素和钾元素的自然伽马放射性标准模型上求取对应的m个特征谱段的计数率其中x表示自然伽马放射性元素钍、铀-镭、钾中的一种。具体方法为：
1)在自然伽马放射性本底标准模型上求取钍系特征谱段的计数率铀-镭系特征谱段的计数率/>钾特征谱段的计数率/>
2)在标称含量为QTh的自然伽马放射性钍元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QTh的自然伽马放射性钍元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QTh的自然伽马放射性钍元素标准模型上求取钾元素所对应的各个特征谱段的计数率
3)在标称含量为QU的自然伽马放射性铀元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QU的自然伽马放射性铀元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QU的自然伽马放射性铀元素标准模型上求取钾元素所对应的各个特征谱段的计数率
4)在标称含量为QK的自然伽马放射性钾元素标准模型上求取钍元素所对应的各个特征谱段的计数率
在标称含量为QK的自然伽马放射性钾元素标准模型上求取铀-镭元素所对应的各个特征谱段的计数率
在标称含量为QK的自然伽马放射性钾元素标准模型上求取钾元素所对应的各个特征谱段的计数率
(3)求取基于能谱测井特征谱段的铀矿定量能谱剥离系数，元素x对元素y在特征谱段m上的剥离系数可表示为：
其中x、y表示任意的自然伽马放射性元素，x/y表示元素x对元素y，Qx表示放射性标准模型元素x的标称含量，m表示各个特征谱段，m＝i+j+k，m∈y。
具体过程为：
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)，求得基于特征谱段的自然伽马能谱剥离系数：
可用公式表示为：
(4)求取基于能谱测井特征谱段的铀矿定量换算系数，元素y在特征谱段m上的换算系数可表示为：
其中y表示任意的自然伽马放射性元素，Qy表示放射性标准模型元素y的标称含量，m表示各个特征谱段，m＝i+j+k，m∈y；
具体步骤为：
1)利用自然伽马放射性本底标准模型、标称含量为QTh的自然伽马放射性钍元素标准模型，求取钍元素对应各个特征谱段的换算系数：
2)利用自然伽马放射性本底标准模型、标称含量为QU的自然伽马放射性铀-镭元素标准模型，求取铀-镭元素对应各个特征谱段的换算系数：
3)利用自然伽马放射性本底标准模型、标称含量为QK的自然伽马放射性钾元素标准模型，求取钾元素对应各个特征谱段的换算系数：
可用矩阵方式可表示为：
(5)求取基于能谱测井特征谱段的铀矿定量刻度系数：
其中n、x、y表示任意的自然伽马放射性元素，n＝x+y，m表示各个特征谱段，m＝i+j+k，m∈y。
可用矩阵方式可表示为：
(6)根据如下公式求取自然伽马放射性元素的含量qx：
可用公式表示为：
其中x、y表示任意的自然伽马放射性元素(分别为钍、铀-镭和钾)，x/y表示元素x对元素y，m表示各个特征谱段，m＝i+j+k，m∈y。
若i、j、k都为1，则钍元素、铀-镭元素、钾元素都选用一个特征谱段，由实施例1的图4所示过程完成自然伽马能谱测井并求得钍、铀-镭、钾的放射性元素含量，得到在同一井孔中不同测井速度情况下的铀-镭放射性元素解释含量对比，效果如图5所示。并通过在模型井中分别进行自然伽马总量测井和基于特征谱段法的自然伽马能谱测井，得到总量测井解释结果和能谱测井解释结果中钍、铀-镭、钾的放射性元素含量对比，效果如图6所示。从图5、图6的对比效果可看出，通过利用基于能谱测井特征谱段求取放射性元素含量的方法，当自然伽马能谱测井速度达到4m/min时，仍可实现钍、铀-镭、钾等放射性元素的准确定量，满足生产应用要求。
表1.天然放射性衰变的伽马核素数据表(仅列出放射几率与能量均较高的特征伽马射线)
注：仅仅放射几率>0.001(指单次放射性衰变的放射几率)、能量>0.4MeV的钍系、铀-镭系和钾放出的特征伽马射线及其伽马核素的数据被列于表中。
Claims (1)
1.一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法，步骤如下：
(1)根据不同放射性元素的特征峰所在位置，将自然伽马能谱曲线其划分为m个特征谱段，分别为：选取钍系的i个特征峰将其扩展成i个对应的钍系特征谱段，选取铀-镭系的j个特征峰将其扩展成j个对应的钍系特征谱段，选取钾的k个特征峰将其扩展成k个对应的钍系特征谱段，其中m＝i+j+k；
(2)在自然伽马放射性本底标准模型上求取m个特征谱段的计数率分别在钍元素、铀-镭元素和钾元素的自然伽马放射性标准模型上求取对应的m个特征谱段的计数率/>其中x表示自然伽马放射性元素钍、铀-镭、钾中的一种；
(3)求取基于能谱测井特征谱段的铀矿定量能谱剥离系数，元素x对元素y在特征谱段m上的剥离系数可表示为：
其中x、y表示任意的自然伽马放射性元素，x/y表示元素x对元素y，Qx表示放射性标准模型元素x的标称含量，m表示各个特征谱段，m＝i+j+k，m∈y；
(4)求取基于能谱测井特征谱段的铀矿定量换算系数，元素y在特征谱段m上的换算系数可表示为：
其中y表示任意的自然伽马放射性元素，Qy表示放射性标准模型元素y的标称含量，m表示各个特征谱段，m＝i+j+k，m∈y；
(5)求取基于能谱测井特征谱段的铀矿定量刻度系数：
其中n、x、y表示任意的自然伽马放射性元素，n＝x+y，m表示各个特征谱段，m＝i+j+k，m∈y；
(6)根据如下公式求取自然伽马放射性元素的含量qx：
可用公式表示为：
其中x、y表示任意的自然伽马放射性元素，分别为钍、铀-镭和钾，x/y表示元素x对元素y，m表示各个特征谱段，m＝i+j+k，m∈y。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011330400.0A CN112523742B (zh) | 2020-11-24 | 2020-11-24 | 一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011330400.0A CN112523742B (zh) | 2020-11-24 | 2020-11-24 | 一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112523742A CN112523742A (zh) | 2021-03-19 |
| CN112523742B true CN112523742B (zh) | 2023-11-17 |
Family
ID=74993382
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011330400.0A Active CN112523742B (zh) | 2020-11-24 | 2020-11-24 | 一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112523742B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113359203B (zh) * | 2021-06-18 | 2024-01-30 | 西安石油大学 | 一种基于自然伽马能谱测井探测深部喷流沉积型矿床的方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3940610A (en) * | 1974-02-15 | 1976-02-24 | Mobil Oil Corporation | Natural gamma radiation borehole logging system |
| US4542292A (en) * | 1982-04-01 | 1985-09-17 | Schlumberger Technology Corporation | Correction of natural gamma radiation logs for the effects of gamma ray emission from and attenuation by the borehole fluid |
| CN103306663A (zh) * | 2012-03-06 | 2013-09-18 | 中国原子能科学研究院 | 铀矿测井方法 |
| CN107533150A (zh) * | 2015-04-27 | 2018-01-02 | 桑德克斯有线有限公司 | 用于以谱的伽马射线测井确定元素浓度的方法 |
| CN108457640A (zh) * | 2018-01-26 | 2018-08-28 | 东华理工大学 | 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 |
| CN108825219A (zh) * | 2018-05-25 | 2018-11-16 | 东华理工大学 | 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法 |
| CN109557120A (zh) * | 2019-01-14 | 2019-04-02 | 东华理工大学 | 伽马能谱与主动式x荧光联合同时测铀钍镭钾的方法 |
| CN111045072A (zh) * | 2019-12-27 | 2020-04-21 | 核工业北京地质研究院 | 一种适用于CeBr3探测器的γ能谱迭代解谱算法 |
| CN111335882A (zh) * | 2020-02-28 | 2020-06-26 | 东华理工大学 | 一种融合自然γ探测器的铀裂变瞬发中子测井控制方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10042082B2 (en) * | 2014-12-04 | 2018-08-07 | Halliburton Energy Services, Inc. | Gamma ray measurement apparatus, systems, and methods |
-
2020
- 2020-11-24 CN CN202011330400.0A patent/CN112523742B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3940610A (en) * | 1974-02-15 | 1976-02-24 | Mobil Oil Corporation | Natural gamma radiation borehole logging system |
| US4542292A (en) * | 1982-04-01 | 1985-09-17 | Schlumberger Technology Corporation | Correction of natural gamma radiation logs for the effects of gamma ray emission from and attenuation by the borehole fluid |
| CN103306663A (zh) * | 2012-03-06 | 2013-09-18 | 中国原子能科学研究院 | 铀矿测井方法 |
| CN107533150A (zh) * | 2015-04-27 | 2018-01-02 | 桑德克斯有线有限公司 | 用于以谱的伽马射线测井确定元素浓度的方法 |
| CN108457640A (zh) * | 2018-01-26 | 2018-08-28 | 东华理工大学 | 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 |
| CN108825219A (zh) * | 2018-05-25 | 2018-11-16 | 东华理工大学 | 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法 |
| CN109557120A (zh) * | 2019-01-14 | 2019-04-02 | 东华理工大学 | 伽马能谱与主动式x荧光联合同时测铀钍镭钾的方法 |
| CN111045072A (zh) * | 2019-12-27 | 2020-04-21 | 核工业北京地质研究院 | 一种适用于CeBr3探测器的γ能谱迭代解谱算法 |
| CN111335882A (zh) * | 2020-02-28 | 2020-06-26 | 东华理工大学 | 一种融合自然γ探测器的铀裂变瞬发中子测井控制方法 |
Non-Patent Citations (3)
| Title |
|---|
| 利用LaBr_3(Ce)伽马谱仪直接测定铀矿体中铀含量的方法;吴永鹏;汤彬;程建平;梁永顺;王哲;;物探与化探;36(03);414-417 * |
| 基于γ谱形差异的定量分析对含铀矿体直接铀定量方法;吴永鹏;汤彬;程建平;王哲;梁永顺;;核电子学与探测技术;32(02);150-155 * |
| 高精度能谱测井与~(234m)Pa特征γ射线的"直接铀定量"技术;汤彬;吴永鹏;张雄杰;王海涛;王仁波;周书民;陈锐;刘庆成;;核技术;35(10);745-750 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112523742A (zh) | 2021-03-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10845502B2 (en) | Optimization of neutron-gamma tools for inelastic gamma-ray logging | |
| US5097123A (en) | Broad energy spectra neutron source for logging and method | |
| CN108457640B (zh) | 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 | |
| US5021653A (en) | Geochemical logging apparatus and method for determining concentrations of formation elements next to a borehole | |
| CN108825219B (zh) | 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法 | |
| CN112904436B (zh) | 一种元素产额与热中子计数比组合的孔隙度测量方法 | |
| CN112523741B (zh) | 一种基于能谱测井交叉谱段的铀矿定量刻度系数求法 | |
| US7112783B2 (en) | Neutron measurement method for determining porosity of a formation surrounding a borehole | |
| US4187908A (en) | Method for ascertaining optimum location for well stimulation and/or perforation | |
| CN112523742B (zh) | 一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法 | |
| US5053620A (en) | Logging apparatus and method for determining concentrations of subsurface formation elements | |
| CN112630848B (zh) | 一种基于能谱测井特征谱段的铀矿定量剥离系数求法 | |
| CN112630849B (zh) | 一种基于能谱测井特征谱峰的铀矿定量剥离系数求法 | |
| CN112649888B (zh) | 一种基于能谱测井特征谱峰的铀矿定量刻度系数求法 | |
| CN112630850A (zh) | 一种基于能谱测井交叉谱段的铀矿定量剥离系数求法 | |
| Killeen | Nuclear techniques for ore grade estimation | |
| Killeen | Borehole logging for uranium by measurement of natural γ-radiation | |
| US2943197A (en) | Method of well logging | |
| CN115032705B (zh) | 一种γ能谱测井铀矿定量五点反褶积解释方法 | |
| CN115614022B (zh) | 铀裂变瞬发中子测井的井眼水层原位修正方法 | |
| Schweitzer et al. | Review of nuclear techniques in subsurface geology | |
| USRE31122E (en) | Method for ascertaining optimum location for well stimulation and/or perforation | |
| RU2624985C1 (ru) | Способ нейтронного каротажа для определения содержания урана в ураново-рудных формациях, пересеченных скважиной | |
| Fan et al. | A New Method of Porosity Determination by DT Neutron Generator and Dual CLYC Detector | |
| Caldwell | Nuclear logging methods |
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 |