首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The sediment saturation recovery process (i.e. the adaptation of suspended sediment concentration [SSC] to local forcing) is the main feature of the non‐equilibrium suspended sediment transport (SST) frequently occurring in fluvial, estuarine and coastal waters. In order to quantitatively describe this phenomenon, a series solution is analytically derived, including the evolution of both vertical SSC profile and near‐bed sediment flux (NBSF), and is verified by net erosion and net deposition experiments, respectively. The results suggest that the sediment saturation recovery process involves vertically varying fluxes that are not represented correctly by depth‐averaging. Consequently, a vertical two‐dimensional (2D) combined scheme is established and applied respectively in to a dredged trench and to a sand wave feature to demonstrate this argument. By analyzing the variations of the calculated depth‐averaged SSC and NBSF we reveal that the equilibrium state presented by the sediment carrying capacity (SCC) form of the NBSF, which is usually applied in depth‐integrated SST models, lags behind the actual dynamic bed equilibrium state. Moreover, the key factor α, the so‐called saturation recovery coefficient within this form, is not only a function of local Rouse number but also is influenced by the local SSC profile. Finally, a three‐dimensional (3D) non‐orthogonal curvilinear body‐fitted SST model is developed and validated in the Yangtze estuary, China, combined with the in situ hourly hydrographic data from August 14–15, 2007 during spring tide in the wet season. Model results confirm that the vertically varying sediment saturation recovery process, the discrepancies between the actual and SCC form of NBSF and non‐constant value of α are significant in actual real geomorphic cases. The quantitative morphological change resulting from variations in environmental conditions may not be correctly represented by uncorrected depth‐integrated SST models if they do not treat the effects of vertical motion on the sediment saturation recovery process. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

2.
AVO investigations of shallow marine sediments   总被引：2，自引：0，他引：2  
Amplitude‐variation‐with‐offset (AVO) analysis is based on the Zoeppritz equations, which enable the computation of reflection and transmission coefficients as a function of offset or angle of incidence. High‐frequency (up to 700 Hz) AVO studies, presented here, have been used to determine the physical properties of sediments in a shallow marine environment (20 m water depth). The properties that can be constrained are P‐ and S‐wave velocities, bulk density and acoustic attenuation. The use of higher frequencies requires special analysis including careful geometry and source and receiver directivity corrections. In the past, marine sediments have been modelled as elastic materials. However, viscoelastic models which include absorption are more realistic. At angles of incidence greater than 40°, AVO functions derived from viscoelastic models differ from those with purely elastic properties in the absence of a critical angle of incidence. The influence of S‐wave velocity on the reflection coefficient is small (especially for low S‐wave velocities encountered at the sea‐floor). Thus, it is difficult to extract the S‐wave parameter from AVO trends. On the other hand, P‐wave velocity and density show a considerably stronger effect. Attenuation (described by the quality factor Q) influences the reflection coefficient but could not be determined uniquely from the AVO functions. In order to measure the reflection coefficient in a seismogram, the amplitudes of the direct wave and the sea‐floor reflection in a common‐midpoint (CMP) gather are determined and corrected for spherical divergence as well as source and streamer directivity. At CMP locations showing the different AVO characteristics of a mud and a boulder clay, the sediment physical properties are determined by using a sequential‐quadratic‐programming (SQP) inversion technique. The inverted sediment physical properties for the mud are: P‐wave velocity α=1450±25 m/s, S‐wave velocity β=90±35 m/s, density ρ=1220±45 kg/m3, quality factor for P‐wave QP=15±200, quality factor for S‐wave QS=10±30. The inverted sediment physical properties for the boulder clay are: α=1620±45 m/s,β=360±200 m/s,ρ=1380±85 kg/m3,QP=790±660,QS=25±10.  相似文献   

