Obtaining the steady‐state drawdown solutions of constant‐head and constant‐flux tests

The solutions of constant‐head and constant‐flux tests are commonly used to predict the temporal or spatial drawdown distribution or to determine aquifer parameters. Theis and Thiem equations, for instance, are well‐known transient and steady‐state drawdown solutions, respectively, of the constant‐flux test. It is known that the Theis equation is not applicable to the case where the aquifer has a finite boundary or the pumping time tends to infinity. On the other hand, the Thiem equation does not apply to the case where the aquifer boundary is infinite. However, the issue of obtaining the Thiem equation from the transient drawdown solution has not previously been addressed. In this paper, the drawdown solutions for constant‐head and constant‐flux tests conducted in finite or infinite confined aquifers with or without consideration of the effect of the well radius are examined comprehensively. Mathematical verification and physical interpretation of the solutions to these two tests converging or not converging to the Thiem equation are presented. The result shows that there are some finite‐domain solutions for these two tests that can converge to the Thiem equation when the time becomes infinitely large. In addition, the time criteria to give a good approximation to the finite‐domain solution by the infinite‐domain solution and the Thiem equation are investigated and presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Surface motion of a cylindrical hill of circular-arc cross-section for incident plane SH waves

A closed-form solution of two-dimensional scattering of plane SH waves by a cylindrical hill of circular-arc cross-section in a half-space is presented using the wave functions expansion method. The solution is reduced to solving a set of infinite linear algebraic equations, using the auxiliary functions and the exterior region form of the Graf's addition theorem. Numerical solutions are obtained by truncation of the infinite equations and their accuracies are demonstrated by convergence of the numerical results and by the extent to which the numerical results fit the exact boundary conditions with increasing the truncation order. The numerical results for some typical cases are then presented for checking accuracies of various numerical methods. The effects of the height-to-width ratio of the hill on surface ground motion are finally illustrated.     

Vibration of vertical rectangular plate in contact with water on one side

In this paper, the vibratory characteristics of a rectangular plate in contact with water on one side are studied. The elastic plate is considered to be a part of a vertical rectangular rigid wall in contact with water, the edges of which are elastically restrained and parallel to those of the rigid wall. The location and size of the plate on the rigid wall may vary arbitrarily. The water with a free surface is in a rectangular domain infinite in the length direction. The effects of free surface waves, compressibility of the water and the hydrostatic water pressure are neglected in the analysis. An analytical‐Ritz method is developed to analyse the interaction of the plate–water system. First of all, by using the method of separation of variables and the method of Fourier series expansion, the exact expression of the motion of water is derived in the form of integral equations including the dynamic deformation of the plate. Then the Rayleigh–Ritz approach is used to derive the eigenfrequency equation of the system via the variational principle of energy. By selecting beam vibrating functions as the admissible functions of the plate, the added virtual mass incremental (AVMI) matrices for plate vibration are obtained. The convergency studies are carried out. The effects of some parameters such as the depth and width of water, the support stiffnesses, location and aspect ratio of the plate and the plate–water size and density ratios on the eigenfrequencies of the plate–water system are investigated. Several numerical examples are given. The validity of AVMI factor approach is also confirmed by comparing the AVMI factor solutions with the analytical‐Ritz solutions. The results show that the approach presented here can also be used as excellent approximate solutions for rectangular plates in contact with water of infinite width and/or infinite depth. Copyright © 2000 John Wiley & Sons, Ltd.     

任意圆弧形凸起地形对平面SH波的散射

本文采用波函数展开方法提出弹性半空间表面一任意圆弧形凸起边界对平面ＳＨ波二维散射的封闭级数解答，利用引入的辅助函数和推广的外域型Ｇｒａｆ加法公式将解答归结为一代无穷代数方程组的求解，解答的数值结果可由无穷方程的截断计算得出，文中从级数项数增加时计算结果的收敛以及边界条件的收敛满足两方面检验了截断计算的精度，并指出了位移解答的收敛速度与主要参数之间的关系，同时，对凸起角点的应力奇异性及其对解答精度的     

地表下圆形夹塞区出平面散射对地面运动的影响

