首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Coastal groundwater discharge (CGD) plays an important role in coastal hydrogeological systems as they are a water resource that needs to be managed, particularly in wetland areas. Despite its importance, identifying and monitoring CGD often presents physical and logistical constraints, restraining the application of more traditional submarine groundwater discharge surveying techniques. Here we investigate the capability of electrical resistivity imaging (ERI) in the Peníscola wetland (Mediterranean coast, Spain). ERI surveying made it possible to identify and delineate an ascending regional groundwater flow of thermal and Ra‐enriched groundwater converging with local flows and seawater intrusion. The continuous inputs of Ra‐rich groundwater have induced high activities of Ra isotopes and 222Rn into the marsh area, becoming among the highest previously reported in wetlands and coastal lagoons. Geoelectrical imaging enabled inferring focused upward discharging areas, leaking from the aquifer roof through a confining unit and culminating as spring pools nourishing the wetland system. Forward modelling over idealized subsurface configurations, borehole datasets, potentiometric records from standpipe piezometers, petrophysical analysis, and four natural and independent tracers (224Ra, 222Rn, temperature and salinity) permitted assessing the geoelectrical model and a derived hydrogeological pattern. The research highlights the potential of ERI to improve hydrogeological characterization of subsurface processes in complex contexts, with different converging flows. Additionally, a hydrogeological conceptual model for a groundwater‐fed coastal wetland was proposed, based on the integration of surveying datasets. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

2.
ABSTRACT

Integrated two-dimensional electrical resistivity imaging (ERI) and hydrochemical surveys were used to investigate the groundwater alluvial aquifer in Kuala Langat, Malaysia. The study in the Langat basin considered the thickness of the aquifer, the depth of the bedrock, the regions influenced by seawater intrusion, and the monitoring of water levels. The resistivity imaging results show that the upper layer consists of clay, while the second layer is an aquifer whose thickness varies mostly in the range of 10–30 m, and in some cases extends to 40 m. The bedrock depth varies from 30 to 65 m. The chemical analyses were carried out on groundwater samples from nine boreholes collected between 2008 and 2012. The analyses indicate that the total dissolved solids (TDS) exceed 1000 mg L-1 near the coastal area and are often less than 500 mg L-1 further inland. The ERI and hydrochemical analyses reveal that groundwater in the study area, especially towards the coast, is a mixture of brackish and fresh waters.
EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR M.D. Fidelibus  相似文献   

3.
Geophysical methods were applied for hydrogeological targets in many countries including Vietnam. This paper presents results of using complex geophysical techniques as well as 2D electrical resistivity imaging (ERI), vertical electrical sounding (VES), very low frequency (VLF), and seismic refraction for geological structure investigation for locating the aquifers and assessing the hydrogeological conditions for groundwater potential in industrial zones of North Hanoi, Vietnam. The locations of two aquifers are determined by their depth and thickness on the basis of resistivity and seismic velocity values which were proved by stratifications of three boreholes to 40–60 m of depth on the study area. There are connections from surface water to shallow aquifer by hydraulic windows, as follows from VLF data. The deeper aquifer can be considered as a potential groundwater supply, but the water level is descending in time, as shown by hydrological monitoring. However, with careful use and by reducing sources of pollution, groundwater can continue to be an important natural resource for future.  相似文献   

