Sensitivity analysis of 2D steady-state shallow water flow. Application to free surface flow model calibration

This paper presents the analytical properties of the solutions of the sensitivity equations for steady-state, two-dimensional shallow water flow. These analytical properties are used to provide guidelines for model calibration and validation. The sensitivity of the water depth/level and that of the longitudinal unit discharge are shown to contain redundant information. Under subcritical conditions, the sensitivities of the flow variables are shown to obey an anisotropic elliptic equation. The main directions of the contour lines for water depth and the longitudinal unit discharge sensitivity are parallel and perpendicular to the flow, while they are diagonal to the flow for the transverse unit discharge sensitivity. Moreover, the sensitivity for all three variables extends farther in the transverse direction than in the longitudinal direction, the anisotropy ratio being a function of the sole Froude number. For supercritical flow, the sensitivity obeys an anisotropic hyperbolic equation. These findings are confirmed by application examples on idealized and real-world simulations. The sensitivities to the geometry, friction coefficient or model boundary conditions are shown to behave in different ways, thus providing different types of information for model calibration and validation.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

An effective scale-dependent dispersivity deduced from a purely convective flow field

Abstract

In the case of straight flow but with hydraulic conductivity varying in a transverse direction, the distribution of hydraulic conductivity has been determined for which the breakthrough curve due to convection only will have the same analytical form as the onedimensional convection/dispersion equation solution at the outlet end of a porous medium. That distribution is found exactly and it is very similar to the lognormal distribution. This result is significant since field evidence indicates that the logarithm of hydraulic conductivity is normally distributed. For the case considered, a simple relation between dispersivity and the coefficient of variation of hydraulic conductivity is found. One can thus determine very simply dispersivity in terms of the parameters of the distribution of hydraulic conductivity. This is particularly useful to estimate dispersivity in various cells of finite difference or finite element models when the distribution of hydraulic conductivity is not stationary, i.e. varies in space.     

Eulerian spatial moments for solute transport in three-dimensional heterogeneous,dual-permeability media

A Eulerian analytical method is developed for nonreactive solute transport in heterogeneous, dual-permeability media where the hydraulic conductivities in fracture and matrix domains are both assumed to be stochastic processes. The analytical solution for the mean concentration is given explicitly in Fourier and Laplace transforms. Instead of using the fast fourier transform method to numerically invert the solution to real space (Hu et al., 2002), we apply the general relationship between spatial moments and concentration (Naff, 1990; Hu et al., 1997) to obtain the analytical solutions for the spatial moments up to the second for a pulse input of the solute. Owing to its accuracy and efficiency, the analytical method can be used to check the semi-analytical and Monte Carlo numerical methods before they are applied to more complicated studies. The analytical method can be also used during screening studies to identify the most significant transport parameters for further analysis. In this study, the analytical results have been compared with those obtained from the semi-analytical method (Hu et al., 2002) and the comparison shows that the semi-analytical method is robust. It is clearly shown from the analytical solution that the three factors, local dispersion, conductivity variation in each domain and velocity convection flow difference in the two domains, play different roles on the solute plume spreading in longitudinal and transverse directions. The calculation results also indicate that when the log-conductivity variance in matrix is 10 times less than its counterpart in fractures, it will hardly influence the solute transport, whether the conductivity field is matrix is treated as a homogeneous or random field.     

Analytical studies on transient groundwater flow induced by land reclamation using different fill materials

Land reclamation may have a significant influence on groundwater regimes. Analytical solutions have been developed in the past to study the impact of land reclamation on a steady‐state groundwater flow and transient flow in fill materials, assuming that the reclamation site consists of a single zone of uniform hydraulic parameters. In this paper, we derive analytical solutions to describe the transient water table change in response to multi‐stage land reclamation where the fill material is uniform in each stage but the hydraulic conductivity of the fill material varies from stage to stage. By introducing the method of separation of variables, we develop a transient analytical solution to study the impact of land reclamation consisting of fill material with different hydraulic properties on groundwater dynamics. The results show that the water table first increases significantly into the reclaimed zone following the fill material deposition, and then the increase gradually propagates into the original aquifer. The change of water table in the original aquifer mainly depends on the value of hydraulic conductivity of the fill materials. Examples in this paper illustrate how the aquifer system experiences a long time unsteady‐state flow as a result of the reclamation, and it takes at least tens of years for the system to approach a new equilibrium. It is suggested that for a large‐scale reclamation project, the response of the groundwater regime to reclamation should be carefully studied. Copyright © 2013 John Wiley & Sons, Ltd.     

