首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
《国际泥沙研究》2020,35(2):115-124
Bed-load transport plays a critical role in river morphological change and has an important impact on river ecology.Although there is good understanding of the role of the variation of river bed grain size on transport dynamics in equilibrium conditions,much less is understood for non-equilibrium conditions when the channel is either aggrading or degrading.In particular,the relative role of different grain sizes in the promotion and hindering of the transport of coarse and fine fractions in a degrading channel has yet to be investigated.The current study attempts to provide new understanding through a series of flume experiments done using uniform and graded sediment particles.The experiments revealed coarser grain-size fractions for a poorly-sorted sediment,relative to uniform-sized sediment,reduced the transport of finer grains and finer fractions enhanced the transport of coarse grains.This hinderingpromotion effect,caused by relative hiding and exposure of finer and coarse fractions,increased with bed slope and decreased with relative submergence.In particular,as relative submergence increased,the graded fractions tended towards behaving more like their unifo rm-sized counterparts.Also,the bed-load parameter of the graded fractions increased more with a rise in bed slope than observed for the uniformsized counterparts.These results revealed,for degrading channel conditions,such as downstream of a dam,bed-load equations developed for uniform bed sediment are inappropriate for use in natural river systems,particularly in mountain streams.Furthermore,changes in river bed composition due to activities that enhance the input of hill-slope sediment,such as fire,logging,and agricultural development,are likely to cause significant changes in river morphology.  相似文献   

