首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
This study is about use of spatially distributed rain in physically based hydrological models. In recent years, spatially distributed radar rainfall data have become available. The distributed radar rain is used to precisely model hydrologic processes and it is more realistic than the past practice of distribution methods like Thiessen polygons. Radar provides a highly accurate spatial distribution of rainfall and greatly improves the basin average rainfall estimates. However, quantification of the exact amount of rainfall from radar observation is relatively difficult. Thus, the fundamental idea of this study is to apply hourly gauge and radar rainfall data in a distributed hydrological model to simulate hydrological parameters. Hence the comparison is made between the outcomes of the WetSpa model from radar rainfall distribution and gauge rainfall distributed by the Thiessen polygon technique. The comparative plots of the hydrograph and the results of hydrological components such as evapotranspiration, surface runoff, soil moisture, recharge and interflow, reflect the spatially distributed radar input performing well for model outflow simulation.
EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR F. Pappenberger  相似文献   

2.
ABSTRACT

Multisource rainfall products can be used to overcome the absence of gauged precipitation data for hydrological applications. This study aims to evaluate rainfall estimates from the Chinese S-band weather radar (CINRAD-SA), operational raingauges, multiple satellites (CMORPH, ERA-Interim, GPM, TRMM-3B42RT) and the merged satellite–gauge rainfall products, CMORPH-GC, as inputs to a calibrated probability distribution model (PDM) on the Qinhuai River Basin in Nanjing, China. The Qinhuai is a middle-sized catchment with an area of 799 km2. All sources used in this study are capable of recording rainfall at high spatial and temporal resolution (3 h). The discrepancies between satellite and radar data are analysed by statistical comparison with raingauge data. The streamflow simulation results from three flood events suggest that rainfall estimates using CMORPH-GC, TRMM-3B42RT and S-band radar are more accurate than those using the other rainfall sources. These findings indicate the potential to use satellite and radar data as alternatives to raingauge data in hydrological applications for ungauged or poorly gauged basins.  相似文献   

3.
Guoqiang Wang  Zongxue Xu 《水文研究》2011,25(16):2506-2517
A grid‐based distributed hydrological model, PDTank model, is used to simulate hydrological processes in the upper Tone River catchment. The Tone River catchment often suffers from heavy rainfall events during the typhoon seasons. The reservoirs located in the catchment play an important role in flood regulation. Through the coupling of the PDTank model and a reservoir module that combines the storage function and operation function, the PDTank model is used for flood forecasting in this study. By comparing the hydrographs simulated using gauging and radar rainfall data, it is found that the spatial variability of rainfall is an important factor for flood simulation and the accuracy of the hydrographs simulated using radar rainfall data is slightly improved. The simulation of the typhoon flood event numbered No. 9 shows that the reservoirs in the catchment attenuate the peak flood discharge by 423·3 m3/s and validates the potential applicability of the distributed hydrological model on the assessment of function of reservoirs for flood control during typhoon seasons. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

4.
The emergence of regional and global satellite‐based rainfall products with high spatial and temporal resolution has opened up new large‐scale hydrological applications in data‐sparse or ungauged catchments. Particularly, distributed hydrological models can benefit from the good spatial coverage and distributed nature of satellite‐based rainfall estimates (SRFE). In this study, five SRFEs with temporal resolution of 24 h and spatial resolution between 8 and 27 km have been evaluated through their predictive capability in a distributed hydrological model of the Senegal River basin in West Africa. The main advantage of this evaluation methodology is the integration of the rainfall model input in time and space when evaluated at the sub‐catchment scale. An initial data analysis revealed significant biases in the SRFE products and large variations in rainfall amounts between SRFEs, although the spatial patterns were similar. The results showed that the Climate Prediction Center/Famine Early Warning System (CPC‐FEWS) and cold cloud duration (CCD) products, which are partly based on rain gauge data and produced specifically for the African continent, performed better in the modelling context than the global SRFEs, Climate Prediction Center MORPHing technique (CMORPH), Tropical Rainfall Measuring Mission (TRMM) and Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN). The best performing SRFE, CPC‐FEWS, produced good results with values of R2NS between 0·84 and 0·87 after bias correction and model recalibration. This was comparable to model simulations based on traditional rain gauge data. The study highlights the need for input specific calibration of hydrological models, since major differences were observed in model performances even when all SRFEs were scaled to the same mean rainfall amounts. This is mainly attributed to differences in temporal dynamics between products. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

