氡是天然放射性铀系中一种惰性气体，具有极强迁移活动性，但其化学性质极不活泼，一般情况下很难与其它元素发生化学反应。通常，有空气的地方，就存有氡及其子体。自然界，氡存在3种天然放射性同位素222Rn、220Rn和219Rn，其分别来源于天然放射系中的铀系、钍系和锕系，半衰期分别为3.825 d、54.5 s和3.92 s。由于220Rn和219Rn半衰期短，222Rn半衰期相对较长，故通常所指的氡为226Ra直接衰变的产物222Rn。氡极易溶于液体(水)和脂肪，通过呼吸道与皮肤侵入机体，其衰变子体常吸附在呼吸道，造成内照射，该照射剂量在体内累积到一定程度，引发机体病变，最终导致肺癌。现代研究表明室内氡浓度长期暴露会增加肺癌的危险度。此外，氡可能引发畸形、基因畸形遗传和白血病等不良后果[1]。因此，室内环境中氡浓度对人体健康的影响已成为理论界亟待解决的重大问题。要研究氡浓度对人体健康的影响，就必须先搞清室内氡浓度与其来源之间的可量性关系。
近5 a对室内氡的研究大多集中于室内氡的检测、水平扩散系数测量与有效防护方法等，如王振基等[2]对标准氡室内氡气体积活度的计算方法与测量验证的研究；连福龙等[3]对主动式活性炭盒法测氡影响因素的探讨；竺云波等[4]对居室氡浓度及其影响因素的分析；张强等[5]对建材制品中测定氡的影响因素及其在防氡建材分析中的应用研究；李先杰等[6]对氡析出率测量仪的检定装置的介绍；王榆元等[7]探讨了含水率对建材砖氡析出率的影响；李晓玲[8]对中国室内氡污染现状及相关控制标准的相关综述；其他国家如UNSCEAR对室内氡的研究发出的最新报告[9]；Ielsch等[10]提出了一种利用岩石和土壤的地球化学和物理性质参数来计算土壤表面的氡析出率的模式TRACHGE；Fazalur等人用固体核径迹法和静电收集法对不同种类建材的析出率进行了测量；EScobar等[11]用活性炭吸附法测量土壤的222Rn析出率[12]。然而，对室内氡浓度分析较零散(与一种氡源的关系分析)。综合定量分析室内氡浓度与氡来源的关系，能使室内氡浓度对人体的危害更具有可衡量性与可操作性，为人类如何降低室内氡提供科学指导。
本文在前人研究成果的基础上，拟通过对室内氡来源现状分析建立室内氡浓度的定量计算模型，从而为室内氡浓度的研究与评价提供参考依据。
现代研究表明：室内氡主要来源于建筑物地基和周围土壤、建筑材料、空气中的氡、家用燃料以及供水[13-14]。一些学者已对室内氡的氡源进入率作了相关研究。
土壤中氡可自由扩散，形成氡气体渗流进入室内，但土壤的氡及其衰变源进入室内，受诸多因素的影响，它不但取决于土壤的湿度、密度、孔隙率、房子的建筑结构以及气象参数中的气温、气压、风向和风速等客观因素，还受人为因素的影响，如换气率。任天山等人认为土壤氡进入率可以表示为[15]：
式中：C为土壤氡的进入率；Xs为土壤气体中的氡浓度，Bq/m3；Qs为渗流率，即单位时间内进入室内的土壤气体的体积占房间体积的份数，h-1。
现代研究发现，建筑材料衰变过程中产生的氡部分能通过扩散进入室内，其进入率主要受建筑材料性状、表面积、建筑材料的氡扩散长度以及建材的厚度等影响。据此，学者们认为建筑材料析氡进入率可以表示为[16]：
式中：N=3 600;v是房间体积, m3;s1和s2分别为墙壁和地板的表面积, m2；r1和r2分别为砖和混凝土的面积析氡率。
某些情况下, 室外空气中的氡是室内氡的一个重要来源，主要进入途径是室内外空气交换。假定室外空气中的氡浓度Aout(Bq/m3), 换气率为λ(/h)，则氡进入率为[15]：
人类在家中使用燃料时，若遇到通风不畅，便产生天然气，据研究显示，天然气中氡的浓度有时极高，可能是室内氡的重要来源之一，相关研究显示，天然气所致氡的进入率计算公式可以表示成[17]：
其中:Dng为天然气中的氡浓度，Bq/m3；Qng为消耗率, m3/h；v为房间体积。
设水中氡浓度为Xw, Bq/m3，水的消耗量为Qw, m3/h，脱气系数为ε，房子体积为v, 则由水氡脱气所致的进氡率E为[18]:
基于对前人研究的探讨，并结合室内氡来源的特性，可看出一般情况下家用燃料产生的天然气、供水、建筑材料以及建筑地基土壤中氡进入室内部分相对较稳定，但室外空气氡进入率最不稳定，易受房间通风性能的影响。不同来源的氡进入室内以后，室内222Rn随时间不断积累；同时氡不断衰变，加之受通风影响，从而使室内氡浓度发生变化。已有研究显示，房间的通风是影响室内氡浓度的关键因子之一[15,20]，根据自然界空气流动的特性，空气流动方向总是从浓度高的流向浓度低的形成一个小空气循环。因此，通风时室外空气流入室内的同时，室内空气流出室外，室内外空气交换随之产生，室内氡浓度随之改变。要合理计算出不同时刻室内氡的浓度，就得分析房间的通风性能；而要分析房间的通风性能必然涉及到房间送风量与房间体积。通风研究显示，房间送风量受房间的结构、位置、送风方式以及室内外空气正负压变差等诸多因子的制约。由于生活中房间结构与位置确定后为静态，一般不轻易改变。因此，暂且把房间的通风性能取决于房间送风方式和室内外空气的压差。房间送风方式受空气通过某介质进入房间接触表面积的大小与当时室内外空气的流动速率等影响，故模型引入通风量系数S 和室内外空气流动交换速率vt 2个变量。通风量系数S 指某时刻空气通过某介质如窗户、门等进入室内时，介质表面积的大小，m2；室内外空气流动速率vt指某时刻通风时，空气通过某介质如窗户、门等进入室内的空气的进入速率，m/h，大小与送风方式、室内空气变差程度以及房间通风道朝向等诸多因子有关。通风状态下，随时间的推移，房间空气的总体积总保持相对一致，即房间流出的空气体积v出等于流进室内的空气体积v进，又因为室内外空气交换时，房间的通风量系数与通风时间相对一定。因此，室内空气流动速率应该与室外空气流动速率一致，否则v进与v出不相等。因此，通风条件下，可把室内外空气流动速率vt近似看成一个常数，不随通风量系数S变化。在上述前提下，我们可近似地把通风性能问题理解为通风时，空气进入室内介质表面积大小的改变，即房间通风量系数S的大小的变化。如当完全不通风时，即S为0，不发生室内外空气交换，室内外空气交换流动速率也为0；当通风时，即S>0, 室内外空气发生交换，产生室内外空气流动速率vt，当完全通风时，即S为房间通风通道的总表面积。其理论实质为换气次数模型的再次演绎，把分子中房间送风量Ma二次变量递推化，即房间通风量系数S与空气流动速率的乘积vt与房间总体积的比，h-1。代入换气次数模型中得：${\lambda _{\text{v}}} = {M_{\text{a}}}/v = S{v_{\text{t}}}/v$, 式中，λv换气次数, h-1；Ma为房间送风量，m3/h；v为房间体积, m3；S为通风量系数，m2；vt为室内外空气流动交换速率，m/h。
综上所述，假定室内外空气氡浓度均一混合，空气中氡的浓度一致；氡均匀地散发到房间的整个空间；送入室内的空气立即与室内空气充分混合，送风量等于排风量，则室内氡的变化可用下式描述：
式中：CRn是t时刻的氡浓度, Bq/m3；λRn是氡的衰变常数, 7.6×10-3h-1；B、C、D、E分别为建筑材料所致氡进入率、建筑物地基土壤氡进入率、家用燃料的氡进入率和供水所致氡的进入率；A为空气氡的浓度；v为房间的体积；v1、v2、v3、v4分别为建筑材料、建筑物地基土壤、家用燃料和用水的体积。如果把B、C、D和E当常数。则通过对式(6)整理，我们可以得到以下微分方程的解：
其中：C0为常数，Bq；t为时间，h。
把${\lambda _{\text{v}}} = S{v_{\text{t}}}/v $代入式(7)我们可以得到：
参照UNSCEAR公布的最新各种室内氡源的相关参数及标准，并利用式(1)~(5)计算与实验中的对比，我们可以得到室外空气氡浓度A、建筑材料所致氡进入率B、建筑物地基土壤氡进入率C、家用燃料的氡进入率D和供水所致氡的进入率E的相对合理值分别为5、10、29、1和0.3 Bq/(m3·h)；氡的衰变常数λRn为7.6×10-3h-1。参照UNSCEAR公布的最新标准，在有地板且没有裂缝的房基土壤和建筑材料所致氡进入室内厚度通常在0.2 m[19]；初始可以直接从室内测的，房间的体积和房间通风道的表面积以及空气流动速率都可以从实际生活中得到。该模型不仅考虑到自然通风，也涉及到机械通风，由于机械通风可以近似看着一个常数，在微分中不影响结果，积分后加上该常数。因此，该模型具有预测的可行性。
方程式(7)、(8)是2个可以用来预测室内氡浓度变化的基本方程，要预测一间房的氡浓度，只要给出房间的长宽高、墙壁的厚度、通风量系数、用水量和家用燃料量以及空气交换流动率，就可以利用上述方程计算出该房间的氡浓度。如果房间在10 m2以上，可以不考虑土壤所致氡的来源[19]。
分析上述方程可以发现，室内氡的浓度除了受不同氡来源的影响外，房间的通风量系数S和空气流动速率vt，时间t和氡衰变常数λRn亦影响着室内氡浓度的估计值的准确性。参照UNSCEAR公布的最新各种室内氡进入率的相关参数及标准和相关文献[20-22]，令室外空气氡浓度A、建筑材料所致氡进入率B、建筑物地基土壤氡进入率C、家用燃料的氡进入率D和供水所致氡的进入率E分别为5、10、29、1、0.3 Bq/(m3·h)；氡的衰变常数λRn为7.6×10-3h-1。假定房子的长宽高分别为为10、10、3.5 m，即房间室内体积v为350 m3；空气的交换流动速率vt为185 m/h；室内氡的初始浓度为40 Bq/m3 ，即t=0时，CRn的值。代入式(8)计算得到理论值为(6.35～40 Bq/m3)与任天山和王玫等文献的实测结果(6～50 Bq/m3)一致，说明了该模型在计算室内氡浓度的可靠性。
预测室内氡的理论模型，几乎都没有把房间的通风量系数S考虑在他们的模型之中，而式(8)充分证明了房间的通风介质表面积大小是影响室内外空气交换的关键因子之一，更是改变室内氡浓度变化的核心因素之一，对室内氡浓度计算的精度影响大。例如，在不考虑室内氡浓度初始值的前提下，即不考虑式(8)常数C0，令t=0, 并且其他参数都不变，只改变房间通风量系数S，具体数据如表 1。则根据式(8)，可得到图 1。结合图 1与表 1，我们发现当室内完全不通风时，不考虑室内氡浓度的初始值的前提下，室内理论氡浓度CRn为475.2 Bq/m3；一旦通风室内氡浓度就快速下降，但当房间通风量系数S的值到达一定值时，室内氡浓度开始趋向于缓和，也可发现当S足够大时，室外氡浓度与室内浓度一致。因此，从模型中可以发现通风量系数S是降低室内氡浓度的重要手段，在生活中多开门窗能有效降低室内氡的浓度。
不同来源的氡进入室内后，室内222Rn随着时间积累, 同时由于放射性衰变和通风使之不断减少。因此，时间t是室内氡浓度变化的一个重要因子。假定其他因子不变，房间的通风量系数S为1 m2, 则根据方程(8)拟合，得到表 2、3和图 2结果如下。
从表 2和图 2中，我们发现当房间通风时，且通风量系数S和其它一定时，室内氡浓度CRn随时间快速降低，但当时间到达一定值(拟合值t为6 h时，室内氡浓度CRn斜率变小，开始趋向于0，随着t增大，室内氡浓度最终和室外空气氡浓度趋向一致，根据方程(8)的运算和表 2数据，当t≥22 h时，CRn趋向稳定值11.66，验证了上述理论。在生活中当房间的门窗面积一定的时候，可通过增加开窗的时间来降低室内氡的浓度。我们已计算出房间室内体积为350m3，通风量系数为1 m2时，通风6 h就可以使室内氡浓度与室外趋向一致；当S为2 m2时，通风4 h就能室内氡浓度与室外趋向一致；当S为5 m2时，不超过2 h就能到达一致。从表 3中，可以发现当通风量系数S不同时，室内氡浓度CRn趋向稳定值不同，如当S=2 m2时，CRn趋向于8.36 Bq/m3；而S为5时，CRn趋向于6.35 Bq/m3。所以，该模型为房子的通风设计提供了合理的依据。
将222Rn衰变常数考虑在模型的建立中十分必要，用表 2的数据对模型(8)拟合，结果显示(图 3)222Rn衰变常数对室内氡浓度的计算有很大的影响，特别是当通风量系数增大时，这种情况下更不容忽视。
前人的模型均是以室内不同氡源的析出面积与表面析出率来表述氡含量，再通过相关分析，获得室内氡的浓度。这样分析可能有一定的缺陷，由于不同的建筑物具有不同的建筑结构，如地板、墙壁的厚度也许就不一致。因此，其结果具有一定的误差，特别是涉及到家用燃料和供水时。而式(6)用不同室内氡源进入率与氡源的体积来构建定量分析，能更好的把房间地板和墙壁等的厚度以及家用燃料和供水充分考虑进去，从而减少计算误差。
在假定室内外氡浓度变化时保持均一混合的基础上，根据室内外空气流动平衡理论与氡自身衰变的特性，建立了一个非线性的室内氡浓度变化定量模型。该模型的建立，能比较好地反映房间的通风量系数S、时间t和室内氡浓度变化的关系，特别是在建筑设计过程中，利用该模型中房间室内体积v、房间通风量系数S以及时间t与室内氡浓度的关系，可设计出一个有效降低室内氡浓度的房间通风量系数。同时，提供了一个利用房间的通风量系数和空气流动速率来分析室内氡浓度变化的技术方法。模型模拟结果显示：室内氡的浓度变化不仅受不同氡来源含量的影响，还受空气流动速率和房间通风量系数、房间的体积、衰变时间等的影响。当房间的通风量系数一定时，可以通过延长通风时间来降低室内氡的浓度，且当通风时间达一定值时，室内外空气氡浓度趋向一致；当时间一定时，可扩大通风量系数S来减小室内氡的浓度，S大于一定值之后，室内外氡浓度趋向一致；并可以结合氡进入率与氡源的体积来进行预测室内氡浓度。模拟结果值与前人的实测结果存在一致性，这不仅说明建立室内氡浓度变化定量模型的可靠性，而且为调控室内氡浓度以及环境评价提供了技术支持。
2013, Vol. 35 Issue (5): 78-84
PDF
2013, Vol. 35 Issue (5): 78-84
PDF