CN101644779A - 一种γ辐射场剂量率快速测量方法 - Google Patents

一种γ辐射场剂量率快速测量方法 Download PDF

Info

Publication number
CN101644779A
CN101644779A CN200910091590A CN200910091590A CN101644779A CN 101644779 A CN101644779 A CN 101644779A CN 200910091590 A CN200910091590 A CN 200910091590A CN 200910091590 A CN200910091590 A CN 200910091590A CN 101644779 A CN101644779 A CN 101644779A
Authority
CN
China
Prior art keywords
data
radiation field
dose rate
gamma radiation
mean value
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
Application number
CN200910091590A
Other languages
English (en)
Other versions
CN101644779B (zh
Inventor
邓长明
任熠
张世让
杜向阳
张佳
程昶
孟丹
刘芸
王志刚
侯磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanxi Zhongfu Technology Co., Ltd.
Original Assignee
China Institute for Radiation Protection
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Institute for Radiation Protection filed Critical China Institute for Radiation Protection
Priority to CN2009100915902A priority Critical patent/CN101644779B/zh
Publication of CN101644779A publication Critical patent/CN101644779A/zh
Application granted granted Critical
Publication of CN101644779B publication Critical patent/CN101644779B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Radiation (AREA)

Abstract

本发明公开了一种γ辐射场剂量率快速测量方法,属于辐射场测量技术领域。现有的方法对辐射场变化后响应的速度较慢,从而影响了测量速度和准确性。本发明所述方法在现有剂量率滑动平均算法的基础上,增加了辐射场水平的判断,将辐射场水平作为本次测量结束的标志,不仅实现了剂量率快速测量,而且加快了对辐射场变化的响应时间,使测量结果更加准确。

Description

一种γ辐射场剂量率快速测量方法
技术领域
本发明属于辐射场测量技术领域,具体涉及一种γ辐射场剂量率快速测量方法。
背景技术
辐射测量仪表类型可分为核素能量测量和辐射场强度测量,辐射场强度的测量一般用计数率或剂量率表示,剂量率是指单位时间内的累积剂量。老式仪表由于电子技术落后,采用定时计数,只能手动测量。随着单片机技术的发展,给剂量率快速智能化测量带来可能,仪器的智能化程度越来越高。剂量率测量仪表适用于核电站、核军工、保健物理、工业测量、放射化学、核医学、核材料运输与储存、环境监测等领域的辐射测量。
在探测器灵敏度确定的情况下,剂量率测量的准确度与测量时间有密切关系。测量时间短,采集的计数少,统计涨落就大,测量精度低;反之测量时间长,采集的计数多,统计涨落小,测量精度高。现有的剂量率测量仪器一般采用滑动平均算法计算剂量率,当测量时间过长时,仪器对辐射剂量场变化的响应时间比较慢,导致测量的结果不准确,而且对人员的防护也是不利的。
发明内容
针对现有技术中存在的缺陷,本发明的目的是提供一种快速、准确地测量γ辐射场剂量率的方法。
为了实现上述目的,本发明采用的技术方案如下:
一种γ辐射场剂量率快速测量方法,该方法针对γ辐射场中的剂量率每秒采集一组数据,计算该组数据平均值,并将该组数据平均值插入到队列尾;当队列中数据的个数超过预设的最大队列个数N时,将队列头数据删除,将新采集的一组数据的平均值插入到队列尾;
在开始下一秒采集数据时,计算队列中所有数据的平均值,并将其显示在测量设备上;
从第2秒开始,依次判断下列条件是否成立:①该秒采集的一组数据的平均值是否大于设定值P,所述P为正整数,4cps≤P≤6cps,②该秒采集的一组数据的变化方向是否一致,③该秒采集的所有数据与当前队列中所有数据的平均值之差是否均大于或均小于K倍的标准偏差;如果上述条件均满足,则结束本次测量,重新开始测量,测量设备上的显示结果为本次测量的最终结果;如果上述条件有一个不满足,则继续采集数据。
如上所述的γ辐射场剂量率快速测量方法,其中,每秒采集的每组数据的个数为2。
如上所述的γ辐射场剂量率快速测量方法,其中,N为正整数,且60≤N≤100。
如上所述的γ辐射场剂量率快速测量方法,其中,P=5cps。
如上所述的γ辐射场剂量率快速测量方法,其中,判断该秒采集的一组数据的变化方向是否一致的方法为:该秒采集的所有数据是否均大于或均小于当前队列中所有数据的平均值,如是,则变化方向一致,否则,变化方向不一致。
如所述的γ辐射场剂量率快速测量方法,其中,标准偏差的计算公式为:
Figure G2009100915902D00021
M1、M2...Mn分别为第n秒采集的所有数据的平均值。
如上所述的γ辐射场剂量率快速测量方法,其中,K为正整数,且3≤K≤5。
本发明所述方法,在剂量率滑动平均算法的基础上,增加了γ辐射场水平的评断,不仅实现了剂量率的快速测量,而且提高了对γ辐射场变化的响应速度,进而提高了剂量率测量的准确性。
附图说明
图1是具体实施方式中测量γ辐射场剂量率的方法流程图;
图2a是具体实施方式中当测量到60秒时队列中数据示意图,图2b是具体实施方式中当测量到61秒时队列中数据示意图。
具体实施方式
本发明的核心思想是:在现有剂量率滑动平均算法的基础上,增加辐射场水平的判断,将辐射场水平作为本次测量结束的标志,实现了剂量率快速测量的基础上,加快对辐射场变化的响应时间,从而使测量结果更加准确。
下面结合具体实施方式和附图对本发明进行详细描述。
图1出示了本发明所述方法一种实施方式的流程,该方法使用GM计数管针对γ辐射场中的剂量率每秒采集一组数据,计算该组数据平均值,并将该组数据平均值插入到队列尾;当队列中数据的个数超过预设的最大队列个数N(N为正整数且60≤N≤100)时,将队列头数据删除,将新采集的一组数据的平均值插入到队列尾。
在开始下一秒采集数据时,计算队列中所有数据的平均值,并将其显示在测量设备上。
从第2秒开始,依次判断下列条件是否成立:①该秒采集的一组数据的平均值是否大于设定值P,P为正整数,且4cps≤P≤6cps,(cps表示每秒的计数个数),②该秒采集的一组数据的变化方向是否一致,③该秒采集的所有数据与当前队列中所有数据的平均值之差是否均大于或均小于K(K为正整数且3≤K≤5)倍的标准偏差。如果上述条件均满足,则结束本次测量,重新开始测量,测量设备上的显示结果为本次测量的最终结果;如果上述条件有一个不满足,则开始下一秒采集数据。
上述方法中,判断该秒采集的一组数据的变化方向是否一致的方法为:该秒采集的所有数据是否均大于或均小于当前队列中所有数据的平均值,如是,则变化方向一致,否则,变化方向不一致。
上述方法中,标准偏差的计算公式为:
Figure G2009100915902D00041
其中M1、M2...Mn分别为第n秒采集的所有数据的平均值。
实施例
下面以每秒采集2个测量数据、N等于60、P等于5cps,K等于4为例,对本发明测量γ辐射场剂量率的过程进行详细描述,如图1所示。
第1秒采集两个测量数据A11、A12,计算平均值M1 M 1 = A 11 + A 12 2 , 将M1插入队列尾,并显示在测量设备上。
第2秒采集两个测量数据A21、A22,计算平均值M2 M 2 = A 21 + A 22 2 , 将M2插入到队列尾,并将M1和M2(队列中所有数据)的平均值显示在测量设备上。
然后,判断平均值M2是否大于5cps。如果大于,则继续判断测量数据的变化方向是否一致,即A21和A22是否均大于或均小于M1和M2(队列中所有数据)的平均值。
也就是判断 A 21 > M 1 + M 2 2 A 22 > M 1 + M 2 2 是否同时成立,或者 A 21 < M 1 + M 2 2 A 22 < M 1 + M 2 2 是否同时成立,如是,则测量数据的变化方向一致。
如果测量数据的变化方向一致,则继续判断A21、A22与M1、M2(队列中所有数据)的平均值之差是否均大于4倍的标准偏差。第2秒采集数据的标准偏差的计算公式为:
Figure G2009100915902D00048
即判断 ( A 21 - M 1 + M 2 2 ) > 4 M 1 + M 2 2 ( A 22 - M 1 + M 2 2 ) > 4 M 1 + M 2 2 是否同时成立。
如果上述条件均满足,则结束本次测量,进行下一次测量,测量设备上显示的数据便是本次测量的结果。如果上述条件有一个不满足,则继续进行数据采集。接下来对采集到的数据处理方式类推第2秒对采集到的数据的处理方式,直到第61秒。第60秒时的队列数据如图2a所示。
第61秒采集2个测量数据A61、A62,计算平均值M61 M 61 = A 61 + A 62 2 . 删除队列头数据M1,将M61插入队列尾,如图2b所示。并计算队列中所有数据M2、M3......M61的平均值显示在测量设备上。
然后,判断M61是否大于5cps。如果大于,则继续判断 A 61 > M 2 + M 3 + . . . + M 61 60 A 62 > M 2 + M 3 + . . . + M 61 60 是否同时成立,或者 A 61 < M 2 + M 3 + . . . + M 61 60 A 62 < M 2 + M 3 + . . . + M 61 60 是否同时成立。如果同时成立,则说明变化方式一致,继续判断 ( A 61 - M 2 + M 3 + . . . + M 61 60 ) > 4 M 2 + M 3 + . . . + M 61 60 ( A 62 - M 2 + M 3 + . . . + M 61 60 ) > 4 M 2 + M 3 + . . . + M 61 60 是否同时成立。
如果上述条件均满足,则结束本次测量,进行下一次测量,测量设备上显示的数据便是本次测量的结果。如果上述条件有一个不满足,则继续进行数据采集。接下来对采集的数据的处理方式类推第61秒对采集的数据的处理方式,直到本次测量结束。
最大队列数N以及标准偏差倍数K的取值都可以根据需要进行调整,以上数值仅作为例子加以参考,本领域的技术人员完全可以根据需要选择不同的数值,一般来说N为正整数且60≤N≤100,K为正整数且3≤K≤5,都是可以选用的数据。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (7)

1.一种γ辐射场剂量率快速测量方法,所述方法针对γ辐射场中的剂量率每秒采集一组数据,计算该组数据平均值,并将该组数据平均值插入到队列尾;当队列中数据的个数超过预设的最大队列个数N时,将队列头数据删除,将新采集的一组数据的平均值插入到队列尾;
在开始下一秒采集数据时,计算队列中所有数据的平均值,并将其显示在测量设备上;
其特征在于:从第2秒开始,依次判断下列条件是否成立:①该秒采集的一组数据的平均值是否大于设定值P,所述P为正整数,4cps≤P≤6cps,②该秒采集的一组数据的变化方向是否一致,③该秒采集的所有数据与当前队列中所有数据的平均值之差是否均大于或均小于K倍的标准偏差;如果上述条件均满足,则结束本次测量,重新开始测量,测量设备上的显示结果为本次测量的最终结果;如果上述条件有一个不满足,则继续采集数据。
2.如权利要求1所述的一种γ辐射场剂量率快速测量方法,其特征在于:每秒采集的每组数据的个数为2。
3.如权利要求1所述的一种γ辐射场剂量率快速测量方法,其特征在于:所述N为正整数,且60≤N≤100。
4.如权利要求1所述的一种γ辐射场剂量率快速测量方法,其特征在于:所述P=5cps。
5.如权利要求1至4之一所述的一种γ辐射场剂量率快速测量方法,其特征在于,所述判断该秒采集的一组数据的变化方向是否一致的方法为:该秒采集的所有数据是否均大于或均小于当前队列中所有数据的平均值,如是,则变化方向一致,否则,变化方向不一致。
6.如权利要求1至4之一所述的一种γ辐射场剂量率快速测量方法,其特征在于,所述标准偏差的计算公式为:
Figure A2009100915900002C1
其中M1、M2…Mn分别为第n秒采集的所有数据的平均值。
7.如权利要求1至4之一所述的一种γ辐射场剂量率快速测量方法,其特征在于,所述K为正整数,且3≤K≤5。
CN2009100915902A 2009-08-28 2009-08-28 一种γ辐射场剂量率快速测量方法 Active CN101644779B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100915902A CN101644779B (zh) 2009-08-28 2009-08-28 一种γ辐射场剂量率快速测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100915902A CN101644779B (zh) 2009-08-28 2009-08-28 一种γ辐射场剂量率快速测量方法

Publications (2)

Publication Number Publication Date
CN101644779A true CN101644779A (zh) 2010-02-10
CN101644779B CN101644779B (zh) 2012-10-31

Family

ID=41656742

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100915902A Active CN101644779B (zh) 2009-08-28 2009-08-28 一种γ辐射场剂量率快速测量方法

Country Status (1)

Country Link
CN (1) CN101644779B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111506857A (zh) * 2020-03-20 2020-08-07 中国辐射防护研究院 一种水路运输放射性物品时公众集体剂量计算方法
CN112330931A (zh) * 2020-09-09 2021-02-05 中国原子能科学研究院 一种优化的临界事故报警方法及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN87211165U (zh) * 1987-07-31 1988-05-11 哈尔滨市东方医疗电子仪器厂 个人射线累积存储剂量仪
JPH0197893A (ja) * 1987-10-12 1989-04-17 Nippon Sci Kk 低エネルギーβ線用自動効率トレーサ法
CN1033720C (zh) * 1993-12-23 1997-01-01 中国辐射防护研究院 环境γ辐射监测探头温度补偿方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111506857A (zh) * 2020-03-20 2020-08-07 中国辐射防护研究院 一种水路运输放射性物品时公众集体剂量计算方法
CN111506857B (zh) * 2020-03-20 2023-05-02 中国辐射防护研究院 一种水路运输放射性物品时公众集体剂量计算方法
CN112330931A (zh) * 2020-09-09 2021-02-05 中国原子能科学研究院 一种优化的临界事故报警方法及系统
CN112330931B (zh) * 2020-09-09 2022-03-11 中国原子能科学研究院 一种优化的临界事故报警方法及系统

Also Published As

Publication number Publication date
CN101644779B (zh) 2012-10-31

Similar Documents

Publication Publication Date Title
CN103323615B (zh) 一种通过加速度传感器计算步行速度的移动终端及方法
CN205175854U (zh) 一种径流与泥沙的自动监测装置
CN104181878B (zh) 对碳排放过程的监测控制方法
CN109461095A (zh) 一种用于非侵入式负荷辨识的负荷事件检测的数据处理方法及设备
CN104678425B (zh) 一种基于液体闪烁探测器的快中子多重性测量分析方法
CN108665096A (zh) 人流量告警方法及装置
CN106597096A (zh) 一种时钟频率监测方法
CN105759216A (zh) 一种软包锂离子电池荷电状态估算方法
CN105788666A (zh) 一种核反应堆临界过程控制方法
CN109613440A (zh) 电池的分级方法、装置、设备和存储介质
CN101644779B (zh) 一种γ辐射场剂量率快速测量方法
CN107728192A (zh) 多探测器的核素识别系统及方法
CN104198811B (zh) 低频信号频率测量方法和装置
CN105387913A (zh) 基于指数关系和支持向量回归的粮仓重量检测方法及装置
CN104237632A (zh) 一种实时双向电能独立计量系统
CN105353371B (zh) 基于ar谱扩展分形的海面雷达目标检测方法
CN109117984A (zh) 稻田径流量预测及氮磷流失估测方法
CN102830128B (zh) α能谱测量空气中222Rn、220Rn子体的方法
CN103235350B (zh) 放射性测井仪器稳定性检测及刻度方法与装置
CN115950495A (zh) 一种超声水表测量频率调节方法
CN104182377A (zh) 一种基于β似然函数的参数估计方法
CN107064808A (zh) 电池的soc估计方法及装置
CN107768754A (zh) 锂离子电池电压挑选容量的方法及其应用
CN108982310A (zh) 一种降低人工放射性气溶胶监测仪虚警率的方法
CN106021855A (zh) 一种反应堆周期计算方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20100210

Assignee: Shanxi Zhongfu Technology Co., Ltd.

Assignor: China Radiation Protection Research Inst.

Contract record no.: 2015140000005

Denomination of invention: Gamma radiation field dose rate rapid measuring method

Granted publication date: 20121031

License type: Exclusive License

Record date: 20150618

LICC Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20181109

Address after: 030006 420 Changzhi Road, hi tech Zone, Taiyuan, Shanxi

Patentee after: Shanxi Zhongfu Technology Co., Ltd.

Address before: 030006 No. 102 Xuefu Street, Taiyuan City, Shanxi Province

Patentee before: China Radiation Protection Research Inst.