岩溶地下河污染物运移模型对比研究

两区模型（TRM）和暂时存储模型（TSM）为两种常用的岩溶地下河（管道）污染物运移模型,为了更好表征岩溶地下河污染物运移过程,开展模型对比研究.采用两种模型拟合室内管道示踪实验穿透曲线,从概念模型框架、模型优缺点、参数校正过程和物理意义、参数敏感性等方面对比分析两种模型的异同.结果表明,OTIS软件提供TSM的数值解,而CXTFIT软件提供TRM的解析解,但两者均能通过下游不同位置的穿透曲线分段表征地下河不同阶段的污染物运移过程.两种模型的物理意义相近,均可近似描述岩溶地下河中溶质的暂时性存储现象,较好地拟合溶质运移曲线.采用TSM可校正得到暂时存储区的横截面积,但采用TRM只能获得非流动区所占的比例.两种模型中弥散系数和传质/交换系数的敏感度较低,并且弥散系数的变化对穿透曲线拖尾无显著影响.研究结果为岩溶地下河污染物运移模型的选取与应用提供了依据.     

胡敏酸在土壤中的迁移

研究了不同条件下胡敏酸在土壤中的迁移过程.结果表明,其迁移过程表现出较大的差异,在溶液pH值、胡敏酸浓度和溶液流速均较高时,胡敏酸迁移的阻滞因子和在土壤中的吸附系数均较低,有利于胡敏酸在土壤中的迁移;在溶液pH值、胡敏酸浓度和溶液流速均较低时,土壤对胡敏酸迁移的阻滞作用较大.胡敏酸的迁移还与土壤性质有关,土壤黏粒含量和阳离子交换量越高,胡敏酸在土壤中吸附的越多,越不利于迁移.在土壤中的吸附是胡敏酸迁移的主要控制因素,在预测胡敏酸在水土环境中迁移时,应该充分考虑不同条件下胡敏酸在不同性质土壤中的吸附.     

旅游和区域大气污染对四川九寨沟气溶胶的贡献

2010年4月~2011年4月,在九寨沟连续监测了总悬浮颗粒物(TSP)和水溶性无机离子的浓度.结果表明:旅游活动显著增加了空气中TSP、SO42-、Ca2+、K+、NH4+和NO3-的含量,而Na+、Mg2+和Cl-则主要来自自然源;旅游强度最高时期(6~10月),降水对空气的清洗作用最强,气溶胶污染程度为全年最低;旅游强度较低时期(1~3、4~5和11~12月),降水量较低且西北沙尘易在春季抵达九寨沟,所以气溶胶污染程度较高;区域燃煤可能是气溶胶[SO42-]:[NO3-]比值较高的原因;生物质燃烧可能是K+在秋末至次年初春较高的原因之一.九寨沟大气环境已明显受当地旅游活动和污染物长距离传输的影响.     

四川盆地近地面风场及污染物输送通道统计分析

利用2017年四川盆地18个城市PM_2.5小时浓度数据、欧洲气象中心ERA-Interim再分析资料对四川盆地污染特征及不同季节近地面风场特征进行统计分析.通过分异指数和相关性分析各城市PM_2.5逐小时浓度，结合近地面风场特征分析出四川盆地污染物变化趋势最相似的城市组合，探索污染物传输特征.研究发现，分异指数和相关性分析得到的传输通道城市组合与地面风场基本相符.四川盆地污染物输送途径可能包括以下3条路径，第1条"川西通道"：污染物随气流沿广元→绵阳东南部→德阳→成都、眉山北部→雅安流动；第2条"川中通道"：污染物随气流沿巴中→南充北部→遂宁北部、绵阳东南部→资阳北部→眉山东部→乐山北部流动；第3条"川东通道"：污染物随气流沿重庆北部→达州、广安→南充南部→遂宁中部→资阳东部→内江、自贡→宜宾、泸州流动.川西通道城市污染排放量大，容易引起连片污染，对应城市群应实施联防联控；川东通道末端在川南城市群和重庆形成风场辐合，造成污染物的滞留和累积，因此，建议在中长期产业布局中减少川南城市群的重污染企业.     

Comparison of Microscopic and Macroscopic Modeling Approaches for Subsurface Contaminant Transport

The practice of contaminant transport and remediation has shown significant progress in recent years. However, despite the significant progress made, remediation efforts are often delayed by extremely long breakthrough curve tails that render efforts to bring the level of contaminants below the regulatory standards inefficient. One hypothesis is that these long tails are due to the reservoir-like slow diffusive processes in soil micropore zones. This study compares the effects of micropores at macroscopic and microscopic levels and establishes a link between these approaches for validation and calibration purposes. The link between macroscopic and microscopic levels is established through comparisons and testing of the two models while incorporating appropriate scale and boundary effects. Despite the differences in conceptual approaches and simulation time, the two approaches rendered meaningful results. The link helps forecast the effects of micropore zone transport processes in the subsurface efficiently and thus allows development of numerical tools that could contribute towards more efficient remediation design.     

