首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 359 毫秒
1.
Salinity difference between terrestrial river discharge and oceanic tidal water plays a role in modifying the local flow field and, as a consequence, estuarine morphodynamics. Although widely recognized, recent numerical studies exploring the long-term morphological evolution of river-influenced estuaries with two-dimensional, depth-averaged models have mostly neglected salinity. Using a three-dimensional morphodynamic model, we aim to gain more insight into the effect of salinity on the morphodynamics of fluvio-deltaic systems. Model results indicate that the resultant estuarine morphology established after 600 years differs remarkably when a salinity gradient is included. A fan-shaped river-mouth delta exhibits less seaward expansion and is cut through by narrower channels when salinity is included. The inclusion of salinity tends to generate estuarine circulation, which favours landward sediment transport and hence limits the growth of the delta while enhancing the development of intertidal areas. The formation of deltaic channel–shoal patterns resulting from morphodynamic evolution tends to strengthen salinity stratification, which is characterized by an increased gradient Richardson number. The direction of the depth-averaged residual sediment transport over a tide may be opposite to the direction of residual velocity, indicating the significant influence of baroclinic effects on the net sediment transport direction (and hence morphological change). The effect of salinity on morphological evolution becomes less profound when the strength of tidal or fluvial forcing is dominant over the other. The effects of sediment type and flocculation, which are particularly important when salinity gradients are present, are also discussed. Overall, this study highlights that neglecting salinity to simulate long-term estuarine morphodynamics requires more careful justification, particularly when the environment is characterized by fine sediment types (favouring suspended transport), and relatively large river discharge and estuarine depth (favouring baroclinic effects). © 2020 John Wiley & Sons, Ltd.  相似文献   

2.
Erosion, sediment transportation and accumulation in rivers   总被引:8,自引:5,他引:3  
The present paper analyses the interrelation between erosion, sediment transportation and accumulation proposed by N. I. Makkaveyev (1908-1983) and its further development in modem studies of river channel processes in Russia. Spatio-temporal linkages between erosion and accumulation are defined considering channel processes at different scales - river longitudinal profile, channel morphological patterns, alluvial bedforms (bars, dunes) and individual sediment particles. Relations between river geomorphic activity, flow transportation capacity and sediment budgets are established (sediment input and output; channel bed erosion and sediment entrainment into flow - termination of sediment transport and its deposition). Channel planforms, floodplain segments separated by the latter and alluvial channel bedforms are shown to be geomorphic expressions of sediment transport process at different spatial and temporal scales. This paper is dedicated to the 100th anniversary of N. I. Makkaveyev, Professor of the Moscow State University, author of the book "River channel and erosion in its basin" (1955). That book is regarded in Russia as the pioneering work which initiated the complex hydrological and geographical studies of channel processes and laid a basis for the theory of unified fluvial erosion-accumulation process.  相似文献   

3.
This work is inspired by the sudden resurgence of the submersed aquatic vegetation (SAV) bed in the Chesapeake Bay (USA). Because the SAV bed occurs at the mouth of the Bay's main tributary (Susquehanna River), it plays a significant role in modulating sediment and nutrient inputs from the Susquehanna to the Bay. Previous model studies on the impact of submersed aquatic vegetation on the development of river mouth bars lacked a complete mechanistic understanding. This study takes advantage of new advances in 3D computational models that include explicit physical-sedimentological feedbacks to obtain this understanding. Specifically, we used Delft3D, a state-of-the-art hydrodynamic model that provides fine-scale computations of three-dimensional flow velocity and bed shear stress, which can be linked to sediment deposition and erosion. Vegetation is modeled using a parameterization of hydraulic roughness that depends on vegetation height, stem density, diameter, and drag coefficient. We evaluate the hydrodynamics, bed shear stresses, and sediment dynamics for different vegetation scenarios under conditions of low and high river discharge. Model runs vary the vegetation height, density, river discharge, and suspended-sediment concentration. Numerical results from the idealized model show that dense SAV on river mouth bars substantially diverts river discharge into adjacent channels and promotes sediment deposition at ridge margins, as well as upstream bar migration. Increasing vegetation height and density forms sandier bars closer to the river mouth and alteration of the bar shape. Thus, this study highlights the important role of SAV in shaping estuarine geomorphology, which is especially relevant for coastal management. © 2019 John Wiley & Sons, Ltd.  相似文献   

