二维河流功能区划中水质目标控制点确定方法研究

提出在水污染控制规划中，明确水质控制目标必须同时明确水质目标控制点的位置及该点对应的水质类别，并以二维河流为对象，总结、分析、比较了原有的控制断面确定方法，进一步以饮用水水源地保护、天然游泳场等５种情况为例，具体阐明了在上述情况下如何确定水质目标控制点的位置。     

利用二维表对活性污泥工艺模拟时的流量进行平衡的方法

多点进水和多重回流成为提高活性污泥工艺去除污染物能力的常用方法,这导致活性污泥工艺流程较为复杂。对工艺进行模拟计算时,由于模拟池体数量比实际池体数量更多,流量平衡较为困难。根据工艺流程特征提出了主流线、副流线等概念,结合流量分配比、循环流量等参数,利用二维表进行流量平衡计算,简捷明了,便于检验,大大提高了进行流量平衡的准确性和效率。     

未确知信息对河流型水源保护区划分的影响

基于未确知数学理论,考虑河流环境系统信息的未确知性,定义二维水质模型主要参数为未确知参数,建立未确知二维水质模型,递推计算出保护区划分的合理长度及对应的污染物浓度取值区间及可信度分布,并对未确知二维水质模型中各参数的影响做了比较分析。结果表明,利用未确知二维水质模型,以COD为污染物计算的晋江干流水源二级保护区水域长度为6 220m,该长度能在0.934 0的主观可信度上保证水质安全,而基于确定性水质模型,利用各参数均值计算的保护区长度5 350 m存在95.22%的水质超标风险率。未确知二维水质模型各参数影响分析表明,模型应用的未确知参数及各参数变异性对模型计算结果有较大关系,未确知参数考虑越多,保护区划分长度越长;参数变异系数越小,其对计算结果影响越小。     

垃圾焚烧炉二次配风优化数值模拟

针对城市生活垃圾焚烧发电的排放问题,以重庆市某垃圾焚烧炉为原始模型,在炉膛两侧炉壁的适当位置设置二次配风口,并采用CFD（计算流体动力学）方法对炉膛内气体的二次燃烧过程进行数值模拟。通过观察对比有无二次配风以及二次配风口位置不同时炉膛内气体的温度场、气体在炉膛内的停留时间分布以及炉膛内气体的混合程度和湍动能等,重点分析了二次风在气体燃烧过程中所起的作用,并对2种不同二次配风口位置时抑制二恶英产生的效果进行了评价。通过对垃圾焚烧炉二次风的优化数值模拟,获得了适合本焚烧炉的比较合理的二次配风条件,可为焚烧炉的设计和改进提供一些有价值的参考。     

因子分析及其在河网水质综合评价中的应用研究

对因子分析方法在天然水体水质评价中的应用进行研究,证明在数据较为缺乏的情况下,使用因子分析方法可以损失较少信息,减少评价指标,针对主因子进行分析,使水质评价结果更为合理和有效。同时,利用因子的得分值还可以将不同地区进行比较分析。应用因子分析对河网地区典型河流的水质状况进行分析,发现了增长潜力较大的主要因子,并对该地区的水质状况做出了符合实际的总体评价,表明因子分析在水质评价中具有重要应用价值。     

数值仿真在钱塘江引水入城工程水质预测中的应用

使用二维浅水波方程及传质模型,对工程前期排人钱塘江的废水中污染物的传播及分布进行预测,其预测结果可为工程排污口附近钱塘江水质保护对策制定提供依据。     

万家寨引黄北干线桑干河道输水水质模拟与控制

主要对采用河道输水方案的山西万家寨引黄北干线桑干河道水质进行模拟与控制研究,依据各水体特性建立河流水质数值模型和水库数值模型,运用系统分解方法将整个研究河道分解为若干个河段,并采用情景分析方法生成3个水质模拟情景,分析模拟结果表明,在山阴取水口和怀仁取水口水质达地表水Ⅲ类标准条件下,册田水库的入水水质不达标。应用反演思想生成2个水质控制情景和3个子情景,在水库出水口水质达标的情况下,分析确定沿河各排污口的污染物排放控制量,给出了保证册田水库水质安全的控制方法。     

电除尘器气流分布试验研究和数值计算

结合660 WM机组电除尘器的设计要求,对电除尘器内气流流动进行了数值模拟与试验研究.计算中选择SIMPLE方法求解三维雷诺平均方程和κ-ε两方程湍流模型.研究结果表明,无论电除尘器两室流量分配还是第一电场进口断面气流分布,计算结果与试验吻合较好,数值计算结果可以指导电除尘器设计与模型试验研究.     

环评中应用二维稳态混合衰减模式预测水质

应用二维稳态混合衰减模式预测水质时,由于模式中本底值也参与衰减,往往造成预测结果低于下游本底值的情况,提出新的解决办法,在环评中应用较为成功。     

内循环厌氧反应器多相流流场的三维数值分析