4.
Marine intrusion is the most serious problem facing the coastal Jorf shallow aquifer, located in south‐eastern Tunisia on the Mediterranean Sea. Jorf Aquifer is intensively exploited to supply the growing needs of agriculture and domestic sectors. This work proposes a multidisciplinary investigation, involving hydro‐geochemical, geoelectrical survey and geostatistical techniques for modelling the saltwater intrusion. For this purpose, 36 water samples were conducted and analysed. Electric conductivity, pH, total dissolved solids and major ions were measured and analysed. Pie and Durov Diagrams, Q‐mode hierarchical cluster and geostatistical analysis were considered to identify the main groundwater mineralization processes. Results revealed that the Na‐Cl‐Ca‐SO4 is the dominant water type suggesting that dissolution of halite and gypsum was the main mineralization source of groundwater in the central and southern part of study area. However, saltwater intrusion was shown to control groundwater quality essentially in coastal areas. Variographic analyses were used to select the variographic model that best fits the spatial development of apparent resistivity. Kriged apparent resistivity profiles showed an abnormal decrease of resistivity values in the coastal zone, implying highly saline water because of seawater intrusion. Apparent resistivity values also decrease considerably in the faulted areas, suggesting a contribution of faults to seawater intrusion. Finally, saltwater mixing ratio was computed for each sample, and a refined seawater intrusion map was developed. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

5.
Pumping wells are common in coastal aquifers affected by tides. Here we present analytical solutions of groundwater table or head variations during a constant rate pumping from a single, fully-penetrating well in coastal aquifer systems comprising an unconfined aquifer, a confined aquifer and semi-permeable layer between them. The unconfined aquifer terminates at the coastline (or river bank) and the other two layers extend under tidal water (sea or tidal river) for a certain distance L. Analytical solutions are derived for 11 reasonable combinations of different situations of the L-value (zero, finite, and infinite), of the middle layer’s permeability (semi-permeable and impermeable), of the boundary condition at the aquifer’s submarine terminal (Dirichlet describing direct connection with seawater and no-flow describing the existence of an impermeable capping), and of the tidal water body (sea and tidal river). Solutions are discussed with application examples in fitting field observations and parameter estimations.  相似文献   

6.
The southern coastal plain of Laizhou Bay, which is the area most seriously affected by salt water intrusion in north China, is a large alluvial depression, which represents one of the most important hydrogeological units in the coastal region of northern China. Chlorofluorocarbons (CFCs, including CFC‐11, CFC‐12 and CFC‐113) and tritium were used together for dating groundwater up to 50 years old in the study area. There are two cones of depression, caused by intensive over‐exploitation of fresh groundwater in the south and brine water in the north. The assigned CFC apparent ages for shallow groundwater range from 8 a to >50 a. A binary mixing model based on CFC‐113 and CFC‐12 concentrations in groundwater was used to estimate fractions of young and pre‐modern water in shallow aquifers and to identify groundwater mixing processes during saltwater intrusion. Discordance between concentrations of different CFC compounds indicate that shallow groundwater around the Changyi cone of depression is vulnerable to contamination. Pumping activities, CFC contamination, mixing and/or a large unsaturated zone thickness (e.g. >20 m) may be reasons for some groundwater containing CFCs without tritium. Saline intrusion mainly occurs because of large head gradients between fresh groundwater in the south and saline water bodies in the north, forming a wedge of saline water below/within fresh aquifer layers. Both CFC and tritium dates indicate that the majority of the saline water is from >50 a, with little or no modern seawater component. Based on the distribution of CFC apparent ages, tritium contents plus chemical and physical data, a conceptual model of groundwater flow along the investigated Changyi‐Xiaying transect has been developed to describe the hydrogeological processes. Three regimes are identified from south to north: (i) fresh groundwater zone, with a mixing fraction of 0.80–0.65 ‘young’ water calculated with the CFC binary mixing model (groundwater ages <34 a) and 1.9–7.8TU of tritium; (ii) mixing zone characterized by a mixing fraction of 0.05–0.65 young groundwater (ages of 23–44 a), accompanied by local vertical recharge and upward leakage of older groundwater; and (iii) salt water zone, mostly comprising waters with ages beyond the dating range of both CFCs and tritium. Some shallow groundwater in the north of the Changyi groundwater depression belongs to the >50a water group (iii), indicating slow velocity of groundwater circulation and possible drawing in of saline or deep groundwater that is tracer‐free. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