5.
This paper provides a comparison of gauge and radar precipitation data sources during an extreme hydrological event. November–December 2006 was selected as a time period of intense rainfall and large river flows for the Severn Uplands, an upland catchment in the United Kingdom. A comparison between gauge and radar precipitation time‐series records for the event indicated discrepancies between data sources, particularly in areas of higher elevation. The HEC‐HMS rainfall‐runoff model was selected to assess the accuracy of the precipitation to simulate river flows for the extreme event. Gauge, radar and gauge‐corrected radar rainfall were used as model inputs. Universal cokriging was used to geostatistically interpolate gauge data with radar and elevation data as covariates. This interpolated layer was used to calculate the mean‐field bias and correct the radar composites. Results indicated that gauge‐ and gauge‐corrected radar‐driven models replicated flows adequately for the extreme event. Gauge‐corrected flow predictions produced an increase in flow prediction accuracy when compared with the raw radar, yet predictions were comparative in accuracy to those using the interpolated gauge network. Subsequent investigations suggested this was due to an adequate spatial and temporal resolution of the precipitation gauge network within the Severn Uplands. Results suggested that the six rain gauges could adequately represent precipitation variability of the Severn Uplands to predict flows at an approximately equal accuracy to that obtained by radar. Temporally, radar produced an increase in flow prediction accuracy in mountainous reaches once the gauge time step was in excessive of an hourly interval. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

6.
Abstract

Abstract After the destructive flood in 1998, the Chinese government planned to build national weather radar networks and to use radar data for real-time flood forecasting. Hence, coupling of weather radar rainfall data and a hydrological (Xinanjiang) model became an important issue. The present study reports on experience in such coupling at the Shiguanhe watershed. After having corrected the radar reflectivity and the attenuation data, the weather radar rainfall was estimated and then corrected in real time using a Kalman filter. In general, the precipitation estimated from weather radar is reasonably accurate in most of the catchment investigated, after corrections as above. Compared to the results simulated by raingauge data, the simulations based on the weather radar data are of similar accuracy. Present research results show that rainfall estimated from the weather radar, the radar data correction method, the method of coupling, and the Xinanjiang model lend themselves well to application in operational real-time flood forecasting.  相似文献   

7.
Images from satellite platforms are a valid aid in order to obtain distributed information about hydrological surface states and parameters needed in calibration and validation of the water balance and flood forecasting. Remotely sensed data are easily available on large areas and with a frequency compatible with land cover changes. In this paper, remotely sensed images from different types of sensor have been utilized as a support to the calibration of the distributed hydrological model MOBIDIC, currently used in the experimental system of flood forecasting of the Arno River Basin Authority. Six radar images from ERS‐2 synthetic aperture radar (SAR) sensors (three for summer 2002 and three for spring–summer 2003) have been utilized and a relationship between soil saturation indexes and backscatter coefficient from SAR images has been investigated. Analysis has been performed only on pixels with meagre or no vegetation cover, in order to legitimize the assumption that water content of the soil is the main variable that influences the backscatter coefficient. Such pixels have been obtained by considering vegetation indexes (NDVI) and land cover maps produced by optical sensors (Landsat‐ETM). In order to calibrate the soil moisture model based on information provided by SAR images, an optimization algorithm has been utilized to minimize the regression error between saturation indexes from model and SAR data and error between measured and modelled discharge flows. Utilizing this procedure, model parameters that rule soil moisture fluxes have been calibrated, obtaining not only a good match with remotely sensed data, but also an enhancement of model performance in flow prediction with respect to a previous calibration with river discharge data only. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

8.
Multi‐step ahead inflow forecasting has a critical role to play in reservoir operation and management in Taiwan during typhoons as statutory legislation requires a minimum of 3‐h warning to be issued before any reservoir releases are made. However, the complex spatial and temporal heterogeneity of typhoon rainfall, coupled with a remote and mountainous physiographic context, makes the development of real‐time rainfall‐runoff models that can accurately predict reservoir inflow several hours ahead of time challenging. Consequently, there is an urgent, operational requirement for models that can enhance reservoir inflow prediction at forecast horizons of more than 3 h. In this paper, we develop a novel semi‐distributed, data‐driven, rainfall‐runoff model for the Shihmen catchment, north Taiwan. A suite of Adaptive Network‐based Fuzzy Inference System solutions is created using various combinations of autoregressive, spatially lumped radar and point‐based rain gauge predictors. Different levels of spatially aggregated radar‐derived rainfall data are used to generate 4, 8 and 12 sub‐catchment input drivers. In general, the semi‐distributed radar rainfall models outperform their less complex counterparts in predictions of reservoir inflow at lead times greater than 3 h. Performance is found to be optimal when spatial aggregation is restricted to four sub‐catchments, with up to 30% improvements in the performance over lumped and point‐based models being evident at 5‐h lead times. The potential benefits of applying semi‐distributed, data‐driven models in reservoir inflow modelling specifically, and hydrological modelling more generally, are thus demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

