首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Yang  Heejun  Tawara  Yasuhiro  Shimada  Jun  Kagabu  Makoto  Okumura  Azusa 《Hydrogeology Journal》2021,29(6):2091-2105

The hydraulic conductivity of an unconfined carbonate aquifer at the uplifted atoll of Minami-Daito, Japan, was evaluated by a combination of cross-spectral analysis, analytical solution, and density-dependent groundwater modeling based on observed groundwater levels in 15 wells and at sea level. The island area was divided into 10 subregions based on island morphology and on inland propagation of ocean tides. The hydraulic conductivity was obtained for each subregion using analytical solutions based on phase lags of M2 constituents of ocean tides at each well by assuming two aquifer thicknesses (300 and 1,800 m) and two effective porosities (0.1 and 0.3). The density-dependent groundwater model evaluated the hydraulic conductivity of the subregions by reproducing observed groundwater levels. The hydraulic conductivity in the subregions was estimated as 3.46?×?10?3 to 6.35?×?10?2 m/s for aquifer thickness of 300 m and effective porosity of 0.1, and as 1.73?×?10?3 to 3.17?×?10?2 m/s for aquifer thickness of 1,800 m and the effective porosity of 0.3. It was higher in southern and northern areas, and higher in interior lowland than in the western and eastern areas. Fissures and dolomite distributions on the island control differences of the omnidirectional ocean tidal propagation and cause these differences in hydraulic conductivity. The method used for this study may also be applicable to other small islands that have few or no data for hydraulic conductivity.

  相似文献   

2.
Sandstone is a common lithology in a number of groundwater reservoirs. Studying the skeleton sandstone architectural units, therefore, lays the basis for characterizing aquifer systems, groundwater quality, and resource evaluation. This comprehensive analysis of Cretaceous aquiferous sandstones in the Ordos basin, China, shows that there exists a basin-scale skeleton sandstone in the Luohe Formation which contains 11 isolated barrier beds, 12 small skeleton sandstone bodies in Huanhe Formation, and 3 in the Luohandong Formation. The spatial structure and superimposed relationship as well as the medium properties of these skeleton sandstones and isolated barrier beds can be shown by 3D visualization models. Simultaneously, resource quantity can be evaluated with the 3D inquiry functions. The comparison between property models and structural models indicates that the salinity of groundwater of the Luohe Formation has a close connection with the locations of isolated barrier beds that contain abundant gypsum. Through quantitative calculation, groundwater resource of the Cretaceous Luohe and Luohandong formations is estimated to be 1.6×10~(12) m~3, and the total groundwater resource of the Cretaceous system in the Ordos basin is more than 2×10~(12) m~3.  相似文献   

3.
Combined geophysical techniques such as multi-electrode resistivity, induced polarization, and borehole geophysical techniques were carried out on volcano-sedimentary rocks in the north of Gemas as part of the groundwater resource’s investigations. The result identifies four resistivity units: the tuffaceous mudstone, tuffaceous sandstone, the tuff bed, and the shale layer. Two types of aquifer systems in terms of storage were identified within the area: one within a fracture system (tuff), which is the leaky area through which vertical flow of groundwater occurs, and an intergranular property of the sandy material of the aquifer which includes sandstone and tuffaceous sandstone. The result also reveals that the aquifer occupies a surface area of about 3,250,555 m2 with a mean depth of 43.71 m and a net volume of 9.798?×?107?m3. From the approximate volume of the porous zone (28 %) and the total aquifer volume, a usable capacity of (274.339?±?30.177)?×?107?m3 of water in the study area can be deduced. This study provides useful information that can be used to develop a much broader understanding of the nature of groundwater potential in the area and their relationship with the local geology.  相似文献   

