共查询到20条相似文献，搜索用时 15 毫秒
1.
A mathematical model describing the hydraulic head distribution for a constant‐head test performed in a well situated at the centre of a patchy aquifer is presented. The analytical solution for the mathematical model is derived by the Laplace transforms and the Bromwich integral method. The solution for the hydraulic head has been shown to satisfy the governing equations, related boundary conditions, and continuity requirements for the hydraulic head and flow rate at the interface of the patch and outer regions. An efficient numerical approach is proposed to evaluate the solution, which has an integral covering an integration range from zero to infinity and an integrand consisting the product and square of the Bessel functions. This solution can be used to produce the curves of dimensionless hydraulic head against dimensionless time for investigating the effect of the contrast of formation properties on the dimensionless hydraulic head distribution. Define the ratio of outer‐region transmissivity to patch‐region transmissivity as α. The dimensionless hydraulic head for α=0.1 case is about 2.72 times to that for α=10 case at dimensionless large time (e.g. τ?106) when the dimensionless distance (ρ) equals 10. The results indicate that the hydraulic head distribution highly depends on the hydraulic properties of two‐zone formations. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
2.
工程建设中当距离抽水井r=rb处水位基本没有变化或不受抽水影响时,或当此处存在止水帷幕时,含水层系统视为侧向有限延伸,rb为有限半径。为此,构建更加符合工程实际的侧向有限延伸的典型弱透水层-承压水层系统中非完整井非稳定流计算模型,同时考虑井径和井储效应的影响,应用Laplace变换和分离变量法得到了水位降深在拉氏空间下的解析解,并应用拉氏数值逆变换Stehfest法得到真实空间下的水位降深。新建立的解析解可以进一步退化为诸多已有解,并进一步将其与已知解和有限元数值解进行对比,验证了所得解的正确性和可靠性。基于新建解重点分析了侧向边界和井的完整性对承压水层水位降深的影响。结果表明:含水层系统的侧向有限边界仅对抽水后期的水位降深影响明显,含水层系统侧向无限延伸情况下的水位降深要大于情形1(在r=rb处为定水头边界)且明显小于情形2(在r=rb处为不透水边界)下的水位降深,rb越小,两者之间的误差越大;抽水井的完整性对整个抽水期间不同情形下的水位降深均有明显的影响,承压含水层顶板处的水位降深随着抽水井滤管的长度和埋深的增加而减小。 相似文献
3.
In this study, non-Darcian flow to a larger-diameter partially penetrating well in a confined aquifer was investigated. The flow in the horizontal direction was assumed to be non-Darcian and described by the Izbash equation, and the flow in the vertical direction was assumed to be Darcian. A linearization procedure was used to approximate the nonlinear governing equation. The Laplace transform associated with the finite cosine Fourier transform was used to solve such non-Darcian flow model. Both the drawdowns inside the well and in the aquifer were analyzed under different conditions. The results indicated that the drawdowns inside the well were generally the same at early times under different conditions, and the features of the drawdowns inside the well at late times were similar to those of the drawdowns in the aquifer. The drawdown in the aquifer for the non-Darcian flow case was larger at early times and smaller at late times than their counterparts of Darcian flow case. The drawdowns for a partially penetrating well were the same as those of a fully penetrating well at early times, and were larger than those for a fully penetrating well at late times. A longer well screen resulted in a smaller drawdown in the aquifer at late times. A larger power index n in the Izbash equation resulted in a larger drawdown in the aquifer at early times and led to a smaller drawdown in the aquifer at late times. A larger well radius led to a smaller drawdown at early times, but it had little impact on the drawdown at late times. The wellbore storage effect disappears earlier when n is larger. 相似文献
4.
This study investigated non‐Darcian flow to a well in a leaky aquifer considering wellbore storage and a finite‐thickness skin. The non‐Darcian flow is described by the Izbash equation. We have used a linearization procedure associated with the Laplace transform to solve such a non‐Darcian flow model. Besides, the Stehfest method has been used to invert the Laplace domain solutions for the drawdowns. We further analyzed the drawdowns inside the well for different cases. The results indicated that a smaller BD results in a smaller drawdown at late times and the leakage has little effect on the drawdown inside the well at early times, where BD is a dimensionless parameter reflecting the leakage. We have also found that the flow for the negative skin case approaches the steady‐state earlier than that for the positive skin. In addition, the drawdown inside the well with a positive skin is larger than that without skin effect at late times, and a larger thickness of the skin results in a greater drawdown inside the well at late times for the positive skin case. A reverse result has been found for the negative skin case. Finally, we have developed a finite‐difference solution for such a non‐Darcian flow model and compared the numerical solution with the approximate analytical solution. It has been shown that the linearization procedure works very well for such a non‐Darcian flow model at late times, and it underestimates the drawdowns at early times. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
5.
Mazda Kompani-Zare Nozar Samani Siavash Behrooz-Koohenjani 《Hydrogeology Journal》2009,17(5):1175-1187
Large diameter fully cased wells that gain water from the bottom are often dug in sandy and collapsible aquifers. They have cylindrical vertical walls lined with brick or concrete. The well bottom is partially filled with aquifer material through which the flow is vertically upward. When the vertical hydraulic gradient reaches a critical value, quicksand occurs and the well structure can be destroyed. Another difficulty encountered is drawdown in the wellbore and the drying up of the well. To overcome these problems, the flow around and beneath these wells is numerically simulated. The simulation results are used to investigate the effect of well and aquifer parameters on quicksand and drawdown. For practical purposes, the dimensionless drawdown-time and dimensionless vertical gradient-time curves are developed. It was found that the ratio of filling material thickness to well radius affects the shape of these type curves. The type curves may be used to predict the time after pumping commences when quicksand occurs and the well dries up. They are also useful to design the safe pumping rate and duration as well as the optimum well radius. These are demonstrated by analyzing the pumping test data from a case study in the arid Chah Kutah region, southern Iran. 相似文献
6.
Non-Darcian flow to a partially penetrating well in a confined aquifer with a finite-thickness skin was investigated. The Izbash equation is used to describe the non-Darcian flow in the horizontal direction, and the vertical flow is described as Darcian. The solution for the newly developed non-Darcian flow model can be obtained by applying the linearization procedure in conjunction with the Laplace transform and the finite Fourier cosine transform. The flow model combines the effects of the non-Darcian flow, partial penetration of the well, and the finite thickness of the well skin. The results show that the depression cone spread is larger for the Darcian flow than for the non-Darcian flow. The drawdowns within the skin zone for a fully penetrating well are smaller than those for the partially penetrating well. The skin type and skin thickness have great impact on the drawdown in the skin zone, while they have little influence on drawdown in the formation zone. The sensitivity analysis indicates that the drawdown in the formation zone is sensitive to the power index (n), the length of well screen (w), the apparent radial hydraulic conductivity of the formation zone (K r2), and the specific storage of the formation zone (S s2) at early times, and it is very sensitive to the parameters n, w and K r2 at late times, especially to n, while it is not sensitive to the skin thickness (r s). 相似文献
7.
Non-Darcian flow to a well in a leaky aquifer was investigated using a finite difference method. Flow in the leaky aquifer is assumed to be non-Darcian and horizontal, while flow in the aquitard is assumed to be Darcian and vertical. The Forchheimer equation was employed to describe the non-Darcian flow in the aquifer. The finite difference solution was compared with the solution of Birpinar and Sen (2004). The latter overestimates the drawdown at early times and underestimates the drawdown at late times; also, the impact of β D on the drawdown depends on the value of B D, where β D is a dimensionless turbulent factor in the Forchheimer equation and B D is the dimensionless leakage parameter. The impact of leakage on drawdown is similar to that of Darcian flow. A sensitivity analysis indicated that the drawdown is very sensitive to the change in the dimensionless well radius r cD when B D is relatively large, while it is sensitive to the change in B D when B D is relatively small. The numerical solution has been applied to analyze the pumping test data in Chaj-Doab area of Pakistan. Birpinar ME, Sen Z (2004) Forchheimer groundwater flow law type curves for leaky aquifers. J Hydrol Eng 9(1):51??9 相似文献
8.
Analytical approach for time‐dependent groundwater inflow into shield tunnel face in confined aquifer 
下载免费PDF全文
下载免费PDF全文 Prediction of time‐dependent groundwater inflow into a shield tunnel is a significant task facing engineers. Published literature shows that there is no available method with which to predict time‐dependent groundwater inflow into a tunnel. This paper presents a prediction approach for time‐dependent groundwater inflow into a tunnel in both anisotropic and isotropic confined aquifers. The proposed solution can predict groundwater inrush from the tunnel cutting face. To obtain the time‐dependent groundwater flow quantity, the concept of a horizontal‐well pumping test based on the theory of a point source is adopted. Multiple factors, eg, drawdown, thickness of aquifer, conductivities, and specific storage, are taken into account. Both groundwater inflow to the cross section of a tunnel face in the y‐z plane and total tunnel inflow are obtained. Based on the proposed approach, the time‐dependent groundwater inflow to a tunnel can be classified as either a uniform or non‐uniform flow. The proposed approach is applied to analyse groundwater inflow of 2 field cases: (1) Metro line No. 7, Guangzhou City and (2) an underground tunnel in Huizhou, Guangdong Province. Results show that the proposed method can predict the measured values, and drawdown‐related curves are also derived. In addition, the calculated results also reveal that the effect of hydraulic conductivity kz on the total groundwater inflow differs from that of hydraulic conductivities kx and ky and the thickness of the aquifer. 相似文献
9.
抽水井附近由于流速过快往往发生非达西流,而远离抽水井随着流速下降又变为达西流.为了描述这些特征,建立了承压含水层中非完整井附近“非达西-达西”两区渗流模型,即距离抽水井较近的区域由于流速较快假设发生非达西渗流,并利用Izbash公式刻画,而距离抽水井较远由于流速较慢假设仍然满足达西定律,含水层中垂向流速较小也利用达西定律描述.通过线性化近似方法结合Laplace变换和有限Fourier余弦变换对模型进行了求解,分析探讨了该两区模型下水位降深曲线特征.结果表明:抽水初期,非达西渗流区域水位降深与全非达西渗流模型结果吻合,而抽水后期两区模型非达西渗流区域的水位降深与全达西模型水位降深基本一致,但大于全非达西渗流模型的水位降深;抽水初期,两区模型中达西渗流区域的水位降深比全达西渗流模型结果大,但比全非达西渗流模型结果小;对不同时间的水位降深随井距变化曲线分析发现非达西渗流区域水位降深随Izbash公式中的幂指数n增大而减小,而在达西渗流区域水位降深基本不受n值的影响.研究成果对非完整井抽水试验参数反演具有重要理论意义. 相似文献
10.
Solutions are presented for the problem of isothermal dessiccation shrinkage in a double‐layer porous partially saturated medium. The rheological model taken into account is linear poroelastic. Hence the analysis is mainly focused on hydromechanical coupling effects and contrasts of mechanical and hydraulic properties between two materials: a low thickness skin comprised between the outer boundary and the reference porous material. Three one‐dimensional ideal structures are taken into account: a wall of finite thickness (cartesian geometry), a thick cylinder and a thick sphere. The solution of the time‐dependent problem is arrived at by applying Laplace transforms to the field variables. Exact solutions are obtained in Laplace transform space using Mathematica© to solve the field equations whilst taking into account the continuity equations at the interface and the boundary conditions. The Talbot's modified algorithm has been performed to invert the Laplace transform solutions. A bibliographical and numerical study shows that this method is remarkably precise, stable and close to the analytical inversion. Results are presented using poroelastic data representative of a concrete material and involve a strong coupling effect between hydraulical and mechanical behaviours. A first approach elastic modelling of degradation process have been presented using a thin outer layer. Apart from emphasising the semi‐explicit solution utility due to accurate speed calculation, this paper deals with more complex problems than those which can be solved using purely analytical solutions. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
11.
An analytical solution is proposed for transient flow and deformation coupling of a fluid‐saturated poroelastic medium within a finite two‐dimensional (2‐D) rectangular domain. In this study, the porous medium is assumed to be isotropic, homogeneous, and compressible. In addition, the point sink can be located at an arbitrary position in the porous medium. The fluid–solid interaction in porous media is governed by the general Biot's consolidation theory. The method of integral transforms is applied in the analytical formulation of closed‐form solutions. The proposed analytical solution is then verified against both exact and numerical results. The analytical solution is first simplified and validated by comparison with an existing exact solution for the uncoupled problem. Then, a case study for pumping from a confined aquifer is performed. The consistency between the numerical solution and the analytical solution confirms the accuracy and reliability of the analytical solution presented in this paper. The proposed analytical solution can help us to obtain in‐depth insights into time‐dependent mechanical behavior due to fluid withdrawal within finite 2‐D porous media. Moreover, it can also be of great significance to calibrate numerical solutions in plane strain poroelasticity and to formulate relevant industry norms and standards. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
12.
Anisotropy and bedding effects on the hydro geological regime in a confined aquifer to design an appropriate dewatering system 总被引：1，自引：1，他引：0
A. Aryafar F. Doulati Ardejani Ph.D. 《International Journal of Environmental Science and Technology》2009,6(4):563-570
Prediction of groundwater inflow into mining excavations is very important in order to design an effective dewatering system to keep the mine workings dry and create prolonged cone of depression. The effects of anisotropy ratio and bedding on the hydraulic head and drawdown curves of a dewatering test carried out in a fully penetrating well in a confined aquifer have been investigated. An existing numerical finite element model has been used to perform the simulations. The results of the numerical model are compared to those from analytical Jacob and Lohman solution for estimating hydraulic heads and drawdown curves. It was found that the anisotropy ratio and bedding should not have a significant effect on drawdown and the quantity of inflow into a confined aquifer. It was further found that taking the simultaneous effects of anisotropy and bedding into account reduces the differences in the results of analytical and numerical methods. Comparison of the field data and model predictions showed that, the modelling results for a three layer anisotropic aquifer fit well to the field data than those results obtained for a single layer aquifer and the relative error decreased from 4.81 % to 2.98 %. 相似文献
13.
The flow rate to fully screened, partially penetrating wells in an unconfined aquifer is numerically simulated using MODFLOW 2000, taking into account the flow from the seepage face and decrease in saturated thickness of the aquifer towards the well. A simple three-step method is developed to find the top of the seepage face and hence the seepage-face length. The method is verified by comparing it with the results of previous predictive methods. The results show that the component of flow through the seepage face can supply a major portion of the total pumping rate. Variations in flow rate as a function of the penetration degree, elevation of the water level in the well and the distance to the far constant head boundary are investigated and expressed in terms of dimensionless curves and equations. These curves and equations can be used to design the degree of penetration for which the allowable steady pumping rate is attained for a given elevation of water level in the well. The designed degree of penetration or flow rate will assure the sustainability of the aquifer storage, and can be used as a management criterion for issuing drilling well permits by groundwater protection authorities. 相似文献
14.
Conventional graphical or computer methods for identifying aquifer parameters have their own inevitable limitations. This paper proposes a computer method based on a drawdown model and a heuristic approach of simulated annealing (SA) to determine the best‐fit aquifer parameters of the confined and unconfined aquifer systems. The drawdown model for the confined aquifer is the Theis solution and the unconfined aquifer is the Neuman solution. The estimated results of proposed method have better accuracy than those of the graphical methods and agree well with those of the computer methods based on the extended Kalman filter and Newton's method. Finally, the sensitivity analyses for the control parameters of SA indicate that the proposed method is very robust and stable in parameter identification procedures. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
15.
This study aims to model temperature distributions in an aquifer thermal extraction (ATE) system that contains a single extraction well in a thin confined aquifer. The aquifer is bounded by hot dry rocks with different thermomechanical properties and thicknesses. Based on the heat convection–conduction equation, a mathematical model is developed to describe the spatial and temporal temperature distributions of the ATE systems. The mechanisms of heat transfer in the model involve horizontal convection and thermal conduction in the aquifer, and vertical thermal conduction in both rocks. A semi‐analytical solution in dimensionless form is developed using the Laplace transform technique and its corresponding time‐domain result is computed by the modified Crump method. In addition, the steady‐state solution is obtained by applying the final value theorem. The simulation results from the semi‐analytical solution indicate that the aquifer temperature distributions are affected by aquifer thickness, the thermomechanical properties of the aquifer and rocks, geothermal gradient, outer boundary temperatures of the rocks, extraction rate, and operating time. The present solution can be used as a preliminary tool for assessing heat extraction efficiency in ATE systems. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
16.
A three‐dimensional finite‐element model of contaminant migration in fissured clays or contaminated sand which includes multiple sources of non‐equilibrium processes is proposed. The conceptual framework can accommodate a regular network of fissures in 1D, 2D or 3D and immobile solutions in the macro‐pores of aggregated topsoils, as well as non‐equilibrium sorption. A Galerkin weighted‐residual statement for the three‐dimensional form of the equations in the Laplace domain is formulated. Equations are discretized using linear and quadratic prism elements. The system of algebraic equations is solved in the Laplace domain and solution is inverted to the time domain numerically. The model is validated and its scope is illustrated through the analysis of three problems: a waste repository deeply buried in fissured clay, a storage tank leaking into sand and a sanitary landfill leaching into fissured clay over a sand aquifer. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
17.
A non-linear distribution of vertical displacement versus aquifer depth is calculated in the case of a partially penertrating well. For a fully penetrating well, however, a linear distribution is observed. The solution exhibits a vertically uniform horizontal displacement in the case of a fully penetrating well and, for a partially penetrating well, the maximum horizontal displacement occurs at the elevation of the well bottom. 相似文献
18.
Analyses for the response of a linear visco‐elastic system subjected to axi‐symmetric vertical circular loading are presented. Hankel transforms with respect to the radial spatial coordinate are used to reduce the three‐dimensional problem to that involving only a single spatial dimension, which is then discretized using the finite element method. Three techniques are employed to handle the time factor in the visco‐elastic material: (i) direct time integration; (ii) Fourier transforms; and (iii) Laplace transforms. These methods are compared and evaluated through their numerical results. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
19.
针对煤层顶板承压含水层涌水模式不清的问题,从煤层回采过程中顶板含水层涌水的时空变化特征入手,提出顶板含水层涌水量由静态储存量和动态补给量构成,认为静态储存量主要受来压步距、顶板垮落和导水裂隙(合称冒裂)影响区含水层厚度、含水层给水度控制,动态补给量主要受冒裂影响区外围含水层厚度、渗透性流场中水力梯度和过水断面面积控制;根据导水裂隙波及含水层情况,将顶板含水层涌水模式划分为井底进水的触及井涌水、井壁及井底进水的非完整井涌水和井壁进水的完整井涌水3种模式,并基于地下水渗流理论给出不同涌水模式下动态补给水量计算公式;针对以往疏放水钻孔数量多及疏放水量大的问题,以实现工作面顶板含水层静态储存量疏放后动态补给量可控为目的,提出冒裂区高度控制钻孔深度、单孔水位影响半径控制钻孔布置间距、钻孔疏放水量稳定时间控制超前疏放时间的疏放水钻孔优化设计理念,对疏放水及疏放钻孔布置进行优化,形成系统的顶板含水层水疏放体系。研究结果丰富了煤层顶板含水层涌水量计算和控制方法,对顶板水害防控具有实际的指导意义。 相似文献
20.
The fundamental solutions were obtained for step‐like point forces acting in three orthogonal directions and an instantaneous fluid point source in a fluid‐saturated, porous, infinite solid of transversely isotropic elasticity and permeability. After expressing the governing equations in the form of matrix in the Laplace space, we employed Kupradze's method together with the triple Fourier transforms. This method reduces the simultaneous partial differential equations with respect to three displacement components and a pore fluid pressure to a differential equation in terms of only one potential scalar function, which can be operationally solved in the transformed space. After the Laplace inversion of the potential, the residue theorem was applied to its Fourier inverse transform with respect to one of the transformation variables. The Fourier transforms with respect to two other variables were rewritten into the Hankel transforms. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献