首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Q. He  D. E. Walling 《水文研究》1998,12(7):1079-1094
River floodplains have been widely recognized as important sinks for storing suspended sediment and associated contaminants transported by river systems. The grain size composition of floodplain deposits exerts an important influence on contaminant concentrations, and commonly exhibits significant spatial variability in response to the dynamic nature of overbank flow and sediment transport. Information on the spatial variability of the grain size composition of overbank deposits is therefore essential for developing an improved understanding of the processes controlling sediment transport on floodplains, and for investigating the fate of sediment-associated contaminants. Such information is also important for validating existing floodplain sedimentation models. This paper reports the results of a study aimed at investigating the spatial variability of the grain size composition of floodplain sediments at different spatial scales, through analysis of surface sediment samples representative of contemporary floodplain deposits collected from frequently inundated floodplain sites on five British lowland rivers. Significant lateral and downstream variations in the grain size composition of the sediment deposits have been identified in the study reaches. An attempt has been made to relate the observed spatial distribution of the grain size composition of the overbank deposits to the local floodplain geometry and topography. The importance of the particle size characteristics of the suspended sediment transported by the rivers in influencing the spatial variability of the grain size composition of the overbank sediments deposited on these floodplains is also considered. © 1998 John Wiley & Sons, Ltd.  相似文献   

2.
This article presents results from an investigation of the hydraulic characteristics of overbank flows on topographically‐complex natural river floodplains. A two‐dimensional hydraulic model that solves the depth‐averaged shallow water form of the Navier–Stokes equations is used to simulate an overbank flow event within a multiple channel reach of the River Culm, Devon, UK. Parameterization of channel and floodplain roughness by the model is evaluated using monitored records of main channel water level and point measurements of floodplain flow depth and unit discharge. Modelled inundation extents and sequences are assessed using maps of actual inundation patterns obtained using a Global Positioning System, observational evidence and ground photographs. Simulation results suggest a two‐phase model of flooding at the site, which seems likely to be representative of natural floodplains in general. Comparison of these results with previous research demonstrates the complexity of overbank flows on natural river floodplains and highlights the limitations of laboratory flumes as an analogue for these environments. Despite this complexity, frequency distributions of simulated depth, velocity and unit discharge data closely follow a simple gamma distribution model, and are described by a shape parameter (α) that exhibits clear systematic trends with changing discharge and floodplain roughness. Such statistical approaches have the potential to provide the basis for computationally efficient flood routing and overbank sedimentation models. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

3.
The aim of this paper is to quantify peakflow attenuation and/or amplification in a river, investigating lateral flow from the intermediate catchment during floods. This is a challenge for the study of the hydrological response of permeable/intermittent streams, and our contribution refers to a modelling framework based on the inverse problem for the diffusive wave model applied in a karst catchment. Knowing the upstream and downstream hydrographs on a reach between two stations, we can model the lateral one, given information on the hydrological processes involved in the intermediate catchment. The model is applied to 33 flood events in the karst reach of the Iton River in French Normandy where peakflow attenuation is observed. The monitored zone consists of a succession of losing and gaining reaches controlled by strong surface‐water/groundwater (SW/GW) interactions. Our results show that despite a high baseflow increase in the reach, peakflow is attenuated. Model application shows that the intensity of lateral outflow for the flood component is linked to upstream discharge. A combination of river loss and overbank flow for highest floods is proposed for explaining the relationships. Our approach differentiates the role of outflow (river loss and overbank flow) and that of wave diffusion on peakflow attenuation. Based on several sets of model parameterization, diffusion is the main attenuation process for most cases, despite high river losses of up to several m3/s (half of peakflow for some parameterization strategies). Finally, this framework gives new insight into the SW/GW interactions during floods in karst basins, and more globally in basins characterized by disconnected river‐aquifer systems.  相似文献   

