微曝气垂直流湿地处理城郊低浓度生活污水模拟实验

本文提出微曝气垂直流湿地的概念,并通过柱状模型实验模拟研究了其对云南滇池流域城郊低浓度生活污水的处理效果.测定了模型的氧利用效率,最高为14.04%.在高水力负荷15.29m3/(m2·d)、气水比1:2的条件下,BOD5和NI4-N去除率分别为74.4%和63.5%,出水平均含量分别为4.96mg/L和272mg/L.结果说明,微曝气垂直流湿地工艺能满足低浓度生活污水高水力负荷的处理要求,微曝气耗能折价O.04-0.05元/m3.     

Density-dependent dispersion in heterogeneous porous media Part I: A numerical study

In this paper, we describe carefully conducted numerical experiments, in which a dense salt solution vertically displaces fresh water in a stable manner. The two-dimensional porous media are weakly heterogeneous at a small scale. The purpose of these simulations, conducted for a range of density differences, is to obtain accurate concentration profiles that can be used to validate nonlinear models for high-concentration-gradient dispersion. In this part we focus on convergence of the computations, in numerical and statistical sense, to ensure that the uncertainty in the results is small enough.Concentration variances are computed, which give estimates of the uncertainty in local concentration values. These local variations decrease with increasing density contrast. For tracer transport, obtained longitudinal dispersivities are in accordance with analytical findings. In the case of high-density contrasts, stabilizing gravity forces counteract the growth of dispersive fingers, decreasing the effective width of the transition zone. For small log-permeability variances, the decrease of the apparent dispersivity that is found is in agreement with laboratory results for homogeneous columns.     

Two-dimensional power series solution for non-axisymmetrical transport in a radially convergent tracer test with scale-dependent dispersion

It has been known for many years that dispersivity increases with solute travel distance in a subsurface environment. The increase of dispersivity with solute travel distance results from the significant variation of hydraulic properties of heterogeneous media and was identified in the literature as scale-dependent dispersion. This study presents an analytical solution for describing two-dimensional non-axisymmetrical solute transport in a radially convergent flow tracer test with scale-dependent dispersion. The power series technique coupling with the Laplace and finite Fourier cosine transform has been applied to yield the analytical solution to the two-dimensional, scale-dependent advection–dispersion equation in cylindrical coordinates with variable-dependent coefficients. Comparison between the breakthrough curves of the power series solution and the numerical solutions shows excellent agreement at different observation points and for various ranges of scale-related transport parameters of interest. The developed power series solution facilitates fast prediction of the breakthrough curves at any observation point.     

Multi-resolution adaptive modeling of groundwater flow and transport problems

Many groundwater flow and transport problems, especially those with sharp fronts, narrow transition zones, layers and fingers, require extensive computational resources. In this paper, we present a novel multi-resolution adaptive Fup approach to solve the above mentioned problems. Our numerical procedure is the Adaptive Fup Collocation Method (AFCM), based on Fup basis functions and designed through a method of lines (MOL). Fup basis functions are localized and infinitely differentiable functions with compact support and are related to more standard choices such as splines or wavelets. This method enables the adaptive multi-resolution approach to solve problems with different spatial and temporal scales with a desired level of accuracy using the entire family of Fup basis functions. In addition, the utilized collocation algorithm enables the mesh free approach with consistent velocity approximation and flux continuity due to properties of the Fup basis functions. The introduced numerical procedure was tested and verified by a few characteristic groundwater flow and transport problems, the Buckley–Leverett multiphase flow problem, the 1-D vertical density driven problem and the standard 2-D seawater intrusion benchmark–Henry problem. The results demonstrate that the method is robust and efficient particularly when describing sharp fronts and narrow transition zones changing in space and time.     

Analytical models of contaminant transport in coastal aquifers

The Henry formulation, which couples subsurface flow and salt transport via a variable-density flow formulation, can be used to evaluate the extent of sea water intrusion into coastal aquifers. The coupling gives rise to nontrivial flow patterns that are very different from those observed in inland aquifers. We investigate the influence of these flow patterns on the transport of conservative contaminants in a coastal aquifer. The flow is characterized by two dimensionless parameters: the Péclet number, which compares the relative effects of advective and dispersive transport mechanisms, and a coupling parameter, which describes the importance of the salt water boundary on the flow. We focus our attention on two regimes – low and intermediate Péclet number flows. Two transport scenarios are solved analytically by means of a perturbation analysis. The first, a natural attenuation scenario, describes the flushing of a contaminant from a coastal aquifer by clean fresh water, while the second, a contaminant spill scenario, considers an isolated point source.     

Coarsening of three-dimensional structured and unstructured grids for subsurface flow

We present a generic, semi-automated algorithm for generating non-uniform coarse grids for modeling subsurface flow. The method is applicable to arbitrary grids and does not impose smoothness constraints on the coarse grid. One therefore avoids conventional smoothing procedures that are commonly used to ensure that the grids obtained with standard coarsening procedures are not too rough. The coarsening algorithm is very simple and essentially involves only two parameters that specify the level of coarsening. Consequently the algorithm allows the user to specify the simulation grid dynamically to fit available computer resources, and, e.g., use the original geomodel as input for flow simulations. This is of great importance since coarse grid-generation is normally the most time-consuming part of an upscaling phase, and therefore the main obstacle that has prevented simulation workflows with user-defined resolution. We apply the coarsening algorithm to a series of two-phase flow problems on both structured (Cartesian) and unstructured grids. The numerical results demonstrate that one consistently obtains significantly more accurate results using the proposed non-uniform coarsening strategy than with corresponding uniform coarse grids with roughly the same number of cells.     