利用数值模拟的方法,引入欧拉双流体多相流模型及标准k-s紊流模型,模拟计算内循环厌氧反应器的三相流三维流场,并通过改变污泥颗粒密度及进水流量,针对固相流速及固含率的变化情况,分析条件的改变对流场的影响。研究结果表明,应用数值模拟方法可以获得内循环厌氧反应器内的流场特征;污泥颗粒密度及进水流量的改变对于反应器内污泥颗粒的流速及分布的均匀性有较为明显的影响。模拟结果对反应器的应用及优化设计具有一定的参考价值。     

The distribution characteristics of pollutants released at different cross-sectional positions of a river

The distribution characteristics of heavier or lighter pollutants released at different cross-sectional positions of a wide river is investigated with a well-tested three-dimensional numerical model of gravity flows based on Reynolds-Averaged Navier-Stokes equations and turbulence k-? model. By focusing on investigating the influences of flow and buoyancy on pollutants, it is found that while carrying by the river flow downstream: i) a heavier pollutant released from the cross-sectional side position, forms transverse oscillation between two banks with decreased amplitude, i.e. forms kind of helical flow pattern along the straight part of channel bed; ii) a heavier pollutant released from the cross-sectional middle position, forms collapse oscillation in the middle of the straight channel part with reduced amplitude; iii) in the downstream sinuous channel, heavier pollutant is of higher concentration on the outer side of channel bends; iv) a light pollutant released from the cross-sectional side position, slips partly to the other side of the river, resulting in higher concentrations on two sides of the channel top; v) a light pollutant released from the cross-sectional middle position, splits into two parts symmetrically along two sides of the channel top; vi) in the downstream sinuous channel, light pollutant presents higher concentration on the inner side of channel bends. These findings may assist in cost-effective scientific countermeasures to be taken for accidental or planned pollutant releases into a river.     

Occurrence and distribution of pharmaceuticals in surface water, suspended solids and sediments of the Ebro river basin, Spain

The occurrence of 43 pharmaceuticals belonging to predominant therapeutic classes and their distribution in surface water, suspended solids and sediments has been investigated in the Ebro river basin in the Northeast of Spain. WWTP effluents were found to be a main source of contamination and the spatial distribution was affected by the river flow at the sampling point and corresponding dilution factor, resulting in higher concentrations and higher loads in small tributary rivers than in the Ebro river. The study showed that some compounds are preferentially found bound to suspended solids and not detected in river water. Generally, compounds with basic characteristics (pKa > 7) showed higher tendency to bind to suspended solids. The sediment samples generally presented lower concentrations than suspended solids.     

Future climate scenarios and rainfall--runoff modelling in the Upper Gallego catchment (Spain)

Global climate change may have large impacts on water supplies, drought or flood frequencies and magnitudes in local and regional hydrologic systems. Water authorities therefore rely on computer models for quantitative impact prediction. In this study we present kernel-based learning machine river flow models for the Upper Gallego catchment of the Ebro basin. Different learning machines were calibrated using daily gauge data. The models posed two major challenges: (1) estimation of the rainfall-runoff transfer function from the available time series is complicated by anthropogenic regulation and mountainous terrain and (2) the river flow model is weak when only climate data are used, but additional antecedent flow data seemed to lead to delayed peak flow estimation. These types of models, together with the presented downscaled climate scenarios, can be used for climate change impact assessment in the Gallego, which is important for the future management of the system.     

Measuring the flux at the interface of coal-tar impacted sediment and river water near a former MGP site

Determining water movement through contaminated sediment is critical for characterizing transport of chemicals from the sediment to the overlying water. Field studies to characterize the water flow across the sediment-water interface within a river adjacent to a former manufactured gas plant site were conducted. For this purpose, a new design of an interfacial flow meter was developed and tested. The in situ components of the system consisted of: a cylinder with an interfacial area of 2342 cm2; a dome attached to the cylinder; and a flow tube that allows water to flow from inside the dome to the river at the rate equal to the specific discharge across the sediment-water boundary. A 'heat-pulse' method was used to measure flow by heating the center of the flow tube for a brief time period and measuring the temperature profile within the tube over time. The system was calibrated to measure volumetric flux in the range 1.5-4.0 cm d(-1), however using a flow-addition method, the measurement of lower velocities also was accomplished, and calibration at higher fluxes is possible. From the groundwater flow at the interface of the coal-tar impacted sediment and information on the sediment pore water concentrations of several PAHs (poly-cyclic aromatic hydrocarbons), the mass flux of these PAHs to the river were estimated. Information on PAH mass flux at the sediment-water interface is useful for site assessment, including the evaluation of remediation alternatives and longer term site characterization.     

An application of the theory of kinematics of mixing to the study of tropospheric dispersion

