首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 390 毫秒
1.
A Potential-Based Inversion of Unconfined Steady-State Hydraulic Tomography   总被引：1，自引：0，他引：1  
The importance of estimating spatially variable aquifer parameters such as transmissivity is widely recognized for studies in resource evaluation and contaminant transport. A useful approach for mapping such parameters is inverse modeling of data from series of pumping tests, that is, via hydraulic tomography. This inversion of field hydraulic tomographic data requires development of numerical forward models that can accurately represent test conditions while maintaining computational efficiency. One issue this presents is specification of boundary and initial conditions, whose location, type, and value may be poorly constrained. To circumvent this issue when modeling unconfined steady-state pumping tests, we present a strategy that analyzes field data using a potential difference method and that uses dipole pumping tests as the aquifer stimulation. By using our potential difference approach, which is similar to modeling drawdown in confined settings, we remove the need for specifying poorly known boundary condition values and natural source/sink terms within the problem domain. Dipole pumping tests are complementary to this strategy in that they can be more realistically modeled than single-well tests due to their conservative nature, quick achievement of steady state, and the insensitivity of near-field response to far-field boundary conditions. After developing the mathematical theory, our approach is first validated through a synthetic example. We then apply our method to the inversion of data from a field campaign at the Boise Hydrogeophysical Research Site. Results from inversion of nine pumping tests show expected geologic features, and uncertainty bounds indicate that hydraulic conductivity is well constrained within the central site area.  相似文献   

2.
Kai‐Yuan Ke 《水文研究》2014,28(3):1409-1421
This research proposes a combination of SWAT and MODFLOW, MD‐SWAT‐MODFLOW, to address the multi‐aquifers condition in Choushui River alluvial fan, Taiwan. The natural recharge and unidentified pumping/recharge are separately estimated. The model identifies the monthly pumping/recharge rates in multi‐aquifers so that the daily streamflow can be simulated correctly. A multi‐aquifers condition means a subsurface formation composed of at least the unconfined aquifer, the confined aquifer, and an in‐between aquitard. In such a case, the variation of groundwater level is related to pumping/recharge activities in vertically adjacent aquifer and the river‐aquifer interaction. Both factors in turn affect the streamflow performance. Results show that MD‐SWAT‐MODFLOW performs better than SWAT alone in terms of simulated streamflow, especially during low flow period, when pumping/recharge rates are properly estimated. A sensitivity analysis of individual parameter suggests that the vertical leakance may be the most sensitive among all investigated MODFLOW parameters in terms of the estimated pumping/recharge among aquifers, and the Latin‐Hypercube‐One‐factor‐At‐a‐Time sensitivity analysis indicates that the hydraulic conductivity of channel is the most sensitive to the model performance. It also points out the necessity to simultaneously estimate pumping/recharge rates in multi‐aquifers. The estimated net pumping rate can be treated as a lower bound of the actual local pumping rate. As a whole, the model provides the spatio‐temporal groundwater use, which gives the authorities insights to manage groundwater resources. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

3.
A new solution of transient confined–unconfined flow driven by a pumping well is developed and compared to previous approximate solutions of Moench and Prickett [Moench AF, Prickett TA. Radial flow in an infinite aquifer undergoing conversion from artesian to water table conditions. Water Resour Res 1972;8:494–9] and Hu and Chen [Hu L, Chen C. Analytical methods for transient flow to a well in a confined–unconfined aquifer. Ground Water 2008;46(4):642–6]. The problem is rewritten in dimensionless form with the Boltzmann transform. The nonlinear equation for flow in the unconfined zone is solved with the Runge–Kutta method. Position of the conversion interface is determined with an iteration scheme. This study shows that the confined–unconfined flow depends on three dimensionless parameters that represent the confined–unconfined storativity ratio (aD), the ratio of the initial hydraulic head over the aquifer thickness (fi), and the dimensionless pumping rate (qD). The rate of expansion of the unconfined zone increases with qD, but decreases with aD and fi. Differences between the two previous approximate solutions and the new solution of this study are observable in the estimated position of the conversion interface and the drawdown–time curves. The new solution can be applied to estimate the time for confined–unconfined conversion to occur (critical conversion time), and the time when the pumping well becomes dry (critical drying time). The critical conversion time is found to be very sensitive to the initial hydraulic head. The critical drying time is often much larger than the critical conversion time and may never be observed during a finite pumping period.  相似文献   