采用波函数展开方法，提出弹性半空间埋藏的圆形弹性夹塞区对平面ＳＨ波二维散射的封闭级数解答．利用内域型Ｇｒａｆ加法公式，将解答归结为一组无穷代数方程组的求解．通过对无穷方程组的截断计算，得到了解答的数值结果；并通过边值条件的收敛，检验了截断计算的精度．最后给出若干典型算例说明夹塞区产生的散射场对地表运动的复杂影响．     

Use of the Fourier-Laplace transform and of diagrammatical methods to interpret pumping tests in heterogeneous reservoirs

Advances in computer power and in reservoir characterization allow simulation of pressure transients in complex reservoirs generated stochastically. Generally, interpretation of these transient gives useful information about the reservoir hydraulic properties: a major goal is to interpret these transients in the stochastic context. First we ensemble average the pressure over all the random permeability field realizations to derive an equation which drives the ensemble averaged pressure. We use the Fourier transform in space and the Laplace transform in time, in conjuction with a perturbation series expansion in successive powers of the permeability fluctuations to obtain an explicit solution. The Nth order term of this series involves the hydrodynamic interaction between N permeability heterogeneities and after averaging we obtain an expansion containing correlation functions of permeability fluctuations of increasing order.Next, Feynman graphs are introduced allowing a more attractive graphical interpretation of the perturbation series. Then series summation techniques are employed to reduce the graph number to be summed at each order of the fluctuation expansion. This in turn gives useful physical insights on the homogenization processes involved. In particular, it is shown that the sum of the so-called ‘one-particle irreducible graphs’ gives the kernel of a linear integro-differential equation obeyed by the ensemble average pressure. All the information about the heterogeneity structure is contained in this renormalized kernel, which is a limited range function.This equation on its own is the starting point of useful asymptotic results and approximations. In particular it is shown that interpretation of pumping tests yields the steady-state equivalent permeability after a sufficiently long time for an infinite reservoir, as expected.     

Analytical solutions for transport of decaying solutes in rivers with transient storage

Analytical solutions are presented for solute transport in rivers including the effects of transient storage and first order decay. The solute transport model considers an advection–dispersion equation for transport in the main channel linked to a first order mass exchange between the main channel and the transient storage zones. In case of a conservative tracer, it is shown that different analytical solutions presented in the literature are mathematically identical. For non-conservative solutes, first order decay reactions are considered with different reaction rate coefficients in the main river channel and in the dead zones. New analytical solutions are presented for different boundary conditions, i.e. instantaneous injection in an infinite river reach, and variable concentration time series input in a semi-infinite river reach. The correctness and accuracy of the analytical solutions is verified by comparison with the OTIS numerical model. The results of analytical and numerical approaches compare favourably and small differences can be attributed to the influence of boundary conditions. It is concluded that the presented analytical solutions for solute transport in rivers with transient storage and solute decay are accurate and correct, and can be usefully applied for analyses of tracer experiments and transport characteristics in rivers with mass exchange in dead zones.     

饱和土中孔列对平面快纵波的隔离效应

在Biot饱和土波动理论的基础上,用单排孔列作为屏障,对全空间饱和土中人射平面快纵波进行了隔离分析.采用波函数展开法,将人射波和散射波的势函数展开成Fourier-Bessel函数的无穷级数的形式,利用一组可以求解多散射问题的圆柱坐标系统和界面处位移连续的条件,得到问题的理论解.对影响隔离效果的孔间距和饱和土的渗透性等...     

Spherical harmonic analysis of the Navier-Stokes equation in magnetofluid dynamics

The equation of motion (Navier-Stokes equation) for a uniformly rotating, compressible, magnetic, viscous fluid is analyzed in terms of infinite series of spherical surface harmonics. Differential equations are obtained for the radial functions of the poloidal and toroidal harmonics of the velocity, corresponding to those obtained by Bullard and Gellman for the magnetic field from the electromagnetic induction equation. This new analysis opens the way for the dynamical problem of electromagnetic induction in the earth's core to be considered by the spherical harmonic method.     