4.
With riverine flooding set to be more frequent in many parts of the world as a result of climate change, the interactions between fluvial morphodynamics and riparian vegetation may depend in part on the sequence of flood events. This paper describes a laboratory study of the geomorphic adjustment of a braided river to sequences of floods across five different strengths of braidplain vegetation. By using alfalfa as a proxy for braidplain vegetation, the differing plant life stages were used to represent the varying strengths of biogeomorphic feedbacks across the floods. Boundary conditions were constrained by sets of experimental runs with both equilibrium sediment loads and deficit loads. Changes in bed topography were monitored and assessed using a detailed digital elevation model, digital imagery and continuous monitoring of the transported sediment. Results demonstrate that in absence of plant colonization, vegetation placed the rivers in a non-equilibrium condition, in which riparian vegetation encouraged the development of new channels, increased the system channel width and enhanced topographic irregularity, these effects being more noticeable during the low-flow periods. The morphodynamics was found to be less sensitive to variations in flood discharges as the vegetation influence (strength) increased from minimum to maximum, until vegetation began to die back and the impacts of flood sequences became yet again evident. Although the overall sediment transport rate was reduced under full-grown vegetation conditions, the presence of the mature plants across the braid bars resulted in the greatest channel scour depths. Results are considered in light of expected changes in flood frequency with climate and likely morphodynamic responses of river systems as a result.  相似文献   

5.
Confluences with low discharge and momentum ratios, where narrow steep tributaries with high sediment load join a wide low‐gradient main channel that provides the main discharge, are often observed in high mountain regions such as in the upper‐Rhone river catchment in Switzerland. Few existing studies have examined the hydro‐morphodynamics of this type of river confluence while considering sediment discharge in both confluent channels. This paper presents the evolution of the bed morphology and hydrodynamics as observed in an experimental facility with a movable bed. For that purpose, one experiment was carried out in a laboratory confluence with low discharge and momentum ratios, where constant sediment rates were supplied to both flumes. During the experiment, bed topography and water surface elevations were systematically recorded. When the bed topography reached a steady state (so‐called equilibrium) and the outgoing sediment rate approximated the incoming rate, flow velocity was measured at 12 different points distributed throughout the confluence, and the grain size distribution of the bed surface was analyzed. Typical morphodynamic features of discordant confluences such as a bank‐attached bar and a flow deflection zone are identified in this study. Nevertheless, the presence of a marked scour hole in the discordant confluence and distinct flow regimes for the tributary and main channel, differ from results obtained in previous studies. Strong acceleration of the flow along the outer bank of the main channel is responsible for the scour hole. This erosion is facilitated by the sediment discharge into the confluence from the main channel which inhibits bed armoring in this region. The supercritical flow regime observed in the tributary is the hydrodynamic response to the imposed sediment rate in the tributary. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