7.
Along the tropical coastline of the Great Barrier Reef (GBR) region, little is known to date about submarine groundwater discharge (SGD) into the near-shore ocean. In an oceanographic sense, SGD consists of freshwater flow from land as well as seawater circulated through sediments. Recent radiochemical and geophysical studies, using the tracer (222)Rn and apparent ground conductivity respectively, provide evidence for SGD to occur in a variety of hydrogeological settings. In this paper, a non-quantitative overview of different settings of SGD in the region is presented: (1) recirculation of seawater through animal burrows in mangrove forests, (2) freshwater SGD from unconfined aquifers as a narrow coastal fringe of freshwater along Wet Tropics beaches, (3) SGD from coastal dune systems in form of localised freshwater springs in the intertidal zone, (4) inner-shelf SGD from confined submarine aquifer systems comprised of riverine paleochannels incised into the shelf.  相似文献   

8.
Intensive pumping in urban coastal areas is a common threat to water resource quality due to seawater intrusion. In those areas where subsurface water resources are not usually used for human consumption or irrigation, intensive pumping is associated with other activities like the lowering of the water table necessary to support underground structures and building foundations. This activity also increases the likelihood of soil settlement that affects building stability and the corrosion of concrete structures due to groundwater salinity. Under these circumstances, the awareness of a certain municipality (Calonge, NE Spain) of the potential effects of groundwater withdrawal upon foundations has led to an integrated approach to anticipate seawater intrusion related to urban development. Geological mapping and correlation of borehole logs, electrical resistivity tomography, and hydrochemical data provide comprehensive knowledge of the geology and hydrogeology of the area and act as screening tools necessary to discern the influence of hydrological processes in coastal areas. Developing Strack's analytical solution, new comprehensive, dimensionless expressions are herein derived to determine the critical pumping rate necessary to prevent seawater intrusion, as well as to reproduce the evolution of the wedge toe and the water table stagnation point under different withdrawal rates. Furthermore, the Dupuit–Forchheimer well discharge formula allows the estimation of the effects of the water table lowering due to such critical pumping in the surrounding building foundations. Field data from the Calonge coastal plain illustrate this approach and provide assessment criteria for future urban development and planning. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

9.
Geophysical surveys have been carried out in two fjord delta aquifers in western Norway. Geophysical methods comprised ground penetrating radar (GPR), shallow reflection seismic and resistivity measurements. The purpose of the investigations was to evaluate hydrogeological properties of the fjord delta aquifers with emphasis on the possibilities of abstracting saline groundwater for use in fish farming. At Sunndalsøra, reflection seismic was helpful in the mapping of the base of the aquifer. Resistivity profiles revealed both the spatial and seasonal variability of seawater intrusion. Penetration depths of ground penetrating radar (GPR-D) revealed the spatial variation of the boundary of the transition zone between saline and fresh groundwater. Maps of GPR-D can be used to indicate direction and magnitude of mean hydraulic gradients and thus to locate optimal sites for both fresh and saline groundwater abstraction.In fjord delta aquifers, both rivers and abandoned river channels constitute important groundwater divides. Optimal locations for saline groundwater abstraction wells are in areas separated from the rivers by abandoned river channels. In areas between the rivers and the abandoned river channels, groundwater abstraction will result in decreasing salinity of the water.  相似文献   

10.
Five years of hydrogeological monitoring and field activities performed in the complex hydrogeological system of the Acque Albule basin (AAB) were conducted to define the hydrogeological setting, the relationship between deep and shallow aquifers and a conceptual groundwater flow model of this exploited area using conventional quantitative techniques. The basin, which is located close to Rome (Italy) on the west side of the Apennine chain and just north of the Colli Albani volcano, subsided after development of a north–south fault system (about 115 000 y bp). The AAB experiences intense hydrothermal activity, which has produced a large travertine deposit (80‐m thick). The travertine deposit constitutes a fractured aquifer that is the final destination of more than 5 m3 s‐1 of water and is strongly dewatered by quarry activities. The complex hydrogeology of this basin was investigated, revealing two main hydraulically connected aquifers, one thermalised and partly confined into the limestone bedrock and one unconfined in the travertine. The two aquifers are separated by a non‐continuous clayey aquiclude. The hydrogeological survey and geological characterisation contributed to the development of the groundwater flow conceptual model. Analysis and comparison of the monitored levels highlighted the pattern of flow between the deep and shallow parts of the flow system. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

