首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 23 毫秒
1.
《国际泥沙研究》2020,35(2):203-216
Accurate evaluation and prediction of bedload transport are crucial in studies of fluvial hydrodynamic characteristics and river morphology.This paper presents a one-dimensional numerical model based on the one-dimensional lateral distribution method(1 D-LDM) and six classic bedload transport formulae that can be used to simulate hydrodynamic characteristics and bedload transport discharge in cross sections.Two gravel-bed rivers,i.e.the Danube River located approximately 70 km downstream from Bratislava in Slovakia and the Tolten River in south of Chile are used as examples.In the 1 D-LDM,gravity,bed shear stress,turbulent diffusion,and secondary flow are included to allow for accurate predictions of flow velocity and consequently bed shear stress in the cross sections.Six classic formulae were applied to evaluate the non-dimensional bedload transport rate,and the bedload transport discharge through a river cross section is obtained by integrating the bedload transport rate over the width of the cross section.The results show that the root mean square error(RMSE) and mean absolute error(MAE) of velocity and water discharge were less than 8% of the observed magnitude,while the correlation coefficient between model predictions and observations was close to unity.The formulae proposed by Ashida and Michiue(1972),in which particle collision with the bed is taken into account,and by Camenen and Larson(2005),which allows for yielding a non-zero bedload transport rate even when the bed shear stress is smaller than the critical bed shear stress value,appeared to be more appropriate for predicting the observed bedload transport rate in the studied cross sections of two gravel-bed rivers.If non-uniform sediment mixtures were considered,the bedload transport discharge through a cross-section could change considerably by up to 22.5% of the observed magnitude.The relations proposed by Ashida and Michiue(1972) and Egiazaroff(1965) for parameterizing the hiding factor yielded more realistic model predictions in comparison with observations for the measured data set collected for the Tolten River,while the one proposed by Wilcock and Crowe(2003) performs the best for the data set measured for the Danube River.  相似文献   

2.
Taking the width/depth ratio of a river channel as an independent variable, a variational analysis of basic flow relationships shows that alluvial‐channel flow adjusts channel geometry to achieve stationary equilibrium when the condition of maximum flow efficiency (MFE) is satisfied. As a test of the veracity of MFE and to examine if this theory of self‐adjusting channel morphodynamics can be practically applied to large river systems, this study examines the degree of correspondence between theoretically determined equilibrium channel geometries and actual measurements along the middle and lower Yangtze River. Using four different forms of the Meyer‐Peter and Müller bedload relation and relations of flow continuity and resistance we show that the Meyer‐Peter and Müller bedload relation modified on the basis of MFE theory predicts channel dimensions most accurately when applied to the middle and lower Yangtze River. This provides convincing evidence supporting MFE equilibrium theory. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

3.
Bagnold developed his formula for bedload transport over several decades, with the final form of the relation given in his 1980 paper. In this formula, bedload transport rate is a function of stream power above some threshold value, depth and grain size. In 1986, he presented a graph which illustrated the strength of his relation. A double‐log graph of bedload transport rate, adjusted for depth and grain size, versus excess stream power was shown to collapse along a line having a slope of 1·5. However, Bagnold based his analyses on limited data. In this paper, the formula is re‐examined using a large data set in order to define the most consistent empirical representation, and dimensional analysis is performed to seek a rationalization of the formula. Functional analysis is performed for the final version of the equation defined by Bagnold to determine if the slope of 1·5 is preserved and to assess the strength of the relation. Finally, relations between excess stream power and bedload transport are examined for a fixed slope of 1·5 to assess the performance of various depth and grain size adjustment factors. The rational scaling is found to provide the best result. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