4.
A physically based inverse method is developed using hybrid formulation and coordinate transform to simultaneously estimate hydraulic conductivity tensors, steady‐state flow field, and boundary conditions for a confined aquifer under ambient flow or pumping condition. Unlike existing indirect inversion techniques, the physically based method does not require forward simulations to assess model‐data misfits. It imposes continuity of hydraulic head and Darcy fluxes in the model domain while incorporating observations (hydraulic heads, Darcy fluxes, or well rates) at measurement locations. Given sufficient measurements, it yields a well‐posed inverse system of equations that can be solved efficiently with coarse grids and nonlinear optimization. When pumping and injection are active, well rates are used as measurements and flux sampling is not needed. The method is successfully tested on synthetic aquifer problems with regular and irregular geometries, different hydrofacies and flow patterns, and increasing conductivity anisotropy ratios. All problems yield stable inverse solutions under increasing head measurement errors. For a given set of observations, inversion accuracy is strongly affected by the conductivity anisotropy ratio. Conductivity estimation is also affected by flow pattern: within a hydrofacies, when Darcy flux component is very small, the corresponding directional conductivity perpendicular to streamlines becomes less identifiable. Finally, inversion is successful even if the location of aquifer boundaries is unknown. In this case, the inversion domain is defined by the location of the measurements.  相似文献   

5.
Seasonal signals of stable isotopes in precipitation, combined with measurements of isotope ratios in soil water, can be used for quantitative estimation of groundwater recharge rates. This study investigates the applicability of using the piston flow principle and the peak shift displacement method to estimate actual groundwater recharge rates in a humid Nordic region located in the province of Quebec, Canada. Two different sites with and without vegetation (C1 and C2) in an unconfined aquifer were tested by measuring soil water isotope ratios (18O/16O and 2H/1H) and volumetric pore water content. Core samples were obtained along the vadose zone down to the groundwater table at the two sites (2.45 m for Site C1 and 4.15 m for Site C2). The peak shift method to estimate groundwater recharge rates was shown to be accurate only in certain specific conditions inherent to the soil properties and the topographical situation of the investigated sites. Indeed, at Site C2, recharge from the snowmelt could not be estimated because of heterogeneity in the lower part of the vadose zone. At this same site the later recharge after the snowmelt (in the period from late spring to early autumn) could be estimated accurately because the upper part of the vadose zone was homogeneous. Furthermore, at site C1, runoff/runon phenomena hampered calculations of actual infiltration and thus produced inaccurate results for recharge. These two different site effects (heterogeneity in the first site and runoff/runon in the other site) were identified as being limiting factors in the accurate assessment of actual recharge. This study therefore recommends the use of the peak shift method for (1) humid Nordic regions, (2) homogeneous and thick vadose zones, and (3) areas with few or limited site effects (runoff/runon).  相似文献   

6.
Following a brief overview of the history and the development of the Surface Wave Method—with a focus on techniques for processing and inverting field data—a Simplified Inversion Method (SIM) is described, which constitutes an improvement of the Satoh et al. (1991) [1] method. The SIM is a direct inversion method of surface wave dispersion data, making use of a penetration depth coefficient, aR, whose value is a function of Poisson's ratio and the overall shape of the dispersion curve. In the present study the coefficient aR has been evaluated using data from (a) an extensive database compiled from the technical literature and containing results of inverted surface wave measurements and nearby cross-hole/down-hole measurements, (b) results of side by side surface wave and cross-hole measurements, performed at five sites in the course of this study, (c) finite element analyses simulating the performance of surface wave measurements and thus providing “virtual” data, and (d) applying a current advanced inversion code, available on the Web. Based on all the above data, optimum values of aR (and of the corresponding uncertainty of the derived Vso vs. depth profile) have been estimated. These values were found to be independent of depth from ground surface. The results of all analyses and comparisons indicate that for the majority of realistic soil profiles (including cases of normal and inverse dispersion conditions) the proposed SIM provides very reliable Vso vs. depth profiles when a value of aR=0.63–0.67 is used in the inversion process. It is concluded that the SIM can be used with confidence as a direct inversion method of surface wave dispersion data.  相似文献   

