共查询到20条相似文献，搜索用时 953 毫秒
1.
Water levels and water quality of open borehole wells in fractured bedrock are flow-weighted averages that are a function of the hydraulic heads and transmissivities of water contributing fractures, properties that are rarely known. Without such knowledge using water levels and water quality data from fractured bedrock wells to assess groundwater flow and contaminant conditions can be highly misleading. This study demonstrates a cost-effective single packer method to determine the hydraulic heads and transmissivities of water contributing fracture zones in crystalline bedrock wells. The method entails inflating a pipe plug to isolate sections of an open borehole at different depths and monitoring changes in the water level with time. At each depth, the change in water level with time was used to determine the sum of fracture transmissivities above the packer and then to solve for individual fracture transmissivity. Steady-state wellbore heads along with the transmissivities were used to determine individual fracture heads using the weighted average head equation. The method was tested in five wells in crystalline bedrock located at the University of Connecticut in Storrs. The single packer head and transmissivity results were found to agree closely with those determined using conventional logging methods and the dissolved oxygen alteration method. The method appears to be a simple and cost-effective alternative in obtaining important information on flow conditions in fractured crystalline bedrock wells. 相似文献
2.
Fransson A 《Ground water》2002,40(2):201-204
The transmissivities of individual fractures along a borehole are difficult to obtain unless each fracture is tested. To estimate a fracture-transmissivity distribution from section transmissivities, a method was developed based on fixed-interval-length transmissivities and the corresponding number of fractures for each interval. The method is nonparametric and iterative, and the fractures are viewed as two-dimensional features, in which the total transmissivity of a borehole is equal to the sum of individual fracture transmissivities. Initially, a linear a priori assumption of the transmissivity distribution is made, and from this a so-called mean transmissivity function is derived. Subsequently, the mean transmissivity of the Nj fractures within a section, j, of the borehole is estimated, and the same value of the mean transmissivity function represents Nj possible fracture transmissivities from the initial distribution. This is repeated for each borehole section, and, eventually, all fracture transmissivities are sorted to give the next iteration's transmissivity distribution and the corresponding mean transmissivity function. Finally, the distributions converge, yielding a possible fracture-transmissivity distribution. The method was verified for a synthetic data sample and then tested on a sample from a borehole at the Asp? Hard Rock Laboratory, Sweden. For the synthetic data, the method gave a distribution that was fairly close to the original one; for the Asp? data, 15% of the fractures had a transmissivity larger than the measurement limit (1 x 10(-9) m2/sec), and these transmissivities follow a log-normal distribution. 相似文献
3.
Magnetic Resonance Sounding (MRS) is nowadays accepted as a new geophysical method that can be used for a reliable determination of the ground water content distribution in the top 150 m. A great effort has also been made in MRS development to deduce the hydraulic transmissivity, based on empiric relationships of the permeability with a factor F which is calculated with NMR parameters measured at laboratory scale. To use this relationship under field conditions a calibration coefficient CT = Tpt / F has to be previously established, which demands the knowledge of the transmissivity Tpt evaluated in the pumping test. The transmissivity can then be calculated at any other site of the same aquifer using the relation Tmrs = CTF. The CT values reported suggest a certain relationship with the lithology, but with a great dispersion and contradictory results. MRS surveys carried out in alluvial aquifers in Spain have shown that the value of CT evaluated at one site may not be valid at another place of the same aquifer, because of the great heterogeneity of this kind of geological environment. The demand of a pumping test at each site where a MRS is measured invalidates the method actually used for MRS transmissivity evaluation. More than 50 MRS have been used to propose a new methodology. The aquifers visited cover a great range of transmissivities (from 2 × 10− 6 to 9 × 10− 3 m2/s). The MRS signal amplitude varies between 20 and 1400 nV, the signal/noise ratio is in the range from 0.6 to 42, and the value of the decay time constant varies from 200 to 800 ms. It has been demonstrated that when the transmissivity increases, the value of F decreases, and CT increases, except for certain groups of MRS taken at the same aquifer or part of one aquifer, for which F increases with Tpt, keeping CT constant. A function CT(F) of the type CT = mF− n has been obtained that allows the transmissivity evaluation without the need of Tpt. Considering that both values of transmissivity, Tpt and Tmrs, are subjected to deviations due to the experimental errors as well as due to evaluation errors, the prediction achieved by the proposed equation is rather good. To perform a better evaluation of the values of the coefficients m and n it is necessary to have a greater number of MR soundings of good quality and with a trustworthy inversion at locations where a really comparable and good performed pumping test is available, covering a sufficient range of transmissivities. Though the data we have used do not always fulfil these conditions, the result is promising. Once a trustable function is available, the forecast of the transmissivity using MRS will not need the existence of any pumping test in the area. The general extension of this methodology demands the availability of MRS taken at all kinds of geological and hydrogeological environments, which is impossible without the existence of a universal MRS data base. 相似文献
4.
《Journal of Applied Geophysics》2009,67(3-4):128-139
Magnetic Resonance Sounding (MRS) is nowadays accepted as a new geophysical method that can be used for a reliable determination of the ground water content distribution in the top 150 m. A great effort has also been made in MRS development to deduce the hydraulic transmissivity, based on empiric relationships of the permeability with a factor F which is calculated with NMR parameters measured at laboratory scale. To use this relationship under field conditions a calibration coefficient CT = Tpt / F has to be previously established, which demands the knowledge of the transmissivity Tpt evaluated in the pumping test. The transmissivity can then be calculated at any other site of the same aquifer using the relation Tmrs = CTF. The CT values reported suggest a certain relationship with the lithology, but with a great dispersion and contradictory results. MRS surveys carried out in alluvial aquifers in Spain have shown that the value of CT evaluated at one site may not be valid at another place of the same aquifer, because of the great heterogeneity of this kind of geological environment. The demand of a pumping test at each site where a MRS is measured invalidates the method actually used for MRS transmissivity evaluation. More than 50 MRS have been used to propose a new methodology. The aquifers visited cover a great range of transmissivities (from 2 × 10− 6 to 9 × 10− 3 m2/s). The MRS signal amplitude varies between 20 and 1400 nV, the signal/noise ratio is in the range from 0.6 to 42, and the value of the decay time constant varies from 200 to 800 ms. It has been demonstrated that when the transmissivity increases, the value of F decreases, and CT increases, except for certain groups of MRS taken at the same aquifer or part of one aquifer, for which F increases with Tpt, keeping CT constant. A function CT(F) of the type CT = mF− n has been obtained that allows the transmissivity evaluation without the need of Tpt. Considering that both values of transmissivity, Tpt and Tmrs, are subjected to deviations due to the experimental errors as well as due to evaluation errors, the prediction achieved by the proposed equation is rather good. To perform a better evaluation of the values of the coefficients m and n it is necessary to have a greater number of MR soundings of good quality and with a trustworthy inversion at locations where a really comparable and good performed pumping test is available, covering a sufficient range of transmissivities. Though the data we have used do not always fulfil these conditions, the result is promising. Once a trustable function is available, the forecast of the transmissivity using MRS will not need the existence of any pumping test in the area. The general extension of this methodology demands the availability of MRS taken at all kinds of geological and hydrogeological environments, which is impossible without the existence of a universal MRS data base. 相似文献
5.
Abstract In order to calculate the transmissivity from the inverse problem corresponding to the groundwater flow in an isotropic horizontal aquifer, a numerical conservative approach is tested. The method deals with triangulation of the domain and applies the conservation of mass to elements of the mesh using the harmonic mean for internodal transmissivities. An optimal sweeping algorithm is used to evaluate nodal transmissivities from one element to another with a minimal relative error accumulation. The practical importance of the method is demonstrated through two synthetic examples representing those experienced in the field, then through application to a Moroccan aquifer. The computed hydraulic head is well fitted to the reference one, which confirms the validity of the identified transmissivity model. 相似文献
6.
R. T. Miller 《Ground water》1984,22(5):532-537
The U.S. Geological Survey is studying the potential for storage of heated water in a sandstone aquifer in St. Paul, Minnesota. The efficiency of the aquifer to store thermal energy is related, in part, to the hydrogeologic characteristics of the aquifer. The movement of heat away from the injection well is directly related to the anisotropy. Aquifer tests were conducted to determine the anisotropy near the heated-water injection well. The maximum and minimum values of transmissivity along the principal directions of the hydraulic conductivity tensors of the Ironton and Galesville Sandstones in St. Paul, Minnesota are approximately 1,090 and 480 feet squared per day. The storage coefficient is 4.5 × 10−5. These values represent the average of four determinations of nonsteady flow to a well in an idealized infinite anisotropic aquifer. Analysis of the values of transmissivity and storage coefficient for hypothetical changes in location of two of the monitoring wells where depth-deviation surveys were not available indicates that computed transmissivities vary less than 5 percent and storage coefficients vary less than ±6 percent. 相似文献
7.
Joseph G. Meert Douglas L. Smith Len Fishkin 《Journal of Volcanology and Geothermal Research》1991,47(3-4)
Heat flow values were calculated from direct measurements of temperature and thermal conductivity at thirteen sites in the Arkansas-Missouri Ozark Plateau region. These thirteen values are augmented by 101 estimates of heat flow, based on thermal conductivity measurements and temperature gradients extrapolated from bottom-hole temperatures. The regional heat flow profile ranges from 9 mW m−2 to over 80 mW m−2, but at least two distinct thermal regimes have been identified. Seven new heat flow determinations are combined with three previously published values for the St. Francois Mountains (SFM), a Precambrian exposure of granitic and rhyolitic basement rocks, average 47 mW m−2. Radioactive heat production of 76 samples of the exposed rocks in the SFM averages 2.4 μW m−2 and a typical continental basement contribution of 14 mW m−2 is implied. Conversely, the sedimentary rock sequence of the plateau is characterized by an anomalously low heat flow, averaging approximately 27 mW m−2. Groundwater transmissivity values that are based on data from 153 wells in deep regional aquifers demonstrate an inverse relationship to the observed heat flow patterns. The areas of high transmissivity that correspond to areas of low total heat flux suggest that the non-conservative vertical heat flow within the Ozark sedimentary sequence can be attributed to the effects of groundwater flow. 相似文献
8.
Interpretation of single-well tests with the Cooper-Jacob method remains more reasonable than most alternatives. Drawdowns from 628 simulated single-well tests where transmissivity was specified were interpreted with the Cooper-Jacob straight-line method to estimate transmissivity. Error and bias as a function of vertical anisotropy, partial penetration, specific yield, and interpretive technique were investigated for transmissivities that ranged from 10 to 10,000 m(2)/d. Cooper-Jacob transmissivity estimates in confined aquifers were affected minimally by partial penetration, vertical anisotropy, or analyst. Cooper-Jacob transmissivity estimates of simulated unconfined aquifers averaged twice the known values. Transmissivity estimates of unconfined aquifers were not improved by interpreting results with an unconfined aquifer solution. Judicious interpretation of late-time data consistently improved estimates where transmissivity exceeded 250 m(2)/d in unconfined aquifers. 相似文献
9.
Assessment of hydrogeological conditions in basement aquifers of the Precambrian Oban massif, southeastern Nigeria 总被引：1，自引：0，他引：1
A combined analysis of lineament length density from radar imagery and surface resistivity data is used to assess the hydrogeological conditions in the Oban massif, Nigeria. The results show that the data guided the qualitative and quantitative estimation of some aquifer parameters. These include resistivity of the water bearing formations (280–740 ω m), thickness (5–140 m), limited hydraulic conductivity (8.53-13.18 m/day) and transmissivity (410.65–725.88 m2/day) data. In addition, the lineament length density for the area ranged between less than 0.2 to slightly more than 0.4. Site evaluation for the location of productive boreholes/wells using a groundwater potential index (GWPI) indicates that areas with a GWPI of greater than 35 are consistent with relatively high yield. 相似文献
10.
Jamal Asfahani 《水文科学杂志》2013,58(2):468-482
Abstract Traditional and polar vertical electrical sounding (VES) techniques were used for computing and characterizing the transmissivity of the Quaternary and Palaeogene aquifers in the Khanasser Valley, northern Syria. The VES technique with its different applied approaches is inexpensive, easy and gives rapid results with high spatial resolution. The Dar-Zarrouk parameters were determined to estimate the groundwater potential of the aquifers. The results show the mean transmissivity of the Quaternary aquifer to be 49 m2/d (minimum: 10 m2/d; maximum: 131 m2/d; standard deviation: 32 m2/d), and that of the Palaeogene aquifer is 0.94 m2/d (minimum: 0.31 m2/d; maximum: 1.54 m2/d; standard deviation: 0.33 m2/d). Editor D. Koutsoyiannis Citation Asfahani, J., 2013. Groundwater potential estimation using vertical electrical sounding measurements in the semi-arid Khanasser Valley region, Syria. Hydrological Sciences Journal, 58 (2), 468--482. 相似文献
11.
Inverse modeling is a useful tool in ground water flow modeling studies. The most frequent difficulties encountered when using this technique are the lack of conditioning information (e.g., heads and transmissivities), the uncertainty in available data, and the nonuniqueness of the solution. These problems can be addressed and quantified through a stochastic Monte Carlo approach. The aim of this work was to compare the applicability of two stochastic inverse modeling approaches in a field-scale application. The multi-scaling (MS) approach uses a downscaling parameterization procedure that is not based on geostatistics. The pilot point (PP) approach uses geostatistical random fields as initial transmissivity values and an experimental variogram to condition the calibration. The studied area (375 km2) is part of a regional aquifer, northwest of Montreal in the St. Lawrence lowlands (southern Québec). It is located in limestone, dolomite, and sandstone formations, and is mostly a fractured porous medium. The MS approach generated small errors on heads, but the calibrated transmissivity fields did not reproduce the variogram of observed transmissivities. The PP approach generated larger errors on heads but better reproduced the spatial structure of observed transmissivities. The PP approach was also less sensitive to uncertainty in head measurements. If reliable heads are available but no transmissivities are measured, the MS approach provides useful results. If reliable transmissivities with a well inferred spatial structure are available, then the PP approach is a better alternative. This approach however must be used with caution if measured transmissivities are not reliable. 相似文献
12.
Andrew C. Mills 《Ground water》2010,48(4):604-608
Parameter-estimation methods, including an exhaustive-search method and PEST (Parameter ESTimation) software, were applied to recovery-test data and slug-test data to obtain best estimates of transmissivity (T) by minimizing the sums of residuals. Each residual represents the difference between the field-measured water-level value and the value calculated by the appropriate non-linear equation. The exhaustive-search method in both cases involves computing the sums of residuals for an array of transmissivity and storativity values selected by the user for testing. Two new Fortran programs are presented that employ the exhaustive-search method. They utilize Picking's method for analyzing recovery-test data and the analytical equation for analyzing slug-test data derived by Cooper, Bredehoeft, and Papadopulos. Picking's method involves application of the Papadopulos and Cooper's equation for drawdown in finite-diameter wells. Utilizing field data reported in the literature, the estimated transmissivity values from the exhaustive-search methods were compared to the literature values obtained by type-curve matching techniques. The exhaustive-search values corresponded closely to the curve-matching values. Estimates for T were also obtained from recovery-test and slug-test data from two sites in southeastern Pennsylvania. For these sites, the PEST program was also applied to the data to evaluate the accuracy of the exhaustive-search methods. The results from the two methods were generally in good agreement. The two new Fortran programs are practical tools for the hydrogeologist, as they require less time compared to type-curve matching and the PEST method, and they yield accurate estimates of transmissivity. 相似文献
13.
Violaine Ponsin Amélie Chablais Julien Dumont Olivier Radakovitch Patrick Höhener 《Ground Water Monitoring & Remediation》2015,35(2):30-38
The objective of this study was to investigate whether 222Rn in groundwater can be used as a tracer for light non‐aqueous phase liquid (LNAPL) quantification at a field site treated by dual‐phase LNAPL removal. After the break of a pipeline, 5 ha of soil in the nature reserve Coussouls de Crau in southern France was contaminated by 5100 m3 of crude oil. Part of this oil seeped into the underlying gravel aquifer and formed a floating oil body of about 3.9 ha. The remediation consists of plume management by hydraulic groundwater barriers and LNAPL extraction in the source zone. 222Rn measurements were performed in 21 wells in and outside the source zone during 15 months. In uncontaminated groundwater, the radon activity was relatively constant and remained always >11 Bq/L. The variability of radon activity measurements in wells affected by the pump‐and‐skim system was consistent with the measurements in wells that were not impacted by the system. The mean activities in wells in the source zone were, in general, significantly lower than in wells upgradient of the source zone, owing to partitioning of 222Rn into the oil phase. The lowest activities were found in zones with high non‐aqueous phase liquid (NAPL) recovery. LNAPL saturations around each recovery well were furthermore calculated during a period of high groundwater level, using a laboratory‐determined crude oil–water partitioning coefficient of 38.5 ± 2.9. This yielded an estimated volume of residual crude oil of 309 ± 93 m3 below the capillary fringe. We find that 222Rn is a useful and cheap groundwater tracer for finding zones of good LNAPL recovery in an aquifer treated by dual‐phase LNAPL removal, but that quantification of NAPL saturation using Rn is highly uncertain. 相似文献
14.
Ground water reservoirs in the Choshuichi alluvial fan, central western Taiwan, were investigated using direct-current (DC) resistivity soundings at 190 locations, combined with hydrogeological measurements from 37 wells. In addition, attempts were made to calculate aquifer transmissivity from both surface DC resistivity measurements and geostatistically derived predictions of aquifer properties. DC resistivity sounding data are highly correlated to the hydraulic parameters in the Choshuichi alluvial fan. By estimating the spatial distribution of hydraulic conductivity from the kriged well data and the cokriged thickness of the correlative aquifer from both resistivity sounding data and well information, the transmissivity of the aquifer at each location can be obtained from the product of kriged hydraulic conductivity and computed thickness of the geoelectric layer. Thus, the spatial variation of the transmissivities in the study area is obtained. Our work is more comparable to Ahmed et al. (1988) than to the work of Niwas and Singhal (1981). The first "constraint" from Niwas and Singhal's work is a result of their use of linear regression. The geostatistical approach taken here (and by Ahmed et al. [1988]) is a natural improvement on the linear regression approach. 相似文献
15.
John P. McDonald 《Ground water》2021,59(5):694-709
High transmissivity aquifers typically have low hydraulic gradients (i.e., a flat water table). Measuring low gradients using water levels can be problematic because measurement error may be greater than the true difference in water levels (i.e., a low signal-to-noise ratio). In this study, the feasibility of measuring a hydraulic gradient in the range of 10−6 to 10−5 m/m was demonstrated. The study was performed at a site where the depth to water from land surface ranged from 40.1 to 94.2 m and the aquifer transmissivity was estimated at 41,300 m2/d (hydraulic conductivity of 18,800 m/d). The goals of the study were to reduce measurement error as much as practicable and assess the importance of factors affecting water level measurement accuracy. Well verticality was the largest source of error (0.000 to 0.168 m; median of 0.014 m), and geodetic survey of casing elevations was the next most important source of error (0.002 to 0.013 m; median of 0.005 m). Variability due to barometric pressure fluctuations was not an important factor at the site. Hydraulic heads were measured to an accuracy of ±0.0065 m, and the average hydraulic gradient was estimated to be 8.0 × 10−6 (±0.9 × 10−6) m/m. The improvement in accuracy allowed for two reversals in the groundwater flow direction to be identified, after which the gradient averaged 2.5 × 10−5 (±0.4 × 10−5) m/m. This study showed it is possible to sufficiently control sources of error to measure hydraulic gradients in the 10−6 to 10−5 m/m range. 相似文献
16.
Hydraulic fracturing has become an important technique for enhancing the permeability of hydrocarbon source rocks and increasing aquifer transmissivity in many hard rock environments where natural fractures are limited, yet little is known about the nature or behaviour of these hydraulically induced fractures as conduits to flow and transport. We propose that these fractures tend to be smooth based on observed hydraulic and transport behaviour. In this investigation a multi‐faceted approach was used to quantify the properties and characteristics of an isolated hydraulically induced fracture in crystalline rocks. Packers were used to isolate the fracture that is penetrated by two separate observation wells located approximately 33 m apart. A series of aquifer tests and an induced gradient tracer test were performed to better understand the nature of this fracture. Aquifer test results indicate that full recovery is slow because of the overall low permeability of the crystalline rocks. Drawdown tests indicate that the fracture has a transmissivity of 1–2 m2/day and a specific storage on the order of 2–9 × 10?7/m. Analysis of a potassium–bromide tracer test break through curve shows classic Fickian behaviour with minimal tailing analogous to parallel plate flow. Virtually all of the tracer was recovered, and the breakthrough curve dilution indicates that the swept area is only about 11% of a radial flow field and the estimated aperture is ≤0.5 mm, which implies a narrow linear flow region. These outcomes suggest that transport within these hydraulically induced ‘smooth’ fractures in crystalline rocks is rapid with minimal mixing, small local velocity fluctuations and no apparent diffusion into the host rock or secondary fractures. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
17.
Klusman K 《Ground water》2004,42(6-7):945-948
A new computer program has been developed to automate analysis of brief single-well pumping tests. Adapted from a procedure developed by Picking (1994) that does not require measurement of the pumping rate, this new program is menu-driven and eliminates one significant source of imprecision in Picking's original method, namely, selection of "well function of u" values by interpolation in a lookup table. This new program has been applied to tests of 25 domestic wells penetrating bedrock, each pumped for <2 min. 相似文献
18.
Evaluating an Analytical Model to Predict Subsurface LNAPL Distributions and Transmissivity from Current and Historic Fluid Levels in Groundwater Wells: Comparing Results to Numerical Simulations 
下载免费PDF全文
下载免费PDF全文 Robert J. Lenhard Kaveh Sookhak Lari John L. Rayner Greg B. Davis 《Ground Water Monitoring & Remediation》2018,38(1):75-84
A recent analytical model predicts free, entrapped, and residual LNAPL saturations and the LNAPL transmissivity in the subsurface from current and historic fluid levels in groundwater wells. As such, the model accounts for effects of fluid level fluctuations in a well. The model was developed to predict LNAPL specific volumes and transmissivities from current fluid level measurements in wells and either recorded historic fluid level fluctuations in wells or estimates. An assumption is made in the model that the predictions are not dependent on whether the historic highest or lowest fluid level elevations in a well occur first. To test the assumption, we conduct two simulations with a modified multiphase flow numerical code TMVOC that incorporates relative permeability‐saturation‐capillary head relations employed in the model. In one simulation, the initial condition is for fluid levels in a well at the historic highest elevations. In the other simulation, the initial condition is for fluid levels in a well at the historic lowest elevations. We change the boundary conditions so both historical conditions occur followed by generating the current condition. Results from the numerical simulations are compared to model predictions and show the assumption in the analytical model is reasonable. The analytical model can be used to develop/refine conceptual site models and for assessing potential LNAPL recovery endpoints, especially on sites with fluctuating fluid levels in wells. 相似文献
19.
The seepage meter and the mini-piezometer were utilized in an attempt to evaluate ground-water reservoir interactions over a 12-month period at Echo Bay in Lake Mead. In conjunction with these techniques three standard piezometers, refraction seismic data, and water chemistry data were utilized to interpret seepage device results. During a four-month period, from December 1979 to March 1980, an 8 ft (2.5 m) rise in reservoir stage, the reservoir contributed water to Echo Wash bank storage at rates of up to 0.29 gpd/ft2 (12 lpd/m2). Ground-water discharge occurred for the remainder of the project, during a stage decline from April 1980 to May 1980, a rise in June 1980, and leveling off and slight decline for the remainder of the year, July 1980 to December 1980. The maximum seepage meter ground-water discharge rate of 3.0 gpd/ft2 (122 lpd/m2) was recorded in December 1980. Seepage meter water chemistry data for June were similar to Lake Mead water chemistry and were interpreted to be previously recharged Lake Mead water. September water chemistry data showed two possible components of ground-water discharge, a high SpC calcium sulfate Echo Wash ground water and a lower SpC Lake Mead recharged bank storage water. December ground-water chemistry data showed discharge to be a high SpC calcium sulfate water similar to Echo Wash ground-water quality which was apparently unaffected by Lake Mead inflow. Mini-piezometer data were collected at each seepage meter site. However, these data usually did not provide correlative results with seepage meter data probably because of suspended sediment in the piezometer water column and plugging of the perforated tip. Seepage meters were successfully utilized to characterize reservoir ground-water interaction in Echo Bay. 相似文献
20.
Abstract This paper describes a study of groundwater flow in a coastal Miliolite limestone aquifer in western India. An examination of field information suggested that the transmissivity of the aquifer varies significantly between high and low groundwater heads. Pumping tests indicate that this is due to the development of major fissures in the upper part of the aquifer. A regional groundwater model with varying transmissivities is used to represent the field behaviour. The model is also used to examine the effect of artificial recharge on the alleviation of saline intrusion problems in the coastal area. 相似文献