6.
An in‐house fully three‐dimensional general‐purpose finite element model is applied to solve the hydrodynamic structure in a periodic Kinoshita‐generated meandering channel. The numerical model solves the incompressible Reynolds‐averaged Navier–Stokes equations for mass and momentum, while solving the k ? ε equations for turbulence. The free surface is described by the rigid‐lid approximation (using measured water surface data) for flat (smooth‐bed) and self‐formed (rough‐bed) conditions. The model results are compared against experimental measurements in the ‘Kinoshita channel’, where three‐dimensional flow velocities and turbulence parameters were measured. This validation was carried out for the upstream‐valley meander bend orientation under smooth (flat bed) conditions. After validation, several simulations were carried out to predict the hydrodynamics in conditions where either it was not possible to perform measurements (e.g. applicability of the laboratory acoustic instruments) and to extrapolate the model to other planform configurations. For the flat smooth‐bed case, a symmetric (no skewness) planform configuration was modeled and compared to the upstream‐skewed case. For the self‐formed rough‐bed case, prediction of the hydrodynamics during the progression of bedforms was performed. It appears that the presence of bedforms on a bend has the following effects: (i) the natural secondary flow of the bend is disrupted by the presence of the bedforms, thus depending on the location of the dune, secondary flows might differ completely from the traditional orientation; (ii) an increment on both the bed and bank shear stresses is induced, having as much as 50% more fluvial erosion, and thus a potential increment on the migration rate of the bend. Implications on sediment transport and bend morphodynamics are also discussed in the paper. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

7.
Sediment often enters rivers in the form of sediment pulses associated with landslides and debris flows. This is particularly so in gravel‐bed rivers in earthquake‐prone mountain regions, such as Southwest China. Under such circumstances, sediment pulses can rapidly change river topography and leave the river in repeated states of gradual recovery. In this paper, we implement a one‐dimensional morphodynamic model of river response to pulsed sediment supply. The model is validated using data from flume experiments, so demonstrating that it can successfully reproduce the overall morphodynamics of experimental pulses. The model is then used to explore the evolution of a gravel‐bed river subject to cycled hydrographs and repeated sediment pulses. These pulses are fed into the channel in a fixed region centered at a point halfway down the calculational domain. The pulsed sediment supply is in addition to a constant sediment supply at the upstream end. Results indicate that the river can reach a mobile‐bed equilibrium in which two regions exist within which bed elevation and surface grain size distribution vary periodically in time. One of these is at the upstream end, where a periodic discharge hydrograph and constant sediment supply are imposed, and the other is in a region about halfway down the channel where periodic sediment pulses are introduced. Outside these two regions, bed elevation and surface grain size distribution reach a mobile‐bed equilibrium that is invariant in time. The zone of fluctuation‐free mobile‐bed equilibrium upstream of the pulse region is not affected by repeated sediment pulses under the scenarios tested, but downstream of the pulse region, the channel reaches different fluctuation‐free mobile‐bed equilibriums under different sediment pulse scenarios. The vertical bed structure predicted by the simulations indicates that the cyclic variation associated with the hydrograph and sediment pulses can affect the substrate stratigraphy to some depth. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

8.
Morphodynamics in sand‐bed braided rivers are associated with simultaneous evolution of mid‐channel bars and channels on the braidplain. Bifurcations around mid‐channel bars are key elements that divide discharge and sediment. This, in turn, may control the evolution of connected branches, with effects propagating to both upstream and downstream bifurcations. Recent works on bifurcation stability and development hypothesize major roles of secondary flow and gradient advantage. However, this has not been tested for channel networks within a fully developed dynamic braided river. A reason for this is a lack of detailed measurements with sufficient temporal and spatial length, covering multiple bifurcations. Therefore we used a physics‐based numerical model to generate a dataset of bathymetry, flow and sediment transport of an 80 km river reach with self‐formed braid bars and bifurcations. The study shows that bar dissection due to local transverse water surface gradients is the dominant bifurcation initiation mechanism, although conversion of unit bars into compound bars dominates in the initial stage of a braided river. Several bifurcation closure mechanisms are equally important. Furthermore, the study showed that nodal point relations for bifurcations are unable to predict short‐term bifurcation evolution in a braided river. This is explained by occurrence of nonlinear processes and non‐uniformity within the branches, in particular migrating bars and larger‐scale backwater‐effects, which are not included in the nodal point relations. Planform morphology, on the other hand, has predictive capacity: bifurcation angle asymmetry and bar‐tail limb shape are indicators for near‐future bifurcation evolution. Remote sensing data has predictive value, for which we developed a conceptual model for interactions between bars, bifurcations and channels in the network. We conducted a preliminary test of the conceptual model on satellite images of the Brahmaputra. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