4.
Only comparatively few experimental studies have been carried out to investigate the performance of the HEC-6 river morphological model. The model was developed by the Hydrologic Engineering Center of the US Army Corps of Engineers. In this study, experiments were carried out in a 20 m long concrete flume 0.6 m wide with varying rectangular cross-sections. The channel bed is paved with uniform sand of D50 = 0.9 mm and D90 = 1.2 mm within the test reach of 12 m. Two types of experiments were carried out with sediment transport, one under steady uniform flow and another under steady non-uniform flow conditions. Nine steady uniform flow experiments were carried out to compare the measured equilibrium relationship of flow and sediment transport rate with two bedload formulae, namely, Du Boys and Meyer–Peter and Muller, and with three total load formulae, namely, Toffaleti, Laursen and Yang. It was found that even though the sediment transport consists of a certain portion of bedload, the total load formulae give satisfactory results and better agreement than the two bedload formulae. Five steady non-uniform flow experiments were carried out under various conditions of varying bed profile and channel width and also with sediment addition and withdrawal. The measured transient water surface and bed profiles are compared with the computed results from the HEC-6 model. It was found that the Toffaleti and Yang total load formulae used in the HEC-6 model give the most satisfactory prediction of actual bed profiles under various conditions of non-uniform flow and sediment transport. The effects of Manning's n, variations of sediment inflow, various sediment transport formulae, sediment grain size and the model numerical parameters, i.e. distance interval Δx and numerical weighting factor, on the computed water surface and bed profiles were determined. It was found that the selection of the sediment transport formulae has the most significant effect on the computed results. It can be concluded that the HEC-6 model can predict satisfactorily a long-term average pattern of local scour and deposition along a channel with either a small abrupt change in geometry or gradually varying cross-sections. However, the accuracy of the model prediction is reduced in the regions where highly non-uniform flow occurs.  相似文献   

5.
Automatic and continuously recording samplers are deployed in a Hertfordshire gravel-bed stream to show that bedload transport is related to stream power. The pattern is similar to that already established for North American channels but, because the record is so detailed, it is possible to identify the cause of the considerable scatter that is normal in such relationships. A major factor is the occurrence of rhythmic pulses in bedload discharge that are not matched by similar fluctuations in hydraulic variables. It is suggested that these pulses reflect downstream differences in the concentration of mobile particles in a slow-moving traction carpet, and that they may be likened to kinematic waves. The record also reveals that the threshold of sediment transport—always presumed hithero to be associated with incipient motion—is related to the cessation of bedload transport in a river flood. Indeed, the mean value of stream power at the finish of bedload transport is only 20 percent of that prevailing at the moment of incipient sediment motion. Because of this, there is an inevitably poor correlation between actual bedload transport rates and those predicted by bedload equations which rely upon a single traction threshold. These new data show that the general inverse relationship between bedload discharge and water-depth : grain-size ratio proposed by Bagnold (1977, 1980) is not universal. Transport efficiency for this gravel-bed stream is typically 0.05 per cent of available stream power, which compares with 1.6 per cent for a river moving both gravel and sand, and 5 per cent for another channel where bedload is composed predominantly of sand-sized particles. It is argued that coarse and fine-grained alluvial channels may need to be considered separately. By allowing for differences in traction threshold at the beginning and end of bedload events, and by averaging bedload discharge flood by flood in order to smooth out the effect of pulses, it is possible to achieve a reasonably good prediction of average bedload transport rate in terms of stream power.  相似文献   

6.
  相似文献   

7.
It is increasingly recognized that effective river management requires a catchment scale approach. Sediment transport processes are relevant to a number of river functions but quantifying sediment fluxes at network scales is hampered by the difficulty of measuring the variables required for most sediment transport equations (e.g. shear stress, velocity, and flow depth). We develop new bedload and total load sediment transport equations based on specific stream power. These equations use data that are relatively easy to collect or estimate throughout stream networks using remote sensing and other available data: slope, discharge, channel width, and grain size. The new equations are parsimonious yet have similar accuracy to other, more established, alternatives. We further confirm previous findings that the dimensionless critical specific stream power for incipient particle motion is generally consistent across datasets, and that the uncertainty in this parameter has only a minor impact on calculated sediment transport rates. Finally, we test the new bedload transport equation by applying it in a simple channel incision model. Our model results are in close agreement to flume observations and can predict incision rates more accurately than a more complicated morphodynamic model. These new sediment transport equations are well suited for use at stream network scales, allowing quantification of this important process for river management applications. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

