共查询到20条相似文献，搜索用时 781 毫秒
1.
The impact of ground water pumping on nearby streams is often estimated using analytic models of the interconnected stream-aquifer system. A common assumption of these models is that the pumped aquifer is underlain by an impermeable formation. A new semianalytic solution for drawdown and stream depletion has been developed that does not require this assumption. This solution shows that pumping-induced flow (leakage) through an underlying aquitard can be an important recharge mechanism in many stream-aquifer systems. The relative importance of this source of recharge increases with the distance between the pumping well and the stream. The distance at which leakage becomes the primary component of the pumping-induced recharge depends on the specific properties of the aquifer, aquitard, and streambed. Even when the aquitard is orders of magnitude less transmissive than the aquifer, leakage can be an important recharge mechanism because of the large surface area over which it occurs. Failure to consider aquitard leakage can lead to large overestimations of both the drawdown produced by pumping and the contribution of stream depletion to the pumping-induced recharge. The ramifications for water resources management and water rights adjudication can be significant. A hypothetical example helps illustrate these points and demonstrates that more attention should be given to estimating the properties of aquitards underlying stream-aquifer systems. The solution presented here should serve as a relatively simple but versatile tool for practical assessments of pumping-induced stream-aquifer interactions. However, this solution should not be used for such assessments without site-specific data that indicate pumping has induced leakage through the aquitard. 相似文献
2.
In the analysis of streamflow depletion, the Hunt (1999) solution has an important advantage because it considers a partially penetrating stream. By extending the Hunt drawdown solution, this paper presents semianalytical solutions for gaining streams that evaluate the induced stream infiltration and base flow reduction separately. Simulation results show that for a given deltah (the initial hydraulic head difference between stream and aquifer beneath the channel), the base flow reduction is in direct proportion to the product of streambed leakage (lambda) and the distance between pumping well and stream (L), and the induced stream infiltration is in inverse proportion to lambdaL. Deltah has a significant effect on the ratio of stream infiltration to base flow reduction. The results from the semianalytical solutions agree well with those from MODFLOW simulations. The semianalytical solutions are useful in the verification of numerical simulations and in the analysis of stream-aquifer interactions where water quantity or quality is concerned. 相似文献
3.
Migration of induced-infiltrated stream water into nearby aquifers due to seasonal ground water withdrawal 总被引：1，自引：0，他引：1
Chen X 《Ground water》2001,39(5):721-728
Analysis of stream-aquifer interaction due to ground water extraction has traditionally focused on the determination of the amount of water depleted in the stream. Less attention has been paid to the movement of infiltrated stream water inside aquifer, particularly for agricultural areas. This paper presents a method of using particle-tracking techniques to evaluate the transport of the leaked stream water in the nearby aquifers. Simple stream-aquifer conditions are used to demonstrate the usefulness of the analysis. Travel times, pathlines, and influence zones of stream water were determined between a stream and nearby pumping wells for seasonal ground water extraction areas. When water quantity is a concern, the analyses provide additional information about stream depletion; when water quality is an issue, they offer information for wellhead protection. Analyses were conducted for transient conditions, and both pumping and nonpumping periods were considered. According to the results from the simulation examples, migration of infiltrated stream water into the nearby aquifers is generally slow and most infiltrated stream water does not arrive at the pumping well at the end of a 90-day irrigation season. Infiltrated stream water may remain in the aquifer for several years before arriving at the pumping well. For aquifers with a regional hydraulic gradient toward streams, part of the infiltrated stream water may discharge back to streams during a recovery period. 相似文献
4.
Ground water pumping from aquifer systems that are hydraulically connected to streams depletes streamflow. The amplitude and timing of stream depletion depend on the stream depletion factor (SDF(i)) of the pumping wells, which is a function of aquifer hydraulic characteristics and the distance from the wells to the stream. Wells located at different locations, but having the same SDF and the same rate and schedule of pumping, will deplete streamflow equally. Wells with small SDF(i) deplete streamflow approximately synchronously with pumping. Wells with large SDF(i) deplete streamflow at approximately a constant rate throughout the year, regardless of the pumping schedule. For large values of SDF(i), artificial recharge that occurs on a different schedule from pumping can offset streamflow depletion effectively. The requirements are (1) that the pumping and recharge wells both have the same SDF(i) and (2) that the annual total quantities of recharge and pumping be equal. At larger SDF(i) values, it takes longer for pumping to impact streamflow in a wide aquifer than it does in a narrow aquifer. In basins that are closed to further withdrawals because streamflow is fully allocated, water-use changes replace new allocations as the source of water for new developments. Ground water recharge can be managed to offset the impacts of new ground water developments, allowing for changes in the timing and source of withdrawals from a basin without injuring existing users or instream flows. 相似文献
5.
Drawdown and stream depletion produced by pumping in the vicinity of a partially penetrating stream 总被引：3，自引：0，他引：3
Commonly used analytical approaches for estimation of pumping-induced drawdown and stream depletion are based on a series of idealistic assumptions about the stream-aquifer system. A new solution has been developed for estimation of drawdown and stream depletion under conditions that are more representative of those in natural systems (finite width stream of shallow penetration adjoining an aquifer of limited lateral extent). This solution shows that the conventional assumption of a fully penetrating stream will lead to a significant overestimation of stream depletion (> 100%) in many practical applications. The degree of overestimation will depend on the value of the stream leakance parameter and the distance from the pumping well to the stream. Although leakance will increase with stream width, a very wide stream will not necessarily be well represented by a model of a fully penetrating stream. The impact of lateral boundaries depends upon the distance from the pumping well to the stream and the stream leakance parameter. In most cases, aquifer width must be on the order of hundreds of stream widths before the assumption of a laterally infinite aquifer is appropriate for stream-depletion calculations. An important assumption underlying this solution is that stream-channel penetration is negligible relative to aquifer thickness. However, an approximate extension to the case of nonnegligible penetration provides reasonable results for the range of relative penetrations found in most natural systems (up to 85%). Since this solution allows consideration of a much wider range of conditions than existing analytical approaches, it could prove to be a valuable new tool for water management design and water rights adjudication purposes. 相似文献
6.
Stream water infiltration, bank storage, and storage zone changes due to stream-stage fluctuations 总被引：3，自引：0，他引：3
During a flood period, stream-stage increases induce infiltration of stream water into an aquifer; subsequent declines in stream stage cause a reverse motion of the infiltrated water. This paper presents the results of the water exchange rate between a stream and aquifer, the storage volume of the infiltrated stream water in the surrounding aquifer (bank storage), and the storage zone. The storage zone is the part of aquifer where groundwater is replaced by stream water during the flood. MODFLOW was used to simulate stream–aquifer interactions and to quantify rates of stream infiltration and return flow. MODPATH was used to trace the pathlines of the infiltrated stream water and to determine the size of the storage zone. Simulations were focused on the analyses of the effects of the stream-stage fluctuation, aquifer properties, the hydraulic conductivity of streambed sediments, regional hydraulic gradients, and recharge and evapotranspiration (ET) rates on stream–aquifer interactions. Generally, for a given stream–aquifer system, larger flow rates result from larger stream-stage fluctuations; larger storage volumes and storage zones are produced by larger and longer-lasting fluctuations. For a given stream-stage hydrograph, a lower-permeable streambed, an aquitard, or an anisotropic aquifer of low vertical hydraulic conductivity can significantly reduce the rate of infiltration and limit the size of the storage zone. The bank storage solely caused by the stage fluctuation differs slightly between gaining and losing streams. Short-term rainfall recharge and ET loss in the shallow groundwater slightly influence on the flow rate, but their effects on bank storage in a larger area for a longer period can be considerable. 相似文献
7.
A semianalytical method commonly used for quantifying stream depletion caused by ground water pumping was reviewed for applicability in narrow alluvial aquifers. This stream depletion factor (SDF) method is based on the analytic Glover model, but uses a numerical model-derived input parameter, called the SDF, to partly account for mathematically nonideal conditions such as variable transmissivity and nearby aquifer boundaries. Using the SDF can improve and simplify depletion estimates. However, the method's approximations introduce error that increases with proximity to the impermeable aquifer boundary. This article reviews the history of the method and its assumptions. New stream depletion response curves are presented as functions of well position within bounded aquifers. A simple modification to modeled SDF values is proposed that allows the impermeable boundary to be accounted for with image wells, but without overaccounting for boundary effects that are already reflected in modeled SDFs. It is shown that SDFs for locations closer to the river than to the aquifer boundary do not reflect impermeable-boundary effects, and thus need no modification, and boundary effects in the other portion of the aquifer follow a predictable removable pattern. This method is verified by comparing response curves using modified SDFs with response curves from an extensively calibrated numerical model of a managed ground water recharge site. The modification improves SDF-based stream depletion estimates in bounded aquifers while still benefiting from the additional information contained in SDF maps and retaining their value as standardized references for water rights administration. 相似文献
8.
Stream-aquifer interaction has been the subject of much research for cases of good hydraulic connection (continuous saturated zone) between a river and an aquifer. Under these conditions, floods do not represent a very large net input to the aquifer because most of the water that enters the aquifer during the flood returns to the river when its stage recedes. The situation is different in disconnected stream-aquifer systems, where the streambed lies above the water level in the aquifer, thus preventing return flow from the aquifer. Under these conditions, floods may represent large, but hard to quantify, water inputs. Here, we present a methodology to estimate recharge from floods for disconnected stream-aquifer systems. Recharge is estimated as the product of a flood time function (dependent on the streamflow) and an unknown factor, which is obtained from calibrating a ground water flow model to aquifer heads. The approach can also benefit from concentration data, which can be very informative when river water concentrations vary over time. This methodology is applied to a field situation where recharge from river flooding is found to amount to nearly 15 million m(3)/year on the average, which represents 40% of the total aquifer inputs. Recharge from flooding helps explain major head recoveries, suggesting that basin water management programs should allow some floods to occur. 相似文献
9.
Modeling Stream-Aquifer Interactions Under Seasonal Groundwater Pumping and Managed Aquifer Recharge
In South Korea, a significant amount of groundwater is used for the heating of water-curtain insulated greenhouses during the winter dry season, which had led to problems of groundwater depletion. A managed aquifer recharge (MAR) project is currently underway with the goal of preventing such groundwater depletion in a typical cultivation area, located on an alluvial aquifer near the Nam River. In the present study, FEFLOW, a three-dimensional finite element model, was used to evaluate different strategies for MAR of the cultivation areas. A conceptual model was developed to simulate the stream-aquifer dynamics under the influence of seasonal groundwater pumping and MAR. The optimal rates and duration of MAR were assessed by analyzing the recovery of the groundwater levels and the change in the groundwater temperature. The simulation results indicate that a MAR rate of 8000 m3/d effectively restores the groundwater level when the injection wells are located inside the groundwater depletion area. It is also demonstrated that starting the MAR before the beginning of the seasonal pumping is more effective. Riverbank filtration is preferable for securing the injection water owing to plentiful source of induced recharge from the river. Locating the pumping wells adjacent to the river where there are thick permeable layers could be a good strategy for minimizing decreases in the groundwater level and temperature. 相似文献
10.
C. M. Schmidt A. T. Fisher A. Racz C. G. Wheat M. Los Huertos B. Lockwood 《水文研究》2012,26(15):2235-2247
Artificial recharge of groundwater is an increasingly important method for augmenting groundwater supply and can have a positive or negative influence on the quality of water resources. We instrumented a managed aquifer recharge (MAR) pond in central coastal California to assess how patterns of infiltration and recharge affect the load of nitrate delivered to the underlying aquifer. The concentration of nitrate in infiltrating water consistently decreased during passage through the first metre of subsurface soils. Enrichment of 18O and 15 N in the residual nitrate in infiltrating water proceeded in a ratio of 1:2, indicating that denitrification plays a significant role in the quantitative reduction of nutrients exported during infiltration through shallow soils. The extent and rate of nitrate removal was spatially and temporally variable across the bottom of the recharge pond, with 30% to 60% of the nitrate load being removed over the first 6 weeks of managed aquifer recharge operation. During the period of highest N loading to the system, when the average infiltration rate was > 1 m/day, the recharge pond achieved a load reduction efficiency of 7 kg NO3?‐N/day/ha, which compares favourably to nitrate load reductions achieved by treatment wetlands. Groundwater mounding and water composition below the recharge pond suggest that recharge and subsequent lateral transport occur heterogeneously in the underlying aquifer. Nitrate concentrations in the aquifer following infiltration were lowered primarily by dilution, with little evidence for additional denitrification occurring in the aquifer in comparison to high rates documented during shallow infiltration. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
11.
A study on flood water infiltration and ground water recharge of a shallow alluvial aquifer was conducted in the hyperarid section of the Kuiseb River, Namibia. The study site was selected to represent a typical desert ephemeral river. An instrumental setup allowed, for the first time, continuous monitoring of infiltration during a flood event through the channel bed and the entire vadose zone. The monitoring system included flexible time domain reflectometry probes that were designed to measure the temporal variation in vadose zone water content and instruments to concurrently measure the levels of flood and ground water. A sequence of five individual floods was monitored during the rainy season in early summer 2006. These newly generated data served to elucidate the dynamics of flood water infiltration. Each flood initiated an infiltration event which was expressed in wetting of the vadose zone followed by a measurable rise in the water table. The data enabled a direct calculation of the infiltration fluxes by various independent methods. The floods varied in their stages, peaks, and initial water contents. However, all floods produced very similar flux rates, suggesting that the recharge rates are less affected by the flood stages but rather controlled by flow duration and available aquifer storage under it. Large floods flood the stream channel terraces and promote the larger transmission losses. These, however, make only a negligible contribution to the recharge of the ground water. It is the flood duration within the active streambed, which may increase with flood magnitude that is important to the recharge process. 相似文献
12.
Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding 总被引：2，自引：0，他引：2
Marc Descloitres Laurent Ruiz M. Sekhar Anatoly Legchenko Jean‐Jacques Braun M. S. Mohan Kumar S. Subramanian 《水文研究》2008,22(3):384-394
A groundwater recharge process of heterogeneous hard rock aquifer in the Moole Hole experimental watershed, south India, is being studied to understand the groundwater flow behaviour. Significant seasonal variations in groundwater level are observed in boreholes located at the outlet area indicating that the recharge process is probably taking place below intermittent streams. In order to localize groundwater recharge zones and to optimize implementation of boreholes, a geophysical survey was carried out during and after the 2004 monsoon across the outlet zone. Magnetic resonance soundings (MRS) have been performed to characterize the aquifer and measure groundwater level depletion. The results of MRS are consistent with the observation in boreholes, but it suffers from degraded lateral resolution. A better resolution of the regolith/bedrock interface is achieved using electrical resistivity tomography (ERT). ERT results are confirmed by resistivity logging in the boreholes. ERT surveys have been carried out twice—before and during the monsoon—across the stream area. The major feature of recharge is revealed below the stream with a decrease by 80% of the calculated resistivity. The time‐lapse ERT also shows unexpected variations at a depth of 20 m below the slopes that could have been interpreted as a consequence of a deep seasonal water flow. However, in this area time‐lapse ERT does not match with borehole data. Numerical modelling shows that in the presence of a shallow water infiltration, an inversion artefact may take place thus limiting the reliability of time‐lapse ERT. A combination of ERT with MRS provides valuable information on structure and aquifer properties respectively, giving a clue for a conceptual model of the recharge process: infiltration takes place in the conductive fractured‐fissured part of the bedrock underlying the stream and clayey material present on both sides slows down its lateral dissipation. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
13.
In Belgium, IWVA uses managed aquifer recharge (MAR) to recharge the aquifer with treated wastewater generated from the communities to sustain the potable water supply on the Belgian coast. This MAR facility is faced with a challenge of reduced infiltration rates during the winter season when pond water temperatures near 4°C. This study involves the identification of the predominant factor influencing the rate of infiltration through the pond bed. Several factors, including pumping rates, natural recharge, tidal influences of the North Sea and pond-water temperature, were identified as potential causes for variation of the recharge rate. Correlation statistics and linear regression analysis were used to determine the sensitivity of the infiltration rate to the aforementioned factors. Two groundwater flow models were developed in visual MODFLOW to simulate the water movement under the pond bed and to obtain the differences in flux to track the effects of variation of hydraulic conductivity during the two seasons. A 32% reduction in vertical hydraulic gradient in the top portion of the aquifer was observed in winter, causing the recharge rates to fluctuate. Results showed that water temperature caused a 30% increase in hydraulic conductivity in summer as compared with winter and has the maximum impact on infiltration rate. Cyclic variations in water viscosity, occurring because of seasonal temperature changes, influence the saturated hydraulic conductivity of the pond bed. Results from the models confirm the impact on infiltration rate by temperature-influenced hydraulic conductivity. 相似文献
14.
William R. Wise Chi-Chung Chang Rick A. Klopp Philip B. Bedient 《Ground Water Monitoring & Remediation》1991,11(2):93-100
At an aviation gasoline spill site in Traverse City, Michigan, historical records indicate a positive correlation between significant rainfall events and increased concentrations of slightly soluble organic compounds in the monitoring wells of the site. To investigate the recharge effect on ground water quality due to infiltrating, water percolating past residual oil and into the saturated zone, an in situ infiltration experiment was performed at the site. Sampling cones were set at various depths below a circular test area, 13 feet (4 meters) in diameter. Rainfall was simulated by sprinkling the test area at a rate sufficiently low to prevent runoff. The sampling cones for soil-gas and ground water quality were installed in the unsaturated and saturated zones to observe the effects of the recharge process. At the time of the test, the water table was below the residual oil layer. The responses of the soil-gas and ground water quality were monitored during the recharge and drainage periods, which resulted from the sprinkling.
Infiltrated water was determined to have transported organic constituents of the residual oil, specifically benzene, toluene, ethylbenzene, and ortho-xylene (BTEX), into the ground water beneath the water table, elevating the aqueous concentrations of these constituents in the saturated zone. Soil-gas concentrations of the organic compounds in the unsaturated zone increased with depth and time after the commencement of infiltration. Reaeration of the unconfined aquifer via the infiltrated water was observed. It is concluded that water quality measurements are directly coupled to recharge events for the sandy type of aquifer with an overlying oil phase, which was studied in this work. Ground water sampling strategies and data analysis need to reflect the effect of recharge from precipitation on shallow, unconfined aquifers where an oil phase may be present. 相似文献
Infiltrated water was determined to have transported organic constituents of the residual oil, specifically benzene, toluene, ethylbenzene, and ortho-xylene (BTEX), into the ground water beneath the water table, elevating the aqueous concentrations of these constituents in the saturated zone. Soil-gas concentrations of the organic compounds in the unsaturated zone increased with depth and time after the commencement of infiltration. Reaeration of the unconfined aquifer via the infiltrated water was observed. It is concluded that water quality measurements are directly coupled to recharge events for the sandy type of aquifer with an overlying oil phase, which was studied in this work. Ground water sampling strategies and data analysis need to reflect the effect of recharge from precipitation on shallow, unconfined aquifers where an oil phase may be present. 相似文献
15.
Pumping test evaluation of stream depletion parameters 总被引：1，自引：0，他引：1
16.
Analysis of the vapor in passive vapor samplers retrieved from a streambed in fractured rock terrain implied that volatile organic carbon (VOC) discharge from ground water to surface water substantially increased following installation of a contaminant recovery well using air rotary drilling. The air rotary technique forced air into the aquifer near the stream. The injection produced an upward hydraulic gradient that appears to have transported water and contaminants from deeper parts of the aquifer through fractures into shallow parts of the aquifer. Once in the shallow flow regime, the contamination was transported to the stream, where it discharged during the next several weeks following well installation. After the recovery well was activated and began continuously pumping contaminated ground water to a treatment facility, the VOC concentrations in the stream bottom passive vapor samplers decreased to below detectable concentrations, suggesting that the withdrawal had captured the contaminated ground water that previously had discharged to the stream. 相似文献
17.
Bredehoeft J 《Ground water》2011,49(4):468-475
An aquifer, in a stream/aquifer system, acts as a storage reservoir for groundwater. Groundwater pumping creates stream depletion that recharges the aquifer. As wells in the aquifer are moved away from the stream, the aquifer acts to filter out annual fluctuations in pumping; with distance the stream depletion tends to become equal to the total pumping averaged as an annual rate, with only a small fluctuation. This is true for both a single well and an ensemble of wells. A typical growing season in much of the western United States is 3 to 4 months. An ensemble of irrigation wells spread more or less uniformly across an aquifer several miles wide, pumping during the growing season, will deplete the stream by approximately one-third of the total amount of water pumped during the growing season. The remaining two-thirds of stream depletion occurs outside the growing season. Furthermore, it takes more than a decade of pumping for an ensemble of wells to reach a steady-state condition in which the impact on the stream is the same in succeeding years. After a decade or more of pumping, the depletion is nearly constant through the year, with only a small seasonal fluctuation: ±10%. Conversely, stream depletion following shutting down the pumping from an ensemble of wells takes more than a decade to fully recover from the prior pumping. Effectively managing a conjunctive groundwater and surface water system requires integrating the entire system into a single management institution with a long-term outlook. 相似文献
18.
In this study, we use borehole temperature data and stable isotopes to delineate the flow system and estimate the effect of urbanization in the Nagaoka area of Japan. Temperature profiles were measured four times in observation wells during the period 2000-2001 and compared with those measured in the same wells during the period 1977-1983 (Taniguchi 1986). Water was sampled in both observation and pumping wells during the same period. The temporal and spatial variability in temperature indicate clearly the effect of urban warming and heavy pumping on the ground water system. Urban warming caused higher temperatures recently as compared to the older values, and pumping caused induced recharge from the river to the ground water. The stable isotope data show the ground water flow system is divided into shallow, intermediate, and deep systems, and that land use and infiltration rate are affecting the shallow flow system. 相似文献
19.
《Advances in water resources》2007,30(4):1016-1026
We have proved that the Hantush’s model [Hantush MS. Wells near streams with semipervious beds. J Geophys Res 1965;70:2829–38] in a half-domain can be extended to a whole-domain and becomes identical to that of Hunt [Hunt B. Unsteady stream depletion from ground water pumping. Ground Water 1999;37(1):98–102] for a shallow and infinitely narrow stream, provided that the Dupuit assumption is adopted. This proof helps correct a false concept that regards the Hantush’s model as less useful because of its fully penetrating stream assumption. This study deals with interaction of an aquifer with two parallel streams based on the Hantush’s model. Semi-analytical solutions are obtained based on rigorous mass conservation requirement by maintaining continuity of flux and head along the aquifer–streambed boundaries. This study shows that the hydraulic conductivity ratio of the two streambeds appears to be the most important factor controlling the stream–aquifer interaction, followed by a less important role played by the thickness ratio of the two streambeds. When the low-permeability streambeds do not exist, the steady-state stream depletion from one stream is linearly proportional to the ratio of the shortest distance from the pumping well to the other stream over the shortest distance between the two streams. When the low-permeability streambeds are presented, similar conclusion can be drawn except that the stream depletion now also strongly depends on the hydraulic conductivity ratio of the two streambeds. When the values of the hydraulic conductivity of the two streambeds are different by an order of magnitude, the location of the pumping well that receives equal flux from two streams can be off the middle-line between the two streams by nearly 90%. 相似文献
20.
Little attention has been given to the role of groundwater in the hydrological cycle of lowland watersheds. Our objective in this study was to estimate total recharge to groundwater by analysing water table response to storm events and the rate at which water was transferred into the shallow aquifer. This was conducted at three sites in a rural watershed in the lower Atlantic coastal plain near Charleston, South Carolina, USA. A novel version of the water table fluctuation method was used to estimate total recharge to the shallow aquifer by comparing hourly data of water table position following storm events and measuring water table recession behavior, rather than subjective graphical analysis methods. Also, shallow aquifer recharge rates (vertical fluxes) were estimated using Darcy's Law by comparing static water levels in a water table well and in a shallow piezometer during dry periods. The total annual recharge estimated ranged from 107 ± 39 mm·yr–1 (5–10% of annual precipitation) at a poorly drained topographic low area to 1140 ± 230 mm·yr–1 (62–94% of annual precipitation) for a moderately well‐drained upland site. The average aquifer recharge rate was 114 ± 60 mm·yr–1, which is similar to previous estimations of base flow for the ephemeral third‐order streams in this watershed. The difference in the two methods may have been caused by processes not accounted for in the Darcy flux method, soil moisture deficits, and average evapotranspiration demand, which is about 1100 mm·yr–1 for this region. Although other factors also can affect partitioning of recharge, an integrated approach to inspecting easily gathered groundwater data can provide information on an often neglected aspect of water budget estimation. We also discuss the effects of land use change on recharge reduction, given a typical development scenario for the region. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献