Stochastic perturbation analysis of groundwater flow. Spatially variable soils,semi-infinite domains and large fluctuations

As is well known, a complete stochastic solution of the stochastic differential equation governing saturated groundwater flow leads to an infinite hierarchy of equations in terms of higher-order moments. Perturbation techniques are commonly used to close this hierarchy, using power-series expansions. These methods are applied by truncating the series after a finite number of terms, and products of random gradients of conductivity and head potential are neglected. Uncertainty regarding the number or terms required to yield a sufficiently accurate result is a significant drawback with the application of power series-based perturbation methods for such problems. Low-order series truncation may be incapable of representing fundamental characteristics of flow and can lead to physically unreasonable and inaccurate solutions of the stochastic flow equation. To support this argument, one-dimensional, steady-state, saturated groundwater flow is examined, for the case of a spatially distributed hydraulic conductivity field. An ordinary power-series perturbation method is used to approximate the mean head, using second-order statistics to characterize the conductivity field. Then an interactive perturbation approach is introduced, which yields improved results compared to low-order, power-series perturbation methods for situations where strong interactions exist between terms in such approximations. The interactive perturbation concept is further developed using Feynman-type diagrams and graph theory, which reduce the original stochastic flow problem to a closed set of equations for the mean and the covariance functions. Both theoretical and practical advantages of diagrammatic solutions are discussed; these include the study of bounded domains and large fluctuations.     

Scattering of plane SH waves by a cylindrical alluvial valley of circular-arc cross-section

A closed-form solution of two-dimensional scattering of plane SH waves by a cylindrical alluvial valley of circular-arc cross-section in a half-space is presented using the wave functions expansion. The solution is reduced to solving a set of infinite linear algebraic equations using the exterior region form of Graf's addition theorem. Numerical solutions are obtained by truncation of the infinite equations and their accuracies are demonstrated by convergence of the numerical results to the exact boundary condition with the increasing of the truncation order. The present solution is compared with the existing one presented by Todorovska and Lee for the same problem and their differences are analysed. Complicated effects of the depth-to-width ratio of the alluvial valley on surface ground motion are finally illustrated.     

Exact solutions ofSH wave equation in transversely isotropic inhomogeneous elastic media

Summary TheSH wave equation in a transversely isotropic inhomogeneous elastic medium, where the elastic parameters and density are functions of vertical coordinate, is considered. A general procedure is given for finding the inhomogeneities for which the equation can be solved in terms of hypergeometric, Whittaker, Bessel and exponential functions. A few simple inhomogeneities and the corresponding solutions in terms of these transcendental functions are presented.     

Scattering of Plane P Waves by Circular-arc Canyon Topography:High-frequency Solution

Through the Fourier-Bessel series expansion of wave functions,the analytical solution to the two-dimensional scattering problem of incidental plane P waves by circular-arc canyon topography with different depth-to-width ratio is deduced.Unlike other existing analytical solutions,in order to ensure that the analytical solution is valid for higher frequency incident waves,the asymptotic properties of cylindrical functions are in this paper introduced to directly determine the unknown coefficients of scattering waves,avoiding the solution of linear equation systems and corresponding numerical issues,which in turn expand the frequency band in which the analytical solution is valid.Comparison with other existing analytical solutions demonstrates that the proposed analytical solution is correct.Furthermore,the scattering effects of a circular-arc canyon on the incident plane P wave are analyzed in a comparatively broad frequency band.     

Explicit expressions and numerical calculations for the Fréchet and second derivatives in 2.5D Helmholtz equation inversion

In order to perform resistivity imaging, seismic waveform tomography or sensitivity analysis of geophysical data, the Fréchet derivatives, and even the second derivatives of the data with respect to the model parameters, may be required. We develop a practical method to compute the relevant derivatives for 2.5D resistivity and 2.5D frequency-domain acoustic velocity inversion. Both geophysical inversions entail the solution of a 2.5D Helmholtz equation. First, using differential calculus and the Green's functions of the 2.5D Helmholtz equation, we strictly formulate the explicit expressions for the Fréchet and second derivatives, then apply the finite-element method to approximate the Green's functions of an arbitrary medium. Finally, we calculate the derivatives using the expressions and the numerical solutions of the Green's functions. Two model parametrization approaches, constant-point and constant-block, are suggested and the computational efficiencies are compared. Numerical examples of the derivatives for various electrode arrays in cross-hole resistivity imaging and for cross-hole seismic surveying are demonstrated. Two synthetic experiments of resistivity and acoustic velocity imaging are used to illustrate the method.     