Shear flow instabilities in a rotating system

Abstract

The stability of a plane parallel shear flow with the profile U(z) = tanh z is considered in a rotating system with the axis of rotation in the z-direction. The establishment of the basic flow requires a baroclinic state, but baroclinic effects are suppressed in the stability analysis by assuming a limit of high thermal conductivity. It is shown that the strongest growing disturbance changes from a purely transverse form in the limit of vanishing rotation rate to a nearly longitudinal form as the angular velocity of rotation increases. An analytical solution of the stability equation is obtained for vanishing growth rates of the transverse form of the instability. But, in general, the solution of the problem requires numerical integrations which demonstrate that the preferred direction of the wave vector of the instability is towards the left of the direction of the mean flow.     

Using data assimilation method to calibrate a heterogeneous conductivity field conditioning on transient flow test data

A data assimilation method is developed to calibrate a heterogeneous hydraulic conductivity field conditioning on transient pumping test data. The ensemble Kalman filter (EnKF) approach is used to update model parameters such as hydraulic conductivity and model variables such as hydraulic head using available data. A synthetical two-dimensional flow case is used to assess the capability of the EnKF method to calibrate a heterogeneous conductivity field by assimilating transient flow data from observation wells under different hydraulic boundary conditions. The study results indicate that the EnKF method will significantly improve the estimation of the hydraulic conductivity field by assimilating continuous hydraulic head measurements and the hydraulic boundary condition will significantly affect the simulation results. For our cases, after a few data assimilation steps, the assimilated conductivity field with four Neumann boundaries matches the real field well while the assimilated conductivity field with mixed Dirichlet and Neumann boundaries does not. We found in our cases that the ensemble size should be 300 or larger for the numerical simulation. The number and the locations of the observation wells will significantly affect the hydraulic conductivity field calibration.     

Numerical modelling of stratified tidal flow over a fjord sill

Field observations of tidally driven stratified flow in the sill area of Knight Inlet (British Columbia) revealed a very complicated structure, which includes solitary waves, upstream bifurcation, hydraulic jump and mixing processes. Recent observations suggest that the flow instabilities on the plunging pycnocline at the lee side of the sill may contribute to solitary wave generation through a subharmonic interaction. The present study reports on a series of numerical experiments of stratified tidal flow in Knight Inlet performed with the help of a fine resolution fully non-linear non-hydrostatic numerical model. The model reproduces all important stages of the baroclinic tidal dynamics observed in Knight Inlet. Results demonstrate that solitary waves are generated apart from the area of hydrodynamic instability. Accelerating tidal flux forms a baroclinic hydraulic jump just above the top of the sill, whereas the bifurcations and zones of shear instabilities are formed downstream of the sill. The first baroclinic mode having the largest velocity escapes from the generation area and propagates upstream, disintegrating further into a packet of solitary waves reviling the classical “non-subharmonic” mechanism of generation. The remaining part of the disturbance (slow baroclinic modes) is arrested by tidal flow and carried away to the lee side of the obstacle, where shear instability, billows and mixing processes are developed. Some sensitivity runs were performed for different value of tidal velocity.     

An efficient tool for accelerating the numerical solution of the stochastic subsurface flow problem using neural networks

 An efficient numerical solution for the two-dimensional groundwater flow problem using artificial neural networks (ANNs) is presented. Under stationary velocity conditions with unidirectional mean flow, the conductivity realizations and the head gradients, obtained by a traditional finite difference solution to the flow equation, are given as input-output pairs to train a neural network. The ANN is trained successfully and a certain level of recognition of the relationship between input conductivity patterns and output head gradients is achieved. The trained network produced velocity realizations that are physically plausible without solving the flow equation for each of the conductivity realizations. This is achieved in a small fraction of the time necessary for solving the flow equations. The prediction accuracy of the ANN reaches 97.5% for the longitudinal head gradient and 94.7% for the transverse gradient. Head-gradient and velocity statistics in terms of the first two moments are obtained with a very high accuracy. The cross covariances between head gradients and the fluctuating log-conductivity (log-K) and between velocity and log-K obtained with the ANN approach match very closely those obtained by a traditional numerical solution. The same is true for the velocity components auto-covariances. The results are also extended to transport simulations with very good accuracy. Spatial moments (up to the fourth) of mean-concentration plumes obtained using ANNs are in very good agreement with the traditional Monte Carlo simulations. Furthermore, the concentration second moment (concentration variance) is very close between the two approaches. Considering the fact that higher moments of concentration need more computational effort in numerical simulations, the advantage of the presented approach in saving long computational times is evident. Another advantage of the ANNs approach is the ability to generalize a trained network to conductivity distributions different from those used in training. However, the accuracy of the approach in cases with higher conductivity variances is being investigated.     