7.
In this study, we attempted to analyse a drawdown pattern around a pumping well in an unconfined sandy gravelly aquifer constructed in a laboratory tank by means of both experimental and numerical modelling of groundwater flow. The physical model consisted of recharge, aquifer and discharge zones. Permeability and specific yield of the aquifer material were determined by Dupuit approximation under steady‐state flow and stepwise gravitational drainage of groundwater, respectively. The drawdown of water table in pumping and neighbouring observation wells was monitored to investigate the effect of no‐flow boundary on the drawdown pattern during pumping for three different boundary conditions: (i) no recharge and no discharge with four no‐flow boundaries (Case 1); (ii) no recharge and reservoir with three no‐flow boundaries (Case 2); (iii) recharge and discharge with two no‐flow boundaries (Case 3). Based on the aquifer parameters, numerical modelling was also performed to compare the simulated drawdown with that observed. Results showed that a large difference existed between the simulated drawdown and that observed in wells for all cases. The reason for the difference could be explained by the formation of a curvilinear type water table between wells rather than a linear one due to a delayed response of water table in the capillary fringe. This phenomenon was also investigated from a mass balance study on the pumping volume. The curvilinear type of water table was further evidenced by measurement of water contents at several positions in the aquifer between wells using time domain reflectometry (TDR). This indicates that the existing groundwater flow model applicable to an unconfined aquifer lacks the capacity to describe a slow response of water table in the aquifer and care should be taken in the interpretation of water table formation in the aquifer during pumping. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

8.
Following a brief overview of the history and the development of the Surface Wave Method—with a focus on techniques for processing and inverting field data—a Simplified Inversion Method (SIM) is described, which constitutes an improvement of the Satoh et al. (1991) [1] method. The SIM is a direct inversion method of surface wave dispersion data, making use of a penetration depth coefficient, aR, whose value is a function of Poisson's ratio and the overall shape of the dispersion curve. In the present study the coefficient aR has been evaluated using data from (a) an extensive database compiled from the technical literature and containing results of inverted surface wave measurements and nearby cross-hole/down-hole measurements, (b) results of side by side surface wave and cross-hole measurements, performed at five sites in the course of this study, (c) finite element analyses simulating the performance of surface wave measurements and thus providing “virtual” data, and (d) applying a current advanced inversion code, available on the Web. Based on all the above data, optimum values of aR (and of the corresponding uncertainty of the derived Vso vs. depth profile) have been estimated. These values were found to be independent of depth from ground surface. The results of all analyses and comparisons indicate that for the majority of realistic soil profiles (including cases of normal and inverse dispersion conditions) the proposed SIM provides very reliable Vso vs. depth profiles when a value of aR=0.63–0.67 is used in the inversion process. It is concluded that the SIM can be used with confidence as a direct inversion method of surface wave dispersion data.  相似文献   