3.
The important effects of bottom sediments on the seismic response of arch dams are studied in this paper. To do so, a three‐dimensional boundary element model is used. It includes the water reservoir as a compressible fluid, the dam and unbounded foundation rock as viscoelastic solids, and the bottom sediment as a two‐phase poroelastic domain with dynamic behaviour described by Biot's equations. Dynamic interaction among all those regions, local topography and travelling wave effects are taken into account. The results obtained show the important influence of sediment compressibility and permeability on the seismic response. The former is associated with a general change of the system response whereas the permeability has a significant influence on damping at resonance peaks. The analysis is carried out in the frequency domain considering time harmonic excitation due to P and S plane waves. The time‐domain results obtained by using the Fourier transform for a given earthquake accelerogram are also shown. The possibility of using simplified models to represent the bottom sediment effects is discussed in the paper. Two alternative models for porous sediment are tested. Simplified models are shown to be able to reproduce the effects of porous sediments except for very high permeability values. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

4.
Mountain ranges are frequently subjected to mass wasting events triggered by storms or earthquakes and supply large volumes of sediment into river networks. Besides altering river dynamics, large sediment deliveries to alluvial fans are known to cause hydro‐sedimentary hazards such as flooding and river avulsion. Here we explore how the sediment supply history affects hydro‐sedimentary river and fan hazards, and how well can it be predicted given the uncertainties on boundary conditions. We use the 2D morphodynamic model Eros with a new 2D hydrodynamic model driven by a sequence of flood, a sediment entrainment/transport/deposition model and a bank erosion law. We first evaluate the model against a natural case: the 1999 Mount Adams rock avalanche and subsequent avulsion on the Poerua river fan (West Coast, New Zealand). By adjusting for the unknown sediment supply history, Eros predicts the evolution of the alluvial riverbed during the first post‐landslide stages within 30 cm. The model is subsequently used to infer how the sediment supply volume and rate control the fan aggradation patterns and associated hazards. Our results show that the total injected volume controls the overall levels of aggradation, but supply rates have a major control on the location of preferential deposition, avulsion and increased flooding risk. Fan re‐incision following exhaustion of the landslide‐derived sediment supply leads to sediment transfer and deposition downstream and poses similar, but delayed, hydro‐sedimentary hazards. Our results demonstrate that 2D morphodynamics models are able to capture the full range of hazards occurring in alluvial fans including river avulsion aggradation and floods. However, only ensemble simulations accounting for uncertainties in boundary conditions (e.g., discharge history, initial topography, grain size) as well as model realization (e.g., non‐linearities in hydro‐sedimentary processes) can be used to produce probabilistic hazards maps relevant for decision making. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

5.
Numerical analyses of liquefiable sand are presented in this paper. Liquefaction phenomenon is an undrained response of saturated sandy soils when they are subjected to static or dynamic loads. A fully coupled dynamic computer code is developed to predict the liquefaction potential of a saturated sandy layer. Coupled dynamic field equations of extended Biot's theory with u–P formulation are used to determine the responses of pore fluid and soil skeleton. Generalized Newmark method is employed for integration in time. The soil behavior is modelled by two constitutive models; a critical state two-surface plasticity model, and a densification model. A class ‘B’ analysis of a centrifuge experiment is performed to simulate the dynamic response of level ground sites. The results of the numerical analyses demonstrate the capability of the critical sate two-surface plasticity model in producing pore pressures that are consistent with observations of the behavior of liquefiable sand in the centrifuge test.  相似文献   

6.
  相似文献   

7.
Distributed watershed models are beneficial tools for the assessment of management practices on runoff and water‐induced erosion. This paper evaluates, by application to an experimental watershed, two promising distributed watershed‐scale sediment models in detail: the Kinematic Runoff and Erosion (KINEROS‐2) model and the Gridded Surface Subsurface Hydrologic Analysis (GSSHA) model. The physics behind each model are to some extent similar, though they have different watershed conceptualizations. KINEROS‐2 was calibrated using three rainfall events and validated over four separate rainfall events. Parameters estimated by this calibration process were adapted to GSSHA. With these parameters, GSSHA generated larger and retarded flow hydrographs. A 30% reduction in both plane and channel roughness in GSSHA along with the assumption of Green‐Ampt conductivity KG‐A = Ks, where Ks is the saturated conductivity, resulted in almost identical hydrographs. Sediment parameters not common in both models were calibrated independently of KINEROS‐2. A comparative discussion of simulation results is presented. Even though GSSHA's flow component slightly overperformed KINEROS‐2, the latter outperformed GSSHA in simulations for sediment transport. In spite of the fact that KINEROS‐2 is not geared toward continuous‐time simulations, simulations performed with both models over a 1 month period generated comparable results. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

