首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
The Green–Ampt infiltration equation is an incomplete governing equation for rainfall infiltration due to the absence of an inertia term. The estimation of the capillary pressure head at the wetting front is difficult to determine. Thus, a major limitation of the Green–Ampt model is the constant, non‐zero surface ponding depth. This paper proposes an integrated rainfall infiltration model based on the Green–Ampt model and the SCS‐CN model. It achieves a complete governing equation for rainfall infiltration by momentum balance and the water budget based on the Green–Ampt assumption, and uses the curve number from the SCS‐CN method to calculate the initial abstraction, which is used as a basic parameter for the governing equation of the intensity of rainfall loss during the runoff period. The integrated rainfall infiltration model resolves the dilemma for capillary pressure head estimation, overcomes the limitation of constant, non‐zero surface ponding depth, and facilitates the calculation of runoff for individual flood simulations. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

2.
A procedure combining the Soil Conservation Service‐Curve Number (SCS‐CN) method and the Green–Ampt (GA) infiltration equation was recently developed to overcome some of the drawbacks of the classic SCS‐CN approach when estimating the volume of surface runoff at a sub‐daily time resolution. The rationale of this mixed procedure, named Curve Number for Green–Ampt (CN4GA), is to use the GA infiltration model to distribute the total volume of the net hyetograph (rainfall excess) provided by the SCS‐CN method over time. The initial abstraction and the total volume of rainfall given by the SCS‐CN method are used to identify the ponding time and to quantify the hydraulic conductivity parameter of the GA equation. In this paper, a sensitivity analysis of the mixed CN4GA parameters is presented with the aim to identify conditions where the mixed procedure can be effectively used within the Prediction in Ungauged Basin perspective. The effects exerted by changes in selected input parameters on the outputs are evaluated using rectangular and triangular synthetic hyetographs as well as 100 maximum annual storms selected from synthetic rainfall time series. When applied to extreme precipitation events, which are characterized by predominant peaks of rainfall, the CN4GA appears to be rather insensitive to the input hydraulic parameters of the soil, which is an interesting feature of the CN4GA approach and makes it an ideal candidate for the rainfall excess estimation at sub‐daily temporal resolution at ungauged sites. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

3.
To investigate the source, flow paths, and chemistry of rich resources of high‐quality, shallow groundwater in the alluvial fan between the Tedori and Sai rivers in central Japan, we analysed stable isotope ratios of H, O, and Sr and concentrations of major dissolved ions and trace elements in groundwater, river water, and paddy water. The 87Sr/86Sr ratios of the groundwater are related to near‐surface geology: groundwater in sediment from the Tedori River has high 87Sr/86Sr ratios (>0.711), whereas that from the Sai River in the north of the fan has low 87Sr/86Sr ratios (<0.711). δ2H and δ18O values and 87Sr/86Sr ratios indicate that groundwater in the central and southern fans is recharged by the Tedori River, whereas recharge in the north is from the Sai River. Mg2+, Ca2+, Sr2+, HCO3?, and SO42? concentrations and δ2H and δ18O values in the groundwater are high in the central fan and, except for the northern area, tend to increase with distance from the Tedori River. There are linear relationships between 87Sr/86Sr ratio and the reciprocal concentrations of Sr2+, Mg2+, and Ca2+. These geochemical characteristics suggest that as groundwater recharged from the Tedori River flows towards the central fan, it mixes with waters derived from precipitation and paddy water that have become enriched in these components during downward infiltration. These results are consistent with our hydrological analysis and numerical simulation of groundwater flow, thus verifying the validity of the model we used in our simulation of groundwater flow. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

4.
Abstract

The potential impacts of future climate change on the evolution of groundwater recharge are examined at a local scale for a 546-km2 watershed in eastern Canada. Recharge is estimated using the infiltration model Hydrologic Evaluation of Landfill Performance (HELP), with inputs derived from five climate runs generated by a regional climate model in combination with the A2 greenhouse gas emissions scenario. The model runs project an increase in annual recharge over the 2041–2070 period. On a seasonal basis, however, a marked decrease in recharge during the summer and a marked increase during the winter are observed. The results suggest that increased evapotranspiration resulting from higher temperatures does not offset the large increase in winter infiltration. In terms of individual water budget components, clear differences are obtained for the different climate change scenarios. Monthly recharge values are also found to be quite variable, even for a given climate scenario. These findings are compared with results from two regional-scale studies.
Editor D. Koutsoyiannis; Associate editor M. Besbes  相似文献   

