Ground water contaminant transport by nondivergence-free, unsteady and nonstationary velocity fields

Pore flow velocity is assumed to be a nondivergence-free, unsteady, and nonstationary random function of space and time for ground water contaminant transport in a heterogeneous medium. The laboratory-scale stochastic contaminant transport equation is up scaled to field scale by taking the ensemble average of the equation by using the cumulant expansion method. A new velocity correction, which is a function of mean pore flow velocity divergence, is obtained due to strict second order cumulant expansion (without omitting any term after the expansion). The field scale transport equations under the divergence-free pore flow velocity field assumption are also derived by simplifying the nondivergence-free field scale equation. The significance of the new velocity correction term is investigated on a two dimensional transport problem driven by a density dependent flow.     

Kriging Regionalized Positive Variables Revisited: Sample Space and Scale Considerations

Frequently, regionalized positive variables are treated by preliminarily applying a logarithm, and kriging estimates are back-transformed using classical formulae for the expectation of a lognormal random variable. This practice has several problems (lack of robustness, non-optimal confidence intervals, etc.), particularly when estimating block averages. Therefore, many practitioners take exponentials of the kriging estimates, although the final estimations are deemed as non-optimal. Another approach arises when the nature of the sample space and the scale of the data are considered. Since these concepts can be suitably captured by an Euclidean space structure, we may define an optimal kriging estimator for positive variables, with all properties analogous to those of linear geostatistical techniques, even for the estimation of block averages. In this particular case, no assumption on preservation of lognormality is needed. From a practical point of view, the proposed method coincides with the median estimator and offers theoretical ground to this extended practice. Thus, existing software and routines remain fully applicable.     

起伏地表条件下地震波场数值模拟有限积分变换有限差分方法

针对起伏地表条件下的地震波数值模拟问题,提出将起伏的地表映射到一个规则的长方形网格坐标系中,并在此基础上推导出变换域中的波动方程。根据导出来的波动方程应用有限余弦变换有限差分方法进行地震波场的数值模拟,以解决起伏地表条件下的数值模拟。     

The solution of the Korn–Lichtenstein equations of conformal mapping: the direct generation of ellipsoidal Gauß–Krüger conformal coordinates or the Transverse Mercator Projection

The differential equations which generate a general conformal mapping of a two-dimensional Riemann manifold found by Korn and Lichtenstein are reviewed. The Korn–Lichtenstein equations subject to the integrability conditions of type vectorial Laplace–Beltrami equations are solved for the geometry of an ellipsoid of revolution (International Reference Ellipsoid), specifically in the function space of bivariate polynomials in terms of surface normal ellipsoidal longitude and ellipsoidal latitude. The related coefficient constraints are collected in two corollaries. We present the constraints to the general solution of the Korn–Lichtenstein equations which directly generates Gau?–Krüger conformal coordinates as well as the Universal Transverse Mercator Projection (UTM) avoiding any intermediate isometric coordinate representation. Namely, the equidistant mapping of a meridian of reference generates the constraints in question. Finally, the detailed computation of the solution is given in terms of bivariate polynomials up to degree five with coefficients listed in closed form. Received: 3 June 1997 / Accepted: 17 November 1997     

方向观测后方交会的点位精度

以虚拟误差方程,经严密推证,给出了以方向作为独立观测值的后方交会点位精度估算公式。以大桥施工控制网为例,介绍了３个方向与４个方向后方交会的精度状况及其应用。     

线性方程组迭代解的随机模型测试研究

本文讨论大型线性方程组迭代解的随机模型测试评价问题。给出了常用迭代解法ＣＧ、ＬＳＱＲ、ＳＩＲＴ、ＳＡＲＴ、ＳＡＳＩＲＴ等的测试结果。结果表明：（１）方程组系数矩阵的特性（条件数）及解结构都对解精度有重要影响。解模型越粗糙,解的精度越低。（２）各种求解算法都有一定的平滑效应,同时各种算法也都会产生误差大于２００％的奇异解,奇异解元素数一般约占１０％。（３）数据的拟合残差一般不能真实反映解的精度。（４）对含误差数据的求解问题,较好的求解算法是ＤＬＳＱＲ与ＳＡＳＩＲＴ。     

考虑权值非负约束的泛克立格算法