8.
Recognizing that soil–structure interaction affects appreciably the earthquake response of highway overcrossings, this paper compares approximate analytical solutions and finite element results to conclude on a simple procedure that allows for the estimation of the kinematic response functions and dynamic stiffnesses of approach embankments. It is shown that the shear‐wedge model yields realistic estimates for the amplification functions of typical embankments and reveals the appropriate levels of dynamic strains which are subsequently used to estimate the stiffness and damping coefficients of embankments. The shear‐wedge model is extended to a two‐dimensional model in order to calculate the transverse static stiffness of an approach embankment loaded at one end. The formulation leads to a sound closed‐form expression for the critical length, Lc, that is the ratio of the transverse static stiffness of an approach embankment and the transverse static stiffness of a unit‐width wedge. It is shown through two case studies that the transverse dynamic stiffness (‘spring’ and ‘dashpot’) of the approach embankment can be estimated with confidence by multiplying the dynamic stiffness of the unit‐width wedge with the critical length, Lc. The paper concludes that the values obtained for the transverse kinematic response function and dynamic stiffness can also be used with confidence to represent the longitudinal kinematic response function and dynamic stiffness, respectively. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

9.
Glaciers are major agents of erosion that increase sediment load to the downstream fluvial system. The Castle Creek Glacier, British Columbia, Canada, has retreated ~1.0 km in the past 70 years. Suspended sediment concentration (SSC) and streamflow (Q) were monitored independently at five sites within its pro‐glacial zone over a 60 day period from July to September 2011, representing part of the ablation season. Meteorological data were collected from two automatic weather stations proximal to the glacier. The time‐series were divided into hydrologic days and the shape and magnitude of the SSC response to hydro‐meteorological conditions (‘cold and wet’, ‘hot and dry’, ‘warm and damp’, and ‘storm’) were categorized using principal component analysis (PCA) and cluster analysis (CA). Suspended sediment load (SSL) was computed and summarized for the categories. The distribution of monitoring sites and results of the multivariate statistical analyses describe the temporal and spatial variability of suspended sediment flux and the relative importance of glacial and para‐glacial sediment sources in the pro‐glacial zone. During the 2011 study period, ~ 60% of the total SSL was derived from the glacial stream and sediment deposits proximal to the terminus of the glacier; during ‘storm’ events, that contribution dropped to ~40% as the contribution from diffuse and point sources of sediment throughout the pro‐glacial zone and within the meltwater channels increased. While ‘storm’ events accounted for just 3% of the study period, SSL was ~600% higher than the average over the monitoring period, and ~20% of the total SSL was generated in that time. Determining how hydro‐meteorological conditions and sediment sources control sediment fluxes will assist attempts to predict how pro‐glacial zones respond to future climate changes. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

10.
Predicting sediment yield at the catchment scale is one of the main challenges in geomorphologic research. The application of both physics‐based models and regression models has until now not provided very satisfying results for prediction of sediment yield for medium to large sized catchments (c. >50 km2). The explanation for this lies in a combination of the large data requirements of most models and a lack of knowledge to describe all processes and process interactions at the catchment scale. In particular, point sources of sediment (e.g. gullies, mass movements), connectivity and sediment transport remain difficult to describe in most models. From reservoir sedimentation data of 44 Italian catchments, it appeared that there was a (non‐significant) positive relation between catchment area and sediment yield. This is in contrast to what is generally expected from the theory of decreasing sediment delivery rates with increasing catchment area. Furthermore, this positive relation suggests that processes other than upland erosion are responsible for catchment sediment yield. Here we explore the potential of the Factorial Scoring Model (FSM) and the Pacific Southwest Interagency Committee (PSIAC) model to predict sediment yield, and indicate the most important sediment sources. In these models different factors are used to characterize a drainage basin in terms of sensitivity to erosion and connectivity. In both models an index is calculated that is related to sediment yield. The FSM explained between 36 and 61 per cent of the variation in sediment yield, and the PSIAC model between 57 and 62 per cent, depending on the factors used to characterize the catchments. The FSM model performed best based on a factor to describe gullies, lithology, landslides, catchment shape and vegetation. Topography and catchment area did not explain additional variance. In particular, the addition of the landslide factor resulted in a significantly increased model performance. The FSM and PSIAC model both performed better than a spatially distributed model describing water erosion and sediment transport, which was applied to the same catchments but explained only between 20 and 51 per cent of the variation in sediment yield. Model results confirmed the hypothesis that processes other than upland erosion are probably responsible for sediment yield in the Italian catchments. A promising future development of the models is by the use of detailed spatially distributed data to determine the scores, decrease model subjectivity and provide spatially distributed output. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