9.
This work deals with the impacts of dams on large gravel -bed rivers in terms of altering coarse transport regimes and the relationship with river morphodynamics. Using data collected by a tracer -based monitoring programme carried out in a 4 -km -long study sector of the Parma River (Italy), located downstream from a relatively recently established dam, we applied a virtual velocity approach to estimate the coarse bed material load at four river cross -sections. Monitoring and calculation results provided new insights into the impacts of the dam on streambed material mobility and the sediment regime over the 17 -month calculation period. A longitudinal gradient of effects was observed along the study sector. Sections located closer to the dam are characterized by more evident impacts due to deficits in coarse sediment input from upstream. Sediment mobility here is strongly altered, especially in the highly armoured main channel, and the overall bed material load is extremely low. A partial recovery of sediment dynamics was observed at the sections located further from the dam, where estimates indicate higher sediment yield. The observed longitudinal trend in the coarse sediment transport regime matches the morphology, as the river shifts downstream from a sinuous configuration with alternate bars to a wandering one. The novel insights into alteration of coarse sediment dynamics and the relationship with river morphodynamics are potentially applicable to many other fluvial contexts affected by similar impoundments. © 2019 John Wiley & Sons, Ltd.  相似文献   

10.
Antidunes are fluvial bedforms that form in rivers with supercritical flows. The water surface over antidunes is strongly in phase with the bed surface, and the water surface is amplified to produce large surface waves. Many experimental studies have addressed antidunes; however, the shapes of three-dimensional antidunes in a wide channel with alternate bars have not yet been appropriately understood. In this study, we experimentally investigated the streamwise and transverse length scales of antidunes under conditions with a large width–depth ratio. Our experimental results provide evidence for the coevolution of antidunes and free alternate bars, and show for the first time that the development of free bars greatly alters the three-dimensional shape of water surface waves over antidunes. In the absence of free bars in a wide channel, multiple longitudinal wave trains form, and the number of wave trains counted in the transverse direction increases with increases in the width–depth ratio. However, the presence of free bars affects the local flow characteristics, resulting in a decrease of the number of wave trains in the transverse direction. Therefore, we propose a simple model for predicting the reduction in the number of wave trains by combining two previous theories for antidunes and free bars. Results obtained by the model were found to largely agree with experimental observations. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd  相似文献   

11.
Occurrence and development of channel bars are major components of the morphodynamics of rivers and their relation to river meandering has been much explored through theory and experimentation. However, field and documentary data of characteristics and evolution over timescales from years to several decades are lacking. Four sets of aerial photographs in the period 1984–2007 were used to map and quantify bar numbers and areas in GIS on an active meandering reach. Bar types were classified. Additional temporal resolution was provided by annual ground photography and mapping for 1981–2010. Analysis was extended backward by use of large scale Ordnance Survey maps from 1873 onwards. As expected, point bars are the most common type but ‘free’ bars of several types are major components of bar deposition. Point bars and attached bars are significantly larger in size than mid‐channel and side bars. Spatial distribution of bars varies down the reach and over time but is related to channel sinuosity, gradient and mobility and to bend evolution. Different types of bar occur in distinctive channel locations, with point and concave‐bend bars in zones of high curvature. Bar activity shows a relation with discharge events and phases and possibly with changing riparian conditions, but superimposed on this is a common sequence of bar evolution from incipient gravel mid‐channel bars to full floodplain integration. This life‐cycle is identified as 7–9 years on average. No evidence for mobility of free bars within the course is found. The results are compared with bar and bend theory; the bars are forced and conform in general to bend theory but detailed variation relates to geomorphic factors and to autogenic sequences of bends and bars. Mid‐channel bars are width induced. Variability of bar occurrence needs to be taken into account in river management and ecological evaluation, including for the EU WFD. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