8.
This paper presents results of recent measurements of sand transport made in Chioggia inlet as part of an extensive monitoring programme in the Venetian inlets. Measurements were made in order: (1) to define a relationship between sand transport magnitude and tidal flow; (2) to derive the thresholds for sand transport; (3) to identify the dominant modes of transport; (4) to evaluate the concentration profiles of sand within the benthic boundary layer; (5) to compare bedload transport observations with model predictions using existent bedload formulae; and (6) to produce yearly estimates of bedload transport across the inlet. The vertical distribution of sand in the water column was sampled using modified Helley–Smith bedload samplers at three sites. Transport was found to vary according to the flow and bed grain size, with considerable temporal and spatial variability. A difference of up to three orders of magnitude in transport was observed through the inlet, with higher transport rates measured on the seaward part. The dominant mode of transport in the central inlet was suspension, while bedload was dominant in the mouths. The measured profiles of sand concentration varied with the tidal stage and seabed grain size according to the Rouse parameter (R). R was high at the inlet mouths (1<R<2), indicative of a well-developed bedload layer. The inverse movability number (Ws/U*) was also higher at these sites and appeared to be grain size dependant. Formulae for bedload transport were tested against field data; stochastic methods such as Einstein–Brown, Engelund–Hansen and Van Rijn produce the best fits. The coupled model SHYFEM-Sedtrans05 appears to simulate well observed transport for most conditions of flow. Long-term bedload predictions indicate a dominant export of sand, with a yearly average of 4500 m3.  相似文献   

9.
A comprehensive monitoring programme focusing on bedload transport behaviour was conducted at a large gravel‐bed river. Innovative monitoring strategies were developed during five years of preconstruction observations accompanying a restoration project. A bedload basket sampler was used to perform 55 cross‐sectional measurements, which cover the entire water discharge spectrum from a 200‐year flood event in 2013 to a rare low flow event. The monitoring activities provide essential knowledge regarding bedload transport processes in large rivers. We have identified the initiation of motion under low flow conditions and a decrease in the rate of bedload discharge with increasing water discharge around bankfull conditions. Bedload flux strongly increases again during high flood events when the entire inundation area is flooded. No bedload hysteresis was observed. The effective discharge for bedload transport was determined to be near mean flow conditions, which is therefore at a lower flow discharge than expected. A numerical sediment transport model was able to reproduce the measured sediment transport patterns. The unique dataset enables the characterisation of bedload transport patterns in a large and regulated gravel‐bed river, evaluation of modern river engineering measures on the Danube, and, as a pilot project has recently been under construction, is able to address ongoing river bed incision, unsatisfactory ecological conditions for the adjacent national park and insufficient water depths for inland navigation. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

10.
We report on bedload transport observations using piezoelectric bedload impact sensors (PBIS), an indirect method of estimating the volume of bedload transport of coarse sediment. The PBIS device registers vibrations produced by bedload (particle diameter >～20 mm) and records the signal as a sum of the number of impulses per time. Sediment transport at the Erlenbach stream has been continuously monitored with a PBIS array starting in 1986. The sensor array spans the width of an entire cross‐section and is mounted flush with the surface of a check dam immediately upstream of a sediment retention basin. We compare PBIS data with long‐term sedimentation records obtained from repeated surveys of material stored in the sediment retention basin, with artificial sediment input under controlled conditions in the field, and also with laboratory experiments. The rate of bedload transport is proportional to the number of impacts on the sensor per unit time. The reliability of the calibration relationship increases with the length of the observation period, e.g. for higher numbers of impacts and larger bedload volumes. Sediment volumes for individual flood events estimated with the PBIS method are in agreement with volumes estimated using an independent empirical method based on the effective runoff volume of water, the peak water discharge, and the critical discharge for the onset of sediment transport. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