4.
An innovative mode of groundwater recharge to a buried esker aquifer is considered. The current conceptual model affords a natural safeguard to underlying aquifers from the overlying muds. A hypothesis of groundwater recharge to a buried esker aquifer via preferential pathways across its overlying muds is tested here by heuristic numerical one-dimensional and two-dimensional modeling simulations. The hypothesis has been tested against two other conventionally accepted scenarios involving: (1) distal esker outcrop areas and (2) remote shallow-bedrock recharge areas. The main evidence comes from documented recharge pressure pulses in the overlying mud aquitard and in the underlying esker hydraulic-head time series for the Vars-Winchester esker aquifer in Eastern Ontario, Canada. These perturbations to the potentiometric surface are believed to be the aquifer response to recharge events. The migration rate of these pressure pulses is directly related to the hydraulic diffusivity of the formation. The measured response time and response amplitude between singular radar precipitation events and well hydrographs constituted the heuristic model calibration targets. The main evidence also includes mud-layering deformation (water escape features) which was observed in seismic surveys of the over-esker muds. These disturbed stratigraphic elements provide a realistic mechanism for migrating water to transit through the muds. The effective hydraulic conductivities of these preferential pathways in the muds were estimated to be between 2?×?10?6 and 7?×?10?6 m/s. The implications of these findings relate to the alleged natural safeguard of these overlying muds.  相似文献   

5.
A fracture-karst aquifer is a karst aquifer with a fractured rock matrix, and its parameters are difficult to determine. Two sequential pumping tests in a fracture-karst aquifer system at the Zhangji well field of China are considered, one carried out before (in 2000) and one after (in 2005) the operation of a pumping station in the well field (2003–2005). The sequential tests serve multiple purposes. First, they provide a cross check of the parameters obtained. Second, they can be used to assess the effect of long-term groundwater exploitation of the aquifer. A three-dimensional finite-element transient flow model has been developed to simulate groundwater flow at the site. Generally good agreement has been found between the simulated and observed hydraulic heads for both tests. The hydraulic parameters obtained from the 2005 test are generally consistent with their counterparts from the 2000 test. However, a small but steady increase of hydraulic conductivities from 2000 to 2005 at the site has been observed. A 10-year prediction of groundwater resources has been made and indicates that the well field can accommodate the proposed 8.0?×?104 m3/day exploitation rate under relative drought conditions without causing a steady decline of groundwater levels.  相似文献   

6.
The present study indicates that the factors controlling the hydraulic relation between surface water and groundwater at the western lake shoreline change from one locality to another. This depends upon the lithological characteristics and the major structures. In the southern sectors, sedimentation at the bottom and sides of the lake prevents the water movement to the Nubian sandstone aquifer. The potentiometric map reveals that the water level altitudes range between 170 m in the vicinity of the lakeshore line and 110 m west of the lake. The groundwater flow lines show that the main recharge to the aquifer comes from the southwest direction, as well as from the lake inland to variable distances (about 30 Km). During the present study, Darcy’s law was applied to calculate the recharge from the western shoreline of Lake Nasser to the adjacent Nubian aquifer. The maximum value of seepage was at Garf Hussein (27.71?×?106 m³/year), which may be related to high permeability and hydraulic gradient. Also, it may be related to the N–S strike faults that cut the area on both sides of the Lake, and the groundwater is expected to have free circulation through the faults of this trend. The minimum value was recorded in Adindan section (0.61?×?106 m³/year). This may be related to the limited recharge from the lake to the aquifer, due to the sedimentation that dislocates this recharge.  相似文献   

7.
The current study aimed to evaluate hydrogeologically the Nubian sandstone aquifer in El-Bahariya Oasis. It represents the main water-bearing horizon in the study area and consists of continental elastic sediments, mainly sandstone alternating with shale and clays. The general flow lines are directed from SW to NE direction, as detected from the constructed potentiometric head contour map. The piezometric surface reaches 149 m in El-Heiz area at the southern part, while it reaches 90 m at the northern, reflecting higher pressure head of the aquifer in the southern part. The map also illustrates that the southern part is considered as the most promising location for development. The structural elements play an important role in the deposition and distribution of the sedimentary succession of the Nubian sandstone sediments. Consequently, this sedimentary pattern affects the occurrences and movements of the groundwater within the aquifer system. Along the structurally high areas, in the study area, the piezometric head increases, while the reverse is recorded along the structurally low areas. The step-drawdown tests data were carried out by calculating the aquifer loss coefficient (B) and the well loss constant (C). The B values are smaller compared with C values, indicating that the aquifer under pressure has a behavior of leaky aquifer; therefore, it shows hydraulic connection with surrounding formation. The values of well efficiency range from 78.50% to 87.76%. Analysis of 12 pumping test data (constant discharge tests) was carried out in order to calculate the Nubian aquifer hydraulic parameters (transmissivity, hydraulic conductivity, and storage coefficient). The transmissivity values decrease from 3,045 m2/day in the southern part (El-Heiz area) to 236 m2/day in the northeastern part (El-Harra area). Accordingly, the aquifer classified as a high to moderate potentiality. Transmissivity contour map observes gradual increase of transmissivity values from the southern to northeastern direction. This may be due to the increase of shale or clay content in the concerned aquifer in that direction. The storage coefficient values range between 1.04 × 10?4 and 5.22 × 10?3, as obtained from the results of pumping test analysis, which ensure that the Nubian sandstone aquifer is classified as semi-confined to confined aquifer type. The S values show a decrease from southwest to northeast direction as detected from S-map. The hydraulic conductivity values vary from to 0.46 m/day in the northern part to 10.88 m/day in the southern part with an average of 5.67 m/day. According to the classification based on K values, the aquifer is mainly composed of coarse sand.  相似文献   

