CN106197544B - 颗粒堆积型射气介质氡迁移物理参数的确定方法及其测量装置 - Google Patents
颗粒堆积型射气介质氡迁移物理参数的确定方法及其测量装置 Download PDFInfo
- Publication number
- CN106197544B CN106197544B CN201610565316.4A CN201610565316A CN106197544B CN 106197544 B CN106197544 B CN 106197544B CN 201610565316 A CN201610565316 A CN 201610565316A CN 106197544 B CN106197544 B CN 106197544B
- Authority
- CN
- China
- Prior art keywords
- mrow
- radon
- msub
- mfrac
- sample
- 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
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 title claims abstract description 275
- 229910052704 radon Inorganic materials 0.000 title claims abstract description 270
- 239000002245 particle Substances 0.000 title claims abstract description 34
- 238000012856 packing Methods 0.000 title claims abstract description 31
- 230000005012 migration Effects 0.000 title claims abstract description 28
- 238000013508 migration Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000009792 diffusion process Methods 0.000 claims abstract description 38
- 238000004519 manufacturing process Methods 0.000 claims abstract description 22
- 238000005259 measurement Methods 0.000 claims abstract description 16
- 230000035699 permeability Effects 0.000 claims abstract description 15
- 239000011148 porous material Substances 0.000 claims abstract description 8
- 238000007789 sealing Methods 0.000 claims abstract description 3
- 239000000523 sample Substances 0.000 claims description 100
- 238000011068 loading method Methods 0.000 claims description 23
- 238000012544 monitoring process Methods 0.000 claims description 19
- 239000006101 laboratory sample Substances 0.000 claims description 18
- 230000006641 stabilisation Effects 0.000 claims description 11
- 238000011105 stabilization Methods 0.000 claims description 11
- 238000001035 drying Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 9
- 238000005325 percolation Methods 0.000 claims description 7
- 238000002474 experimental method Methods 0.000 claims description 5
- 238000002360 preparation method Methods 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 5
- 238000009530 blood pressure measurement Methods 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 230000007480 spreading Effects 0.000 claims description 3
- 238000003892 spreading Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 210000005239 tubule Anatomy 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 230000001133 acceleration Effects 0.000 claims description 2
- 238000007664 blowing Methods 0.000 claims description 2
- 239000012530 fluid Substances 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000009423 ventilation Methods 0.000 claims description 2
- 239000000243 solution Substances 0.000 claims 3
- 238000010521 absorption reaction Methods 0.000 claims 1
- 238000005422 blasting Methods 0.000 claims 1
- 239000003610 charcoal Substances 0.000 claims 1
- 238000010790 dilution Methods 0.000 claims 1
- 239000012895 dilution Substances 0.000 claims 1
- 238000012423 maintenance Methods 0.000 abstract 1
- 238000001556 precipitation Methods 0.000 description 5
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 229910052770 Uranium Inorganic materials 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 206010058467 Lung neoplasm malignant Diseases 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 201000005202 lung cancer Diseases 0.000 description 2
- 208000020816 lung neoplasm Diseases 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 229910052705 radium Inorganic materials 0.000 description 2
- HCWPIIXVSYCSAN-UHFFFAOYSA-N radium atom Chemical compound [Ra] HCWPIIXVSYCSAN-UHFFFAOYSA-N 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 241000282376 Panthera tigris Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000711 cancerogenic effect Effects 0.000 description 1
- 231100000357 carcinogen Toxicity 0.000 description 1
- 239000003183 carcinogenic agent Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 235000019504 cigarettes Nutrition 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007865 diluting Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000002345 respiratory system Anatomy 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 230000031068 symbiosis, encompassing mutualism through parasitism Effects 0.000 description 1
- 229910000442 triuranium octoxide Inorganic materials 0.000 description 1
- 239000010863 uranium mill tailing Substances 0.000 description 1
- 229910001727 uranium mineral Inorganic materials 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
-
- 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/08—Investigating permeability, pore-volume, or surface area of porous materials
-
- 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/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明涉及一种颗粒堆积型射气介质氡迁移物理参数的确定方法及测量装置，可同时确定此类介质的氡扩散系数、可运移氡产生率、介质渗透率以及孔隙率。测量采用本发明提供的实验装置，确定纯扩散平衡状态下介质底部孔隙氡浓度和表面环境空气的氡浓度、两种密闭循环体积下不同时刻循环气体的氡浓度值和介质的气体压力梯度；然后利用得到多孔介质的孔隙率ε、单位体积孔隙介质产生的运移氡量α、多孔介质的渗透率k以及扩散系数D。本发明方法可同时确定多个物理参数，装置结构简单、制造和维修成本低、易操作。
Description
技术领域
本发明属于多孔射气介质氡迁移模拟与控制领域，具体涉及一种砂壤、废石堆、矿石堆和尾矿堆等颗粒堆积型射气介质氡迁移的物理参数的确定方法及测量装置，特别是对砂壤、废石堆、矿石堆和尾矿堆等颗粒堆积型射气介质的氡扩散系数、可运移氡产生率、渗透率以及孔隙率的确定方法及测量装置。
背景技术
氡是放射性气体，当人们吸入后，氡衰变产生的α粒子将对呼吸***造成辐射损伤，诱发肺癌。氡是世界卫生组织(World Health Organization简称WHO)公布的19种主要致癌物质之一，是仅次于香烟导致人类肺癌的第二大元凶。
砂壤、废石堆、矿石堆和尾矿堆等含颗粒堆积射气介质是铀矿山大气氡的主要来源。这类介质产生的氡气将在浓度梯度和压力梯度的作用下沿着颗粒间连通的孔隙通道运移而析出。由于颗粒堆积型射气介质的氡扩散系数、可运移氡产生率、渗透率以及孔隙率受到颗粒粒径分布和堆积方式等的影响较大，依据渗流-扩散理论，采用室内实验的手段研究这类介质氡的迁移和析出规律，以及进行氡析出率估算时，为了研究介质内氡的迁移规律和估算表面氡析出率，必须确定介质的氡扩散系数、可运移氡产生率、渗透率以及孔隙率等四个物理参数，以确保实验获得的氡迁移和析出规律的准确性。
目前，这四个物理参数的测定方法较复杂，如可运移氡产生率的确定过程还需要测定其他中间参数：镭含量、射气系数以及介质的密度；同时，每一个物理参数均需要独自的实验装置。因此，若逐个测定这四个物理参数需要耗费大量的人力、物力和时间。
因此，若提供一种同时确定砂壤、废石堆、矿石堆和铀尾矿堆等颗粒堆积型射气介质的氡扩散系数、可运移氡产生率、渗透率以及孔隙率的简易方法和测量装置，不仅能减少大量的人力、物力和时间，而且也对准确开展氡在该类介质中的迁移和析出实验研究具有重要的实用价值。
发明内容
本发明的目的在于公开一套新的***：颗粒堆积型射气介质氡迁移的物理参数的确定方法及测量装置。利用该颗粒堆积型射气介质氡迁移模拟物理参数的确定方法及测量装置，由公式(7)、(9)和(15)可得出颗粒堆积型射气介质内可运移氡产生率α、孔隙率ε、介质渗透率k以及氡扩散系数D等参数。
本发明的技术方案是：
一种颗粒堆积型射气介质氡迁移物理参数的测量装置，该测量装置包括氡扩散装置、外部集氡筒体、闭式循环氡浓度测量装置以及压差测量装置，
所述氡扩散装置包括底板(1)、有纱网的孔板(2)、橡胶垫圈(3)、多孔取样管(6)以及装样筒体(7)；底板(1)、有纱网的孔板(2)、橡胶垫圈(3)和装样筒体(7)依次连接；多孔取样管(6)用纱布包裹住后***到装样筒体(7)底部，多孔取样管(6)中的圆孔朝下，测压管(12)管口用橡胶塞(5)密封。
所述外部集氡筒体包括带外接管嘴(14)的上部集氡筒体(13)和下部集氡筒体(8)；所述闭式循环氡浓度测量装置包括测氡仪(17)、干燥装置(18)、流量计(19)、循环气泵(21)、流量调节阀(20)、连接管(15)以及三岔连接管(16)；
所述压差测量装置包括测压管(12)、U型倾斜式压差计(10)、支架(9)以及连接管(11)；所述氡扩散装置下部的底板(1)去掉，氡扩散装置的下部与下部集氡筒体(8)连接，上部与上部集氡筒体(13)连接；装样筒体(7)通过连接管(11)和测压管(12)与设在支架(9)上的U型倾斜式压差计(10)连接；所述闭式循环氡浓度测量装置的上端通过外接管嘴(14)与上部集氡筒体(13)连接，下端通过外接管嘴(14)与下部集氡筒体(8)连接；
所述法兰间夹有橡胶垫圈(3)密封，以免气体泄露；
所述上部集氡筒体(13)可更换不同体积的筒体。
一种颗粒堆积型射气介质氡迁移模拟物理参数的确定方法，利用该颗粒堆积型射气介质氡迁移物理参数的确定方法及测量装置，由公式(7)、(9)和(15)得出介质内氡扩散系数D和可运移氡产生率α、介质渗透率k以及孔隙率ε参数，具体测量和确定方法如下：
(1)纯扩散平衡状态下介质内孔隙氡浓度和表面环境空气氡浓度的测量
1)氡扩散装置的准备
将底板(1)、有纱网的孔板(2)、橡胶垫圈(3)和装样筒体(7)连接；用纱布包裹住多孔取样小管(6)后***到装样筒体(7)底部，取样小管(6)中的圆孔朝下，测压管(12)管口用橡胶塞(5)密封，安装好氡扩散装置，注满水检查装置的气密性；
2)实验样品准备
确定测定对象，并采集实验样品，之后将其在温度为105℃的烘干炉中烘干；
3)装样
按照要求将烘干的实验样品分层均匀装入装样筒体(7)中，直至样品表面与上部法兰面持平；
4)静置平衡和孔隙氡浓度的测定
让装有颗粒堆积型射气介质的氡扩散装置敞开在通风良好的大气环境中，使其内产生的氡气处于纯扩散状态。静止3天后，每6小时采用针孔注射器抽取取样管(6)内的气体样，注入已准备好的闪烁室内。静置3～3.5小时后测定样品的氡浓度。孔隙氡浓度按下式计算；
式(1)中：C(x)—测定样品底部孔隙氡浓度，单位Bq·m-3；NC—本底加样品计数率的平均值，单位cpm；NB—本底计数率的平均值，单位cpm；NL—室外空气计数率，单位cpm；VS—闪烁室体积，单位m3；V—针孔注射器的取样体积，单位m3；L—刻度系数，单位Bq·m-3/cpm；
若连续3次测得的孔隙氡浓度相近，实验样品达到了扩散平衡状态；
5)环境空气氡浓度和底部孔隙扩散平衡氡浓度的确定
在此期间，采用连续测氡仪监测室内环境空气的氡浓度。由于室内空气氡浓度变化不大，取监测氡浓度的平均值为实验样品表面环境空气的氡浓度Ca。取步骤4)中确定扩散平衡状态的3次孔隙氡浓度的平均值为装样筒体底部x＝0m的稳定孔隙氡浓度Cd；
(2)两种封闭体积下装置内循环气体的稳定氡浓度和介质的气体压力梯度的测量
1)装置内循环气体的稳定氡浓度和介质的气体压力梯度测量装置的准备
如附图2所示，先将氡扩散装置下部的底板(1)去掉；接着将其下部与下部集氡筒体(8)连接，上部与上部集氡筒体(13)连接；之后将装样筒体(7)、支架(9)、U型倾斜式压差计(10)、连接管(11)和测压管(12)连接；最后将闭式循环氡浓度测量装置分别与上部集氡筒体(13)和下部集氡筒体(8)连接；
测量装置安装好后，检查整个装置的气密性。具体方法为：通过装样筒体(7)底部的取样管(6)管口向装置内鼓入一定压力的气体后，关闭鼓气管道的阀门，观察倾斜式压差计示数是否稳定，若稳定，则表明装置的气密性良好；反之，则检查装置的连接情况再做气密性测试；
气密性检查完毕后，泄压至大气压力状态后，用橡胶塞(5)密封取样管(6)管口，准备后续实验；
2)不同时刻封闭装置内循环气体氡浓度的测定
先将上部集氡筒体(13)的外接管嘴与闭式循环氡浓度测量装置断开，启动连续测氡仪(17)和循环气泵(21)，调节流量调节阀(20)使循环流量远大于连续测氡仪(17)的循环流量，通过上部集氡筒体的管嘴吸入外界气体对装置内的气体进行稀释，出口用活性炭吸附；
稀释一段时间后，将上部集氡筒体的外接管嘴(14)立即与闭式循环氡浓度测量装置连接，并设置好连续测氡仪(17)的连续监测工作模式后，将第一次监测得到的循环气体氡浓度记为C0。此后不同时刻监测得到的循环气体氡浓度依次记为C1、C2、……、Cn(n≥5)；
将上部集氡筒体(13)更换为大体积的上部集氡筒体。采用同样的方法测定另一封闭体积下封闭装置内循环气体氡浓度C'0、C′1、……、C'n(n≥5)；
3)试样内气体压力梯度的测定
在测定装置内循环气体的稳定氡浓度C2时，记录流量计(19)的读数Q、室内的温度T和U型倾斜式压差计(10)的长度差△l，两测压管(12)之间的高度差△H。试样内气体的压力梯度按下式计算：
式(2)中：△P为试样内的气体压力梯度，单位Pa/m；△H为两测压管(12)的高度差，单位m；△l为U型倾斜式压差计(10)管内液面的长度差，单位m；η为U型倾斜式压差计(10)的倾斜系数；ρ为U型倾斜式压差计(10)管内液体的密度，单位kg/m3；g为重力加速度，g＝9.8m/s2；
试样内气体的渗流速度按下式计算：
式(3)中：v为试样内气体的渗流速度，单位m/s；Q为循环气体的流量，单位m3/s；S为试样的横截面积，单位m2；
(3)颗粒堆积型射气介质内可运移氡产生率α、孔隙率ε、介质渗透率k以及氡扩散系数D的确定
1)试样的可运移氡产生率α和孔隙率ε的确定
当实验样品产生的氡气在封闭空间内循环流动。依据核素衰变理论和质量守恒定律可得平衡状态下氡浓度的微分方程：
式(4)中：V0为样品的体积，单位m3；V外为外部循环空间的体积，单位m3；λ为衰变系数，2.1×10-6/s；ε为多孔介质的孔隙率；α为样品的可运移氡产生率，单位Bq/(m3·s)；
若t＝0,C＝C0，则式(4)的解析解为：
令封闭空间循环气体的氡浓度随时间增长的解析式(5)可简化为通式C(t)＝A+(C0-A)e-λt；
当装置外接普通的上部集氡筒体时，通式中A的计为A1，当外接大体积的上部集氡筒体时，通式中A的计为A2。则：
式(6)中：V1为普通的上部集氡筒体的体积，单位m3；V2为大体积的上部集氡筒体的体积，单位m3；V3为下部集氡筒体(8)的体积，单位m3；V4为闭式循环氡浓度测量装置的体积，单位m3；
运用最小二乘法，将外接上部集氡筒体(13)时测得的不同时刻氡浓度值C0、C1、C2、……、Cn(n≥5)按C(t)＝A+(C0-A)e-λt进行拟合，得到A1的值；同理将外接上部集氡筒体(13)时测得的不同时刻氡浓度值C'0、C′1、……、C'n(n≥5)按C(t)＝A+(C0-A)e-λt进行拟合，得到A2的值；
将最小二乘法拟合得到A1和A2代入公式(7)，得到实验样品的可运移氡产生率α和孔隙率ε；
2)试样的渗透率k的确定
由达西定律可以得出渗流速度与压力梯度变化关系为：
式(8)中：k为渗透率，单位m2；μ为空气粘滞系数，单位Pa·S；
由式(2)、(3)、(8)可以得出：
3)试样的氡扩散系数D的确定
在纯扩散稳定情况下，厚度为l0颗粒堆积型射气介质226Ra衰变产生的氡在自身衰变和扩散作用下而最终趋于稳定。根据核素衰变理论、气体扩散理论和质量守恒定律建立氡气的扩散迁移微分方程为：
在稳定状态，式(10)变为：
在通风状态良好的环境中，式(11)及其边界条件为：
式(12)的解析解为：
从式(13)中可以看出，纯扩散条件下一维颗粒堆积射气介质的氡浓度分布与扩散系数、介质本身的物理特征(可运移氡产生率α和孔隙率ε)以及氡的衰变系数λ有关。由于可运移氡产生率α和孔隙率ε公式(5)已确定，λ为已知常量。故可以将x＝0,C(0)＝Cd代入解析式得：
式(14)中：D为扩散系数，单位m2/s；λ为衰变系数，2.1x10-6/s；ε为多孔介质的孔隙率；α为单位体积孔隙介质产生的运移氡量，单位Bq/(m3·s)；由式(14)可得到扩散系数D的解为：
本发明提供的颗粒堆积型射气介质内氡的可运移氡产生率α、孔隙率ε、介质渗透率k以及氡扩散系数D的计算方法及其测量装置的优点如下：
1、通过改变外部集氡体的体积，并测量稳定时期筒内氡的浓度计算出介质内氡的扩散系数D、可运移氡产生率α、渗透率k以及孔隙率ε，大大提高了计算结果的准确性，可以避免查找资料的繁琐和中间测定参数的步骤，这样既可以保证参数的准确度，又可以节省人力，物力，财力；
2、此计算方法和装置也可以适用于其它材料参数的测定；
3、测量装置具有制作简单，成本低，适用性强，可重复使用。
附图说明
图1为氡扩散装置示意图；
图2为装置内循环气体的氡浓度和介质内气体压力梯度测量装置一；
图3为装置内循环气体的氡浓度和介质内气体压力梯度测量装置二。
具体实施方式
一种颗粒堆积型射气介质氡迁移的物理参数的确定方法，这里以两种样品为例，说明颗粒堆积型射气介质氡迁移的物理参数的确定。
样品1：来自新上堆浸矿堆的矿石，矿石中的矿物由热液矿物和围岩残留矿物组成，矿山提供的矿物分析资料表明主要铀矿物为沥青铀矿及少量铀的次生矿物，并与黄铁矿、石英等共生。本实验用样品密度为2.61×103kg/m3，粒径0mm～9mm，铀品位约为0.145％，铀镭平衡系数约为1.04；
样品2：实验用砂壤土来自南华大学卧虎山。样品的镭含量为15Bq/kg，粒径为0mm～2mm，松散堆积密度约为1.6g/cm3；
1、实验步骤
(1)纯扩散平衡状态下介质内孔隙氡浓度和表面环境空气氡浓度的测量
测量装置尺寸以及其它参数如下：
装样筒体(7)高度为0.5m，内径为0.1m，厚度为6mm，则样品横截面积A为A＝7.8500×10-3m2，装样筒体体积(含橡胶垫圈的体积)为V＝3.925×10-3m3；样品1和2的体积为V0＝V＝3.925×10-3m3；
下部集氡筒体(8)高度为0.125m，内径为0.1m，厚度为6mm，体积为V3＝0.98125×10-3m3；外部测量管均采用的是外径为6mm，内径为4mm的橡皮管道，总长度为1.5m，总体积为V4＝0.72L(包含测氡仪内部体积0.7L和管道的体积0.0188L)；
普通上部集氡筒体体积(含橡胶垫圈的体积)为V1＝0.98125×10-3m3。大体积的上部集氡筒体积(含橡胶垫圈的体积)为V2＝1.9625×10-3m3，室内温度T约为T＝20℃，则空气粘滞系数μ为μ＝1.8107×10-5Pa·S，U型倾斜式压差计(10)管内液体的密度为ρ＝1×103kg/m3，流量计(19)为微型流量计，其测量范围为0.5-8L/min，U型倾斜式压差计(10)的测量范围为0～150mm，其倾斜系数η为：0.1、0.2、0.5、0.7，循环气泵(21)为微型气泵，最大压力为40kpa，最大循环流量为11L/min；
如附图1所示安装好氡扩散装置后，注满水检查装置的气密性；
将实验样品在维持温度为105℃的烘干炉中烘6小时。将烘干的实验样品分层均匀装入装样筒体(7)中，直至样品表面与上部法兰面持平。让装有颗粒堆积型射气介质的氡扩散装置敞开在通风良好的大气环境中，使其内产生的氡气处于纯扩散状态。静止3天后，每6小时采用20ml的针孔注射器抽取取样管(6)内的气体样，注入已准备好的闪烁室内。静置3～3.5小时后测定样品的氡浓度，孔隙氡浓度按式(1)计算。若连续3次测得的孔隙氡浓度相近，实验样品达到了扩散平衡状态；
在此期间，采用连续测氡仪监测室内环境空气的氡浓度。取监测氡浓度的平均值为实验样品表面环境空气的氡浓度。取扩散平衡状态的3次孔隙氡浓度的平均值为装样筒体底部x＝0m的稳定孔隙氡浓度。样品1和样品2底部(x＝0)孔隙氡浓度与表面环境空气氡浓度分别见表1及表2；
表1样品1内底部(x＝0)孔隙氡浓度Cd与表面环境空气氡浓度Ca
表2样品2底部(x＝0)孔隙氡浓度Cd与表面环境空气氡浓度Ca
(2)两种封闭体积下装置内循环气体的稳定氡浓度和介质的气体压力梯度的测量
先将上部集氡筒体(13)的管嘴与闭式循环氡浓度测量装置断开，启动连续测氡仪(17)和循环气泵(21)，调节流量调节阀(20)使循环流量远大于连续测氡仪(17)的循环流量，通过上部集氡筒体的管嘴吸入外界气体对装置内的气体进行稀释，出口用活性炭吸附；
稀释30min之后将装置连接封闭，测量循环气体氡浓度记为C0。此后按监测周期为30min监测得到的循环气体氡浓度依次记为C1、C2、……、Cn(n≥5)。将普通上部集氡筒体更换为大体积的上部集氡筒体，采用同样的方法测定另一封闭体积下装置内循环气体氡浓度C'0、C′1、……、C'n；
测定样品2的监测周期为60min，其他同样品1；
当测试对象为样品1和样品2时不同封闭体积装置内循环气体的氡浓度监测结果分别见表3及表4；
表3测试对象为样品1时不同封闭体积装置内循环气体的氡浓度监测结果
表4测试对象为样品2时不同封闭体积装置内循环气体的氡浓度监测结果
| 监测序号 | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
| C1(t)/(Bq·m3) | 57 | 146 | 211 | 301 | 426 | 475 | 558 | 691 | 767 | 826 |
| C2(t)/(Bq·m3) | 36 | 111 | 157 | 238 | 326 | 380 | 471 | 549 | 611 | 648 |
表3中相邻监测点时间间隔为30min，即相隔1800s；表4中相邻监测点时间间隔为60min，即相隔1800s。将表3和表4中的C1(t)和C2(t)分别按C(t)＝A+(C0-A)e-λt进行拟合，拟合结果见表5；
表5拟合结果
3)试样内气体压力梯度的测定
在测定装置内循环气体的稳定氡浓度C2时，记录流量计(19)的读数Q、室内的温度T和U型倾斜式压差计(10)的长度差△l。测试样品1时倾斜系数选取为0.2，测试样品2时倾斜系数选取为0.7；两测压管(12)之间的高度差△H。试样内气体的压力梯度按式子(2)计算；试样内气体的渗流速度按式子(3)计算。试样内气体渗流相关参数监测结果见表5；
表6试样内气体渗流相关参数监测结果
(3)据此可以得出样品氡迁移的物理参数：
将表1-表5中的监测数据代入本发明提供的公式(7)、(9)和(15)得到多孔介质的孔隙率ε、单位体积孔隙介质产生的运移氡量α、多孔介质的渗透率k以及扩散系数D；
表7两个试样氡迁移的物理参数
Claims (5)
1.颗粒堆积型射气介质氡迁移物理参数的测量装置，包括氡扩散装置、外部集氡筒体、闭式循环氡浓度测量装置以及压差测量装置，其特征在于，
所述氡扩散装置包括底板(1)、有纱网的孔板(2)、橡胶垫圈(3)、多孔取样管(6)以及装样筒体(7)；底板(1)、有纱网的孔板(2)、橡胶垫圈(3)和装样筒体(7)依次连接；多孔取样管(6)用纱布包裹住后***到装样筒体(7)底部，多孔取样管(6)中的圆孔朝下，测压管(12)管口用橡胶塞(5)密封；
所述外部集氡筒体包括带外接管嘴(14)的上部集氡筒体(13)和下部集氡筒体(8)；所述闭式循环氡浓度测量装置包括测氡仪(17)、干燥装置(18)、流量计(19)、循环气泵(21)、流量调节阀(20)、连接管(15)以及三岔连接管(16)；测氡仪(17)、干燥装置(18)连接在三岔连接管(16)上，三岔连接管(16)下端依次连接流量计(19)和循环气泵(21)，循环气泵(21)通过流量调节阀(20)连接到下部集氡筒体(8)，三岔连接管(16)上端连接连接管(15)，连接管(15)连接到上部集氡筒体(13)的外接管嘴(14)；
所述压差测量装置包括测压管(12)、U型倾斜式压差计(10)、支架(9)以及连接管(11)；
所述氡扩散装置下部的底板(1)去掉，氡扩散装置的下部与下部集氡筒体(8)连接，上部与上部集氡筒体(13)连接；装样筒体(7)通过连接管(11)和测压管(12)与设在支架(9)上的U型倾斜式压差计(10)连接；所述闭式循环氡浓度测量装置的上端通过外接管嘴(14)与上部集氡筒体(13)连接，下端通过外接管嘴(14)与下部集氡筒体(8)连接。
2.根据权利要求1所述的颗粒堆积型射气介质氡迁移物理参数的测量装置，其特征在于，所述法兰间夹有橡胶垫圈(3)密封，以免气体泄露。
3.根据权利要求1所述的颗粒堆积型射气介质氡迁移物理参数的测量装置，其特征在于，所述上部集氡筒体(13)可更换不同体积的筒体。
4.一种颗粒堆积型射气介质氡迁移模拟物理参数的确定方法，其特征在于，具体测量和确定方法如下：
(1)纯扩散平衡状态下介质内孔隙氡浓度和表面环境空气氡浓度的测量，具体包括：
1)氡扩散装置的准备
将底板(1)、有纱网的孔板(2)、橡胶垫圈(3)和装样筒体(7)连接；用纱布包裹住多孔取样小管(6)后***到装样筒体(7)底部，多孔取样管(6)中的圆孔朝下，测压管(12)管口用橡胶塞(5)密封，安装好氡扩散装置，注满水检查装置的气密性；
2)实验样品准备
确定测定对象，并采集实验样品，将其在温度为105℃的烘干炉中烘干，
3)装样
将烘干的实验样品分层均匀装入装样筒体(7)中，直至样品表面与上部法兰面持平；
4)静置平衡和孔隙氡浓度的测定
让装有颗粒堆积型射气介质的氡扩散装置敞开在通风良好的大气环境中，使其内产生的氡气处于纯扩散状态，静止3天后，每6小时采用针孔注射器抽取多孔取样管(6)内的气体样，注入已准备好的闪烁室内，静置3～3.5小时后测定样品的氡浓度，孔隙氡浓度按下式计算
<mrow>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mi>S</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>C</mi>
</msub>
<mo>-</mo>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>S</mi>
</msub>
<mo>-</mo>
<mi>V</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>N</mi>
<mi>L</mi>
</msub>
</mrow>
<mi>V</mi>
</mfrac>
<mo>&times;</mo>
<mi>L</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式(1)中：C(x)—测定样品底部孔隙氡浓度，单位Bq·m-3；NC—本底加样品计数率的平均值，单位cpm；NB—本底计数率的平均值，单位cpm；NL—室外空气计数率，单位cpm；VS—闪烁室体积，单位m3；V—针孔注射器的取样体积，单位m3；L—刻度系数，单位Bq·m-3/cpm；
若连续3次测得的孔隙氡浓度相近，实验样品达到了扩散平衡状态；
5)环境空气氡浓度和底部孔隙扩散平衡氡浓度的确定
采用连续测氡仪监测室内环境空气的氡浓度，由于室内空气氡浓度变化不大，取监测氡浓度的平均值为实验样品表面环境空气的氡浓度Ca，取步骤4)中确定扩散平衡状态的3次孔隙氡浓度的平均值为装样筒体底部x＝0m的稳定孔隙氡浓度Cd。；
(2)两种封闭体积下装置内循环气体的稳定氡浓度和介质的气体压力梯度的测量，具体包括：
2.1)装置内循环气体的稳定氡浓度和介质的气体压力梯度测量装置的准备
先将氡扩散装置下部的底板(1)去掉；接着将其底部与下部集氡筒体(8)连接，上部与上部集氡筒体(13)连接；之后将装样筒体(7)、支架(9)、U型倾斜式压差计(10)、连接管(11)和测压管(12)连接；最后将闭式循环氡浓度测量装置分别与上部集氡筒体(13)和下部集氡筒体(8)连接；
测量装置安装好后，检查整个装置的气密性；
气密性检查完毕后，泄压至大气压力状态后，用橡胶塞(5)密封取样管(6)管口，准备后续实验；
2.2)不同时刻封闭装置内循环气体氡浓度的测定
先将上部集氡筒体(13)的外接管嘴与闭式循环氡浓度测量装置断开，启动连续测氡仪(17)和循环气泵(21)，调节流量调节阀(20)使循环流量远大于连续测氡仪(17)的循环流量，通过上部集氡筒体(13)的外接管嘴吸入外界气体对装置内的气体进行稀释，出口用活性炭吸附；
稀释后，将上部集氡筒体的外接管嘴(14)立即与闭式循环氡浓度测量装置连接，并设置好连续测氡仪(17)的连续监测工作模式后，将第一次监测得到的循环气体氡浓度记为C0，此后不同时刻监测得到的循环气体氡浓度依次记为C1、C2、……、Cn(n≥5)；
将普通的上部集氡筒体更换为大体积的上部集氡筒体，采用同样的方法测定该封闭体积下封闭装置内循环气体氡浓度C'0、C1'、……、C'n(n≥5)；
2.3)试样内气体压力梯度的测定
在测定装置内循环气体的稳定氡浓度C2时，记录流量计(19)的读数Q、室内的温度T和U型倾斜式压差计(10)的长度差△l，两测压管(12)之间的高度差△H，试样内气体的压力梯度按下式计算：
<mrow>
<mi>&Delta;</mi>
<mi>P</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&rho;</mi>
<mi>g</mi>
<mi>&Delta;</mi>
<mi>l</mi>
<mi>&eta;</mi>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>H</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式(2)中：△P为试样内的气体压力梯度，单位Pa/m；△H为两测压管(12)的高度差，单位m；△l为U型倾斜式压差计(10)管内液面的长度差，单位m；η为U型倾斜式压差计(10)的倾斜系数；ρ为U型倾斜式压差计(10)管内液体的密度，单位kg/m3；g为重力加速度，g＝9.8m/s2；
试样内气体的渗流速度按下式计算：
<mrow>
<mi>v</mi>
<mo>=</mo>
<mfrac>
<mi>Q</mi>
<mi>S</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式(3)中：v为试样内气体的渗流速度，单位m/s；Q为循环气体的流量，单位m3/s；S为试样的横截面积，单位m2；
(3)颗粒堆积型射气介质内可运移氡产生率α、孔隙率ε、介质渗透率k以及氡扩散系数D的确定，具体包括：
3.1)试样的可运移氡产生率α和孔隙率ε的确定
当实验样品产生的氡气在封闭空间内循环流动，依据核素衰变理论和质量守恒定律得氡浓度的微分方程：
式(4)中：V0为样品的体积，单位m3；V外为外部循环空间的体积，单位m3；λ为衰变系数，2.1×10-6/s；ε为多孔介质的孔隙率；α为样品的可运移氡产生率，单位Bq/(m3·s)；
若t＝0,C＝C0，则式(4)的解析解为：
令封闭空间循环气体的氡浓度随时间增长的解析式(5)简化为通式C(t)＝A+(C0-A)e-λt，
当装置外接普通的上部集氡筒体(13)时，通式中A的记为A1，外接大体积的上部集氡筒体时，通式中A的记为A2，则：
式(6)中：V1为普通的上部集氡筒体的体积，单位m3；V2为大体积的上部集氡筒体的体积，单位m3；V3为下部集氡筒体(8)的体积，单位m3；V4为闭式循环氡浓度测量装置的体积，单位m3。
运用最小二乘法，将外接普通的上部集氡筒体时测得的不同时刻循环气体的氡浓度值C0、C1、C2、……、Cn(n≥5)按C(t)＝A+(C0-A)e-λt进行拟合，得到A1的值；同理将外接大体积的上部集氡筒体时测得的不同时刻循环气体的氡浓度值C'0、C1'、……、C'n(n≥5)按C(t)＝A+(C0-A)e-λt进行拟合，得到A2的值；
将最小二乘法拟合得到A1和A2代入公式(6)，得到实验样品的孔隙率ε和可运移氡产生率α，
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&epsiv;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mn>4</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mn>4</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&alpha;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mn>4</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mn>4</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
3.2)试样的渗透率k的确定
由达西定律得出渗流速度与压力梯度变化关系为：
<mrow>
<mi>v</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&kappa;</mi>
<mi>&Delta;</mi>
<mi>P</mi>
</mrow>
<mi>&mu;</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式(8)中：k为渗透率，单位m2；μ为空气粘滞系数，单位Pa·S，
由式(2)、(3)和(8)得出：
<mrow>
<mi>&kappa;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>Q</mi>
<mi>&mu;</mi>
<mi>&Delta;</mi>
<mi>H</mi>
</mrow>
<mrow>
<mi>S</mi>
<mi>&rho;</mi>
<mi>g</mi>
<mi>&Delta;</mi>
<mi>l</mi>
<mi>&eta;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
3.3)试样的氡扩散系数D的确定
在纯扩散稳定情况下，厚度为l0颗粒堆积型射气介质226Ra衰变产生的氡在自身衰变和扩散作用下而最终趋于稳定，根据核素衰变理论、气体扩散理论和质量守恒定律建立氡气的扩散迁移微分方程为：
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>C</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>D</mi>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>C</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mi>&lambda;</mi>
<mi>C</mi>
<mo>+</mo>
<mfrac>
<mi>&alpha;</mi>
<mi>&epsiv;</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
在稳定状态，式(10)变为：
<mrow>
<mn>0</mn>
<mo>=</mo>
<mi>D</mi>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>C</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mi>&lambda;</mi>
<mi>C</mi>
<mo>+</mo>
<mfrac>
<mi>&alpha;</mi>
<mi>&epsiv;</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
在通风状态良好的环境中，式(11)及其边界条件为：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>D</mi>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>C</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mi>&lambda;</mi>
<mi>C</mi>
<mo>+</mo>
<mfrac>
<mi>&alpha;</mi>
<mi>&epsiv;</mi>
</mfrac>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>C</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mo>=</mo>
<msub>
<mi>l</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<mi>C</mi>
<mo>=</mo>
<msub>
<mi>C</mi>
<mi>a</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
式(12)的解析解为：
<mrow>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>&alpha;</mi>
<mrow>
<mi>&lambda;</mi>
<mi>&epsiv;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mrow>
<mi>x</mi>
<msqrt>
<mfrac>
<mi>&lambda;</mi>
<mi>D</mi>
</mfrac>
</msqrt>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>x</mi>
<msqrt>
<mfrac>
<mi>&lambda;</mi>
<mi>D</mi>
</mfrac>
</msqrt>
</mrow>
</msup>
<mo>)</mo>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>a</mi>
</msub>
<mi>&lambda;</mi>
<mi>&epsiv;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mi>&epsiv;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>l</mi>
<mn>0</mn>
</msub>
<msqrt>
<mfrac>
<mi>&lambda;</mi>
<mi>D</mi>
</mfrac>
</msqrt>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>0</mn>
</msub>
<msqrt>
<mfrac>
<mi>&lambda;</mi>
<mi>D</mi>
</mfrac>
</msqrt>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
从式(13)中看出，纯扩散条件下一维颗粒堆积射气介质的氡浓度分布与扩散系数、可运移氡产生率α、孔隙率ε，及氡的衰变系数λ有关，由于可运移氡产生率α和孔隙率ε公式(5)已确定，λ为已知常量，故将x＝0,C(0)＝Cd代入解析式得：
<mrow>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>&alpha;</mi>
<mrow>
<mi>&lambda;</mi>
<mi>&epsiv;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>a</mi>
</msub>
<mi>&lambda;</mi>
<mi>&epsiv;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mi>&epsiv;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>l</mi>
<mn>0</mn>
</msub>
<msqrt>
<mfrac>
<mi>&lambda;</mi>
<mi>D</mi>
</mfrac>
</msqrt>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>0</mn>
</msub>
<msqrt>
<mfrac>
<mi>&lambda;</mi>
<mi>D</mi>
</mfrac>
</msqrt>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
式(14)中：D为扩散系数，单位m2/s；λ为衰变系数，2.1x10-6/s；ε为多孔介质的孔隙率；α为单位体积孔隙介质产生的运移氡量，单位Bq/(m3·s)；由式(14)得到扩散系数D的解为：
<mrow>
<mi>D</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>l</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mi>&lambda;</mi>
</mrow>
<mrow>
<mi>I</mi>
<mi>n</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>a</mi>
</msub>
<mi>&lambda;</mi>
<mi>&epsiv;</mi>
<mo>+</mo>
<msqrt>
<mrow>
<mo>-</mo>
<msubsup>
<mi>C</mi>
<mi>d</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>C</mi>
<mi>a</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&alpha;C</mi>
<mi>d</mi>
</msub>
<mi>&lambda;</mi>
<mi>&epsiv;</mi>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&alpha;C</mi>
<mi>a</mi>
</msub>
<mi>&lambda;</mi>
<mi>&epsiv;</mi>
</mrow>
</msqrt>
<mo>-</mo>
<mi>&alpha;</mi>
</mrow>
<mrow>
<mi>&alpha;</mi>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<mi>&lambda;</mi>
<mi>&epsiv;</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
5.根据权利要求4所述的一种颗粒堆积型射气介质氡迁移模拟物理参数的确定方法，其特征在于，检查整个装置的气密性，具体方法为：通过装样筒体(7)底部的取样管(6)管口向装置内鼓入气体后，关闭鼓气管道的阀门，观察倾斜式压差计示数是否稳定，若稳定，则表明装置的气密性良好；反之，则检查装置的连接情况再做气密性测试。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610565316.4A CN106197544B (zh) | 2016-07-18 | 2016-07-18 | 颗粒堆积型射气介质氡迁移物理参数的确定方法及其测量装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610565316.4A CN106197544B (zh) | 2016-07-18 | 2016-07-18 | 颗粒堆积型射气介质氡迁移物理参数的确定方法及其测量装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106197544A CN106197544A (zh) | 2016-12-07 |
| CN106197544B true CN106197544B (zh) | 2018-05-18 |
Family
ID=57494377
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610565316.4A Active CN106197544B (zh) | 2016-07-18 | 2016-07-18 | 颗粒堆积型射气介质氡迁移物理参数的确定方法及其测量装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106197544B (zh) |
Families Citing this family (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106644823B (zh) * | 2016-12-21 | 2023-03-28 | 山西大学 | 煤壁涂料的气密性测试装置及其气密性测试方法 |
| CN106886041B (zh) * | 2017-03-07 | 2019-01-04 | 中南大学 | 压差可调式射气岩石介质氡析出率测量装置及其使用方法 |
| CN106959464B (zh) * | 2017-03-29 | 2023-06-09 | 南华大学 | 一种氡析出率的测量装置和测量方法 |
| CN106842275B (zh) * | 2017-04-10 | 2023-08-01 | 福建师范大学 | 自然环境条件下表面氡析出率测量装置及方法 |
| CN106873023B (zh) * | 2017-04-20 | 2019-02-01 | 福建师范大学 | 一种利用活性炭盒测量空气氡浓度的方法 |
| CN107704682B (zh) * | 2017-09-30 | 2021-08-10 | 西南科技大学 | 基于概率用于核素近场、远场迁移评估的空间域描述方法 |
| CN108931465B (zh) * | 2018-05-21 | 2020-07-14 | 南华大学 | 测定多孔射气介质中氡扩散系数和可运移氡产生率的方法 |
| CN108897031B (zh) * | 2018-05-21 | 2020-02-21 | 南华大学 | 一种开环式测量射气介质氡析出率的方法和装置 |
| CN108680467B (zh) * | 2018-05-21 | 2020-07-10 | 南华大学 | 地下工程原位测定氡扩散系数和可运移氡产生率的方法 |
| CN108896467B (zh) * | 2018-08-28 | 2020-12-29 | 衡阳师范学院 | 一种动态测量土壤孔隙率的方法及装置 |
| CN109030785B (zh) * | 2018-08-28 | 2021-05-11 | 衡阳师范学院 | 一种利用表层土壤动态测量土壤潜势氡浓度的方法及装置 |
| CN108918395B (zh) * | 2018-09-07 | 2021-03-16 | 衡阳师范学院 | 一种利用静电收集测氡仪测量土壤孔隙率的方法及装置 |
| CN108827859B (zh) * | 2018-09-11 | 2024-04-02 | 东华理工大学 | 一种氡气扩散式岩石有效孔隙度测量装置及方法 |
| CN109060633B (zh) * | 2018-09-14 | 2019-06-28 | 衡阳师范学院 | 一种静态测量土壤孔隙率的方法及装置 |
| CN108872544B (zh) * | 2018-09-14 | 2020-12-29 | 衡阳师范学院 | 一种利用表层土壤静态测量土壤潜势氡浓度的方法及装置 |
| CN109283115B (zh) * | 2018-11-09 | 2021-08-31 | 福州大学 | 一种基于溶质挥发体系的微胶囊内部空隙体积的计算方法 |
| CN109323962A (zh) * | 2018-11-16 | 2019-02-12 | 核工业北京化工冶金研究院 | 一种破碎射气介质多功能氡扩散渗流试验装置 |
| CN109269943B (zh) * | 2018-11-27 | 2021-04-20 | 衡阳师范学院 | 快速测量氡在薄膜中有效扩散系数的方法 |
| CN109655856B (zh) * | 2018-12-14 | 2022-02-18 | 衡阳师范学院 | 静电收集法快速测量氡在薄膜中有效扩散系数的装置及方法 |
| CN110927017A (zh) * | 2019-11-11 | 2020-03-27 | 南华大学 | 一种确定多孔射气介质中自由220Rn迁移参数的装置及方法 |
| CN114264585B (zh) * | 2020-09-16 | 2023-11-14 | 宝山钢铁股份有限公司 | 生产用碳铁复合炉料的透气性模拟测定方法 |
| CN116429663B (zh) * | 2023-06-08 | 2023-09-12 | 太原理工大学 | 一种测量煤岩介质中氡气渗流速率的装置及方法 |
| CN117110172B (zh) * | 2023-10-24 | 2024-01-19 | 中国矿业大学 | 一种气体渗流启动压力梯度测试装置及应用方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101620274A (zh) * | 2009-08-12 | 2010-01-06 | 成都理工大学 | 基于壤中氡静态扩散静电吸附原理壤中氡α能谱测量方法 |
| CN103439398A (zh) * | 2013-09-06 | 2013-12-11 | 北华航天工业学院 | 一种基于模糊逻辑的氡检测***和方法 |
| CN104181283A (zh) * | 2014-09-05 | 2014-12-03 | 河南理工大学 | 受载煤破裂过程中氡气析出测定方法及装置 |
| CN104863631A (zh) * | 2015-04-22 | 2015-08-26 | 南华大学 | 降低铀矿气态流出物对周边环境危害的优化方法 |
| CN105067789A (zh) * | 2015-10-01 | 2015-11-18 | 南华大学 | 一种开环式原位测量水体氡析出率的方法和装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140061117A1 (en) * | 2012-04-23 | 2014-03-06 | Aqua Ip | Water filter with multiple internal filtration modules |
-
2016
- 2016-07-18 CN CN201610565316.4A patent/CN106197544B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101620274A (zh) * | 2009-08-12 | 2010-01-06 | 成都理工大学 | 基于壤中氡静态扩散静电吸附原理壤中氡α能谱测量方法 |
| CN103439398A (zh) * | 2013-09-06 | 2013-12-11 | 北华航天工业学院 | 一种基于模糊逻辑的氡检测***和方法 |
| CN104181283A (zh) * | 2014-09-05 | 2014-12-03 | 河南理工大学 | 受载煤破裂过程中氡气析出测定方法及装置 |
| CN104863631A (zh) * | 2015-04-22 | 2015-08-26 | 南华大学 | 降低铀矿气态流出物对周边环境危害的优化方法 |
| CN105067789A (zh) * | 2015-10-01 | 2015-11-18 | 南华大学 | 一种开环式原位测量水体氡析出率的方法和装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106197544A (zh) | 2016-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106197544B (zh) | 颗粒堆积型射气介质氡迁移物理参数的确定方法及其测量装置 | |
| CN108919329B (zh) | 一种闭环式测量射气介质氡析出率的方法和装置 | |
| CN105158788B (zh) | 开环式同步测量有效衰变常数及氡析出率的方法 | |
| CN106886041B (zh) | 压差可调式射气岩石介质氡析出率测量装置及其使用方法 | |
| CN106124354B (zh) | 高含硫气藏的硫溶解度在线测试装置及方法 | |
| CN109490146B (zh) | 陆地-大气界面气体交换通量的标定***及方法 | |
| CN105425267A (zh) | 一种氡析出率的测量***及方法 | |
| CN104502163B (zh) | 一种可调深度土壤二氧化碳采集装置及其使用方法 | |
| CN105067789B (zh) | 一种开环式原位测量水体氡析出率的方法和装置 | |
| CN106707325B (zh) | 用于测量一定时间段内介质表面氡释放总量的方法与*** | |
| ITPD20080204A1 (it) | Apparato metodo e apparato di rilevamento della frazione di origine remota di radon presente in un sito di misura. | |
| WO2023240817A1 (zh) | 一种连续测试变吸力下非饱和土体气体渗透系数的装置和方法 | |
| CN105929443A (zh) | 主动任意调节氡析出率及有效衰变常数的标准装置及方法 | |
| Ye et al. | Simultaneous determination of the radon diffusion coefficient and the free radon production rate from compact porous emanation media | |
| CN108732329B (zh) | 一种煤层瓦斯压力测定相似模拟实验装置及方法 | |
| CN108614001A (zh) | 测定半挥发性有机污染物释放量的微池热脱附装置及应用 | |
| CN102608001A (zh) | 氢扩散渗透特性测量装置 | |
| CN110320072A (zh) | 一种用于水中总α、总β测量的样品源制备装置及方法 | |
| CN108680467B (zh) | 地下工程原位测定氡扩散系数和可运移氡产生率的方法 | |
| CN109814144A (zh) | 一种氡和氡子体同步测量及剂量率评价方法 | |
| CN108897031B (zh) | 一种开环式测量射气介质氡析出率的方法和装置 | |
| CN115753543B (zh) | 考虑概率分布的页岩支撑裂缝相对渗透率测定装置及方法 | |
| CN108489771B (zh) | 流气式土壤氡气源收集装置及收集方法 | |
| CN105973337A (zh) | 一种气体产物计量与干燥取样的实验装置及实验方法 | |
| CN105301201A (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |