首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
High regulation of dams or sluices disturbs flow regimes and pollutant transformation process significantly in most basins over the world. The water quality and quantity simulation in highly regulated river basins is always a very complicated task. The Huai River Basin is a typical area in China with the highest density of water projects and serious pollution problems simultaneously. In this paper, a procedure based on Soil and Water Assessment Tool (SWAT) was proposed to carry out the water quantity and quality simulation of Huai River Basin by incorporating the operation rules of dams or sluices into the reservoir regulation module. The water quality module in SWAT is extended to suit the actual situation of river basin in China. The results show that: for runoff simulation, 28 of all the 38 stations (73.7%) have the acceptable performance, with the average correlation coefficient and efficiency coefficient of 0.83 and 0.66, respectively in the calibration period. In the validation period, 17 of all the 34 stations (50%) have the acceptable performance, with the average correlation coefficient and efficiency coefficient of 0.77 and 0.54, respectively. The model performs the worst for reservoirs, little better for sluices and the best for unregulated stations. Comparing to the low flow and high flow simulation of original SWAT model, the low flow simulation of reservoirs and the high flow simulation of sluices are much better because of their different main purposes. For water quality simulation, the standard reaching rates of NH3–N and CODMn are 55.9% (19/34) and 67.6% (23/34) respectively for all the stations, and the average relation coefficients are 0.46 and 0.48 respectively. Comparing with the results of original SWAT model, the improved model better reproduces the long-term water quantity and quality processes in the Huai River Basin of China. This study provides a new approach and reference to understand the variation of water quantity and quality in highly regulated river basin, and is expected as technical support for the environment restoration and integrated management in the basins, especially in China.  相似文献   

2.
The soil and water assessment tool (SWAT) has been widely used and thoroughly tested in many places in the world. The application of the SWAT model has pointed out that 2 of the major weaknesses of SWAT are related to the nonspatial reference of the hydrologic response unit concept and to the simplified groundwater concept, which contribute to its low performance in baseflow simulation and its inability to simulate regional groundwater flow. This study modified the groundwater module of SWAT to overcome the above limitations. The modified groundwater module has 2 aquifers. The local aquifer, which is the shallow aquifer in the original SWAT, represents a local groundwater flow system. The regional aquifer, which replaces the deep aquifer of the original SWAT, represents intermediate and regional groundwater flow systems. Groundwater recharge is partitioned into local and regional aquifer recharges. The regional aquifer is represented by a multicell aquifer (MCA) model. The regional aquifer is discretized into cells using the Thiessen polygon method, where centres of the cells are locations of groundwater observation wells. Groundwater flow between cells is modelled using Darcy's law. Return flow from cell to stream is conceptualized using a non‐linear storage–discharge relationship. The SWAT model with the modified aquifer module, the so‐called SWAT‐MCA, was tested in 2 basins (Wipperau and Neetze) with porous aquifers in a lowland area in Lower Saxony, Germany. Results from the Wipperau basin show that the SWAT‐MCA model is able (a) to simulate baseflow in a lowland area (where baseflow is a dominant source of streamflow) better than the original model and (b) to simulate regional groundwater flow, shown by the simulated groundwater levels in cells, quite well.  相似文献   

3.
Using a mass balance algorithm, this study develops an extension module that can be embedded in the commonly used Soil and Water Assessment Tool (SWAT). This module makes it possible to assess effects of riparian wetlands on runoff and sediment yields at a watershed scale, which is very important for aquatic ecosystem management but rarely documented in the literature. In addition to delineating boundaries of a watershed and its subwatersheds, the module groups riparian wetlands within a subwatershed into an equivalent wetland for modelling purposes. Further, the module has functions to compute upland drainage area and other parameters (e.g. maximum volume) for the equivalent wetland based on digital elevation model, stream network, land use, soil and wetland distribution GIS datasets. SWAT is used to estimate and route runoff and sediment generated from upland drainage area. The lateral exchange processes between riparian wetlands and their hydraulically connected streams are simulated by the extension module. The developed module is empirically applied to the 53 km2 Upper Canagagigue Creek watershed located in Southern Ontario of Canada. The simulation results indicate that the module can make SWAT more reasonably predict flow and sediment loads at the outlet of the watershed and better represent the hydrologic processes within it. The simulation is sensitive to errors of wetland parameters and channel geometry. The approach of embedding the module into SWAT enables simulation of hydrologic processes in riparian wetlands, evaluation of wetland effects on regulating stream flow and sediment loading and assessment of various wetland restoration scenarios. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

4.
Stream flow predictions in ungauged basins are one of the most challenging tasks in surface water hydrology because of nonavailability of data and system heterogeneity. This study proposes a method to quantify stream flow predictive uncertainty of distributed hydrologic models for ungauged basins. The method is based on the concepts of deriving probability distribution of model's sensitive parameters by using measured data from a gauged basin and transferring the distribution to hydrologically similar ungauged basins for stream flow predictions. A Monte Carlo simulation of the hydrologic model using sampled parameter sets with assumed probability distribution is conducted. The posterior probability distributions of the sensitive parameters are then computed using a Bayesian approach. In addition, preselected threshold values of likelihood measure of simulations are employed for sizing the parameter range, which helps reduce the predictive uncertainty. The proposed method is illustrated through two case studies using two hydrologically independent sub‐basins in the Cedar Creek watershed located in Texas, USA, using the Soil and Water Assessment Tool (SWAT) model. The probability distribution of the SWAT parameters is derived from the data from one of the sub‐basins and is applied for simulation in the other sub‐basin considered as pseudo‐ungauged. In order to assess the robustness of the method, the numerical exercise is repeated by reversing the gauged and pseudo‐ungauged basins. The results are subsequently compared with the measured stream flow from the sub‐basins. It is observed that the measured stream flow in the pseudo‐ungauged basin lies well within the estimated confidence band of predicted stream flow. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

5.
针对现有的河道水流洪水演算模型只能模拟单一变量(流量或水位)的问题,以水流连续方程和河段蓄水量的两种不同表达形式(蓄水量等于平均过水断面面积与河段长乘积,蓄水量等于河段平均流量与传播时间的乘积)为基础,对马斯京根模型进行了通用性改进,提出了双变量耦合通用演算模型.选取了四大水系(包括内陆河流和入海河流)的16个河段汛期洪水资料进行模型检验,模型验证考虑了地理范围、不同的河段特征和水力特征、洪水量级等因素,全面地检验了模型结构的合理性和模拟实际洪水的有效性.将双变量耦合通用演算模型与传统的马斯京根法进行了效果比较,结果表明双变量耦合通用演算模型的模拟精度高于马斯京根法,模拟效果比马斯京根法稳定一些,而且具有较好的通用性.  相似文献   

6.
Streamflow measurements provide information about the flow generation characteristics of land surfaces as well as the flow transferring nature of the channel network. In this study, such flow transferring properties of the channel network that were obtained from downstream flow observations were used for predicting flow in ungauged basins. A temporally averaged transfer function (ATF) of the channel segments of Kentucky River Basin (KRB) in Kentucky, USA, was extracted from observed hydrographs in a time‐invariant system as a function of drainage area. The ATF was regionalized through multiple regression analysis for 194 combinations of drainage areas that differ in topography, terrain, and geology. The application of ATF for flow prediction in ungauged basins was performed for Goose Creek, a subbasin of KRB by integrating ATF with the TOPMODEL. In addition, the ATF was shown to be capable of providing calibration and validation data for ungauged basins in a backward technique from a measured stream gauge downstream, with minimal data requirement of drainage area. The applicability of ATF was illustrated across a range of streamflow conditions from watersheds that varied greatly in their terrain and geology. Nash–Sutcliffe efficiency of the proposed method, as a function of drainage areas of the corresponding basins, to predict daily streamflow from ungauged basins ranged from 0.83 to 0.92. The results of the study concluded that the ATF obtained from measured streamflow thus proved to be a quick and simple tool for assessment of streamflow in both operational and modeling hydrology. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