2.
LIN~DUCTIONTheinvenigationsofrivermorphologyandOfriverbeddebilityresultalwaysinidalsabouttheinterrelationshipbetweenthehydrauliccondihons,sc(hant~rtandbeddynndcs.ThedeCOuragingfactforsuchinveedgationsarethecomplexityoftheseSystemsandthehugenumberofimportantvariables.IfwefollowSimonsetal(1965)andall~eneeandtheconcentrationsofSUSPendedlOadandbedloadthefollowingparametersarethepnncapleVarialllesf-av~velocityofnow,diStributionofnowvelociticsacdeadofturbulence-meandepthacthedistributionofde…  相似文献   

3.
A comprehensive understanding of the dynamics of erosion and sedimentation in reservoirs under different management conditions is required to anticipate sedimentation issues and implement effective sediment management strategies. This paper describes a unique approach combining fluvial geomorphology tools and morphodynamic modeling for analyzing the sediment dynamics of an elongated hydropower reservoir subjected to management operations: the Génissiat Reservoir on the Rhône River. Functional sub‐reaches representative of the reservoir morphodynamics were delineated by adapting natural river segmentation methods to elongated reservoirs. The segmentation revealed the link between the spatial and temporal reservoir changes and the variability of longitudinal flow conditions during reservoir management operations. An innovative modeling strategy, incorporating the reservoir segmentation into two sediment transport codes, was implemented to simulate the dynamics of erosion and sedimentation at the reach scale during historic events. One code used a bedload approach, based on the Exner equation with a transport capacity formula, and the other used a suspended load approach based on the advection–dispersion equation. This strategy provided a fair quantification of the dynamics of erosion and sedimentation at the reach scale during different management operations. This study showed that the reservoir morphodynamics is controlled by bedload transport in upper reaches, graded suspended load transport of sand in middle reaches and suspended load transport of fine sediments in lower reaches. Eventually, it allowed a better understanding of the impact of dam management on sediment dynamics. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

4.
Fluvial sediment transport in the high mountain Partnach River (Reintal Valley, Bavarian Alps) was investigated during a 10‐year observation period (2001–2010). During this period, the downstream reach of the Partnach River was decoupled from upstream sediment throughput by a rockslide deposit until 2005. In August 2005, the dam was partially breached during a flood event resulting in renewed sediment coupling between the upstream and downstream reaches. A comparison of pre‐ and post‐dambreak river sediment load data showed that the dissolved load dominated sediment transport prior to August 2005 with a switch to the dominance of bedload transport, post‐dambreak. The higher post‐dambreak bedload rates were particularly evident during the first years after the dam failure due to significant coarse material coupling between active sediment sources (undercut banks/talus cones) and the Partnach River. In the last years of the observation period (2009 and 2010) the dominance of dissolved load transport was re‐established. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

5.
Riverine sediments have played an important role in the morphological evolution of river channels and river deltas. However, the sediment regime in the many world's rivers has been altered in the context of global changes. In this study, temporal changes in the sediment regime of the Pearl River were examined at different time scales, that is, annual, seasonal, and monthly time scales, using the Mann–Kendall test. The results revealed that precipitation variability was responsible for monthly and seasonal distribution patterns of the sediment regime and the long‐term changes in the water discharge; however, dam operation has smoothed the seasonal distribution of water discharge and resulted in decreasing trends in the annual, wet‐season, and dry‐season sediment load series since the 1950s. Due to the different regulation magnitudes of dam operation, differences were observed in sediment regime changes among the three tributaries. In addition, human activities have altered the hysteresis of seasonal rating curves and affected hysteresis differences between increasing and decreasing water discharge stages. Sediment supply is an important factor controlling river channel dynamics, affecting channel morphology. From the 1950s to the 1980s, siltation was dominant in river channels across the West River and North River deltas in response to the sediment increases; however, scouring occurred in the East River deltas due to sediment reduction. Significant erosion occurred in river channels in the 1990s, which was mostly due to downcutting of the river bed caused by sand excavations and partly because of the reduced sediment load from upstream. Although sand excavations have been banned and controlled by authority agencies since 2000, the erosion of cross sections was still observed in the 2000s because of reduced sediment caused by dam construction. Our study examines the different effects of human activities on the sediment regime and downstream channel morphology, which is of substantial scientific importance for river management.  相似文献   

6.
We investigate the use of the short‐lived fallout radionuclide beryllium‐7 (7Be; t1/2 = 53·4 days) as a tracer of medium and coarse sand (0·25–2 mm), which transitions between transport in suspension and as bed load, and evaluate the effects of impoundment on seasonal and spatial variations in bed sedimentation. We measure 7Be activities in approximately monthly samples from point bar and streambed sediments in one unregulated and one regulated stream. In the regulated stream our sampling spanned an array of flow and management conditions during the annual transition from flood control in the winter and early spring to run‐of‐the‐river operation from late spring to autumn. Sediment stored behind the dam during the winter quickly became depleted in 7Be activity. This resulted in a pulse of ‘dead’ sediment released when the dam gates were opened in the spring which could be tracked as it moved downstream. Measured average sediment transport velocities (30–80 metres per day (m d?1)) exceed those typically reported for bulk bed load transport and are remarkably constant across varied flow regimes, possibly due to corresponding changes in bed sand fraction. Results also show that the length scale of the downstream impact of dam management on sediment transport is short (c. 1 km); beyond this distance the sediment trapped by the dam is replaced by new sediment from tributaries and other downstream sources. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

7.
A sediment mass balance constructed for a 16‐km reach of the Snake River downstream from Jackson Lake Dam (JLD) indicates that river regulation has reduced the magnitude of sediment mass balance deficit that would naturally exist in the absence of the dam. The sediment budget was constructed from calibrated bed load transport relations, which were used to model sediment flux into and through the study reach. Calibration of the transport relations was based on bed load transport data collected over a wide range of flows on the Snake River and its two major tributaries within the study area in 2006 and 2007. Comparison of actual flows with unregulated flows for the period since 1957 shows that operations of JLD have reduced annual peak flows and increased late summer flows. Painted tracer stones placed at five locations during the 2005 spring flood demonstrate that despite the reduction in flood magnitudes, common floods are capable of mobilizing the bed material. The sediment mass balance demonstrates that more sediment exits the study reach than is being supplied by tributaries. However, the volume of sediment exported using estimated unregulated hydrology indicates that the magnitude of the deficit would be greater in the absence of JLD. Calculations suggest that the Snake River was not in equilibrium before construction of JLD, but was naturally in sediment deficit. The conclusion that impoundment lessened a natural sediment deficit condition rather than causing sediment surplus could not have been predicted in the absence of sediment transport data, and highlights the value of transport data and calculation of sediment mass balance in informing dam operations. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

8.
Bo Wang  Yi-Jun Xu 《水文研究》2020,34(13):2864-2877
Bed material transport at river bifurcations is crucial for channel stability and downstream geomorphic dynamics. However, measurements of bed material transport at bifurcations of large alluvial rivers are difficult to make, and standard estimates based on the assumption of proportional partitioning of flow and bedload transport at bifurcations may be erroneous. In this study, we employed a combined approach based on observed topographic change (erosion/deposition) and bed material transport predicted from a one-dimensional model to investigate bed material fluxes near the engineering-controlled Mississippi-Atchafalaya River diversion, which is of great importance to sediment distribution and delivery to Louisiana's coast. Yang's (1973) sediment transport equation was utilized to estimate daily bed material loads upstream, downstream, and through the diversion during 2004–2013. Bathymetric changes in these channels were assessed with single beam data collected in 2004 and 2013. Results show that over the study period, 24% of the Mississippi River flow was diverted into the Atchafalaya River, while the rest remained in the mainstem Mississippi. Upstream of the diversion, the bed material yield was predicted to be 201 million metric tons (MT), of which approximately 35 MT (i.e., 17%) passed through the bifurcation channel to the Atchafalaya River. The findings from this study reveal that in the mainstem Mississippi, the percentage of bed material diversion (83%) is larger than the percentage of flow diversion (76%); Conversely, the diversion channel receives a disproportionate amount of flow (24%) relative to bed material supply (17%). Consequently, severe bed scouring occurred in the controlled Outflow Channel to the Atchafalaya River, while riverbed aggradation progressed in the mainstem Mississippi downstream of the diversion structures, implying reduced flow capacity and potential risk of a high backwater during megafloods. The study demonstrates that Yang's sediment transport equation provides plausible results of bed material fluxes for a highly complicated large river diversion, and that integration of the sediment transport equation with observed morphological changes in riverbed is a valuable approach to investigate sediment dynamics at controlled river bifurcations.  相似文献   

9.
The bed of estuaries is often characterized by ripples and dunes of varying size. Whereas smaller bedforms adapt their morphological shape to the oscillating tidal currents, large compound dunes (here: asymmetric tidal dunes) remain stable for periods longer than a tidal cycle. Bedforms constitute a form roughness, that is, hydraulic flow resistance, which has a large-scale effect on tidal asymmetry and, hence, on hydrodynamics, sediment transport, and morphodynamics of estuaries and coastal seas. Flow separation behind the dune crest and recirculation on the steep downstream side result in turbulence and energy loss. Since the energy dissipation can be related to the dune lee slope angle, asymmetric dune shapes induce variable flow resistance during ebb and flood phases. Here, a noncalibrated numerical model has been applied to analyze the large-scale effect of symmetric and asymmetric dune shapes on estuarine tidal asymmetry evaluated by residual bed load sediment transport at the Weser estuary, Germany. Scenario simulations were performed with parameterized bed roughness of symmetric and asymmetric dune shapes and without dune roughness. The spatiotemporal interaction of distinct dune shapes with the main drivers of estuarine sediment and morphodynamics, that is, river discharge and tidal energy, is shown to be complex but substantial. The contrasting effects of flood- and ebb-oriented asymmetric dunes on residual bed load transport rates and directions are estimated to be of a similar importance as the controls of seasonal changes of discharge on these net sediment fluxes at the Lower Weser estuary. This corroborates the need to consider dune-induced directional bed roughness in numerical models of estuarine and tidal environments.  相似文献   

10.
1INTRODUCTIONNon-equilibriumsedimenttransportina"at'Uralstreamwithnon-uniformbedmaterialisasubjectofilltensiveresearch.ThesedimentconcentrationmaybenotequaltothetransportcapacityofflowinanalluvialrivedItgraduallyapproachesequilibriumbydegradationoraggradationinalongdistance.usuallyover100kilometers.Thisproblemisofgreatimportancetothepredictionofthedistributionoferosionordepositionalongachannel,especiallyforlargeriversinChina.Scientistsstudiednonequilibriumsedimenttransportandthemainresul…  相似文献   

11.
1 INTRODUCTION Increasing attention is being given to sedimentation hazards downstream from reservoirs as dams built during the past century accumulate progressively greater volumes of sediment. The sediment storage both decreases reservoir capacity and operating efficiency of the dam, and creates a 搒ediment-shadow?downstream where sediment-starved flows commonly erode channel boundaries and create long-term channel instabilities. Numerous studies have documented downstream channel change…  相似文献   

12.
The morphological active width, defined as the lateral extent of bed material displacement over time, is a fundamental parameter in multi‐threaded gravel‐bed rivers, linking complex channel dynamics to bedload transport. Here, results are presented from five constant discharge experiments, and three event hydrographs, covering a range of flow strengths and channel configurations for which morphological change, bedload transport rates, and stream power were measured in a physical model. Changes in channel morphology were determined via differencing of photogrammetrically‐derived digital elevation models (DEMs) of the model surface generated at regular intervals over the course of ~115 h of experimental runs. Independent measures of total bedload output were made using downstream sediment baskets. Results indicate that the morphological active width increases with total and dimensionless stream power and is strongly and positively correlated with bulk change (total volume of bed material displaced over time) and active braiding intensity (ABI). Although there is considerable scatter due to the inherent variability in braided river morphodynamics, the active width is positively correlated with independent measurements of bedload transport rate. Active width, bulk change, and bedload transport rates were all negligible below a dimensionless stream power threshold value of ~ 0.09, above which all increase with flow strength. Therefore, the active width could be used as a general predictor of bulk change and bedload transport rates, which in turn could be approximated from total and dimensionless stream power or ABI in gravel‐bed braided rivers. Furthermore, results highlight the importance of the active width, rather than the morphological active depth, in predicting volumes of change and bedload transport rates. The results contribute to the larger goals of better understanding braided river morphodynamics, creating large high‐resolution datasets of channel change for model calibration and validation, and developing morphological methods for predicting bedload transport rates in braiding river systems. Copyright © 2018 John Wiley & Sons, Ltd.  相似文献   

13.
Dam break flows and resulting river bed erosion can have disastrous impacts on human safety,infrastructure,and environmental quality.However,there is a lack of research on the mobility of non-uniform sediment mixtures resulting from dam break flows and how these differ from uniform sized sediment.In this paper,laboratory flume experiments revealed that coarse and fine fractions in non-uniform sediment had a higher and a lower bed-load parameter,respectively,than uniform sediments of the same size.Thus,the finer fractions were more stable and the coarser fractions were more erodible in a nonuniform bed compared to a uniform-grained bed.These differences can be explained by the hiding and protrusion of these fractions,respectively.By investigating changes in mobility of the mixed-size fractions with reservoir water levels,the results revealed that at low water levels,when the coarser fractions were only just mobile,the bed-load parameter of the finer fractions was higher than the coarser fractions.The opposite was observed at a higher water level,when a significant proportion of the coarsest fractions was mobilized.The higher protrusion of these grains had an important effect on their mobility relative to the finer grains.The transported sediment on these mixed-sized beds was coarser than the initial bed sediment,and became coarser with an increase in reservoir water level.  相似文献   

14.
The Tongariro Power Development Scheme (TPDS) is used to regulate flow in the headwaters of the largest catchment on the North Island of New Zealand (the Waikato). Two small dams, the Rangipo Dam and the Poutu Intake Dam, were constructed in 1973 and 1983. The flow regime of the river is managed to divert freshes into the power scheme, but allows flows larger than 100 m3 s?1 to be released, to rework and transport sediment through the catchment. Analysis of aerial photos and maps spanning 1928 to 2007, alongside field measurements, show that there have been few hydrogeomorphic adjustments since dam construction. This includes limited changes to channel geometry, channel planform and bed material organization immediately downstream of the dams. In addition, offsite effects are minimal, both 500 m downstream of each dam, and in the more sensitive, less confined reaches in the lower catchment (11 km downstream of the Poutu Intake dam). The limited changes can be attributed to the locations of the dams within reaches characterised by bedrock gorges and confined within terraces. These locations act to flush sediments and impose margins that allow minimal adjustment of the channel. Bed material within this reach is characterised by the presence of a boulder lag. This is sourced from long-term incision into lahar deposits, and acts to limit the rate of incision, creating a steep and stable base upon which active fractions are transported. Just as importantly, significant storage in the low-relief volcanic plateau located in the upper catchment acts to disconnect and store the high sediment yields generated by active volcanic cones in the western sub-catchment upstream of the dams. This limits the rate of sediment supply to regulated reaches. Findings from this study show that analysis of reach-scale controls is essential in framing dam site locations in relation to the distribution of reaches and landscape units across the catchment. In this instance, tributary inputs downstream of the dams do not replenish the sediment and flow removed at the dam locations, as has been observed in other regulated systems. Rather, the river itself is resilient to change and flow variability is well managed allowing geomorphically effective floods to occur. Landscape setting is a key consideration in determining the hydrogeomorphic impact of flow regulation.  相似文献   

15.
We explore the link between channel‐bed texture and river basin concavity in equilibrium catchments using a numerical landscape evolution model. Theory from homogeneous sediment transport predicts that river basin concavity directly increases with bed sediment size. If the effective grain size on a river bed governs its concavity, then natural phenomena such as grain‐size sorting and channel armouring should be linked to concavity. We examine this hypothesis by allowing the bed sediment texture to evolve in a transport‐limited regime using a two grain‐size mixture of sand and gravel. Downstream ?ning through selective particle erosion is produced in equilibrium. As the channel‐bed texture adjusts downstream so does the local slope. Our model predicts that it is not the texture of the original sediment mixture that governs basin concavity. Rather, concavity is linked to the texture of the sorted surface layer. Two different textural regimes are produced in the experiments: a transitional regime where the mobility of sand and gravel changes with channel‐bed texture, and a sand‐dominated region where the mobility of sand and gravel is constant. The concavity of these regions varies depending on the median gravel‐ or sand‐grain size, erosion rate, and precipitation rate. The results highlight the importance of adjustments in both surface texture and slope in natural rivers in response to changes in ?uvial and sediment inputs throughout a drainage network. This adjustment can only be captured numerically using multiple grain sizes or empirical downstream ?ning rules. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

16.
Changes in river regime after the construction of upstream reservoirs   总被引：1，自引：0，他引：1  
This article presents and analyses many years of investigations in China on the fluvial processes downstream of impounding and detention reservoirs. The study covers the change in hydrograph, the recovering of sediment concentration along the river course, the degradation of stream bed, the adjustment of longitudinal profile, the coarsening of bed material, the change in channel width, and the trend of channel pattern variation for alluvial streams downstream of impounding reservoirs. Without confluence of major tributaries, the degradation may extend to a great distance below the dam. In the process of reducing the sediment carrying capacity of the flow to match the diminished sediment supply, the coarsening of bed material is a factor of equal, if not greater, importance as compared with the flattening of channel gradient. In places where the flow has not been sufficiently cut down and the bank is erosive non-resistant, a receding of banklines may take place in concurrence with the deepening of the river bed. Below detention reservoirs, even if the total runoff and sediment supply remain essentially unchanged, the modification of the hydrograph is sufficient to enhance the deterioration of the downstream channel.  相似文献   

17.
Tributaries may either ameliorate or exacerbate the geomorphic and ecologic impacts of flow regulation by altering the flux of water and sediment into the flow‐regulated mainstem. To capture the effects of tributary influences on a flow regulated river, long‐term discharge and cross‐sectional data are used to assess the geomorphic and hydrologic impacts of impoundment. In addition, the use of the short‐lived cosmogenic radioisotope 7Be (half‐life 53·4 days) to link sediment transport dynamics to benthic macroinvertebrate community structure is evaluated. It is found that the 7Be activity of transitional bed load sediment is highly seasonal and reflects both variations in activity of sediment sources and limited sediment residence time within the junction. Benthic communities also exhibit a strong seasonal variability. In the spring, neither the 7Be activity of the sediment, nor benthic communities exhibit clear relationships with sample site location. In contrast, during the late summer the ratio of Ephemeroptera (mayflies)/Trichoptera (caddisflies) decreased significantly below tributary junctions. This decrease in benthic community ratio was driven by increases in caddisfly abundance and was strongly correlated with the presence of recently 7Be tagged transitional bedload sediment. These observations are probably associated with the presence of coarse, stable, and unembedded substrate downstream of tributaries and the rapid turnover of sediment that may also be associated with a rapid flux in nutrients or seston. The results show that tributaries are impacting the flow‐regulated mainstem and that these impacts are reflected in the benthic community structure and in the 7Be activity of transitional bed load sediment. Moreover, the observed reduction in competence and capacity of the mainstem following flood control suggests that these spatial discontinuities may be a consequence of impoundment. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

18.
The 3D numerical model, ECOMSED (open source code), was used to simulate flow and sediment transport in rivers. The model has a long history of successful applications to oceanic, coastal and estuarine waters. Improvements in the advection scheme, treatment of river roughness parameterization and shear stress partitioning were necessary to reproduce realistic and comparable results in a river application. To account for the dynamics of the mobile bed boundary, a model for the bed load transport was included in the code. The model reproduced observed secondary currents, bed shear stress distribution and erosion-deposition patterns on a curved channel. The model also successfully predicted the general flow patterns and sediment transport characteristics of a 1-km long reach of the River Klar?lven, located in the north of the county of V?rmland, Sweden.  相似文献   

19.
In mixed bedrock–alluvial rivers, the response of the system to a flood event can be affected by a number of factors, including coarse sediment availability in the channel, sediment supply from the hillslopes and upstream, flood sequencing and coarse sediment grain size distribution. However, the impact of along-stream changes in channel width on bedload transport dynamics remains largely unexplored. We combine field data, theory and numerical modelling to address this gap. First, we present observations from the Daan River gorge in western Taiwan, where the river flows through a 1 km long 20–50 m wide bedrock gorge bounded upstream and downstream by wide braidplains. We documented two flood events during which coarse sediment evacuation and redeposition appear to cause changes of up to several metres in channel bed elevation. Motivated by this case study, we examined the relationships between discharge, channel width and bedload transport capacity, and show that for a given slope narrow channels transport bedload more efficiently than wide ones at low discharges, whereas wider channels are more efficient at high discharges. We used the model sedFlow to explore this effect, running a random sequence of floods through a channel with a narrow gorge section bounded upstream and downstream by wider reaches. Channel response to imposed floods is complex, as high and low discharges drive different spatial patterns of erosion and deposition, and the channel may experience both of these regimes during the peak and recession periods of each flood. Our modelling suggests that width differences alone can drive substantial variations in sediment flux and bed response, without the need for variations in sediment supply or mobility. The fluctuations in sediment transport rates that result from width variations can lead to intermittent bed exposure, driving incision in different segments of the channel during different portions of the hydrograph. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd  相似文献   

20.
Cohesive sediment dynamics in mountainous rivers is poorly understood even though these rivers are the main providers of fine particles to the oceans. Complex interactions exist between the coarse matrix of cobble bed rivers and fine sediments. Given that fine sediment load in such environments can be very high due to intense natural rainfall or snowmelt events and to man‐induced reservoir or dam flushing, a better understanding of the deposition and sedimentation processes is needed in order to reduce ecohydrological downstream impacts. We tested a field‐based approach on the Arc and Isère alpine rivers combining measurements of erosion and settling properties of river bed deposits before and after a dam flushing, with the U‐GEMS (Gust Erosion Microcosm System) and SCAF (System Characterizing Aggregates and Flocs), respectively. These measurements highlight that critical shears, rates of erosion, settling velocities and propensity of particles to flocculate are highly variable in time and space. This is reflective of the heterogeneity of the hydrodynamic conditions during particle settling, local bed roughness, and nature and size of particles. Generally the deposits were found to be stable relative to what is measured in lowland rivers. It was, however, not possible to make a conclusive assessment of the extent to which the dynamics of deposits after reservoir flushing were different from those settled after natural events. The absence of any relationships between erosion and deposition variables, making it impossible to predict one from another, underlined the need to measure all of them to have a full assessment of the fine sediment dynamics and to obtain representative input variables for numerical models. While the SCAF was found to be effective, an alternative to the U‐GEMS device will have to be found for the erodibility assessment in cobble bed rivers, in order to make more rapid measurements at higher shears. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

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

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