8.
Phu tok Aquifer is the most productive water bearing rock in the Northeastern region of Thailand. Generally, well drilled in this aquifer yields approximately 20–50 m3/hr with some wells yield over 150 m3/hr. This aquifer characterized both unconfined and confined aquifer. Data collected from core sampling and TV bore hole show water bearing fractures in the fine grained massive sandstone intercalated with thin shale layers. These intermittently bedding plane fracture zones were found vertically at 60-65 and 95-110 m below ground surface. Three directions of vertically fractures were also found at 103-104 m in the test well with underlain densely sandstone until 120 m of depth. Wells drilled to 60-70 m fractures at the area where surface elevation between 150 and 170 m above mean sea level were artesian well with rising water about 1-3 m. Natural discharge rate from 2-5 inches of diameter casing is 5-10 m3/hr. Fractures at 90-110 m yield 150 m3/hr through 6 inch casing well with 5-6 m of risen water above the ground. It is found from pumping test of Phu Tok Aquifer that hydraulic conductivity of unconfined aquifer at shallower than 50 m is 0.005-17 m/d. Transmissivity and storage coefficient are 0.05-20 m2/d and 7×10-3-0.725 respectively. The confined aquifer at depth not exceed 90 m has hydraulic conductivity value of 0.2-10 m/d while transmissivity and storage coefficient are 3.19-150 m2/d and 1×10-10-1.6×10-2. Another confined aquifer at 90-120 m of depth have hydraulic conductivity value of 0.08-15 m/d and transmissivity and storage coefficient values of 1.7-178 m2/d and 4×10-7-4.5×10-3 respectively.  相似文献   

9.
The Khor Arbaat basin is the main source of potable water supply for the more than 750,000 inhabitants of Port Sudan, eastern Sudan. The variation in hydraulic conductivity and storage capacity is due to the heterogeneity of the sediments, which range from clay and silt to gravely sand and boulders. The water table rises during the summer and winter rainy seasons; it reaches its lowest level in the dry season. The storage capacity of the Khor Arbaat aquifer is estimated to be 21.75?×?106 m3. The annual recharge through the infiltration of flood water is about 1.93?×?106 m3. The groundwater recharge, calculated as underground inflow at the ‘upper gate’, is 1.33?×?105 m3/year. The total annual groundwater recharge is 2.06?×?106 m3. The annual discharge through underground outflow at the ‘lower gate’ (through which groundwater flows onto the coastal plain) is 3.29?×?105 m3/year. Groundwater discharge due to pumping from Khor Arbaat basin is 4.38?×?106 m3/year on average. The total annual groundwater discharge is about 4.7?×?106 m3. A deficit of 2.6?×?106 m3/year is calculated. Although the total annual discharge is twice the estimated annual recharge, additional groundwater flow from the fractured basement probably balances the annual groundwater budget since no decline is observed in the piezometric levels.  相似文献   