The most common technique used for numerical simulations of tracer mixing is that of the numerical solution of the advection–diffusion equation with the unresolved fluxes parameterized using the similarity theory. Despite correct predictions of the overall directions of transport, models based on a numerical solution of the advection–diffusion equation lack sufficient accuracy to correctly reproduce the coupling of mixing with small scale processes which are sensitive to the microstructure of the tracer distribution. The objective of this paper is to revisit the basic formalism employed in numerical models used to investigate atmospheric tracers. The main mathematical method proposed here is the theory of kinematics of mixing which could be applied effectively for simulations of atmospheric transport processes. At the beginning of the paper, we introduce simple mathematical transformations in order to demonstrate how complex topological structures are created by mixing processes. These idealistic flow systems are essential to explain transport properties of much more complex three-dimensional geophysical flows. An example of the application of the kinematics of mixing to the analysis of tracer transport on a planetary scale is presented in the following sections. The complex filamentary structures simulated in the numerical experiment are evaluated using some commonly applied statistical measures in order to compare the results with the data published in the literature. The results of the experiment are also analysed with the help of simple conceptual models of fluid filaments. The microstructure of the tracer distribution introduced in the paper is essential to increase our understanding of atmospheric transport and to develop more realistic parameterizations of small-scale mixing. The presented results could also be used to improve calculations of the coupling between microphysical processes and tracer mixing.     

Modelling radionuclide transport for time varying flow in a channel network

Water flowrates and flow directions may change over time in the subsurface for a number of reasons. In fractured rocks flow takes place in channels within fractures. Solutes are carried by the advective flow. In addition, solutes may diffuse in and out of stagnant waters in the rock matrix and other stagnant water regions. Sorbing species may sorb on fracture surfaces and on the micropore surfaces in the rock matrix. We present a method by which solute particles can be traced in flowing water undergoing changes in flowrate and direction in a complex channel network where the solutes can also interact with the rock by diffusion in the rock matrix. The novelty of this paper is handling of diffusion in the rock matrix under transient flow conditions. The diffusive processes are stochastic and it is not possible to follow a particle deterministically. The method therefore utilises the properties of a probability distribution function for a tracer moving in a fracture where matrix diffusion is active. The method is incorporated in a three dimensional channel network model. Particle tracking is used to trace out a multitude of flowpaths, each of which consists of a large number of channels within fractures. Along each channel the aperture and velocity as well as the matrix sorption properties can vary. An efficient method is presented whereby a particle can be followed along the variable property flowpath. For stationary flow conditions and a network of channels with advective flow and matrix diffusion, a simple analytical solution for the residence time distribution along each pathway can be used. Only two parameter groups need to be integrated along each path. For transient flow conditions, a time stepping procedure that incorporates a stochastic Monte-Carlo like method to follow the particles along the paths when flow conditions change is used. The method is fast and an example is used for illustrative purposes. It is exemplified by a case where land rises due to glacial rebound. It is shown that the effects of changing flowrates and directions can be considerable and that the diffusive migration in the matrix can have a dominating effect on the results.     

An interpretation of potential scale dependence of the effective matrix diffusion coefficient

Matrix diffusion is an important process for solute transport in fractured rock, and the matrix diffusion coefficient is a key parameter for describing this process. Previous studies have indicated that the effective matrix diffusion coefficient values, obtained from a large number of field tracer tests, are enhanced in comparison with local values and may increase with test scale. In this study, we have performed numerical experiments to investigate potential mechanisms behind possible scale-dependent behavior. The focus of the experiments is on solute transport in flow paths having geometries consistent with percolation theories and characterized by multiple local flow loops formed mainly by small-scale fractures. The water velocity distribution through a flow path was determined using discrete fracture network flow simulations, and solute transport was calculated using a previously derived impulse-response function and a particle-tracking scheme. Values for effective (or up-scaled) transport parameters were obtained by matching breakthrough curves from numerical experiments with an analytical solution for solute transport along a single fracture. Results indicate that a combination of local flow loops and the associated matrix diffusion process, together with scaling properties in flow path geometry, seems to be an important mechanism causing the observed scale dependence of the effective matrix diffusion coefficient (at a range of scales).     

Use of tritium to predict soluble pollutants transport in Ebro River waters (Spain)

The Ebro River, in Northeast Spain, discharges into the Mediterranean Sea after flowing through several large cities and agricultural, mining and industrial areas. The Ascó nuclear power plant (NPP) is located in its lower section and comprises two pressurised water reactor units, from which low-level liquid radioactive waste is released to river waters under authority control. Tritium routinely released by the NPP was used as a radiotracer to determine the longitudinal dispersion coefficient and velocity of the river waters. Several field experiments, in co-ordination with the NPP, were carried out during 1991 and 1992. During each field experiment, the flow rate was kept constant by dams located upstream from the NPP. After each tritium release, water was sampled downstream at periodic intervals over several hours and tritium was measured with a low-background liquid scintillation counter. Velocity and dispersion coefficient were determined in river waters for several river discharges using an analytical, box-type and numerical approach to solve the one-dimensional advection-diffusion equation. The set of calibrated parameters was used to predict the displacement and dispersion of soluble pollutants in river waters. Velocity was determined as a function of river discharge and river slope, and dispersion coefficient was determined as a function of distance. Finally, sensitivity of the model predictions was studied and uncertainties of the fitted parameters were estimated.