共查询到20条相似文献，搜索用时 88 毫秒
1.
地下水的补排主要包括垂向的地面入渗补给、蒸发排泄（蒸发可视为入渗的负值）及侧向的地表水补给、排泄。水文地质学最基本的问题之一——地下水可持续开釆量的评价准则,涉及补给的增量与排泄的减量,因此地下水开采的预测模型必须包含上述两类的补给、排泄因素,否则不能满足要求。然而,经典的Theis不稳定井流模型,即使在傍河抽水,也只有侧边界的补给、排泄作用,而不涉及上边界的地面入渗补给。这样一来,这个解析模型基本上不能够用于预测,而只能在旱季用于井流试验求取含水系统的参数。为此,文章的目标是发展具地面入渗补给的Theis不稳定潜水井流模型。对于潜水流问题,不能再用承压水流的以水头为应变量的方程来建立,应采用第二类线性化方法的势函数来建立潜水流问题。对于既有降雨入渗补给,又有抽水井作用的复杂的水文地质问题所概化数学模型的求解,采取的方法是把它分解成若干个简单的子模型问题求解,然后将其合成为原来复杂数学模型的解。基于质量守恒原理,假定渗流服从Darcy定律并满足Dupuit徦定建立了水流基本微分方程。然后对于两平行河流及一河流平行一隔水边界形成的两类条形区域,具地面均匀稳定入渗补给的井流问题,获得通用水位方程和几类常见的特定条件水位方程及其流量方程。此外,提出并采用“边界对边界的反映法”用以求解一河流平行一隔水边界条形区域的同一问题,减少了许多推导过程。最后,作为上述理论成果的初步应用,也是一个重要的应用,即在河水水质不能满足要求的河流附近,设有一口抽水井,计算该抽水井在不汲取河水的前提下的临界流量方程,获得具重要意义的结构简洁的关系式。该方程也可以用于滨海区的抽水井,在不发生海水入侵前提下的临界抽水流量计算。给出了上述条件不稳定井流过程某时刻的地下水流网图,其流网与文献中常见的傍河井流的流网相比,具显明的特征。 相似文献
2.
We present a second-order analytic solution [in terms of a heterogeneous log-transmissivity Y(r) = ln T(r)] for the hydraulic head field in a finite 2D confined heterogeneous aquifer under steady radial flow conditions assuming fixed head boundary conditions at the well and at a circular exterior boundary. The solution may be used to obtain the gradient used in calculation of solute transport to a well in a heterogeneous transmissivity field. The solution, obtained using perturbation methods coupled with Green's function techniques, leads us to postulate a more general form of the head for arbitrarily large-variance fields and may be used to obtain moment relations between the log-transmissivity and head under convergent flow conditions when Y(r) is expressed as a random space function. We present expressions for the mean head field when the log-transmissivity is Gaussian and conditioned on the transmissivity value at the well for an arbitrary ln T covariance. Finally, we look at the effect of parameter variations on the mean head behavior and present numerical simulations verifying the second-order mean head expressions. 相似文献
3.
4.
潜水水流的动态随机模拟是一个复杂而难解决的问题.通过建立二维潜水非稳定流模拟的摄动随机有限元模型, 把控制方程的主要参数渗透系数和给水度随机变量、及源汇项和边界条件看作随机变量.在充分考虑4种随机因素的条件下, 推导出求解潜水二维非稳定流均值和方差的9个方程; 重点介绍了不同方程数值离散的特殊处理方法.通过设定理想例子对模拟结果进行了分析, 表明随机变量中边界条件值方差、渗透系数方差变化对水头方差变化的影响很小, 给水度方差的变化对水头方差的变化影响很大.本模型考虑因素全面, 对一般的潜水非稳定流随机模拟都可应用.本研究给出了边界、渗透系数、给水度的随机因素对潜水动态模拟的影响, 丰富和补充了地下水运动的随机理论. 相似文献
5.
Darcian flow under/through a leaky cutoff wall: Terzaghi–Anderson's seepage problem revisited 
下载免费PDF全文
下载免费PDF全文 An analytical solution is obtained for 2‐D steady Darcian flow under and through a cutoff wall partially obstructing a homogeneous isotropic foundation of a dam. The wall is leaky; that is, flow across it depends on the ratio of hydraulic conductivity of the wall and the wall thickness that results in the third‐type (Robin) boundary condition along the wall, as compared with the Terzaghi problem for an impermeable wall. The Laplace equation for the hydraulic head is meshlessly solved in a non‐standard flow tube. A Fredholm equation of the second kind is obtained for the intensity of leakage across the wall. The equation is tackled numerically, by adjusted successive iterations. Flow characteristics (total Darcian discharge and its components through the wall and the window between the wall top and horizontal bedrock, stream function, head distribution, and Darcian velocity along the wall and tailwater bed) are obtained for various conductivity ratios, head drops across the structure, thicknesses of the foundation, and the degree of its blockage by the wall. Comparisons with the Terzaghi limit of an impermeable wall show that for common wall materials and thicknesses, the leakage may constitute tens of percent of the discharge under the dam. The through‐flow hydraulic gradients on a vertical wall face (Robin's boundary condition) as well as the exit gradients along a horizontal tailwater boundary (Dirichlet's boundary condition) acting for decades have deleterious impacts on dam stability because of potential heaving, piping, and mechanical–chemical suffusion. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
6.
Peter Indelman 《Computational Geosciences》2000,4(4):351-381
Average nonuniform flows in heterogeneous formations are modeled with the aid of the nonlocal effective Darcy's law. The mean head for flow toward source of instantaneous discharge in a heterogeneous medium of given statistics represents the fundamental solution of the average flow equation and is called the Mean Green Function (MGF). The general representation of the MGF is obtained for weakly heterogeneous formations as a functional of the logconductivity correlation function. For Gaussian logconductivity correlation, the MGF is derived in terms of one quadrature in time t and it is analyzed for isotropic media of any dimensionality d and for 3D axisymmetric formations. The MGF is further applied to determining the mean head distribution for flow driven by a continuous source of constant discharge. The large time asymptotic of the mean head is analyzed in details. 相似文献
7.
L. Ya. Aranovich I. V. Zakirov N. G. Sretenskaya T. V. Gerya 《Geochemistry International》2010,48(5):446-455
New experimental data on the solubility of NaCl in gaseous CO2 were obtained at pressures (P) of 30–70 MPa and temperatures of 623 and 673 K on experimental equipment making possible to sample a portion of the gas
in the course of the experiment. The new measures have demonstrated that the NaCl solubility increases with increasing temperature
(T) and pressure and is approximately four to five orders of magnitude higher than the saturated vapor pressure of NaCl at the
corresponding temperature. The paper also reports newly obtained experimental data on the equilibrium conditions of the reaction
of talc decomposition into enstatite and quartz at a variable H2O/NaCl ratio in the fluid. The results of the experiments validate the empirical equations previously suggested for H2O and NaCl activities in concentrated aqueous salt solutions that can be used in describing silica-saturated fluids at high
T-P parameters. A new empirical equation is suggested for the Gibbs free mixing energy in the H2O-CO2-NaCl ternary system, with the parameters of the equation calibrated against experimental data on phase equilibria in marginal
binary systems and on the location of the boundary of the region of homogeneous three-component fluid according to data on
synthetic fluid inclusions in quartz. 相似文献
8.
A new methodology is developed in assessing environmental impacts of desalination plants discharging brine into the ground.
The main environmental problem of the desalination of seawater is the brine disposal. The brine is commonly discharged into
the sea or injected into a saline aquifer. In the case of injection into the ground, it is necessary to design a disposal
system in a way that respects the environment and is sustainable. Laboratory and computational methods have been utilized
to simulate the unsteady three-dimensional (3D) phenomena of subsurface brine disposal. The computational software used is
SEAWAT, which is a 3D unsteady variable-density flow simulation model. The model is first used to simulate the laboratory
results, and good agreement is achieved. Then, hypothetical problems are designed and simulated of groundwater extraction
and brine disposal by desalination stations. The major purpose of these hypothetical problems is to delineate a methodology
and to create design charts for design and management of production and injection well fields for coastal desalination plants.
Several design charts have been developed with 36 scenarios for two well configurations created by four design parameters:
relative salt concentration (RSC), production and injection rates (Q
d
, Q
r
), well spacing (S), and simulation period (T). 相似文献
9.
By extending Darcy's law to the dry domain above the free surface and specifying the boundary condition on the potential seepage surfaces as Signorini's type, a partial differential equation (PDE) defined in the entire domain of interest is formulated for non‐steady seepage flow problems with free surfaces. A new parabolic variational inequality (PVI) formulation equivalent to the PDE formulation is then proposed, in which the flux part of the complementary condition of Signorini's type in the PDE formulation is transformed into natural boundary condition. Consequently, the singularity at the seepage points is eliminated and the difficulty in selecting the trial functions is significantly reduced. By introducing an adaptive penalized Heaviside function in the finite element analysis, the numerical stability of the discrete PVI formulation is well guaranteed. The proposed approach is validated by the existing laboratory tests with sudden rise and dropdown of water heads, and then applied to capture the non‐steady seepage flow behaviors in a homogeneous rectangular dam with five drainage tunnels during a linear dropdown of upstream water head. The non‐steady seepage flow in the surrounding rocks of the underground powerhouse in the Shuibuya Hydropower Project is further modeled, in which a complex seepage control system is involved. Comparisons with the in situ monitoring data show that the calculation results well illustrate the non‐steady seepage flow process during impounding and the operation of the reservoir as well as the seepage control effects of the drainage hole arrays and drainage tunnels. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
10.
Guang Y. Sheu 《国际地质力学数值与分析法杂志》2004,28(5):437-463
Earlier solutions (Bull. Seismol. Soc. Amer. 1985; 75 :1135–1154; Bull. Seismol. Soc. Amer. 1992; 82 :1018–1040) of deformations caused by the movements of a shear or tensile fault in an isotropic half-space for finite rectangular sources of strain nucleus have been extended for a transversely isotropic half-space. Results of integrating previous solutions (Int. J. Numer. Anal. Meth. Geomech. 2001; 25 (10): 1175–1193) of deformations due to a shear or tensile fault in a transversely isotropic half-space for point sources of strain nucleus over the fault plane are presented. In addition, a boundary element (BEM) model (POLY3D:A three-dimensional, polygonal element, displacement discontinuity boundary element computer program with applications to fractures, faults, and cavities in the Earth's crust. M.S. Thesis, Stanford University, Department of Geology, 1993; 62) is given. Different from similar researches (e.g. Thomas), the Akaike's view on Bayesian statistics (Akaike Information Criterion Statistics. D. Reidel Publication: Dordrecht, 1986) is applied for inverting deformations due to a fault to obtain displacement discontinuities on the fault plane. An example is given for checking displacements predicted by proposed analytical expressions. Another example is generated for the use of proposed BEM model. It demonstrates the effectiveness of this model in exploring displacement behaviours of a fault. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
11.
Kusuma G. Rao V. N. Lykossov A. Prabhu S. Sridhar E. Tonkacheyev 《Journal of Earth System Science》1996,105(3):227-260
An attempt has been made here to study the sensitivity of the mean and the turbulence structure of the monsoon trough boundary
layer to the choice of the constants in the dissipation equation for two stations Delhi and Calcutta, using one-dimensional
atmospheric boundary layer model withe-ε turbulence closure. An analytical discussion of the problems associated with the constants of the dissipation equation is
presented. It is shown here that the choice of the constants in the dissipation equation is quite crucial and the turbulence
structure is very sensitive to these constants. The modification of the dissipation equation adopted by earlier studies, that
is, approximating the Tke generation (due to shear and buoyancy production) in theε-equation by max (shear production, shear + buoyancy production), can be avoided by a suitable choice of the constants suggested
here. The observed turbulence structure is better simulated with these constants. The turbulence structure simulation with
the constants recommended by Aupoixet al (1989) (which are interactive in time) for the monsoon region is shown to be qualitatively similar to the simulation obtained
with the constants suggested here, thus implying that no universal constants exist to regulate dissipation rate.
Simulations of the mean structure show little sensitivity to the type of the closure parameterization betweene-l ande-ε closures. However the turbulence structure simulation withe-ε. closure is far better compared to thee-l model simulations. The model simulations of temperature profiles compare quite well with the observations whenever the boundary
layer is well mixed (neutral) or unstable. However the models are not able to simulate the nocturnal boundary layer (stable)
temperature profiles. Moisture profiles are simulated reasonably better. With one-dimensional models, capturing observed wind
variations is not up to the mark. 相似文献
12.
泥石流作为非牛顿体,屈服应力大,运动过程通常不稳定。前人建立了许多模型来研究沟床揭底和堰塞体溃决对泥石流不稳定动力过程的影响,沟岸侧蚀对泥石流不稳定动力过程的影响研究较少。通过侧蚀为主的模型和完全底蚀的模型两种水槽实验的对比,针对泥石流的动力过程展开研究。实验发现两种工况条件下泥石流正应力和孔隙水压力随着龙头高度沿程波动性的增长而相应地波动性增大,但侧蚀作用使得这种波动特征更加明显。通过力学分析,证明侧蚀作用导致泥石流龙头的阻力更大,但是龙身颗粒和龙头颗粒的速度差更大,使得龙头附加坡降更大,因此,侧蚀作用使得泥石流龙头的平均速度更快。泥石流龙头浓度和容重的不断增大,使得阻力不断增大,阻力和动力的动态平衡关系是泥石流不稳定运动的原因之一。 相似文献
13.
提出了一种改进的时间差分有限分析算法,它能够避免在求解水头剧烈变化时用一般有限分析法计算精度较差的缺陷.算法的稳定性好. 相似文献
14.
Hugo A. Loáiciga 《Mathematical Geosciences》2008,40(2):199-211
The equation of groundwater flow in marine island aquifers in which there is time-independent, spatially-variable recharge
and pumping is solved in closed form for rectangular, circular, and elliptical island geometries. The solution of the groundwater
flow equation is expressed in terms of the elevation of the phreatic surface within the flow domain. The depth of the seawater-freshwater
interface below mean sea level follows from the Dupuit–Ghyben–Herzberg relation. The method of solution presented in this
work relies on expanding the hydraulic head and forcing function (recharge and groundwater extraction) as Fourier series that
transforms the two-dimensional Poisson-type flow equations into second-order ordinary differential equations solvable using
classical theory. The important case of constant recharge (without groundwater extraction) leads to solutions in which the
hydraulic head is expressible as the product of a flow factor equal to the squared root of the ratio of recharge over hydraulic
conductivity times a geometric factor involving island shape parameters and flow boundary conditions. Estimability conditions
for the hydraulic conductivity are derived for the cases of constant recharge and spatially variable recharge with pumping. 相似文献
15.
边界元区域分解算法在地下水中的应用 总被引：3，自引：0，他引：3
将边界元方法与区域分解算法进行耦合，结合两种方法的优点，并将此方法应用到地下水中，并计算了承压含水层稳定流问题。对非稳定问题则需通过变量替换和Laplace变换，将含时间变量的问题化成只含空间变量的问题，再针对Laplace变换解求其逆变换即可。 相似文献
16.
17.
18.
This paper presents a boundary element method (BEM) procedure for a linear elastic fracture mechanics analysis in two‐dimensional anisotropic bimaterials. In this formulation, a displacement integral equation is only collocated on the uncracked boundary, and a traction integral equation is only collocated on one side of the crack surface. A fundamental solution (Green's function) for anisotropic bimaterials is also derived and implemented into the boundary integral formulation so that except for the interfacial crack part, the discretization along the interface can be avoided. A special crack‐tip element is introduced to capture the exact crack‐tip behavior. A computer program using FORTRAN has been developed to effectively calculate the stress intensity factors of an anisotropic bimaterial. This BEM program has been verified to have a good accuracy with previous studies. In addition, a central cracked bimaterial Brazilian specimen constituting cement and gypsum is prepared to conduct the Brazilian test under diametral loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
19.
The properties of the solution of the kinematic equation (Kepler’s equation) presented as a series in powers of a function
of time are considered for the case of unperturbed, rectilinear elliptical and rectilinear hyperbolic motion. Kepler’s equation
for unperturbed, rectilinear elliptical motion has the form E − sinE = z
3/6, where E is the eccentric anomaly determining the position in the orbit, z
3/6 is the mean anomaly, which is proportional to the time measured from an encounter, and its solution can be represented
as a series in powers of z. It is established that the coefficients of the series are positive. The asymptotic for the coefficients in the region of
convergence of the series is found, which covers the entire orbit. The series continues to converge over the entire boundary
of the circle of convergence. The kinematic equation for unperturbed, rectilinear hyperbolic motion has the form sinhH − H = ζ
3/6, where, as before, ζ
3 is proportional to time. The substitution E = iH, z = iζ reduces one equation to the other. The series for the solution in the hyperbolic case differs from the series for the elliptical
solution only in its alternating-sign coefficients. However, the region of convergence covers only part of the orbit in the
hyperbolic case. 相似文献
20.