4.
Part 1 of this study assessed the effectiveness of dispersion models based around random particle tracking (RPT) applied an overbank flow in the Flood Channel Facility at HR Wallingford. Several of the RPT models that were used predicted the mixing behaviour of the tracer in the complex channel flow to within an acceptable accuracy. If there is no uniqueness in the model structure which can be used to represent the true system to within the limitations of the available observations, then this implies an inherent degree of ambiguity in our knowledge of the physically based model structure. This suggests that we should be less forthright in the optimisation of each individual model structure, and perhaps investigate more of the parameter combinations for each model which yield feasible simulations of the system. An alternative fuzzy calibration technique is introduced which avoids the optimisation process and takes account of uncertainties in the model structure, parameter sets and observed data in prediction.  相似文献   

5.
Relative to those at sub‐bankfull flow, hydraulic conditions at overbank flow, whether in the channel or on the floodplain, are poorly understood. Here, velocity conditions are analysed over an unusually wide range of flows in the arid zone river of Cooper Creek with its complex system of anastomosing channels and large fluctuations in floodplain width. At‐a‐station hydraulic geometry relationships reveal sharp discontinuities in velocity at the inbank–overbank transition, the nature of the discontinuity varying with the degree of flow confinement and the level of channel–floodplain interaction. However, despite inter‐sectional differences, velocities remain modest throughout the flow range in this low‐gradient river, and the large increases in at‐a‐station discharge are principally accommodated by changes in cross‐sectional area. Velocity distribution plots suggest that within‐channel conditions during overbank flow are characterized by a central band of high velocity which penetrates far toward the bed, helping to maintain already deep cross‐sections. Floodplain resistance along Cooper Creek is concentrated at channel bank tops where vegetation density is highest, and the subsequent flow retardation is transmitted across the surface of the channels over distances as large as 50–70 m. The rough floodplain surface affects flood wave transmission, producing significant decreases in wave speeds downstream. The character of the wave‐speed–discharge relationship also changes longitudinally, from log–linear in the upper reaches to nonlinear where the floodplain broadens appreciably. The nonlinear form is similar in several respects to relationships proposed for more humid rivers, with flood wave speed reaching an intermediate maximum at about four‐fifths bankfull discharge before decreasing to a minimum at approximately Q2·33. It does not regain the value at the intermediate maximum until the 10 year flood, by which time floodplain depths have become relatively large and broad floodways more active. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

6.
The channel boundary conditions along the Lower Yellow River (LYR) have been altered significantly since the 1950s with the continual reinforcement and construction of both main and secondary dykes and river training works. To evaluate how the confined complex channel–floodplain system of the LYR responds to floods, this study presents a detailed investigation of the relationship between the tempo‐spatial distribution of sedimentation/erosion and overbank floods occurred in the LYR. For large overbank floods, we found that when the sediment transport coefficient (ratio of sediment concentration of flow to flow discharge) is less than 0.034, the bankfull channel is subject to significant erosion, whereas the main and secondary floodplains both accumulate sediment. The amount of sediment deposited on the main and secondary floodplains is closely related to the ratio of peak discharge to bankfull discharge, volume of water flowing over the floodplains, and sediment concentration of overbank flow, whereas the degree of erosion in the bankfull channel is related to the amount of sediment deposited on the main and secondary floodplains, water volume, and sediment load in flood season. The significant increase in erosion in the bankfull channel is due to the construction of the main and secondary dykes and river training works, which are largely in a wide and narrow alternated pattern along the LYR such that the water flowing over wider floodplains returns to the channel downstream after it drops sediment. For small overbank floods, the bankfull channel is subject to erosion when the sediment transport coefficient is less than 0.028, whereas the amount of sediment deposited on the secondary floodplain is associated closely with the sediment concentration of flow. Over the entire length of the LYR, the situation of erosion in the bankfull channel and sediment deposition on the main and secondary floodplains occurred mainly in the upper reach of the LYR, in which a channel wandering in planform has been well developed.  相似文献   

