梯度模式原理及其应用

本文推导梯度模式，引用大气边界层湍流系数进行实际运算，与现场大气环境监测结果对照；同时与高斯模式计算结果比较，得出该模式在城市复杂流场的实用性．     

烟囱排烟扩散的影响因素初步研究

为减少和控制工业锅炉烟囱排放中的有害气体对自然生态环境的污染，设计单位的专业工程设计人员在锅炉房设计布置中，除了应根据锅炉的蒸发量来确定风，烟道的截面外，还须根据不同地区的气象条件，局部地形等实际情况，来决定锅炉房烟囱的设计。文章简妻介绍了风速、风向，湍流强度及大气稳定度对烟囱排烟扩散的影响。     

Estimation of far-field horizontal and vertical turbulent diffusion coefficients from the concentration field of a wastewater plume near the Akashi Strait

We carried out a field study of the plume discharged by a near-shore wastewater outfall near the Akashi Strait, Japan. Using an Acoustic Doppler Current Profiler and a tow-body CTD, we measured the near-surface salinity and temperature fields in the region throughout an M₂ tidal cycle. We filtered the data in T–S space to remove water masses other than the wastewater, and then used the adiabatic mixing assumption to calculate the concentration of wastewater in the far field of this plume. Averaging the T–S fields of repeated surveys over a time period during which the tidal regime did not change substantially, allowed comparison of the time-averaged plume with the analytical solution for a plume diffusing in both the horizontal and vertical dimensions. The resulting vertical turbulent diffusion coefficients agreed well with those resulting from Thorpe scales determined via a vertically-profiling CTD, as well as with the canonical value for open channel flow of D _z = 0.067hu _*. The corresponding horizontal turbulent diffusion coefficients, however, were two orders of magnitude larger than those typically observed in straight channels, and an order of magnitude larger than those observed in meandering rivers. This is likely a result of enhanced horizontal mixing due to barotropic eddies generated by the interaction of strong tidal flow with headlands and levees, as well as due to the time-varying nature of tidal flow, and baroclinic spreading of the buoyant wastewater plume.     

Field study of TCE diffusion profiles below DNAPL to assess aquitard integrity

An area where a free-product accumulation of trichloroethylene (TCE) dense non-aqueous phase liquid (DNAPL) occurs at the bottom of a 10-m-thick surficial sand aquifer was studied to determine the integrity of the underlying, 20-m-thick, clayey silt aquitard formed of glaciolacustrine sediment. TCE concentration-versus-depth profiles determined from aquitard cores collected at five locations indicated penetration of detectable TCE 2.5 to 3.0 m into the aquitard. Two of the profiles show persistent DNAPL at the aquitard interface, while two others indicate that DNAPL, present initially, was completely dissolved away producing concentration declines at the aquitard interface. The fifth profile suggests shallow DNAPL penetration (<0.5 m) into the aquitard, however, this penetration, which was likely caused by cross-contamination during core collection or cone penetrometry (CPT) of the aquitard interface, did not increase the maximum depth of TCE penetration. Combining the field profiles with one-dimensional model simulations, downward migration of the aqueous TCE front, defined as the EPA MCL of 5 microg/l, which was below the analytical detection limit, was projected to a distance between 4 and 5 m below the top of the aquitard. Using a single set of estimated aquitard parameter values, simulations of aqueous TCE migration into the aquitard provided a good fit to four of the field profiles with a migration time of 35 to 45 years, consistent with the history of TCE use at the site. These simulations indicate aqueous TCE migration is diffusion-dominated with only small advective influence by the downward groundwater velocity of 2 to 3 cm/year or less in the aquitard due to pumping of the underlying aquifer to supply water to the facility in the past 50 years. The applicability of the parameter values was confirmed by in situ diffusion experiments of 1-year duration, in which stainless steel cylinders containing DNAPL were inserted into the aquitard. The diffusion-dominated nature of the profiles indicates that the aquitard provides long-term protection of the underlying aquifer from contamination from this DNAPL zone. Simulations of long-term migration of the TCE solute front indicate breakthrough to the lower aquifer at 1200 years for the no advection scenario and at 500 years if the strong downward hydraulic gradient persists. However, even after breakthrough, the mass flux through the aquitard to the underlying aquifer remains relatively low, and when considered in terms of potential impacts to pumping wells, concentrations are not expected to increase significantly above present-day MCLs. The use of contaminant profiles of different time and distance scales, in addition to hydraulic data, dramatically improves the ability to assess aquitard integrity, and provides improved transport parameter values for estimating contaminant arrival times and fluxes. The apparent lack of deep preferential pathways for TCE migration, such as open fractures, is probably due to the softness of the silty aquitard deposit and minimal physical or chemical weathering of the aquitard provides long-term protection of the underlying aquifer from contamination from this DNAPL zone. Simulations of long-term migration of the TCE solute front indicate breakthrough to the lower aquifer at 1200 years for the no advection scenario and at 500 years if the strong downward hydraulic gradient persists. However, even after the breakthrough, the mass flux through the aquitard to the underlying aquifer remains relatively low, and when considered in terms of potential impacts to pumping wells , concentrations are not expected to increase significantly above present-day MCLs. The use of contaminant profiles of different time and distance scales, in addition to hydraulic data, dramatically improves the ability to assess aquitard integrity, and provides improved transport parameter values for estimating contaminant arrival times and fluxes. The apparent lack of deep preferential pathways for TCE migration, such as open fractures, is probably due to the softness of the silty aquitard deposit and minimal physical or chemical weathering of the aquitard.     