Thermo-elastic interactions in an infinite elastic solid due to a concentrated transient heat source

Summary The present paper is concerned with the determination of thermo-elastic stress and temperature distribution in an infinite elastic solid when it is subjected to a concentrated transient heat source taking into account the effect of coupling. The fundamental partial differential equation for the thermo-elastic potential is solved by means of the operational method. The solutions are obtained for small values of time in terms of known functions such as complementary error functions and associated complementary error functions and have the form of power series with respect to the parameter of coupling . Since this parameter is generally very small for most of the metals, the expressions for stress and temperature distribution have been obtained by considering terms upto linear .     

Nonlinear aspects of phreatic water filtration

Horizontal gravity filtration of groundwater in soil is considered. Under Boussinesq approximation, the problem is reduced to a one-dimensional nonlinear parabolic equation in phreatic water level. The problem of linearizing the original equation is discussed. The comparison of gravity-filtration problem solutions in the nonlinear and linearized formulations shows considerable discrepancies to exist between the solutions, especially, for boundary problems with mixed boundary conditions, when the value of the function is not fixed on the right boundary. An analytical solution is obtained for steady-state flow from a water body into the soil with subsequent leakage into underlying beds. Two regimes are shown to exist: one with an infinite exponential tail, and another in the form of a finite groundwater mound. A new approach is proposed to the linearization problem—quasilinearization with the use of the Burgers equation.     

GENERALIZED KUNETZ-TYPE EQUATIONS*

A Kunetz equation is often used as the starting point in the development of solutions for the inversion of one-dimensional, noise free, normal incident seismograms, for which |r_o|= 1. In this paper we demonstrate a need for a Kunetz-type equation in which filtered signals can be used, so that noise effects can be reduced. We then show that an infinite number of Kunetz-type equations exist for the lossless wave equation in layered media. Finally, we show that it is indeed valid to formulate and solve the inverse problem using filtered signals.     

Stochastic space transforms in subsurface hydrology — Part 2: Generalized spectral decompositions and plancherel representations

In earlier publications, certain applications of space transformation operators in subsurface hydrology were considered. These operators reduce the original multi-dimensional problem to the one-dimensional space, and can be used to study stochastic partial differential equations governing groundwater flow and solute transport processes. In the present work we discuss developments in the theoretical formulation of flow models with space-dependent coefficients in terms of space transformations. The formulation is based on stochastic Radon operator representations of generalized functions. A generalized spectral decomposition of the flow parameters is introduced, which leads to analytically tractable expressions of the space transformed flow equation. A Plancherel representation of the space transformation product of the head potential and the log-conductivity is also obtained. A test problem is first considered in detail and the solutions obtained by means of the proposed approach are compared with the exact solutions obtained by standard partial differential equation methods. Then, solutions of three-dimensional groundwater flow are derived starting from solutions of a one-dimensional model along various directions in space. A step-by-step numerical formulation of the approach to the flow problem is also discussed, which is useful for practical applications. Finally, the space transformation solutions are compared with local solutions obtained by means of series expansions of the log-conductivity gradient.     

Diffraction of anti-plane SH waves by a semi-circular cylindrical hill with an inside concentric semi-circular tunnel

A closed-form analytic solution of two-dimensional scattering and diffraction of plane SH waves by a semicylindrical hill with a semi-cylindrical concentric tunnel inside an elastic half-space is presented using the cylindrical wave functions expansion method. The solution is reduced to solving a set of infinite linear algebraic equations. Fourier expansion theorem with the form of complex exponential function and cosine function is used. Numerical solutions are obtained by truncation of the infinite equations. The accuracy of the presented numerical results is carefully verified.