7.
Field studies suggest that a cohesive floodplain is a necessary condition for meandering in contrast to braided rivers. However, it is only partly understood how the balance between floodplain construction by overbank deposition and removal by bank erosion and chutes leads to meandering. This is needed because only then does a dynamic equilibrium exist and channels maintain meandering with low width–depth ratios. Our objective is to understand how different styles of floodplain formation such as overbank deposition and lateral accretion cause narrower channels and prevent chute cutoffs that lead to meandering. In this study we present two experiments with a self‐forming channel in identical conditions, but to one we added cohesive silt at the upstream boundary. The effect of cohesive silt on bank stability was tested in auxiliary bank erosion experiments and showed that an increase in silt reduced erosion rates by a factor of 2. The experiment without silt developed to a braided river by continuous and extensive shifting of multiple channels. In contrast, in the meandering river silt deposits increased bank stability of the cohesive floodplain and resulted in a reduction of chute cutoffs and increased sinuosity by continuous lateral migration of a single channel. Overbank flow led to deposition of the silt and two styles of cohesive floodplain were observed: first, overbank vertical‐accretion of silt, e.g. levee, overbank sedimentation or splays; and second, lateral point bar accretion with silt on the scrolls and in the swales. The first style led to a reduction in bank erosion, while the second style reduced excavation of chutes. We conclude that sedimentation of fine cohesive material on the floodplain by discharge exceeding bankfull is a necessary condition for meandering. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

8.
A palynological approach was used to estimate overbank deposition rates in a forested catchment affected by logging. The palynological approach uses downcore variations in total fossil pollen and fossil pollen assemblage to calculate rates of overbank deposition and has a distinct advantage over radioisotopic approaches in that it is not limited by radioactive decay. Using this approach, we determined that overbank deposition rates increased over 400 per cent within years of logging events and that the increased rates persisted for less than 4 years. After logging‐induced deposition peaked, overbank deposition decreased over 60 per cent relative to the pre‐logging background values. The decreased deposition rates persisted for over 40 years. The immediate effect of logging in this catchment was to induce mass‐wasting events in hollows that produced rapidly travelling sediment pulses. In the subsequent recovery period, reduced sediment loading occurred as a result of a reduction in the volume of sediment available for transport. The reduction in sediment load led to a reduction in overbank deposition rates until subsequent logging disturbances destabilized and emptied other hollows. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

9.
Abstract

Mathematical modelling of overbank inundation and flows faces many problems and is still in its infancy. Work to date has generally been restricted to small reaches. Large-scale models based on longer reaches of river channel are likely to be of greater value for engineering and flood plain management purposes, but the problems associated with the transition from small to large scales need to be assessed. A large-scale finite element model, RMA-2, has been applied to the flood plain of the lower reaches of the River Culm in southeast Devon, UK. Patterns of radiocaesium accumulation by overbank accretion during flood water inundation were used to assess the potential of using such models for explaining sedimentation rates and patterns. A strong correlation was found between values of the 137Cs inventory and surface concentration and the predicted flood water patterns derived using the RMA-2 model. Except where recession pondage occurs, an inverse relationship existed between 137Cs deposition and water depth. However, the discretization model developed cannot presently cope with large-scale compartmentalization of flows by barriers to flow and small-scale local features, such as ditches crossing the flood plain and the microtopography of the flood plain. This study appraises the potential for using the RMA-2 model to predict patterns of overbank deposition and represents an initial stage in the development of an integrated model of hydraulic and sediment dynamics.  相似文献   

10.
This paper outlines a numerical model for the prediction of floodplain inundation sequences, overbank deposition rates and deposit grain size distributions. The model has two main components: first, a simplified hydraulic scheme which predicts floodwater flow depths and velocities, and second, a sediment transport element which employs a mass balance relation describing suspended sediment dispersion by convective and diffusive processes and sediment deposition as a function of particle settling rates. These relationships are solved numerically on a finite difference grid that accurately replicates the complex topographic features typical of natural river floodplains. The model is applied to a 600 m reach of the River Culm, Devon, U.K. using data derived from a range of field and laboratory techniques. Continuous records of river stage and suspended sediment concentration provide the model's upstream boundary input requirements. These are supplemented by measurements of the in situ settling characteristics of the suspended sediment load. The model's sediment transport component is calibrated with the aid of a dataset of measured overbank deposition amounts derived from flood events over a 16 month period. The model is shown to predict complicated floodwater inundation sequences and patterns of suspended sediment dispersion and deposition, which are largely a product of the complex topography of the floodplain. These results compare favourably with observations of overbank processes and are an improvement over those of previous models which have employed relatively simple representations of floodplain geometry. © 1997 by John Wiley & Sons, Ltd.  相似文献   