11.
Vertical 2D slice laboratory experiments were carried out in homogenous and layered sand tanks to elucidate the effects of a highly permeable (coarse‐grained sand) interlayer on seawater intrusion and transport of contaminants to a coastal sea. Tidal fluctuations produced oscillations in the seawater–freshwater transition zone, fluctuations of the contaminant infiltration rate and a zigzag contaminant plume outline. The seawater wedge became discontinuous at the (vertical) edges of the interlayer because of increased lateral movement of the seawater–freshwater interface within the interlayer. The contaminant plume formed a tail within the interlayer depending on the tidal stage, and similar to the wedge, its movement was accentuated. A simple analytical model that neglected vertical flow reliably predicted steady‐state seawater intrusion into the coastal aquifer. Numerical modeling was used to gain insight into the groundwater hydrodynamics and contaminant migration. The numerical results confirmed the experimental findings, i.e. that a highly permeable interlayer can provide a rapid transit path for contaminants to reach the seaward boundary and that the interlayer amplifies the effects of tidal fluctuations, resulting in wider transition zones for the seawater wedge and contaminant plume. Numerical simulations further showed that, with increasing interlayer hydraulic conductivity, the maximum seawater intrusion distance inside the interlayer increases approximately linearly. For the fixed‐head contaminant injection condition used, the model showed that contaminant infiltration increases approximately logarithmically with increasing interlayer hydraulic conductivity (other factors held fixed). Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

12.
The concentrations of chlorofluorocarbons (CFC‐11, CFC‐12 and CFC‐113) and tritium (3H) content in groundwater were used to date groundwater age, delineate groundwater flow systems and estimate flow velocity in the Hohhot basin. The estimated young groundwater age is fallen in the bracket of 21 ~ 50 a and indicates the presence of two different age profiles and flow systems in the shallow groundwater system. Older age waters occur under the topographically low areas, where the aquifer is double‐layer aquifer system consisting of shallow unconfined‐semi‐confined aquifer and deep confined aquifer. This reflects long flow paths associated with regional flow. Groundwater (range from 21 to 34 years) in the north piedmont and east hilly areas, where the aquifer is a single‐layer aquifer consisting of alluvial fans, are typically younger than those in the low areas. The combination of CFCs dating with hydrogeological information indicates that both local and regional flow systems are present at the basin. The regional groundwater flow mainly flows from the north and east to the southwest, the local groundwater flow system occurs nearby the Hohhot city. The mean regional groundwater flow velocity of the shallow groundwater is estimated about 0.73 km/a. These findings can aid in refining hydrogeological conceptual model of the study area. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

13.
Studies investigating the effects of inland recharge on coastal groundwater dynamics were carried out typically in unconfined aquifers, with few in confined aquifers. This study focused on the groundwater dynamics in confined aquifers with seasonally sinusoidally fluctuated inland groundwater head and constant sea level by numerical simulations. It is known that the mixing zone (MZ) of saltwater wedge in response to the seasonal oscillations of inland groundwater head swings around the steady-state MZ. However, our simulation results indicate that even the most landward freshwater-saltwater interface over a year is seaward from the steady-state location when the hydraulic conductivity K is ≤10−4 m/s under certain boundary conditions with given parameter values. That is, seasonal oscillations of inland groundwater head may reduce seawater intrusion in confined coastal aquifers when K ≤ 10−4 m/s. Sensitivity analysis indicates that for aquifers of K ≤ 10−4 m/s, the larger the inland head fluctuation amplitude is, the less the seawater intrudes. This is probably due to the reason that the seawater intrusion time decreases with the increase of fluctuation amplitude when K ≤ 10−4 m/s. Numerical simulations demonstrate that seasonal inland groundwater head oscillations promote the annual averaged recirculated seawater discharge across the seaward boundary.  相似文献   