9.
Australian arid zone ephemeral rivers are typically unregulated and maintain a high level of biodiversity and ecological health. Understanding the ecosystem functions of these rivers requires an understanding of their hydrology. These rivers are typified by highly variable hydrological regimes and a paucity, often a complete absence, of hydrological data to describe these flow regimes. A daily time‐step, grid‐based, conceptual rainfall–runoff model was developed for the previously uninstrumented Neales River in the arid zone of northern South Australia. Hourly, logged stage data provided a record of stream‐flow events in the river system. In conjunction with opportunistic gaugings of stream‐flow events, these data were used in the calibration of the model. The poorly constrained spatial variability of rainfall distribution and catchment characteristics (e.g. storage depths) limited the accuracy of the model in replicating the absolute magnitudes and volumes of stream‐flow events. In particular, small but ecologically important flow events were poorly modelled. Model performance was improved by the application of catchment‐wide processes replicating quick runoff from high intensity rainfall and improving the area inundated versus discharge relationship in the channel sections of the model. Representing areas of high and low soil moisture storage depths in the hillslope areas of the catchment also improved the model performance. The need for some explicit representation of the spatial variability of catchment characteristics (e.g. channel/floodplain, low storage hillslope and high storage hillslope) to effectively model the range of stream‐flow events makes the development of relatively complex rainfall–runoff models necessary for multisite ecological studies in large, ungauged arid zone catchments. Grid‐based conceptual models provide a good balance between providing the capacity to easily define land types with differing rainfall–runoff responses, flexibility in defining data output points and a parsimonious water‐balance–routing model. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

10.
Rivers draining the southern Himalaya provide most of the water supply for the densely populated Indo-Gangetic plains. Despite the importance of water resources in light of climate change, the relative contributions of rainfall, snow and glacier melt to discharge are not well understood, due to the scarcity of ground-based data in this complex terrain. Here, we quantify discharge sources in the Sutlej Valley, western Himalaya, from 2000 to 2012 with a distributed hydrological model that is based on daily, ground-calibrated remote-sensing observation. Based on the consistently good model performance, we analyzed the spatiotemporal distribution of hydrologic components and quantified their contribution to river discharge. Our results indicate that the Sutlej River's annual discharge at the mountain front is sourced to 55% by effective rainfall (rainfall reduced by evapotranspiration), 35% by snow melt and 10% by glacier melt. In the high-elevation orogenic interior glacial runoff contributes ∼30% to annual river discharge. These glacier melt contributions are especially important during years with substantially reduced rainfall and snowmelt runoff, as during 2004, to compensate for low river discharge and ensure sustained water supply and hydropower generation. In 2004, discharge of the Sutlej River totaled only half the maximum annual discharge; with 17.3% being sourced by glacier melt. Our findings underscore the importance of calibrating remote-sensing data with ground-based data to constrain hydrological models with reasonable accuracy. For instance, we found that TRMM (Tropical Rainfall Measuring Mission) product 3B42 V7 systematically overestimates rainfall in arid regions of our study area by a factor of up to 5. By quantifying the spatiotemporal distribution of water resources we provide an important assessment of the potential impact of global warming on river discharge in the western Himalaya. Given the near-global coverage of the utilized remote-sensing datasets this hydrological modeling approach can be readily transferred to other data-sparse regions.  相似文献   

11.
This paper analyses the skills of fuzzy computing based rainfall–runoff model in real time flood forecasting. The potential of fuzzy computing has been demonstrated by developing a model for forecasting the river flow of Narmada basin in India. This work has demonstrated that fuzzy models can take advantage of their capability to simulate the unknown relationships between a set of relevant hydrological data such as rainfall and river flow. Many combinations of input variables were presented to the model with varying structures as a sensitivity study to verify the conclusions about the coherence between precipitation, upstream runoff and total watershed runoff. The most appropriate set of input variables was determined, and the study suggests that the river flow of Narmada behaves more like an autoregressive process. As the precipitation is weighted only a little by the model, the last time‐steps of measured runoff are dominating the forecast. Thus a forecast based on expected rainfall becomes very inaccurate. Although good results for one‐step‐ahead forecasts are received, the accuracy deteriorates as the lead time increases. Using the one‐step‐ahead forecast model recursively to predict flows at higher lead time, however, produces better results as opposed to different independent fuzzy models to forecast flows at various lead times. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