7.
Integrated river basin models should provide a spatially distributed representation of basin hydrology and transport processes to allow for spatially implementing specific management and conservation measures. To accomplish this, the Soil and Water Assessment Tool (SWAT) was modified by integrating a landscape routing model to simulate water flow across discretized routing units. This paper presents a grid‐based version of the SWAT landscape model that has been developed to enhance the spatial representation of hydrology and transport processes. The modified model uses a new flow separation index that considers topographic features and soil properties to capture channel and landscape flow processes related to specific landscape positions. The resulting model is spatially fully distributed and includes surface, lateral and groundwater fluxes in each grid cell of the watershed. Furthermore, it more closely represents the spatially heterogeneous distributed flow and transport processes in a watershed. The model was calibrated and validated for the Little River Watershed (LRW) near Tifton, Georgia (USA). Water balance simulations as well as the spatial distribution of surface runoff, subsurface flow and evapotranspiration are examined. Model results indicate that groundwater flow is the dominant landscape process in the LRW. Results are promising, and satisfactory output was obtained with the presented grid‐based SWAT landscape model. Nash–Sutcliffe model efficiencies for daily stream flow were 0.59 and 0.63 for calibration and validation periods, and the model reasonably simulates the impact of the landscape position on surface runoff, subsurface flow and evapotranspiration. Additional revision of the model will likely be necessary to adequately represent temporal variations of transport and flow processes in a watershed. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

8.
This paper presents the results of a study that compares the sediment routing of the Simiyu River using the hydrologic model, Soil and Water Assessment Tool (SWAT) and the 1D hydrodynamic simulation software for Rivers and Estuaries (SOBEK-RE) model. Routing in SWAT is completed using the simplified Bagnold’s equation and in the SOBEK-RE model is undertaken using the Saint Venant equation. The upstream boundary conditions for the routing modules were derived from the subcatchments sediment yields that were estimated by SWAT using the Modified Universal Soil Loss Equation (MUSLE). The sediment loads extrapolated or interpolated from the sediment rating curve for the catchment outlet were used for calibration and validation purposes. The SWAT model predicted an erosion rate of 2.09 Mt/yr. The total sediment load transported to the main outlet of the catchment simulated by the SWAT and SOBEK-RE models was equal to 2.94 and 2.72 Mt/yr, respectively. Thus the models computed a net erosion in the channels of 0.84 Mt/yr (SWAT) and 0.63 Mt/yr (SOBEK-RE). When comparing the results of the models for the different reaches of the main channel and main tributaries, the models showed different results both in magnitude and in sign (erosion/deposition). However, in a situation where data is scarce (such as grain size, channel geometry), the more complex hydrodynamic model does not necessarily lead to more reliable results.  相似文献   

9.
Bank erosion is the main source of suspended sediment (SS) and diffuse total phosphorus (TP) in many lowland catchments. This study compared a physically based sediment routing method (Physical method), which distinguishes between stream bed and bank erosion, with the original sediment routing method (Original method) within the Soil and Water Assessment Tool (SWAT) version 2009, for simulating SS and TP losses from a lowland catchment. A SWAT model was set up for the lowland River Odense catchment in Denmark and calibrated against observed stream flow and phosphate (PO4) loads. On the basis of an initial calibration of hydrological and PO4 parameters, the SWAT model with the Original method (Original model) and the SWAT model with the Physical method (Physical model) were calibrated separately against observed SS and TP loads. The SWAT model simulated daily stream flow well but underestimated PO4 loads. The Physical model simulated daily SS and TP better than the Original model. The simulated contribution of bank erosion to SS in the Physical model (99%) was close to the estimated contribution from in situ erosion measurements (90–94%). Compared with the Original method, the Physical method is not only more conceptually correct but also improves model performance. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

10.
Despite the significant role of precipitation in the hydrological cycle, few studies have been conducted to evaluate the impacts of the temporal resolution of rainfall inputs on the performance of SWAT (soil and water assessment tool) models in large-sized river basins. In this study, both daily and hourly rainfall observations at 28 rainfall stations were used as inputs to SWAT for daily streamflow simulation in the Upper Huai River Basin. Study results have demonstrated that the SWAT model with hourly rainfall inputs performed better than the model with daily rainfall inputs in daily streamflow simulation, primarily due to its better capability of simulating peak flows during the flood season. The sub-daily SWAT model estimated that 58 % of streamflow was contributed by baseflow compared to 34 % estimated by the daily model. Using the future daily and 3-h precipitation projections under the RCP (Representative Concentration Pathways) 4.5 scenario as inputs, the sub-daily SWAT model predicted a larger amount of monthly maximum daily flow during the wet years than the daily model. The differences between the daily and sub-daily SWAT model simulation results indicated that temporal rainfall resolution could have much impact on the simulation of hydrological process, streamflow, and consequently pollutant transport by SWAT models. There is an imperative need for more studies to examine the effects of temporal rainfall resolution on the simulation of hydrological and water pollutant transport processes by SWAT in river basins of different environmental conditions.  相似文献   

11.
River corridors exhibit landforms nested within landforms repeatedly down spatial scales. In Pasternack et al. ( 2018 ), a new, scale‐independent, hierarchical river classification was developed that uses five landform types to map the domains of a single fluvial process – flow convergence routing – at each of three–five spatial scales. Given those methods, this study investigated the details of how flow convergence routing organizes nested landform sequences. The method involved analyzing landform abundance, sequencing, and hierarchical nesting along the 35 km gravel/cobble lower Yuba River in California. Independent testing of flow convergence routing found that hydraulic patterns at every flow matched the essential predictions from classification, substantiating the process–morphology link. River width and bed elevation sequences exhibit large, nonrandom, and linked oscillations structured to preferentially yield wide bars and constricted pools at base flow and bankfull flow. At a flow of 8.44 times bankfull, there is still an abundance of wide bar and constricted pool landforms, but larger topographic drivers also yield an abundance of nozzle and oversized landforms. The nested structure of flow convergence routing landforms reveals that base flow and bankfull landforms are nested together within specific floodprone valley landform types, and these landform types control channel morphodynamics during moderate to large floods. As a result, this study calls into question the prevailing theory that the bankfull channel of a gravel/cobble river is controlled by in‐channel, bankfull, and/or small flood flows. Such flows may initiate sediment transport, but they are too small to control landform organization in a gravel/cobble river with topographic complexity. Copyright © 2018 John Wiley & Sons, Ltd.  相似文献   

12.
Lei Wang  Jaehyung Yu 《水文研究》2012,26(19):2973-2984
The construction of stormwater detention basins is a best management practice to effectively control floods, to provide additional surface storage for excess floodwater and to compensate for the adverse effects of urban development. Traditional field‐based levelling survey methods are very time consuming and subject to human‐induced arbitrariness and error. This article presents an approach to modelling detention basins measured from light detection and ranging remote sensing data. A case study is illustrated by using the White Oak Bayou watershed of Harris County, Texas. The storage–stage curve obtained from the volumetric analysis is used in a modified detention basins routing model, which was developed by adding the weir structure control to the traditional hydrologic reservoir routing equations. The model simulation showed that the peak flow of the synthetic 100‐year reoccurrence event was effectively reduced and delayed by the detention basins. The comparison with the simulation results from the traditional reservoir routing model suggested that previous studies using the reservoir routing model were likely to underestimate the flood reduction effect of detention basins. The sensitivity analysis of the parameters showed that the detention basin design and evaluation should pay more attention on the weir height and river channel's roughness. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

13.
Hydrological models are useful tools to analyze present and future conditions of water quantity and quality. The integrated modelling of water and nutrients needs an adequate representation of the different discharge components. In common with many lowlands, groundwater contribution to the discharge in the North German lowlands is a key factor for a reasonable representation of the water balance, especially in low flow periods. Several studies revealed that the widely used Soil and Water Assessment Tool (SWAT) model performs poorly for low flow periods. This paper deals with the extension of the groundwater module of the SWAT model to enhance low flow representation. The current two‐storage concept of SWAT was further developed to a three‐storage concept. This was realized due to modification of the groundwater module by splitting the active groundwater storage into a fast and a slow contributing aquifer. The results of this study show that the groundwater module with three storages leads to a good prediction of the overall discharge especially for the recession limbs and the low flow periods. The improved performance is reflected in the signature measures for the mid‐segment (percent bias ?2.4% vs ?15.9%) and the low segment (percent bias 14.8% vs 46.8%) of the flow duration curve. The three‐storage groundwater module is more process oriented than the original version due to the introduction of a fast and a slow groundwater flow component. The three‐storage version includes a modular approach, because groundwater storages can be activated or deactivated independently for subbasin and hydrological response unit level. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
Abstract

In a typical reservoir routing problem, the givens are the inflow hydrograph and reservoir characteristic functions. Flood attenuation investigations can be easily accomplished using a hydrological or hydraulic routing of the inflow hydrograph to obtain the reservoir outflow hydrograph, unless the inflow hydrograph is unavailable. Although attempts for runoff simulation have been made in ungauged basins, there is only a limited degree of success in special cases. Those approaches are, in general, not suitable for basins with a reservoir. The objective of this study is to propose a procedure for flood attenuation estimation in ungauged reservoir basins. In this study, a kinematic-wave based geomorphic IUH model was adopted. The reservoir inflow hydrograph was generated through convolution integration using the rainfall excess and basin geomorphic information. Consequently, a fourth-order Runge-Kutta method was used to route the inflow hydrograph to obtain the reservoir outflow hydrograph without the aid of recorded flow data. Flood attenuation was estimated through the analysis of the inflow and outflow hydrographs of the reservoir. An ungauged reservoir basin in southern Taiwan is presented as an example to show the applicability of the proposed analytical procedure. The analytical results provide valuable information for downstream flood control work for different return periods.  相似文献   

15.
Drainage area and the variation of channel geometry downstream   总被引：1，自引：0，他引：1  
Recent geomorphological studies tend to deal with small basins. The understanding of small basin dynamics provides important information for the understanding of large basin dynamics assuming that the extrapolation of small basin data to larger basins is valid. This work tests the validity of this extrapolation of data with reference to channel geometry. An analysis of the variation of channel width downstream reveals that the value b =0.5 (W = aQb) is a ‘good’ average. However, the use of a one-line model consisting of a simple power function incurs a loss of a considerable amount of relevant information concerning the channel form and hence the channel processes. It has been shown that the –b– value for small basins and very big basins is lower than the one for the intermediate basins.  相似文献   

16.
This paper investigates the specific contributions of river network geomorphology, hillslope flow dynamics and channel routing to the scaling behavior of the hydrologic response as function of drainage area. Scaling relationships emerged from the observations of geomorphological and hydrological data and were reproduced in previous works through mathematical models, for both idealized self-similar networks and natural basins. Recent literature highlighted that scale invariance of hydrological quantities depends not only on the metrics of the drainage catchment but also on effective flow routing. In this study we employ a geomorphological width function scheme to test the simple scaling hypothesis adopting more realistic dynamic conditions than in previous approaches, specifically taking into account the role of hillslopes. The analysis is based on the derivation of the characteristic distributions of path lengths and travel times, inferred from DEM processing and measurements of rainfall and runoff data. The study area is located in the Tiber River region (central Italy).Results indicate that, while scaling properties clearly emerge when the hydrologic response is defined on the basis of the sole geomorphology, scale invariance is broken when less idealized flow dynamics are taken into account. Lack of scaling appears in particular as a consequence of the catchment to catchment variability of hillslope velocities.  相似文献   

17.
The peakedness of a basin and its variation with drainage area were analyzed for three areas. Peakedness of a basin is calculated as mean flow as a percentage of highest flow. A fitted power curve relating ‘Peakedness index’ (PI) to drainage area for each of the three areas indicates a break point in a basin of about 300 km2. This break point divides the basins into small basins which are more peaky and large basins which are less peaky. The break point is an outcome of a difference in order of magnitude between channel flow velocity from the headwater sources and hillside flow velocity. When the basin responds to hillslope flow the runoff from the head water sources has already flown about 30 km downstream.  相似文献   