11.
In a compound meandering channel, patterns of flow structures and bed variations change with increasing water depth owing to complex momentum exchange between high-velocity flow in a main channel and low-velocity flows in flood plains. We have developed a new quasi-three-dimensional model without the shallow water assumption, i.e., hydrostatic pressure distribution; our method is known as the general bottom velocity computation (BVC) method. In this method, a set of depth-integrated equations, including depth-integrated momentum and vorticity equations, are prepared for evaluating bottom velocity and vertical velocity distributions. The objective of this study is to develop a bed variation calculation method for both single and compound meandering channels by using the BVC method coupled with a sediment transport model. This paper shows that the BVC method can reproduce the pattern change of bed variation in a compound meandering channel flow with increasing relative depth. The variation in sediment transport rate due to overbank flow is explained by experimental and computational results.  相似文献   

12.
Particulate organic matter (POM) transiting through rivers could be lost to overbank storage, stored in‐channel, added to by erosion or autochthonous production, or turned over to release greenhouse gases to the atmosphere (either while in the water column or while stored in the channel). In the UK, a net loss of POM across catchments has been recorded, and the aim here was to investigate the balances of processes acting on the POM. This study considered records of suspended sediment and POM flux in comparison to stream flow, velocity, stream power, and residence time for the River Trent (English Midlands, 8,231 km2). We show that for the lower two thirds (106 km) of the River Trent, 2% is lost to overbank storage; 10% is lost to the atmosphere in the water column; and 31% is turned over while in temporary storage. Permanent in‐channel storage is negligible, and for the lower course of the river, material stored in‐channel will have a residence time of the order of hundreds of days between the last flood hydrograph of one winter and the first winter storm of the next winter (usually in the same calendar year). When considered at the scale of the UK, 1% POM in transit would be lost to overbank sedimentation; 5% turned over in the water column, and 14% turned over while in temporary storage. In the upper third of the study river channel, there is insufficient stream power to transport sediment and so in‐channel storage or in‐channel turnover over to the atmosphere dominate. The in‐channel processes of the River Trent do not conform to that expected for river channels as the headwaters are not eroding or transporting sediment. Therefore, the source of sediment must be lower down the channel network.  相似文献   

13.
Valley-fill ignimbrites of the Middle Laacher See Tephra (MLST) in the proximal depositional fan south of Laacher See volcano are laterally continuous with fine-grained overbank-facies deposits, 0.5–1.5 m thick, covering higher elevations and interfluves between the paleovalleys. The overbank deposits consist of up to 12 ash layers, each 4–10 cm thick, which show internal structures typical of ash-flow transport, such as poor sorting, reverse size-grading of pumice, local normal grading of coarse ash-sized lithics above a fine-grained basal layer, cross-stratification behind obstacles, and erosional unconformities. Thickness, median grain-size, and number of individual layers decrease systematically with distance from the vent. Overbank ash layers thicken at the valley slopes and form discrete valley-fill ignimbrite flow units in the paleochannels with median grain size increasing from 63 m to 350 m. Toward the center of paleochannels, however, the well-defined overbank facies is obscured by mutual erosion of individual flow units. Overlapping data fields in ternary grain-size variation diagrams indicate the overbank facies to have evolved from ash flows chiefly through depletion of lapilli and coarse ash. Overbank-facies ash layers do not represent dust layers resulting from elutriation clouds of ash flows (co-ignimbrite ash) or surge deposits developed on higher ground due to low concentration of solids. They are similar in some parameters to Taupo-type ignimbrite veneer deposits, but are interpreted differently. The thin, fine-grained, Laacher See ash layers are thought to have been deposited from diluted portions of the flow proceeding directly from the eruption column while the main pyroclastic flows were confined to the paleovalleys radiating away from the eruptive center. The wide distribution of the thin ash layers is attributed to the balance of deposition from various flow parts and turbulent entraining and heating of ambient air that sustained sufficient mobility of the diluted flows to spread across hills and level ground.  相似文献   