12.
The development of alternate bars in channelized rivers can be explained theoretically as an instability of the riverbed when the active channel width to depth ratio exceeds a threshold. However, the development of a vegetation cover on the alternate bars of some channelized rivers and its interactions with bar morphology have not been investigated in detail. Our study focused on the co‐evolution of alternate bars and vegetation along a 33 km reach of the Isère River, France. We analysed historical information to investigate the development of alternate bars and their colonization by vegetation within a straightened, embanked river subject to flow regulation, sediment mining, and vegetation management. Over an 80 year period, bar density decreased, bar length increased, and bar mobility slowed. Vegetation encroachment across bar surfaces accompanied these temporal changes and, once established, vegetation cover persisted, shifting the overall system from an unvegetated to a vegetated dynamic equilibrium state. The unvegetated morphodynamics of the impressively regular sequence of alternate bars that developed in the Isère following channelization is consistent with previous theoretical morphodynamic work. However, the apparent triggering dynamics of vegetation colonization needs to be investigated, based on complex biophysical instability processes. If instability related to vegetation colonization is confirmed, further work needs to focus on the relevance of initial conditions for this instability, and on related feedback effects such as how the morphodynamics of bare‐sediment alternate bars may have affected vegetation development and, in turn, how vegetation has created a new dynamic equilibrium state. Copyright © 2018 John Wiley & Sons, Ltd.  相似文献   

13.
1INTRODUCTIONHumanactivitiesmayinducesignificantchangesinsedimenttransportcapacityofaluvialriverchannelsandcausechangesinthe...  相似文献   

14.
As economic development upstream in the Yangtze River basin has progressed in recent decades,the demand for sediment has rapidly increased and contributed to an expansion in sediment excavation that may affect the river’s stability and navigation safety.In the current study,the distribution of gravel mining in the upstream reach of the Yangtze River was investigated using field measurements obtained from2008 to 2017.An experimental investigation was then done to analyze the bed load behavior in ...  相似文献   

15.
In the middle and lower reaches of alluvial rivers, various kinds of river regulation projects affecting natural channel evolution often are distributed due to the requirements of flood control, navigation,and channel stability. However, the influence of large-scale river regulation works on fluvial processes is not fully known. Therefore, a two-dimensional(2D) morphodynamic model has been improved to address this problem. The new detailed procedure is presented in this paper:(i) First, each nod...  相似文献   

16.
Channel bars and banks strongly affect the morphology of both braided and meandering rivers. Accordingly, bar formation and bank erosion processes have been greatly explored. There is, however, a lack of investigations addressing the interactions between bed and bank morphodynamics, especially over short timescales. One major implication of this gap is that the processes leading to the repeated accretion of mid‐channel bars and associated widenings remain unsolved. In a restored section of the Drau River, a gravel‐bed river in Austria, mid‐channel bars have developed in a widening channel. During mean flow conditions, the bars divert the flow towards the banks. One channel section exhibited both an actively retreating bank and an expanding mid‐channel bar, and was selected to investigate the morphodynamic processes involved in bar accretion and channel widening at the intra‐event timescale. We repeatedly surveyed riverbed and riverbank topography, monitored riverbank hydrology and mounted a time‐lapse camera for continuous observation of riverbank erosion processes during four flow events. The mid‐channel bar was shown to accrete when it was submerged during flood events, which at the subsequent flow diversion during lower discharges narrowed the branch along the bank and increased the water surface elevation upstream from the riffle, which constituted the inlet into the branch. These changes of bed topography accelerated the flow along the bank and triggered bank failures up to 20 days after the flood events. Four analysed flow events exhibited a total bar expansion from initially 126 m2 to 295 m2, while bank retreat was 6 m at the apex of the branch. The results revealed the forcing role of bar accretion in channel widening and highlighted the importance of intra‐event scale bed morphodynamics for bank erosion, which were summarized in a conceptual model of the observed bar–bank interactions. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

