负压条件下地层中VOCs污染的数值模拟分析.pdf

6环境科学导刊2023，42（4）CN 53-1205/XISSN 1673-9655负压条件下地层中 VOCs 污染的数值模拟分析程 昕1,3，赵胜豪2，钱 瑜1,3（1.南京大学污染控制与资源化研究国家重点实验室，江苏 南京 210023；2.南京国环科技股份有限公司，江苏 南京 210042；3.南京大学环境学院，江苏 南京 210023）摘 要：以ASU House为场地案例进行三维数值模拟，使用达西定律与多孔介质中稀物质传递定律，模拟土壤气速度场与进入建筑物的三氯乙烯通量，从而获得负压条件下地层中气态三氯乙烯的时空分布规律。结果显示，与基线状态相比，负压条件下土壤气中的三氯乙烯浓度上升，蒸气入侵效应显著增强；距负压中心与地下污染源的距离越近，三氯乙烯浓度越高。关键词：VOCs；蒸气入侵；负压控制；数值模拟中图分类号：X14 文献标志码：A 文章编号：1673-9655（2023）04-0006-050 引言蒸气入侵指VOCs从地下污染源通过挥发释放出来后，以气态形式经由土壤包气带孔隙运移到建筑物地基附近，最终进入建筑物室内，并产生人体暴露的过程1。在原址改造类的棕地再开发项目中，原场地部分或全部建筑物将被保留2，蒸气入侵成为VOCs暴露的重要途径，其风险评估结果尤为重要3-4。影响蒸气入侵风险评估准确性的最大因素是VOCs浓度具有时间和空间异质性，仅依据单次的地下水、土壤、土壤气等采样结果，难以准确评估5。有学者提出，通过人工调控场地建筑物为负压状态后再开展污染物采样工作，可以消除时空异质性的干扰6。在实际场地调研方面，McHugh等7在美国六栋建筑中进行了压力控制方法的验证，在基线、负压和正压控制条件下分别进行VOCs采样。研究结果表明通过在窗户或门口安装风扇可以有效调控建筑物的气压状态，其中负压控制可以解决污染物的时空异质性以及背景干扰问题，该方法可以代替需要多次进行的土壤气采样监测。Holton等8针对ASU House开展长期实验，在超过一年的时间内执行负压控制，将各监测指标与自然状态进行对比分析。研究结果表明负压条件下污染物浓度的时空异质性降低，且监测结果比自然状态高12个数量级，因此该方法可以用于最坏情形下蒸气入侵的评估。Guo等9-10在蒸气入侵污染场地多次执行负压控制实验，研究排风扇位置、室内外压差、压差控制持续时间、室内空气混合、室内空气取样位置与规范等因素对污染物采样结果的影响。在模型模拟研究方面，Song等11将现有的空气入渗建筑物模型与蒸气入侵模型相结合，构建烟囱效应（负压）、风效应与土壤气进入速率之间的定量关系。Yao等12通过数值模型，分析建筑物的室内压力波动对蒸气入侵的影响。研究表明当建筑物室内压力波动超过5 Pa/h时，室内污染物浓度将显著提高。目前针对负压控制下蒸气入侵风险评估的研究，大多为场地实验执行结果的分析与讨论，而对负压条件下地层中VOCs时空分布规律的机理研究尚存不足。因此，本研究以污染场地ASU House为案例，使用达西定律与多孔介质中稀物质传递定律，构建负压条件下三氯乙烯蒸气在地层中迁移的三维数值模型，模拟负压条件下气态三氯乙烯的时空分布，分析其时空分布规律，指导蒸气入侵污染场地的调查与风险评估。1 负压条件下地层中 VOCs 迁移的模型1.1 模拟场景本研究所模拟的场景为典型蒸气入侵场景下，负压控制时VOCs蒸气在地层中的迁移过程。蒸气入侵概念模型如图1所示。负压控制通过利用建筑物现有的暖通空调系统或单独安装风扇风机来实现，其原理在于建筑物室内和土壤气之间的气压差决定了穿越地板的气流流向13。负压控制执行时，建筑物内部形成负压中心，在压力的驱收稿日期：2022-11-25基金项目：国家重点研发计划项目（No.2016YFC0207603）。作者简介：程昕（1998-），女，硕士研究生，研究方向：土壤污染状况调查与风险评估。负压条件下地层中 VOCs 污染的数值模拟分析程昕7动下，地层中的VOCs蒸气向建筑物地基处迁移，并通过地基处的裂缝进入室内14。而本研究主要关注负压条件下，VOCs蒸气在地层中的迁移过程。模拟场景主要由以受污染地下水为主的挥发性气体释放源、土壤介质以及目标建筑物地基构成，且污染源与建筑物室内空气、地表大气之间存在明显的浓度差。图1 概念模型图1.2 模拟方法及基础原理本研究使用的模拟软件为Comsol Multiphysics。构建模型的主要过程有：构建几何体、设置物理场、划分有限元网格、求解计算和结果后处理。模型涉及的物理场为达西定律和多孔介质中稀物质传递两个模块。其中，达西定律模块用于模拟流体在多孔介质中通过间隙的流动，主要适用于低速流动介质或渗透率和孔隙率较小的介质。对于这些介质，压力梯度是流体流动的主要驱动力，流动主要受孔隙内摩擦阻力的影响。而多孔介质中稀物质传递模块用于计算多孔介质中的物质浓度和输运，可以模拟物质在多孔介质中的扩散、对流、迁移、分散、吸附、挥发等过程，适用于一个或多个稀物质，在饱和或不饱和多孔介质中的移动。在达西定律模块中，根据压力分布求解土壤气体流场：（1）式中：u土壤气的流速，m/s；k土壤的空气渗透系数，m2；土壤气的动态黏滞系数，g/(ms)；土壤气的气压梯度，g/(m2s2)。在多孔介质中稀物质传递模块中，土壤气体浓度的时间和空间分布通过将土壤气体流速u代入以下公式求解：()gigeffgcaucDct=+（2）式中：ai总孔隙率，m3-soil gas/m3-soil；cg污染物在土壤气中的浓度，g/m3；t时间，s；u土壤气的流速，m/s；Deff 土壤中污染物的有效扩散系数15。iocwiocbgiiKfaHH=+（3）式中：ai总孔隙率，m3-soil gas/m3-soil；土壤中的气体孔隙度，m3-soil gas/m3-soil；土壤中的水分孔隙度，m3-H2O/m3-soil；Hi污染物i的亨利常数，无量纲；土壤中有机碳对污染物i的吸附系数，L/g；foc土壤中的有机碳质量分数，无量纲；b土壤密度，g/L。10/310/322ggweffgtitDDDH=+（4）式中：Deff 土壤中污染物的有效扩散系数，m2/s；Dg气相中污染物的扩散系数，m2/s；Dw水相中污染物的扩散系数，m2/s；土壤中的气体孔隙度，m3-soil gas/m3-soil；土壤中的水分孔隙度，m3-H2O/m3-soil；土壤的气体孔隙度和水分孔隙度的总和。1.3 初始及边界条件设置模型设置了多层土壤结构。假设污染源为无限大的地下水源，假设建筑物地基处的裂缝是污染物从土壤进入室内的唯一途径，不考虑污染物蒸气入侵过程中的生物降解作用，渗流区除地表外的边界表面的边界条件均考虑为无流量。通常地表都是直接与大气接触，因此该表面上的压强等于大气压。裂缝设置：根据场地的实地调研，获取地基处裂隙的位置、长度、宽度、深度等数值，以压力井的边界条件模拟裂隙，设置相应的井宽与压力参数。与过往学者研究中假设地基四周处存在裂缝相比，提高了模拟的准确性。污染源通量设置：基于负压是VOCs从地层迁移至建筑物内部的主要驱动力，模型将污染物的质量传递速度设置为裂缝处达西速度的垂向分量，从而获得污染源的通量大小。具体计算过程如下：假设VOCs在水相-气相之间达到了分配平衡，则地下水中VOCs的浓度与土壤气中VOCs的浓度符合亨利定律：iigiwcHc=（5）式中：igc土壤气中污染物的浓度；Hi污染物的亨利常数；iwc地下水中污染物的浓度。根据环境科学导刊第 42 卷第 4 期2023 年 8 月8亨利定律，求解土壤气中VOCs的源浓度。根据达西定律模块模拟获得的土壤气体速度场，计算地基裂缝处土壤气达西速度在垂直方向的平均速度，作为进入建筑物室内的土壤气中VOCs的质量传递速度。2 案例研究2.1 场地概况ASU House是一栋两层的分体式房屋，位于朝南的斜坡上，占地面积约为85 m2。该建筑物下方2.7 m处存在受氯化溶剂污染的地下水羽流，主要污染物包括：1,1,1-三氯乙烷、三氯乙烯（TCE）以及相关的降解子产物（如1,1-二氯乙烯）。本研究关注的目标污染物为三氯乙烯。根据过往四年的研究，地下水中TCE的浓度范围为1050 g/L，平均值为249 g/L。ASU House内部结构见图2，模型参数见表1，其中裂缝位置及土层厚度为场地实际调研获得，裂缝长度、裂缝深度、土层渗透率及孔隙率为文献经验值。ABCDA 车库C 洗手间裂缝D 卧室B 洗衣房图2 ASU House结构示意图2.2 模拟过程地下水中TCE浓度为24 g/L，根据文献调研，25时TCE的亨利常数为0.050916，因此热力学平衡时，土壤气中TCE蒸气的浓度为9.30E-06 mol/m3。根据模拟土壤气体流场模拟结果（图3），地基底部土壤气达西速度在垂直方向的平均速度为 9.12E-04 m/s，最终计算得到进入建筑物室内的TCE通量为7.63E-12 mol/（m2s）。污染物浓度模拟最终结果见图4。为确保模拟结果准确性，使用污染物浓度实测数据对结果进行验证。根据11个土壤气监测点位处的模拟浓度与实测浓度对比，计算得到两组数据的MAPE为0.42，MAE为1.38E-06，MSE为6.17E-12，RMSE为2.48E-06。结果（见表2）表明，模型预测值与场地实测数据吻合度良好，其中靠近地基处的监测点位具有更高的吻合度，因此认为该模型可以用于负压控制条件下地层中VOCs蒸气的分布模拟。而造成模拟差异的主要原因有：场地实验数据量不足；降解等过程的简化；裂缝等相关参数缺少实际数据而取自文献经验值；模型计算的迭代限制17等。2.3 VOCs时空分布规律分析根据模型模拟结果，以点位2为例，各时间点的TCE浓度模拟结果见表3，在负压施加的过程中，TCE浓度整体呈平稳上升趋势，增速先上升后下降。在t=76 h，浓度达到1.92E-06 mol/m3，与基线状态相比，污染物浓度波动性下降，蒸气入表1 模型参数参数数值模拟域长度/m50模拟域宽度/m50模拟域深度/m2.7地基长度/m11.3地基宽度/m7.3地基深度/m0.2污染源浓度/（g/L）24裂缝长度/m0.3裂缝深度/m0.001裂缝处压强差/Pa-2.3-4.2-4.4有效扩散系数/m2/s4.5E-7土层厚度(a)/m0.9渗透率(a)/m21.02E-7孔隙率(a)0.375土层厚度(b)0.9渗透率(b)1.79E-9孔隙率(b)0.481土层厚度(c)/m0.7渗透率(c)/m21.79E-9孔隙率(c)0.481流线：达西速度场4040202000yzx1mmm10-42018161412108642图3 达西速度场模拟结果负压条件下地层中 VOCs 污染的数值模拟分析程昕9侵效应显著增强。污染物浓度随时间的变化规律如图5所示。因此，建筑物压力控制技术可以用于污染场地蒸气入侵最坏情形的调查评估。点位1,2,3（a组）与点位9,10,11（b组）是两组平面位置相同而深度不同的采样点，横向纵向对比该两组点位的污染物浓度模拟结果（见表4与图6），分析污染物浓度的空间分布规律。水平方向上，a组点位浓度大于b组，浅层土壤处略有反常。对比两组点位的空间位置，与b组相比，a组与裂缝（即负压中心）之间的距离更小，推测这是a组浓度整体高于b组的原因。浅层土壤处，土壤气流速加大，土壤气中的污染物通过裂缝进入建筑物室内空间，这可能是点位1污染物浓度小于点位9的原因。因此，在负压控制下的蒸气入侵污染场地调查中，应重点关注裂缝附近的污染情况，在裂缝处加密布置采样点位，以获取最坏情形下的评估结果。在垂直方向上，三层采样深度分别为0.7 m、1.6 m和2.5 m，随着深度的增加，污染物浓度逐渐提高。分析其原因为：流线：总通量 流线 Color：浓度/（mol/m3）10-6987654321040402020001mmmyzx图4 污染物浓度模拟结果表2 模型验证结果点位 模拟浓度/（mol/m3）实测平均浓度/（mol/m3）百分误差11.76E-081.64E-080.08 21.92E-061.17E-060.64 37.16E-064.80E-060.49 41.29E-081.62E-08-0.20 53.52E-094.07E-09-0.13 64.97E-081.01E-07-0.51 74.36E-086.26E-08-0.30 85.56E-064.11E-060.35 94.45E-088.71E-08-0.49 108.33E-074.57E-06-0.82 115.14E-061.19E-05-0.57 MAPE（平均绝对百分比误差）0.42MAE（平均绝对值误差）1.38E-06MSE（均方误差）6.17E-12RMSE（均方根误差）2.48E-06表3 点位2处TCE浓度模拟结果时间/h浓度/（mol/m3）06.00E-1046.46E-0982.11E-08121.83E-08162.15E-08202.51E-08244.81E-08288.84E-08321.35E-07362.28E-07403.21E-07444.50E-07485.93E-07527.48E-07569.28E-07601.11E-06641.31E-06681.51E-06721.71E-06761.92E-06时间/hR2=0.99632.50E-062.00E-061.50E-061.00E-065.00E-070.00E+00TCE浓度/（mol/m3）04812 16 20 24 28 32 36 40 44 48 52 56 60 64 68 72 76图5 点位2污染物浓度模拟结果表4 a组与b组TCE浓度模拟结果点位（a组）浓度/（mol/m3）点位（b组）浓度/（mol/m3）11.76E-0894.45E-0821.92E-06108.33E-0737.16E-06115.14E-06R2=0.9048点位点位点位点位点位点位R2=0.9821浓度浅层点位中层点位深层点位点位图6 a组与b组污染物浓度模拟结果深层土壤更靠近地下水污染源，导致土壤气中TCE浓度较高，且靠近地基处的污染物通过裂缝环境科学导刊第 42 卷第 4 期2023 年 8 月10进入室内空气，导致地基处土壤气中污染物浓度降低。因此在蒸气入侵污染场地调查中，距地下污染源垂向距离与裂缝位置是需要重点关注评价的指标。3 结论与讨论3.1 结论（1）负压条件下地层中VOCs的时间分布规律：在施加负压的过程中，土壤气中污染物浓度平稳上升，增速先上升后下降，最终模拟浓度为1.92E-06 mol/m3（t=76 h），与基线状态相比，污染物浓度波动性下降，蒸气入侵效应增强。因此，建筑物压力控制技术可以用于污染场地蒸气入侵最坏情形的调查评估。（2）负压条件下地层中VOCs的空间分布规律：水平方向上，与裂缝（即负压中心）之间的距离越小，污染物浓度越高。浅层土壤处，土壤气流速加大，污染物通过裂缝进入建筑物室内空间，浓度略有下降；垂直方向上，深层土壤更靠近地下水污染源，深度越大，土壤气中的污染物浓度越高。靠近地基处的污染物通过裂缝进入室内空气，导致地基处土壤气中的污染物浓度降低。因此，在负压控制下的蒸气入侵污染场地调查中，应重点关注裂缝位置与距地下污染源的垂向距离，在裂缝附近加密布置采样点位，以获取最坏情形下的评估结果。（3）指导蒸气入侵污染场地调查与风险评估：根据11个点位处的实际采样数据可知，本研究建立的模型与实测数据吻合度良好；分析模拟结果可知，通过负压控制，能够消除VOCs的时空异质性，获取最坏情况下的污染数据。因此，实际场调时，可以在负压条件下进行VOCs采样，能够更高效、准确地评估蒸气入侵的风险。同时在蒸气入侵污染场地调查中，距地下污染源垂向距离与裂缝位置是需要重点关注评价的指标。3.2 不足与展望本研究构建的模型适用于氯代烃污染场地，在模拟的过程中，没有考虑生物降解、吸附解析等作用。同时模型的应用需要获取孔隙率、渗透率、地基裂缝等相关参数，若缺少实际数据而参考文献经验值则会带来一定的误差。未来可以进一步研究场地特征对负压执行效果的影响，进而更精确地分析该技术的有效性，总结相应的执行指南（包括压力控制时间、压力控制方式、污染物采样点位、污染物采样方式、蒸气入侵风险评估指标等），更好地为该技术在污染场地调查与风险评估中的实际应用提供技术指导。参考文献：1 Usepa.OSWER technical guide for assessing and mitigating the vapor intrusion to indoor air pathway from groundwater and soils(subsurface vapor intrusion guidance)R.2015.2 常健铭.基于景观都市主义理论的工业遗址改造与更新研究D.长沙:湖南大学,2013.3 中国生态环境部.建设用地土壤污染风险评估技术导则:HJ 25.32019S.2019.4 中国生态环境部.建设用地土壤污染状况调查技术导则:HJ 25.12019S.2019.5 马杰.污染场地VOCs蒸气入侵风险评估与管控M.北京:科学出版社,2020.6 Mchugh T,Loll P,Eklund B.Recent advances in vapor intru-sion site investigationsJ.Journal of Environmental Manage-ment,2017(204):783-792.7 Mchugh T E,Beckley L,Bailey D,et al.Evaluation of vapor intru-sion using controlled building pressureJ.Environmental Science&Technology,2012,46(9):4792-4799.8 Holton C,Guo Y,Luo H,et al.Long-term evaluation of the controlled pressure method for assessment of the vapor intrusion pathwayJ.Environmental Science&Technology,2015,49(4):2091-8.9 Guo Y,Holton C,Luo H,et al.Identification of alternative vapor intrusion pathways using controlled pressure testing,Soil Gas Monitoring,and Screening Model CalculationsJ.Environmental Science&Technology,2015,49(22):13472-82.10 Guo Y,Dahlen P,Johnson P C.Development and validation of a con-trolled pressure method test protocol for vapor intrusion pathway assessmentJ.Environmental Science&Technology,2020,54(12):7117-7125.11 Song S,Schnorr B A,Ramacciotti F C.Quantifying the influence of stack and wind effects on vapor intrusionJ.Human and Ecologi-cal Risk Assessment,2014,20(5):1345-1358.12 Yao Y,Xiao Y,Luo J,et al.High-frequency fluctuations of indoor pressure:A potential driving force for vapor intrusion in urban areasJ.Science of the Total Environment,2020(710).13 Millington R,Quirk J P.Permeability of porous solidsJ.Transac-tions of the Faraday Society,1961,57(8):1200-1205.14 Lutes C C,Holton C W,Truesdale R,et al.Key design elements of building pressure cycling for evaluating vapor intrusion-a literature reviewJ.Ground Water Monitoring and Remedia-tion,2019,39(1):66-72.15 Millington R,Quirk J P.Permeability of porous solidsJ.Transac-tions of the Faraday Society,1961,57(8):1200-1205.16 吕燕,张巍,黄流雅,等.挥发性有机污染物水中亨利常数的简易测定方法J.华东理工大学学报(自然科学版),2009,35(6):860-865.17 刘妍秋.建筑物室内压强控制技术在蒸气入侵调查中的应用研究D.杭州:浙江大学,2021.（下转第76页）环境科学导刊第 42 卷第 4 期2023 年 8 月7614 白洁,王欢欢,刘世存,等.流域水环境承载力评价以白洋淀流域为例J.农业环境科学学报,2020,39(5):1070-1076.15 傅聪颖,赖昭豪,郭熙.基于熵权TOPSIS模型的区域资源环境承载力评价及障碍因素诊断J.生态经济,2020,36(1):198-204.16 Stephen C.Mangi,Callum M.Roberts,Lynda D.Rodwell.Reef fisheries management in Kenya:Preliminary approach using the driver-pressure-state-impacts-response(DPSIR)scheme of indicatorsJ.Ocean&coastal management,2007,50(5/6).Evaluation and Prediction of Carrying Capacity of Water Resources and Environment in Xian CitySI Xun-lian,ZHAO Pei-yu,DONG Zhong-hui(School of Economics and Management,Xian Shiyou University,Xian Shanxi 710065,China)Abstract:Based on the connotation of carrying capacity of water resources and environment system,the index system of carrying capacity of water resources and environment in Xi an was constructed from four aspects of water resources,water environment,society and economy.The TOPSIS method and fuzzy comprehensive evaluation method were adopted to evaluate the carrying capacity of water resources and environment in Xi an from 2006 to 2020.Based on the evaluation results,GM(1,1)model was used to predict the change trend of carrying capacity of water resources and environment the next two years.The results showed that the carrying capacity of water resources and environment in Xi an fluctuated and rised,and the comprehensive scores of each subsystem all rised on the whole from 2006 to 2020.In the next two years,the carrying capacity of water resources and environment would increase year by year,and both of them have reached a high-quality level.The water environment and economic subsystem were green lights,and the social subsystem would rise to a green light in 2022,however the level of water resources subsystem was still worrying.Key words:carrying capacity of water resources and environment;TOPSIS method;fuzzy comprehensive evaluation;GM（1,1）modelNumerical Simulation of the VOCs Underground in Negative Pressure ConditionCHENG Xin1,3,ZHAO Sheng-hao2,QIAN Yu1,3(1.State Key Laboratory of Pollution Control and Resource Reuse,Nanjing Jiangsu 210023,China)Abstract:A three-dimensional numerical simulation was conducted using ASU House as the site case,and the Darcys law and the law of dilute matter transfer in porous media were used to simulate the soil gas velocity field and the flux of VOCs into the building,so as to obtain the spatial and temporal distribution law of VOCs in the stratum under negative pressure conditions.The results showed that the pollutant concentration increases significantly under the negative pressure condition and the vapor intrusion process was enhanced compared with the baseline condition.The closer the distance from the negative pressure center to the pollution source,the higher the pollutant concentration was.Key words:VOCs;vapor intrusion;negative pressure control;numerical simulation（上接第 10 页）