共查询到20条相似文献，搜索用时 238 毫秒
1.
Aquifer information carried by aquifer test data may be affected by the presence of a finite thickness skin around the wellbore. The mathematical treatment for an aquifer accounting for the skin zone can be characterized by five parameters, that is, the outer radius of the skin zone and the transmissivity and storativity for each of the skin and aquifer zones. Sensitivity analysis was performed to examine the ground water flow behavior in the skin and aquifer zones in terms of the constant-head test (CHT) data. The simulated annealing procedure was applied to simultaneously determine the skin and aquifer parameters from the analysis of CHT data. Toward the previously mentioned goals, four suites of CHT data were analyzed in this article. The analyses of wellbore flow rate at the test well and the specific drawdown at the observation well gave accurate estimates for the skin and aquifer parameters, respectively. Only the skin thickness and both the skin and the aquifer diffusivities could be accurately estimated from the analysis of drawdown data in the observation well. The estimates for all skin and aquifer parameters from the composite analysis of flow rate and drawdown data were the most accurate. The results of sensitivity analyses and parameter estimations provide instructive references in the analysis of the skin-affected CHT data. 相似文献
2.
Aquifer parameter estimation using an incremental area method 总被引：2，自引：0，他引：2
Theoretical well functions have been derived over the years to predict ground water level behaviour in aquifer systems under stress owing to groundwater extraction. The drawdown data collected during pump tests are typically analysed using graphical curve‐matching procedures to estimate aquifer parameters based on these well functions. Difficulty in aquifer characteristic identification and parameter estimation may arise when the field data do not perfectly match the drawdown curves obtained from the well functions. The present study provides a new method for the interpretation of aquifer pump tests which supplements the existing curve‐matching procedures in case ideal conditions do not exist; the proposed method provides a greater degree of flexibility in the data analysis for diagnostic tool purposes. The method, referred to as the Incremental Area Method (IAM) is based on integrating the logarithmic‐based drawdown curves within a discrete time and matching the results with a corresponding time integral of the Theis ( 1935 ) Well Function which governs ideal confined aquifers. The application of the proposed method to synthetically generated data and field data showed that IAM represents a viable method which yields information on potential non‐idealness of the aquifer and provides aquifer parameter estimates thus potentially overcoming drawdown data curve‐matching difficulties. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
3.
将Kp指数与国际地磁学与高空物理学协会(International Association of Geomagnetism and Aeronomy,下文简称IAGA)推荐的FMI方法(V1.0)计算的K值对比分析后发现,二者的差值在中国区域内出现"不同台站之间K1~K8两头离散、K值整体高出K_p约0.5个单位、个别日期S_R计算错误、相邻日期S_R频繁出现跨天台阶"现象.针对这些现象,分析了可能产生的原因,并对原FMI算法进行了四处改进.用改进后的FMI算法计算了全球8个台站、中国5个台站在2003-2005年和2009年共4年的K指数数据,并将计算结果分别与Kp指数以及这些台站用其他方法计算的K值进行一致性比较,发现改进后的计算结果与二者更为接近,说明本文对FMI方法的改进是有效的. 相似文献
4.
This study presents analytical solutions of the three‐dimensional groundwater flow to a well in leaky confined and leaky water table wedge‐shaped aquifers. Leaky wedge‐shaped aquifers with and without storage in the aquitard are considered, and both transient and steady‐state drawdown solutions are derived. Unlike the previous solutions of the wedge‐shaped aquifers, the leakages from aquitard are considered in these solutions and unlike similar previous work for leaky aquifers, leakage from aquitards and from the water table are treated as the lower and upper boundary conditions. A special form of finite Fourier transforms is used to transform the z‐coordinate in deriving the solutions. The leakage induced by a partially penetrating pumping well in a wedge‐shaped aquifer depends on aquitard hydraulic parameters, the wedge‐shaped aquifer parameters, as well as the pumping well parameters. We calculate lateral boundary dimensionless flux at a representative line and investigate its sensitivity to the aquitard hydraulic parameters. We also investigate the effects of wedge angle, partial penetration, screen location and piezometer location on the steady‐state dimensionless drawdown for different leakage parameters. Results of our study are presented in the form of dimensionless flux‐dimensionless time and dimensionless drawdown‐leakage parameter type curves. The results are useful for evaluating the relative role of lateral wedge boundaries and leakage source on flow in wedge‐shaped aquifers. This is very useful for water management problems and for assessing groundwater pollution. The presented analytical solutions can also be used in parameter identification and in calculating stream depletion rate and volume. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
5.
Data from a large-scale canal-drawdown test were used to estimate the specific yield (sy) of the Biscayne Aquifer, an unconfined limestone aquifer in southeast Florida. The drawdown test involved dropping the water level in a canal by about 30 cm and monitoring the response of hydraulic head in the surrounding aquifer. Specific yield was determined by analyzing data from the unsteady portion of the drawdown test using an analytical stream-aquifer interaction model (Zlotnik and Huang 1999). Specific yield values computed from drawdown at individual piezometers ranged from 0.050 to 0.57, most likely indicating heterogeneity of specific yield within the aquifer (small-scale variation in hydraulic conductivity may also have contributed to the differences in sy among piezometers). A value of 0.15 (our best estimate) was computed based on all drawdown data from all piezometers. We incorporated our best estimate of specific yield into a large-scale two-dimensional numerical MODFLOW-based ground water flow model and made predictions of head during a 183-day period at four wells located 337 to 2546 m from the canal. We found good agreement between observed and predicted heads, indicating our estimate of specific yield is representative of the large portion of the Biscayne Aquifer studied here. This work represents a practical and novel approach to the determination of a key hydrogeological parameter (the storage parameter needed for simulation and calculation of transient unconfined ground water flow), at a large spatial scale (a common scale for water resource modeling), for a highly transmissive limestone aquifer (in which execution of a traditional pump test would be impractical and would likely yield ambiguous results). Accurate estimates of specific yield and other hydrogeological parameters are critical for management of water supply, Everglades environmental restoration, flood control, and other issues related to the ground water hydrology of the Biscayne Aquifer. 相似文献
6.
Abstract Unconfined aquifer parameters, viz. transmissivity, storage coefficient, specific yield and delay index from a pumping test are estimated using the genetic algorithm optimization (GA) technique. The parameter estimation problem is formulated as a least-squares optimization, in which the parameters are optimized by minimizing the deviations between the field-observed and the model-predicted time–drawdown data. Boulton's convolution integral for the determination of drawdown is coupled with the GA optimization technique. The bias induced by three different objective functions: (a) the sum of squares of absolute deviations between the observed and computed drawdown; (b) the sum of squares of normalized deviations with respect to the observed drawdown; and (c) the sum of squares of normalized deviations with respect to the computed drawdown, is statistically analysed. It is observed that, when the time–drawdown data contain no errors, the objective functions do not induce any bias in the parameter estimates and the true parameters are uniquely identified. However, in the presence of noise, these objective functions induce bias in the parameter estimates. For the case considered, defining the objective function as the sum of the squares of absolute deviations between the observed and simulated drawdowns resulted in the best possible estimates. A comparison of the GA technique with the curve-matching procedure and a conventional optimization technique, such as the sequential unconstrained minimization technique (SUMT), is made in estimating the aquifer parameters from a reported field pumping test in an unconfined aquifer. For the case considered, the GA technique performed better than the other two techniques in parameter estimation, with the sum-of-squares errors obtained from the GA about one fourth of those obtained by the curve matching procedure, and about half of those obtained by SUMT. Citation Rajesh, M., Kashyap, D. & Hari Prasad, K. S. (2010) Estimation of unconfined aquifer parameters by genetic algorithms. Hydrol. Sci. J. 55(3), 403–413. 相似文献
7.
A new type of vertical circulation well (VCW) is used for groundwater dewatering at construction sites. This type of VCW consists of an abstraction screen in the upper part and an injection screen in the lower part of a borehole, whereby drawdown is achieved without net withdrawal of groundwater from the aquifer. The objective of this study is to evaluate the operation of such wells including the identification of relevant factors and parameters based on field data of a test site and comprehensive numerical simulations. The numerical model is able to delineate the drawdown of groundwater table, defined as free‐surface, by coupling the arbitrary Lagrangian–Eulerian algorithm with the groundwater flow equation. Model validation is achieved by comparing the field observations with the model results. Eventually, the influences of selected well operation and aquifer parameters on drawdown and on the groundwater flow field are investigated by means of parameter sensitivity analysis. The results show that the drawdown is proportional to the flow rate, inversely proportional to the aquifer conductivity, and almost independent of the aquifer anisotropy in the direct vicinity of the well. The position of the abstraction screen has a stronger effect on drawdown than the position of the injection screen. The streamline pattern depends strongly on the separation length of the screens and on the aquifer anisotropy, but not on the flow rate and the horizontal hydraulic conductivity. 相似文献
8.
In the past, graphical or computer methods were usually employed to determine the aquifer parameters of the observed data obtained from field pumping tests. Since we employed the computer methods to determine the aquifer parameters, an analytical aquifer model was required to estimate the predicted drawdown. Following this, the gradient‐type approach was used to solve the nonlinear least‐squares equations to obtain the aquifer parameters. This paper proposes a novel approach based on a drawdown model and a global optimization method of simulated annealing (SA) or a genetic algorithm (GA) to determine the best‐fit aquifer parameters for leaky aquifer systems. The aquifer parameters obtained from SA and the GA almost agree with those obtained from the extended Kalman filter and gradient‐type method. Moreover, all results indicate that the SA and GA are robust and yield consistent results when dealing with the parameter identification problems. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
9.
The value of subsidence data in ground water model calibration 总被引：2，自引：0，他引：2
The accurate estimation of aquifer parameters such as transmissivity and specific storage is often an important objective during a ground water modeling investigation or aquifer resource evaluation. Parameter estimation is often accomplished with changes in hydraulic head data as the key and most abundant type of observation. The availability and accessibility of global positioning system and interferometric synthetic aperture radar data in heavily pumped alluvial basins can provide important subsidence observations that can greatly aid parameter estimation. The aim of this investigation is to evaluate the value of spatial and temporal subsidence data for automatically estimating parameters with and without observation error using UCODE-2005 and MODFLOW-2000. A synthetic conceptual model (24 separate cases) containing seven transmissivity zones and three zones each for elastic and inelastic skeletal specific storage was used to simulate subsidence and drawdown in an aquifer with variably thick interbeds with delayed drainage. Five pumping wells of variable rates were used to stress the system for up to 15 years. Calibration results indicate that (1) the inverse of the square of the observation values is a reasonable way to weight the observations, (2) spatially abundant subsidence data typically produce superior parameter estimates under constant pumping even with observation error, (3) only a small number of subsidence observations are required to achieve accurate parameter estimates, and (4) for seasonal pumping, accurate parameter estimates for elastic skeletal specific storage values are largely dependent on the quantity of temporal observational data and less on the quantity of available spatial data. 相似文献
10.
本文基于考虑人工场源辐射的可控源地面电磁勘探三维数值模拟技术,以目前常用的标量可控源音频大地电磁勘探测量方式为依据,讨论地面可控源三维电磁勘探中,为达到最佳的探测效果,如何选择合适的反演数据类型.地面可控源数据的三维反演采用有限内存拟牛顿方法.反演过程中,三维可控源频率域响应数值计算采用交错网格有限差分法,求解基于二次电场的Helmholtz方程.发射装置采用长度为1000m的有限长直导线源,测量频率为10Hz;测点个数200个,分布在10条剖面上.在异常体分布区分别观测(1)电场Ex分量的振幅和相位,(2)Ex振幅,(3)Ex与Ey分量的振幅和相位,(4)Ex与Hy的振幅和相位.反演的数据类型分别为上述4种数据以及导出的阻抗Zxy振幅和相位.反演模型由30×30×20个网格组成,测区内水平方向大小为50m×50m,垂直方向厚度为50m,最后5层厚度倍增.反演都从均匀半空间开始,迭代120次结束.数值模拟结果发现,(1)单个电场分量Ex,其相位信息对异常体信息提取非常重要,若只反演该电场振幅,深部电阻率分辨率低;(2)观测正交电场的效果比单个电场分量效果好,其浅部异常的边缘效应明显减弱,深部异常形态完整;(3)同时观测正交的电场和磁场,反演效果与只观测电场振幅和相位的相当;(4)从阻抗数据反演得到的异常位置和电阻率分布均有较大的改善,反演模型分辨率最好.因此,在理论上,建议在野外数据采集时,最好测量正交电磁场;次之,同时测量正交的电场;最次的,可以只观测电场Ex的振幅和相位,亦可取得较足够的信息.但如果只能获得电场的振幅信息,反演结果的深部将出现较大的不确定性. 相似文献
11.
Aquifer Properties Determined from Two Analytical Solutions 总被引：3，自引：0，他引：3
In the analysis of pumping test data, the quality of the determined aquifer parameters can be greatly improved by using a proper model of the aquifer system. Moench (1995) provided an analytical solution for flow to a well partially penetrating an unconfined aquifer. His solution, in contrast to the Neuman solution (1974), accounts for the noninstantaneous decline of the water table (delayed yield). Consequently, the calculated drawdown in these two solutions is different under certain circumstances, and this difference may therefore affect the computation of aquifer properties from pumping test data. This paper uses an inverse computational method to calculate four aquifer parameters as well as a delayed yield parameter, α1 from pumping test data using both the Neuman (1974) and Moench (1995) solutions. Time-drawdown data sets from a pumping test in an unconfined alluvial aquifer near Grand Island, Nebraska, were analyzed. In single-well analyses, horizontal hydraulic conductivity values derived from the Moench solution are lower, but vertical hydraulic conductivity values are higher than those calculated from the Neuman solution. However, the hydraulic conductivity values in composite-well analyses from both solutions become very close. Furthermore, the Neuman solution produces similar hydraulic conductivity values in the single-well and composite-well analyses, but the Moench solution does not. While variable α1, seems to play a role in affecting the computation of aquifer parameters in the single-well analysis, a much smaller effect was observed in the composite-well analysis. In general, specific yield determined using the Moench solution could be slightly higher than the values from the Neuman solution; however, they are still lower than the realistic values for sand and gravel aquifers. 相似文献
12.
Stephen J. Baker 《Ground water》1984,22(5):609-616
Abstract. Two useful programs have been developed for the Hewlett Packard HP41CV programmable calculator. The THEIS program is designed to simulate a well pumping from a confined or unconfined aquifer. Drawdown, residual drawdown, t/t1 and t/r2 are calculated. The BOUN program is designed to solve for drawdown in a well pumping from an aquifer bounded by two parallel impermeable barriers. The programs can be used in aquifer pumping test design, pumping test analysis, and aquifer response predictions. 相似文献
13.
Edwards DA 《Ground water》2012,50(4):554-561
The ability to manipulate analytical expressions for aquifer drawdown can provide insights into groundwater flow processes and assist with assessing strengths and weaknesses of aquifer parameter estimation methods. In the Cooper and Jacob (1946) parameter estimation method, the antilog of the horizontal-axis intercept in a plot of drawdown vs. log(time) is referred to as time naught (t(0)), which is used for estimating storativity. This article briefly reviews traditional uses of the time-naught concept and then spends time introducing new insights and applications involving (1) time-naught/distance relationships, including ways to compensate for certain missing data; (2) use of time naught in a simple method providing a quick visual check of which data in a Cooper-Jacob plot are suitable for use in linear regression; (3) application of time naught, as determined for one well, in estimating the later minimum time for which data from a distant well can be used in the Cooper-Jacob method; (4) development of relationships between drawdown and time naught; (5) use of time naught in a simple algebraic equation to estimate drawdown at smaller times than feasible using the Cooper-Jacob method; and (6) employment of time naught and a vertical-axis intercept on a plot of drawdown vs. log(time) for evaluating storativity. This information may be useful to new hydrogeologists or others interested in further developing their analytical well hydraulics skills. 相似文献
14.
An analytical model of stream/aquifer interaction is proposed that predicts drawdown in an aquifer with leakage from a finite-width stream induced by pumping from a well. The model is formulated based on the assumptions of stream partial penetration, a semipervious streambed, and distributed recharge across a finite-width stream. Advantages of the analytical solution include its simple structure, consisting of the Theis well function with integral modifications. The solution is derived for the semi-infinite domain between the stream and pumping well, which is of primary interest to hydrogeologists. Previous stream/aquifer analytical models are compared to the analytical solution based on dimensionless drawdown profiles. Drawdown in the aquifer near a wide stream was found to be less than that predicted by a solution that ignored stream width. Deviations between the proposed analytical solutions and previous solutions increase as stream width increases. For a hypothetical stream/aquifer system, the proposed analytical solution was equivalent to prior solutions when the ratio of the distance between the stream and aquifer to the stream width was greater than 25. This analytical solution may provide improved estimates of aquifer and streambed leakage parameters by curve fitting experimental field drawdown data. 相似文献
15.
Sang-Heon Shim Thomas S. Duffy Guoyin Shen 《Physics of the Earth and Planetary Interiors》2000,120(4):327-338
Pressure–volume measurements have been performed for CaSiO3 perovskite to 108 GPa at 300 K using NaCl and argon pressure-transmitting media, and energy dispersive X-ray diffraction (EDXD) in a diamond-anvil cell (DAC). By determining a parameter that is the product of the elastic anisotropy, S, and the uniaxial stress component, t, for each data point, we define the stress condition of the sample. For different points at the same pressure in a temperature-quenched sample, the St value can differ by as much as a factor of 5, indicating heterogeneity in the stress condition. This may be responsible for the large scatter of earlier P–V measurements in the DAC which in general used a large diameter X-ray beam. Also, the St value provides insight into the elastic anisotropy, S, of CaSiO3 perovskite and platinum. The sign of S (positive) for CaSiO3 perovskite agrees with first principles calculations but the magnitude may be inconsistent. A new compression curve at 300 K was obtained for CaSiO3 perovskite by using those data points which represent the most nearly hydrostatic conditions. It is observed that the data points with high St values yield larger volumes than the points with small St values at a given pressure. By selecting the data points having low St values (St≤0.005), combining with lower pressure large volume press (LVP) measurements and fitting to third order Birch–Murnaghan equation of state (EOS), we find that CaSiO3 perovskite is more compressible (V0=45.58±0.05 Å3, KT0=236±4 GPa, and KT0′=3.9±0.2 GPa) than suggested by previous studies. The density and bulk modulus of CaSiO3 perovskite at lower mantle pressures and 300 K are 1–3% greater and 5–15% smaller, respectively, than found in previous studies. This study demonstrates that defining the stress state of the sample is crucial to obtain an accurate 300 K compression curve for unquenchable high-pressure phases. 相似文献
16.
《Advances in water resources》2004,27(10):981-990
The Theis equation has been widely used to study the transient movement of groundwater as a result of pumping in a confined aquifer. It is well known that the observed drawdown at early times has an obvious departure from the theoretical drawdown based on the Theis equation. The Theis equation was derived under the assumption that total stress in the aquifer was constant and the mechanical behavior of the confining unit was neglected. However, most geological formations, especially those which are well consolidated, have rigidity and therefore may bend like a plate to a certain extent. The increase in the effective stress in the aquifer due to pumping may not contribute entirely to the compression of the aquifer, but may be partially cancelled out by bending of the overlying aquitard. This means only a part of the total stress is used to compact the aquifer, or the aquifer cannot produce as much water as estimated from the Theis equation. This paper investigated the impact of the bending effect of the confining unit on drawdown. An analytical model which couples flow in the aquifer and bending of the confining unit was presented. The theory is based on elastic plates and solutions were given to the drawdown of groundwater level and deflection of the overlying formation. The drawdown estimated from the new equation was compared with that from the Theis equation. It can be concluded that drawdown from the Theis equation is less than the drawdown predicted by including the bending effect of the confining unit. Both a hypothetical example and a field pumping test in Shandong Province, China, were used to demonstrate the bending effect of the confining unit in the analysis of pumping test data. This paper demonstrated that the initial disagreement between observed drawdown and the Theis solution could be caused by the bending effect of the confining unit, a phenomenon not well addressed in traditional pumping test analysis. A quantitative understanding of this phenomenon can provide improved guidelines for analyzing drawdown data in a confined aquifer. 相似文献
17.
选取易县地震台YRY-4型分量式钻孔应变仪2012—2020年观测数据,进行Venedikov调和分析,并对2015—2020年四分量数据S1、S2、S3、S4进行自洽检验,发现应变观测数据连续性、可靠性较高。选取该台2018—2019年钻孔应变原始数据S1、S2、S3、S4,利用面应变S1 + S3、S2 + S4相关系数,分析河北省3次ML≥4.0地震前数据异常变化,并利用剪应变S1 - S3、S2 - S4的大小,分析受抽水干扰的NE分量主应力方位角,判断干扰源方位,结果发现：①面应变与剪应变在地震发生前均会出现异常变化,且剪应变异常变化幅度较大,应与震中相对于观测井的位置有关;②对于受抽水干扰的钻孔应变NE分量,利用剪应变固体潮数据推算主应力方位角,可得干扰源在钻孔应变观测室NW80.96°方位,与实际位置存在一定偏差。 相似文献
18.
Tomas Perina 《Ground water》2020,58(6):993-999
Hydraulic testing for aquifer characterization at contaminated sites often includes tests of short duration and of different types, such as slug tests and pumping tests, conducted at different phases of investigation. Tests conducted on a well cluster installed in a single aquifer can be combined in aggregate inverse analysis using an analytical model for groundwater flow near a test well. A genetic algorithm performs parallel search of the parameter space and provides starting parameter values for a Markov chain Monte Carlo simulation to estimate the parameter distribution. This sequence of inverse methods avoids guessing of the initial parameter vector and the often encountered difficult convergence of gradient-based methods and estimates the parameter covariance matrix from a distribution rather than from a single point in the parameter space. Combination of different tests improves the resolution of the estimated aquifer properties and allows an assessment of the uniformity of the aquifer. Estimated parameter correlations and standard deviations are used as relative metrics to distinguish well resolved and poorly resolved parameters. The methodology is demonstrated on example field tests in unconfined and leaky aquifers. 相似文献
19.
John J. Hickey 《Ground water》1984,22(5):544-547
The acceptability of the hypothesis of Darcian flow through a semiconfined carbonate aquifer was tested prior to running a multiple-day aquifer test in Pinellas County, Florida. The approach used to test the hypothesis was to run a number of hour-long aquifer tests at different discharges with drawdown measured at the same time during each test in two observation wells, one at 35 feet and the other at 733 feet from the pumped well. If the hypothesis were acceptable, a plot of drawdown versus discharge should describe a straight line. A linear relation was described in 15-minute and 1-hour data plots of drawdown versus discharge for both observation wells. Thus, the hypothesis of Darcian flow through the semiconfined carbonate aquifer was deemed acceptable. Drawdown data from the distant observation well collected during the following multiple-day aquifer test were then analyzed for aquifer-property values of transmissivity, storage coefficient, and leakance coefficient using standard methods. Discharge for the multiple-day test was within the discharge range of the hour-long tests. 相似文献
20.
In this study, we use a linearization procedure and a finite difference method to solve non-Darcian flow to a well in an aquifer–aquitard system. The leakage effect is considered. Flow in the aquifer is assumed to be non-Darcian and horizontal, whereas flow in the aquitard is assumed to be Darcian and vertical. The Izbash equation [Izbash SV. O filtracii V Kropnozernstom Materiale. USSR: Leningrad; 1931 [in Russian]] is employed to describe the non-Darcian flow. The wellbore storage is also considered in this study. An approximate semi-analytical solution has been obtained by the linearization procedure, and a numerical solution has been obtained by using a finite difference method. The previous solutions for Darcian flow case and non-Darcian flow case without leakage can be described as special cases of the new solutions. The error caused by the linearization procedure has also been analyzed. The relative error caused by the linearization procedure is nearly 100% at early times, and decreases to zero at late times. We have also compared the results in this study with Wen et al. [Wen Z, Huang G, Zhan H. A numerical solution for non-Darcian flow to a well in a confined aquifer using the power law function. J Hydrol, 2008d [in revision]] in which the leakage effect is not considered, and Hantush and Jacob [Hantush MS, Jacob CE. Non-steady radial flow in an infinite leaky aquifer. Trans Am Geophys Union 1955;36(1):95–100] who investigated a similar problem in Darcian flow case. The comparison of this study and Wen et al. (2008d) indicates the dimensionless drawdown in the aquifer with leakage is less than that without leakage, and the leakage has little effect at early times. The comparison between the results of this study and that of Hantush and Jacob (1955) indicates that the dimensionless drawdown in the aquifer for non-Darcian flow is larger at early times and smaller at late times, than their counterparts for Darcian flow. A larger dimensionless non-Darcian conductivity kD results in a smaller dimensionless drawdown in the aquifer at late times, and leads to a larger dimensionless drawdown in the aquifer at early times. A smaller dimensionless leakage parameter BD results in a smaller drawdown at late times, and the leakage does not affect the early-time drawdown. The analysis of the dimensionless drawdown inside the well has also been included in this study when the wellbore storage is considered. 相似文献