11.
Laboratory flume experiments were carried out to evaluate the effect of particle density on bedload transport of sand‐sized particles and the effect of a suspended load of clay particles (kaolinite) on bedload transport of sand‐sized particles in rill flow conditions. Three materials in the range 400–600 µm were selected to simulate bedload transport of primary particles and aggregates: sand (2650 kg/m3), crushed brick (2450 kg/m3) and anthracite (1300–1700 kg/m3). In the two first experiments, two different methods were applied to determine bedload transport capacity of coarse particles for various conditions of flow discharge (from 2 to 15 L/min) and slope (2.2, 3 and 4%). In the third experiment, clear water was replaced with kaolinite–water mixture and bedload transport capacity of crushed brick particles was determined for a 4% slope and different concentrations of kaolinite (0, 7, 41 and 84 g/L). The results showed that bedload transport increased significantly with the decrease in particle density but the effect of particle density on transport rates was much less important than flow discharge. Velocity measurements of clear flow, flow mixed with coarse particles and coarse particles confirmed the existence of a differentiation between suspended load and bedload. In these experimental conditions, suspended load of kaolinite did not affect bedload rates of crushed brick particles. Three transport capacity formulae were tested against observed bedload rates. A calibration of the Foster formula revealed that the shear stress exponent should be greater than 1.5. The Low and the Govers unit stream power (USP) equations were then evaluated. The Low equation was preferred for the prediction of bedload rates of primary particles but it was not recommended in the case of aggregates of low density because of the limited experimental conditions applied to derive this equation. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

12.
Aeolian sand transport was studied at the Lanphere Dunes, a coastal dune complex in northern California, by comparing slipface advance rates with transport predicted based on local wind data. The slipfaces of a 2·5 m high transverse ridge and 10 m high parabolic dune were monitored over a period of three months to estimate sand discharge. The study was performed during the dry season, which has the maximum sand‐driving potential. Over the three month study period, average sand discharge was 12·5 m3 per m width per year at the transverse ridge and 8·8 m3 per m width per year at the parabolic dune. A method was developed for modelling slipfaces that are sinuous and where sediment transport rates are not constant across the width of the slipface. Field measurements were used to generate three‐dimensional representations of dune slipfaces. Periodic measurements over the course of three months were used to compute the volume of displaced sediment. Theoretical sand transport was computed from local wind data using the Bagnold model and compared with the observed transport rates. Predicted rates were substantially lower than observed rates. Wind velocities rarely exceeded the threshold velocity. Discrepancies between the observed and predicted values appear to be caused by a combination of wind data recording procedures and differences between wind velocities at the anemometer location and the site where sand transport was measured. Wind data collected by weather bureaux have been utilized in numerous studies for modelling sediment transport. Such data typically have sample intervals of one hour or greater and are often averaged prior to reporting. The effect of averaging was investigated by comparing sand transport estimates based on daily average wind velocities with those based on the original hourly observations. The daily average data were depleted of high velocity winds and sand transport estimates were accordingly much lower than those based on the hourly data. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