The diffusion and sorption of volatile organic compounds through kaolinitic clayey soils

Laboratory experiments to estimate the effective molecular diffusion coefficient (D(e)) and sorption coefficient (K(d)) for volatile organic compounds through natural clayey soils were conducted using diffusion testing apparatus. The compounds tested were methyl ethyl ketone (MEK), toluene and trichloroethylene (TCE). The D(e) and K(d) values were determined by a curve fitting procedure. The compound losses, and the effects of porous disks used in the apparatus were significant. The transport of MEK was faster than that of TCE and toluene because of the lower sorption to the soils. The D(e) values of all the compounds were of the order of 10(-10) m(2)/s and smaller than the diffusion coefficient in pure aqueous solution at infinite dilution (D(0)), due to the tortuosity of the samples. The effects of the sample thickness on the parameter determination were not significant. Comparison to the K(d) values estimated from batch sorption tests and from organic carbon content (f(oc))-based predictions showed that the diffusion test results were intermediate between those from the other two methods. The diffusion tests use compacted soil samples and should be more relevant to in situ conditions, but the reliability of the tests is affected by large compound losses that cause uncertainties in their interpretation. It is recommended that more than one method be used to assess K(d) values.     

Characterisation of HTO diffusion properties by an in situ tracer experiment in Opalinus clay at Mont Terri

A long-term single borehole diffusion experiment using tritiated water as tracer was carried out in Opalinus clay, an argillaceous rock formation that is accessible at the Mont Terri Underground Research Laboratory, situated in the Swiss Jura. The tracer was diluted in reconstituted formation water and introduced into a packed-off section of a borehole located in saturated rock. Pressure in this interval was maintained equal to the pore pressure of the surrounding rock in order to prevent any hydraulic gradient around the borehole and to avoid advective transport processes. The evolution of the tracer concentration in the injection system was monitored over time. After 1 year of diffusion, the claystone surrounding the interval was retrieved by overcoring the whole borehole and packer system, and by an adjacent oblique borehole. Compressed air was used as drilling fluid to reduce rock disturbances. The recovered overcore was sampled along profiles perpendicular to the borehole wall with a view to determining the tracer-concentration profiles in the rock. To avoid further evaporation of tritiated water, subsamples were immediately transferred into polyethylene bottles and disaggregated by adding a known amount of tracer-free water. Fifteen profiles were determined and showed a decreasing tracer concentration with distance into the rock. The pore-water contents were constant along those profiles, confirming that only very little water was lost during overcoring operations. The evolution of tritium-tracer concentration in the injection system over time and in situ profiles were interpreted with a 3-D numerical simulation of the experiment. That allowed for the identification of the transport parameters (orthotropic diffusion tensor and porosity) by minimising the relative quadratic error between the experimental and simulated data. The fitting is good and the results are consistent with data obtained on drill-core samples. The result of tritiated water is discussed regarding (1) the potential effect of mechanical and/or chemical disturbances around the injection borehole and (2) the specific behaviour of tritiated water.     

Experimental verification for real-time environmental emergency response system: WSPEEDI by European tracer experiment

