CN108457640B - 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 - Google Patents
融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 Download PDFInfo
- Publication number
- CN108457640B CN108457640B CN201810076976.5A CN201810076976A CN108457640B CN 108457640 B CN108457640 B CN 108457640B CN 201810076976 A CN201810076976 A CN 201810076976A CN 108457640 B CN108457640 B CN 108457640B
- Authority
- CN
- China
- Prior art keywords
- uranium
- logging
- ore
- radium
- natural gamma
- 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
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 title claims abstract description 241
- 229910052770 Uranium Inorganic materials 0.000 title claims abstract description 237
- 238000001228 spectrum Methods 0.000 title claims abstract description 53
- 238000004445 quantitative analysis Methods 0.000 title description 2
- 239000011435 rock Substances 0.000 claims abstract description 56
- 238000000034 method Methods 0.000 claims abstract description 53
- 238000011002 quantification Methods 0.000 claims abstract description 42
- 230000004044 response Effects 0.000 claims abstract description 20
- ZKQISJYCQJRICV-UHFFFAOYSA-N [Rn].[Ra].[U] Chemical compound [Rn].[Ra].[U] ZKQISJYCQJRICV-UHFFFAOYSA-N 0.000 claims abstract description 18
- 238000005553 drilling Methods 0.000 claims abstract description 18
- 238000006243 chemical reaction Methods 0.000 claims abstract description 16
- 230000005477 standard model Effects 0.000 claims abstract description 15
- 238000005070 sampling Methods 0.000 claims abstract description 13
- 238000012937 correction Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 230000004927 fusion Effects 0.000 claims abstract description 4
- HCWPIIXVSYCSAN-UHFFFAOYSA-N radium atom Chemical compound [Ra] HCWPIIXVSYCSAN-UHFFFAOYSA-N 0.000 claims description 62
- 229910052705 radium Inorganic materials 0.000 claims description 57
- 229920006395 saturated elastomer Polymers 0.000 claims description 37
- 230000002285 radioactive effect Effects 0.000 claims description 15
- ZSLUVFAKFWKJRC-IGMARMGPSA-N 232Th Chemical compound [232Th] ZSLUVFAKFWKJRC-IGMARMGPSA-N 0.000 claims description 14
- 229910052776 Thorium Inorganic materials 0.000 claims description 14
- 230000014509 gene expression Effects 0.000 claims description 13
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 13
- 239000011707 mineral Substances 0.000 claims description 13
- 239000000126 substance Substances 0.000 claims description 10
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 claims description 7
- 239000011591 potassium Substances 0.000 claims description 7
- 229910052700 potassium Inorganic materials 0.000 claims description 7
- 238000004519 manufacturing process Methods 0.000 claims description 6
- 150000001875 compounds Chemical class 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 150000001224 Uranium Chemical class 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000009826 distribution Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 claims description 2
- 230000003993 interaction Effects 0.000 claims description 2
- 238000005316 response function Methods 0.000 claims description 2
- 238000011065 in-situ storage Methods 0.000 abstract description 8
- 238000002386 leaching Methods 0.000 abstract description 8
- 238000005065 mining Methods 0.000 abstract description 8
- 238000004458 analytical method Methods 0.000 abstract description 7
- 238000009614 chemical analysis method Methods 0.000 abstract description 5
- 230000004992 fission Effects 0.000 description 12
- 230000005251 gamma ray Effects 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 239000012530 fluid Substances 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 229910001220 stainless steel Inorganic materials 0.000 description 2
- 239000010935 stainless steel Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect 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
- 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
Landscapes
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (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)
- Analysing Materials By The Use Of Radiation (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种用于铀矿勘探和地浸采铀领域的融合瞬发中子时间谱测井修正自然γ总量测井的铀矿定量方法。具体是指:针对钻孔中的含铀岩矿层,利用自然γ总量测井和瞬发中子时间谱测井的两类测井数据,构建出基于自然γ总量测井且有“铀‑镭‑氡”平衡修正的铀矿定量解释方法;并利用标准模型井的测井数据,构建出铀矿定量所需的换算系数、本底响应等刻度参数求取方法。相比单一瞬发中子时间谱测井或自然γ总量测井的铀矿定量方法,本发明无需岩芯取样和化学分析方法,既能求取“铀‑镭‑氡”平衡系数,还能大大提高测井速度,进而具有提高钻探效率、降低勘探成本、缩短铀定量周期等优点,并能通过计算机编程实现铀矿定量的现场解析。
Description
技术领域
本发明涉及用于铀矿勘探和地浸采铀领域的核测井铀矿定量方法,尤其是一种融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法。
背景技术
固体矿产定量主要依赖于钻孔岩芯取样的化学分析方法,该法是最直接的岩矿层元素定量的传统方法,但也带来了取芯钻探效率低、岩芯取样成本高、化学分析周期长等缺点。然而,我国铀矿普遍具有含量偏低、矿厚偏薄与矿体分散等特点,例如砂岩型铀矿的最低可采品位仅为0.01%铀含量,可采厚度仅需70cm以上,因此岩芯取样的化学分析方法更是加剧了上述缺点。
地壳中的铀系、钍系和钾元素贡献了约占99%的自然γ射线,因此铀矿勘探和地浸采铀几乎离不开自然γ总量或γ能谱测井。通常以自然γ总量测井的定量解释方法为主,岩芯取样的化学分析方法为辅,共同实现铀矿定量。鉴于我国铀、钍、钾等放射性混合的矿产不足5%,仅需自然γ总量测井就能初步确定矿体的铀或镭的含量。
铀、钍作为铀系、钍系的母体元素,母体自身几乎不放射γ射线,而是某些衰变子体放射γ射线,称为γ核素。当母体与子体的衰变率相同时,称其达到放射性平衡,处于放射性平衡的铀系、钍系就能按各自γ核素放射的γ射线推断出铀、钍及子体含量。通常,钍系极易达到放射性平衡,按钍系γ射线确定钍含量并无悬念;然而铀系中的“铀-镭-氡”平衡很难保证,但将铀系作为铀组、镭组两个子系看待,子系放射性平衡更易保证。又因镭组占整个铀系γ射线的份额高达97%以上,依据铀系γ射线确定镭含量并非难事,但依此确定铀含量的风险很大。也就是,在镭组γ射线绝对占优的前提下,难以确保测井解释的铀含量具有可靠性,因此我国《铀矿地质勘查规范》要求不低于30%的岩芯化学分析求取“铀-镭-氡”平衡系数,并对测井解释的铀含量进行修正,才能求得真实的铀含量。可见,该法难以彻底规避取芯钻探效率低、岩芯取样成本高、化学分析周期长等缺点。
从理论上说,脉冲中子引发铀裂变的瞬发中子测井方法不受铀以外的放射性元素干扰,是用于铀矿定量的理想方法。但相比自然γ总量测井,铀裂变的瞬发中子数很少,使得瞬发中子测井的速度很慢,其测井效率甚至低到难以实用。为此,将两种测井方法有机融合并形成新测井方法,即使钻孔中单个测点的瞬发中子测井计数率很低,但它在整个岩矿层中的计数率之和并不低,只要采用自然γ总量测井求取岩矿层各测点的镭或镭组元素含量,再利用瞬发中子测井的计数率之和求取整个岩矿层的平均铀含量,进而求得矿体的“铀-镭-氡”平衡系数,就能将各测点的镭含量修正到铀含量,因而无须岩芯取样和化学分析,也能快速、准确地确定各测点的真实铀含量。
可见,本发明公开的融合瞬发中子时间谱测井修正自然γ总量测井的铀矿定量方法适用于铀矿勘探和地浸采铀,特别是在地浸采铀领域,当铀不断从岩矿层中采走之后,镭及γ核素仍然留在原地,有效区分镭和剩余铀的测井方法和刻度参数求取方法已成为亟待攻克的科技难题。因此,本发明具有重要的科学意义和实用价值。
发明内容
本发明的目的在于研发出一种用于铀矿勘探和地浸采铀领域的融合瞬发中子时间谱测井修正自然γ总量测井的铀矿定量方法。
本发明的技术方案为:一种融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法,其特征在于:该方法以自然γ总量测井为主、瞬发中子时间谱测井为辅,融合这两类测井数据,在钻孔中求得含铀岩矿层的“铀-镭-氡”平衡系数,并修正自然γ总量测井的铀矿定量解释结果,依此求得各个测点的真实铀含量;以及利用饱和矿层构建的标准模型井及测井数据,依此建立铀矿定量解释所需的换算系数、本底响应等刻度参数的求取方法。
所述融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法,鉴于自然γ射线很难对3He正比计数管等中子探测器产生计数率贡献,但中子能激发钻孔及岩矿层物质产生非弹和俘获等γ射线,如果γ探测器和中子探测器之距不足1.5m,非弹和俘获等γ射线就将对γ探测器产生计数率贡献;因此,瞬发中子测井仪附带的自然γ测井功能难以独立实现铀矿定量,以自然γ测井解释结果实现最终的铀矿定量,并借助瞬发中子时间谱测井来修正自然γ总量测井的铀矿定量结果,就能大大提高测井速度,并实现低成本的快速铀矿定量;
由于瞬发中子时间谱测井的计数率仅与含铀岩矿层的铀(235U)含量存在正比关系;对于忽略钍系和钾,甚至忽略铀系铀组γ射线的铀矿定量而言,自然γ总量的测井计数率主要与岩矿层的镭子体含量存在正比关系;如果沿钻孔方向将含铀岩矿层分解为无数个铀矿薄层,并将任一铀矿薄层的铀或镭含量记为qk(z)函数,与之对应的自然γ总量与瞬发中子时间谱的测井计数率记为Nk(Z)函数,由此构建的qk(z)与Nk(Z)之间的函数关系为:
上式称为铀矿定量正演方程的积分表达式,与之相应的微分表达式为:
式中,
Z、z分别表示钻孔中任一测点、任一铀矿薄层的深度坐标,且|Z-z|为任一测点Z到任一铀矿薄层dz的距离;
H1、H2表示当沿钻孔方向以z点将含铀岩矿层分为上下两部分时,该矿层上部和下部各自所占的厚度,即H=H1+H2为含铀岩矿层的总厚度;
qk(z)表示深度坐标z处的铀矿薄层所含铀或镭含量,由下标k(=1,2)来区分铀或镭,即q1(z)为铀含量,q2(z)为镭含量;
Nk(Z)表示含铀岩矿层在深度坐标Z处形成的测井计数率,由下标k(=1,2)来区分与铀定量或镭定量对应的测井计数率,即N1(Z)为瞬发中子时间谱测井的计数率,N2(Z)为自然γ总量测井的计数率;此外,任一铀矿薄层在深度坐标Z处形成的测井计数率采用微分元dNk(Z)表示;
称为地质脉冲响应函数近似表达式;其中用于表达任一铀矿薄层在任一测点所响应的测井计数率,与两者之距|Z-z|密切相关,近似为负指数衰减,且衰减速率称为特征参数αk,取值与射线种类、射线与岩矿层的相互作用等因素相关;
Ak称为换算系数,表示单位含量的铀或镭元素在饱和矿层中心部位所能形成的测井计数率;其中,饱和矿层是指“无限厚”均匀矿层,实指厚度H中的H1和H2至少均达0.6m以上,此时的矿层中心部位基本满足饱和矿层定义;应当注意,镭含量通常采用平衡铀含量表示,对于平衡铀系构建的饱和矿层,铀与镭含量在数值上必定相等,即q1(z)=q2(z),但两种的计数率并不相等;可见,Ak是一个与岩矿层物质成分和结构相关的参量;
Bk称为本底响应;包括测井仪、钻孔与岩矿层等环境因素产生的本底计数率,当环境因素确定后,Bk取值为常数;通常,自然γ总量的本底计数率包含非弹与俘获等γ射线产生的测井计数率,特别是中子与γ探测器之距不足1.5m时,中子激发岩矿层和钻孔物质产生非弹与俘获γ射线的效应尤为明显;上述(1)式和(2)式分别是定量解释铀、镭含量沿钻孔方向分布的积分表达式和微分表达式,统称为铀矿定量正演方程,且积分表达式最为常用。
所述融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法,由铀矿定量正演方程求解岩矿层内任一点铀、镭含量的前提是:必须事先确定换算系数Ak、本底响应Bk等刻度参数;
选取一个饱和矿层的模型井,因该饱和矿层内任一点的铀、镭含量qk(z)处处相等,则可将深度坐标z处的铀、镭含量作为常数看待,并采用该岩矿层中心点的含量表示该常数,即qk(z)=qk(Z0),其中Z0为该岩矿层中心点的深度坐标;由此求得深度坐标为Z的任一测点的测井计数率为:
由此可知,饱和矿层任一测点的测井计数率也处处相等,也可采用饱和矿层中心点表示,即Nk(Z)=Nk(Z0),由此求得换算系数Ak的计算公式为:
采用不含铀、镭等任何放射性元素的饱和矿层求得本底响应Bk,也就是选择零值饱和矿层的模型井,即按qk(Z0)=0代入(3)式,由此求得:
Bk=Nk(Z0)(其中k=1,2) (4)在野外生产测井及其铀矿定量时,为了选取合适的换算系数Ak和本底响应Bk,通常必须建造出岩矿层物质成分、钻孔结构及测井条件与野外生产井相同或相近的饱和矿层刻度模型井,这样的刻度模型井称为标准模型井,简称标准井,包括铀系平衡或不平衡的标准铀模型井、标准零值模型井等等,并至少分为砂岩型标准模型井、硬岩型标准模型井两大系列。
所述融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法,在确定合适的换算系数Ak和本底响应Bk等刻度参数之后,按照铀矿定量正演方程的积分表达式,即(1)式,就可求解出含铀岩矿层内任一点的铀或镭含量qk(z);当镭含量采用平衡铀含量表示时,“铀-镭-氡”平衡系数KP可按单个含铀岩矿层甚至整个铀矿区的铀矿定量结果求得,即:
式中,
Zj是对钻孔中连续测量点进行离散取样后的测点深度坐标,其中取样间距满足Zj=Zj-1+ΔZ,ΔZ为等间距取样,铀矿勘探一般选择ΔZ=10或20cm;
实验表明,瞬发中子时间谱测井与自然γ总量测井的计数率几乎相差2个量级,即瞬发中子时间谱测井的计数率很低,统计涨落误差很大,使得该法的最高测井速度仅被建议为0.3~0.6m/min,与我国《铀矿地质勘查规范》和《铀矿γ测井规范》要求的2~4m/min的测井速度有很大差距;为提高测井速度并降低统计涨落误差,将自然γ总量测井作为铀矿定量解释的主要测井方法,在定量解释铀、镭含量qk(z)之前,直接采用上式对单个含铀岩矿层或整个铀矿区段的测井计数率进行求和,并按其比值求得“铀-镭-氡”平衡系数KP,然后对自然γ总量测井的定量解释结果进行“铀-镭-氡”平衡修正,由此求得各测点的真实铀含量必定具有较小的统计涨落误差;
应当注意,通常还需将铀矿定量正演方程(1)式或(2)式转化为相应的铀矿定量反演方程,也就是《铀矿γ测井规范》要求的“平均含量法”、“五点式反褶积法”等铀矿定量解释方法,并根据野外生产井的钻孔井径大小、井液成分、套管材料与厚度,及其岩矿成分等影响因素,对采用标准模型井求得的换算系数Ak和本底响应Bk等刻度参数进行影响因素修正,依此由自然γ总量测井的铀矿定量解释结果及“铀-镭-氡”平衡修正才能求得矿体的真实铀含量,并最终用于铀矿储量估算。
本发明的优点在于:相比单一瞬发中子时间谱测井或自然γ总量测井的铀矿定量方法,本发明无需岩芯取样和化学分析方法,既能求取“铀-镭-氡”平衡系数,还能大大提高测井速度,进而具有提高钻探效率、降低勘探成本、缩短铀定量周期等优点,并能通过计算机编程实现铀矿定量的现场解析。
附图说明
图1、自然γ总量测井与瞬发中子时间谱测井的组合测井仪基本结构;
图2(a)、实测标准饱和模型Nb-4(本底模型)的自然γ总量曲线;
图2(b)、实测标准饱和模型Nb-4(本底模型)的瞬发中子时间谱曲线;
图3(a)、实测标准饱和模型Nu-3(铀含量为0.1%)的自然γ总量曲线;
图3(b)、实测标准饱和模型Nu-3(铀含量为0.1%)的瞬发中子时间谱曲线;
图4(a)、实测标准饱和模型Nu-1(铀含量为0.03%)的自然γ总量曲线;
图4(b)、实测标准饱和模型Nu-1(铀含量为0.03%)的瞬发中子时间谱曲线;
图5(a)、实测标准饱和模型Nu-2(铀含量为0.07%)的自然γ总量曲线;
图5(b)、实测标准饱和模型Nu-2(铀含量为0.07%)的瞬发中子时间谱曲线;
图6、利用饱和矿层刻度模型井求取关键测井参数的步骤流程图。
图中:脉冲中子发生器1、双中子时间谱测量系统2、自然γ总量测量系统3、马笼头4、不锈钢探管外壳5、测井电缆6、矿层或岩层7。
具体实施方式
本发明公开了用于铀矿探采(铀矿勘探和地浸采铀)领域的自然γ总量测井与铀裂变瞬发中子时间谱测井融合在一起的组合测井仪基本结构,以及利用该组合测井确定岩矿层铀、镭、钍、钾等天然放射性元素含量的刻度参数求法。包括该组合测井仪基本结构(其框图如图1所示)、实测砂岩饱和矿层刻度模型井获得的自然γ总量测井数据(饱和矿层中心点的自然γ总量曲线示例结果如图2(a)、图3(a)、图4(a)、图5(a)所示)、铀裂变瞬发中子时间谱测井数据(饱和矿层中心点的超热中子与热中子时间谱曲线示例如图2(b)、图3(b)、图4(b)、图5(b)所示)、求取该组合测井仪换算系数、本底响应等关键测井参数的刻度公式,以及刻度数据(如表1)等内容。本发明公开的组合测井仪基本结构与刻度参数求法的示例及说明如下:
表1、标准饱和模型井的测井数据及其刻度参数求取结果:
表1.标准饱和模型井的测井数据及其刻度参数求取结果注2
表2、标准饱和模型井的测井数据及其解释结果:
表2.标准饱和模型井的测井数据及其解释结果注2
注1:镭含量用平衡时的铀含量表示;
注2:上表数据为2017年在核工业放射性勘查计量站(石家庄)标准模型井上测得的计数率,以及求得的解释结果;
注3:铀对应的特征峰表示经热中子修正后的超热中子计数率,中子时间谱范围为128~800us,镭对应的特征能区为0.40~3.00MeV。
1、组合测井仪基本结构与特征能区计数率
因铀系、钍系和钾元素含有的γ核素并不多,然而它们放出上百种能量的特征γ射线。现有自然γ总量测井仅仅能够分辨出为数不多的十几种特征γ射线,即放射几率>0.01、能量>0.4MeV且无重叠峰的特征γ射线。其中,铀组放出的特征γ射线因其放射几率和能量都太低,因而几乎无法分辨出铀的特征峰,也很难采用单一的自然γ总量测井方法确定天然岩矿层的铀含量。因而本发明将自然γ总量测井与铀裂变瞬发中子时间谱测井融合在一起,设计出了一种组合测井仪,它的基本结构如图1所示。从图1可知,该组合测井仪由马笼头4、不锈钢探管外壳5及探管内的脉冲中子发生器1、双中子时间谱测量系统2、自然γ总量测量系统3组成,通过测井电缆6控制脉冲中子发生器、读取自然γ总量数据和瞬发中子时间谱数据;
根据该类组合测井仪的基本结构及其对放射性元素的定量能力,以及该组合测井仪在砂岩饱和矿层刻度模型井中实测的示例数据,如图2、图3、图4、图5所示,选定放射性元素定量的特征能区计数率如下:
对于铀定量,选定铀裂变瞬发中子时间谱测井的铀裂变时间段在128us~800us计超热中子与累计热中子的比值为铀特征计数率N1(Z),它在饱和矿层刻度模型井中实测的示例数据分别如图2(b)、图3(b)、图4(b)和图5(b)所示。
对于镭定量,选择γ射线能量范围在0.4MeV~3.0MeV的γ总量计数,依此确定镭特征计数率N2(Z),它在饱和矿层刻度模型井中实测的示例数据分别如图2(a)、图3(a)、图4(a)和图5(a)所示。
应当注意,本发明公开的组合测井仪可以同步开展自然γ总量测井与铀裂变瞬发中子时间谱测井,但要以自然γ总量测井的探测器中心点为基准,并根据铀裂变瞬发中子时间谱测井的探测器中心点的偏差大小,将两者的测井数据调整为同一测点的测井数据。
2、组合测井仪用于铀矿定量时的刻度参数求法与示例
在开展野外生产测井之前,需在饱和矿层刻度模型井中对该组合测井仪进行刻度,并依此求得该组合测井仪的换算系数、本底响应等关键测井参数。利用饱和矿层刻度模型井求取关键测井参数的步骤流程如图6所示。
也就是在饱和矿层刻度模型井中部(深度坐标记为Z0),由铀裂变瞬发中子时间谱测井的累计超热中子与累计热中子比值作为铀特征计数率N1(Z0);由自然γ总量测井的特征计数率作为镭特征计数率N2(Z0)。
按照饱和矿层刻度模型井的放射性元素含量qk(Z0)为已知,将其代入(4)式和(5)式,依此求得该组合测井仪的换算系数、本底响应等关键测井参数,其结果示例如表2所示。
应当注意,在铀矿测井中,主要由镭元素为主的自然γ射线对中子探测器基本没有产生计数率贡献,则A12=0;因中子激发钻孔物质产生非弹和俘获等γ射线,比自然γ总量的测井计数率要低1~2个量级,因而在一定程度上也可以认为,以镭元素为主的自然γ射线对铀特征计数率没有贡献,则A21=0。
可见,对自然γ总量测井与铀裂变瞬发中子时间谱测井的测井曲线融合,就能找到对应于铀、镭等放射性元素(对应元素编号k=1,2)定量的特征峰或等效的瞬发中子时间谱测井数据(对应特征γ射线编号i=1,2)。
还应注意,如果刻度模型井与野外钻孔的情况相差较大,还应当制作各类修正量版,包括孔径修正量版、井液修正量版、套管修正量版等,以便将换算系数Aki与本底相应Bi修正为野外钻孔的实际情况中。
在本发明建立的瞬发中子时间谱与自然γ总量的组合测井仪的刻度参数求法中,虽然涉及的公式较多,具体实施也较复杂,但可以将它们编制为计算机软件,以此形成铀矿定量解释系统。
为对上述方法进行验证,根据图6所示步骤,利用两个饱和矿层刻度模型井求取了关键测井参数,并通过组合测井仪测量了两个不平衡饱和砂岩模型对算法进行了验证,测量得到的自然伽马能谱和瞬发中子时间谱图如图2、图3、图4、图5所示,模型各特征计数率及刻度参数求解结果见表1、各特征计数率及解释结果见表2。
从表1、表2可以看出,对于钍、钾含量很低的砂岩型铀矿(标准模型井),通过本文构建融合了自然γ总量测井和瞬发中子时间谱测井的组合测井方法,任选一个饱和砂岩铀标准模型井和一个零值标准模型井刻度出来的换算系数Ak和本底响应Bk等刻度参数,并通过自然γ总量测井和瞬发中子时间谱测井的组合测井方法,能够求准同类铀矿(例如砂岩型铀矿)的铀、镭含量,其含量绝对误差≤±0.002%eU,相对误差≤±2.5%,完全能够满足《铀矿地质勘查规范》的要求。
表1中,镭、铀定量的换算系数分别为A2=271.36(s-1/10-4eU),A1=3.55(s-1/10- 4eU)。即0.01%铀含量的饱和矿层产生的γ总计数率为271.36s-1,超热中子计数率为3.35s-1,两者几乎相差2个量级。表明瞬发中子时间谱测井的计数率太小,为降低统计涨落误差,使得最高测井速度仅能达到0.3m/min左右,与《铀矿地质勘查规范》要求的2-4m/min的测井速度有很大差距。本发明提出的在整个矿段内分别按两种测井方法求得各自的测井计数率之和(该累计值的统计涨落误差必定很小),再求得该矿段内的平均镭、铀含量,按平均镭、铀含量之比求得该矿段的铀-镭-氡平衡系数;然后,采用统计涨落误差很小的自然γ总量测井定量解释各测点的镭含量,最后通过铀-镭-氡平衡修正求出各测点的真实铀含量。
根据以上文本,本发明的优点在于:利用该组合测井仪确定岩矿层铀、镭等天然放射性元素含量的刻度参数求法,便可在实际测井过程中无需借助岩芯取样和化学分析的基础上,直接求出放射性矿层或岩矿层中铀、镭等放射性元素含量,以及铀-镭-氡平衡系数,对于铀矿勘探和地浸采铀领域的铀矿储量估算,具有提高钻探效率、降低勘探成本、缩短铀定量周期等技术优势。
Claims (4)
1.一种融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法,其特征在于:该方法以自然γ总量测井为主、瞬发中子时间谱测井为辅,融合这两类测井数据:在钻孔中求得含铀岩矿层的“铀-镭-氡”平衡系数,并修正自然γ总量测井的铀矿定量解释结果,依此求得各个测点的真实铀含量;以及利用饱和矿层构建的标准模型井及测井数据,依此建立铀矿定量解释所需的换算系数、本底响应刻度参数的求取方法;具体为:将自然γ总量测井作为铀矿定量解释的主要测井方法,因为单个矿层或整个铀矿区段的“铀-镭-氡”平衡系数不会有急剧变化,在定量解释铀、镭含量之前,求得单个含铀岩矿层或整个铀矿区段的分别代表铀、镭含量的测井计数率之和,并按其比值求得单个含铀岩矿层或整个铀矿区段的平均“铀-镭-氡”平衡系数,然后对各测点由自然γ总量测井求得的镭含量进行“铀-镭-氡”平衡修正,得到各测点的真实铀含量。
2.根据权利要求1所述融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法,其特征在于:
鉴于自然γ射线很难对3He正比计数管中子探测器产生计数率贡献,但中子能激发钻孔及岩矿层物质产生非弹和俘获γ射线,如果γ探测器和中子探测器之距不足1.5m,非弹和俘获γ射线就将对γ探测器产生计数率贡献;因此,瞬发中子测井仪附带的自然γ测井功能难以独立实现铀矿定量,以自然γ测井解释结果实现最终的铀矿定量,并借助瞬发中子时间谱测井来修正自然γ总量测井的铀矿定量结果,就能大大提高测井速度,并实现低成本的快速铀矿定量;
由于瞬发中子时间谱测井的计数率仅与含铀岩矿层的铀含量存在正比关系;对于忽略钍系和钾,甚至忽略铀系铀组γ射线的铀矿定量而言,自然γ总量的测井计数率主要与岩矿层的镭子体含量存在正比关系;如果沿钻孔方向将含铀岩矿层分解为无数个铀矿薄层,并将任一铀矿薄层的铀或镭含量记为qk(z)函数,与之对应的自然γ总量与瞬发中子时间谱的测井计数率记为Nk(Z)函数,由此构建的qk(z)与Nk(Z)之间的函数关系为:
上式称为铀矿定量正演方程的积分表达式,与之相应的微分表达式为:
式中,
Z、z分别表示钻孔中任一测点、任一铀矿薄层的深度坐标,且|Z-z|为任一测点Z到任一铀矿薄层dz的距离;
H1、H2表示当沿钻孔方向以z点将含铀岩矿层分为上下两部分时,该矿层上部和下部各自所占的厚度,即H=H1+H2为含铀岩矿层的总厚度;
qk(z)表示深度坐标z处的铀矿薄层所含铀或镭含量,由下标k来区分铀或镭,即q1(z)为铀含量,q2(z)为镭含量;
Nk(Z)表示含铀岩矿层在深度坐标Z处形成的测井计数率,由下标k来区分与铀定量或镭定量对应的测井计数率,即N1(Z)为瞬发中子时间谱测井的计数率,N2(Z)为自然γ总量测井的计数率;此外,任一铀矿薄层在深度坐标Z处形成的测井计数率采用微分元dNk(Z)表示;
称为地质脉冲响应函数近似表达式;其中用于表达任一铀矿薄层在任一测点所响应的测井计数率,与两者之距|Z-z|密切相关,近似为负指数衰减,且衰减速率称为特征参数αk,取值与射线种类、射线与岩矿层的相互作用因素相关;
Ak称为换算系数,表示单位含量的铀或镭元素在饱和矿层中心部位所能形成的测井计数率;其中,饱和矿层是指“无限厚”均匀矿层,实指厚度H中的H1和H2至少均达0.6m以上,此时的矿层中心部位基本满足饱和矿层定义;应当注意,镭含量采用平衡铀含量表示,对于平衡铀系构建的饱和矿层,铀与镭含量在数值上必定相等,即q1(z)=q2(z),但两种的计数率并不相等;可见,Ak是一个与岩矿层物质成分和结构相关的参量;
Bk称为本底响应;包括测井仪、钻孔与岩矿层环境因素产生的本底计数率,当环境因素确定后,Bk取值为常数;自然γ总量的本底计数率包含非弹与俘获γ射线产生的测井计数率,特别是中子与γ探测器之距不足1.5m时,中子激发岩矿层和钻孔物质产生非弹与俘获γ射线的效应尤为明显;
上述(1)式和(2)式分别是定量解释铀、镭含量沿钻孔方向分布的积分表达式和微分表达式,统称为铀矿定量正演方程。
3.根据权利要求2所述融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法,其特征在于:
由铀矿定量正演方程求解岩矿层内任一点铀、镭含量的前提是:必须事先确定换算系数Ak、本底响应Bk刻度参数;
选取一个饱和矿层的模型井,因该饱和矿层内任一点的铀、镭含量qk(z)处处相等,则可将深度坐标z处的铀、镭含量作为常数看待,并采用该岩矿层中心点的含量表示该常数,即qk(z)=qk(Z0),其中Z0为该岩矿层中心点的深度坐标;由此求得深度坐标为Z的任一测点的测井计数率为:
由此可知,饱和矿层任一测点的测井计数率也处处相等,采用饱和矿层中心点表示,即Nk(Z)=Nk(Z0),由此求得换算系数Ak的计算公式为:
采用不含铀、镭任何放射性元素的饱和矿层求得本底响应Bk,也就是选择零值饱和矿层的模型井,即按qk(Z0)=0代入(3)式,由此求得:
Bk=Nk(Z0)其中k=1,2 (4)
在野外生产测井及其铀矿定量时,为了选取合适的换算系数Ak和本底响应Bk,必须建造出岩矿层物质成分、钻孔结构及测井条件与野外生产井相同或相近的饱和矿层刻度模型井,这样的刻度模型井称为标准模型井,简称标准井,包括铀系平衡或不平衡的标准铀模型井、标准零值模型井,并至少分为砂岩型标准模型井、硬岩型标准模型井两大系列。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810076976.5A CN108457640B (zh) | 2018-01-26 | 2018-01-26 | 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810076976.5A CN108457640B (zh) | 2018-01-26 | 2018-01-26 | 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108457640A CN108457640A (zh) | 2018-08-28 |
CN108457640B true CN108457640B (zh) | 2021-08-06 |
Family
ID=63239027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810076976.5A Active CN108457640B (zh) | 2018-01-26 | 2018-01-26 | 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108457640B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108825219B (zh) * | 2018-05-25 | 2021-08-10 | 东华理工大学 | 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法 |
CN111413741B (zh) * | 2019-01-04 | 2022-12-02 | 中国石油天然气股份有限公司 | 一种砂岩型铀矿资源量计算方法和装置 |
CN111335882B (zh) * | 2020-02-28 | 2022-11-01 | 东华理工大学 | 一种融合自然γ探测器的铀裂变瞬发中子测井控制方法 |
CN112523742B (zh) * | 2020-11-24 | 2023-11-17 | 东华理工大学 | 一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法 |
CN112630849B (zh) * | 2020-11-24 | 2023-10-10 | 东华理工大学 | 一种基于能谱测井特征谱峰的铀矿定量剥离系数求法 |
CN112523741B (zh) * | 2020-11-24 | 2023-04-14 | 东华理工大学 | 一种基于能谱测井交叉谱段的铀矿定量刻度系数求法 |
CN113671556B (zh) * | 2021-07-30 | 2024-02-20 | 中国原子能科学研究院 | 一种铀裂变电离室铀靶铀同位素原子核数量定值方法 |
CN115032705B (zh) * | 2022-01-18 | 2023-02-17 | 东华理工大学 | 一种γ能谱测井铀矿定量五点反褶积解释方法 |
CN115614022B (zh) * | 2022-11-10 | 2023-04-28 | 东华理工大学 | 铀裂变瞬发中子测井的井眼水层原位修正方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3845299A (en) * | 1972-11-29 | 1974-10-29 | Mobil Oil Corp | Method for uranium exploration employing radioactive reconnaissance and assay logging |
US4201912A (en) * | 1978-10-03 | 1980-05-06 | The United States Of America As Represented By The United States Department Of Energy | Subthreshold neutron interrogator for detection of radioactive materials |
US4326129A (en) * | 1978-05-10 | 1982-04-20 | Jacob Neufeld | Evaluation of the interaction with radiant energy of substances traversed by a bore hole |
US4378498A (en) * | 1980-10-07 | 1983-03-29 | Mobil Oil Corporation | Simulation of subsurface formations for the testing and calibration of radioactive well logging tools |
US20060192096A1 (en) * | 2005-02-28 | 2006-08-31 | Schlumberger Technology Corporation | Apparatus and methods for interlaced density and neutron measurements |
CN101078776A (zh) * | 2007-06-29 | 2007-11-28 | 西安奥华电子仪器有限责任公司 | 核测井数据多尺度处理方法 |
CN101078775A (zh) * | 2007-06-29 | 2007-11-28 | 西安奥华电子仪器有限责任公司 | 脉冲中子双谱饱和度测井方法 |
CN102175704A (zh) * | 2011-02-18 | 2011-09-07 | 中国原子能科学研究院 | 铀同位素丰度分析方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA1095635A (en) * | 1977-01-17 | 1981-02-10 | Wyatt W. Givens | Prompt fission neutron logging for uranium |
US4698499A (en) * | 1983-12-19 | 1987-10-06 | Halliburton Company | Quantitative evaluation of uranium ore zones |
CN1115571C (zh) * | 1999-08-18 | 2003-07-23 | 上海伽玛射线技术有限公司 | 放射性比活度测定方法及装置 |
CN104074501B (zh) * | 2013-03-29 | 2016-08-10 | 核工业航测遥感中心 | 一种裂变中子测井修正方法 |
CN103470252B (zh) * | 2013-09-13 | 2017-01-04 | 东华理工大学 | 基于超热中子时间谱的瞬发中子测井及铀矿定量方法 |
US20150226875A1 (en) * | 2014-02-12 | 2015-08-13 | Ge Oil & Gas Logging Services, Inc. | Single Sensor for Detecting Neutrons and Gamma Rays |
RU2579822C1 (ru) * | 2015-02-05 | 2016-04-10 | Открытое акционерное общество Опытно-демонстрационный центр вывода из эксплуатации уран-графитовых ядерных реакторов | Способ контроля стабильности внутренних барьеров безопасности в пункте консервации уран-графитового реактора |
-
2018
- 2018-01-26 CN CN201810076976.5A patent/CN108457640B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3845299A (en) * | 1972-11-29 | 1974-10-29 | Mobil Oil Corp | Method for uranium exploration employing radioactive reconnaissance and assay logging |
US4326129A (en) * | 1978-05-10 | 1982-04-20 | Jacob Neufeld | Evaluation of the interaction with radiant energy of substances traversed by a bore hole |
US4201912A (en) * | 1978-10-03 | 1980-05-06 | The United States Of America As Represented By The United States Department Of Energy | Subthreshold neutron interrogator for detection of radioactive materials |
US4378498A (en) * | 1980-10-07 | 1983-03-29 | Mobil Oil Corporation | Simulation of subsurface formations for the testing and calibration of radioactive well logging tools |
US20060192096A1 (en) * | 2005-02-28 | 2006-08-31 | Schlumberger Technology Corporation | Apparatus and methods for interlaced density and neutron measurements |
CN101078776A (zh) * | 2007-06-29 | 2007-11-28 | 西安奥华电子仪器有限责任公司 | 核测井数据多尺度处理方法 |
CN101078775A (zh) * | 2007-06-29 | 2007-11-28 | 西安奥华电子仪器有限责任公司 | 脉冲中子双谱饱和度测井方法 |
CN102175704A (zh) * | 2011-02-18 | 2011-09-07 | 中国原子能科学研究院 | 铀同位素丰度分析方法 |
Non-Patent Citations (1)
Title |
---|
钻孔γ场理论与核测井分层解释方法研究与应用;汤彬;《中国博士学位论文全文数据库.基础科学辑》;20090315(第3期);第41、46-47、61、66-67、71、107-111页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108457640A (zh) | 2018-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108457640B (zh) | 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 | |
US10845502B2 (en) | Optimization of neutron-gamma tools for inelastic gamma-ray logging | |
CN108825220B (zh) | 融合自然γ能谱与中子时间谱的铀矿测井仪及铀定量方程 | |
US5097123A (en) | Broad energy spectra neutron source for logging and method | |
CN108825219B (zh) | 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法 | |
CN105874354A (zh) | 用于进行井下测量的系统和方法 | |
EP0473467A2 (en) | Gamma-ray spectroscopy method and apparatus for determining concentration of elements in an unknown material | |
US20140197324A1 (en) | Methods to enhance nuclear spectroscopy analysis | |
CN112987117B (zh) | 基于自然γ能谱测井多特征峰组合的铀矿定量换算系数求法 | |
CN112523741B (zh) | 一种基于能谱测井交叉谱段的铀矿定量刻度系数求法 | |
CN112630849B (zh) | 一种基于能谱测井特征谱峰的铀矿定量剥离系数求法 | |
CN112523742B (zh) | 一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法 | |
CN111859675B (zh) | 一种基于铀裂变瞬发中子测井数据泥浆密度修正方法 | |
CN112630848B (zh) | 一种基于能谱测井特征谱段的铀矿定量剥离系数求法 | |
Moake | Characterizing natural-gamma-ray tools without the API calibration formation | |
CN112649888B (zh) | 一种基于能谱测井特征谱峰的铀矿定量刻度系数求法 | |
CN115614022B (zh) | 铀裂变瞬发中子测井的井眼水层原位修正方法 | |
Keys et al. | Well logging with californium-252 | |
Zhang | Development of a surrogate model for elemental analysis using a natural gamma ray spectroscopy tool | |
Schweitzer et al. | Review of nuclear techniques in subsurface geology | |
Aksoy et al. | Elemental analysis using natural gamma-ray spectroscopy | |
Martynenko et al. | The Prospects for the Application of the New Generation Neutron Logging System Developed by Rusburmash Inc.(Russian Federation) for Handling Geotechnical Issues at Sandstone-Hosted Uranium Deposits | |
Yang et al. | Design and performance of an epithermal neutron detector based on PFN | |
Tanner et al. | Assay for uranium and determination of disequilibrium by means of in situ high resolution gamma-ray spectrometry | |
CN112630850A (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 |