Influence of transient flow on contaminant biodegradation

The rate of biodegradation in contaminated aquifers depends to a large extent on dispersive mixing processes that are now generally accepted to result from spatial variations in the velocity field. It has been shown, however, that transient flow fields can also contribute to dispersive mixing. The influence of transient flow on biodegrading contaminants is particularly important since it can enhance mixing with electron acceptors, further promoting the reactive process. Using numerical simulations, the effect of transient flow on the behavior of a biodegradable contaminant is evaluated here both with respect to the development of apparently large horizontal transverse dispersion and also with respect to enhanced mixing between the substrate (electron donor) and electron acceptor. The numerical model BIO3D, which solves for advective-dispersive transport coupled with Monod-type biodegradation of substrates in the presence of an electron acceptor, was used for the simulations. The model was applied in a two-dimensional plan view mode considering a single substrate. Transient flow fields were found to yield larger apparent transverse dispersion because the longitudinal dispersivity also acts transverse to the mean flow direction. In the reactive case, the transient flow field increases substrate-oxygen mixing, which in turn enhances the overall rate of biodegradation. The results suggest that in the case of moderate changes of flow directions, a steady-state flow field can be justified, thereby avoiding the higher computational costs of a fully transient simulation. The use of a higher transverse horizontal dispersivity in a steady flow field can, under these conditions, adequately forecast plume development.     

Flow in unsaturated random porous media,nonlinear numerical analysis and comparison to analytical stochastic models

This work presents a rigorous numerical validation of analytical stochastic models of steady state unsaturated flow in heterogeneous porous media. It also provides a crucial link between stochastic theory based on simplifying assumptions and empirical field and simulation evidence of variably saturated flow in actual or realistic hypothetical heterogeneous porous media. Statistical properties of unsaturated hydraulic conductivity, soil water tension, and soil water flux in heterogeneous soils are investigated through high resolution Monte Carlo simulations of a wide range of steady state flow problems in a quasi-unbounded domain. In agreement with assumptions in analytical stochastic models of unsaturated flow, hydraulic conductivity and soil water tension are found to be lognormally and normally distributed, respectively. In contrast, simulations indicate that in moderate to strong variable conductivity fields, longitudinal flux is highly skewed. Transverse flux distributions are leptokurtic. the moments of the probability distributions obtained from Monte Carlo simulations are compared to modified first-order analytical models. Under moderate to strong heterogeneous soil flux conditions (σ²_y≥1), analytical solutions overestimate variability in soil water tension by up to 40% as soil heterogeneity increases, and underestimate variability of both flux components by up to a factor 5. Theoretically predicted model (cross-)covariance agree well with the numerical sample (cross-)covarianaces. Statistical moments are shown to be consistent with observed physical characteristics of unsaturated flow in heterogeneous soils.©1998 Elsevier Science Limited. All rights reserved     

Two-region non-Darcian flow toward a well in a confined aquifer

We have derived an analytical solution for two-region flow toward a well in a confined aquifer based on a linearization method. The two-region flow includes Izbash non-Darcian flow near the well and Darcian flow in the rest of the aquifer. The wellbore storage is also considered. The type curves in the non-Darcian and Darcian flow domains are obtained by a numerical Laplace inversion method incorporated in MATLAB programs. We have compared our results with the one-region Darcian flow model (Theis). Our solutions agree with those of Sen [Sen Z. Type curves for two-region well flow. J Hydr Eng 1988;114(12):1461–84] which were obtained using the Boltzmann transform at late times for fully turbulent flow, while some difference has been found at early and moderate times. We have defined a dimensionless non-Darcian hydraulic conductivity term which is shown to be a key parameter for analyzing the two-region flow. A smaller dimensionless non-Darcian hydraulic conductivity results in a larger drawdown in the non-Darcian flow region at late times. However, the dimensionless non-Darcian hydraulic conductivity does not affect the slope of the dimensionless drawdown versus the logarithmic dimensionless time in the non-Darcian flow region at late times. The dimensionless non-Darcian hydraulic conductivity does not affect the late time drawdown in the Darcian flow region.     