Flow and deposition of pyroclastic granular flows: A type example from the 1975 Ngauruhoe eruption,New Zealand

Small-volume pyroclastic density currents (PDCs) are generated frequently during explosive eruptions with little warning. Assessing their hazard requires a physical understanding of their transport and sedimentation processes which is best achieved by the testing of experimental and numerical models of geophysical mass flows against natural flows and/or deposits. To this end we report on one of the most detailed sedimentological studies ever carried out on a series of pristine small-volume PDC deposits from the 1975 eruption of Ngauruhoe volcano, whose emplacement were also witnessed during eruption. Using high-resolution GPS surveys, a series of lateral excavations across the deposits, and bulk sedimentological analysis we constrained the geomorphology, internal structure and texture of the deposits with respect to laterally varying modes of deposition.     

Magmatic structures in the Krkonoše–Jizera Plutonic Complex,Bohemian Massif: evidence for localized multiphase flow and small-scale thermal–mechanical instabilities in a granitic magma chamber

The present paper examines magmatic structures in the Jizera and Liberec granites of the Krkonoše–Jizera Plutonic Complex, Bohemian Massif. The magmatic structures are here interpreted to preserve direct field evidence for highly localized magma flow and other processes in crystal-rich mushes, and to capture the evolution of physical processes in an ancient granitic magma chamber. We propose that after chamber-wide mixing and hybridization, as suggested by recent petrological studies, laminar magma flow became highly localized to weaker channel-like domains within the higher-strength crystal framework. Mafic schlieren formed at flow rims, and their formation presumably involved gravitational settling and velocity gradient flow sorting coupled with interstitial melt escape. Local thermal or compositional convection may have resulted in the formation of vertical schlieren tubes and ladder dikes whereas subhorizontal tubes or channels formed during flow driven by lateral gradients in magma pressure. After the cessation or deceleration of channel flow, gravity-driven processes (settling of crystals and enclaves, gravitational differentiation, development of downward dripping instabilities), accompanied by compaction, filter pressing and melt segregation, dominated in the crystal mush within the flow channels. Subsequently, magmatic folds developed in schlieren layers and the magma chamber recorded complex, late magmatic strains at high magma crystallinities. Late-stage magma pulsing into localized submagmatic cracks represents the latest events of magmatic history of the chamber prior to its final crystallization. We emphasize that the most favorable environments for the formation and preservation of magmatic structures, such as those hosted in the Jizera and Liberec granites, are slowly cooling crystal-rich mushes. Therefore, where preserved in plutons, these structures may lend strong support for a “mush model” of magmatic systems.     

Magnetic petrology of the 1991–1995 dacite lava of Unzen volcano,Japan: Degree of oxidation and implications for the growth of lava domes

To understand the oxidation state and process of oxidation of lava domes, we carried out magnetic petrological analyses of lava samples obtained from domes and block-and-ash-flow deposits associated with the 1991–1995 eruption of Unzen volcano, Japan. As a result, we recognize three different types of magnetic petrology, each related to deuteric high-temperature oxidation during initial cooling. Type A oxides are characterized by homogenous titanomagnetite and titanohematite, indicating a low oxidation state and high titanomagnetite concentrations. Type B oxides are weakly exsolved and contain titanohematite laths and rutile lenses, indicating a higher oxidation state. Type C oxides, which represent the highest oxidation state, are completely exsolved and composed of Ti-poor titanomagnetite, titanohematite, rutile, and pseudobrookite, indicating high hematite concentrations. Some grains in Types A and B show indications of reduction, which was related to interaction with volcanic gases subsequent to high-temperature oxidation. In terms of geological occurrence, the oxidation processes probably differed for endogenous and exogenous domes. Endogenous dome lavas are oxidized concentrically and are classified into the three types according to their location within the dome: samples from the surface are strongly oxidized and classified as Type C, while the inner part is unoxidized and classified as Type A. Exogenous dome lavas are unoxidized and assigned to Type A. Some samples show signs of reduction, which may have occurred around fumaroles. We propose that location within the dome and the process of dome growth are the factors that control oxidation.     

Nonlinear determinism in river flow: prediction as a possible indicator

Whether or not river flow exhibits nonlinear determinism remains an unresolved question. While studies on the use of nonlinear deterministic methods for modeling and prediction of river flow series are on the rise and the outcomes are encouraging, suspicions and criticisms of such studies continue to exist as well. An important reason for this situation is that the correlation dimension method, used as a nonlinear determinism identification tool in most of those studies, may possess certain limitations when applied to real river flow series, which are always finite and often short and also contaminated with noise (e.g. measurement error). In view of this, the present study addresses the issue of nonlinear determinism in river flow series using prediction as a possible indicator. This is done by (1) reviewing studies that have employed nonlinear deterministic methods (coupling phase‐space reconstruction and local approximation techniques) for river flow predictions and (2) identifying nonlinear determinism (or linear stochasticity) based on the level of prediction accuracy in general, and on the prediction accuracy against the phase‐space reconstruction parameters in particular (termed as the ‘inverse approach’). The results not only provide possible indications to the presence of nonlinear determinism in the river flow series studied, but also support, both qualitatively and quantitatively, the low correlation dimensions reported for such. Therefore, nonlinear deterministic methods are a viable complement to linear stochastic ones for studying river flow dynamics, if sufficient caution is exercised in their applications and in interpreting the outcomes. Copyright © 2006 John Wiley & Sons, Ltd.