11.
The implicit dissipative generalized‐ α method is analyzed using discrete control theory. Based on this analysis, a one‐parameter family of explicit direct integration algorithms with controllable numerical energy dissipation, referred to as the explicit KR‐α method, is developed for linear and nonlinear structural dynamic numerical analysis applications. Stability, numerical dispersion, and energy dissipation characteristics of the proposed algorithms are studied. It is shown that the algorithms are unconditionally stable for linear elastic and stiffness softening‐type nonlinear systems, where the latter indicates a reduction in post yield stiffness in the force–deformation response. The amount of numerical damping is controlled by a single parameter, which provides a measure of the numerical energy dissipation at higher frequencies. Thus, for a specific value of this parameter, the resulting algorithm is shown to produce no numerical energy dissipation. Furthermore, it is shown that the influence of the numerical damping on the lower mode response is negligible. It is further shown that the numerical dispersion and energy dissipation characteristics of the proposed explicit algorithms are the same as that of the implicit generalized‐ α method. A numerical example is presented to demonstrate the potential of the proposed algorithms in reducing participation of undesired higher modes by using numerical energy dissipation to damp out these modes. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

12.
A 2D depth‐averaged hydrodynamic, sediment transport and bed morphology model named STREMR HySeD is presented. The depth‐averaged sediment transport equations are derived from the 3D dilute, multiphase, flow equations and are incorporated into the hydrodynamic model STREMR. The hydrodynamic model includes a two‐equation turbulence model and a correction for the mean flow due to secondary flows. The suspended sediment load can be subdivided into different size classes using the continuum (two‐fluid) approach; however, only one bed sediment size is used herein. The validation of the model is presented by comparing the suspended sediment transport module against experimental measurements and analytical solutions for the case of equilibrium sediment‐laden in a transition from a rigid bed to a porous bed where re‐suspension of sediment is prevented. On the other hand, the bed‐load sediment transport and bed evolution numerical results are compared against bed equilibrium experimental results for the case of a meander bend. A sensitivity analysis based on the correction for secondary flow on the mean flow including the effect of secondary flow on bed shear stresses direction as well as the downward acceleration effect due to gravity on transverse bed slopes is performed and discussed. In general, acceptable agreement is found when comparing the numerical results obtained with STREMR HySeD against experimental measurements and analytical solutions. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