14.
Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.  相似文献   

15.
Effects of Recharge Wells and Flow Barriers on Seawater Intrusion   总被引：2，自引：0，他引：2  
The installation of recharge wells and subsurface flow barriers are among several strategies proposed to control seawater intrusion on coastal groundwater systems. In this study, we performed laboratory‐scale experiments and numerical simulations to determine the effects of the location and application of recharge wells, and of the location and penetration depth of flow barriers, on controlling seawater intrusion in unconfined coastal aquifers. We also compared the experimental results with existing analytical solutions. Our results showed that more effective saltwater repulsion is achieved when the recharge water is injected at the toe of the saltwater wedge. Point injection yields about the same repulsion compared with line injection from a screened well for the same recharge rate. Results for flow barriers showed that more effective saltwater repulsion is achieved with deeper barrier penetration and with barriers located closer to the coast. When the flow barrier is installed inland from the original toe position however, saltwater intrusion increases with deeper barrier penetration. Saltwater repulsion due to flow barrier installation was found to be linearly related to horizontal barrier location and a polynomial function of the barrier penetration depth.  相似文献   

16.
Subsurface dams are rather effective and used for the prevention of saltwater intrusion in coastal regions around the world. We carried out the laboratory experiments to investigate the elevation of saltwater wedge after the construction of subsurface dams. The elevation of saltwater wedge refers to the upward movement of the downstream saltwater wedge because the subsurface dams obstruct the regional groundwater flow and reduce the freshwater discharge. Consequently, the saltwater wedge cannot further extend in the longitudinal direction but rises in the vertical profile resulting in significant downstream aquifer salinization. In order to quantitatively address this issue, field-scale numerical simulations were conducted to explore the influence of various dam heights, distances, and hydraulic gradients on the elevation of saltwater wedge. Our investigation shows that the upward movement of the saltwater wedge and its areal extension in the vertical domain of the downstream aquifer become more severe with a higher dam and performed a great dependence on the freshwater discharge. Furthermore, the increase of the hydraulic gradient and the dam distance from the sea boundary leads to a more pronounced wedge elevation. This phenomenon comes from the variation of the freshwater discharge due to the modification of dam height, location, and hydraulic gradient. Large freshwater discharge can generate greater repulsive force to restrain the elevation of saltwater wedge. These conclusions provide theoretical references for the behaviour of the freshwater–seawater interface after the construction of subsurface dams and help optimize the design strategy to better utilize the coastal groundwater resources.  相似文献   

17.
The aquifer of Nador has suffered significant salinization due to seawater intrusion. It was strongly exploited during the 1980s and 1990s. A piezometric analysis in April 2012 showed the piezometric level to lie at 0 m a.s.l. over the plain; as a result, this aquifer is highly sensitive to the marine intrusion with an electrical conductivity of the groundwater in of exceeds 2500 μS/cm and so there are no abstractions for irrigation or drinking purpose from these sectors. The geoelectric study also showed the lateral variation in the electrical resistivity for two moments separated in time by more than 45 years. The fall in resistivity may be due to the encroachment of seawater into previously freshwater zones and/or infiltration during the era of pumped abstractions downstream. The resistivity surveys reveal two distinct sectors: the saturated aquifer in brackish and saltwater having resistivity values to 36-10 Ωm, which extends nearly 1600 m inland.  相似文献   