5.
The groundwater in shallow loess aquifers in high mountain–hills in the western Loess Plateau in China is almost the sole water resource for local residents. However, the question about how the loess groundwater naturally circulates in these high mountain–hills, characterized by low precipitation and high potential evaporation, remains unclear. The objectives of this study are to evaluate the application of hydrogen and oxygen isotopes to (1) examine temporal variations of the isotopic composition of precipitation and shallow groundwater and (2) uncover the mechanism of groundwater recharge in high mountain–hills. Results from 2 years of monitoring data show a difference in the stable isotopes for groundwater and local precipitation between the winter and summer periods. Similar to precipitation, stable isotopes in groundwater are observed to be depleted in winter and enriched in summer, particularly in oxygen isotope. A prominent characteristic is that H and O isotopes of groundwater show a very clear response to strong precipitation in the rainy season in 2013. The results highlight that local precipitation is the likely recharge source for groundwater in shallow loess aquifers. Annual recharge from local precipitation maintains the groundwater resource in the shallower loess aquifer. The mechanisms governing shallow loess groundwater recharge in high mountain–hills were evaluated. In addition to possible vertical slow percolation of soil water through the unsaturated zone, rapid groundwater recharge mechanisms have been identified as temporal preferential infiltration through sinkholes, slip surface or landslide surface and through the interface of loess layer and palaeo‐soils. Most groundwater can be recharged after a heavy rainy season. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

6.
It is often assumed that the net groundwater flow direction is towards the channel in headwater streams in humid climates, with magnitudes dependent on flow state. However, studies that characterize stream–groundwater interactions in ephemeral and intermittent streams in humid landscapes remain sparse. Here, we examined seasonally driven stream–groundwater interactions in response to temporary streamflow on the basis of field observations of streamflow and groundwater on an adjacent hillslope. The direction of hydraulic head gradients between the stream and groundwater shifted seasonally. The stream gained water (head gradients were towards the stream) when storage state was high. During this period, streamflow was persistent. The stream lost water to the groundwater system (head gradients were away from the stream) when storage state was low. During this period, streamflow only occurred in response to precipitation events, and head gradients remained predominantly away from the stream during events. This suggested that mechanisms other than deep groundwater contributions produced run‐off when storage was low, such as surface and perched subsurface flowpaths above the water table. Analysis of the annual water balance for the study period showed that the residual between precipitation inputs and streamflow and evapotranspiration outputs, which were attributed to the loss of water to the deeper, regional groundwater system, was similar in magnitude to streamflow. This, coupled with results that showed bidirectionality in stream–groundwater head gradients, indicated that headwaters composed of temporary (e.g., ephemeral and intermittent) streams can be important focal areas for regional groundwater recharge, and both contribute to and receive water, solutes, and materials from the groundwater system.  相似文献   

7.
ABSTRACT

Many oases are experiencing severe groundwater depletion due to increased population, expanding agriculture and economic development. For sustainable development, quantifying groundwater recharge resources are fundamentally important. In this study, stable isotope techniques were employed to identify recharge sources of groundwater and quantitatively evaluate their contribution ratios in the Dunhuang Oasis, northwest China. Our findings indicate that heavy isotopes in shallow groundwater are more negative than those in deep groundwater, which is attributed to shallow groundwater that was modern and deep groundwater that was old. Irrigated return water and lateral groundwater flow from the Qilian Mountains are considered as the two main sources of shallow groundwater, accounting for 35% and 65% of the total recharge, respectively. Thus, as the main groundwater source of the Dunhuang Oasis, the Qilian Mountain Front should be protected against over-exploitation. Our results provide not only fundamental knowledge for groundwater management of aquifers of the Oasis, but also valuable water management information for other similar arid oases worldwide.  相似文献   

