典型土壤中NAPLs态甲苯运移特性剖析与模拟研究

典型土壤中NAPLs态甲苯运移特性剖析与模拟研究摘要甲苯属于轻非水相液体（LNAPLs，NAPLs大类），密度小于水、难溶于水、挥发性强，是加油站、石化场地地下水土典型污染物。当发生泄漏后，NAPL相甲苯在包气带、毛细带、潜水面附近发生多相渗流+相间传质，同时伴随固相吸附、微生物降解，形成长期“持续污染源”。本文剖析NAPL甲苯在典型土壤（砂土、粉土、黏土）的运移机制、主控影响因素，建立多相流-溶质耦合数学模型，开展数值模拟，揭示其迁移分布规律，为污染场地风险评估、修复方案设计提供理论支撑。关键词：LNAPLs；甲苯；包气带；多相渗流；运移特性；TMVOC/GMS/COMSOL；数值模拟1绪论1.1研究背景与意义地下储油罐、油品输送管道泄漏会释放含甲苯的NAPL。NAPL相甲苯不会快速溶于地下水，会在土壤孔隙中形成自由相，在重力、毛细力共同作用向下迁移；到达潜水面后漂浮于水面，长期缓慢溶解、挥发，持续污染地下水、土壤气，存在呼吸暴露、地下水饮用风险。

传统单一溶质对流弥散模型无法描述自由NAPL相流动，必须采用水-NAPL-气三相多相流耦合模型刻画全过程。1.2国内外研究现状现有研究多利用土柱、砂箱试验观测LNAPL运移；数值模拟常采用TOUGH2-TMVOC、GMS、COMSOL等，研究水位波动、介质非均质性对LNAPL分布的影响。现有研究对粉土、黏土夹层、透镜体条件下NAPL甲苯滞留、残余饱和度、相间传质的定量刻画仍有待完善。1.3研究内容与技术路线测定典型土壤（中砂、粉土、黏土）基础物性参数与甲苯物化参数；剖析NAPL甲苯在包气带→毛细带→含水层全过程运移机理；构建水-NAPL-气三相耦合数学模型；数值模拟，分析介质类型、泄漏量、水位波动对甲苯分布的影响；敏感性分析，识别主控因子，提出场地管控建议。2NAPL态甲苯理化性质与典型土壤介质特征2.1甲苯基础物化参数（25℃）密度：0.867g/cm³＜水，属于LNAPL，浮于潜水面水中溶解度：约526mg/L（低溶解度）黏度：0.59mPa・s，流动性较好亨利常数：0.27，挥发性较强有机碳吸附系数Koc2.2典型土壤介质特征砂土：大孔隙、高渗透系数、毛细压力小；NAPL向下迁移快，易形成连续自由相，残余NAPL饱和度低；粉土：中等孔隙，毛细作用显著，NAPL易在粉土层大量滞留，形成残余油；黏土：极低渗透性，毛细压力大；NAPL难以穿透，多在黏土层上方积聚，横向铺展。介质非均质性（透镜体、夹层）会改变NAPL运移路径，优先沿高渗通道迁移，在低渗界面大量截留。3NAPL态甲苯在土壤中的运移特性剖析NAPL甲苯进入地下多孔介质后，存在4种赋存相：自由NAPL相、水溶解相、土壤气相、固相吸附相，多过程耦合。3.1运移物理机制自由NAPL相多相渗流

泄漏初期，在重力主导下垂向入渗；孔隙毛细阻力阻碍NAPL进入微小孔隙。当NAPL饱和度低于残余饱和度，甲苯不再流动，固定残留在孔隙，成为长期污染源。到达潜水面后，漂浮在水面之上，沿水力梯度侧向铺展，形成浮油层。相间传质过程溶解：自由NAPL相缓慢溶解进入孔隙水，形成溶解羽，随地下水对流弥散，是地下水污染的核心来源；挥发：NAPL甲苯挥发进入土壤气相，气相扩散，可向上迁移至地表形成土壤气污染；气压波动、水位升降会显著强化挥发传质；吸附/解吸：水相中甲苯吸附于土壤有机质表面；土壤含水率升高会降低固相吸附容量。衰减转化