13.
River incision is fundamental in shaping the Earth's surface. In mountainous regions with steep river beds, fluvial bedrock erosion by bedload transport is an important mechanism forming channels. However, there are only a few complete field datasets that can be used to improve process understanding and evaluate erosion models, especially at the process scale. To provide a simultaneous dataset of hydraulics, bedload transport and bedrock erosion at high temporal and spatial resolution, a new measuring device has been installed in the Erlenbach, a gauged stream in the Swiss Pre‐Alps. In this stream, bedload transport rates can be calculated from surveying deposits and from geophone plate sensors and bedload transport samples can be taken directly by an automated moving basket system. To measure bedrock erosion rates simultaneously, two natural stone slabs were mounted flush with the channel bed in a steel frame hosting various measurement devices. Force sensors below the slabs record normal stress and shear stress. At‐a‐point erosion rates on the slab surfaces are continuously measured at sub‐millimetre precision at three locations on each slab. In addition, the slab topography is monitored following erosive flood events. In this article (i) the ‘erosion scale’ device is described, (ii) data resolution and data quality is assessed by means of tests and event data, and (iii) the first transport event is discussed. The erosion scales are confirmed to provide data at high spatio‐temporal resolution for process analysis. The preliminary data show evidence for the tools effect in bedrock erosion. The bedrock slabs can be exchanged to obtain measurements for catchments with different lithologies for comparison. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

14.
It is widely recognized nowadays that there are at least two different phases of bedload sediment transport in gravel‐bed rivers. However, the transition between these phases is still poorly or subjectively defined, especially at bends in rivers, where cross‐stream sediment transport can strongly influence changes in the texture of the transported sediment. In this paper, we use piecewise models to identify objectively, at two points in the cross‐section of a river bend, the discharge at which the transition between bedload transport phases occurs. Piecewise models were applied to a new bedload data set collected during a wide range of discharges while analysing the associated changes in sediment texture. Results allowed the identification of two well‐differentiated phases of sediment transport (phase I and phase II), with a breakpoint located around bankfull discharge. Associated with each phase there was a change in bedload texture. In phase I there was non‐dominance in the transport of fine or coarse fractions at a particular sampling point; but in phase II bedload texture was strongly linked to the position of the sampling point across the channel. In this phase, fine particles tended to be transported to the inner bank, while coarse sizes were transferred throughout the middle parts of the channel. Moreover, bedload texture at the inner sampling point became bimodal while the transport of pebble‐sized particles was increasing in the central parts of the river channel. It is suggested that this general pattern may be related both to secondary currents, which transfer finer particles from the outer to the inner bank, and to the progressive dismantling of the riverbed surface layer. Copyright © 2018 John Wiley & Sons, Ltd.  相似文献   

15.
Previous analyses have identified the active width of braided rivers, the bed area over which bed load flux and short‐term morphological change occurs, as an important element of braiding dynamics and predictions of bed load flux. Here we compare theoretical predictions of active width in gravel‐bed braided rivers with observations from Sunwapta River, and from a generic physical model of gravel braided rivers, to provide general observations of the variation in active width, and to develop an understanding of the causes of variation. Bed topography was surveyed daily along a 150 m reach of the pro‐glacial Sunwapta River for a total of four weeks during summer when flow was above threshold for morphological activity. In the laboratory, detailed digital elevation models (DEMs) were derived from photogrammetric survey at regular intervals during a constant discharge run. From the field and flume observations there is considerable local and circumstantial variation in active width, but also a general trend in average active width with increasing discharge. There is also a clear relationship of active width with active braiding index (number of active branches in the braided channel network), and with dimensionless stream power, which appears to be consistent across the range of data from field and physical models. Thus there is a link between active width and the river morphology and dynamics, and the possibility of a general relationship for estimating active width from channel pattern properties or reach‐scale stream power values, from which approximate bedload flux calculations may be made. The analysis also raises questions about differences between hydraulically‐based numerical model computations of instantaneous active width and observation of time‐integrated morphological active width. Understanding these differences can give insight into the nature of bedload transport in braided rivers and the relationship to morphological processes of braiding. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