9.
The success of modeling groundwater is strongly influenced by the accuracy of the model parameters that are used to characterize the subsurface system. However, the presence of uncertainty and possibly bias in groundwater model source/sink terms may lead to biased estimates of model parameters and model predictions when the standard regression‐based inverse modeling techniques are used. This study first quantifies the levels of bias in groundwater model parameters and predictions due to the presence of errors in irrigation data. Then, a new inverse modeling technique called input uncertainty weighted least‐squares (IUWLS) is presented for unbiased estimation of the parameters when pumping and other source/sink data are uncertain. The approach uses the concept of generalized least‐squares method with the weight of the objective function depending on the level of pumping uncertainty and iteratively adjusted during the parameter optimization process. We have conducted both analytical and numerical experiments, using irrigation pumping data from the Republican River Basin in Nebraska, to evaluate the performance of ordinary least‐squares (OLS) and IUWLS calibration methods under different levels of uncertainty of irrigation data and calibration conditions. The result from the OLS method shows the presence of statistically significant (p < 0.05) bias in estimated parameters and model predictions that persist despite calibrating the models to different calibration data and sample sizes. However, by directly accounting for the irrigation pumping uncertainties during the calibration procedures, the proposed IUWLS is able to minimize the bias effectively without adding significant computational burden to the calibration processes.  相似文献   

10.
The identification of groundwater parameters in heterogeneous systems is a major challenge in groundwater modeling. Flexible parameterization methods are needed to assess the complexity of the spatial distributions of these parameters in real aquifers. In this article, we introduce an adaptative parameterization to identify the distribution of hydraulic conductivity within the large‐scale (4400 km2) Upper Rhine aquifer. The method is based on adaptative multiscale triangulation (AMT) coupled with an inverse problem procedure that identifies the parameters' distributions by reducing the error between measured and simulated heads. The AMT method has the advantage of combining both zonation and interpolation approaches. The AMT method uses area‐based interpolation rather than an interpolation based on stochastic features. The method is applied to a standard 2D groundwater model that takes into account the interactions between the aquifer and surface water bodies, groundwater recharge, and pumping wells. The simulation period covers 204 months, from January 1986 to December 2002. Recordings at 109 piezometers are used for model calibration. The simulated heads are globally quite accurate and reproduce the main dynamics of the system. The local hydraulic conductivities resulting from the AMT method agree qualitatively with existing local experimental observations across the Rhine aquifer.  相似文献   

11.
This study investigates the inverse solution on a buried and polarized sphere-shaped body using the self-potential method via multilayer perceptron neural networks (MLPNN). The polarization angle (α), depth to the centre of sphere (h), electrical dipole moment (K) and the zero distance from the origin (x 0) were estimated. For testing the success of the MLPNN for sphere model, parameters were also estimated by the traditional Damped Least Squares (Levenberg–Marquardt) inversion technique (DLS). The MLPNN was first tested on a synthetic example. The performance of method was also tested for two S/N ratios (5 % and 10 %) by adding noise to the same synthetic data, the estimated model parameters with MLPNN and DLS method are satisfactory. The MLPNN also applied for the field data example in ?zmir, Urla district, Turkey, with two cross-section data evaluated by MLPNN and DLS, and the two methods showed good agreement.  相似文献   

12.
ABSTRACT

This study presents a systematic illustration quantifying how misleading the calibration results of a groundwater simulation model can be when recharge rates are considered as the model parameters to be estimated by inverse modelling. Three approaches to recharge estimation are compared: autocalibration (Model 1), the empirical return coefficient method (Model 2), and distributed hydrological modelling using the Soil and Water Assessment Tool, SWAT (Model 3). The methodology was applied in the Dehloran Plain, western Iran, using the MODFLOW modular flow simulator and the PEST method for autocalibration. The results indicate that, although Model 1 performed the best in simulating water levels at observation wells in the calibration stage, it did not perform satisfactorily in real future scenarios. Model 3, with SWAT-based recharge rates, performed better than the other models in the validation stage. By not evaluating the model performance solely on calibration results, we demonstrate the relative significance of using more accurate recharge estimates when calibrating groundwater simulation models.
EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR M. Besbes  相似文献   

13.
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.  相似文献   