13.
The stream power incision model (SPIM) is a cornerstone of quantitative geomorphology. It states that river incision rate is the product of drainage area and channel slope raised to the power exponents m and n, respectively. It is widely used to predict patterns of deformation from channel long profile inversion or to model knickpoint migration and landscape evolution. Numerous studies have attempted to test its applicability with mixed results prompting the question of its validity. This paper synthesizes these results, highlights the SPIM deficiencies, and offers new insights into the role of incision thresholds and channel width. By reviewing quantitative data on incising rivers, I first propose six sets of field evidence that any long‐term incision model should be able to predict. This analysis highlights several inconsistencies of the standard SPIM. Next, I discuss the methods used to construct physics‐based long‐term incision laws. I demonstrate that all published incising river datasets away from knickpoints or knickzones are in a regime dominated by threshold effects requiring an explicit upscaling of flood stochasticity neglected in the standard SPIM and other incision models. Using threshold‐stochastic simulations with dynamic width, I document the existence of composite transient dynamics where knickpoint propagation locally obeys a linear SPIM (n=1) while other part of the river obey a non‐linear SPIM (n>1). The threshold‐stochastic SPIM resolves some inconsistencies of the standard SPIM and matches steady‐state field evidence when width is not sensitive to incision rate. However it fails to predict the scaling of slope with incision rate for cases where width decreases with incision rate. Recent proposed models of dynamic width cannot resolve these deficiencies. An explicit upscaling of sediment flux and threshold‐stochastic effects combined with dynamic width should take us beyond the SPIM which is shown here to have a narrow range of validity. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
It seems that the explicit KR‐α method (KRM), which was developed by Kolay and Ricles, is promising for the step‐by‐step integration because it simultaneously integrates unconditional stability, explicit formulation, and numerical dissipation together. It was shown that KRM can inherit the numerical dispersion and energy dissipation properties of the generalized‐α method [1] for a linear elastic system, and it reduces to CR method (CRM), which was developed by Chen and Ricles [2] if ρ∞ = 1 is adopted, where ρ∞ is the spectral radius of the amplification matrix of KRM as the product of the natural frequency and the step size tends to infinity. However, two unusual properties were found for KRM and CRM, and they might limit their application to solve either linear elastic or nonlinear systems. One is the lack of capability to capture the structural nonlinearity, and the other is that it is unable to realistically reflect the dynamic loading. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

15.
Flow diversion terraces (FDT) are commonly used beneficial management practice (BMP) for soil conservation on sloped terrain susceptible to water erosion. A simple GIS‐based soil erosion model was designed to assess the effectiveness of the FDT system under different climatic, topographic, and soil conditions at a sub‐basin level. The model was used to estimate the soil conservation support practice factor (P‐factor), which inherently considered two major outcomes with its implementation, namely (1) reduced slope length, and (2) sediment deposition in terraced channels. A benchmark site, the agriculture‐dominated watershed in northwestern New Brunswick (NB), was selected to test the performance of the model and estimated P‐factors. The estimated P‐factors ranged from 0·38–1·0 for soil conservation planning objectives and ranged from 0·001 to 0·45 in sediment yield calculations for water‐quality assessment. The model estimated that the average annual sediment yield was 773 kg ha?1 yr ?1 compared with a measured value of 641 kg ha?1 yr?1. The P‐factors estimated in this study were comparable with predicted values obtained with the revised universal soil loss equation (RUSLE2). The P‐factors from this study have the potential to be directly used as input in hydrological models, such as the soil and water assessment tool (SWAT), or in soil conservation planning where only conventional digital elevation models (DEMs) are available. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

16.
Incremental dynamic analysis (IDA) leads to curves expressed in terms of structural response versus intensity, commonly known as the IDA curves. It is known that implementation of IDA usually involves significant computational effort and most often significant scaling of the original records to various intensity levels. Employing as the performance variable the critical demand to capacity ratio (DCR) throughout the structure, which is equal to unity at the onset of the limit state, facilitates the identification of the intensity values at the onset of a desired limit state and hence the implementation of the IDA procedure. Employing the structural response to un‐scaled records and the corresponding regression‐based response predictions (a.k.a., the “Cloud Analysis”) helps in identifying the range of intensity values corresponding to demand to capacity ratio values in the vicinity of unity. The Cloud to IDA procedure for structural fragility assessment is proposed on the premise of exploiting the Cloud Analysis results to obtain the IDA curves both with minimum number of analyses and minimum amount of scaling. The transverse frame of a shear‐critical 7‐story older RC building in Van Nuys, CA, which is modeled in Opensees with fiber‐section considering the flexural‐shear‐axial interactions and the bar slip, is employed as a case study. It is demonstrated, by comparing the results to IDA and other state of the art non‐linear dynamic procedures based on no scaling or spectral‐shape‐compatible scaling, that the Cloud to IDA procedure leads to reliable results in terms of structural fragility and risk for the prescribed limit state.  相似文献   