17.
If increased sediment supply to a river channel exceeds its transport capacity, deposition necessarily occurs as the bed adjusts to accommodate the increased supply. Both the mean and spatial patterns in bed elevation and grain size may change and an ability to understand their relative importance is needed to predict bed response. We report on an experiment in a field‐scale flume in which sediment supply is increased to a gravel bed with alternate bars. Sediment was recirculated in the experiments, but augmented in two steps, after which the bed was allowed to reach a new steady state. The transport rate at the end of the experiment was three times larger than at the start. High‐resolution sediment flux and topographic measurements, grain size derived from photographs, and hydrodynamic modeling allow us to document the topographic and textural response of the bed to increased sediment supply. The spatial patterns of bed topography and texture were forced by the flume setup and the initial and final steady states included long stationary alternate bars with associated grain size sorting. The transient bed contained several scales of shorter wavelength migrating bedforms superimposed on, and temporarily replacing the stationary alternate bars. Bed topography and textural patterns adjusted to increased sediment supply over different timescales. Bed slope and mean stress increased directly with sediment supply rate to produce a new transport steady state in a time about 2.5 times the minimum needed to deposit the required sediment wedge, indicating a trap efficiency of about 40% for the aggrading wedge. Adjustments in local topography and sorting, primarily in the form of smaller, migrating bars, continued for a period approximately equal to that required to initially reach transport steady state. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

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

19.
The paper describes the derivation of a depth-averaged, two-dimensional form of the sediment balance equation, suitable to study the morphodynamics of movable sediment beds even when the flow depth attains values comparable to bed irregularities. This equation is derived by double-averaging in time and in space the instantaneous three-dimensional sediment balance equation. For this, a proper phase function is introduced, which depends on the statistics of bed topography. The structure of the macroscopic volumetric sediment discharge vector resulting from the averaging procedure is discussed for the case of dominant bedload transport. The theoretical framework developed within the paper sets the stage for a proper parametrization of the physical processes acting at spatial scales smaller than those usually resolved by depthaveraged numerical models.  相似文献   

20.
Laboratory observations and computational results for the response of bedform fields to rapid variations in discharge are compared and discussed. The simple case considered here begins with a relatively low discharge over a flat bed on which bedforms are initiated, followed by a short high‐flow period with double the original discharge, during which the morphology of the bedforms adjusts, followed in turn by a relatively long period of the original low discharge. For the grain size and hydraulic conditions selected, the Froude number remains subcritical during the experiment, and sediment moves predominantly as bedload. Observations show rapid development of quasi‐two‐dimensional bedforms during the initial period of low flow with increasing wavelength and height over the initial low‐flow period. When the flow increases, the bedforms rapidly increase in wavelength and height, as expected from other empirical results. When the flow decreases back to the original discharge, the height of the bedforms quickly decreases in response, but the wavelength decreases much more slowly. Computational results using an unsteady two‐dimensional flow model coupled to a disequilibrium bedload transport model for the same conditions simulate the formation and initial growth of the bedforms fairly accurately and also predict an increase in dimensions during the high‐flow period. However, the computational model predicts a much slower rate of wavelength increase, and also performs less accurately during the final low‐flow period, where the wavelength remains essentially constant, rather than decreasing. In addition, the numerical results show less variability in bedform wavelength and height than the measured values; the bedform shape is also somewhat different. Based on observations, these discrepancies may result from the simplified model for sediment particle step lengths used in the computational approach. Experiments show that the particle step length varies spatially and temporally over the bedforms during the evolution process. Assuming a constant value for the step length neglects the role of flow alterations in the bedload sediment‐transport process, which appears to result in predicted bedform wavelength changes smaller than those observed. However, observations also suggest that three‐dimensional effects play at least some role in the decrease of bedform wavelength, so incorporating better models for particle hop lengths alone may not be sufficient to improve model predictions. Published in 2011. This article is a US Government work and is in the public domain in the USA.  相似文献   

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

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