14.
Slow earth sliding is pervasive along the concave side of Red River meanders that impinge on Lake Agassiz glaciolacustrine deposits. These failures form elongated, low‐angled (c. 6 to 10°) landslide zones along the valleysides. Silty overbank deposits that accumulated during the 1999 spring freshet extend continuously along the landslide zones over hundreds of metres and aggraded the lower slopes over a distance 50 to 80 m from the channel margin. The aggradation is not obviously related to meander curvature or location within a meander. Along seven slope profiles surveyed in 1999 near Letellier, Manitoba, the deposits locally are up to 21 cm thick and generally thin with increasing distance from, and height above, the river. Local deposit thickness relates to distance from the channel, duration of inundation of the landslide surface, mesotopography, and variations in vegetation cover. Immediately adjacent to the river, accumulated overbank deposits are up to 4 m thick. The 1999 overbank deposits also were present along the moderately sloped (c. 23 to 27°) concave banks eroding into the floodplain, but the deposits are thinner (locally up to c. 7 cm thick) and cover a narrower area (10 to 30 m wide) than the deposits within the landslide zones. Concave overbank deposition is part of a sediment reworking process that consists of overbank aggradation on the landslide zones, subsequent gradual downslope displacement from earth sliding, and eventually reworking by the river at the toe of the landslide. The presence of the deposits dampens the outward migration of the meanders and contributes to a low rate of contemporary lateral channel migration. Concave overbank sedimentation occurs along most Red River meanders between at least Emerson and St. Adolphe, Manitoba. © Her Majesty the Queen in right of Canada.  相似文献   

15.
Intensive field monitoring of a reach of upland gravel‐bed river illustrates the temporal and spatial variability of in‐channel sedimentation. Over the six‐year monitoring period, the mean bed level in the channel has risen by 0·17 m with a maximum bed level rise of 0·5 m noted at one location over a five month winter period. These rapid levels of aggradation have a profound impact on the number and duration of overbank flows with flood frequency increasing on average 2·6 times and overbank flow time increasing by 12·8 hours. This work raises the profile of coarse sediment transfer in the design and operation of river management, specifically engineering schemes. It emphasizes the need for the implementation of strategic monitoring programmes before engineering work occurs to identify zones where aggradation is likely to be problematic. Exploration of the sediment supply and transfer system can explain patterns of channel sedimentation. The complex spatial, seasonal and annual variability in sediment supply and transfer raise uncertainties into the system's response to potential changes in climate and land‐use. Thus, there is a demand for schemes that monitor coarse sediment transfer and channel response. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

16.
This paper investigates the impact of a 1000‐year flood in August 2002 on floodplains and valley morphology of an Austrian mixed alluvial bed rock river. Discharges with a recurrence interval between 500 and 2000 years caused distinctive overbank scouring and material deposition in the floodplains. After the 1000‐year flood, those morphologically affected areas were at random intervals documented over the whole longitudinal profile. In addition to overbank erosion in curved sections (cut‐offs), the river bed locally widened, floodplain stripping occurred and local overbank scours were documented along straight parts of the river. A hydrodynamic‐numerical model, combined with field measurements, was used to analyse the cause of these erosional landforms. Based on the modelled hydraulic conditions for a one‐year flood (30–78 m3 s–1) and the catastrophic 2002 event (700–800 m3 s–1), the numerical results allowed a cause‐effect study with 19 parameters. Deterministic and statistical analysis (ANOVA, discriminant analysis) showed that the morphodynamic effects of the 2002 flood were influenced by the variability of valley morphology of the Kamp River, which led partially to supercritical flow during flood constriction. These processes were in some cases also anthropogenically influenced. Lateral constriction and expansion of the valley geometry over short distances led to scouring and aggradation within the inundated areas during the event. These morphological features were therefore responsible for the elongated scour holes in the floodplains. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