水力压裂对速度场及微地震定位的影响

水力压裂是页岩气开发过程中的核心增产技术,微地震则广泛用于压裂分析、水驱前缘监测和储层描述.微地震反演过程中,用于反演的速度模型往往基于测井、地震或标定炮资料构建,忽略了压裂过程中裂缝及孔隙流体压力变化对地层速度的影响.本文首先基于物质守恒、渗流理论和断裂力学模拟三维水力压裂过程,得到地下裂缝发育特征和孔隙压力分布.继而根据Coates-Schoenberg方法和裂缝柔量参数计算裂缝和孔隙压力对速度场的影响,得到压裂过程中的实时速度模型.最后利用三维射线追踪方法正演微地震走时和方位信息,并采用常规微地震定位方法反演震源位置及进行误差分析.数值模拟结果表明,检波器空间分布影响定位精度,常规方法的定位误差随射线路径在压裂带中传播距离增加而变大,且不同压裂阶段的多点反演法与单点极化法精度相当.     

A closed form slug test theory for high permeability aquifers

We incorporate a linear estimate of casing friction into the analytical slug test theory of Springer and Gelhar (1991) for high permeability aquifers. The modified theory elucidates the influence of inertia and casing friction on consistent, closed form equations for the free surface, pressure, and velocity fluctuations for overdamped and underdamped conditions. A consistent, but small, correction for kinetic energy is included as well. A characteristic velocity linearizes the turbulent casing shear stress so that an analytical solution for attenuated, phase shifted pressure fluctuations fits a single parameter (damping frequency) to transducer data from any depth in the casing. Underdamped slug tests of 0.3, 0.6, and 1 m amplitudes at five transducer depths in a 5.1 cm diameter PVC well 21 m deep in the Plymouth-Carver Aquifer yield a consistent hydraulic conductivity of 1.5 x 10(-3) m/s. The Springer and Gelhar (1991) model underestimates the hydraulic conductivity for these tests by as much as 25% by improperly ascribing smooth turbulent casing friction to the aquifer. The match point normalization of Butler (1998) agrees with our fitted hydraulic conductivity, however, when friction is included in the damping frequency. Zurbuchen et al. (2002) use a numerical model to establish a similar sensitivity of hydraulic conductivity to nonlinear casing friction.     

A geostatistically based inverse model for three-dimensional variably saturated flow

We present a geostatistically based inverse model for characterizing heterogeneity in parameters of unsaturated hydraulic conductivity for three-dimensional flow. Pressure and moisture content are related to perturbations in hydraulic parameters through cross-covariances, which are calculated to first-order. Sensitivities needed for covariance calculations are derived using the adjoint state sensitivity method. Approximations of the conditional mean parameter fields are then obtained from the cokriging estimator. Correlation between parameters and pressure – moisture content perturbations is seen to be strongly dependent on mean pressure or moisture content. High correlation between parameters and pressure data was obtained under saturated or near saturated flow conditions, providing accurate estimation of saturated hydraulic conductivity, while moisture content measurements provided accurate estimation of the pore size distribution parameter under unsaturated flow conditions.     

A geostatistically based inverse model for three-dimensional variably saturated flow

We present a geostatistically based inverse model for characterizing heterogeneity in parameters of unsaturated hydraulic conductivity for three-dimensional flow. Pressure and moisture content are related to perturbations in hydraulic parameters through cross-covariances, which are calculated to first-order. Sensitivities needed for covariance calculations are derived using the adjoint state sensitivity method. Approximations of the conditional mean parameter fields are then obtained from the cokriging estimator. Correlation between parameters and pressure – moisture content perturbations is seen to be strongly dependent on mean pressure or moisture content. High correlation between parameters and pressure data was obtained under saturated or near saturated flow conditions, providing accurate estimation of saturated hydraulic conductivity, while moisture content measurements provided accurate estimation of the pore size distribution parameter under unsaturated flow conditions.     

The effect of non divergence-free velocity fields on field scale ground water solute transport