泛克立格法是地质统计学的一种重要估值方法,当区域化变量在空间变异几何域内非平稳时常用泛克立格法来估值。然而,用泛克立格方程组计算估计权值时由于未对权的符号作任何限制,从而使得计算出的权值经常出现负权现象,而负权的存在有许多弊端,所以在许多应用中有必要对权值作非负要求。本文基于线性规划方法提出了一种考虑权值非负约束的泛克立格算法,该算法既考虑到了权值的非负约束条件,又利用了线性规划方法求解简便快捷的优点。     

At‐a‐station hydraulic geometry relations, 1: theoretical development

This paper, the first of two, hypothesizes that: (1) the temporal variation of stream power of a river channel at a given station with varying discharge is accomplished by the temporal variation in channel form (flow depth and channel width) and hydraulic variables, including energy slope, flow velocity and friction; (2) the change in stream power is distributed among the changes in flow depth, channel width, flow velocity, slope, and friction, depending on the boundary conditions that the channels has to satisfy. The second hypothesis is a result of the principle of maximum entropy and the theory of minimum energy dissipation or its simplified minimum stream power. These two hypotheses lead to families of at‐a‐station hydraulic geometry relations. The conditions under which these families of relations can occur in the field are discussed. Copyright © 2007 John Wiley & Sons, Ltd.     

A mass balance approach to estimate the dilution and removal of the pollutants in stream water polluted by acid mine drainage

 A few simple mass balance equations were developed to simultaneously estimate how much the pollutants from acid mine drainage (AMD) in stream water are diluted and removed during their migration. The application of the equations requires knowledge of the variations in the concentrations of the dissolved pollutants and the stoichiometry of the precipitation reaction of the pollutants when none of the pollutant shows a conservative behavior along the stream path. The calculation should be restricted to the pollutants showing much higher concentrations in the polluted main stream water than in the combining or diluting water of the same target area. The mass balance equations were applied to estimate the dilution factor and precipitation fractions of pollutants in Imgok Creek such as Fe, SO₄ and Al from the AMD of Yeongdong mine. The results show that the estimation, especially for SO₄ and Al, significantly depends on the kinds of the precipitates. When FeOHSO₄ and AlOHSO₄ are assumed to precipitate, the maximum removal fractions of SO₄ and Al by precipitation are respectively 34% and 46% of the original input, which is much higher than the values estimated when SO₄ is considered to be perfectly conservative. It indicates that the stoichiometry of precipitation reaction is very important in the interpretation of the pollutant dilution and migration and assessment of environmental impacts of AMD. The applicability of the mass balance equations may still need to be verified. However, examining the calculated dilution factor and precipitation fractions with the equations can provide invaluable information on not only the behavior but also unexpected input of the pollutants in the stream water polluted by AMD and other point sources. Received: 12 November 1997 · Accepted: 30 March 1998     

A fuzzy dynamic flood routing model for natural channels

A fuzzy dynamic flood routing model (FDFRM) for natural channels is presented, wherein the flood wave can be approximated to a monoclinal wave. This study is based on modification of an earlier published work by the same authors, where the nature of the wave was of gravity type. Momentum equation of the dynamic wave model is replaced by a fuzzy rule based model, while retaining the continuity equation in its complete form. Hence, the FDFRM gets rid of the assumptions associated with the momentum equation. Also, it overcomes the necessity of calculating friction slope (S_f) in flood routing and hence the associated uncertainties are eliminated. The fuzzy rule based model is developed on an equation for wave velocity, which is obtained in terms of discontinuities in the gradient of flow parameters. The channel reach is divided into a number of approximately uniform sub‐reaches. Training set required for development of the fuzzy rule based model for each sub‐reach is obtained from discharge‐area relationship at its mean section. For highly heterogeneous sub‐reaches, optimized fuzzy rule based models are obtained by means of a neuro‐fuzzy algorithm. For demonstration, the FDFRM is applied to flood routing problems in a fictitious channel with single uniform reach, in a fictitious channel with two uniform sub‐reaches and also in a natural channel with a number of approximately uniform sub‐reaches. It is observed that in cases of the fictitious channels, the FDFRM outputs match well with those of an implicit numerical model (INM), which solves the dynamic wave equations using an implicit numerical scheme. For the natural channel, the FDFRM outputs are comparable to those of the HEC‐RAS model. Copyright © 2009 John Wiley & Sons, Ltd.