17.
Measurements from a fixed‐bed, Froude‐scaled hydraulic model of a stream in northeastern Vermont demonstrate the importance of forested riparian vegetation effects on near‐bank turbulence during overbank flows. Sections of the prototype stream, a tributary to Sleepers River, have increased in channel width within the last 40 years in response to passive reforestation of its riparian zone. Previous research found that reaches of small streams with forested riparian zones are commonly wider than adjacent reaches with non‐forested, or grassy, vegetation; however, driving mechanisms for this morphologic difference are not fully explained. Flume experiments were performed with a 1:5 scale, simplified model of half a channel and its floodplain, mimicking the typical non‐forested channel size. Two types of riparian vegetation were placed on the constructed floodplain: non‐forested, with synthetic grass carpeting; and forested, where rigid, randomly distributed, wooden dowels were added. Three‐dimensional velocities were measured with an acoustic Doppler velocimeter at 41 locations within the channel and floodplain at near‐bed and 0·6‐depth elevations. Observations of velocity components and calculations of turbulent kinetic energy (TKE), Reynolds shear stress and boundary shear stress showed significant differences between forested and non‐forested runs. Generally, forested runs exhibited a narrow band of high turbulence between the floodplain and main channel, where TKE was roughly two times greater than TKE in non‐forested runs. Compared to non‐forested runs, the hydraulic characteristics of forested runs appear to create an environment with higher erosion potential. Given that sediment entrainment and transport can be amplified in flows with high turbulence intensity and given that mature forested stream reaches are wider than comparable non‐forested reaches, our results demonstrated a possible driving mechanism for channel widening during overbank flow events in stream reaches with recently reforested riparian zones. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

18.
The basaltic Kaupulehu 1800–1801 lava flow of Hualalai Volcano, Hawaii contains abundant ultramafix xenoliths. Many of these xenoliths occur as bedded layers of semi-rounded nodules, each thinly coated with a veneer (typically 1 mm thick) of lava. The nodule beds are analogous to cobble deposits of fluvial sedimentary systems. Although several mechanisms have been proposed for the formation of the nodule beds, it was found that, at more than one locality, the nodule beds are overbank levee deposits. The geological occurrence of the nodules, certain diagnostic aspects of the flow morphology and consideration of the inferred emplacement process indicate that the Kaupulehu flow had an exceptionally low viscosity on eruption and that the flow of the lava stream was extremely rapid, with flow velocities of at least 10 m s-1 (more than 40 km h-1). This flow is the youngest on Hualalai Volcano and future eruptions of a similar type would pose considerable hazard to life as well as property.  相似文献   

19.
Dominant discharge may be defined as that discharge which transports most bed sediment in a stream that is close to steady-state conditions. The concept is examined in relation to two single thread gravel-bedded streams. One stream is alluvial and free to adjust its geometry whilst in the other, channel capacity and form are partially constrained by cohesive till-banks and a heavily compacted bed. The total quantity of bedload and its calibre were measured for every flood over a six year period, so that the relationship between the grain-size of bedload and the most effective discharge could be examined in the context of thresholds for channel change. The dominant discharge concept was applicable to the alluvial stream in that the bankfull value is an effective discharge for maintaining channel capacity. The concept applied less well to the ‘non-alluvial’ stream. Although in both streams the bankfull value was exceeded for less than 0.34 per cent of the time, overbank flows are important in instigating channel change. It is only during overbank flows that the largest bed material is entrained in quantity. For within-channel flows a threshold separates flows which winnow fine matrix from those which entrain the finer bed gravels. This threshold occurred at 60 per cent bankfull. It was concluded that the dominant discharge concept can be applied to streams close to steady-state which are alluvial, competent, and free to adjust their boundaries. An important proviso is that two channel-stability domains can be recognized. These domains represent channel maintenance and channel adjustment and are defined by intrinsic thresholds in the bed material entrainment function.  相似文献   

20.
In a special opportunity, detailed measurements of the flow in an overbank flow in the Flood Channel Facility at HR Wallingford were used in conjunction with tracer test data to assess the effectiveness of dispersion models based around random particle tracking (RPT). Ten different RPT models based on different assumptions and levels of information about the nature of the Lagrangian velocity field were investigated. Multiple simulations were used to calibrate variable parameters controlling the average magnitude of the perturbations for each model by comparison with observed concentrations at one cross-section. The calibrated models were then used to predict concentration distributions further downstream. Several of the calibrated models showed close agreement between observed and predicted concentration distributions. The most complex models using the most information about the velocity distributions were no better (and in some cases worse) in prediction than the simplest models investigated. It would appear that our knowledge of the system, despite the quality of the experiments, is too uncertain to infer a precise model structure.  相似文献   

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

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