10.
The Guarani aquifer system (GAS) represents one of the biggest aquifers in the world and is the most relevant groundwater resource in South America. For the first time, by combining field and laboratory measurements, a high-resolution aquifer analog model of fluvial–aeolian sediments of the GAS in São Paulo State (Brazil) is constructed. Three parallel sections of frontal outcrops, 28 m × 5.8 m, and two parallel sections of lateral outcrops, 7 m × 5.8 m, are recorded during open-pit mining of sandy sediments and describe in detail the three-dimensional distribution of the local lithofacies and hydrofacies. Variations of hydraulic conductivity, K, and porosity, n, are resolved on the centimeter scale, and the most permeable units of the fluvial–aeolian facies association are identified. The constructed aquifer analog model shows moderate hydraulic heterogeneity and a mean K value of 1.36 × 10?4 m/s, which is greater than the reported range of K values for the entire GAS in São Paulo State. The results suggest that the examined sedimentary unit constitutes a relevant portion of the GAS in São Paulo State in the context of groundwater extraction and pollution. Moreover, the constructed aquifer analog is considered an ideal basis for future numerical model experiments, aiming at in-depth understanding of the groundwater flow and contaminant transport patterns at this GAS portion or at comparable fluvial–aeolian facies associations.  相似文献   

11.
Physical and hydraulic properties of sediment from two karst aquifers were measured to determine (1) the similarity of sediment between karst aquifer systems and (2) the importance of sediment in modeling flow through karst aquifers. The sediment from the two systems was similar in size and composition. Within both aquifers, the silt-sized sediment was composed primarily of quartz, with minor amounts of plagioclase and clays. Hydraulic conductivity of the sediment measured directly (falling-head test) ranged from 1.61×10−7 to 1.33×10−6 m s–1 and estimated using the Campbell equation ranged from 8.30×10−8 to 8.98×10−7 m s–1. These values of hydraulic conductivity fall within the span of hydraulic conductivities for carbonate rocks, indicating that the sediment and carbonate matrix could be represented as one mathematical unit in modeling flow through karst aquifers. Statistical agreement in the hydraulic conductivity values generated by the two methods indicates that the estimation technique could be used to calculate hydraulic conductivities; thus allowing karst scientist to collect bulk sediment samples instead of having to collect cores from within karst aquifers. Electronic Publication  相似文献   

12.
The River Gash Basin is filled by the Quaternary alluvial deposits, unconformably overlying the basement rocks. The alluvial deposits are composed mainly of unconsolidated layers of gravel, sand, silt, and clays. The aquifer is unconfined and is laterally bounded by the impermeable Neogene clays. The methods used in this study include the carry out of pumping tests and the analysis of well inventory data in addition to the river discharge rates and other meteorological data. The average annual discharge of the River Gash is estimated to be 1,056?×?106 m3 at El Gera gage station (upstream) and 587?×?106 m3 at Salam-Alikum gage station (downstream). The annual loss mounts up to 40% of the total discharge. The water loss is attributed to infiltration and evapotranspiration. The present study proofs that the hydraulic conductivity ranges from 36 to 105 m/day, whereas the transmissivity ranges from 328 to 1,677 m2/day. The monitoring of groundwater level measurements indicates that the water table rises during the rainy season by 9 m in the upstream and 6 m in the midstream areas. The storage capacity of the upper and middle parts of the River Gash Basin is calculated as 502?×?106 m3. The groundwater input reach 386.11?×?106 m3/year, while the groundwater output is calculated as 365.98?×?106 m3/year. The estimated difference between the input and output water quantities in the upper and middle parts of the River Gash Basin demonstrates a positive groundwater budget by about 20?×?106 m3/year  相似文献   

13.
The hydrogeologic and hydraulic characteristics of a lateritic terrain in West Bengal, India, were investigated. Test drilling was conducted at ten sites and grain-size distribution curves (GSDCs) were prepared for 275 geologic samples. Performance evaluation of eight grain-size-analysis (GSA) methods was carried out to estimate the hydraulic conductivity (K) of subsurface formations. Finally, the GSA results were validated against pumping-test data. The GSDCs indicated that shallow aquifer layers are coarser than the deeper aquifer layers (uniformity coefficient 0.19–11.4). Stratigraphy analysis revealed that both shallow and deep aquifers of varying thickness exist at depths 9–40 and 40–79 m, respectively. The mean K estimates by the GSA methods are 3.62–292.86 m/day for shallow aquifer layers and 0.97–209.93 m/day for the deeper aquifer layers, suggesting significant aquifer heterogeneity. Pumping-test data indicated that the deeper aquifers are leaky confined with transmissivity 122.69–693.79 m2/day, storage coefficient 1.01?×?10?7–2.13?×?10?4 and leakance 2.01?×?10?7–34.56?×?10?2 day?1. Although the K values yielded by the GSA methods are generally larger than those obtained from the pumping tests, the Slichter, Harleman and US Bureau Reclamation (USBR) GSA methods yielded reasonable values at most of the sites (1–3 times higher than K estimates by the pumping-test method). In conclusion, more reliable aquifers exist at deeper depths that can be tapped for dependable water supply. GSA methods such as Slichter, Harleman and USBR can be used for the preliminary assessment of K in lateritic terrains in the absence of reliable field methods.  相似文献   