16.
Indirect, passive approaches for monitoring coarse bedload transport could allow cheaper, safer, higher‐resolution, longer‐term data that revolutionises bedload understanding and informs river management. Here, insights provided by seismic impact plates in a downstream reach of a flashy gravel‐bed river (River Avon, Devon, UK) are explored in the context of plate performance. Monitoring of a centrally‐situated plate (IP1) during an extremely wet 12‐month period demonstrated that impacts were related to discharge as a measure of transport potential (R2 = 0.38) but that factors other than transport limitations are important. Analysis of discrete flow events revealed consistent rising‐limb and falling‐limb impact spikes biased toward the latter for larger events. Such patterns may result from disruption of the upstream armour layer (rising limb) and supply enhancements related to both upstream mass bank failures and/or flood routing of non‐local sediment sources (falling limb). Installation of additional impact plates indicated that plate IP1 was indeed dominantly related to instantaneous discharge, that a three‐plate lateral array somewhat better explained impact variability (R2 = 0.49), and that the bedload track shifts laterally with discharge. Aggregating event‐total IP1 impacts against volumetric discharge further increases explanation as intra‐event and stochastic bedload factors are subsumed but left 26% unexplained variance related to the unsampled bedload mass, inter‐event supply differences, and attributes of plate performance. Annualising the data created an impact‐based 'effective discharge’ for this extremely wet year that was closer to morphological bar‐full in magnitude than bankfull, but the preceding results imply this outcome is related as much to supply limitations as transport limitations. Overall, passive approaches offer a liberating prospect for bedload monitoring, capable of producing insights only achievable through high resolution, extended time periods. Such results could potentially inform threshold conditions and geomorphological effectiveness of flows for future river management strategies. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

17.
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.  相似文献   

18.
Understanding bedload transport fluctuations in rivers is crucial for complementing the existing knowledge on sediment transport theory. In this contribution, we use a natural-scale laboratory flume to analyse bedload transport fluctuations in non-uniform sand under normal flow conditions. Based on the significance of downward seepage, we incorporate the seepage effect on bedload transport over a non-uniform sand bed channel. The weight of the dry material was measured, and the volumetric transport rate per unit width (bedload transport rate) was estimated. An important observation is that the bedload transport rate initially rapidly increases with time and reaches a maximum value. Based on experimental data, we propose an empirical expression to estimate temporal bedload transport. In addition, an empirical model for bedload transport is proposed by incorporating downward seepage among other variables. The performance of several existing bedload transport formulae was also taken into account by the experimental datasets.  相似文献   

19.
ＵＮＣＥＲＴＡＩＮＴＹＡＮＤＳＥＮＳＩＴＩＶＩＴＹＡＮＡＬＹＳＥＳＯＦＳＥＤＩＭＥＮＴＴＲＡＮＳＰＯＲＴＦＯＲＭＵＬＡＳＫｅｈＣｈｉａＹＥＨ１ａｎｄＳｅｎＬｏｎｇＤＥＮＧ２ＡＢＳＴＲＡＣＴＩｎｖｉｅｗｏｆｔｈｅｒａｎｄｏｍｃｈａｒａｃｔｅｒｉｓｔ...  相似文献   

20.
Recalibrating aeolian sand transport models   总被引：1，自引：0，他引：1  
A quality‐controlled data set comprising measurements of aeolian sand transport rates obtained at three disparate field sites is used to evaluate six commonly employed transport rate models (those of Bagnold, Kawamura, Zingg, Owen, Hsu, and Lettau and Lettau) and to recalibrate the empirical constants in those models. Shear velocity estimates were obtained using the von Kármán constant and an apparent von Kármán parameter. Models were recalibrated using non‐linear regression and non‐linear regression with least‐squares lines forced through axes origins. Recalibration using the apparent von Kármán parameter and forced regression reduced the empirical constants for all models. The disparity between the predictions from the different models is reduced from about an order of magnitude to about a quarter of an order of magnitude. The recalibrated Lettau and Lettau model provided the greatest statistical agreement between observed and predicted transport rates, with a coefficient of determination of 0·77. Evaluation of the results suggests that our estimations of threshold shear velocity may be too slow, causing errors in predicted transport rates. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

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

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