CN104765060B - 一种测氡仪刻度因子的定值方法 - Google Patents
一种测氡仪刻度因子的定值方法 Download PDFInfo
- Publication number
- CN104765060B CN104765060B CN201510137825.2A CN201510137825A CN104765060B CN 104765060 B CN104765060 B CN 104765060B CN 201510137825 A CN201510137825 A CN 201510137825A CN 104765060 B CN104765060 B CN 104765060B
- Authority
- CN
- China
- Prior art keywords
- measure
- radon
- cycle
- calibration factor
- emanometer
- 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.)
- Expired - Fee Related
Links
Landscapes
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种测氡仪刻度因子的定值方法,先确定测量周期为3600s时的刻度因子K0,然后根据下式确定测量周期不为3600s时的测氡仪刻度因子Kx:该方法能快速地对测氡仪刻度因子进行测定,易于实施。
Description
技术领域
本发明涉及一种测氡仪刻度因子的定值方法。
背景技术
环境中氡污染水平评价、氡的特性研究及地震预报和资源勘测等示踪应用中需要对氡进行连续可靠的测量。当环境氡浓度起伏较大时,目前现有测氡仪难以准确地测量其变化规律,无法满足快速可靠测量的要求。为实现快速追踪环境氡浓度的变化需要采用静电收集法收集222Rn衰变产生的第一代子体218Po,并对其衰变产生的高能α粒子能量进行识别并分别计数,由此反推得到氡的浓度,各个不同测量周期反推得到氡的浓度需要进行多次反复测量才能得到较的刻度因子,该过程操作繁琐且测量误差大。氡测量时短测量周期测量当次测量受前一次测量结果影响较大,虽然可以采用迭代修正进行修正,但是对刻度因子的刻度带来较大误差,参见文献:何正忠,肖德涛,赵桂芝,等,吴喜军.连续测氡方法迭代修正因子的理论计算与实验测定[J].原子能科学技术,2013,47(6):1040-1043。
因此,有必要设计一种测氡仪刻度因子的定值方法。
发明内容
本发明所要解决的技术问题是提供一种测氡仪刻度因子的定值方法,该测氡仪刻度因子的定值方法能快速地对测氡仪刻度因子进行测定,易于实施。
发明的技术解决方案如下:
一种测氡仪刻度因子的定值方法,先确定测量周期为3600s时的刻度因子K0,然后根据下式确定测量周期不为3600s时的测氡仪刻度因子Kx:
其中,λR和λP分别为氡222Rn和218Po的衰变常量;
TM0表示测量周期T=3600s时的氡测量时间;【TM0=T-Ts是3300s】;
TMx表示任意测量周期T=3600s外的其它测量周期的氡测量时间;【换言之,TMx表示测量周期T≠3600s的取样完成后氡测量的时间】,有TMx=T-Ts,Ts为取样时间,60s≤Ts≤300s;
其中,为被测环境氡浓度;【是已知量】;
ΔNP(T0)′是半导体探测器【如PIPS探测器】在时间TM0收到的218Po衰变产生的α粒子数(是直接测量得到的计数);
Kx对应的测量周期为900s和1800s。
被测环境氡浓度由标准氡室的标准测量仪器PQ2000测量得到。
本发明的计算公式推导过程说明如下:
静电收集法测氡是在测量室与半导体探测器之间建立高压静电场,采用抽气泵取样时长为Ts,使含氡空气以一定的流速经高效过滤器后进入测量室,由222Rn衰变产生的第一代子体218Po带正电荷会在静电场的作用下吸附到半导体探测器表面,当其子体进一步衰变时,半导体探测器对其衰变产生的高能α粒子能量进行识别并在时间T-Ts内计数,最后根据α粒子计数与氡浓度的关系(即刻度因子)确定氡浓度。
静电收集法快速测222Rn浓度原理公式为:
为被测环境氡浓度,R表示温湿度修正因子[6],ΔN表示氡衰变的第一代子体218Po进一步衰变时产生的6.00MeV的α粒子计数,K为测氡刻度因子,(1)式中的η表示探测效率,V表示测量室体积,R表示不同温湿度条件下的温湿度因子。
探测器收到的218Po衰变产生的α粒子数ΔNP(T0)′=ηVRΔNP(T0),(1)式可化为:
设氡222Rn、218Po的衰变常量分别为λR、λP,放射性活度分别为AR、AP;在时刻t,,222Rn、218Po的原子核数分别为NR(t)、NP(t);在t=0时,只有222Rn,即NP(t)=0。
222Rn的放射性活度为
关于218Po单位时间核数目的变化一方面以速率λRNR(t)从222Rn中产生,另一方面又以速率λPNP(t)衰变为214Bi,即
利用(3)式对(4)式微分方程求解,得
由(5)式得218Po的放射性活度为
从时刻Ts~T内218Po衰变的原子数,由(6)式积分得
化简(8)式得
式中:T为测量周期,TS为取样时间,TM为测量时间,当TS=300second时,取T=600s测量周期(TM=3300s)刻度因子为K0,由(2)式得:
Kx表示其它不同测量周期Tx刻度因子,由(1)、(9)式可得
将(8)式代入(11)式得
其中,设测量周期为3600s的刻度因子为K0,
为被测环境氡浓度,R表示温湿度修正因子[6],ΔNP(T0)′是探测器收到的218Po衰变产生的α粒子数,是直接测量得到的计数。
λR和λP分别为氡222Rn和218Po的衰变常量;
TM0表示测量周期T=3600s时的氡测量时间;TMx表示测量周期T=3600s外的其它测量周期的氡测量时间;
(12)式表时同台测氡仪其刻度因子仅与测量周期和取样时间有关,与收集效率及探测效率无关。当同台测氡仪通过实验获得60min(3600s)测量周期的刻度因子K0后,通过(12)式可计算出该测氡仪在任意测量周期和相应抽气取样时间下的刻度因子,实现对氡浓度的快速标定。计算不同测量周期下的刻度因子,取Ts=300s,λR=2.1×10-6s-1,λP=3.7×10-3s-1,当测量周期T分别为900s、1800s时,则TM分别为600s、1500s时其刻度因子K15、K30值为
K15=8.42K0 (13)
K30=2.46K0 (14)
从(13)、(14)式可看出,固定取样时间的情况下,测量周期越长,刻度因子会越小。
有益效果:
本发明的测氡仪刻度因子的定值方法,该方法是通过实验测得60min测量周期的刻度因子后,再依据该测量周期收集到218Po衰变产生的α粒子数与其它不同测量周期收集到218Po衰变产生的α粒子数之比与60min测量周期刻度因子的乘积来快速定值测氡仪刻度因子。60min测量周期内218Po衰变产生的α粒子数与其它不同测量周期218Po衰变产生的α粒子数之比即等价于公式中的通过对3台改进型NRL-II型测氡仪进行浓度刻度,并对测量结果的影响因素进行了讨论。结果表明,采用该方法快速定值的刻度因子与实验测量的刻度因子符合很好,同台测氡仪的刻度因子仅与测量时间有关。
采用该方法,同一台测氡仪只需要刻度1个测量周期,就能得到所有测量周期的刻度因子,极大的节药了测氡仪的刻度时间,同时还减少了统计涨落带来的误差。
附图说明
图1为测氡仪测量示意图。
具体实施方式
以下将结合附图和具体实施例对本发明做进一步详细说明:
一种测氡仪刻度因子的定值方法,先确定测量周期为3600s时的刻度因子K0,然后根据下式确定测量周期不为3600s时的测氡仪刻度因子Kx:
其中,λR和λP分别为氡222Rn和218Po的衰变常量;
TM0表示测量周期T=3600s时的氡测量时间;【TM0=T-Ts是3300s】;
TMx表示任意测量周期T=3600s外的其它测量周期的氡测量时间;【换言之,TMx表示测量周期T≠3600s的取样完成后氡测量的时间】,有TMx=T-Ts,Ts为循环取样时间,60s≤Ts≤300s;
其中,为被测环境氡浓度;【是已知量】;
ΔNP(T0)′是半导体探测器【如PIPS探测器】在时间TM0收到的218Po衰变产生的α粒子数(是直接测量得到的计数);
Kx对应的测量周期为900s和1800s。
被测环境氡浓度由标准氡室的标准测量仪器PQ2000测量得到。
实验方案
测氡仪刻度因子快速刻度实验验证装置由标准氡室、1台PQ2000、3台改进型NRL-II氡连续测量仪(改进型NRL-II氡连续测量仪的介绍:为了获得不同条件,将a号测氡仪将更换为的探测器,使得探测效率η有所区别,b号测氡仪将原来的测量室由原来950ml更换为785ml,使得收集体测量室体积有所改变。后面有说明)构成。实验测量得到刻度因子K值过程如下:整点开启3台测氡仪与PQ2000同步进行测量,首先测量仪器自身和空气的本底计数nb,再从标准氡室抽取含氡空气,并测量该仪器测量室内218Po粒子发射α粒子的总计数ΔNP(T)′,将PQ2000测量的氡浓度的平均值作为氡室标准浓度依据(10)式可得该测量周期的刻度因子K值。为了得到稳定室氡浓度采用南华大学总体积为25m3的大体积标准氡室,整个系统的漏气率小于0.1%,内部还设置了温湿度调节装置,可实现氡活度浓度和温湿度的调节和维持动态稳定。为减少放射性的统计误差影响,使氡室氡浓度控制在2000~2500Bq·m-3范围内;为了保证温湿度修正因子R不变,温、湿度分别稳定在(20±2)℃、(45±5)%。将a号测氡仪将更换为的探测器,使得探测效率η有所区别,b号测氡仪将原来的测量室由原来950ml更换为785ml,使得收集体测量室体积有所改变。
3结果分析与讨论
分别对3台仪器在周期为60min(3600s)的情况下进行3次测量,再与PQ2000测量的氡浓度的平均值作为氡室标准浓度,对3次测量的结果求平均再依据(10)式得到每台仪器测量周期为60min(3600s)的刻度因子,3台测氡仪通过实验刻度得到的刻度因子K0a=1.62Bq.m-3.cpm-1、K0b=2.01Bq.m-3.cpm-1、
K0c=1.57Bq.m-3.cpm-1,再将K0a、K0b、K0c代入(12)、(13)式得K15、K30的理论计算值,再对3台仪器在测量周期为15min(900s)、30min(1800s)的分别进行3次测量,与PQ2000测量的氡浓度的平均值作为氡室标准浓度依据(2)式可得该测量周期的刻度因子K值。计算及测量如果如表1所示。
表1 3台NRL-II测氡仪理论计算刻度因子与实验刻度因子比较
表1测量结果显示,其实验测量值与理论值符合较好,说明实验结果是可靠的。根据计算得到的刻度因子看出,在900s和1800s测量周期的刻度因子是3600s测量周期刻度因子的8.60和2.47倍,为方便测氡仪的刻度需要采用该方法进行快速刻度是很有必要的,以提高仪器连续测量的准确性。
对同台NRL-II氡连续测量仪的刻度因子是60min测量周期内218Po衰变产生的α粒子数与其它不同测量周期218Po衰变产生的α粒子数之比与60min测量周期刻度因子的乘积来快速定值测氡仪刻度因子,60min测量周期内218Po衰变产生的α粒子数与其它不同测量周期218Po衰变产生的α粒子数之比即等价于公式中的通过对NRL-II型测氡仪进行浓度刻度,并对测量结果的影响因素进行了讨论。结果表明,通过方法来快速刻度测氡仪在不同测量周期时的刻度因子与实验获得的刻度因子符合很好。
Claims (3)
1.一种测氡仪刻度因子的定值方法,其特征在于,先确定测量周期为3600s时的刻度因子K0,然后根据下式确定测量周期不为3600s时的测氡仪刻度因子Kx:
其中,λR和λP分别为氡222Rn和218Po的衰变常量;
TM0表示测量周期T=3600s时的氡测量时间;
TMx表示任意测量周期T=3600s外的其它测量周期的氡测量时间;有TMx=T-Ts,Ts为取样时间,60s≤TS≤300s;
其中,为被测环境氡浓度;
ΔNP(T0)′是半导体探测器在时间TM0收到的218Po衰变产生的α粒子数;
2.根据权利要求1所述的测氡仪刻度因子的定值方法,其特征在于,Kx对应的测量周期为900s和1800s。
3.根据权利要求1所述的测氡仪刻度因子的定值方法,其特征在于,被测环境氡浓度由标准氡室的标准测量仪器PQ2000测量得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510137825.2A CN104765060B (zh) | 2015-03-26 | 2015-03-26 | 一种测氡仪刻度因子的定值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510137825.2A CN104765060B (zh) | 2015-03-26 | 2015-03-26 | 一种测氡仪刻度因子的定值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104765060A CN104765060A (zh) | 2015-07-08 |
CN104765060B true CN104765060B (zh) | 2017-05-17 |
Family
ID=53647007
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510137825.2A Expired - Fee Related CN104765060B (zh) | 2015-03-26 | 2015-03-26 | 一种测氡仪刻度因子的定值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104765060B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105954789B (zh) * | 2016-04-20 | 2018-05-01 | 南华大学 | 一种氡子体测量仪α能谱峰重叠修正因子的准确刻度方法 |
CN108363089B (zh) * | 2018-02-07 | 2019-02-22 | 衡阳师范学院 | 氡析出率快速定值方法和装置 |
CN108646286B (zh) * | 2018-06-01 | 2019-12-24 | 衡阳师范学院 | 一种连续可靠调控氡析出率的参考装置和方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644700A (zh) * | 2008-12-10 | 2010-02-10 | 南华大学 | 测量介质表面氡析出率的方法及装置 |
CN101957328A (zh) * | 2010-08-14 | 2011-01-26 | 南华大学 | 连续可靠测量空气中氡浓度的方法及装置 |
CN102043159A (zh) * | 2010-10-22 | 2011-05-04 | 衡阳师范学院 | 一种连续快速跟踪氡浓度变化的测量方法 |
CN102141527A (zh) * | 2010-12-28 | 2011-08-03 | 衡阳师范学院 | 采用0阶近似快速测量氡浓度的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100734222B1 (ko) * | 2006-06-12 | 2007-07-06 | 박영웅 | 라돈농도 측정기의 상대교정 방법 및 장치 |
-
2015
- 2015-03-26 CN CN201510137825.2A patent/CN104765060B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644700A (zh) * | 2008-12-10 | 2010-02-10 | 南华大学 | 测量介质表面氡析出率的方法及装置 |
CN101957328A (zh) * | 2010-08-14 | 2011-01-26 | 南华大学 | 连续可靠测量空气中氡浓度的方法及装置 |
CN102043159A (zh) * | 2010-10-22 | 2011-05-04 | 衡阳师范学院 | 一种连续快速跟踪氡浓度变化的测量方法 |
CN102141527A (zh) * | 2010-12-28 | 2011-08-03 | 衡阳师范学院 | 采用0阶近似快速测量氡浓度的方法 |
Non-Patent Citations (4)
Title |
---|
222Rn-220Rn分辨探测器的刻度与比对;武云云 等;《核技术》;20081231;第31卷(第12期);920-924 * |
a novel algorithm for quick and continuous tracing the change of radon concentration in environment;yanliang tan et al.;《Review of scientific instruments》;20110430;第82卷(第4期);043503-1至043503-4 * |
EXPERIMENTAL AND THEORETICAL CONSIDERATIONS ON THE CALIBRATION FACTOR K BETWEEN a-ACTIVITY CONCENTRATION AND TRACK DENSITY FOR APPLICATION IN RADON DOSIMETRY;A. F. MAGED et al.;《Journal of Radioanalytical and Nuclear Chemistry》;19931231;第170卷(第2期);423-431 * |
连续测氡方法迭代修正因子的理论计算与实验测定;何正忠 等;《原子能科学技术》;20130630;第47卷(第6期);1040-1043 * |
Also Published As
Publication number | Publication date |
---|---|
CN104765060A (zh) | 2015-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bronson | Validation of the accuracy of the LabSOCS software for mathematical efficiency calibration of Ge detectors for typical laboratory samples | |
Bartley et al. | Respirable aerosol sampler performance testing | |
McHenry Jr et al. | Axial mixing of binary gas mixtures flowing in a random bed of spheres | |
Mishra et al. | Superhorizon fluctuations and acoustic oscillations in relativistic heavy-ion collisions | |
Townsend | The diffusion behind a line source in homogeneous turbulence | |
Soo | Instrumentation for fluid particle flow | |
CN104765060B (zh) | 一种测氡仪刻度因子的定值方法 | |
CN102680999B (zh) | 静电收集测氡探测效率的温湿度自动补偿方法 | |
CN105425267A (zh) | 一种氡析出率的测量系统及方法 | |
CN106021710A (zh) | 基于大气电离层参数的震前卫星轨道异常识别方法 | |
CN104614753B (zh) | 一种连续测量介质表面氡析出率的方法和装置 | |
CN106707325A (zh) | 用于测量一定时间段内介质表面氡释放总量的方法与系统 | |
Megantoro et al. | Validation method for digital flow meter for fuel vendors | |
CN103900666B (zh) | 膜式燃气表长期运行的计量精度控制方法 | |
CN101957328A (zh) | 连续可靠测量空气中氡浓度的方法及装置 | |
CN108572387A (zh) | 一种体源测量探测器的校准方法 | |
CN105806756B (zh) | 一种利用气象参数修正室内外pm2.5监测结果的方法 | |
CN109061304A (zh) | 一种钯在极稀薄氢气环境中电导率变化量计算方法 | |
CN104360370B (zh) | 一种连续测量钍射气子体浓度的方法 | |
Zhang et al. | The geometric factor of high energy protons detector on FY-3 satellite | |
CN108363089A (zh) | 氡析出率快速定值方法和装置 | |
CN106940448A (zh) | 一种铀浓缩厂分离功率在线监测方法 | |
Yang et al. | Optimized design of the restrictor for a MEMS thermal gas flowmeter | |
Mamedov et al. | Development of an ultra-sensitive radon detector for the SuperNEMO experiment | |
Kelly et al. | Turbulent interchange in triangular array bare rod bundles |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170517 Termination date: 20180326 |
|
CF01 | Termination of patent right due to non-payment of annual fee |