Japan Atomic Energy Research Institute has developed an emergency response system WSPEEDI to forecast long-range atmospheric dispersions of radionuclides discharged into the atmosphere. The latest version of WSPEEDI consists of an atmospheric dynamic model MM5 for calculating meteorological fields and a particle random-walk model for atmospheric dispersion. The performance of WSPEEDI was evaluated by data obtained from a field tracer experiment over Europe (ETEX) in this paper. The model validation was done with respect to the following points: (1) the dependence of model accuracy on the temporal and spatial resolutions of the meteorological fields and (2) the superiority of an atmospheric dynamic model over a mass-consistent wind model. Regarding (1), it was shown that the calculation accuracy of the new version with high temporal resolution was improved, especially at the edge of the plume. Moreover, although the increase in horizontal spatial resolution of the old version had no substantial effect on the model performance, increase in horizontal resolution of the new version contributed to the significant improvement of the calculation accuracy. These results showed that the dynamically calculated meteorological field with the spatial resolution of the meso-β–γ scale greatly improved calculation accuracy.     

Diffusion of HTO, Cl and I in Opalinus Clay samples from Mont Terri: Effect of confining pressure

Diffusion coefficients (T=23±2 °C) and accessible porosities for HTO, ³⁶Cl⁻ and ¹²⁵I⁻ were measured on Opalinus Clay (OPA) samples from the Mont Terri Underground Rock Laboratory (URL) using the through-diffusion technique. The direction of transport (diffusion) was perpendicular to bedding. Special cells that allowed the application of confining pressure were designed and constructed. The pressures ranged from 1 to 5 MPa, the latter value simulating the overburden at the Mont Terri URL (about 200 m). The test solution used in the experiments was a synthetic version of the Opalinus Clay pore water, which has Na⁺ and Cl⁻ as the main components (I=0.42 M).The measured values of the effective diffusion coefficients (D_e) and rock capacity factors (α) are: D_e=1.2–1.5×10⁻¹¹ m² s⁻¹ and α=0.09–0.11 for HTO, D_e=4.0–5.5×10⁻¹² m² s⁻¹ and α=0.05 for ³⁶Cl⁻ and D_e=3.2–4.6×10⁻¹² m² s⁻¹ and α=0.07–0.10 for ¹²⁵I⁻. For non-sorbing tracers (HTO, ³⁶Cl) the rock capacity factor α is equal to the diffusion-accessible porosity . The experimental results showed that pressure only had a small effect on the value of the diffusion coefficients. Increasing the pressure from 1 to 5 MPa resulted in a decrease of the diffusion coefficient of 17% for HTO, 28% for ³⁶Cl⁻ and 30% for ¹²⁵I⁻. Moreover, the diffusion coefficients for ³⁶Cl⁻ and ¹²⁵I⁻ are smaller than for HTO, which is consistent with an effect arising from anion exclusion.The diffusion coefficients of HTO and ¹²⁵I⁻ measured in this study are in good agreement with recent measurements at three other laboratories performed within the framework of a laboratory comparison exercise. The values of the diffusion-accessible porosities show a larger degree of scatter.     

Diffusion of tritiated water and Na through non-degraded hardened cement pastes

Diffusion experiments through hardened cement pastes (HCP) using tritiated water (HTO) and ²²Na⁺, considered to be conservative tracers, have been carried out in triplicates in a glove box under a controlled nitrogen atmosphere. Each experiment consisted of a through-diffusion test followed by an out-diffusion test.The experimental data were inversely modelled applying an automated Marquardt–Levenberg procedure. The analysis of the through-diffusion data allowed the extraction of values for the effective diffusion coefficients, D_e, and the rock capacity factor, α. Good agreement between measured and calculated tracer breakthrough curves was achieved using both a simple diffusion model without sorption and a diffusion/linear sorption model. The best-fit K_d-values were found to be consistent with R_d-values measured in previous batch-sorption experiments.The best-fit values from the through-diffusion tests were then used to predict the results of subsequent out-diffusion experiments. Good agreement between experimental data and predictions was achieved only for the case of linear sorption.Isotopic exchange can only partially account for both the amount of tracer taken up in the batch-sorption tests and the measured retardation in the diffusion experiments and, hence, additional mechanisms have to be invoked to explain the data.