8.
In the shallow groundwater areas of the North China Plain (NCP), precipitation infiltration and evapotranspiration in the vertical direction are the main processes of the water cycle, in which the unsaturated zone plays an important role in the transformation process between precipitation and groundwater. In this paper, two typical sites in Cangzhou (CZ) and Hengshui (HS) of Hebei province with shallow water tables were selected to analyse the relationship among precipitation, soil water and groundwater. At each site, precipitation, soil water at depths 10, 20, 30, 50, 70, 100, 150, 200, 300 cm, and groundwater were sampled to analyse the stable isotope compositions of hydrogen and oxygen. The soil water potentials at the corresponding depths were observed. Although the climates at the two sites are similar, there are some differences in the infiltration process, soil water movement and groundwater recharge sources. Evaporation occurred at the upper depths, which led to the decrease of soil potential and the enrichment of heavy isotopes. At the CZ site, precipitation infiltrated with piston mode, and an obvious mixture effect existed during the infiltration process. Preferential flow may exist in the soil above 100 cm depth. However, at the HS site soil water moved in piston mode, and groundwater was mainly recharged by precipitation. When precipitation recharged the groundwater it experienced a strong evaporation effect. The results of the soil water movement mechanism provides the transformation relationship among precipitation, soil water and groundwater in the middle and eastern NCP. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

9.
Uruguay has stimulated the development of its forest sector since the promulgation of Forest Law N° 15 939 in December of 1987. Nevertheless, the substitution of natural grasslands with forest plantations for industrial use has raised concerns regarding hydrological processes of groundwater recharge and water consumption involving evapotranspiration. The purpose of this study is to assess the effects of this substitution approach on water resources. Input data were collected from two small experimental watersheds of roughly 100–200 hectares located in western Uruguay. The watersheds are characterized by Eucalyptus Globulus ssp. Maidenni and natural grasslands for cattle use. Total rainfall, stream discharge, rainfall redistribution, soil water content and groundwater level data were collected. Groundwater recharge was estimated from water table fluctuations and from groundwater contributions to base flows. Seasonal and annual water budgets were computed from October of 2006 to September of 2014 to evaluate changes in the hydrological processes. The data show a decrease in annual specific discharge of roughly 17% for mean hydrological years and no conclusive effects on annual groundwater recharge in the forested watershed relative to the reference pasture watershed. Reduced annual specific discharge is equivalent to the mean annual interception. The computed actual annual evapotranspiration is consistent with international catchment measurements. Reduction rates vary seasonally and according to accumulated rainfall and its temporary distribution. The degree of specific discharge decline is particularly high for drier autumns and winters (32 to 28%) when the corresponding rainfall varies from 275 to 400 mm. These results are of relevance for water resources management efforts, as water uses downstream can be affected. These findings, based on a study period dominated by anomalous wet springs and summers and by dry autumns and winters, oppose earlier results based on 34 years of rainfall and discharge data drawn from Uruguayan large basins. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

10.
The tenets of sustainable development are often ignored by growing cities in developing countries. Salt Lake City, located on the eastern fringes of Kolkata, despite being a planned township, has failed to manage its water resources in a sustainable way. Results from the study show that the water supplied at the households is hard in nature, saline, and replete with dissolved solids. The problem is further aggravated due to over-extraction of groundwater resources, resulting in an intense fall of the groundwater level in the area. This paper attempts to suggest measures for sustainable water resource management, with particular emphasis on rainwater harvesting.  相似文献   

