CN106556856B - 长期密封矿样氡射气系数快速测量的方法 - Google Patents
长期密封矿样氡射气系数快速测量的方法 Download PDFInfo
- Publication number
- CN106556856B CN106556856B CN201611094699.8A CN201611094699A CN106556856B CN 106556856 B CN106556856 B CN 106556856B CN 201611094699 A CN201611094699 A CN 201611094699A CN 106556856 B CN106556856 B CN 106556856B
- Authority
- CN
- China
- Prior art keywords
- radon
- sample
- atomicity
- formula
- coefficient
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/167—Measuring radioactive content of objects, e.g. contamination
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measurement Of Radiation (AREA)
Abstract
长期密封矿样氡射气系数快速测量的方法,首先用高纯锗γ谱仪测得长期密封矿样内平衡状态的氡原子数N1,然后拆开密封样品,并在开放环境中以搅拌的方式排出样品中的自由氡,再重新密封后测得样品内未平衡状态氡原子数N2,最后利用N1、N2及建立的长期密封矿样氡射气系数理论计算公式计算出样品的氡射气系数η。所述的长期密封矿样氡射气系数理论计算公式如下:其中:
Description
技术领域
本发明涉及辐射测量技术领域,特别是一种对长期密封矿样的氡射气系数进行快速测量的方法。
背景技术
铀资源是核行业发展的基础,而铀矿的安全开发和利用更是核安全的第一环节和重中之重。铀矿开采过程中,矿石内由镭衰变产生的氡会析出到空气中对人体产生辐射危害,因此矿山对于氡的防护是日常辐射安全防护工作的重要内容。由矿石内析出到空气中的氡被称为自由氡,而另一部分被束缚在矿石晶格内的氡为束缚氡,氡射气系数是指单位时间内从矿石析出到空气中的氡占该段时间内由镭衰变产生总氡量的份额,因此氡射气系数是表征矿石内氡析出到空气中多少最直观的系数,也是评价矿石氡辐射危害的重要参数。
目前,测量样品氡射气系数有射气法和能谱法。
采用射气法测量样品氡射气系数时,须在测定样品镭含量之后将样品放入集气瓶中密封3-4周进行氡累积,然后根据镭含量及达到镭氡长期平衡时的稳定氡浓度计算出样品的氡射气系数,采用该方法测量矿样氡射气系数需要耗费较长的时间。
采用常规能谱法测量矿样氡射气系数时,只需分别测量同一样品密封3h和3天的氡原子数,并根据两次测得的氡原子数即可计算出样品的氡射气系数,采用能谱法不须测定样品的镭含量,且操作更简便。
但是在对长期密封矿样氡射气系数的实际测量中,样品盒内的氡原子数由于镭氡建立长期平衡而处于稳定状态,因此无法利用常规能谱法直接对长期密封样品进行测量,须先将样品拆封并排除样品内的自由氡,然后再按照能谱法的常规步骤进行测量,测定一个样品至少需要4-6天,测量效率较低。
发明内容
本发明的目的是针对常规能谱法在测量长期密封矿样氡射气系数耗时长的问题而提出一种对长期密封矿样的氡射气系数进行快速测量的方法。
本发明的技术方案是:长期密封矿样氡射气系数快速测量的方法,首先用高纯锗谱仪测得长期密封矿样内平衡状态的氡原子数N1,然后拆开密封样品,并在开放环境中以搅拌的方式排出样品中的自由氡,再重新密封后测得样品内未平衡状态氡原子数N2,最后利用N1、N2及建立的长期密封矿样氡射气系数理论计算公式计算出样品的氡射气系数η。
所述的长期密封矿样氡射气系数理论计算公式建立过程如下:
已知样品平衡状态和未平衡状态氡原子数分别为N1、N2,镭的衰变系数为λ1,氡的衰变系数为λ2,由于样品中的镭始终以固定比例衰变产生自由氡和束缚氡,为了便于计算认为密封后的样品中产生自由氡的那部分镭的原子数为NRa’,镭的总原子数为NRa,常规情况下密封样品氡原子数由N2累积至N1所需的时间为T,氡射气系数为η;
将长期密封矿样氡射气系数理论计算公式建立过程分为以下三个步骤:
A.密封样品氡原子数由N2累积至N1所需时间T的理论计算公式:设样品在累积的t时刻由镭原子数NRa新产生的氡为N总(t),则:
(1)
解微分方程得:
(2)
根据样品中氡原子数达到平衡状态的时间T得到等式:
(3)
将(2)式代入(3)式中,为了简便计算,令,,并化简得
(4)
即能求得:
(5)
由于,则
(6)
B.镭原子数NRa’的理论计算公式:
根据测量知道样品平衡状态氡原子数为N1,未平衡状态氡原子数为N2,因此T时间段内由NRa’产生的自由氡不仅包含自由氡的净增量(N1-N2),还包括在时间T内产生并衰变的自由氡量,设T时间段内每个时刻的自由氡原子数为N(t),则T时间段内自由氡的变化存在以下等式:
(7)
N(t)表示由N2开始累积了时间t时刻样品中的氡原子数,则在时间t内新产生的自由氡可由公式(8)计算得到:
(8)
解(8)式微分方程得:
(9)
将(9)式代入(7)式中,得:
(10)
C.氡射气系数η的理论计算公式:
在样品密封的时间T内产生的自由氡量之比即为氡射气系数,即:
(11)
根据上述理论,时间T内密封样品产生的自由氡均由镭NRa’衰变产生,因此自由氡量为:
(12)
同理,密封样品在时间T内产生的总氡量为:
(13)
因此密封样品氡射气系数η为:
(14)
公式(14)化简为:
(15)
已知样品中镭氡建立长期平衡时,存在等式:
(16)
将(16)式代入(15)式得到:
(17)
将(10)式代入到(17)式中即可得到氡射气系数的计算公式:
(18)
由公式(18)得知,只要利用N1、N2即能够计算得到长期密封矿样的氡射气系数。
本发明与现有技术相比具有如下特点:
1、在样品氡原子数的两次测量期间不需要长时间等待氡重新累积,测量周期短,可实现对长期密封矿样的氡射气系数进行快速简便测量。
2、由于测量过程中考虑了密封样品中自由氡的自身衰变,相对于常规方法,本方法的测量结果更加准确可靠。
以下结合具体实施方式对本发明作进一步描述。
具体实施方式
长期密封矿样氡射气系数快速测量的方法,首先用高纯锗谱仪测得长期密封矿样内平衡状态的氡原子数N1,然后拆开密封样品,并在开放环境中以搅拌的方式排出样品中的自由氡,再重新密封后测得样品内未平衡状态氡原子数N2,最后利用N1、N2及建立的长期密封矿样氡射气系数理论计算公式计算出样品的氡射气系数η。
所述的长期密封矿样氡射气系数理论计算公式建立过程如下:
已知样品平衡状态和未平衡状态氡原子数分别为N1、N2,镭的衰变系数为λ1,氡的衰变系数为λ2,由于样品中的镭始终以固定比例衰变产生自由氡和束缚氡,为了便于计算认为密封后的样品中产生自由氡的那部分镭的原子数为NRa’,镭的总原子数为NRa,常规情况下密封样品氡原子数由N2累积至N1所需的时间为T,氡射气系数为η;
将长期密封矿样氡射气系数理论计算公式建立过程分为以下三个步骤:
A.密封样品氡原子数由N2累积至N1所需时间T的理论计算公式:设样品在累积的t时刻由镭原子数NRa新产生的氡为N总(t),则:
(1)
解微分方程得:
(2)
根据样品中氡原子数达到平衡状态的时间T得到等式:
(3)
将(2)式代入(3)式中,为了简便计算,令,,并化简得
(4)
即能求得:
(5)
由于,则
(6)
B.镭原子数NRa’的理论计算公式:
根据测量知道样品平衡状态氡原子数为N1,未平衡状态氡原子数为N2,因此T时间段内由NRa’产生的自由氡不仅包含自由氡的净增量(N1-N2),还包括在时间T内产生并衰变的自由氡量,设T时间段内每个时刻的自由氡原子数为N(t),则T时间段内自由氡的变化存在以下等式:
(7)
N(t)表示由N2开始累积了时间t时刻样品中的氡原子数,则在时间t内新产生的自由氡可由公式(8)计算得到:
(8)
解(8)式微分方程得:
(9)
将(9)式代入(7)式中,得:
(10)
C.氡射气系数η的理论计算公式:
在样品密封的时间T内产生的自由氡量之比即为氡射气系数,即:
(11)
根据上述理论,时间T内密封样品产生的自由氡均由镭NRa’衰变产生,因此自由氡量为:
(12)
同理,密封样品在时间T内产生的总氡量为:
(13)
因此密封样品氡射气系数η为:
(14)
公式(14)化简为:
(15)
已知样品中镭氡建立长期平衡时,存在等式:
(16)
将(16)式代入(15)式得到:
(17)
将(10)式代入到(17)式中即可得到氡射气系数的计算公式:
(18)
由公式(18)得知,只要利用N1、N2即能够计算得到长期密封矿样的氡射气系数。
在本实施例中,利用上述测量方法对A、B、C、D、E五个长期密封矿样的氡射气系数进行测量,具体测量过程如下:
首先利用高纯锗谱仪测量得到上述五个样品平衡状态的氡原子数分别为N1A=10800、N1B=10037、N1C=23043、N1D=3078、N1E=8537;然后将密封样品拆开,并在开放环境中以搅拌的方式排出样品中的自由氡;接着将样品重新密封,等待3h后测得A、B、C、D、E五个样品的未平衡状态氡原子数分别为N2A=10000、N2B=9725、N2C=20487、N2D=2980、N2E=7391,再将已测得的N1A、N1B、N1C、N1D、N1E、N2A、N2B、N2C、N2D、N2E及镭、氡的衰变常数λ1和λ2分别代入到公式(5)中计算出参数KA、KB、KC、KD、KE,其中λ1=1.37×10-11,λ2=2.1×10-6;
最后将计算出的KA、KB、KC、KD、KE 、及N1A、N1B、N1C、N1D、N1E、N2A、N2B、N2C、N2D、N2E分别代入到氡射气系数的理论计算公式(18)中:
从而计算出ηA、ηB、ηC、ηD、ηE,得到五个样品的长期密封矿样氡射气系数。
此外采用常规方法对A、B、C、D、E五个样品进行氡射气系数的测量,测量结果记为ηA’、ηB’、ηC’、ηD’、ηE’,两种方法的测量结果对比见表1。
表1 两种方法测量结果的比较
样品编号 | A | B | C | D | E |
η | 24.15% | 10.14% | 36.15% | 10.38% | 43.75% |
η’ | 26.24% | 9.10% | 33.41% | 9.56% | 47.94% |
相对偏差 | 8.29% | 10.48% | 7.88% | 8.22% | 9.14% |
从表1可知,A、B、C、D、E五个样品的两种测量结果相对偏差均在10%左右,此外相关文献表明氡射气系数能谱测量方法的不确定度通常约为10%,因此可以说上述两种方法具有较好的一致性,采用本方法能够对长期密封矿样的氡射气系数进行准确可靠、快速简便的测量。
Claims (2)
1.长期密封矿样氡射气系数快速测量的方法,其特征是:首先用高纯锗谱仪测得长期密封矿样内平衡状态的氡原子数N1,然后拆开密封样品,并在开放环境中以搅拌的方式排出样品中的自由氡,再重新密封后测得样品内未平衡状态氡原子数N2,最后利用N1、N2及建立的长期密封矿样氡射气系数理论计算公式计算出样品的氡射气系数η;
所述的长期密封矿样氡射气系数理论计算公式如下:
其中:
λ1为镭的衰变系数,λ2为氡的衰变系数。
2.如权利要求1所述的长期密封矿样氡射气系数快速测量的方法,其特征是:所述的长期密封矿样氡射气系数理论计算公式建立过程如下:
已知样品平衡状态和未平衡状态氡原子数分别为N1、N2,镭的衰变系数为λ1,氡的衰变系数为λ2,由于样品中的镭始终以固定比例衰变产生自由氡和束缚氡,为了便于计算认为密封后的样品中产生自由氡的那部分镭的原子数为NRa’,镭的总原子数为NRa,常规情况下密封样品氡原子数由N2累积至N1所需的时间为T,氡射气系数为η;
将长期密封矿样氡射气系数理论计算公式建立过程分为以下三个步骤:
A.密封样品氡原子数由N2累积至N1所需时间T的理论计算公式:设样品在累积的t时刻由镭原子数NRa新产生的氡为N总(t),则:
(1)
解微分方程得:
(2)
根据样品中氡原子数达到平衡状态的时间T得到等式:
(3)
将(2)式代入(3)式中,为了简便计算,令,,并化简得
(4)
即能求得:
(5)
由于,则
(6)
B.镭原子数NRa’的理论计算公式:
根据测量知道样品平衡状态氡原子数为N1,未平衡状态氡原子数为N2,因此T时间段内由NRa’产生的自由氡不仅包含自由氡的净增量(N1-N2),还包括在时间T内产生并衰变的自由氡量,设T时间段内每个时刻的自由氡原子数为N(t),则T时间段内自由氡的变化存在以下等式:
(7)
N(t)表示由N2开始累积了时间t时刻样品中的氡原子数,则在时间t内新产生的自由氡可由公式(8)计算得到:
(8)
解(8)式微分方程得:
(9)
将(9)式代入(7)式中,得:
(10)
C.氡射气系数η的理论计算公式:
在样品密封的时间T内产生的自由氡量之比即为氡射气系数,即:
(11)
根据上述理论,时间T内密封样品产生的自由氡均由镭NRa’衰变产生,因此自由氡量为:
(12)
同理,密封样品在时间T内产生的总氡量为:
(13)
因此密封样品氡射气系数η为:
(14)
公式(14)化简为:
(15)
已知样品中镭氡建立长期平衡时,存在等式:
(16)
将(16)式代入(15)式得到:
(17)
将(10)式代入到(17)式中即可得到氡射气系数的计算公式:
(18)
由公式(18)得知,只要利用N1、N2即能够计算得到长期密封矿样的氡射气系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611094699.8A CN106556856B (zh) | 2016-12-02 | 2016-12-02 | 长期密封矿样氡射气系数快速测量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611094699.8A CN106556856B (zh) | 2016-12-02 | 2016-12-02 | 长期密封矿样氡射气系数快速测量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106556856A CN106556856A (zh) | 2017-04-05 |
CN106556856B true CN106556856B (zh) | 2019-03-01 |
Family
ID=58446072
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611094699.8A Active CN106556856B (zh) | 2016-12-02 | 2016-12-02 | 长期密封矿样氡射气系数快速测量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106556856B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108931465B (zh) * | 2018-05-21 | 2020-07-14 | 南华大学 | 测定多孔射气介质中氡扩散系数和可运移氡产生率的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3275468B2 (ja) * | 1993-08-18 | 2002-04-15 | 春視 荒木 | ラドン・エマネーションの測定方法及び地中探査方法 |
CN103091472A (zh) * | 2013-01-25 | 2013-05-08 | 南华大学 | 一种快速测量土壤潜势氡浓度的方法 |
CN105954789A (zh) * | 2016-04-20 | 2016-09-21 | 南华大学 | 一种氡子体测量仪α能谱峰重叠修正因子的准确刻度方法 |
-
2016
- 2016-12-02 CN CN201611094699.8A patent/CN106556856B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3275468B2 (ja) * | 1993-08-18 | 2002-04-15 | 春視 荒木 | ラドン・エマネーションの測定方法及び地中探査方法 |
CN103091472A (zh) * | 2013-01-25 | 2013-05-08 | 南华大学 | 一种快速测量土壤潜势氡浓度的方法 |
CN105954789A (zh) * | 2016-04-20 | 2016-09-21 | 南华大学 | 一种氡子体测量仪α能谱峰重叠修正因子的准确刻度方法 |
Non-Patent Citations (3)
Title |
---|
Determination of the emanation coefficient and the Henry"s law constant for the groundwater radon;Kil Yong Lee et al.;《Journal of radioanalytical and nuclear chemistry》;20101130;第286卷(第2期);第381-385页 * |
Estimation the radon emanation coefficient from crystalline rocks into groundwater;Tadeusz et al.;《Applied Radiation and Isotopes》;20000915;第53卷(第3期);第473-479页 * |
不需积累氡而同时测定铀矿样中的铀、镭和射气系数;张书达;《核技术》;19810531(第5期);第33-36页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106556856A (zh) | 2017-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Roshani et al. | Combination of X-ray tube and GMDH neural network as a nondestructive and potential technique for measuring characteristics of gas-oil–water three phase flows | |
Garcia-Rios et al. | Influence of the flow rate on dissolution and precipitation features during percolation of CO2-rich sulfate solutions through fractured limestone samples | |
Thorn et al. | Three-phase flow measurement in the petroleum industry | |
Huyskens et al. | Evaluation of temperature-time conditions for the chemical abrasion treatment of single zircons for U–Pb geochronology | |
Palmerini et al. | The RGB and AGB Star Nucleosynthesis in Light of the Recent 17O (p, α) 14N and 18O (p, α) 15N Reaction-rate Determinations | |
Polyansky et al. | Calculation of rotation-vibration energy levels of the ammonia molecule based on an ab initio potential energy surface | |
Thompson et al. | Chemical properties of Berkelium | |
EP3667304B1 (en) | Crystal-phase quantitative analysis device, crystal-phase quantitative analysis method, and crystal-phase quantitative analysis program | |
CN103954596A (zh) | 原子荧光光谱法间接测定土壤中铀的含量的方法 | |
CN106556856B (zh) | 长期密封矿样氡射气系数快速测量的方法 | |
Alabidi et al. | Generating primordial black holes via hilltop-type inflation models | |
ITPD20080204A1 (it) | Apparato metodo e apparato di rilevamento della frazione di origine remota di radon presente in un sito di misura. | |
Kirkegaard et al. | Elucidation of the structure and vibrational spectroscopy of synthetic metaschoepite and its dehydration product | |
Tranter et al. | Measuring flow profiles in heated miniature reactors with X-ray fluorescence spectroscopy | |
Blumentritt et al. | Unravelling the origin of the yellow-orange luminescence in natural and synthetic scapolites | |
Oher et al. | Influence of Alkaline Earth Metal Ions on Structures and Luminescent Properties of Na m M n UO2 (CO3) 3 (4–m–2 n)–(M= Mg, Ca; m, n= 0–2): Time-Resolved Fluorescence Spectroscopy and Ab Initio Studies | |
Elliott et al. | The gallium anomaly | |
Zamzamian et al. | Optimization of the marinelli beaker dimensions using genetic algorithm | |
Cooper et al. | Minimum detectable concentration and concentration calculations | |
Li et al. | Effect of temperature on the protonation of N-(2-hydroxyethyl) ethylenediamine-N, N′, N′-triacetic acid in aqueous solutions: Potentiometric and calorimetric studies | |
Suerfu et al. | Zone refining of ultrahigh-purity sodium iodide for low-background detectors | |
CN102692385B (zh) | 一种测定化妆品中滑石粉的方法 | |
CN201335814Y (zh) | 油气田中多相流指标的真双能χ射线测量装置 | |
Hickam et al. | Neptunyl peroxide chemistry: synthesis and spectroscopic characterization of a neptunyl triperoxide compound, Ca2 [NpO2 (O2) 3]· 9H2O | |
Kahoul et al. | New procedure calculation of photon-induced Kβ/Kα intensity ratios for elements 16S to 92U |
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 |