好氧条件下，甲苯可被土著微生物降解；厌氧降解速率慢，在污染核心区降解作用有限。3.2主控影响因素介质属性：渗透系数、孔隙分布、有机质含量、毛细压力-饱和度关系（vanGenuchten模型）；NAPL自身属性：黏度、密度、溶解度、挥发性；水力条件：地下水水位波动、降雨入渗、水力梯度；水位升降会反复冲刷残余NAPL，促进溶解与再迁移；泄漏条件：泄漏总量、泄漏速率、泄漏持续时间；少量泄漏会全部转化为残余NAPL，不会到达潜水面；大量泄漏才能形成连续自由相。3.3不同土层中运移行为差异包气带（非饱和带）：重力向下迁移，一部分挥发，一部分被土壤吸附；细颗粒土层截留效应强；毛细上升带：毛细水饱和，NAPL在毛细带顶部大量聚集，是LNAPL污染的核心滞留区；饱和含水层：自由NAPL漂浮在潜水面之上；溶解相在水面以下形成污染羽，随地下水侧向扩散。4数学模型构建4.1基本假设①等温条件；②多孔介质均质/分层，骨架不可压缩；③三相（水w、NAPLn、气相g）渗流满足达西定律；④相间传质采用一级动力学；吸附采用线性分配模型；⑤微生物降解为一级衰减。4.2多相流控制方程（质量守恒）qβ=−kkrβ(Sβ)μβ(∇P4.3组分溶质运移方程（甲苯组分）源汇项包含：NAPL-水溶解、NAPL-气相挥发、固相吸附、微生物降解。4.4模型参数介质参数：孔隙度、渗透系数、vanGenuchten参数；甲苯参数：密度、黏度、溶解度、亨利常数、Kd动力学参数：相间传质系数、降解速率常数。4.5数值求解方案采用积分有限差分/有限元求解，推荐软件：TOUGH2-TMVOC（多相NAPL首选）、GMS（场地尺度）、COMSOL多物理场。

边界设置：上边界大气开放；底部隔水；侧向定水头；泄漏源设置为恒质量源。5模拟结果与分析5.1典型运移演化过程入渗阶段：NAPL向下快速垂向迁移，在粗砂层形成连续油相；遇到粉土/黏土夹层，垂向受阻，侧向铺展；聚集阶段：NAPL抵达毛细带，在潜水面上方累积，浮油范围持续扩大；孔隙中形成大量残余NAPL；长期缓释阶段：自由NAPL逐渐消耗，但残余NAPL持续溶解、挥发，溶解污染羽缓慢向下游扩散；即使停止泄漏，污染仍可持续数年～数十年。5.2敏感性分析结论土壤渗透系数是最敏感参数，高渗砂土迁移距离显著大于黏土；水位波动会周期性改变NAPL饱和度，增强溶解，扩大地下水污染羽范围；有机质含量升高，固相吸附增强，降低水相甲苯浓度，但不改变自由NAPL渗流路径；残余饱和度决定污染源存续时间，残余饱和度越高，污染持续时间越长。6结论与建议6.1主要结论NAPL态甲苯作为LNAPL，泄漏后在包气带以重力垂向迁移为主，到达潜水面漂浮聚集；土壤孔隙内残留大量残余NAPL，是长期缓释污染源；介质分层/透镜体强烈改变运移路径；砂土利于NAPL迁移，粉黏土截留作用显著；水位波动、降雨入渗强化相间传质，促进甲苯溶解、挥发，加剧土壤气与地下水污染；多相流耦合模型可较好刻画自由NAPL流动+溶解+挥发+吸附全过程；单一水相溶质模型会低估污染风险与污染持续时间。6.2工程建议场地调查不能仅检测地下水，必须识别包气带残余NAPL分布，重点探查毛细带；修复方案：自由相优先抽提（LNAPL回收）；包气带可采用SVE土壤气相抽提；地下水污染采用PRB、生物修复；风险模拟推荐采用TMVOC，必须考虑残余饱和度、相间传质，不可简化为单相溶质模型。6.3不足与展望本研究可进一步考虑介质空间随机非均质性、非等温效应、多组分BTEX共溶效应，开展原位场地校准。参考文献（示例）[1]王洪涛。多孔介质多相流与污染物运移[M