国民经济的发展和能源的短缺促使我国以铀资源为基础的核电工业迅速崛起,在此大环境下,我国对铀矿的开采活动不断增强。而在开采过程中,铀矿废渣、采矿废水和工艺废水在雨水淋滤及扩散作用下不断迁移、浓缩、富集,并进入水体、土壤、大气等环境中,以废气、废水及废渣等形式对矿区及其周边地区产生具有长期性、持续性的放射性污染,危害人体健康[1-7]。
目前,国内外诸多专家学者对场地放射性污染调查与评价主要以放射性核素(226Ra、232Th、40K、238U等)、辐射贯穿剂量率以及总α、总β为研究重点。如Noguera等[8]对乌拉圭海岸带的天然放射性核素226Ra、232Th和40K进行调查;蒋文波等[9]对某退役铀矿周边土壤不同深度226Ra和238U特征进行研究;林武辉等[10]探究大气中γ辐射空气吸收剂量率的波动机制。场地放射性污染评价方法主要是辐射危害指数法[11]、镭当量活度[12]、美国环境保护局(US EPA)推荐的放射性污染场地风险评估模型、地累积指数法等。Asad等[13]运用辐射危害指数法和外部危险指数法对美国2个造船厂土壤的元素和放射性污染进行评估;余根锌[14]估算了我国福建省天然γ辐射剂量率,并量化福建省人居环境安全性标准;哈日巴拉等[15]运用US EPA建议的放射性健康风险模型评估我国巴彦乌拉铀矿周围饮用水中放射性和健康风险。
我国部分铀矿矿井、水冶生产设施于21世纪初开始相继停运,并对在线设备设施均进行了封存处理,但尾渣库、堆浸场、废石场等遗留场地有大量未中和的废渣和废石,还有周边被污染的公路、工业建(构)筑物及各类污染设备、管线等,均存在着一定的辐射环境污染风险。因此,为查清退役铀矿山遗留物对周边环境产生的污染影响及范围,本文以铀矿区及周边环境放射性源项为研究对象,对研究区内各污染源项进行调查及统计性分析,查明各放射性污染之间的相互联系;采用潜在生态危害指数法、US EPA推荐的放射性污染场地风险评估模型以及国际辐射防护委员会(ICRP)提供的年有效剂量模型[16]对研究区放射性污染强度范围进行调查;并针对特征放射性核素对人群可能造成的健康风险进行评价。
研究区位于广西壮族自治区某退役铀矿,地处亚热带季风带,具有典型的山地立体气候,年平均降雨量为1 773.1 mm,属中高山地貌,平均海拔1 500~1 600 m,岩性主要为中粗粒、中细粒似斑状黑云母花岗岩及二云母花岗岩,均为燕山期,矿区内无大型地表水体,年均流量约0.5 m3·s-1。
本次研究调查范围主要有坑井口、废石堆场、工业场地、堆浸渣场、尾矿库、建(构)筑物、设备器材和管线以及可能受到污染的地表水体和道路等。按照《环境核辐射监测规定》(GB 12379—1990)[17],按设施周围东北、西南、东南、西北方位布设,共布设采点19个。每个点在10 m×10 m范围内,用铁锹,移植镘刀等物去除其外围富含有机质的杂填土,采取0~10 cm的表层土,每次采取约1.0~2.0 kg并用塑封袋密封装入布制取样袋,排出空气压实在1 000 mL样品盒中,保证所取土壤样品接近场地样品源。采样点位置如表1所示。
表1 研究区采样点位置
Table 1 Location of sampling points in study area
序号Serial number名称Name区域Region位置(经纬度)Location (latitude and longitude)1老1堆废石场Old 1 pile waste rock quarry2老2堆废石场Old 2 pile waste rock quarry31392废石场1392 waste stone quarry41470废石场1470 waste stone quarry51490废石场1490 waste stone quarry61440废石场1440 waste stone quarry71480废石场1480 waste stone quarry81520废石场1520 waste stone quarry91560废石场1560 waste stone quarry废石场Waste rock quarry110°31'1.93″,26°9'4.39″110°30'59.46″,26°9'5.92″110°31'0.54″,26°9'2.24″110°30'58.46″,26°9'2.87″110°30'55.91″,26°9'3.28″110°31'23.56″,26°9'39.48″110°31'18.77″,26°9'37.74″110°31'18.70″,26°9'36.15″110°31'17.61″,26°9'33.93″10水冶工业场地Hydrometallurgical industrial site11水冶辅助工业场地Hydrometallurgical auxiliary industrial sites12破碎工业场地Crushing industrial sites13旧破碎工业场地Old broken industrial sites14炸药库Explosives depot15西部工棚区Western work shed area16生活办公区Living office area17三期库坝下工业场地Industrial site below the third phase of the reservoir dam工业场地Industrial sites110°31'7.50″,26°8'54.27″110°31'3.44″,26°8'55.03″110°31'5.41″,26°8'58.08″110°31'7.46″,26°9'0.44″110°31'10.16″,26°9'3.63″110°31'11.01″,26°9'1.62″110°31'12.36″,26°9'0.13″110°30'56.99″,26°8'41.72″18一期尾矿库Phase Ⅰ tailings storage19二期尾矿库Phase Ⅱ tailings storage尾矿库Tailings pond110°31'3.02″,26°8'51.70″110°30'59.85″,26°8'47.33″
按照《岩石样品中226Ra的分析方法射气法》(GB/T 13073—1991) [18]和《土壤中放射性核素的γ能谱分析方法》(GB/T 11743—2013)[19]中规定的相对比较法求解样品中238U、226Ra放射性核素比活度,并分析扩展不确定度。将土壤样品经100 ℃烘干、磨碎、过60目孔径筛。每个样品称取100 g放入γ能谱分析专用聚乙烯样品盒,用蜡密封后存放3周后,使放射系中铀、镭和它们的短寿命子体基本达到平衡至低本底γ能谱仪经刻度后,在相同条件下依次测量空样品盒、土壤粉末标准源和土壤样品能谱各8 h,利用仪器自带解谱软件分别求出238U、226Ra子体中γ射线发射几率比较大的92.6 KeV(238U)、609.4 KeV(226Ra)特征峰的全能峰净面积(将空样品盒测量值作为本底计数),根据测量时间,通过相对比活度法得到土壤238U、226Ra比活度,根据扩展不确定度公式求得高纯锗探测器测量的土壤样品238U、226Ra的扩展不确定度均分别<4.15%、<6.15%,满足GB/T 11743—2013中扩展不确定度铀<20%,镭<10%的要求。
按照《表面氡析出率测定》(EJ/T 979—1995)[20]采用真空法取样,取样时记录气温、气压、湿度、风度、风向等数据,用以对测量结果进行分析比较。采用积累法对研究区表面氡析出率进行计算,测量结果符合铀尾矿、废石场表面氡析出率测量差不超过10%的要求。研究区γ辐射吸收剂量率使用R-PD型γ辐射剂量率仪(福州智元仪器设备有限公司,中国)测量,依据《环境地表γ辐射剂量率测定规范》(GB/T 14583—1993)[21],测量时探头在离地面1 m高处进行现场定点测量,连续测量10次并读数,取测量结果算术平均值作为该点γ辐射剂量率测量值,γ辐射剂量率测量值总不确定度应不超过20%。
1.2.1 潜在生态风险指数法
潜在生态风险指数法(potential ecological risk index)最早由瑞典科学家Hakanson提出。以各种重金属在环境中所具有性质及行为特点为基础,加入沉积学理论权重比例,对土壤或沉积物中的重金属污染进行定量评价。该方法不仅兼顾重金属在土壤中浓度,而且涵盖了环境中多种元素协同作用、毒理学特性、污染范围及浓度以及环境对重金属污染的响应效应等因素,因此在环境风险评价中得到了广泛应用[22-23]。计算公式如下:
(1)
(2)
(3)
式中:为单个污染物的污染系数;
为单个污染物潜在生态风险系数;RI为多种污染物的综合潜在生态风险系数;Ci为污染物在0~1 m表层土壤中的实测值
为污染本底值(Bq·kg-1);Tri为污染物i的毒性响应系数。潜在生态风险评价等级如表2所示。
表2 潜在生态风险评价等级
Table 2 Potential ecological risk evaluation level
单个重金属潜在生态风险系数Individual heavy metal potential ecological risk factor单个重金属潜在生态风险程度Degree of potential ecological risk of individual heavy metal多种重金属潜在生态风险指数Multiple heavy metals potential ecological risk index风险等级Risk levelEir<40轻微MinorRI<150轻微Minor40≤Eir<80中等Moderate150≤RI<300中等Moderate80≤Eir<160强Strong300≤RI<600强Strong160≤Eir<320很强Very strong600≤RI很强Very strong320≤Eir极强Extremely strong
1.2.2 放射性风险评价
226Ra、238U作为典型的放射性物质,其衰变过程中产生大量电离辐射,生物体长期暴露在高强度电离辐射中,会导致细胞突变,进而引发癌变以及死亡。故对研究区内226Ra、238U进行风险评估,采用US EPA推荐的放射性污染场地风险评估模型进行健康风险评估[24],其基本公式为:
R=∑(RE×SF)
(4)
式中:R为健康风险值,RE为某一暴露途径的暴露量,SF为某一暴露途径所对应的斜率因子。因此,通过计算放射性物质每种暴露途径的暴露量,再结合相应暴露途径的斜率因子即可得到该污染物质所造成的健康风险值。本文调查对象为某退役铀矿场地,考虑到调查结果和风险评估主要为污染退役治理提供依据,以外照射、污染土壤摄入和吸入3种暴露途径计算其健康风险,其暴露量计算公式为:
RE外照射=C土壤×EF×ED×ET×GSF外照射×
(5)
式中:RE外照射为外照射暴露量(Bq·a·g-1);C土壤为放射性污染物在表层土中的比活度(Bq·g-1);EF为暴露频率(d·a-1);ED为暴露期(a);ET为暴露时间(h·d-1);t为持续时间(a);λ为衰变常数(a -1);GSF外照射为γ屏蔽系数(室外取1.0);ACF外照射为面积修正系数。
(6)
式中:RE摄入为土壤摄入暴露量(Bq);IRS为每日摄入土壤量(mg·d-1)。
(7)
式中:RE吸入为土壤吸入暴露量(Bq);IRA为每日呼吸量(m3·d-1);PEF为颗粒物排放系数(m3·kg-1)。健康风险计算时,土壤中放射性污染物的比活度取研究区中各点位226 Ra、238U的表层土实测值,衰变常数λ以λ=0.693/T1/2(T1/2为半衰期)计算。上述各参数的取值列于表3。
表3 健康风险评价参数取值[25-26]
Table 3 The values of health risk evaluation parameters [25-26]
参数Parameters单位Units数值Numerical value暴露期(ED)Exposure period (ED)a25暴露频率(EF)Exposure frequency (EF)d·a-1250暴露时间(ET)Exposure time (ET)h·d-18成人每日摄入土壤量(IRS)Daily ingestion of soil for adults (IRS)mg·d-1100成人每日呼吸量(IRA)Adult daily respiratory volume (IRA)m3·d-114.5γ屏蔽系数(GSF外照射)γ shielding factor (GSFext)量纲为1Dimensionless1.0(表层Surface)颗粒物排放系数(PEF)Particulate emission factor (PEF)m3·kg-12.97×1011238U面积修正系数(ACF外照射)238U area correction factor (ACFext)量纲为1Dimensionless0.99226Ra面积修正系数(ACF外照射)226Ra area correction factor (ACFext)量纲为1Dimensionless0.53238U衰变常数(λ)238U decay constant (λ)a-11.55×10-10226Ra衰变常数(λ)226Ra decay constant (λ)a-14.33×10-4238U外照射斜率因子(SF外照射)238U external irradiation slope factor (SFext)g· Bq-1·a-13.35×10-9226Ra外照射斜率因子(SF外照射)226Ra external irradiation slope factor (SFext)g· Bq-1·a-16.75×10-7238U摄入斜率因子(SF摄入)238U ingestion slope factor (SFin)Bq-11.26×10-9226Ra摄入斜率因子(SF摄入)226Ra ingestion slope factor (SFin)Bq-17.96×10-9238U吸入射斜率因子(SF吸入)238U inhalation shot slope factor (SFinh)Bq-16.39×10-7226Ra吸入射斜率因子(SF吸入)226Ra inhalation shot slope factor (SFinh)Bq-17.61×10-7
1.2.3 辐射年有效剂量与终生癌症风险
使用空气中的吸收剂量率与有效剂量的转换系数(0.7 Sv)可计算出室外(废石场、工业场地、尾矿库及周边区域)和室内(住房、厂房)年有效剂量当量。室外和室内的占用系数分别为0.2和0.8[27-28]。辐射年有效剂量当量(E)是由以下关系计算出来的:
E室外(mSv·a-1)=D(nGy·h-1)×24×365×0.2×0.7×10-6
(8)
E室内(mSv·a-1)=D(nGy·h-1)×24×365×0.8×0.7×10-6
(9)
式中:D为测量的空气中的吸收剂量率平均值(nGy·h-1)。空气吸收剂量率对人口造成的年有效剂量当量是通过将室外和室内的年有效剂量当量相加来计算,而其对人口的影响,通常用集体有效剂量(SC)来表征,SC由公式(10)[29]来估计:
SC=E×N(P)
(10)
式中:N(P)是研究区周边人口。
一生中的终身癌症风险是指在一定的暴露水平下,一生中发生癌症的概率。它的计算是假设剂量和随机效应之间存在线性关系,它以一个数值表示,在一定数量的人以一定的剂量暴露于致癌物中,预计额外的癌症发生概率。终身癌症风险(ELCR)是用公式(11)来计算的:
ELCR=E×DL×RF
(11)
式中:E为年总有效剂量当量,DL为人均寿命(据第七次全国人口普查,中国人均寿命为77.3 a),RF为风险系数(Sv-1),即单位剂量当量的致命癌症风险。对于随机效应,ICRP使用的数值为0.05[30]。
根据铀尾矿库氡析出率平均值,可估算铀矿各区域的氡年释放量。计算公式如下:
(12)
式中:RTotal表示待测区域的氡年释放量(Bq);JAV平均值表示待测区域的氡析出率平均值(Bq·m-2·s-1);S为待测区域的表面积(m2);t为估算时间,估算1年释放量时间t=3.1536×107 s。
在对研究区放射性源项调查中,238U比活度595.5~2 335 Bq·kg-1,226Ra比活度637~16 551 Bq·kg-1,均高于全国、世界平均水平[31],分别为全国、世界本底值的15.08倍~59.12倍、17.45倍~453.45倍和14.89倍~58.38倍、15.93倍~413.78倍;γ辐射空气吸收剂量率分别为全国、世界平均水平[32]的11.91倍、14.05倍;氡析出率范围处于全国范围的0.34倍~6.62倍(表4)。由图1可知,238U的样本特征基本服从正态分布,无异常值,平均值1 454.2 Bq·kg-1位于中位数上方,属右偏态分布,变异系数0.38为中等变异,表明样本梯度离散程度较大。226Ra比活度、γ辐射空气吸收剂量率、氡析出率样本均不服从正态分布,均值分别为2 962.3 Bq·kg-1、1 179.9 nGy·h-1和1.34 Bq·m-2·s-1,与中位数相差较大,变异系数分别为1.51、0.82和0.90,样本离散程度较大,偏度及峰度均较高,表现出明显的不对称性,分布存在异常值,分别位于一期、二期尾矿库。由于尾矿库表面覆盖厚度较小,氡析出率可随温度、气压变化或降水而有很大的变化[34],当大气压力减少时,氡析出率增加;降雨会增加尾矿库的湿度,在一定的湿度范围内尾矿库表面的氡析出率随着温度的增加而增加。氡作为镭衰变的产物之一,在其衰变过程中会产生电离辐射,使空气辐射剂量率测量值升高,故而这三者在样本特征分布形式上有一定相似性。
图1 放射性源项柱状图和Q-Q图
注:A为老1堆废石场,B为老2堆废石场,C为1392废石场,D为1470废石场,E为1490废石场,F为1440废石场,G为1480废石场,H为1520废石场,I为1560废石场,J为水冶工业场地,K为水冶辅助工业场地,L为破碎工业场地,M为旧破碎工业场地,N为炸药库,O为西部工棚区,P为生活办公区,Q为三期库坝下工业场地,R为一期尾矿库,S为二期尾矿库。
Fig. 1 Histogram of radioactive source term and Q-Q diagram
Note: A represents old 1-pile waste rock site; B represents old 2-pile waste rock site; C represents 1392 waste rock site; D represents 1470 waste rock site; E represents 1490 waste rock site; F represents 1440 waste rock site; G represents 1480 waste rock site; H represents 1520 waste rock site; I represents 1560 waste rock site; J represents hydrometallurgical industrial site; K represents hydrometallurgical auxiliary industrial site; L represents crushing industrial site; M represents old crushing industrial site; N represents explosives depot; O represents western shed area; P represents living office area; Q represents industrial site below the Phase Ⅲ dam; R represents Phase Ⅰ tailings pond; S represents Phase Ⅱ tailings pond.
表4 放射性源项背景值
Table 4 Background values of radioactive source terms
区域Region238U/(Bq·kg-1)[31]226Ra/(Bq·kg-1)[31]γ辐射空气吸收剂量率/(nGy·h-1)[32]Air absorbed dose rate of γradiation/(nGy·h-1)[32]氡析出率/(Bq·m-2·s-1)[31]Radon precipitation rate /(Bq·m-2·s-1)[31]范围Scope均值Mean范围Scope均值Mean范围Scope均值Mean范围Scope均值Mean全国 National1.80~520.039.502.40~134.6036.50-99.100.74世界World40.0040.00-84.00-研究区Research area595.50~2 335.101 454.20637~16 5512 962.30225~4 0221 179.900.25~4.901.34
注:-代表尚无统一标准值。
Note: - means there is no uniform standard value yet.
放射性核素随人类活动在环境中的迁移范围也越来越广,因此区域放射性背景值也随之增高,造成环境放射性污染,影响人类的健康安全[35]。所以,研究各放射性源项之间相互联系,对矿区退役治理工程建设具有重要意义。分别用238U、226Ra比活度对氡析出率和γ辐射空气吸收剂量率作线性拟合图(图2),如图2所示,用238U比活度对氡析出率和γ辐射空气吸收剂量率进行拟合,可决系数分别为0.27和0.19,用226Ra比活度对氡析出率、γ辐射空气吸收剂量率拟合的可决系数为0.79和0.85,可以看出226Ra比活度对放射性源项的拟合度高于238U。在放射性物质的递次衰变中,母体原子数目不断减少,子体原子数目不断增加,在相当长的时间后,母体与子体的放射性比活度比例趋于稳定且不随时间的推移而产生变化,核素达到长期平衡的状态。238U、226Ra的长期平衡往往需要300万年才能建立[36],故而氡的析出很大程度上受它的母体226Ra的比活度控制,与238U比活度相关性不大。而226Ra和氡的比活度与γ辐射空气吸收剂量率成正比,所以这3种放射性源项表现出较好的线性关系。
图2 放射性源项线性拟合
Fig. 2 Linear fit of the radioactive source term
场地放射性源项相关性如图3所示,可进一步说明研究区各放射性源项相互联系。在0.01显著水平下,238U比活度分别与226Ra比活度、氡析出率、空气吸收辐射剂量率的相关系数为0.35、0.52和0.44,呈现出正相关、弱相关性;226Ra比活度分别与氡析出率、空气吸收辐射剂量率相关系数为0.89、0.92,氡析出率与空气吸收辐射剂量率相关系数为0.93,均达到高度相关,再次说明相比于238U比活度的影响,226Ra比活度对氡析出和空气吸收辐射剂量的产生具有重要作用。
图3 放射性相关系数
Fig. 3 Radiological correlation coefficient
潜在生态风险指数法作为环境土壤评价的重要方法,许多学者将其引入放射性土壤的评价体系中[37]。根据《电离辐射防护与辐射源安全基本标准》[38](GB 18871—2002)226Ra、238U分别为极毒、低毒组别对其毒性响应系数进行赋值。研究区场地各区域采样点范围内238U与226Ra所造成的潜在生态风险评价结果列于表5。潜在生态风险指数法评价结果表明,研究区各采样点均<40,属于轻微潜在生态风险;研究区226Ra潜在生态风险程度较高,不存在轻微潜在生态风险点位,中等潜在生态风险占比42.11%,强潜在生态风险占比36.83%,很强潜在生态风险占比10.53%,一期尾矿库、二期尾矿库
分别为1 190.12、1 025.01,属于极强潜在风险程度,占比10.53%。RI总风险结果表明,老1堆废石场、老2堆废石场、1470废石场、1440废石场、水冶辅助工业场地、炸药库、西部工棚区、生活办公区、三期库坝下工业场地点位属于轻微潜在生态风险,占比47.37%;1392废石场、1480废石场、1480废石场、1520废石场、1560废石场、水冶工业场地、破碎工业场地、旧破碎工业场地号点位属于中等潜在生态风险,占比42.10%;位于尾矿库区域一期尾矿库、二期尾矿库点属于很强潜在生态风险,占比10.53%。从评价结果来看,研究区中等潜在生态风险地区占比最大,且
与RI结果具有相似性,表明226Ra对总风险的贡献率较高。
表5 场地潜在生态风险评价结果
Table 5 Results of potential ecological risk evaluation of the site
采样点Sampling pointsEir(238U)Eir(226Ra)RI风险等级Risk level老1堆废石场Old 1 pile waste rock quarry16.47 74.11 90.59 轻微Minor老2堆废石场Old 2 pile waste rock quarry9.78 44.33 54.11 轻微Minor1392废石场1392 waste stone quarry27.59 122.41 150.00 中等Moderate1470废石场1470 waste stone quarry17.06 76.83 93.89 轻微Minor1490废石场1490 waste stone quarry36.11 161.86 197.98 中等Moderate1440废石场1440 waste stone quarry16.07 71.12 87.19 轻微Minor1480废石场1480 waste stone quarry33.68 147.88 181.55 中等Moderate1520废石场1520 waste stone quarry37.13 167.64 204.77 中等Moderate1560废石场1560 waste stone quarry33.47 146.00 179.47 中等Moderate水冶工业场地Hydrometallurgical industrial site27.79 125.19 152.98 中等Moderate水冶辅助工业场地Hydrometallurgical auxiliary industrial sites16.53 74.88 91.41 轻微Minor破碎工业场地Crushing industrial sites24.55 110.09 134.64 中等Moderate旧破碎工业场地Old broken industrial sites32.26 143.28 175.54 中等Moderate炸药库Explosives depot15.52 68.27 83.79 轻微Minor西部工棚区Western work shed area13.41 60.47 73.88 轻微Minor生活办公区Living office area15.13 67.36 82.50 轻微Minor三期库坝下工业场地Industrial site below the third phase of the reservoir dam22.11 98.96 121.07 轻微Minor一期尾矿库Phase Ⅰ tailings storage38.34 1 151.77 1 190.12 很强Very strong二期尾矿库Phase Ⅱ tailings storage20.69 1 004.31 1 025.01 很强Very strong
研究表明,放射性辐射可对生物体产生不可逆转的伤害,一但人体暴露在放射性环境中并达到一定辐射剂量,轻则引起恶心、呕吐、脱发等症状,严重则会导致细胞癌变等,最终造成休克甚至死亡。本文使用US EPA推荐的健康风险模型对场地放射性核素对人体所造成的危害进行评估,评价结果列于表6。结果表明,场地内226Ra造成的健康风险值相较于238U高出1~2个数量级,说明该场地主要风险来源于226Ra,最大风险点位位于尾矿库范围,最大风险值为1.16×10-4,是US EPA建议最大风险(10-4)[39]的1.16倍,238U最大风险值为1.87×10-6,处于可接受范围内。其中,以外照射途径造成的226Ra最大风险值为3.36×10-5,占比28.97%,摄入途径最大风险值为8.19×10-5,占比70.6%,吸入途径最大风险值低于10-10,可忽略,故而摄入是最主要的暴露途径。矿区工作人员作为场地内主要潜在暴露人群,在工作时间内均受到不同途径的辐射暴露,以每天工作8 h条件下估算每年潜在暴露剂量。结果表明,场地226Ra对工作人员以外照射、摄入、吸入途径造成最大年暴露剂量分别为49.78、10 292.65和1.68×10-3 Bq(外照射转化系数取辐射权重因子值1.0),最大总暴露量10 342.43 Bq;238U的3种途径最大暴露量分别为13.12、1 542.11和2.36×10-4 Bq,最大总暴露量1 465.23 Bq,进一步表明226Ra是场地内主要的风险来源。在铀矿开采和挖掘隧道等工程中,将原本在地下,对人类生活环境没有影响的放射性核素如238U、226Ra提升至地面,进入人的生活环境。在矿区人员进行工作时,采取了一定的辐射防护,但依旧会以各种途径收到放射性辐射污染,它不仅影响到矿工工人,也为矿区周围的公众环境带来额外的放射性辐射照射。矿区周边经放射性核素污染的花岗岩被居民用作建筑耗材,加大了放射性污染的风险程度。
表6 场地放射性核素健康风险评价结果
Table 6 Results of health risk evaluation of radionuclides at the site
核素Nuclide比活度范围/(Bq·kg-1)Activity range/(Bq·kg-1)暴露途径Exposure pathways暴露量/BqExposure/Bq最大总暴露量/BqMaximum total exposure/Bq风险值Risk value最大总风险Maximum total risk226Ra595.5~2 335.1外照射External irradiation1.92~49.78摄入Ingestion396.13~10 292.65吸入Inhalation6.45×10-5~1.68×10-310 342.431.29×10-6~3.36×10-53.15×10-6~8.19×10-54.91×10-13~1.28×10-111.16×10-4238U637~16 551外照射External irradiation3.35~13.12摄入Ingestion370.33~1 452.11吸入Inhalation6.03×10-5~2.36×10-41 465.231.12×10-8~4.39×10-84.67×10-7~1.83×10-63.85×10-11~1.51×10-101.87×10-6
本文采用ICRP推荐的估算模型对研究区空气辐射年有效剂量进行计算并评估其终身癌症风险(ELCR),对场地氡析出率做量化处理,计算其年释放量,为进一步治理铀矿区退役工程提供参考。表7、表8分别为研究区年有效剂量和氡年释放量评价结果。结果显示,研究区室外有效剂量高于室内,分别为1.10~19.73、0.28~4.93 mSv,远高于世界平均值(室外0.072 mSv,室内0.412 mSv),研究区最大年有效剂量达24.66 mSv,是世界平均水平年有效剂量(0.484 mSv)的50.9倍,最大集体有效剂量为4 476.32 人·Sv。终生癌症风险是指在一定的暴露水平下终生患癌的概率,研究区最大终生癌症风险为9.53×10-2,是世界平均水平(1.694×10-2)的5.64倍。据计算,研究区废石场区域面积约为32 796 m2,氡年释放量1.35×1012 Bq,工业场地面积约为47 044 m2,氡年释放量1.07×1012 Bq,尾矿库面积约36 503 m2,氡年释放量4.57×1012 Bq。从计算结果来看,尾矿库依旧是空气辐射与氡释放的主要污染场地。研究表明,大陆气团影响下的γ辐射空气吸收剂量率普遍高于海洋气团影响下的结果,大陆气团与陆地下垫面长时间接触,进而拥有更高浓度的陆地土壤中226Ra衰变而来的222Rn及其子体,而海洋气团具有较高的湿度,相应的222Rn及其子体浓度较低[38]。部分研究将γ辐射空气吸收剂量率夏低冬高的季节特征归因为夏季更强的降雨清除效应[40]。并且在小时尺度下白天温度上升,由于蒸腾作用,有利于土壤中水汽携带222Rn上升至大气并释放,导致γ辐射空气吸收剂量率随之增加;温度变化也与大气稳定度密切相关,即白天温度上升促进垂向对流,而夜晚温度下降不利于垂向对流。强对流有利于222Rn及其子体核素的向上稀释扩散,并导致空气中放射性核素浓度降低,以及氡释放量的下降。
表7 年有效剂量评价结果
Table 7 Annual effective dose evaluation results
点位Point location室外年有效剂量/mSvOutdoor annual effective dose/mSv室内年有效剂量/mSvIndoor annual effective dose/mSv年有效剂量/mSvAnnual effective dose/mSv集体有效剂量 /(人·Sv)Collective effective dose/(man·Sv)终生风险(×10-2)Lifetime risk (×10-2)老1堆废石场Old 1 pile waste rock quarry1.244.966.201 125.202.40老2堆废石场Old 2 pile waste rock quarry0.712.833.54642.181.371392废石场1392 waste stone quarry1.355.426.771 228.712.621470废石场1470 waste stone quarry0.913.644.54824.701.761490废石场1490 waste stone quarry0.923.694.62838.061.781440废石场1440 waste stone quarry1.305.206.501 179.742.511480废石场1480 waste stone quarry1.255.016.261 136.332.421520废石场1520 waste stone quarry2.6510.5913.242 402.885.121560废石场1560 waste stone quarry1.817.229.031 638.273.49水冶工业场地Hydrometallurgical industrial site0.873.494.37792.431.69水冶辅助工业场地Hydrometallurgical auxiliary industrial sites0.391.551.93350.580.75破碎工业场地Crushing industrial sites1.465.837.281 322.192.82旧破碎工业场地Old broken industrial sites1.154.605.751 042.842.22炸药库Explosives depot0.702.803.50634.391.35西部工棚区Western work shed area0.552.212.76500.831.07生活办公区Living office area0.281.101.38250.420.53三期库坝下工业场地Industrial site below the third phase of the reservoir dam0.351.421.77321.640.68一期尾矿库Phase Ⅰ tailings storage4.6818.7123.384 243.719.04二期尾矿库Phase Ⅱ tailings storage4.9319.7324.664 476.329.53
表8 氡年释放量评价结果
Table 8 Evaluation results of annual radon release
区域Region面积/m2Area/m2氡析出率平均值/(Bq·m-2·s-1)Average radon precipitation rate/(Bq·m-2·s-1)氡年释放量/BqAnnual release of radon/Bq废石堆场Waste rock dumps32 7961.311.35×1012工业场地Industrial sites47 0440.721.07×1012尾矿库Tailings pond36 5033.974.57×1012
研究区放射性源项调查数据显示,238U比活度595.5~2 335.1 Bq·kg-1,226Ra比活度637~16 551 Bq·kg-1,分别为全国、世界本底值的15.08倍~59.12倍、17.45倍~453.45倍和14.89倍~58.38倍、15.93倍~413.78倍;γ辐射空气吸收剂量率分别为全国、世界平均水平的11.91倍、14.05倍;氡析出率范围处于全国范围的0.34倍~6.62倍,最大超标点位均位于尾矿库区域。
238U比活度对氡析出率和γ辐射空气吸收剂量率拟合可决系数分别为0.27、0.19,226Ra比活度对氡析出率、γ辐射空气吸收剂量率可决系数为0.79和0.85,226Ra比活度对放射性源项的拟合度高于238U。238U比活度与226Ra比活度、氡析出率、空气吸收辐射剂量率呈现弱相关性(相关系数为0.35、0.52和0.44),226Ra比活度、氡析出率、空气吸收辐射剂量率三者之间达到高度相关(相关系数为0.89、0.92和0.93)。
潜在生态风险指数法结果表明强潜在生态风险点位占比最大,226Ra对总风险的贡献率较高;US EPA健康风险评价结果显示226Ra是主要的健康风险来源,摄入途径是最主要的暴露途径。研究区空气辐射最大年有效剂量达24.66 mSv,是世界平均水平年有效剂量的50.9倍,最大集体有效剂量为4 476.32 人·Sv,最大终生癌症风险为9.53×10-2,是世界平均水平的5.64倍。研究区最大氡年释放量为4.57×1012 Bq,影响面积约36 503 m2,各评价方法均表明尾矿库周边区域是污染最为严重地区。
[1] 王津, 李红春, 刘娟, 等. 铀矿开采利用过程中的放射性污染研究进展[J]. 环境与健康杂志, 2013, 30(11): 1033-1036
Wang J, Li H C, Liu J, et al. Progress on radioactive pollution and related problems caused by uranium mining [J]. Journal of Environment and Health, 2013, 30(11): 1033-1036 (in Chinese)
[2] 刘娟, 李红春, 王津, 等. 华南某铀矿开采利用对地表水环境质量的影响[J]. 环境化学, 2012, 31(7): 981-989
Liu J, Li H C, Wang J, et al. Environmental quality of surface waters in a uranium industrial site in South China [J]. Environmental Chemistry, 2012, 31(7): 981-989 (in Chinese)
[3] 朱琳, 王斐. 环境中放射性污染的来源、特点、危害和防护[J]. 河南科技, 2014(13): 185
Zhu L, Wang F. Sources, characteristics, hazards and protection of radioactive pollution in the environment [J]. Journal of Henan Science and Technology, 2014(13): 185 (in Chinese)
[4] 徐园, 孔海宇, 王希涛, 等. 水体环境中的放射性污染及测量[J]. 核安全, 2020, 19(1): 75-79
Xu Y, Kong H Y, Wang X T, et al. Radioactive pollution in water environment and its measurement [J]. Nuclear Safety, 2020, 19(1): 75-79 (in Chinese)
[5] 汪江英, 朱建林, 莫子奋, 等. 土壤中铀的赋存形态及放射性污染治理研究进展[J]. 能源研究与管理, 2020(4): 65-69
Wang J Y, Zhu J L, Mo Z F, et al. Advances in research on occurrence of uranium in soil and radioactive pollution control [J]. Energy Research and Management, 2020(4): 65-69 (in Chinese)
[6] 刘畅荣, 王志勇, 崔娜. 铀矿区大气放射性污染的健康风险评价研究展望[J]. 工业设计, 2016(8): 180-181
Liu C R, Wang Z Y, Cui N. Prospect of health risk assessment of atmospheric radioactive pollution in uranium mining areas [J]. Industrial Design, 2016(8): 180-181 (in Chinese)
[7] Hossain K A, Kumar G A, Sumon R M, et al. An investigation on the possible radioactive contamination of environment during a steam-line break accident in a VVER-1200 nuclear power plant [J]. Current World Environment, 2019, 14(2): 299-311
[8] Noguera A, Pereira H B, Fornaro L. Natural radionuclide survey in the coastal strip of the 290 Ramsar site, Uruguay [J]. Environmental Earth Sciences, 2018, 77(22): 755
[9] 蒋文波, 高柏, 张海阳, 等. 某铀矿区周边土壤238U和226Ra分布特征及污染评价[J]. 中国环境科学, 2021, 41(4): 1799-1805
Jiang W B, Gao B, Zhang H Y, et al. Distribution characteristics and pollution assessment of 238U and 226Ra in soils surrounding a uranium mining area [J]. China Environmental Science, 2021, 41(4): 1799-1805 (in Chinese)
[10] 林武辉, 王诗玥, 黄亚萍, 等. 大气中γ辐射空气吸收剂量率的波动机制[J]. 中国环境科学, 2022, 42(3): 1097-1103
Lin W H, Wang S Y, Huang Y P, et al. Mechanism of the variable γ radiation air absorbed dose rate in the atmosphere [J]. China Environmental Science, 2022, 42(3): 1097-1103 (in Chinese)
[11] Suresh S, Rangaswamy D R, Sannappa J, et al. Estimation of natural radioactivity and assessment of radiation hazard indices in soil samples of Uttara Kannada District, Karnataka, India [J]. Journal of Radioanalytical and Nuclear Chemistry, 2022, 331(4): 1869-1879
[12] Beretka J, Matthew P J. Natural radioactivity of Australian building materials, industrial wastes and by-products [J]. Health Physics, 1985, 48(1): 87-95
[13] Asad S. Investigation of elemental and radiological contamination of soils in two shipyards in Chittagong, Bangladesh [J]. Radiochimica Acta International Journal for Chemical Aspects of Nuclear Science &Technology, 2014, 102(8): 741-749
[14] 余根锌. 福建省环境地表γ辐射剂量率估算与人居环境安全性评价[J]. 物探与化探, 2021, 45(1): 192-199
Yu G X. The environmental terrestrial gamma-radiation dose rate estimation and the living environment safety evaluation in Fujian Province [J]. Geophysical and Geochemical Exploration, 2021, 45(1): 192-199 (in Chinese)
[15] 哈日巴拉, 拓飞, 胡碧涛, 等. 巴彦乌拉铀矿周围饮用水中放射性及健康风险[J]. 中国环境科学, 2017, 37(9): 3583-3590
Haribala, Tuo F, Hu B T, et al. Radioactivity and health risk of drinking water around Bayanwula uranium mining area [J]. China Environmental Science, 2017, 37(9): 3583-3590 (in Chinese)
[16] Chandrasekaran A, Senthilkumar G, Thillaivelavan K, et al. Spatial distribution and lifetime cancer risk due to gamma radioactivity in Yelagiri Hills, Tamilnadu, India [J]. Egyptian Journal of Basic and Applied Sciences, 2014, 1: 38-48
[17] 国家技术监督局. 环境核辐射监测规定: GB 12379—1990[S]. 北京: 中国标准出版社, 1990
[18] 国家技术监督局. 岩石样品中226Ra的分析方法射气法: GB/T 13073—1991[S]. 北京: 中国标准出版社, 1992
[19] 国家质量监督检验检疫总局, 中国国家标准化管理委员会. 土壤中放射性核素的γ能谱分析方法: GB/T 11743—2013[S]. 北京: 中国标准出版社, 2014
[20] 中国核工业总公司. 表面氡析出率测定: EJ/T 979—1995[S]. 北京: 中国核工业总公司, 1995
[21] 国家环境保护局, 国家技术监督局. 环境地表γ辐射剂量率测定规范: GB/T 14583—1993[S]. 北京: 中国标准出版社, 1993
[22] 王漫漫. 太湖流域典型河流重金属风险评估及来源解析[D]. 南京: 南京大学, 2016: 6-7
Wang M M. Risk assessment and source apportionment of heavy metals in typical rivers of Taihu Basin [D]. Nanjing: Nanjing University, 2016: 6-7 (in Chinese)
[23] 阿卜杜萨拉木·阿布都加帕尔, 王宏卫, 杨胜天, 等. 准东地区土壤重金属污染特征及来源分析[J]. 中国矿业, 2019, 28(11): 168-174
Abdusalam A, Wang H W, Yang S T, et al. Pollution characteristics and source analysis of soil heavy metals in East Junggar region [J]. China Mining Magazine, 2019, 28(11): 168-174 (in Chinese)
[24] 李战国, 朱勇兵, 李军, 等. 226Ra污染场地调查与健康风险评估[J]. 安全与环境学报, 2020, 20(6): 2408-2415
Li Z G, Zhu Y B, Li J, et al. Investigation and health risk assessment of 226Ra in deserted contaminated site [J]. Journal of Safety and Environment, 2020, 20(6): 2408-2415 (in Chinese)
[25] U.S. Environmental Protection Agency (US EPA). RAIS risk exposure models for radionuclides user’s guide [EB/OL]. (2019-07-15) [2022-10-12]. https://rais.Ornl.gov/tools/rais_rad _risk_guide.html.
[26] Oak Ridge National Laboratory. Area correction factors for contaminated soil for use in risk and dose assessment models [EB/OL]. [2022-10-12]. https://rais.ornl.gov /tools/ACF_FINAL.Pdf,2014.
[27] United Nations Scientific Committee on the Effects of Atomic Radiations (UNSCEAR). 1993 Sources, effects and risks of ionizing radiation [R]. New York: UNSCEAR, 1993
[28] United Nations Scientific Committee on the Effects of Atomic Radiations (UNSCEAR). 2000 Sources and effects of ionizing radiation. UNSCEAR Report on the general assembly, with scientific annex [R]. New York: UNSCEAR, 2000
[29] Abdurabu W A, Ramli A T, Saleh M A, et al. Terrestrial gamma dose rate, radioactivity and radiological hazards in the rocks of an elevated radiation background in Juban District, Ad Dali’ Governorate, Yemen [J]. Journal of Radiological Protection: Official Journal of the Society for Radiological Protection, 2016, 36(1): 163-177
[30] International Commission on Radiological Protection (ICRP). ICRP Publication 60: 1990 Recommendations of the International Commission on Radiological Protection. Annual of ICRP Reference [R]. Ottawa: ICRP, 1991
[31] 全国环境天然放射性水平调查总结报告编写小组. 全国土壤中天然放射性核素含量调查研究(1983—1990年)[J]. 辐射防护, 1992, 12(2): 122-142
The Writing Group for the Summary Report on Nationwide Survey of Environmental Radioactivity Level in China. Investigation of natural radionuclide contents in soil in China [J]. Radialization Protection, 1992, 12(2): 122-142 (in Chinese)
[32] 潘自强, 修炳林. 联合国原子辐射效应科学委员会及其第49次会议[J]. 辐射防护, 2000, 20(5): 312-317, 320
Pan Z Q, Xiu B L. United nations scientific committee on the effects of atomic radiation (UNSCEAR) and its forty-ninth session [J]. Radialization Protection, 2000, 20(5): 312-317, 320 (in Chinese)
[33] Wang Z Y. Natural radiation environment in China [J]. International Congress Series, 2002, 1225: 39-46
[34] 肖高平. 铀尾矿库氡析出率变化规律的研究及应用[D]. 衡阳: 南华大学, 2016: 31-45
Xiao G P. Study on the variation characteristics of radon exhalation rate and its application in uranium tailings reservoir [D]. Hengyang: University of South China, 2016: 31-45 (in Chinese)
[35] 朱国兴. 西南某磷矿周边环境放射性水平研究[D]. 北京: 中国地质大学(北京), 2018: 21-30
Zhu G X. Study on the radioactivity level of the surrounding environment of a phosphate rock in southwest China [D]. Beijing: China University of Geosciences, 2018: 21-30 (in Chinese)
[36] 刘成, 杨亚新, 彭道峰, 等. 铀矿山废石堆表面氡析出率的变化规律[J]. 东华理工大学学报(自然科学版), 2013, 36(4): 385-389
Liu C, Yang Y X, Peng D F, et al. Change regulation of the radon exhalation rate in the waste dump of uranium mine [J]. Journal of East China Institute of Technology (Natural Science), 2013, 36(4): 385-389 (in Chinese)
[37] 李冠超, 蒋文波, 钟丽艳, 等. 某钽铌厂及周边土壤放射性水平调查与健康风险评估[J]. 有色金属(冶炼部分), 2021(12): 112-119
Li G C, Jiang W B, Zhong L Y, et al. Investigation and health risk assessment of soil radioactivity level in tantalum-niobium plant and its surroundings [J]. Nonferrous Metals (Extractive Metallurgy), 2021(12): 112-119 (in Chinese)
[38] 国家质量监督检验检疫总局. 电离辐射防护与辐射源安全基本标准: GB 18871—2002[S]. 北京: 中国标准出版社, 2004
[39] Pawel D, Boyd M. Studies of radiation health effects inform EPA actions [J]. Journal of Radiological Protection: Official Journal of the Society for Radiological Protection, 2019, 39(4): S40-S57
[40] 谢志辉, 刘春琼, 史凯. 成都市PM10对大气辐射环境的影响[J]. 环境科学研究, 2016, 29(7): 972-977
Xie Z H, Liu C Q, Shi K. Impacts of PM10 on the radiation environments of the atmosphere in Chengdu City [J]. Research of Environmental Sciences, 2016, 29(7): 972-977 (in Chinese)