18.
Hydrological interaction between surface and subsurface water systems has a significant impact on water quality, ecosystems and biogeochemistry cycling of both systems. Distributed models have been developed to simulate this function, but they require detailed spatial inputs and extensive computation time. The soil and water assessment tool (SWAT) model is a semi‐distributed model that has been successfully applied around the world. However, it has not been able to simulate the two‐way exchanges between surface water and groundwater. In this study, the SWAT‐landscape unit (LU) model – based on a catena method that routes flow across three LUs (the divide, the hillslope and the valley) – was modified and applied in the floodplain of the Garonne River. The modified model was called SWAT‐LUD. Darcy's equation was applied to simulate groundwater flow. The algorithm for surface water‐level simulation during flooding periods was modified, and the influence of flooding on groundwater levels was added to the model. Chloride was chosen as a conservative tracer to test simulated water exchanges. The simulated water exchange quantity from SWAT‐LUD was compared with the output of a two‐dimensional distributed model, surface–subsurface water exchange model. The results showed that simulated groundwater levels in the LU adjoining the river matched the observed data very well. Additionally, SWAT‐LUD model was able to reflect the actual water exchange between the river and the aquifer. It showed that river water discharge has a significant influence on the surface–groundwater exchanges. The main water flow direction in the river/groundwater interface was from groundwater to river; water that flowed in this direction accounted for 65% of the total exchanged water volume. The water mixing occurs mainly during high hydraulic periods. Flooded water was important for the surface–subsurface water exchange process; it accounted for 69% of total water that flowed from the river to the aquifer. The new module also provides the option of simulating pollution transfer occurring at the river/groundwater interface at the catchment scale. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

19.
The delicate balance between human utilization and sustaining its pristine biodiversity in the Mara River basin (MRB) is being threatened because of the expansion of agriculture, deforestation, human settlement, erosion and sedimentation and extreme flow events. This study assessed the applicability of the Soil and Water Assessment Tool (SWAT) model for long‐term rainfall–runoff simulation in MRB. The possibilities of combining/extending gage rainfall data with satellite rainfall estimates were investigated. Monthly satellite rainfall estimates not only overestimated but also lacked the variability of observed rainfall to substitute gage rainfall in model simulation. Uncertainties related to the quality and availability of input data were addressed. Sensitivity and uncertainty analysis was reported for alternative model components and hydrologic parameters used in SWAT. Mean sensitivity indices of SWAT parameters in MRB varied with and without observed discharge data. The manual assessment of individual parameters indicated heterogeneous response among sub‐basins of MRB. SWAT was calibrated and validated with 10 years of discharge data at Bomet (Nyangores River), Mulot (Amala River) and Mara Mines (Mara River) stations. Model performance varied from satisfactory at Mara Mines to fair at Bomet and weak at Mulot. The (Nash–Sutcliff efficiency, coefficient of determination) results of calibration and validation at Mara Mines were (0.68, 0.69) and (0.43, 0.44), respectively. Two years of moving time window and flow frequency analysis showed that SWAT performance in MRB heavily relied on quality and abundance of discharge data. Given the 5.5% area contribution of Amala sub‐basin as well as uncertainty and scarcity of input data, SWAT has the potential to simulate the rainfall runoff process in the MRB. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

20.
The objective of this study is to incorporate a time‐dependent Soil Moisture Accounting (SMA) based Curve Number method (SMA_CN) in Soil and Water Assessment Tool (SWAT) and compare its performance with the existing CN method in SWAT by simulating the hydrology of two agricultural watersheds in Indiana, USA. Results show that fusion of the SMA_CN method causes decrease in runoff volume and increase in profile soil moisture content, associated with larger groundwater contribution to the streamflow. In addition, the higher amount of moisture in the soil profile slightly elevates the actual evapotranspiration. The SMA‐based SWAT configuration consistently produces improved goodness‐of‐fit scores and less uncertain outputs with respect to streamflow during both calibration and validation. The SMA_CN method exhibits a better match with the observed data for all flow regimes, thereby addressing issues related to peak and low flow predictions by SWAT in many past studies. Comparison of the calibrated model outputs with field‐scale soil moisture observations reveals that the SMA overhauling enables SWAT to represent soil moisture condition more accurately, with better response to the incident rainfall dynamics. While the results from the modification of the CN method in SWAT are promising, more studies including watersheds with various physical and climatic settings are needed to validate the proposed approach. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

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

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