12.
《水文科学杂志》2013,58(4):567-584
Abstract

Reliable, real-time river flow forecasting in Africa on a time scale of days can provide enormous humanitarian and economic benefits. This study investigates the feasibility of using daily rainfall estimates based on cold cloud duration (CCD) derived from Meteosat thermal infrared imagery as input to a simple rainfall—runoff model and also whether such estimates can be improved by the inclusion of information from numerical weather prediction (NWP) analysis models. The Bakoye catchment in Mali, West Africa has been used as a test area. The data available for the study covered the main months of the rainy season for three years. The rainfall estimates were initially validated against gauge data. Improvements in quality were observed when information relating to African Easterly Wave phase and storm type was included in a multiple linear regression (MR) algorithm. The estimates were also tested by using them as input to a rainfall—runoff model. When contemporaneous calibrations from raingauges were available for calibration, both CCD-only and MR rainfall estimates gave more accurate river flow forecasts than when using raingauge data alone. In the absence of contemporaneous calibrations, the performance was reduced but the MR did better than the CCDonly input in all years. The use of satellite-derived vegetation index did not improve the quality of the river flow forecasts.  相似文献   

13.
Abstract

One of the main challenges faced by hydrologists and water engineers is the estimation of variables needed for water resources planning and management in ungauged river basins. To this end, techniques for transposing information, such as hydrological regional analyses, are widely employed. A method is presented for regionalizing flow-duration curves (FDCs) in perennial, intermittent and ephemeral rivers, based on the extended Burr XII probability distribution. This distribution shows great flexibility to fit data, with accurate reproduction of flow extremes. The performance analysis showed that, in general, the regional models are able to synthesize FDCs in ungauged basins, with a few possible drawbacks in the application of the method to intermittent and ephemeral rivers. In addition to the regional models, we summarize the experience of using synthetic FDCs for the indirect calibration of the Rio Grande rainfall–runoff model parameters in ungauged basins.

Editor D. Koutsoyiannis

Citation Costa, V., Fernandes, W., and Naghettini, M., 2013. Regional models of flow-duration curves of perennial and intermittent streams and their use for calibrating the parameters of a rainfall–runoff model. Hydrological Sciences Journal, 59 (2), 262–277.  相似文献   

14.
C. Fleurant  B. Kartiwa  B. Roland 《水文研究》2006,20(18):3879-3895
The rainfall‐runoff modelling of a river basin can be divided into two processes: the production function and the transfer function. The production function determines the proportion of gross rainfall actually involved in the runoff. The transfer function spreads the net rainfall over time and space in the river basin. Such a transfer function can be modelled using the approach of the geomorphological instantaneous unit hydrograph (GIUH). The effectiveness of geomorphological models is actually revealed in rainfall‐runoff modelling, where hydrologic data are desperately lacking, just as in ungauged basins. These models make it possible to forecast the hydrograph shape and runoff variation versus time at the basin outlet. This article is an introduction to a new GIUH model that proves to be simple and analytical. Its geomorphological parameters are easily available on a map or from a digital elevation model. This model is based on general hypotheses on symmetry that provide it with multiscale versatile characteristics. After having validated the model in river basins of very different nature and size, we present an application of this model for rainfall‐runoff modelling. Since parameters are determined relying on real geomorphological data, no calibration is necessary, and it is then possible to carry out rainfall‐runoff simulations in ungauged river basins. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

15.
ABSTRACT

Streamflow data are important for river management and the calibration of hydrological models. However, such data are only available for gauged catchments. Citizen science offers an alternative data source, and can be used to estimate streamflow at ungauged sites. We evaluated the accuracy of crowdsourced streamflow estimates for 10 streams in Switzerland by asking citizens to estimate streamflow either directly, or based on the estimated width, depth and velocity of the stream. Additionally, we asked them to estimate the stream level class by comparing the current stream level with a picture that included a virtual staff gauge. To compare the different estimates, the stream level class estimates were converted into streamflow. The results indicate that stream level classes were estimated more accurately than streamflow, and more accurately represented high and low flow conditions. Based on this result, we suggest that citizen science projects focus on stream level class estimates instead of streamflow estimates.  相似文献   