18.
Analytical models have been exhaustively used to study simple seawater intrusion problems and the sustainable management of groundwater resources in coastal aquifers because of its simplicity, easy implementation, and low computational cost. Most of these models are based on the sharp‐interface approximation and the Ghyben–Herzberg relation, and their governing equations are expressed in terms of a single potential theory to calculate critical pumping rates in a coastal pumping scenario. The Ghyben–Herzberg approach neglects mixing of fresh water and seawater and implicitly assumes that salt water remains static. Therefore, the results of the analytical solutions may be inaccurate and unacceptable for some real‐complex case studies. This paper provides insight into the validity of sharp‐interface models to deal with seawater intrusion in coastal aquifers, i.e. when they can be applied to obtain accurate enough results. For that purpose, this work compares sharp‐interface solutions, based on the Ghyben–Herzberg approach, with numerical three‐dimensional variable‐density flow simulations for a set of heterogeneous groundwater flow and mass transport parameters, and different scenarios of spatially distributed recharge values and spatial wells placement. The numerical experiment has been carried out in a 3D unconfined synthetic aquifer using the finite difference numerical code SEAWAT for solving the coupled partial differential equations of flow and density‐dependent transport. This paper finds under which situations the sharp‐interface solution gives good predictions in terms of seawater penetration, transition zone width and critical pumping rates. Additionally, the simulation runs indicate to which parameters and scenarios the results are more sensitive. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

19.
1D resistivity sounding and 2D resistivity imaging surveys were integrated with geological and hydrochemical data to assess the aquifer vulnerability and saltwater intrusion in the north of Nile Delta, Egypt. In the present study, the El-Gharbyia main drain was considered as a case study to map the sand bodies within the upper silt and clay aquitard. Twenty Schlumberger soundings and six 2D dipole-dipole profiles were executed along one profile close to the western side of the main drain. In addition, 14 groundwater samples and 4 surface water samples from the main drain were chemically analyzed to obtain the major and trace elements concentrations.The results from the resistivity and hydrochemical data were used to assess the protection of the groundwater aquifer and the potential risk of groundwater pollution. The inverted resistivities and thicknesses of the layers above the aquifer layer were used to estimate the integrated electrical conductivity (IEC) that can be used for quantification of aquifer vulnerability. According to the aquifer vulnerability assessment of an underlying sand aquifer, the southern part of the area is characterized by high vulnerability zone with slightly fresh to brackish groundwater and resistivity values of 11-23 Ω.m below the clay cap. The resistivity sections exhibit some sand bodies within the clay cap that lead to increase the recharging of surface waste water (650 mg/l salinity) and flushing the upper part of underlying saltwater aquifer. The region in the north has saltwater with resistivity less than 6 Ω.m and local vulnerable zones within the clay cap. The inverted 2D dipole-dipole profiles in the vulnerable zones, in combination with drilling information have allowed the identification of subsoil structure around the main drain that is highly affected by waste water.  相似文献   

20.
Abstract

Abstract Geophysical results obtained in the Rharb basin, Morocco are reported. Correlations between hydrogeological well logs reveal several water-bearing Plio-Quaternary units resting on a substrate of blue marls. Geo-electrical borehole analyses were interpreted using bi-logarithmic diagrams which indicate the permeable layers of the aquifer and also its basement. Resistivity data from NE–SW and NW–SE electrical sections allow definition of the permeable/impermeable levels, and identification of ?ditches? that may be favourable sectors for hydrogeological exploitation. Resistivity anomalies were investigated by analysing maps of resistivity at 400 and 1000 m AB. Anomalies identified in the Rharb basin are related to the thickening of the permeable bodies (sand, limestone, sandstone deposits). In the coastal zone (AB = 1000 m), there is a strong decrease of the resistivity gradient (35–10 Ω m), which is probably linked to marine intrusion. Electrical anomalies allow detection of the water-bearing zones notably in the western and southwestern parts of the Rharb basin.  相似文献   

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

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