11.
The objectives of this study are (1) to understand the subsurface hydrology in the Aurku area, Chiayi County, southern Taiwan, and (2) to determine the interaction between the manmade lake and groundwater level through the recharge produced by infiltration by on‐site investigation and laboratory sand tank simulation. The manmade lake was selected as the field site for groundwater recharge effect so as to assess the role of infiltration from the aquaculture ponds in this area. These results can be used as reference for future application of constructing a series of manmade lakes. The field experiment was performed to measure the infiltration rate of the manmade lake by using the water balance method and double‐ring infiltration test. The results demonstrated that the manmade lake had helped the recharge of the groundwater. Raising or maintaining a higher water level of the manmade lake can promote higher infiltration. When the groundwater level is equal to or higher than the bottom of the manmade lake, infiltration will slow or cease. The field experiment and laboratory sand tank simulation demonstrated that the infiltration rate increased with the higher storage depth of the manmade lake. The laboratory simulation also indicated that while the groundwater level was lower than the bottom of manmade lake (i.e. the reference level) and the initial water depth (3 cm) was equal to or greater than 50% of the full water storage depth, the infiltration depth increased with time. However, the infiltration depth would be very small or nearly zero when the groundwater level was higher than the bottom of the manmade lake. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

12.
The main objective of this research was to analyse the effect of soil management on soil sealing and on soil water content under contrasting tillage practices and its influence on corn yield. The experimental research was carried out in a field cultivated with irrigated corn differentiated into three zones representing a gradient of soil texture (Z1, Z2, and Z3, i.e., increasingly coarser). Two plots under different soil management practices (conventional intensive tillage, CT, and no‐tillage, NT) were selected in each zone. The susceptibility to sealing of each soil and the steady infiltration rates were evaluated in the laboratory subjecting the soils to rainfall simulation applied at an intensity of 25 mm h?1. In addition, soil porosity under each treatment was quantified. Soil water content (0–90 cm depth) was determined gravimetrically at the beginning and the end of the growing cycle and at the surface (0–5 cm) during three growing seasons and continuously at two depths (5–15 and 50–60 cm) during the last growing cycle. Soil water content was simulated using the SIMPEL model, which was calibrated for the experimental conditions. Corn yield and above‐ground biomass were also analysed. Significant differences in soil sealing among zones, with decreasing soil sealing for coarser textures, and treatments were observed with infiltration rates that were near twice in NT than in CT, being the effect of soil cover significant in the reduction of soil detachment and soil losses. NT showed higher soil water content than CT, especially in the surface layers. Above‐ground biomass production was smaller in CT than in NT, and in the areas with higher sealing susceptibility was 30% to 45% smaller than in other zones, reaching the smallest values in Z1. A similar reduction in corn yield was observed between treatments being smaller in CT than in NT. No‐tillage has been confirmed as an effective technique that benefits soil physical properties as well as crop yields in relation to CT, being its impact greater in soils susceptible to sealing.  相似文献   

13.
Marathwada Agricultural University, Pharbani, has developed about 560 hectares of Wagarwadi watershed in Pharbani district since 1987. Groundwater monitoring on 16 observations wells at weekly intervals commenced in January 1992, and rainfall and pan evaporation has been measured daily at a hydrometeorological station situated in the nearby university campus. Aquifer parameters, namely, transmissivity and specific yield, have been estimated by carrying out a pumping test on a large diameter well. Groundwater recharge resulting from rainfall has been estimated using a water balance model of the soil moisture zone considering soil zone thickness and crops grown. The SCS (Soil Conservation Service) curve number method was used for surface runoff estimation. The groundwater flow model has been constructed using the nested squares, finite difference method. Nested square meshes of sizes 160 m×160 m and 80 m×80 m have been used and the steady-state condition of aquifer system was simulated in the model assuming the June 1992 water level configuration under equilibrium conditions. The model has been calibrated for transient conditions incorporating additional seepage from the water harvesting structures and their contribution to the groundwater regime has been assessed. © 1998 John Wiley & Sons, Ltd.  相似文献   

14.
15.
Environmental isotopes (2H or D, 18O, 3H), along with geology, hydrochemistry and in situ physicochemical parameters (EC, T, DO, pH) were employed to study surface water (reservoir, lake)–groundwater (spring) relationships at (1) Nagewadi, a minor irrigation project in the State of Maharashtra, Western India; (2) Kanhirapuzha reservoir in the State of Kerala, Southern India and (3) Ghatghar Pumped Storage Hydroelectric Project in the State of Maharashtra, Western India for the purpose of understanding the seepage/leakage and its associated problems. The studies concluded that the springs found downstream of the Nagewadi project originate from the reservoir and not from the abutments or shallow aquifers. The Kanhirapuzha reservoir receives a substantial base‐flow component compared to riverine inputs. The reason for the water‐logging problem at a nearby downstream village during the non‐summer periods is due to the change in the upstream groundwater flow direction under reservoir filling conditions and is not due to reservoir leakage. Most of the springs in the approach tunnel to the underground power house of the Ghatghar Project originate from the lower reservoir and not from the upper reservoir or the overburden rock matrix. The above case studies illustrate the diversity of environmental isotope applications in surface water (reservoir, lake)–groundwater (spring) relationships related to sustainability of hydro‐projects. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

