Distribution and Health Risk Assessment of Cr6+ in Soil and Groundwater of a Chemical Plant
-
摘要:
为掌握某化工厂Cr6+污染状况及风险水平,在研究区布设土壤与地下水采样点各19个,测定Cr6+含量,分析污染特征与成因,并开展健康风险评估。结果表明:表层土壤(0~0.5 m)Cr6+超标率为42.11%,其浓度随深度增加总体上降低,但在地下水位附近(15~20 m)回升,主要因长期淋溶作用使污染物向下迁移至含水层位富集。浅层地下水Cr6+超标率为73.68%,深层地下水超标率为37.50%。整体而言该区土壤和地下水Cr6+污染形势不容乐观,应开展水土协同治理。空间上表层土壤Cr6+分布受人类活动影响,与投产期厂内功能区的划分相关性较低。地下水Cr6+受水动力场影响,表现为西北高东南低,污染羽中心向下游方向迁移。土壤致癌风险均值为1.85×10−6,介于10−6~10−4,风险中等,应引起必要重视。非致癌风险低于1,无慢性毒害影响。经口摄入土壤为主要暴露途径。地下水致癌风险高达10−2数量级,远超过10−4,风险不可接受。非致癌风险均值51.62,远大于1,极可能引发慢性毒害效应。另外需进一步考虑Cr6+通过牛羊肉食物链进入人体的危害。
Abstract:In order to clarify the Cr6+ pollution status and risk level of a chemical plant, 19 soil and groundwater sampling points were laid in the study area, and the content of Cr6+ was measured, to analyze the pollution characteristics and causes, and to carry out health risk assessment. The results showed that the excess rate of Cr6+ in topsoil (0~0.5 m) was 42.11%. In general, the Cr6+ concentration in soil decreased with the increase of depth, but it rose up near the groundwater level (15~20 m), which was mainly due to the long-term leaching effect that the pollutants migrated down to the aquifer and enriched. The excess rate of Cr6+ in shallow and deep groundwater was 73.68% and 37.50%, respectively. Overall, the Cr6+ pollution is not optimistic, and comprehensive water-soil treatment should be carried out. At the spatial scale the distribution of Cr6+ in topsoil was affected by human activities, and had a low correlation with the division of functional zones. And under the influence of hydrodynamic field, Cr6+ in groundwater is higher in the northwest than the southeast, and the center of the pollution plume has migrated downstream. In addition, the mean value of soil carcinogenic risk is 1.85×10−6, which is between 10−6 and 10−4, and the risk is moderate, which should be paid more attention to. Soil non-carcinogenic risk is less than 1 and doesn’t produce chronic toxic effects. Oral ingestion of soil was the main exposure route. The carcinogenic risk of groundwater is up to the order of 10−2, much more than 10−4, which is unacceptable. The average non-carcinogenic risk is 51.62, much higher than 1, which may cause chronic toxic effects. In addition, the health hazards caused by Cr6+ entering the human body through the beef and mutton food chain should be further studied.
-
Keywords:
- heavy metal /
- chromium /
- pollution characteristics /
- health risk assessment
-
高分辨率实景三维建模技术因其多尺度、多学科属性,忠实映射、高保真度等特性而广泛应用于基础测绘及地质调查等技术领域(胡运海,2012;万荧等,2016;李煜,2017;印森林等,2018;崔玉福等,2021;马建雄等,2022;蔡小超等,2022;李林等,2022;贾俊等,2023)。梁京涛等(2020)开展了基于贴近摄影测量技术的高位崩塌早期识别技术方法研究;闫博等(2020)在地质调查技术领域总结无人机工作流程基础上,对现代沉积、古代构造特征、古代岩体切割关系等方面进行了一些尝试;郑明等(2022)在青藏高原难进入地区开展基于虚拟现实技术的多维度、高精度地质调查实验;刘帅等(2022)基于无人机倾斜摄影技术对露头点实景三维模型建模并绘制地质剖面;常晓艳(2021)在地灾监测技术领域利用无人机倾斜摄影测量技术开展应用研究;刘立等(2022)构建地质灾害隐患三维图,为灾后地质灾害隐患管理、风险分析等提供保障。前人在无人机技术应用领域进行较多研究,但在基于多源数据如智能手机、数字相机、消费级无人机、无人机倾斜摄影等地质剖面建模、剖面测量、剖面投影、剖面制图综合等技术领域涉及较少。笔者在前人研究基础上,基于近景摄影测量和倾斜摄影测量等方法获取高分辨率地质剖面实景三维,结合剖面投影技术和制图综合技术对地质剖面测量和制图进行实验讨论。结果表明,该方法对提高地质剖面作业效率和测量精度有一定现实意义。
1. 实景三维建模技术方法
1.1 技术流程
剖面实景三维建模及处理流程(图1)主要包括获取影像、实景三维建模、实景三维剖面测量、剖面制图综合、成果提交等。
便携式实景三维采集设备包括普通智能手机、GNSS数字相机、消费级无人机等(图2);大面积作业设备多采用多镜头多旋翼或固定翼无人机。依据工作范围、工作目的、作业方式、工作区状况、剖面制图比例等制定实景三维获取方案。如智能手机、GNSS数字相机、消费级无人机满足地质队员野外典型地质现象观察和建模;多镜头无人机等实现大面积长距离地质剖面摄影及建模任务;高分卫星影像建模实现中小比例尺剖面测量和制图。
实景三维建模采用Smart3D等软件进行无人机航片、地面数字相机单独建模或者联合建模。建模坐标信息常采用POS或像控点刺点方式。不同类型设备获取实景三维可利用公共点进行精度互检。实景三维剖面模型可采用通用数据格式实现模型相互转化,常见三维交换格式为如OSGB可实现Smart3D、3Dmax、EPS等软件数据交换。
实景三维剖面测量主要包括剖面要素测量、剖面投影面选择、剖面投影。剖面要素测量主要包括产状测量、接触界线解译、断层绘制等;剖面投影面选择是指依据剖面主要产状、剖面走向等要素确定剖面的投影面;依据剖面投影结果进行立面测量,完成剖面要素采集。
在剖面投影及矢量化采集基础上,依照地图制图综合原理对投影至投影面的剖面各种信息进行选取、简化、夸大、位移等进行综合制图。
基于实景三维的地质剖面成果包括传统纸质图件和新型数字产品,可满足不同比例剖面图的生产,新型数字产品可满足剖面3D打印、数字科普展示、地质三维档案、互联网三维共享等多方面的应用需求,并具有良好的可复用特征。
1.2 实验条件
实景三维地质剖面测量结果与三维影像分辨率关系密切,便携式摄影测量方法获取实景三维建模质量主要取决于近景摄影分辨率、相片重叠度、近景摄影POS定位精度等。便携式近景设备GNSS基线较短,在GNSS定位精度不高情况下易导致剖面模型要素如方位、长度在内的形变。主要通过以下方法预防和改善:①增加足够密度像控点。②适当增加近景摄影GNSS基线网精度。③增加近景摄影重叠度及相片数量。④增加人工定向、测距标志、抽检要素等。
无人机作业应采用航线规划并尽量满足贴近摄影要求,获取垂向和倾斜航片。部分实验区高分辨率卫星三维影像可以用于地质剖面宏观要素的解译,从实验精度和效果看,近景摄影测量技术与卫星遥感技术可以形成技术互补和相互验证。
2. 实景三维剖面测量技术
2.1 产状的测量方法
实景三维是一种精细化建模技术,满足高分辨率地质剖面测量。该技术对岩层产状、接触界线、构造裂隙的观察、解译、统计具有一定优势。依据建模结果可进行产状采集及计算、接触界线解译、裂隙等密图及玫瑰图制作等。三维地质剖面产状测量宜采用三点法直接测量和计算(王俊锋等,2014)。图3为笔者在陕西蓝田县采用外置GNSS数字相机获取的实景三维影像,在三维软件(如EPS)中对构造面理进行实景坐标测量,获得三维坐标并输出,在Struckit软件中完成产状各要素计算,完成倾向玫瑰图、裂隙等密图等绘制。点1-2-3、4-5-6、7-8-9为岩石产状面理,10-11-12为沉积地层与侵入岩接触界线、13-14为岩石中裂隙(图3)。
产状测量应避免遮挡、影像拉花、产状面过小等部位。表1为实景三维及罗盘面理测量结果对比,倾向之差≤12°,倾角之差≤10°。依据卢立吉等(2016)方法检查本次实景三维及罗盘产状测量结果,未发现存在显著精度差异。
表 1 实景三维与罗盘产状测量对比Table 1. Comparison of real 3D and compass attitude measurement序号 建模方法 实景三维 罗盘测量 差值 倾向 倾角 倾向 倾角 倾向 倾角 1 单反相机
近景摄影
建模113° 44° 118° 41° −5° 3° 2 134° 74° 133° 77° 1° −3° 3 121° 45° 124° 50° −3° −5° 4 106° 57° 113° 51° −7° 6° 5 109° 59° 102° 61° 7° −2° 6 125° 41° 114° 47° 11° −6° 1 无人机摄影
建模(图5f)189° 88° 187° 84° 2° 4° 2 182° 87° 181° 79° −1° 8° 3 358° 84° 1° 81° −3° 3° 4 189° 87° 192° 86° −3° 1° 5 275° 54° 279° 54° −4° 0° 2.2 实景三维剖面投影原理
空间点在平面的投影可以看成过经特定点的直线与平面的交点。空间点在平面上的投影问题可概括为直线和平面交点求解(朱鹏先等,2021):
设非零向量为:
$$ g = \left( {m,n,p} \right) $$ (1) 设实景三维体一点为:
$$ w = \left( {{x_0},{y_0},{z_0}} \right) $$ (2) 设投影平面S方程为:
$$ Ax + By + Cz + D = 0 $$ (3) 设水平面为H。
经过w点的向量g方程为:
$$ \frac{{x - {x_0}}}{m} = \frac{{y - {y_0}}}{n} = \frac{{z - {z_0}}}{p} $$ (4) 当Am+Bn+Cp≠0时,非零向量与投影面有交点,交点为投影点坐标:
$$ x = {x_0} - \frac{{m(A{x_0} + B{y_0} + C{z_0} + D)}}{{Am + Bn + Cp}} $$ (5) $$ y = {y_0} - \frac{{n(A{x_0} + B{y_0} + C{z_0} + D)}}{{Am + Bn + Cp}} $$ (6) $$ z = {z_0} - \frac{{p(A{x_0} + B{y_0} + C{z_0} + D)}}{{Am + Bn + Cp}} $$ (7) 实景三维与投影面位置关系包括相割、相离、相切等关系(图4)。投影面可依据地质剖面的产状、剖面走向、接触界线等关系设置一个或多个投影面。多投影面可采用分段方式分段投影,不同投影面之间坐标系统存差异。
实景三维剖面投影主要存在几何体的长度、面积、角度变形,如剖面中岩层厚度、共轭节理夹角、表面积等几何要素都可能存在变形。剖面中长度变形对岩层厚度影响较大,设l 为岩层真厚度, l’ 为岩层投影后厚度,g⊥S、S⊥H,θ为向量g与岩层面走向夹角,g为平行投影时满足:l’ =l ×cosθ,当θ=0°时长度变形最小,当θ=30°时,长度变形约大于13%。剖面岩层总体走向、断层走向、接触接线走向等与g夹角不建议大于30°。
当S⊥H时实景三维投影垂向变形最小;当g为平行向量且满足g⊥S且S⊥H时为正投影,正投影时实景三维剖面在投影面S中横轴整体变形最小。因此,实景三维剖面投影面宜选择正投影。剖面走向方位的选择与岩层走向、接触界线走向、断层走向相协调,避免剖面的长度、角度、面积变形过大。
2.3 剖面测量方法
实景三维常见测量方法主要为直接测量法和立面图绘制法。直接测量法从实景三维影像中直接测量坐标位置、岩层厚度等利用制图软件直接进行剖面绘制,该方法适合实景三维直接绘制柱状剖面。如黄土高原中生界水平地层(图5b、图5e、图5g)、探槽、土壤样坑、拨土等工程。
立面图是一种平行正投影的结果,投影面与水平面垂直,常用于建筑物外立面绘制(王莫,2015)。常用软件如EPS等提供立面图绘制功能。利用立面图功能,采集三维剖面要素。主要步骤:①依据剖面测量需要,裁剪实景三维。②依据剖面位置、岩层走向、剖面走向等绘制剖面线并构造垂向投影面。③立面模式下进行各要素选取、采集等,依照采集结果进行地质剖面制图综合。
2.4 测量精度检验方法
在现有技术条件下,实景三维剖面精度检验常见方法有像控点(特征点)检验法、标志抽样检验法、激光水平标检验法等。
像控点(特征点)检验法:在获取实景三维剖面附近按照一定几何间距布设一定数量或一定密度的像控点,利用传统测量仪器如全站仪或CORS等进行测量并获取空间三维坐标,通过空间三维坐标与像控点或特征点的坐标残差分析,以评价实景三维采集精度的方法。该方法在大面积、长距离剖面等较常用。
标志抽样检验法指依照钢尺、罗盘等传统仪器,对实景三维剖面采集的同位置岩层产状、岩层厚度、裂隙走向、构造面理等进行测量检核,以评价实景三维建模过程层中的整体变形及局部细节建模精度。
激光水平标检验法:在狭小的空间如土壤探坑、地质探槽等地区,利用激光水平装置制造激光水平线或激光垂向交叉线,激光水平线满足高程相等、满足激光线走向一致等特性,在获取实景三维后通过检验实景三维影像中激光水平线的高程,或者激光标志线的走向特性。
实景三维影像获取时可利用测量标尺、直角三角形等特定的几何模型构造长度、角度、面积等检验条件进行精度检验。
3. 地质剖面制图综合
与传统遥感解译方式不同,高分辨率实景三维地质剖面分辨率可达毫米级,其地质单元色调和纹理更接近目视效果,依据剖面地质单元几何特征、色调、纹理等通过三维模型直接测量可完成地质剖面制图。剖面制图比例尺理论精度为0.1 mm,表2为实景三维地质剖面制图比例与实景三维分辨率对应表。
表 2 制图比例与分辨率Table 2. Mapping scale and resolution制图比例 分辨率(m) 制图比例 分辨率(m) 1∶50 0.005 1∶500 0.05 1∶100 0.01 1∶ 1000 0.1 1∶250 0.025 1∶ 2000 0.2 地质剖面制图综合是指将大比例尺剖面图缩编成中小比例尺图的过程,依据用图性质、比例尺等通过简化、概括、抽象等方法舍去次要信息,使得图面层次协调、对象清晰的制图过程(曲平等,2012)。常见地质剖面制图综合方法有逐级缩编法、自动综合法等。低分辨率实景三维不能绘制高于同等比例尺精度地质剖面,高分辨率实景三维影像可以通过制图综合方法绘制中小比例尺剖面图。
3.1 逐级缩编法
实景三维地质剖面测量应进行最大比例尺测图,依据制图综合原理将大比例尺向中小比例尺逐级缩编。通过制图综合过程中的简化、位移、夸大、合并、分割等方法完成剖面制图综合。
例如,获取0.01 m分辨率三维影像,应采集1∶100数字剖面,依照制图综合方法逐级完成1∶250、1∶500、1∶
1000 、1∶2000 剖面绘制。剖面绘制过程中,对重要部位进行夸大,非重要部位简化、合并、位移等概括表示。岩层产状、断层走向等要素应在最大分辨率下完成采集,制图综合过程中对密集或次要信息进行抽稀或概略。3.2 自动综合法
自动综合法是指通过技术手段自动生成任意中小比例尺数据,信息量随着比例尺变化自动增减,信息的压缩与复原和比例尺是自动适应的(赵丽宁等,2001)。随地质剖面显示比例尺变化,断层、产状、岩性花纹、裂隙等剖面信息自适应增减,并呈现信息承载的最优化特征。常见的自动综合算法有优先级法(刘波等,2008)、结点自动删除法(李靖涵,2018)、弯曲简化法(段修梅,2012)、移动平均法等方法对大比例地质剖面进行自动选取、简化等。
4. 典型案例分析
4.1 黄土高原中生界垂直剖面
试验区位于陕西省铜川市,出露岩石地层有中三叠统纸坊组(T2z)、上三叠统铜川组(T3t)、胡家组(T3h)、永坪组(T3y)、瓦窑堡组(T3w)、下侏罗统延安组(J1y)、直罗组(J1z)、下白垩统宜君组(K1y)、洛河组(K1l)等。地层产状发育较缓,地质剖面常见为沿沟系发育的垂直剖面,主要岩性为长石石英砂岩、砾岩、泥岩、油页岩、碳质页岩、煤线等。
本区域实景三维资料获取主要采用大疆精灵4无人机平台,采用垂直及倾斜航线飞行,航高约10~30 m,航向重叠度设置80%,旁向重叠度设置60%。便携式近景GNSS数字相机路线摄影测量补充。完成三叠系、侏罗系、白垩系地质剖面总计约2.6 km。实验对比实景三维影(图5b、图5e、图5g)及地质剖面特征。
(1)通过GNSS数字相机获取实景三维影像分辨率可达1~5 mm,基本可以识别10~30 mm砾石等要素,如宜君组中大量砾岩。砂岩、泥岩、煤线等要素几何表达清晰可见,岩层产状、岩层厚度等要素测量结果与检查结果吻合较好。实景三维建模应用效果优于同位置肉眼观察,部分人员难以涉足陡崖地形区有显著的测量优势。
(2)无人机低空摄影可满足大面积长距离剖面的实景三维建模。分辨率可达0.5~5 cm,无人机建模分辨率略低于数字相机,但砂岩、泥岩、煤线等要素较易区分,砂岩中砾石、岩石中垂向节理及小型断裂亦可清晰分辨。实验表明,无人机剖面应用从效率及测量精度上与传统罗盘、钢尺方法无显著区别,与RTK方法(李玲等,2018)地质剖面相比具有显著的效率优势。
(3)GNSS数字相机近景摄影较无人机摄影具有一定高分辨率优势,其缺点也较明显如:摄影基线网较短且在摄影基线布置不合理情况下易导致三维模型垂向变形较大。无人机建模后,通过与数字影像公共点采集一定密度像控点、无人机影像和数字影像联合建模等方案可减小数字相机或者智能手机影像建模的局部误差。
4.2 土壤样坑剖面建模
试验区位于陕西关中平原区农田,为研究农田表层土壤及植物根系状况,进行土壤取样采样坑实验,实验坑长宽高约1 m×1 m×1 m。采用大疆御无人机及GNSS数字相机进行对比实验(图5d)。无人机航高5~15 m,航向重叠度80%,旁向重叠度60%;GNSS数字相机地表拍摄,相机与目标距离约2~5 m,从不同方向完成采样坑相片150余张。
通过实验对比,GNSS数字相机与无人机获取探坑几何要素精度基本相当,直线定向、长度变形差异可忽略不计,满足传统采样坑的三维建模精度要求。
GNSS数字相机实景三维分辨率优于1 mm,实景三维中可见清晰的小麦根须以及杂草根茎。无人机获取实景三维优于5 mm,采样坑土壤含水、土壤成分变化等色调信息实景三维中能较好体现。实景三维土壤采样坑建模结果优于同等距离条件下肉眼观察要素。
4.3 陕西宁陕白钨矿探槽
地质探槽、拨土等是地质施工中重要的工作方法和实物依据,三维可视化后具有较好的工程应用价值。实景三维探槽能够较好的复原探槽及矿体或矿化体的赋存位置及其几何形态。可采用普通智能手机或者GNSS数字相机完成建模数据获取。普通智能手机应布设相控或获得GNSS三维坐标。
试验区位于陕西省宁陕县花岗岩与地层接触部位,因本区域高大植被覆盖较好,无人机完成航摄较为困难。从东、南、西、北、东北、西北、西南、东南等不同方向获取建模相片222张,并从垂向获取探槽壁、探槽底部不同方位的相片,相片重叠度优于90%。数字相机或智能手机与探槽距离约2~10 m。
综合实验结果显示:①智能手机相片建模分辨率优于1~5 mm,GNSS数字相机分辨率优于1~5 mm。实景三维模型分辨率优于同等条件下肉眼所见。②探槽测量中采用激光水平光源制造激光水平线,三维建模后利用激光水平线验证探槽垂向变形状况(图5c)。实景三维产状测量结果与传统罗盘无显著差异,探槽长度等几何形变可忽略不计。③采用激光水平装置制造垂直或水平相交激光线及交点,或采用激光反光标志等制造相控点,并利用全站仪或CORS测量。标志点作为建模像控点或检验点,GNSS定位困难时,将提高三维建模的几何精度。④白钨矿在紫外光下呈现出浅蓝色萤光,紫外光源激发状态下获取近景相片,重建实景三维模型,有利于立体环境下矿体或矿化体观察和三维空间角度矿体的赋存状态研究。
4.4 中高山区中小比例剖面
航空摄影测量存在困难时可采用高分辨率卫星遥感影像获取三维模型,或采用DEM与正射影像形成三维模型,也可采用三维模型与近景摄影测量等综合方式形成三维模型,完成相应的地质剖面测量、制图综合等。
试验区位于新疆乌苏博罗克鲁山腹地,出露地层为东图津河组(C2dt)。剖面沿南北向断裂发育,断裂东西两侧崖壁60°~90°近直立状,崖壁高差约400 m。通过高分卫片、DEM、近景摄影测量等综合方法建模获取实景三维(图5)。近景相片分辨率优于20 cm,实景三维中一般的地层界线、构造界线、垂向节理等可清晰解译。东图津河组一段(C2dt1)主要为火山碎屑岩,近景相片中呈深灰色或灰黑色,二段(C2dt2)主要岩性灰岩,近景相片呈乳白色。东图津河组地质剖面总长度约2.8 km,剖面北侧灰岩段逆冲于火山碎屑岩之上;南测灰岩段与火山碎屑岩段之间发育大量的互层或条带。灰岩中穿插大量火山碎屑岩条带,火山碎屑岩条带随高度变化呈现厚度变薄、频次减弱特征(图6)。依实景三维剖面及测量结果推测,东图津河组海陆交互相特征明显(张天继等,2006)。
实验表明,中高山区或者造山带地质剖面实景三维测量还应参考如下几个步骤和原则:①剖面宏观走向依据岩层的主要产状、接触界线等整体投影形变最小的原则。②实景三维测量先整体后局部。先整体指先绘制剖面的整体界线、地层界线、断层界线等;后局部是指在宏观界线基础上,绘制岩层、裂隙、测量产状等。③剖面的绘制应该按照大比例尺、中比例尺、小比例尺逐级缩编,剖面逐级缩编需要按照制图综合的基本原则完成,对具有特殊意义的地层或非正式填图单元按照一定制图综合原则合理放大、平移、合并等。
5. 结论
(1)测图精度高、效率高。普通数字相机或智能相机像素目前已达到1亿像素,对应地质剖面分辨率部分优于毫米级。通过测试,部分实景三维剖面整体建模效果优于同位置肉眼观察。
(2)设备多样,采集方案多样。地质剖面实景三维获取方式多样,长距离、大面积地质剖面可采用多镜头长航时无人机低空采集,地质队员单人携带便携式摄影器材如数字相机、智能手机等可以完成剖面地质调查或路线地质调查典型地质现象的实景三维的获取工作。
(3)数据共享。实景三维数据基于通用三维建模格式,模型具有良好的几何特征和数据共享特征。实景三维模型存储格式具有精确的几何特征和丰富的纹理特征,可以作为重要工程施工的佐证档案。实景三维与传统的存档相片相比具有精确的可测量性、信息多样性等特征。
(4)基于多光谱相机的实景三维地质剖面,将为剖面测量和矿物诊断提供良好的制图依据。
致谢:本次实验过程感谢中国地质调查局西安地质调查中心朱涛高级工程师、西安地质矿产勘查开发院苏英明等工程师参与研究并提出宝贵意见。
-
表 1 健康风险评价模型相关参数表
Table 1 Parameters of health risk assessment model
类型 参数 单位 取值 数据来源 暴露量评估 ABSd / 0.001 王珊等,2019 Rs mg/d 100 《建设用地土壤污染风险
评估技术导则》(HJ 25.3-2019)ED a 25 EF d/a 250 BW kg 61.8 ABSo / 1 AT d 致癌:27740;
非致癌:9125SSAR mg/cm2 0.2 Ev 次/d 1 H cm 161.5 SER / 0.18 PM10 mg/m3 0.119 Ra m3/d 14.5 PIAF / 0.75 fspi / 0.8 fspo / 0.5 EFi d/a 187.5 Rw L/d 1 Efo d/a 62.5 SAE cm2 765.97 公式(3)计算 毒性评估 IUR m3/mg 0.12 《建设用地土壤污染风险
评估技术导则》(HJ 25.3-2019)RfC mg/m3 0.0001 ABSgi / 0.025 风险评估 SFo [mg/(kg·d)]−1 0.5 王珊等,2019 SFd [mg/(kg·d)]−1 20 公式(10)和(12)计算 SFp [mg/(kg·d)]−1 0.51 RfDo mg/(kg·d) 0.003 《建设用地土壤污染风险
评估技术导则》(HJ 25.3-2019)RfDd mg/(kg·d) 7.50E-05 公式(11)和(13)计算 RfDp mg/(kg·d) 2.35E-05 SAF / 0.5 《建设用地土壤污染风险
评估技术导则》(HJ 25.3-2019)WAF / 0.5 表 2 土壤Cr6+污染浓度统计特征
Table 2 Statistical characteristics of Cr6+ pollution concentration in soil
深度(m) 采样数(个) 最大值(mg/kg) 最小值(mg/kg) 平均值(mg/kg) 方差 C·V 0~0.5 19 37.60 0.00 9.56 11.80 1.24 0.5~2.0 19 27.50 0.00 2.39 6.19 2.58 2.0~4.0 19 5.50 0.00 0.81 1.60 1.97 4.0~6.0 19 1.20 0.00 0.17 0.40 2.32 6.0~8.0 19 4.00 0.00 0.62 1.42 2.31 8.0~10.0 19 4.20 0.00 0.39 1.04 2.62 10.0~15.0 19 16.10 0.00 2.00 3.83 1.92 15.0~20.0 19 20.70 0.00 3.95 5.86 1.48 20.0~25.0 19 12.50 0.00 2.01 3.60 1.80 25.0~30.0 16 2.00 0.00 0.18 0.52 2.85 30.0~孔底 19 1.60 0.00 0.08 0.36 4.24 注:超标率计算以《土壤环境质量 建设用地土壤污染风险管控标准》(GB36600-2018)中第二类用地Cr6+风险筛选值(5.7 mg/kg)为依据。 表 3 地下水Cr6+污染浓度统计特征(mg/kg)
Table 3 Statistical characteristics of Cr6+ pollution concentration in groundwater (mg/kg)
取样位置 采样数 Max Min Ave SD C·V 浅层 19 49.90 0.00 6.58 11.29 1.71 深层 8 1.48 0.00 0.25 0.48 1.95 注:超标率计算以《地下水质量标准》(GB/T 14848-2017)III类标准Cr6+限值(0.05 mg/L)为依据。 表 4 土壤及地下水重金属Cr6+的健康风险表
Table 4 Health risks of heavy metal Cr6+ in soil and groundwater
环境介质 暴露途径 致癌风险 非致癌风险 风险范围与均值 贡献率 风险范围与均值 贡献率 土壤 经口摄入 0~6.854×10−6 94.01% 0~0.028 72.80% 1.742×10−6 0.0071 皮肤接触 0~4.200×10−7 5.77% 0~0.0017 4.46% 1.070×10−7 0.0004 呼吸吸入 0~1.668×10−8 0.23% 0~0.0085 22.25% 4.255×10−9 0.0022 3种途径
健康总风险0~7.291×10−6 100.00% 0~0.0380 100.00% 1.853×10−6 0.0097 地下水 饮用地下水 0~9.096×10−2 100.00% 0~368.70 100.00% 1.274×10−2 51.62 -
白福高. 西北某铬污染场地地下水Cr(VI)污染的抽出处理数值模拟研究[D]. 北京: 中国地质大学(北京), 2017 BAI Fugao. Study on Pumping Strategy for Cr(VI) Contaminated Groundwater of a Northwest of China Chromium Contaminated site on Numerical Simulation[D]. Beijing: China University of Geosciences, 2017.
方晴, 冼萍, 蒙政成. 基于蒙特卡罗模拟的农用地土壤健康风险评价[J]. 环境工程, 2021, 39(2): 147-152 FANG Qing, XIAN Ping, MENG Zhengcheng. Environmental health risk assessment model of agricultural land based on Monte Carlo simulation and it’s Application[J]. Environmental Engineering, 2021, 39(2): 147-152.
郭媛媛. 铬在地下含水层中的迁移转化特征[D]. 长春: 吉林大学, 2008 GUO Yuanyuan. Characteristics of movement and transformation of chromium in underground aquifer[D]. Changchun: Jilin University, 2008.
顾小凡, 党学亚, 杨炳超, 等. 延安吴起县地下水中Cr6+分布规律及来源探讨[J]. 西北地质, 2015, 48(4): 190-203 doi: 10.3969/j.issn.1009-6248.2015.04.019 GU Xiaofan, DANG Xueya, YANG Bingchao, et al. The Distribution Regularity and Sources of Six-valence Chromium in Groundwater from Wuqi County, Yanan City[J]. Northwestern Geology, 2015, 48(4): 190-203. doi: 10.3969/j.issn.1009-6248.2015.04.019
高文武, 姜燕, 赵晋陵. 基于协同克里金插值法的土壤锰元素含量预测[J]. 地理与地理信息科学, 2018, 34(3): 119-124 doi: 10.3969/j.issn.1672-0504.2018.03.019 GAO Wenwu, JIANG Yan, ZHAO Jinling. Predicting and Mapping the Mn Content in Soil Based on Cokriging[J]. Geography and Geo-Information Science, 2018, 34(3): 119-124. doi: 10.3969/j.issn.1672-0504.2018.03.019
高瑞忠, 秦子元, 张生, 等. 吉兰泰盐湖盆地地下水Cr6+、As、Hg健康风险评价[J]. 中国环境科学, 2018, 38(6): 2353-2362 doi: 10.3969/j.issn.1000-6923.2018.06.040 GAO Ruizhong, QIN Ziyuan, ZHANG Sheng, et al. Health risk assessment of Cr6+, As and Hg in groundwater of Jilantai salt lake basin, China[J]. China Environmental Science, 2018, 38(6): 2353-2362. doi: 10.3969/j.issn.1000-6923.2018.06.040
高雅, 胡晨, 张春雷, 等. 安徽石台地区富硒土壤分布及硒的富集迁移规律探讨[J]. 西北地质, 2022, 55(2): 284-291 GAO Ya, HU Chen, ZHANG Chunlei, et al. Study on the distribution of selenium-rich soil and the regularity of selenium enrichment-migration in Shitai area, Anhui, China[J]. Northwestern Geology, 2022, 55(2): 284-291.
韩琳, 徐夕博. 基于PMF模型及地统计的土壤重金属健康风险定量评价[J]. 环境科学, 2020, 41(11): 5114-5124 HAN Lin, XU Xibo. Quantitative Evaluation of Human Health Risk of Heavy Metals in Soils Based on Positive Matrix Factorization Model and Geo-statistics[J]. 2020, 41(11): 5114-5124.
蒋兴超, 许静, 李如意, 等. 广东省汕头市土壤铬的空间分布特征、来源解析及影响因素研究[J]. 地学前缘, 2023, 30(2): 514-525 JIANG Xingchao, XU Jing, LI Ruyi, et al. Soil chromium in Shantou City, Guangdong Province: Spatial distribution characteristics, source apportionment and influencing factors[J]. Earth Science Frontiers, 2023, 30(2): 514-525.
李晶晶, 彭恩泽. 综述铬在土壤和植物中的赋存形式及迁移规律[J]. 工业安全与环保, 2005, 31(3): 31-33 LI Jingjing, PENG Enze. Summarization on the existing from and transferring rules of chroming in soil[J]. Industrial Safety and Environmental Protection, 2005, 31(3): 31-33.
刘瑞平, 朱桦, 亢明仲, 等. 大荔县地下水环境质量评价及成因浅析[J]. 西北地质, 2009, 42(2): 116-125 LIU Ruiping, ZHU Hua, KANG Mingzhong, et al. Assessment of water environment quality and pollution factors for DAli county[J]. Northwestern Geology, 2009, 42(2): 116-125.
廉晶晶, 罗泽娇, 靳孟贵. 某厂电镀车间场地土壤与地下水污染特征[J]. 地质科技情报, 2013, 32(2): 150-155 LIAN Jingjing, LUO Zejiao, JIN Menggui. Contamination Characteristics of Soil and Groundwater in Electroplating Plant[J]. Bulletin of Geological Science and Technology, 2013, 32(2): 150-155.
刘伟江, 陈坚, 刘锐, 等. 郯城某化工厂周边地下水污染现状调查与评价[J]. 安全与环境工程, 2018, 25(6): 67-75 LIU Weijiang, CHEN Jian, LIU Rui, et al. Investigation and evaluation of groundwater pollution around a chemical plant in Tancheng country[J]. Safety and Environmental Engineering, 2018, 25(6): 67-75.
刘柱光, 方樟, 丁小凡. 燃煤电厂贮灰场土壤重金属污染及健康风险评价[J]. 生态环境学报, 2021, 30(9): 1916-1922 LIU Zhuguang, FANG Zhang, DING Xiaofan. Heavy metal pollution and health risk assessment of soil in ash yard of coal-fired power plant [J]. Ecology and Environmental Sciences, 2021, 30(9): 1916-1922.
马海珍. 白洋淀流域平原区地下水环境健康风险评价及预测[D]. 西安: 长安大学, 2021 MA Haizhen. Health risk assessment and prediction of groundwater in plain area of Baiyangdian Basin[D]. Xi’an: Chang’an University, 2021.
彭叶棉, 杨阳, 侯素霞, 等. 外源六价铬在土壤中的有效性及其小麦毒性效应[J]. 生态环境学报, 2020, 29(2): 369-377 PENG Yemian, YANG Yang, HOU Suxia, et al. The bio-availability of exogenous Cr(VI) in soils and its toxic effect on wheat [J]. Ecology and Environmental Sciences, 2020, 29(2): 369-377.
青海省海北州海晏县地下水污染防治试点项目实施方案[R]. 海北藏族自治州生态环境局, 2021. 史锐, 岳荣, 张红. 有色金属采选冶基地周边土壤中重金属纵向分层研究[J]. 土壤通报, 2016, 47(1): 186-191 SHI Rui, YUE Rong, ZHANG Hong. Research on Vertical Distribution of Heavy Metal in Soil around Non-ferrous Metal Industry Area[J]. Chinese Journal of Soil Science, 2016, 47(1): 186-191.
吴敦敖, 鲁文毓. 铬在土壤-地下水系统中的污染研究[J]. 环境科学学报, 1991, 11(3): 276-283 doi: 10.13671/j.hjkxxb.1991.03.004 WU Dun’ao, LU Wenyu. Study on chromium contamination in soil-groundwater system[J]. Acta Scientiae Circumstantiae, 1991, 11(3): 276-283. doi: 10.13671/j.hjkxxb.1991.03.004
王珊, 魏海春. 2018年我国中东部局部地区农田土壤典型重金属健康风险评估[J]. 环境与健康杂志, 2019, 36(9): 807-810 WANG Shan, WEI Haichun. Health risk assessment of typical heavy metals in farmland soils in parts of central and eastern China in 2018[J]. Journal of Environment and Health, 2019, 36(9): 807-810.
王露艳, 刘干斌, 周晔, 等. 电镀场地重金属铬污染土固化率及稳定性研究[J]. 水文地质工程地质, 2022, 49(4): 183-189 WANG Luyan, LIU Ganbin, ZHOU Ye, et al. A study of the curing rate and stability of heavy metal chromium contaminated soil at electroplating sites[J]. Hydrogeology & Engineering Geology, 2022, 49(4): 183-189.
王蕊, 陈楠, 张二喜. 基于总量与形态的矿区周边土壤重金属生态风险与健康风险评估[J]. 环境科学, 2022, 43(3): 1546-1557 WANG Rui, CHEN Nan, ZHANG Erxi. Ecological and Health Risks Assessment Based on the Total Amount and Speciation of Heavy Metals in Soils Around Mining Areas[J]. Environmental Science, 2022, 43(3): 1546-1557.
徐腾, 南丰, 蒋晓锋, 等. 制革场地土壤和地下水中铬污染来源及污染特征研究进展[J]. 土壤学报, 2020, 57(6): 1341-1352 XU Teng, NAN Feng, JIANG Xiaofeng, et al. Advances on Sources and Characteristics of Chromium Pollution in Soils and Groundwater of Tannery Sites[J]. Acta Pedologica Sinica, 2020, 57(6): 1341-1352.
余飞, 张永文, 严明书, 等. 重庆汞矿区耕地土壤和农作物重金属污染状况及健康风险评价[J]. 环境化学, 2022, 41(2): 1-13 YU Fei, ZHANG Yongwen, YAN Mmingshu, et al. Heavy metal pollution and human health risks assessment of soil and crops near the mercuryore in Chongqing[J]. Environmental Chemistry, 2022, 41 (2): 1-13.
周文武, 陈冠益, 穷达卓玛, 等. 拉萨市垃圾填埋场地下水水质的居民健康风险评价[J]. 环境化学, 2020, 39(6): 1513-1522 doi: 10.7524/j.issn.0254-6108.2019041101 ZHOU Wenwu, CHEN Guanyi, QIONG Dazhuoma, et al. Health risk assessment of groundwater quality in Lhasa landfill[J]. Environmental Chemistry, 2020, 39(6): 1513-1522. doi: 10.7524/j.issn.0254-6108.2019041101
Costa Max, Klein Catherine B. Toxicity and Carcinogenicity of Chromium Compounds in Humans[J]. Critical Reviews in Toxicology, 2006, 36(2): 155-163. doi: 10.1080/10408440500534032
Dilek G. Turer, Barry J. Maynard. Heavy metal contamination in highway soils: comparison of Corpus Christi, Texas and Cincinnati, Ohio shows organic matter is key to mobility[J]. Clean technologies and Environmental policy, 2002, 4(4): 235-245.
Eziz Mamattursun, MOHAMMAD Anwar, MAMUT Ajigul, et al. A human health risk assessment of heavy metals in agricultural soils of Yanqi Basin, Silk Road Economic Belt, China [J]. Human and Ecological Risk Assessment, 2018, 24(5-6): 1352-1366.
Kong Jing, Guo Qingjun, Wei Rongfei, et al. Contamination of heavy metal sand isotopic tracing of Pb in surface and profile soils in a polluted farmland from a typical karst area in sounthern China[J]. Science of the Total Environment, 2018, 637/638: 1035-1045.
Zeng Yanyan, Zhou Jinlong, Zhou Yinzhu, et al. Assessment and causes of groundwater organic pollution in typical plain areas in Xinjiang, China[J]. Exposure and Health, 2016, 8(3): 401-417. doi: 10.1007/s12403-016-0211-0