14.
The spatial moments of a contaminant plume undergoing bio-attenuation are coupled to the moments of microbial populations effecting that attenuation. In this paper, a scalable inverse method is developed for estimating field-scale Monod parameters such as the maximum microbial growth rate (μmax), the contaminant half saturation coefficient (Ks), and the contaminant yield coefficient (Ys). The method uses spatial moments that characterize the distribution of dissolved contaminant and active microbial biomass in the aquifer. A finite element model is used to generate hypothetical field-scale data to test the method under both homogeneous and heterogeneous aquifer conditions. Two general cases are examined. In the first, Monod parameters are estimated where it is assumed a microbial population comprised of a single bacterial species is attenuating one contaminant (e.g., an electron donor and an electron acceptor). In a second case, contaminant attenuation is attributed to a microbial consortium comprised of two microbial species, and Monod parameters for both species are estimated. Results indicate the inverse method is only slightly sensitive to aquifer heterogeneity and that estimation errors decrease as the sampling time interval decreases with respect to the groundwater travel time between sample locations. Optimum conditions for applying the scalable inverse method in both space and time are investigated under both homogeneous and heterogeneous aquifer conditions.  相似文献   

15.
We examined the fire‐induced changes in groundwater recharge rate. This aspect is particularly important in the case of large forested areas growing over a coastal aquifer affected by saltwater intrusion. In the Ravenna coastal area (Italy), pine forests grow on coastal dune belts, overlying a sandy unconfined aquifer, which is strongly affected by marine ingression. Three groundwater profiles across the forest and perpendicular to the coastline were monitored for groundwater level, physical, and chemical parameters. The aims were to define groundwater quality, recharge rate, freshwater volume, and highlight change, which occurred after a forest fire with reference to pre‐fire conditions. Analytical solutions based on Darcy Law and the Dupuit Equation were applied to calculate unconfined flow and compare recharge rates among the profiles. The estimated recharge rates increased in the partially and completely burnt areas (219 and 511 mm year?1, respectively) compared with the pristine pine forest area (73 mm year?1). Although pre‐fire conditions were similar in all monitored profiles, a post‐fire decrease in salinity was observed across the burnt forest, along with an increase in infiltration and freshwater lens thickness. This was attributed to decrease canopy interception and evapotranspiration caused by vegetation absence after the fire. This research provided an example of positive forest fire feedback on the quantity and quality of fresh groundwater resources in a lowland coastal aquifer affected by saltwater intrusion, with limited availability of freshwater resources. The fire provided an opportunity to evaluate a new forest management approach and consider the restoration and promotion of native dune herbaceous vegetation.  相似文献   

16.
We jointly invert field data of flowmeter and multiple pumping tests in fully screened wells to estimate hydraulic conductivity using a geostatistical method. We use the steady-state drawdowns of pumping tests and the discharge profiles of flowmeter tests as our data in the inference. The discharge profiles need not be converted to absolute hydraulic conductivities. Consequently, we do not need measurements of depth-averaged hydraulic conductivity at well locations. The flowmeter profiles contain information about relative vertical distributions of hydraulic conductivity, while drawdown measurements of pumping tests provide information about horizontal fluctuation of the depth-averaged hydraulic conductivity. We apply the method to data obtained at the Krauthausen test site of the Forschungszentrum Jülich, Germany. The resulting estimate of our joint three-dimensional (3D) geostatistical inversion shows an improved 3D structure in comparison to the inversion of pumping test data only.  相似文献   

17.
The article presents semi‐analytical mathematical models to asses (1) enhancements of seepage from a canal and (2) induced flow from a partially penetrating river in an unconfined aquifer consequent to groundwater withdrawal in a well field in the vicinity of the river and canal. The nonlinear exponential relation between seepage from a canal reach and hydraulic head in the aquifer beneath the canal reach is used for quantifying seepage from the canal reach. Hantush's (1967) basic solution for water table rise due to recharge from a rectangular spreading basin in absence of pumping well is used for generating unit pulse response function coefficients for water table rise in the aquifer. Duhamel's convolution theory and method of superposition are applied to obtain water table position due to pumping and recharge from different canal reaches. Hunt's (1999) basic solution for river depletion due to constant pumping from a well in the vicinity of a partially penetrating river is used to generate unit pulse response function coefficients. Applying convolution technique and superposition, treating the recharge from canal reaches as recharge through conceptual injection wells, river depletion consequent to variable pumping and recharge is quantified. The integrated model is applied to a case study in Haridwar (India). The well field consists of 22 pumping wells located in the vicinity of a perennial river and a canal network. The river bank filtrate portion consequent to pumping is quantified.  相似文献   