14.
A conceptual groundwater flow model was developed for the crystalline aquifers in southeastern part of the Eastern region, Ghana. The objective was to determine approximate levels of groundwater recharge, estimate aquifer hydraulic parameters, and then test various scenarios of groundwater extraction under the current conditions of recharge. A steady state groundwater flow model has been calibrated against measured water levels of 19 wells in the area. The resulting recharge is estimated to range from 8.97 × 10?5 m/d to 7.14 × 10?4 m/d resulting in a basin wide average recharge of about 9.6% of total annual precipitation, which results in a basin wide quantitative recharge of about 2.4 million m3/d in the area. This compares to recharge estimated from the chloride mass balance of 7.6% of precipitation determined in this study. The general groundwater flow in the area has also been determined to conform to the general northeast–southwest structural grain of the country. The implication is that the general hydrogeology is controlled by post genetic structural entities imposed on the rocks to create ingresses for sufficient groundwater storage and transport. Calibrated aquifer hydraulic conductivities range between 0.99 m/d and over 19.4 m/d. There is a significant contribution of groundwater discharge to stream flow in the study area. Increasing groundwater extraction will have an effect on stream flow. This study finds that the current groundwater extraction levels represent only 0.17% of the annual recharge from precipitation, and that groundwater can sustain future increased groundwater demands from population growth and industrialization.  相似文献   

15.
A multi-layered aquifer, typical of riverbank alluvial deposits in Korea, was studied to determine the hydrologic properties. The geologic logging showed that the subsurface of the study site was comprised of four distinctive hydrogeologic units: silt, sand, highly weathered and fresh bedrock layers. The electrical resistivity survey supplied information on lateral extension of hydrogeologic strata only partially identified by a limited number of the geologic loggings. The laboratory column tracer test for the recovered core of the sand layer resulted in a hydraulic conductivity of 5.00×10−2 cm/s. The slug tests performed in the weathered rock layer yielded hydraulic conductivities of 4.32–7.72×10−4 cm/s. Hydraulic conductivities for the sand layer calculated from the breakthrough curves of bromide ranged between 2.08×10−3 and 2.44×10−2 cm/s with a geometric mean of 6.89×10−3 cm/s, which is 7 times smaller than that from the laboratory column experiment. The trend of increasing hydraulic conductivity with an increase in tracer travel length is likely a result of the increased likelihood of encountering a high conductivity zone as more of the aquifer is tested. The combined hydrogeologic site characterization using hydraulic tests, tracer tests, and column test with geologic loggings and geophysical survey greatly enhanced the understanding of the hydrologic properties of the multi-layered alluvial aquifer.  相似文献   

16.
The aquifer of the semi-arid Kairouan plain has been exploited for decades to supply the growing irrigated agriculture and the need of drinking water. In parallel, the major hydraulic works drastically changed the natural groundwater recharge processes. The continuous groundwater level drop observed since the 1970s naturally raises the question of groundwater storage sustainability. To date, hydrogeological studies focused on groundwater fluxes, but the total amount of groundwater stored in the aquifer system has never been fully estimated. This is the purpose of the present paper. A complete database of all available geological, hydrogeological and geophysical data was created to build a 3D lithology model. Then, the lithological units were combined with the hydraulic properties to estimate the groundwater storage. Over the 700 km2 of the modelled area, the estimated storage in 2013 was around 18?×?109 m3 (equivalent to 80 times the annual consumption of 2010) with a highly variable spatial distribution. In 45 years (1968–2013), 12% of the amount of groundwater stored in the aquifer has been depleted. According to these results, individual farms will face strong regional disparities for their access to groundwater in the near future.  相似文献   