16.
Recharge varies spatially and temporally as it depends on a wide variety of factors (e.g. vegetation, precipitation, climate, topography, geology, and soil type), making it one of the most difficult, complex, and uncertain hydrologic parameters to quantify. Despite its inherent variability, groundwater modellers, planners, and policy makers often ignore recharge variability and assume a single average recharge value for an entire watershed. Relatively few attempts have been made to quantify or incorporate spatial and temporal recharge variability into water resource planning or groundwater modelling efforts. In this study, a simple, daily soil–water balance model was developed and used to estimate the spatial and temporal distribution of groundwater recharge of the Trout Lake basin of northern Wisconsin for 1996–2000 as a means to quantify recharge variability. For the 5 years of study, annual recharge varied spatially by as much as 18 cm across the basin; vegetation was the predominant control on this variability. Recharge also varied temporally with a threefold annual difference over the 5‐year period. Intra‐annually, recharge was limited to a few isolated events each year and exhibited a distinct seasonal pattern. The results suggest that ignoring recharge variability may not only be inappropriate, but also, depending on the application, may invalidate model results and predictions for regional and local water budget calculations, water resource management, nutrient cycling, and contaminant transport studies. Recharge is spatially and temporally variable, and should be modelled as such. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

17.
Abstract

Semi-arid coastal zones often suffer water-stress, as water demand is high and markedly seasonal, due to agriculture and tourism. Driven by scarcity of surface water, the communities in semi-arid coastal regions turn to aquifers as prime water source; but intensive exploitation of coastal aquifers causes seawater intrusion, which degrades the quality of groundwater. The cost-efficient and sustainable development of coastal aquifers can be achieved through a holistic management scheme which combines two non-traditional water sources: (a) saltwater, to be treated to the desired quality, and (b) wastewater, to be re-claimed to augment aquifer recharge for control of seawater intrusion, and also to meet certain demands. This management scheme is based on the idea that it is cost-advantageous to: (i) desalt brackish groundwater, instead of seawater, as the former requires far less energy, and (ii) to re‐use wastewater at only the differential cost to any treatment already practiced. In this paper, we present the general framework of the proposed management scheme, and a decision aid tool (DAT) which has been developed to assist decision makers to explore the scheme's decision space. The DAT uses cost as optimization criterion to screen various management scenarios, via modelling of the dynamic natural-engineered system behaviour, and identifies those cost-efficient ones that meet the water demand and achieve aquifer protection.

Citation Koussis, A. D., Georgopoulou, E., Kotronarou, A., Lalas, D. P., Restrepo, P., Destouni, G., Prieto, C., Rodriguez, J. J., Rodriguez-Mirasol, J., Cordero, T. & Gomez-Gotor, A. (2010 Koussis, A. D., Georgopoulou, E., Kotronarou, A., Mazi, K., Restrepo, P., Destouni, G., Prieto, C., Rodriguez, J. J., Rodriguez-Mirasol, J., Cordero, T., Schwartz, J., Ioannou, C., Georgiou, A. and Zacharias, I. 2010. Cost-efficient management of coastal aquifers in water-stressed regions via recharge with treated wastewater and desalination of brackish groundwater: application to the Akrotiri basin and aquifer, Cyprus. Hydrol. Sci. J, 55(7): 12341245. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]) Cost-efficient management of coastal aquifers via recharge with treated wastewater and desalination of brackish groundwater: general framework. Hydrol. Sci. J. 55(7),1217–1233.  相似文献   