Experimental Study of the Criteria of Flow Laminarization in Two-Dimensional Dense Gas Plumes

Laminarization of flow in a two-dimensional dense gas plume was experimentally investigated in this study. The plume was created by releasing CO₂ through a ground-level line source into a simulated turbulent boundary layer over an aerodynamically rough surface in a meteorological wind tunnel. The bulk Richardson number (Ri_*), based on negative plume buoyancy, plume thickness, and friction velocity, was varied over a wide range so that the effects of stable stratification on plume laminarization could be observed. A variety of ambient wind speeds as well as three different sizes of roughness arrays were used so that possible effects of roughness Reynolds number (Re_*) on plume laminarization could also be identified. Both flow visualization methods and quantitative measurements of velocity and intermittency of turbulence were used to provide quantitative assessments of plume laminarization.Flow visualization provided an overall picture of how the plume was affected by the negative buoyancy. With increasing Ri_*, both the plume depth and the vertical mixing were significantly suppressed, while upstream propagation of the plume from the source was enhanced. The most important feature of the flow revealed by visualization was the laminarization of flow in the lower part of the plume, which appeared to be closely related to both Ri_* and Re_*.Measurements within the simulated dense gas plumes revealed the influence of the stable stratification on mean velocity and turbulence intensity profiles. Both the mean velocity and turbulence intensity were significantly reduced near the surface; and these reductions systematically depended on Ri_*. The roughness Reynolds number also had considerable influence on the mean flow and turbulence structure of the dense gas plumes.An intermittency analysis technique was developed and applied to the digitized instantaneous velocity signals. It not only confirmed the general flow picture within the dense plume indicated by the flow visualization, but also clearly demonstrated the changes of flow regime with variations in Ri_* and Re_*. Most importantly, based on this intermittency analysis, simple criteria for characterizing different flow regimes are formulated; these may be useful in predicting when plume laminarization might occur.     

Validation of the operational emergency response model at the Swedish Meteorological and Hydrological Institute using data from etex and the chernobyl accident

The Eulerian atmospheric tracer transport model MATCH (Multiscale Atmospheric Transport and Chemistry model) has been extended with a Lagrangian particle model treating the initial dispersion of pollutants from point sources. The model has been implemented at the Swedish Meteorological and Hydrological Institute in an emergency response system for nuclear accidents and can be activated on short notice to provide forecast concentration and deposition fields.The model has been used to simulate the transport of the inert tracer released during the ETEX experiment and the transport and deposition of ¹³⁷Cs from the Chernobyl accident. Visual inspection of the results as well as statistical analysis shows that the extent, time of arrival and duration of the tracer cloud, is in good agreement with the observations for both cases, with a tendency towards over-prediction for the first ETEX release. For the Chernobyl case the simulated deposition pattern over Scandinavia and over Europe as a whole agrees with observations when observed precipitation is used in the simulation. When model calculated precipitation is used, the quality of the simulation is reduced significantly and the model fails to predict major features of the observed deposition field.     

应用TOXI模型进行重金属在河流中迁移转化规律的研究

ＴＯＸＩ模型是美国环保局的ＷＡＳＰ（水质分析模拟程序）中的一部分。由于该模型考虑了污染物质在河流中溶解态，颗粒态和底泥态的分布关系，对于模拟重金属在河流中的浓度变化效果较好。为保证计算结果的可靠性，首先应对重金属浓度的衰减变化特性进行研究，本文重点讨论的是重金属在河流中迁移转化的过程，并如何将这一过程用数学公式来进行模拟。     

土壤中水动力弥散系数的研究进展

从溶质运移的基本方程出发 ,系统论述了国内外对土壤中水动力弥散系数测定方法的研究进展 ,并指出了上述各种方法的适用条件及存在的问题     

Transport of volatile compounds in porous media in the presence of a trapped gas phase

The presence of an immobile gaseous phase in an otherwise-saturated porous medium affects the transport of volatile compounds. The linear theory of partitioning tracers suggests that a volatile tracer introduced into such a system should be retarded with a constant retardation factor. Using high concentrations, however, the saturation of the gaseous phase will change as an effect of the tracer test itself. Competitive gas transfer among all volatile compounds and the change of saturation may lead to tracer concentrations that are temporarily higher than those injected. We analyze the system in the framework of the coherence theory by Helfferich [Soc. Pet. Eng. J. 21 (1) (1981) 51]. The governing equations are formulated as functions of total concentration, i.e., the mass of solutes in all phases per pore volume. Neglecting dispersion and mass-transfer kinetics, we derive the characteristic form of the resulting system of hyperbolic equations. In a system with N volatile compounds, a variation of the concentrations splits up into N waves, each traveling with its own characteristic velocity. If the presence of a gaseous phase is sustained, one wave will be a standing one. We perform numerical model calculations for tracers with various Henry's law coefficients and show that the results agree with the semi-analytical solution obtained by coherence theory.