17.
Groundwater in shallow unconsolidated sedimentary aquifers close to the Bornheim fault in the Lower Rhine Embayment (LRE), Germany, has relatively low δ2H and δ18O values in comparison to regional modern groundwater recharge, and 4He concentrations up to 1.7?×?10?4 cm3 (STP) g–1?±?2.2 % which is approximately four orders of magnitude higher than expected due to solubility equilibrium with the atmosphere. Groundwater age dating based on estimated in situ production and terrigenic flux of helium provides a groundwater residence time of ~107 years. Although fluid exchange between the deep basal aquifer system and the upper aquifer layers is generally impeded by confining clay layers and lignite, this study’s geochemical data suggest, for the first time, that deep circulating fluids penetrate shallow aquifers in the locality of fault zones, implying  that sub-vertical fluid flow occurs along faults in the LRE. However, large hydraulic-head gradients observed across many faults suggest that they act as barriers to lateral groundwater flow. Therefore, the geochemical data reported here also substantiate a conduit-barrier model of fault-zone hydrogeology in unconsolidated sedimentary deposits, as well as corroborating the concept that faults in unconsolidated aquifer systems can act as loci for hydraulic connectivity between deep and shallow aquifers. The implications of fluid flow along faults in sedimentary basins worldwide are far reaching and of particular concern for carbon capture and storage (CCS) programmes, impacts of deep shale gas recovery for shallow groundwater aquifers, and nuclear waste storage sites where fault zones could act as potential leakage pathways for hazardous fluids.  相似文献   

18.
Water-table elevation measurements and aquifer parameter estimates are rare in alpine settings because few wells exist in these environments. Alpine groundwater systems may be a primary source of recharge to regional groundwater flow systems. Handcart Gulch is an alpine watershed in Colorado, USA comprised of highly fractured Proterozoic metamorphic and igneous rocks with wells completed to various depths. Primary study objectives include determining hydrologic properties of shallow bedrock and surficial materials, developing a watershed water budget, and testing the consistency of measured hydrologic properties and water budget by constructing a simple model incorporating groundwater and surface water for water year 2005. Water enters the study area as precipitation and exits as discharge in the trunk stream or potential recharge for the deeper aquifer. Surficial infiltration rates ranged from 0.1–6.2×10?5 m/s. Discharge was estimated at 1.28×10?3 km3. Numerical modeling analysis of single-well aquifer tests predicted lower specific storage in crystalline bedrock than in ferricrete and colluvial material (6.7×10?5–2.0×10?3 l/m). Hydraulic conductivity in crystalline bedrock was significantly lower than in colluvial and alluvial material (4.3×10?9–2.0×10?4 m/s). Water budget results suggest that during normal precipitation and temperatures water is available to recharge the deeper groundwater flow system.  相似文献   

19.
20.

Hydraulic properties of coastal, urban aquifers vary spatially and temporally with the complex dynamics of their hydrogeology and the heterogeneity of ocean-influenced hydraulic processes. Traditional aquifer characterisation methods are expensive, time-consuming and represent a snapshot in time. Tidal subsurface analysis (TSA) can passively characterise subsurface processes and establish hydro-geomechanical properties from groundwater head time-series but is typically applied to individual wells inland. Presented here, TSA is applied to a network of 116 groundwater boreholes to spatially characterise confinement and specific storage across a coastal aquifer at city-scale in Cardiff (UK) using a 23-year high-frequency time-series dataset. The dataset comprises Earth, atmospheric and oceanic signals, with the analysis conducted in the time domain, by calculating barometric response functions (BRFs), and in the frequency domain (TSA). By examining the damping and attenuation of groundwater response to ocean tides (OT) with distance from the coast/rivers, a multi-borehole comparison of TSA with BRF shows this combination of analyses facilitates disentangling the influence of tidal signals and estimation of spatially distributed aquifer properties for non-OT-influenced boreholes. The time-series analysed covers a period pre- and post-impoundment of Cardiff’s rivers by a barrage, revealing the consequent reduction in subsurface OT signal propagation post-construction. The results indicate that a much higher degree of confined conditions exist across the aquifer than previously thought (specific storage?=?2.3 × 10?6 to 7.9 × 10?5 m?1), with implications for understanding aquifer recharge, and informing the best strategies for utilising groundwater and shallow geothermal resources.

  相似文献   

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

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