CN113670783A - 基于β射线法的气体中颗粒物浓度直读测算方法及装置 - Google Patents
基于β射线法的气体中颗粒物浓度直读测算方法及装置 Download PDFInfo
- Publication number
- CN113670783A CN113670783A CN202111079440.7A CN202111079440A CN113670783A CN 113670783 A CN113670783 A CN 113670783A CN 202111079440 A CN202111079440 A CN 202111079440A CN 113670783 A CN113670783 A CN 113670783A
- Authority
- CN
- China
- Prior art keywords
- gas
- beta
- ray
- concentration
- particulate matter
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 125
- 230000005250 beta ray Effects 0.000 title claims abstract description 113
- 239000013618 particulate matter Substances 0.000 claims abstract description 82
- 238000004364 calculation method Methods 0.000 claims abstract description 56
- 238000005259 measurement Methods 0.000 claims abstract description 48
- 239000002245 particle Substances 0.000 claims abstract description 39
- 239000000523 sample Substances 0.000 claims abstract description 20
- 238000009423 ventilation Methods 0.000 claims abstract description 8
- 230000008859 change Effects 0.000 claims description 30
- 230000005855 radiation Effects 0.000 claims description 21
- 238000012360 testing method Methods 0.000 claims description 18
- 238000000691 measurement method Methods 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 7
- 230000003068 static effect Effects 0.000 claims description 2
- 238000006467 substitution reaction Methods 0.000 claims description 2
- 239000000428 dust Substances 0.000 abstract description 16
- 239000003344 environmental pollutant Substances 0.000 abstract description 9
- 231100000719 pollutant Toxicity 0.000 abstract description 9
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 239000007789 gas Substances 0.000 description 178
- 239000000126 substance Substances 0.000 description 48
- 239000007787 solid Substances 0.000 description 25
- 238000002474 experimental method Methods 0.000 description 20
- 239000012528 membrane Substances 0.000 description 14
- 239000000779 smoke Substances 0.000 description 13
- 238000010561 standard procedure Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 6
- 230000004044 response Effects 0.000 description 5
- 238000005303 weighing Methods 0.000 description 5
- 230000000052 comparative effect Effects 0.000 description 4
- 238000011065 in-situ storage Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000001678 irradiating effect Effects 0.000 description 3
- 230000009347 mechanical transmission Effects 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000011282 treatment Methods 0.000 description 3
- 238000009795 derivation Methods 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000002912 waste gas Substances 0.000 description 2
- 238000012951 Remeasurement Methods 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000008277 atmospheric particulate matter Substances 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 235000019504 cigarettes Nutrition 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 239000000356 contaminant Substances 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000008016 vaporization Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
Description
技术领域
本发明属于气体中颗粒物浓度测算技术领域,涉及β射线质量衰减理论,直接应用于工况条件下的颗粒物浓度检测,具体地说是基于β射线法的气体中颗粒物浓度直读测算方法及装置。
背景技术
目前,用于测量气体中颗粒物浓度的方法主要有以下三种:
一、重量法。重量法测算气体中颗粒物浓度时,利用采样管抽取一定体积的含颗粒物气体,用滤筒吸附颗粒物后,再将含颗粒物的滤筒拿到实验室中,经恒温恒湿处理后,用高精度天平进行手工称重,计算出颗粒物重量,由颗粒物重量除以抽气体积得出颗粒物浓度。我国于1996年颁布国家标准GB/T 16157《固定污染源排气中颗粒物测定与气态污染物采样方法》,将这种办法作为法定通行办法在全国施行。为解决颗粒物的超低排放问题,环境保护部于2017年颁布国家环境保护标准HJ 836《固定污染源废气低浓度颗粒物的测定重量法》,将滤筒改为滤膜,增大了抽气体积,用电子天平手工称重,依旧用重量法来计算气体中的颗粒物浓度。
在工况条件下利用重量法测颗粒物浓度时,存有以下不足:其一是测算用时长;由于现场采样后需要到实验室进行恒温恒湿处理,干燥、称重处理,通常需要3-5天才可获得测算结果,造成测算用时较长。其二是具有手工操作误差;由于该种测量方式中需要操作者手工操作,容易产生人为误差,影响测算结果。其三是测算数据精准度低;由于使用滤膜过滤颗粒物,滤膜材质不稳定时影响颗粒物重量的计算,从而影响测算数据的精准度。
二、β射线滤膜法。该种方法通过β射线滤膜法颗粒物测定仪器抽取一定体积的气体,在抽气过程中,用滤膜或者滤带吸附气体中的固态颗粒物,由测定仪器中的机械传送部件完成附尘前和附尘后的滤膜的机械转换,使β射线前后两次分别照射附尘前和附尘后的滤膜或者滤带,根据射线量的变化,计算出颗粒物的质量,除以所抽取气体的体积,得到颗粒物浓度。
该种方法相对重量法虽然在一定程度上解决了测算时间长、手工操作容易出现误差的问题,但仍存有以下不足:其一是测算时间较长;由于需要通过气泵抽取一定体积的气体,抽气过程耗费时间在30分钟以上,测算时间相对较长。其二是测算精准度低;由于该种方法是在测定仪器内对附尘前和附尘后的射线量进行测量,因而需要利用机械传送部件完成附尘前和附尘后滤膜的转换,使得该方法的测定仪器机械结构较为复杂,机械操作步骤较多,容易产生机械误差,影响测算结果的精准度。其三是该种方法同样使用滤膜过滤颗粒物,滤膜材质不稳定时影响颗粒物重量的计算,造成计算结果不准确。
三、β射线电离室法。2016年河北省发布地方标准DB13/T 2376《固定污染源废气低浓度颗粒物的测定β射线法》,该方法把采用电离室结构的β射线传感器放入烟道内,其原理可参考说明书附图1所示,用β射线源发出的起始射线量为I0的β射线直接照射被测气体,β射线在通过物质时会与物质中的原子或原子核相互作用,产生能量衰减后射线量为I,通过式6计算物质的质量,再通过物质的质量和电离室的体积计算气体中颗粒物的浓度C,则C的计算公式为式7:
该种方法虽然在一定程度上解决了测定时间长、容易出现机械以及人工误差的问题,但由于该种方法认为β射线在气态物质中不产生衰减,没有考虑气态物质的质量衰减对计算数据的影响,认为该计算公式可直接获得气体中的颗粒物浓度。但实际上,β射线不仅在固态物质中产生衰减,在气态物质中也产生衰减,按照该标准提供的式7计算出的是被测气体的气态物质和固态颗粒物的混合浓度,即C=C气+C固,而不是气体中的颗粒物浓度,致使计算结果与重量法比对偏差较大,影响测量结果的精准性。
发明内容
为解决现有技术中存在的以上不足,本发明旨在提供基于β射线法的气体中颗粒物浓度直读测算方法及装置,以达到缩短颗粒物浓度测算时间、提高测算结果精准性、直接应用于颗粒物工况检测的目的。
为实现上述目的,本发明所采用的技术方案如下:基于β射线法的气体中颗粒物浓度直读测算方法,设不含颗粒物气体的射线量I1、含颗粒物气体的射线量I2、颗粒物的质量衰减系数μm及射线辐照气体体积V为已知,则气体中颗粒物浓度C为:
式中,C的单位为mg/m3,V的单位为m3,μm的单位为mg-1,I1与I2的单位一致。
作为本发明的限定,含颗粒物气体的射线量I2通过测量获值,并通过计算获得I1;确定μm值,计算V,最终由式1计算出C。
作为本发明的进一步限定,该方法包括以下步骤:
步骤一:由β射线颗粒物测定仪直接测量含颗粒物气体的射线量I2,同时测得被测气体的压力P、温度T;根据P、T,计算I1;
步骤二:确定μm,计算V;
步骤三:根据式1计算C。
作为本发明的另一种限定,实时动态测量时,不含颗粒物气体的射线量I1动态、含颗粒物气体的射线量I2、颗粒物的质量衰减系数μm、β射线照射气体的测试时长t、气流通过辐照区域的横截面积S、气体的流动速度v、β射线的辐照体积V为已知,则气体中颗粒物浓度c动态为:
式中,c动态的单位为mg/m3,μm的单位为mg-1,v的单位为m/s,V的单位为m3,t的单位为s,S的单位为m2,I1动态与I2的单位一致。
作为本发明的进一步限定,设β射线穿过不含颗粒物气体辐射量的静态测量数值I1静,则I1动态为:
式中,I1动态与I1静的单位一致。
作为本发明的限定,计算I1静的步骤包括:
A.在标准实验条件下进行测量和计算
A1.在设定的温度T定和压力P定的条件下,测量出不含颗粒物气体的射线量I定,测算出β射线的辐照体积V和气流通过辐照区域的横截面积S;
A2.计算出β射线照射不含颗粒物气体在温度每变化1摄氏度溢出或溢入辐照体积V的距离LT以及压力每变化1Pa溢出或溢入辐照体积V的距离Lp;
B.在工况条件下实地测量
实地测量被测气体的温度T况和压力P况的数值;
C.计算I1静:
作为本发明的进一步限定,LT为不同温度T'条件下LT'数值的平均数,Lp为不同压力p'条件下Lp'数值的平均数,LT'、Lp'的计算公式为:
式中,I1″在实验室测量若干组不同温度T'或不同压力p'条件下的β射线穿过不含颗粒物气体的射线量数值。
作为本发明的进一步限定,测试时长t的测算为:
测量出零流速条件下β射线穿过不含颗粒物气体的射线量I1″′;
测量若干组不同流速v条件下的β射线穿过不含颗粒物气体的射线量I1动态',代入公式:
求出不同流速v条件下各自t'的数值,其平均数即为β射线照射不含颗粒物气体的测试时长t。
本发明还提供了一种基于β射线法的气体中颗粒物浓度直读测算装置,应用于所述的基于β射线法的气体中颗粒物浓度直读测算方法中使用,其技术方案如下:包括壳体,壳体上开设有进气口及出气口,进气口与出气口之间的壳体内连通形成通气通道,通气通道内固装有β射线发射极以及固装有与β射线发射极对应设置的β射线接收极,通气通道内还固装有温度传感探头以及压力传感探头,温度传感探头、压力传感探头及β射线接收极分别与显示仪器的信号输入端相连。
作为本发明的其它限定,壳体上固装有气体流速传感器,壳体固装有增温恒温装置。
由于采用了上述的技术方案,本发明通过创新β射线质量衰减理论和改进计算方法,依据β射线不仅在固态物质中产生衰减,在气态物质中也产生衰减的理论,利用β射线同时照射气、固两种状态的物质,将射线量分步骤测量,分离出气态物质中固态颗粒物的射线衰减量,从而获得了气体中颗粒物浓度值,颠覆了传统的理论和计算方法。本发明与现有技术相比,所取得的有益效果是:
(1)本发明的基于β射线法的气体中颗粒物浓度直读测算方法在测量含颗粒物气体的射线量I2时,可同时测得被测气体的压力P、温度T,根据P、T,可直接计算I1。本发明日常维护简单,没有日常使用的耗材,大大降低生产和维护成本。相对β射线滤膜法所采用的具有机械传送部件的测定仪器,本发明可有效避免机械部件产生的机械误差。
(2)提高工作效率。本发明相对重量法及β射线滤膜法,无需进行抽取气体的操作,可大大缩短现场的测算时间。本发明通过测量现场的含颗粒物气体的射线量I2,并于现场直接计算后得出I1,再将已知的μm、V代入式1中,可快速得到颗粒物的浓度。本发明可缩短现场的测算时间和整体的测算周期,在现场同时测试多个点位,并分析测试结果。
(3)消除了β射线在气态物质中的衰减量对固态颗粒物浓度计算的影响。相对β射线电离室法中认为β射线在气态物质中不产生衰减,按照公式7计算出的是被测气体的气态物质和固态颗粒物的混合浓度,即C=C气+C固,而不是气体中的颗粒物浓度。本发明依据β射线不仅在固态物质中产生衰减,在气态物质中也产生衰减的理论,分离了在气态物质中颗粒物的射线衰减量,直接获得了气体中的颗粒物浓度的计算方法。
(4)提高测量精准度,保证数据真实性。本发明相对重量法及β射线滤膜法,减少了人工和机械的操作步骤,从而减少了人为和机械误差。计算方法中也消除了β射线在气态物质中的衰减量对固态颗粒物浓度计算的影响,具有较高的数据精准度;同时,省去了中间环节,消除了数据造假机会,有效保证了数据的真实性、可靠性和有效性。
(5)简化现场操作环节。本发明通过在实验室测出I定数值,在工况条件下测得I2的同时,也测量气体的温度、压力值,便于现场计算得到工况条件下的I1值,省掉了现场实测I1值的环节,方便使用本装置者现场测试操作。
(6)提供了便捷可行的校准方法。在对基于β射线法的气体中颗粒物浓度直读测算装置定期校准中,只需在设定温度T定和设定压力P定条件下,对I定数值进行重新测量和标定。不必如β射线滤膜法及β射线电离室法一样,不仅需要复杂的校准程序,还需要经常与手工称重法进行比对试验,来不断修正偏离正常值太多的检测数据。
(7)实现了实时动态测算颗粒物的浓度值。本发明能够迅速反映颗粒物浓度值的实时变化,实现了气体中颗粒物排放浓度的实时在线监测,便于迅速采取措施加强生产过程干预,减少污染物排放。
本发明适用于在工况条件下的烟道污染物排放、粉尘污染工作环境、大气颗粒物监测、在线污染物实时监测等方面测算气体中颗粒物的浓度。
附图说明
下面结合附图及具体实施例对本发明作更进一步详细说明。
图1为本发明背景技术利用β射线颗粒物测定仪测定颗粒物浓度的原理图;
图2为本发明实施例1的β射线穿过不含颗粒物的气体后的射线量I1原理图;
图3为本发明实施例1的β射线穿过含颗粒物的气体后的射线量I2原理图;
图4为本发明实施例3的结构示意图;
图5为本发明实施例6的结构示意图;
图6为实施例7中实验室测量若干组不同温度条件下的I1;
图7为实施例7中实验室测量若干组不同流速条件下的I1;
图8为实施例7中流动气体中的颗粒物质量衰减系数。
图中:1、β射线发射极;2、β射线接收极;3、温度传感探头;4、压力传感探头;5、进气口;6、出气口;7、壳体;8、气体流速传感器;9、增温恒温装置。
具体实施方式
以下结合附图对本发明的优选实施例进行说明。应当理解,此处所描述的优选实施例仅用于说明和理解本发明,并不用于限定本发明。
实施例1气体中颗粒物浓度计算公式的推导
本实施例提出了一种新的可直接计算气体中颗粒物浓度计算公式,即依据β射线不仅在固态物质中产生衰减、在气态物质中也产生衰减,利用β射线同时照射气、固两种状态的物质,将固-气混合物质中固态颗粒物的β射线衰减量分离出来后,推导出气体中颗粒物浓度C计算公式,具体推导过程如下:
(1)考虑了β射线在气态物质、固态物质中的衰减量
β射线不仅在固态物质中产生衰减,在气态物质中也产生衰减,β射线对这种含颗粒物气体物质的衰减,存在两种质量衰减系数,即:气态物质衰减系数和固态物质衰减系数。β射线在气、固两种状态物质中的质量衰减系数相差很大,成千上万倍。在气态物质中的质量衰减系数极小,一般可以忽略不计,认为不衰减。但在此种测量和计算中,气态物质相对固态物质的质量,数据差距极大。分析单位体积气体中的颗粒物含量(毫克每立方米),标况条件下,每立方米气态物质的质量在1.29千克以上,固态颗粒物却只有几或几十毫克,数据相差万倍以上。β射线在气态物质中的衰减量,足以影响颗粒物浓度的计算,不能忽略不计,必须予以考虑。
(2)将固-气混合物质中固态颗粒物的β射线衰减量分离
用β射线照射被测气体时,β射线同时照射了气态物质和悬浮其中的固态物质(颗粒物),这两种物质共用一个体积(由于颗粒物体积较小,忽略不计),混合在一起。β射线对这种混合物质的衰减,存在两个衰减量,即:在气态物质中的衰减量I气衰和在固态物质中的衰减量I固衰。测量这种固-气混合物质中颗粒物浓度时,需要对射线衰减量进行分离。通过首先测出β射线穿过不含颗粒物气体后的射线量I1,再测出β射线穿过含颗粒物气体后的射线量I2,即可分离出固态颗粒物的衰减量见式8,其对数关系式见式9:
I固衰=I1-I2--------------------------(式8),
(3)推导可直接计算颗粒物浓度的公式
把测量过程分成两个操作步骤来进行:第一步测量β射线穿过不含颗粒物的气体后的射线量I1,参考说明书附图2;第二步测量β射线穿过含颗粒物的气体后的射线量I2,参考说明书附图3。由公式9和公式7,推导出可直接计算颗粒物浓度的公式1:
式中,C——气体中颗粒物浓度,单位为mg/m3;
I1——不含颗粒物气体的射线量;
I2——含颗粒物气体的射线量,I1与I2的单位保持一致;
μm——颗粒物的质量衰减系数,单位为mg-1;
V——射线辐照气体体积,单位为m3。
通过测量含颗粒物气体的射线量I2,计算不含颗粒物气体的射线量I1、确定颗粒物的质量衰减系数μm及计算射线辐照气体体积V,由式1计算出C。
实施例2基于β射线法的气体中颗粒物浓度直读测算方法及装置
本实施例包括以下步骤,其中步骤二也可先于步骤一进行:
步骤一:由基于β射线法的气体中颗粒物浓度直读测算方法的装置于现场直接测量含颗粒物气体的射线量I2,同时测得被测气体的压力P、温度T。根据P、T,计算I1。I1、I2的单位为V。
如图4所示,基于β射线法的气体中颗粒物浓度直读测算方法装置包括壳体7,壳体7上开设有进气口5及出气口6,进气口5与出气口6之间的壳体7内连通形成通气通道,通气通道内固装有β射线发射极1、β射线接收极2、温度传感探头3以及压力传感探头4,保证β射线发射极1、β射线接收极2与温度传感探头3和压力传感探头4处于同一测量环境条件中。
β射线发射极1作为β射线源用于发射β射线,β射线接收极2用于接收β射线并检测出含颗粒物气体的射线量I2。β射线接收极2与β射线发射极1对应设置,即β射线发射极1与β射线接收极2于通气通道内对侧相对设置。压力传感探头4用于检测被测气体的压力P,温度传感探头3用于检测被测气体的温度T。温度传感探头3、压力传感探头4及β射线接收极2分别通过信号转化器与显示仪器的信号输入端相连,在显示仪器上能够显示出被测气体的温度T、被测气体的压力P以及含颗粒物气体的射线量I2。显示仪器可为显示器,也可为具有显示屏的计算机。
步骤二:确定μm。计算V,即β射线辐照气体的体积。因固体颗粒物占用的体积太小,可以忽略不计。
其中,确定μm与计算V二者之间的顺序并无特定关系,μm及V只需在进行步骤三前获得即可。
步骤三:根据式1计算C。
实施例3基于β射线法的气体中颗粒物浓度直读测算方法
本实施例详细说明计算不含颗粒物气体的射线量I1的方法,还包括以下两种方式:
第一种:现场计算法
在现场测量得到压力P、温度T后,直接通过现场计算得到I1值。具体包括以下步骤:
S11.在标准实验条件下,并在设定温度T定和设定压力P定条件下,测量出不含颗粒物气体的射线量I定,还可以针对不同检测环境下的被测气体,在烟道污染物排放、粉尘污染工作环境、大气颗粒物监测、在线污染物实时监测等方面测定各自设定温度、压力等条件下的I定数值。例如:在烟道污染物排放方面设定温度为140℃、压力为常压;在粉尘污染工作环境方面设定温度为25℃、压力为常压。
S12.在标准实验条件下:在压力不变的情况下升温N度,测量升温后的不含颗粒物气体的射线量I升温;在温度不变的情况下增加N个单位大气压力,测量升压后的不含颗粒物气体的射线量I升压。
S13.计算出不含颗粒物气体射线量的温度变化系数I温系数和不含颗粒物气体射线量的压力变化系数I压系数,其中:
I温系数=(I升温-I定)÷N------------------------(式2);
I压系数=(I定-I升压)÷N—----------------------(式3);
S14.计算I1:
I1=I定+(T-T定)×I温系数+(P-P定)×I压系数-----------(式4)。
此种方式为两种方式中的优选方式,此种方式的前期计算量较小,且具有较高的精准度,便于现场直接计算获得I1值。
第二种:预先储存法
在进行气体中颗粒物浓度测算前,预先于实验室中测算出不同温度与不同压力下对应的不含颗粒物气体的射线量数值I1中。在现场测得的压力P、温度T,由预先测定的多组I1中找到与该压力P、温度相对应的I1值。此种方式能够在一定程度上减少现场测算工作量,但前期所需储存的数据较多,工作量较大。
实施例4基于β射线法的气体中颗粒物浓度直读测算方法
本实施例详细说明实施例2的步骤二中在洁净环境的实验室里计算μm的步骤,具体包括:
S21.在洁净环境的实验室里,在常温、常压、正常湿度条件下,测量出不含颗粒物气体的射线量I1'。
S22.在洁净环境的实验室里,在常温、常压、正常湿度条件下,在步骤S21的不含颗粒物气体中加入一定浓度的颗粒物,测量出气体加颗粒物后的射线量I2'。
S23.用重量法按照标准GB/T 16157-1996或标准HJ 836-2017计算出气体中的颗粒物浓度C'。
S24.测算射线辐照气体体积V',由以下式5计算得到μm:
实施例5基于β射线法的气体中颗粒物浓度直读测算方法
本实施例为实施例3中现场计算法的校准步骤,步骤S11前包括在基于β射线法的气体中颗粒物浓度直读测算装置使用过程中,定期进行的校准步骤:
在每年定期校准中,只需在设定温度T定和设定压力P定条件下,对I定数值进行重新测量和标定。标物为设定温度T定和设定压力P定条件下的不含颗粒物气体,其质量是稳定的,β射线穿过其后的I定数值变化反映的就是β射线发射极和β射线接收极在装置使用一段时间后产生的信号失真变化,重新测量和标定I定数值,等于重新校准了仪器的β射线发射极到接收极的射线量,方法科学,简单可行。温度和压力的校准,按照国家相应的规程进行即可。
实施例6一种基于β射线法的气体中颗粒物浓度直读测算装置及方法
由于β射线在气态物质中的质量衰减系数的数值本身极小而且是个基本不变的常量,在气态物质中进行比较时,气体中各组成物质的成分和比例对I1数值的影响就微乎其微,可以忽略不计。只有气体的温度、压力变化时,对其密度产生的影响较大,能够明显影响其质量的变化,进而对I1的数值产生影响。
本实施例为动态实时测算流动气体I1的测算,利用基于β射线法的气体中颗粒物浓度直读测算装置来实现,本实施例中基于β射线法的气体中颗粒物浓度直读测算装置是在实施例2中基于β射线法的气体中颗粒物浓度直读测算装置的基础上,壳体上固装有气体流速传感器8,壳体外表面固装有增温恒温装置9,增温恒温装置采用现有技术中能够使液态水汽化的结构皆可,如电加热丝,以保证通过检测的气流恒定保持在100度以上。
实时动态测量时,基于β射线法的气体中颗粒物浓度直读测算方法,包括以下步骤:
A.在标准实验条件下进行测量和计算
A1.在设定的温度T定和压力P定的条件下,测量出不含颗粒物气体的射线量I定,测算出β射线的辐照体积V和气流通过辐照区域的横截面积S。
A2.计算出β射线照射不含颗粒物气体在温度每变化1摄氏度溢出或溢入辐照体积V的距离LT(亦为温度与体积的变化系数)以及压力每变化1Pa溢出或溢入辐照体积V的距离Lp(亦为压力与体积的变化系数)。具体地:
LT为不同温度T'条件下LT'数值的平均数,Lp为不同压力p'条件下Lp'数值的平均数,LT'、Lp'的计算公式为:
式中,I1″在实验室测量若干组不同温度T'或不同压力p'条件下的β射线穿过不含颗粒物气体的射线量数值。
B.在工况条件下实地测量
实地测量被测气体的温度T况和压力P况的数值。
C.计算β射线穿过不含颗粒物气体辐射量的I1静态测量数值I1静
D.测算I1动态
D1.在实验室的标准实验条件下,测算出β射线辐照气体的测试时长t(亦为β射线接收极每次测量的响应时间)。具体测算方法是:在实验室先测量出零流速条件下β射线穿过不含颗粒物气体的射线量I1″′,再测量若干组不同流速v条件下的β射线穿过不含颗粒物气体的射线量I1动态',代入公式求出不同流速v条件下各自t'的数值,其平均数即为β射线照射的测试时长t。
D2.现场测量出气体的流动速度v。
D3.计算实时动态测量时不含颗粒物气体的射线量I1动态:
式中,v的单位为m/s,V的单位为m3,t的单位为s,S的单位为m2,I1动态与I1静的单位一致。
E.在流动气体中测量颗粒物质量衰减系数值时,具体包括:
现场测量出流动气体的流速v',β射线的辐照体积V'和气流通过辐照区域的横截面积s',β射线照射气体的测试时长t'后,计算出I1静值,按照公式91计算出I1动态'值。用重量法按照标准GB/T 16157-1996或标准HJ 836-2017计算出流动气体中的颗粒物浓度c动态'值,现场直接测量含颗粒物气体的射线量I2',工况条件下流动气体中的颗粒物质量衰减系数值可以按照下列公式计算:
F.现场直接测量含颗粒物气体的射线量I2,计算实时动态测量时气体中颗粒物浓度c动态为:
式中,c动态的单位为mg/m3,μm的单位为mg-1,v的单位为m/s,V的单位为m3,t的单位为s,S的单位为m2,I1动态与I2的单位一致。
本实施例的技术原理是:
一定量的气体,温度的变化可使体积产生膨胀或缩小的变化,其体积变化与显示温度高低刻度的变化之间存在一个固定系数,参照水银温度计测量温度变化引起刻度升降变化的原理,温度每升降1度引起刻度升降一个固定的数值。反映到气体在β射线辐照区域V内随温度的变化而膨胀或缩小,在溢出或溢入辐照区域的变化就是:气体温度T每变化1度溢出或溢入辐照体积V的距离LT是个固定数,其溢出或溢入的体积数值是s×LT,体积变化系数就是气体随温度变化溢出或溢入的体积s×(T-T定)×LT。在其密度随温度变化后的体积中的占比,即:根据同一密度的物质,体积和质量成正比,其体积的比值可以看作是质量的比值。β射线穿过物质时产生衰减,β射线量的增减与物质的质量增减相关,气体中β射线量的温度变化系数就是:在I定基础上,气体的体积变化系数与I定的乘积,即:
压力变化原理与此一致,气体压力P每变化1帕溢出或溢入辐照体积V的距离LP是个固定数,其溢出或溢入的体积数值是s×LP,体积变化系数就是气体随压力变化溢出或溢入的体积s×(P定-P况)×LP,在其密度随压力变化后的体积中的占比,即:根据同一密度的物质,体积和质量成正比,体积的比值可以看作是质量的比值。β射线穿过物质时产生衰减,β射线量的增减与物质的质量增减相关,气体中β射线量的压力变化系数就是:在I定基础上,气体的体积变化系数与I定的乘积,即:
β射线接收极接收β射线信号不是一般意义上理解的像光速一样瞬时接收,而是受制造接收极材料的特性影响,有一段接收时间,亦即测试时长或响应时间。在这段时间内不同流速的气体进入β射线辐照区域内的体积是不同的,而同一种材料接收极的响应时间t是固定的,与气体流速v有关的体积数值是s×v×t,体积变化系数就是气体随流速变化增加的体积s×v×t,在其变化后的体积V+s×v×t中的占比,即:根据同一密度的物质,体积和质量成正比,体积的比值可以看作是质量的比值,β射线穿过物质时产生衰减,β射线量的增减与物质的质量增减相关,在I1静数值基础上,气体中β射线量的流速变化系数就是气体的体积变化系数与I1静的乘积,计算流动气体中β射线量的公式为:
实施例7实验例
利用实施例6中装置及方法进行实验。
(1)用公式93计算LT
在实验室温度25℃的室温条件下,在9组不同增温恒温装置9温度条件下,测量并记录样机的I1″数值,设141℃时为I定值,再按照公式93计算其它8组的LT'值,计算其平均数。使用一台实施例6中装置进行9次测算,具体计算过程数据见下表1:
表1
LT'取平均值后,得LT为0.0000150674米/度,分析其相对误差:
按公式93计算不同温度下LT'的数值都很接近,根据LT'计算出β射线I1的计算数值,与实测β射线I1的误差最大为-0.256391577%(见下表2和说明书附图6),相对误差极低。考虑实测时极易产生测量误差,取平均数的计算值应该更准确可信。
(2)利用公式95计算t
实验初始条件:在实验室温度25度的室温条件下,将实施例6中装置到风洞实验装置里,先记录零风速下I1实测值,再在不同风速条件下测量并记录样机的I1动态数值。做9组实验,按公式分别计算t'的数值,取其平均数。
使用一台实施例6中装置进行10次测算,按照公式计算t数值是0.001257553秒,具体计算过程数据见下表3:
表3
分析其相对误差:
按公式计算t'的数值都很接近,说明同一材料β射线接收极的响应时间是一致的。根据响应时间t'的平均值计算的I1动态数值与实测值的误差最大值是4.855673013%(见上表3和说明书附图7),考虑实测时极易产生测量误差,取平均数的计算值应该更准确可信。
(3)利用公式52计算μm
实验设备:实施例6中装置、风洞实验装置、崂应3012H型自动烟尘(气)测试仪、梅特勒-托利多仪器(上海)有限公司生产的十万分之一电子天平。
实验初始条件:在实验室温度25度的室温条件下,将样机的测量传感部分放置到风洞实验装置里,将风洞里的气体温度稳定在43.0度,气体流速稳定在10.8米每秒(不能完全稳定,上下有浮动),先记录无颗粒物尘的条件下I1动态实测值,再向风洞里加不同浓度的颗粒物尘,按国标法测量颗粒物的浓度c动态。做6组不同浓度条件下的实验,测量并记录样机的I2数值。由上例已知t的平均数是0.001257553秒,按公式计算各次的质量衰减系数值。
实验过程:在一定的流速下,进行多组实验,测量已知不同浓度颗粒物的气体,用重量法测知气体的颗粒物浓度。同时用本方法测量其不同浓度下的I2数值,按照公式计算各自的μm的数值,以多组数据的平均值来确定μm的数值,按照公式计算出这台样机的μm数值是726.0820426毫克-1,具体计算过程数据见下表4:
表4
分析相对误差:
因为每次的测量都是以国标法(重量法)来测量颗粒物的浓度c动态,以此计算出样机测量颗粒物μm的6组数值都很接近,说明其相对误差极小,而且始终与国标法保持一致(见上表4和说明书附图8),其平均值就更接近颗粒物质量衰减系数的真实数值。
对比实验例8在实验室中采用实施例2与其它测量方法的比对实验
实验温度:23℃
实验湿度:50%RH
实验设施:30m3环境实验舱、恒温恒湿室
污染物:红塔山牌,焦油量8mg产生颗粒物70~80mg/根
主要实验设备:BDY-315型电离室法颗粒物浓度检测仪,实施例2中装置、TSI粉尘测试仪、崂应3012H型自动烟尘(气)测试仪、十万分之一电子天平。
实验设计:在23℃,50%RH恒温恒湿条件下,向30m3环境实验舱中打入不少于20根香烟,同时开动扰动风扇,10min后同时用BDY-315型电离室法颗粒物浓度检测仪(第一组)、实施例2中装置(第二组)、TSI粉尘仪(第三组)、重量法(第四组)来分别测试颗粒物浓度值。下表为各自的实验测试数据:
分析以上数据:
用本方法测量的平均数据(第二组)最接近国标法重量法(第四组)测试颗粒物浓度的数值,而且不同时间段的检测数据,与重量法数值偏离不多。其它仪器的检测数据,无论平均值还是单个不同时段的数值,都与重量法的数据偏差较多。说明实施例2中装置测量的数据与重量法一致性最高。
对比实验例9在工况条件下采用实施例6与按国家标准重量法的比对实验
实验1
实验地点:某电厂锅炉烟气排放出口
烟温:39℃
水分含量:33%
主要实验设备:实施例6中装置、崂应3012H型自动烟尘(气)测试仪、梅特勒-托利多仪器(上海)有限公司生产的电子天平。
实验设计:实施例6中装置和崂应国标法烟尘采样器同时放入排放出口;实施例6中动态实时显示数据,每2分钟读一次平均数据;崂应国标法烟尘采样器采样30分钟,并用电子天平进行实验室称重测量。
排放出口实验数据:
实验2
实验地点:某钢厂烟气排放出口
烟温:68℃
水分含量:9%
主要实验设备:实施例6中装置、崂应3012H型自动烟尘(气)测试仪、梅特勒-托利多仪器(上海)有限公司生产的电子天平。
实验设计:将实施例6中装置和崂应国标法烟尘采样器同时放入排放出口;实施例6中动态实时显示数据,每2分钟读一次平均数据;崂应国标法烟尘采样器采样30分钟,并进行实验室称重测量。
排放出口实验数据:
对以上两次工况实验的数据结果分析,可以得知:
实施例6中测试数据与崂应国标法烟尘采样器的重量法测试数据偏差很小,在误差允许范围内,与重量法基本保持一致。说明在工况条件下按本方法测量颗粒物浓度与按国家标准方法测量的数据高度一致。
需要说明的是,以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照上述实施例对本发明进行了详细的说明,对于本领域技术人员来说,其依然可以对上述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
2.根据权利要求1所述的基于β射线法的气体中颗粒物浓度直读测算方法,其特征在于:含颗粒物气体的射线量I2通过测量获值,并通过计算获得I1;确定μm值,计算V,最终由式1计算出C。
3.根据权利要求2所述的基于β射线法的气体中颗粒物浓度直读测算方法,其特征在于该方法包括以下步骤:
步骤一:由β射线法颗粒物测定仪直接测量含颗粒物气体的射线量I2,同时测得被测气体的压力P、温度T;根据P、T,计算I1;
步骤二:确定μm,计算V;
步骤三:根据式1计算C。
9.一种基于β射线法的气体中颗粒物浓度直读测算装置,其特征在于:应用于权利要求1-8中所述的基于β射线法的气体中颗粒物浓度直读测算方法中使用,包括壳体,壳体上开设有进气口及出气口,进气口与出气口之间的壳体内连通形成通气通道,通气通道内固装有β射线发射极以及固装有与β射线发射极对应设置的β射线接收极,壳体上还固装有温度传感探头以及压力传感探头,温度传感探头、压力传感探头及β射线接收极分别与显示仪器的信号输入端相连。
10.根据权利要求9所述的基于β射线法的气体中颗粒物浓度直读测算装置,其特征在于:壳体上固装有气体流速传感器,壳体加装有增温恒温装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2021103002649 | 2021-03-22 | ||
CN202110300264 | 2021-03-22 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113670783A true CN113670783A (zh) | 2021-11-19 |
Family
ID=78549504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111079440.7A Withdrawn CN113670783A (zh) | 2021-03-22 | 2021-09-15 | 基于β射线法的气体中颗粒物浓度直读测算方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113670783A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114186838A (zh) * | 2021-12-08 | 2022-03-15 | 河北建投能源科学技术研究院有限公司 | 一种对气固两相流测量装置进行精确度评价的设备和方法 |
CN115372220A (zh) * | 2022-10-26 | 2022-11-22 | 杭州泽天春来科技有限公司 | 颗粒物浓度检测装置 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002357532A (ja) * | 2001-06-01 | 2002-12-13 | Shimadzu Corp | 浮遊粒子状物質測定装置 |
JP2003139725A (ja) * | 2001-10-31 | 2003-05-14 | Shimadzu Corp | 浮遊粒子状物質測定装置 |
US20090321635A1 (en) * | 2007-07-03 | 2009-12-31 | Wuhan Tianhong Instruments Co., Ltd. | Beta-ray soot concentration direct readout monitor and method for determining effective sample |
CN202814896U (zh) * | 2011-06-23 | 2013-03-20 | 北京地海云天环境科技有限公司 | 低能β射线快速烟尘检测仪 |
CN204789249U (zh) * | 2015-07-15 | 2015-11-18 | 中国工程物理研究院材料研究所 | 一种高准确度β射线法在线测量大气颗粒物浓度的装置 |
CN105334147A (zh) * | 2015-12-04 | 2016-02-17 | 深圳睿境环保科技有限公司 | 基于β射线法和光散射法的颗粒物在线监测系统及方法 |
CN207336493U (zh) * | 2017-11-02 | 2018-05-08 | 中广核久源(成都)科技有限公司 | 一种抽真空测量放射性气溶胶装置 |
CN108152178A (zh) * | 2016-12-05 | 2018-06-12 | 天津普瑞爱尔科技有限公司 | β射线颗粒物浓度检测传感器及其检测方法 |
CN108195854A (zh) * | 2017-12-26 | 2018-06-22 | 中国计量科学研究院 | 一种x射线空气衰减系数检测方法 |
CN109507074A (zh) * | 2019-01-21 | 2019-03-22 | 上海理工大学 | 一种超低排放烟尘浓度监测装置和监测方法 |
CN109540754A (zh) * | 2018-12-28 | 2019-03-29 | 安徽蓝盾光电子股份有限公司 | 一种基于β射线法的大气颗粒物在线监测装置及方法 |
-
2021
- 2021-09-15 CN CN202111079440.7A patent/CN113670783A/zh not_active Withdrawn
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002357532A (ja) * | 2001-06-01 | 2002-12-13 | Shimadzu Corp | 浮遊粒子状物質測定装置 |
JP2003139725A (ja) * | 2001-10-31 | 2003-05-14 | Shimadzu Corp | 浮遊粒子状物質測定装置 |
US20090321635A1 (en) * | 2007-07-03 | 2009-12-31 | Wuhan Tianhong Instruments Co., Ltd. | Beta-ray soot concentration direct readout monitor and method for determining effective sample |
CN202814896U (zh) * | 2011-06-23 | 2013-03-20 | 北京地海云天环境科技有限公司 | 低能β射线快速烟尘检测仪 |
CN204789249U (zh) * | 2015-07-15 | 2015-11-18 | 中国工程物理研究院材料研究所 | 一种高准确度β射线法在线测量大气颗粒物浓度的装置 |
CN105334147A (zh) * | 2015-12-04 | 2016-02-17 | 深圳睿境环保科技有限公司 | 基于β射线法和光散射法的颗粒物在线监测系统及方法 |
CN108152178A (zh) * | 2016-12-05 | 2018-06-12 | 天津普瑞爱尔科技有限公司 | β射线颗粒物浓度检测传感器及其检测方法 |
CN207336493U (zh) * | 2017-11-02 | 2018-05-08 | 中广核久源(成都)科技有限公司 | 一种抽真空测量放射性气溶胶装置 |
CN108195854A (zh) * | 2017-12-26 | 2018-06-22 | 中国计量科学研究院 | 一种x射线空气衰减系数检测方法 |
CN109540754A (zh) * | 2018-12-28 | 2019-03-29 | 安徽蓝盾光电子股份有限公司 | 一种基于β射线法的大气颗粒物在线监测装置及方法 |
CN109507074A (zh) * | 2019-01-21 | 2019-03-22 | 上海理工大学 | 一种超低排放烟尘浓度监测装置和监测方法 |
Non-Patent Citations (2)
Title |
---|
KATSUYUKI TAKAHASHI ET AL.: ""Examination of discrepancies between beta-attenuation and gravimetric methods for the monitoring of particulate matter"", 《ATMOSPHERIC ENVIRONMENT》 * |
李德文 等: ""基于β 射线法的粉尘质量浓度检测算法研究"", 《矿业安全与环保》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114186838A (zh) * | 2021-12-08 | 2022-03-15 | 河北建投能源科学技术研究院有限公司 | 一种对气固两相流测量装置进行精确度评价的设备和方法 |
CN115372220A (zh) * | 2022-10-26 | 2022-11-22 | 杭州泽天春来科技有限公司 | 颗粒物浓度检测装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Meyer et al. | Two decades of water vapor measurements with the FISH fluorescence hygrometer: a review | |
CN113670783A (zh) | 基于β射线法的气体中颗粒物浓度直读测算方法及装置 | |
KR101408513B1 (ko) | 베타선 매연 농도 직접 판독 모니터 및 유효 샘플 결정법 | |
CN105334147B (zh) | 基于β射线法和光散射法的颗粒物在线监测系统及方法 | |
Hinz et al. | A comprehensive experimental study of aerial pollutants in and emissions from livestock buildings. Part 1: Methods | |
US20110197649A1 (en) | Self-calibrating gas sensor | |
CN109870395B (zh) | 基于滤膜自动称重测量空气中颗粒物浓度的不确定度分析方法 | |
US10495573B1 (en) | Instrument for measuring airborne particulate matter | |
MARPLE et al. | An aerosol chamber for instrument evaluation and calibration | |
Wedding et al. | An automatic particle sampler with beta gauging | |
Speer et al. | An instrument for measuring the liquid water content of aerosols | |
US3700333A (en) | Method and apparatus for making an in-situ determination of the chemical properties of atmospheric aerosols | |
CN105954789B (zh) | 一种氡子体测量仪α能谱峰重叠修正因子的准确刻度方法 | |
CN106707325A (zh) | 用于测量一定时间段内介质表面氡释放总量的方法与系统 | |
L'Orange et al. | A high-throughput, robotic system for analysis of aerosol sampling filters | |
CN109814144A (zh) | 一种氡和氡子体同步测量及剂量率评价方法 | |
CN109655870B (zh) | 流气式氡源任意调节氡析出率及有效衰变常数的装置及方法 | |
CN104406932B (zh) | 固定污染源废气二氧化硫的紫外吸收测定方法 | |
CN216747021U (zh) | 固定污染源废气颗粒物直接测定所用的等速采样系统 | |
CN109253953B (zh) | 一种可测量多种颗粒物质量浓度的方法和系统 | |
CN106018401A (zh) | 一种水中氯离子含量的测定方法 | |
Robson et al. | Evaluation of air particulate sampling equipment | |
CN109752751B (zh) | 空气中131i动态平衡监测装置及监测方法 | |
RU2330308C1 (ru) | Способ контроля массовой доли изотопа уран-235 в газовой фазе гексафторида урана и система измерения для его реализации | |
CN110632642A (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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20211119 |