18.
Aquifer hydraulic parameters are commonly inferred from constant-rate pumping tests, while variable pumping rates are frequently encountered in actual field conditions. In this study, we propose a generally applicable dimensionless form of the analytical solution for variable-rate pumping tests in confined aquifers. In particular, we adopt a piecewise-linear fitting of variable pumping rates and propose a new type-curve method for estimating the hydraulic conductivity (K ) and specific storage (Ss ) of the investigated confined aquifer. For each test, a series of type curves, which depend on the variable pumping rates, the location of observation wells and the introduced first dimensionless inflection time, need to be provided for matching the observed drawdown data on a log-log graph. We first demonstrate the applicability and robustness of this method through a synthetic pumping test. Subsequently, we apply this method to analyze drawdown data from four pumping tests conducted within a multilayered aquifer/aquitard system in Wuxi city, Jiangsu Province, China. The parameter estimates are then compared with those reported by PEST. The K and Ss values estimated by the new type-curve method are found to be quite close to PEST-based estimates. Parameter estimation results demonstrate the difference in K and Ss values between observation wells. The difference could be attributed to the spatial heterogeneity in K and Ss . A future research topic may focus on the characterization of K and Ss heterogeneity with the currently available drawdown data from variable-rate pumping tests.  相似文献   

19.
This paper presents the Kalman Filtered Double Constraint Method (DCM‐KF) as a technique to estimate the hydraulic conductivities in the grid blocks of a groundwater flow model. The DCM is based on two forward runs with the same initial grid block conductivities, but with alternating flux‐head conditions specified on parts of the boundary and the wells. These two runs are defined as: (1) the flux run, with specified fluxes (recharge and well abstractions), and (2) the head run, with specified heads (measured in piezometers). Conductivities are then estimated as the initial conductivities multiplied by the fluxes obtained from the flux run and divided by the fluxes obtained from the head run. The DCM is easy to implement in combination with existing models (e.g., MODFLOW). Sufficiently accurate conductivities are obtained after a few iterations. Because of errors in the specified head‐flux couples, repeated estimation under varying hydrological conditions results in different conductivities. A time‐independent estimate of the conductivities and their inaccuracy can be obtained by a simple linear KF with modest computational requirements. For the Kleine Nete catchment, Belgium, the DCM‐KF yields sufficiently accurate calibrated conductivities. The method also results in distinguishing regions where the head‐flux observations influence the calibration from areas where it is not able to influence the hydraulic conductivity.  相似文献   

20.
We use particle tracking to determine contributing areas (CAs) to wells for transient flow models that simulate cyclic domestic pumping and extreme recharge events in a small synthetic watershed underlain by dipping sedimentary rocks. The CAs consist of strike-oriented bands at locations where the water table intersects high-hydraulic conductivity beds, and from which groundwater flows to the pumping well. Factors that affect the size and location of the CAs include topographic flow directions, rock dip direction, cross-bed fracture density, and position of the well relative to streams. For an effective fracture porosity (ne) of 10−4, the fastest advective travel times from CAs to wells are only a few hours. These results indicate that wells in this type of geologic setting can be highly vulnerable to contaminants or pathogens flushed into the subsurface during extreme recharge events. Increasing ne to 10−3 results in modestly smaller CAs and delayed well vulnerability due to slower travel times. CAs determined for steady-state models of the same setting, but with long-term average recharge and pumping rates, are smaller than CAs in the models with extreme recharge. Also, the earliest-arriving particles arrive at the wells later in the steady-state models than in the extreme-recharge models. The results highlight the importance of characterizing geologic structure, simulating plausible effective porosities, and simulating pumping and recharge transience when determining CAs in fractured rock aquifers to assess well vulnerability under extreme precipitation events.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号