18.
In semi‐arid and arid river basins, understanding the connectivity between rivers and alluvial aquifers is one of the key challenges for the management of groundwater resources. The type of connection present (gaining, losing‐connected, transitional and losing‐disconnected) was assessed at 12 sites along six Murray–Darling Basin river reaches. The assessments were made by measuring the hydraulic head in the riparian zone near the rivers to evaluate if the water tables intersected the riverbeds and by measuring fluid pressure (ψ) in the riverbeds. The rationale for the latter was that ψ will always be greater than or equal to zero under connected conditions (either losing or gaining) and always lesser than or equal to zero under losing‐disconnected conditions. A mixture of losing‐disconnected, losing‐connected and gaining conditions was found among the 12 sites. The losing‐disconnected sites all had a riverbed with a lower hydraulic conductivity than the underlying aquifer, usually in the form of a silty clay or clay unit 0.5–2 m in thickness. The riparian water tables were 6 to 25 m below riverbed level at the losing‐disconnected sites but never lower than 1 m below riverbed level at the losing‐connected ones. The contrast in water table depth between connected and disconnected sites was attributed to the conditions at the time of the study, when a severe regional drought had generated a widespread decline in regional water tables. This decline was apparently compensated near losing‐connected rivers by increased infiltration rates, while the decline could not be compensated at the losing‐disconnected rivers because the infiltration rates were already maximal there. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

19.
A simple conceptual semi‐distributed modelling approach for assessing the impacts of climate change on direct groundwater recharge in a humid tropical river basin is investigated. The study area is the Chaliyar river basin in the state of Kerala, India. Many factors affecting future groundwater recharge include decrease or increase in precipitation and temperature regimes, coastal flooding, urbanization and changes in land use. The model is based on the water‐balance concept and links the atmospheric and hydrogeologic parameters to different hydrologic processes. It estimates daily water‐table fluctuation and is calibrated and validated using 10 years of data. Data for the first 6 years (2000 to 2005) is used for model calibration, and data for the remaining four years (2006 to 2009) is used for validation. For assessing the impact of predicted climate change on groundwater recharge during the period 2071–2100, temperature and precipitation data in two post climate change scenarios, A2 and B2, were predicted using the Regional Climate Model (RCM), PRECIS (Providing Regional Climates for Impact Studies). These data were then corrected for biases and used in a hydrologic model to predict groundwater recharge in the post climate change scenario. Due to lack of reliable data and proper knowledge as to the magnitude and extent of future climatic changes, it may not be possible to include all the possible effects quantitatively in groundwater recharge modelling. However, the study presents a scientific method to assess the impact of predicted climate change on groundwater recharge and would help engineers, hydrologists, administrators and planners to devise strategies for the efficient use as well as conservation of freshwater resources. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

20.
Artificial recharge of groundwater is an increasingly important method for augmenting groundwater supply and can have a positive or negative influence on the quality of water resources. We instrumented a managed aquifer recharge (MAR) pond in central coastal California to assess how patterns of infiltration and recharge affect the load of nitrate delivered to the underlying aquifer. The concentration of nitrate in infiltrating water consistently decreased during passage through the first metre of subsurface soils. Enrichment of 18O and 15 N in the residual nitrate in infiltrating water proceeded in a ratio of 1:2, indicating that denitrification plays a significant role in the quantitative reduction of nutrients exported during infiltration through shallow soils. The extent and rate of nitrate removal was spatially and temporally variable across the bottom of the recharge pond, with 30% to 60% of the nitrate load being removed over the first 6 weeks of managed aquifer recharge operation. During the period of highest N loading to the system, when the average infiltration rate was > 1 m/day, the recharge pond achieved a load reduction efficiency of 7 kg NO3?‐N/day/ha, which compares favourably to nitrate load reductions achieved by treatment wetlands. Groundwater mounding and water composition below the recharge pond suggest that recharge and subsequent lateral transport occur heterogeneously in the underlying aquifer. Nitrate concentrations in the aquifer following infiltration were lowered primarily by dilution, with little evidence for additional denitrification occurring in the aquifer in comparison to high rates documented during shallow infiltration. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

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

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