Solute plume subjected to field scale hydraulic conductivity heterogeneity shows a large dispersion/macrodispersion, which is the manifestation of existing fields scale heterogeneity on the solute plume. On the other hand, due to the scarcity of hydraulic conductivity measurements at field scale, hydraulic conductivity heterogeneity can only be defined statistically, which makes the hydraulic conductivity a random variable/function. Random hydraulic conductivity as a parameter in flow equation makes the pore flow velocity also random and the ground water solute transport equation is a stochastic differential equation now. In this study, the ensemble average of stochastic ground water solute transport equation is taken by the cumulant expansion method in order to upscale the laboratory scale transport equation to field scale by assuming pore flow velocity is a non stationary, non divergence-free and unsteady random function of space and time. Besides the stochastic explanation of macrodispersion and the velocity correction term obtained by Kavvas and Karakas (J Hydrol 179:321–351, 1996) before a new velocity correction term, which is a function of mean pore flow velocity divergence, is obtained in this study due to strict second order cumulant expansion (without omitting any term after the expansion) performed. The significance of the new velocity correction term is investigated on a one dimensional transport problem driven by a density dependent flow field.     

Anisotropic dispersive Henry problem

The Henry problem has played a key role in our understanding of seawater intrusion into coastal aquifers and in benchmarking density dependent flow codes. This paper seeks to modify Henry’s problem to ensure sensitivity to density variations and vertical salinity profiles that resemble field observations. In the proposed problem, the “dispersive Henry problem”, mixing is represented by means of the traditional Scheidegger dispersion tensor (dispersivity times water flux). Anisotropy in the hydraulic conductivity is acknowledged and Henry’s seaside boundary condition of prescribed salt concentration is replaced by a flux dependent boundary condition, which represents more realistically salt transport across the seaside boundary. This problem turns out to be very sensitive to density variations and its solution gets closer to reality. However, an improvement in the traditional Henry problem (gain in sensitivity and realism) can be also achieved if the value of the Peclet number is significantly reduced.Although the dispersive problem lacks an analytical solution, it can shed light on flow in coastal aquifers. It provides significant information about the factors controlling seawater penetration, width of the mixing zone and influx of seawater. The width of the mixing zone depends basically on dispersion with longitudinal and transverse dispersion controlling different parts of the mixing zone but displaying similar overall effects. Toe penetration is mainly controlled by the horizontal permeability and by the geometric mean of the dispersivities. Finally, transverse dispersivity and the geometric mean of the hydraulic conductivity are the leading parameters controlling the amount of saltwater that enters the aquifer.     

Coordinate mapping of analytical contaminant transport solutions to non-uniform flow fields

Existing analytical solutions to 2D and 3D contaminant transport problems are limited by the mathematically convenient assumption of uniform flow. An approximate method is developed herein for coordinate mapping of 2D (vertically-averaged) transport solutions to non-uniform steady-state irrotational and divergence-free flow fields in single-layer aquifers. The method enables existing analytical transport solutions to be applied to aquifer systems with wells, non-uniform saturated thickness, surface water features, and (to a limited degree) heterogeneous hydraulic conductivity and recharge. This mass-conservative coordinate mapping approach is inexact in its approximation of the dispersion process but is still sufficiently accurate for many simple flow systems. The degree of model error is directly proportional to the variation of velocity magnitude within the domain. These mapped analytical solutions are compared to numerical simulation results and the coordinate mapping errors are investigated. The methods described herein may be used in the traditional capacity of analytical transport models, i.e., screening and preliminary site assessment, without sacrificing accuracy by assuming locally uniform flow conditions or applying an ad-hoc coordinate transformation. The solutions benefit from the traditional advantages of analytical methods, particularly the removal of artifacts due to spatial and temporal discretization: no time-stepping or numerical discretization is required.     

Steady state flow with hydraulic conductivity change around large diameter wells

Darcian flow law in aquifers assumes that the aquifer hydraulic conductivity is constant and the groundwater movement is due only to the piezometric level changes through hydraulic gradient. In practice, after the well development the aquifer just around the well has comparatively larger hydraulic conductivity and gradient. Patchy aquifer solutions in the literature consider sudden hydraulic conductivity changes with distance for the steady state flow. The change of transmissivity is demonstrated by the application of slope‐matching procedure to actual field data. It is the main purpose of this paper to derive simple analytical expressions for aquifer parameter evaluations with steadily decreasing hydraulic conductivity around the well. Spatial nonlinear hydraulic conductivity changes around a large‐diameter well within the depression cone of a confined aquifer are considered as exponentially decreasing functions of the radial distance. Copyright © 2010 John Wiley & Sons, Ltd.