17.
In this paper, we derive analytical expressions for one‐way and two‐way kinematical parameters in elliptical tilted transverse isotropy media. We show that the homogeneous elliptical tilted transverse isotropy models result in hyperbolic moveout with a reflection point sideslip x0, which can be considered as an additional traveltime parameter for one‐way wave propagation. For homogeneous elliptical tilted transverse isotropy models we show that the inversion of one‐way traveltime parameters suffers from the ambiguity for large tilts. It is shown that the accuracy of the inversion is sensitive to the error in x0. We also derive and invert the traveltime parameters for a vertically heterogeneous elliptical tilted transverse isotropy model with a tilt gradually changing with depth. The a priori knowledge of parameter δ is very important for inversion. The wrong choise of this parameter results in significant errors in inverted model parameters.  相似文献   

18.
A 487‐year annually laminated (varved) sediment record from Nicolay Lake, Cornwall Island, in the Canadian High Arctic was evaluated to determine the impact that years with high sediment yields had on sediment yields in subsequent years. All of the 40 largest years showed evidence for increased sediment yield in the subsequent 10–30 years. The positive anomalies in lagging years were approximately scaled according to the size of the initiating year, although many intermediate years (25‐ to 100‐year recurrence) showed weak or variable responses. The smallest events considered (10‐ to 25‐year recurrence) showed a consistent, but low‐amplitude response. Additionally the 10‐year events revealed frequent negative sediment yield anomalies in the preceding decade. This behaviour was interpreted as a frequent sediment activation cycle initiated by the modest year, and leading to sediment yield hysteresis lasting 15–25 years. The largest years (greater than 50‐year recurrence) showed consistently above‐average sediment yields in the preceding decade, in part due to the frequent occurrence of moderate (Q10) years. It is hypothesized that temporary storage of sediment and previous initiation of erosion sites resulted in extraordinary sediment yields during intense summer rainfall events. This study demonstrates the potential use of varved lake sediment records to improve our understanding of long‐term sediment dynamics. These records present an opportunity to further develop and test sediment dynamic and routing models to gain insight into the interaction of time and space in fluvial and sediment delivery processes. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

19.
Simulation approaches employed in suspended sediment processes are important in the areas of water resources and environmental engineering. In the current study, neuro‐fuzzy (NF), a combination of wavelet transform and neuro‐fuzzy (WNF), multi linear regression (MLR), and the conventional sediment rating curve (SRC) models were considered for suspended sediment load (S) modeling in a gauging station in the USA. In the proposed WNF model, the discrete wavelet analysis was linked to a NF approach. To achieve this aim, the observed time series of river flow discharge (Q) and S were decomposed to sub time series at different scales by discrete wavelet transform. Afterwards, the effective sub time series were added together to obtain a useful Q and S time series for prediction. Eventually, the obtained total time series were imposed as inputs to the NF method for daily S prediction. The results illustrated that the predicted values by the proposed WNF model were in good agreement with the observed S values and gave better results than other models. Furthermore, the WNF model satisfactorily estimated the cumulative suspended sediment load and produced relatively reasonable predictions for extreme values of S, while NF, MLR, and SRC models provided unacceptable predictions.  相似文献   

20.
Output‐only system identification is developed here towards assessing current modal dynamic properties of buildings under seismic excitation. Earthquake‐induced structural response signals are adopted as input channels for two different Operational Modal Analysis (OMA) techniques, namely, a refined Frequency Domain Decomposition (rFDD) algorithm and an improved Data‐Driven Stochastic Subspace Identification (SSI‐DATA) procedure. Despite that short‐duration, non‐stationary, earthquake‐induced structural response signals shall not fulfil traditional OMA assumptions, these implementations are specifically formulated to operate with seismic responses and simultaneous heavy damping (in terms of identification challenge), for a consistent estimation of natural frequencies, mode shapes, and modal damping ratios. A linear ten‐storey frame structure under a set of ten selected earthquake base‐excitation instances is numerically simulated, by comparing the results from the two identification methods. According to this study, best up‐to‐date, reinterpreted OMA techniques may effectively be used to characterize the current dynamic behaviour of buildings, thus allowing for potential Structural Health Monitoring approaches in the Earthquake Engineering range.  相似文献   

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

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