16.
Ragab Ragab  John Bromley 《水文研究》2010,24(19):2663-2680
A newly Integrated Hydrological Modelling System (IHMS) has been developed to study the impact of changes in climate, land use and water management on groundwater and seawater intrusion (SWI) into coastal areas. The system represents the combination of three models, which can, if required, be run separately. It has been designed to assess the combined impact of climate, land use and groundwater abstraction changes on river, drainage and groundwater flows, groundwater levels and, where appropriate, SWI. The approach is interdisciplinary and reflects an integrated water management approach. The system comprises three packages: the Distributed Catchment Scale Model (DiCaSM), MODFLOW (96 and 2000) and SWI models. In addition to estimating all water balance components, DiCaSM, produces the recharge data that are used as input to the groundwater flow model of the US Geological Survey, MODFLOW. The latter subsequently generates the head distribution and groundwater flows that are used as input to the SWI model, SWI. Thus, any changes in land use, rainfall, water management, abstraction, etc. at the surface are first handled by DiCaSM, then by MODFLOW and finally by the SWI. The three models operate at different spatial and temporal scales and a facility (interface utilities between models) to aggregate/disaggregate input/output data to meet a desired spatial and temporal scale was developed allowing smooth and easy communication between the three models. As MODFLOW and SWI are published and in the public domain, this article focuses on DiCaSM, the newly developed unsaturated zone DiCaSM and equally important the interfacing utilities between the three models. DiCaSM simulates a number of hydrological processes: rainfall interception, evapotranspiration, surface runoff, infiltration, soil water movement in the root zone, plant water uptake, crop growth, stream flow and groundwater recharge. Input requirements include distributed data sets of rainfall, land use, soil types and digital terrain; climate data input can be either distributed or non‐distributed. The model produces distributed and time series output of all water balance components including potential evapotranspiration, actual evapotranspiration, rainfall interception, infiltration, plant water uptake, transpiration, soil water content, soil moisture (SM) deficit, groundwater recharge rate, stream flow and surface runoff. This article focuses on details of the hydrological processes and the various equations used in DiCaSM, as well as the nature of the interface to the MODFLOW and SWI models. Furthermore, the results of preliminary tests of DiCaSM are reported; these include tests related to the ability of the model to predict the SM content of surface and subsurface soil layers, as well as groundwater levels. The latter demonstrates how the groundwater recharge calculated from DiCaSM can be used as input into the groundwater model MODFLOW using aggregation and disaggregation algorithms (built into the interface utility). SWI has also been run successfully with hypothetical examples and was able to reproduce the results of some of the original examples of Bakker and Schaars ( 2005 ). In the subsequent articles, the results of applications to different catchments will be reported. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

17.
This paper reports the results of an investigation into flood simulation by areal rainfall estimated from the combination of gauged and radar rainfalls and a rainfall–runoff model on the Anseong‐cheon basin in the southern part of Korea. The spatial and temporal characteristics and behaviour of rainfall are analysed using various approaches combining radar and rain gauges: (1) using kriging of the rain gauge alone; (2) using radar data alone; (3) using mean field bias (MFB) of both radar and rain gauges; and (4) using conditional merging technique (CM) of both radar and rain gauges. To evaluate these methods, statistics and hyetograph for rain gauges and radar rainfalls were compared using hourly radar rainfall data from the Imjin‐river, Gangwha, rainfall radar site, Korea. Then, in order to evaluate the performance of flood estimates using different rainfall estimation methods, rainfall–runoff simulation was conducted using the physics‐based distributed hydrologic model, Vflo?. The flood runoff hydrograph was used to compare the calculated hydrographs with the observed one. Results show that the rainfall field estimated by CM methods improved flood estimates, because it optimally combines rainfall fields representing actual spatial and temporal characteristics of rainfall. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

18.
Abstract

Modelling of the rainfall–runoff transformation process and routing of river flows in the Kilombero River basin and its five sub-catchments within the Rufiji River basin in Tanzania was undertaken using three system (black-box) models—a simple linear model, a linear perturbation model and a linear varying gain factor model—in their linear transfer function forms. A lumped conceptual model—the soil moisture accounting and routing model—was also applied to the sub-catchments and the basin. The HEC-HMS model, which is a distributed model, was applied only to the entire Kilombero River basin. River discharge, rainfall and potential evaporation data were used as inputs to the appropriate models and it was observed that sometimes the system models performed better than complex hydrological models, especially in large catchments, illustrating the usefulness of using simple black-box models in datascarce situations.  相似文献   

19.
This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

20.
We explore the impact of uncertainties in the spatial–temporal distribution of rainfall on the prediction of peak discharge in a typical mountain basin. To this end, we use a stochastic generator previously developed for rainfall downscaling, and we estimate the basin response by adopting a semi-distributed hydrological model. The results of the analysis provide information on the minimum rainfall resolution needed for operational flood forecasting, and confirm the sensitivity of peak discharge estimates to errors in the determination of the power spectrum of the precipitation field.  相似文献   

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

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