首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 78 毫秒
1.
Predicting the geometry of channels and alluvial rivers is of primary importance in river engineering science. Appropriately designing channels and predicting stable river cross‐sections can decrease costs and prevent the destruction of installations and agricultural land by rivers. Consequently, researchers have applied different empirical and regression methods to achieve relations for predicting stable channel and river geometry. In this study, Group Method of Data Handling ]GMDH) models are used to predict three geometric variables of stable channels, namely width (w), depth (h) and slope (s). The effect of different input parameters, such discharge (Q), median grain size (d50) and the Shields parameter (τ*) on the GMDH models is assessed with regard to predicting stable channel geometry. The results indicate that the GMDH model with mean absolute percentage error (MAPE) of 5.53%, 4.05% and 4.89% for channel width, depth and slope prediction respectively, exhibits good accuracy. Moreover, a comparison of the GMDH models with previous theoretical equations (based on regression analysis) indicates the superiority of GMDH model performance, with error reductions of one‐fifth, one‐eighth and one‐sixth compared with the regression equations for channel width, depth and slope prediction, respectively. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

2.
It has generally been assumed that diffusive sediment transport on soil‐mantled hillslopes is linearly dependent on hillslope gradient. Fieldwork was done near Santa Barbara, California, to develop a sediment transport equation for bioturbation by the pocket gopher (Thomomys bottae) and to determine whether it supports linear diffusion. The route taken by the sediment is divided into two parts, a subsurface path followed by a surface path. The first is the transport of soil through the burrow to the burrow opening. The second is the discharge of sediment from the burrow opening onto the hillslope surface. The total volumetric sediment flux, as a function of hillslope gradient, is found to be: qs (cm3 cm−1 a−1) = 176(dz/dx)3 − 189(dz/dx)2 + 68(dz/dx) + 34(dz/dx)0·4. This result does not support the use of linear diffusion for hillslopes where gopher bioturbation is the dominant mode of sediment transport. A one‐dimensional hillslope evolution program was used to evolve hillslope profiles according to non‐linear and linear diffusion and to compare them to a typical hillslope. The non‐linear case more closely resembles the actual profile with a convex cap at the divide leading into a straight midslope section. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

3.
The local reach gradient of small gravel bed rivers (drainage area 0-8-110 km2) in the Eifel, West Germany, is adjusted to transport the river bed sediments. Transport of gravel becomes possible under high flow conditions (Shields entrainment factor ≈-03). Mean bed material size for riffle sections increases with distance downstream. For small drainage areas channel slope is a negative exponential function of drainage area, while for the larger region the additional influence of bedload size has to be considered. Good agreement with Hack's data (1957) for Virginia and Maryland, U.S.A., is achieved (S = 0.0066 (D50/A)- 40., r = 0.67).  相似文献   

4.
Since 1986, with a sharp decrease in water dis-charges, the Yellow River has entered a period charac-terized by low discharges and seasonally occurring dry-ups[1,2]. Since 1999, more strict management of water diversion has been imposed, and therefore the dry-ups have been well under control. However, the lower reaches of the Yellow River is still predominated by low-discharges, and has become a man-induced shrinking river. In the past 40 years, significant effect of soil and water conservat…  相似文献   

5.
Previous studies of alluvial rivers have shown that channel patterns form a continuum controlled by interactions among factors such as gradient, discharge, sediment size, and bank strength. Data from channels in the permanent wetlands of the Okavango Delta add to these ?ndings by focusing on pattern transitions in channels with banks formed by sedges and grasses that are rooted in peat and underlain by unconsolidated sediment. Channels are well de?ned, and transport ?ne–medium sand as bedload between the vertical, vegetation‐lined banks. Water depths, velocities, grain sizes, and bankline vegetation do not vary signi?cantly or systematically downstream, but the permeable banks allow water to leak from the channels, contributing to an overall downstream decrease in discharge and width. In addition, as the Okavango River ?ows from the <12 km wide ‘Panhandle’ and splits into distributaries in the broader ‘Fan’, valley gradient steepens by c. 60 per cent. These downstream changes result in channel pattern adjustments. In the Panhandle, the Okavango River is a relatively wide (c. 30–100 m), actively meandering, sinuous channel (P > 2·0), but further downstream in the Fan, the narrower (<40 m) distributaries follow laterally stable, less sinuous (‘straight’) courses (P < 1·75). Some channel pattern discrimination diagrams based on simple indices of gradient, discharge, sediment size or stream power are inadequate for analysing the meandering–straight transition in the Okavango but Parker's (1976) approach, based on ratios of depth–width and slope–Froude number, accurately characterizes the transition. Our ?eld observations, combined with the results from previous experimental studies, suggest that in relatively wide channels (w/d > 10), thalweg meandering results in scour of the unconsolidated sediment at the bank base, leading to undermining and collapse of the vegetation, and to slow meander migration. However, as channels narrow downstream (w/d < 10) with discharge losses, proportionally increasing sidewall drag exerted by bankline vegetation suppresses thalweg meandering and bank scour, and channels follow stable, less sinuous courses. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

6.
Groundwaters feeding travertine‐depositing rivers of the northeastern segment of the Barkly karst (NW Queensland, Australia) are of comparable chemical composition, allowing a detailed investigation of how the rate of downstream chemical evolution varies from river to river. The discharge, pH, temperature, conductivity and major‐ion concentrations of five rivers were determined by standard field and laboratory techniques. The results show that each river experiences similar patterns of downstream chemical evolution, with CO2 outgassing driving the waters to high levels of calcite supersaturation, which in turn leads to widespread calcium carbonate deposition. However, the rate at which the waters evolve, measured as the loss of CaCO3 per kilometre, varies from river to river, and depends primarily upon discharge at the time of sampling and stream gradient. For example, Louie Creek (Q = 0·11 m3 s?1) and Carl Creek (Q = 0·50 m3 s?1) have identical stream gradients, but the loss of CaCO3 per kilometre for Louie Creek is twice that of Carl Creek. The Gregory River (Q = 3·07 m3 s?1), O'Shanassy River (Q = 0·57 m3 s?1) and Lawn Hill Creek (Q = 0·72 m3 s?1) have very similar gradients, but the rate of hydrochemical evolution of the Gregory River is significantly less than either of the other two systems. The results have major implications for travertine deposition: the stream reach required for waters to evolve to critical levels of calcite supersaturation will, all others things being equal, increase with increasing discharge, and the length of reach over which travertine is deposited will also increase with increasing discharge. This implies that fossil travertine deposits preserved well downstream of modern deposition limits are likely to have been formed under higher discharge regimes. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

7.
Fluvial seed dispersal considers both the transport and deposition of seeds where channel geomorphic structures, hydrology and seed dispersal traits contribute to transport times and depositional locations. This study examines the influence of stream flow patterns on fluvial seed dispersal of buoyant white alder (Alnus rhombifolia) seeds by applying a one‐dimensional transport model. Conceptually, the model separates the stream into two components: (i) the main channel where the seeds are transported downstream; and (ii) the transient storage zone where seeds are temporarily detained or deposited on the river bank. Transport processes are characterized by an advection–dispersion equation which is coupled to a transient storage model using an exponential decay term. The model parameters: longitudinal dispersion (DL), exchange coefficient (α), main channel area (A) and storage zone (As) are estimated based on field experiments conducted in a confined, bedrock‐gravel bed river with pool‐riffle morphology located in coastal northern California. The riparian zone is inhabited by Alnus rhombifolia that disperse buoyant seeds in mid‐spring coinciding with the end of the wet, Mediterranean season. Artificial seeds, with similar traits of buoyancy and density to alder seeds, were used to quantify transport times and depositional locations. Preferential deposition resulted in stream reaches with larger As, high As/A ratios, and faster exchange coefficients corresponding to divergent stream flow (back‐eddies, re‐circulating flow, flow expansions) caused by geomorphic structures such as the ends of bar/riffle features and bends in the stream. The results demonstrate the importance of transient storage for seed transport and depositional processes. Morphological features that increase a channel's complexity create complex flow structures that detain seeds and provide a greater opportunity for deposition to occur. The model provides a simplification of river hydraulics to represent dispersal dynamics and lends itself to further understanding of hydrochory processes and associated population structure. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

8.
Sediment supply (Qs) is often overlooked in modelling studies of landscape evolution, despite sediment playing a key role in the physical processes that drive erosion and sedimentation in river channels. Here, we show the direct impact of the supply of coarse-grained, hard sediment on the geometry of bedrock channels from the Rangitikei River, New Zealand. Channels receiving a coarse bedload sediment supply are systematically (up to an order of magnitude) wider than channels with no bedload sediment input for a given discharge. We also present physical model experiments of a bedrock river channel with a fixed water discharge (1.5 l min−1) under different Qs (between 0 and 20 g l−1) that allow the quantification of the role of sediment in setting the width and slope of channels and the distribution of shear stress within channels. The addition of bedload sediment increases the width, slope and width-to-depth ratio of the channels, and increasing sediment loads promote emerging complexity in channel morphology and shear stress distributions. Channels with low Qs are characterized by simple in-channel morphologies with a uniform distribution of shear stress within the channel while channels with high Qs are characterized by dynamic channels with multiple active threads and a non-uniform distribution of shear stress. We compare bedrock channel geometries from the Rangitikei and the experiments to alluvial channels and demonstrate that the behaviour is similar, with a transition from single-thread and uniform channels to multiple threads occurring when bedload sediment is present. In the experimental bedrock channels, this threshold Qs is when the input sediment supply exceeds the transport capacity of the channel. Caution is required when using the channel geometry to reconstruct past environmental conditions or to invert for tectonic uplift rates, because multiple configurations of channel geometry can exist for a given discharge, solely due to input Qs. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd  相似文献   

9.
This paper assesses the mechanisms and pathways by which peat blocks are eroded and transported in upland fluvial systems. Observations and experiments from the north Pennines (UK) have been carried out on two contrasting river systems. Mapping of peat block distributions and appraisal of reach‐based sediment budgets clearly demonstrates that macro‐size peat is an important stream load component. In small streams block sizes can approximate the channel width and much of the peat is transported overbank. Local ‘peat jams’ and associated mineral deposition may provide an important component of channel storage. In larger systems peat blocks rapidly move down‐channel and undergo frequent exchanges between bed and bank storage. Results of peat block tracing using painted blocks indicate that once submerged, blocks of all sizes are easily transported and blocks break down rapidly by abrasion. Vegetation and bars play an important role in trapping mobile peat. In smaller streams large block transport is limited by channel jams. Smaller blocks are transported overbank but exhibit little evidence of downstream fining. In larger rivers peat blocks are more actively sorted and show downstream reduction in size from source. A simple model relating peat block diameter (Dp) to average flow depth (d) suggests three limiting transport conditions: flotation (Dp < d), rolling (d < Dp > d/2) and deposition (Dp > d/2). Experiments demonstrate that peat block transport occurs largely by rolling and floating and the transport mechanism is probably controlled by relative flow depth (d/Dp ratio). Transport velocity varies with transport mechanism (rolling is the slowest mode) and transport lengths increase as flow depth increases. Abrasion rates vary with the transport mechanism. Rolling produces greater abrasion rates and more rounded blocks. Abrasion rates vary from 0 to 10 g m?1 for blocks ranging in mass from 10 to 6000 g. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

10.
It was found in previous studies that groundwater levels may fluctuate as a temporal fractal. In this study numerical simulations of groundwater level fluctuations in an unconfined aquifer near a river were conducted to investigate the effects of aquifer heterogeneity and river stage variations on the fractal behavior of the water levels, h(t). Groundwater recharge was taken to be a white-noise process. The aquifer heterogeneity was simulated with a second-order stationary field of hydraulic conductivity (K) with an exponential variogram model. The results showed that groundwater levels fluctuate as a temporal fractal in both homogeneous and heterogeneous aquifers as long as K is less than 10 m/d. Most aquifers may indeed act as a fractal filter which takes a random non-fractal recharge inputs and produces a fractal responses of groundwater level fluctuations. A crossover in temporal scaling of h(t) may appear in more permeable aquifers. Fluctuations of the groundwater level in a homogeneous aquifer are dominated by the recharge process when the river stage is constant or by the river stage variations when the river stage varies in highly permeable aquifers. Heterogeneity plays an important role in the temporal scaling of h(t) in more permeable aquifers: the stronger the heterogeneity, the stronger the temporal scaling of h(t).  相似文献   

11.
Spatial analysis of stream power using GIS: SLk anomaly maps   总被引：2，自引：0，他引：2  
The stream length‐gradient index (SL) shows the variation in stream power along river reaches. This index is very sensitive to changes in channel slope, thus allowing the evaluation of recent tectonic activity and/or rock resistance. Nevertheless, the comparison of SL values from rivers of different length is biased due to the manner in which the index is formulated, thus making correlations of SL anomalies along different rivers difficult. Therefore, when undertaking a comparison of SL values of rivers of different lengths, a normalization factor must be used. The graded river gradient (K) has already been used in some studies to normalize the SL index. In this work, we explore the relationships between the graded river gradient (K), the SL index and the stream power, proposing the use of a re‐named SLk index, which enables the comparison of variable‐length rivers, as well as the drawing of SLk anomaly maps. We present here a GIS‐based procedure to generate SLk maps and to identify SLk anomalies. In order to verify the advantages of this methodology, we compared an SLk map of the NE border of the Granada basin with both simple river profile–knickpoint identification and with an SL map. The results show that the SLk map supplies good results with defined anomalies and suitably reflects the main tectonic and lithological features of the study area. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

12.
Two reaches of Aguapeí River, a left‐bank tributary of the Paraná River in western São Paulo state, Brazil, were studied with the objective of assessing the role of bend curvature on channel migration in this wet‐tropical system and examining if land‐use changes or ENSO (El Niño Southern Oscillation) driven climate anomalies over nearly half a century have changed migration behaviour and planform geometry. Meander‐bend migration rates and morphometric parameters including meander‐bend curvature, sinuosity, meander wavelength and channel width, were measured and the frequency of bend cutoffs was analysed in order to determine the rate of change of channel adjustment over a 48 year period to 2010. Results show that maximum average channel migration rates occur in bends with curvatures of about 2–3 rc/w, similar to other previously studied temperate and subarctic freely meandering rivers although not as pronounced and with a tendency to favour tighter curvature. From 1962 to 2010 the Aguapeí River has undergone a significant reduction in sinuosity, a shift from tightly curving to more open bends, an overall decline in channel migration rates, an associated decrease in the frequency of neck‐cutoffs and an overall increase in channel width. As the majority of the drainage basin (96%) was already deforested in 1962, channel form and process changes were, unlike an interpretation for an adjacent river system, not attributed to altered land‐use but rather to a sharp ENSO‐driven increase in the magnitude of peak flow‐discharges of some 32% since 1972. In summary, this research revealed that recent climate and associated flow regime changes are having a pronounced effect on river channel behaviour in the Aguapeí River investigated here. Copyright © 2018 John Wiley & Sons, Ltd.  相似文献   

13.
The capability of peatland ecosystems to regulate evapotranspiration (ET) following wildfire is a key control on the resilience of their globally important carbon stocks under future climatic conditions. Evaporation dominates post-fire ET, with canopy and sub-canopy removal restricting transpiration and increasing evaporation potential. Therefore, in order to project the hydrology and associated stability of peatlands to a diverse range of post-fire weather conditions and future climates the regulation of evaporation must be accurately parameterised in peatland ecohydrological models. To achieve this, we measure the surface resistance (rs) to evaporation over the growing season one year post-fire within four zones of a boreal peatland that burned to differing depths, relating rs to near surface soil tensions. We show that the magnitude and temporal variability in rs varies with burn severity. At the peatland scale, rs and near-surface tension correlates non-linearly. However, at the point scale no relationship was evident between temporal variations in rs and near-surface tension across all burn severities; in part due to the limited fluctuation in near-surface tensions and the precision of rs measurements. Where automated measurements enabled averaging of errors, the relationship between near-surface tension and rs switched between periods of strong and weak correlation within a burned peat hummock. This relationship, when strong, deviated from that obtained under steady state laboratory conditions; increases in rs were more sensitive to fluctuations in near-surface tension under dynamic field conditions. Calculating soil vapour densities directly from near-surface tensions is shown to require calibration between peat types and provides little if any benefit beyond the derivation of empirical relationships between rs and measured soil tension. Thus, we demonstrate important spatiotemporal fluctuations in post-fire rs that will be key to regulating post-fire peatland hydrology, but highlight the complex challenges in effectively parameterising this important underlying control of near-surface tensions within hydrological simulations.  相似文献   

14.
  相似文献   

15.
In meandering rivers, the local channel migration rate increases with increasing bend sharpness until it reaches a maximum at a certain critical value of the bend sharpness. Beyond this critical value, the migration rate decreases if bend sharpness increases. Similarly, reach‐averaged migration rates attain a maximum at a certain river sinuosity. This work investigates the physics of these phenomena by comparing the results of two physics‐based models of different complexity, in which the migration rates are proportional to the near‐bank flow velocity excess. In the computational tests the river was allowed to meander progressively, starting from an almost straight planimetry. Both models reproduced the observed peak in the curve describing the local migration rate as a function of the ratio radius of curvature‐channel width (R/B), with a rising limb at lower R/B values and a falling limb at higher R/B values. The rising limb can be explained by the decrease in relative lag distance between near‐bank flow velocity and forcing curvature as R/B increases. The falling limb results from the decrease in local channel curvature and near‐bank flow velocity excess. Since the models do not include flow separation, the results indicate that this phenomenon is not needed to explain the decrease of channel migration rates in sharp bends. The models reproduced also the peak in the curve describing the reach‐averaged migration rates as a function of river sinuosity The increase and then decrease of reach‐averaged migration rates as sinuosity increases appears to be mainly caused by the variation of the reach‐averaged value of the ratio R/B. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

16.
Abstract

Rivers have been channelized, deepened and constrained by embankments for centuries to increase agricultural productivity and improve flood defences. This has decreased the hydrological connectivity between rivers and their floodplains. We quantified the hydrological regime of a wet grassland meadow prior to and after the removal of river embankments. River and groundwater chemistry were also monitored to examine hydrological controls on floodplain nutrient status. Prior to restoration, the highest river flows (～2 m3 s?1) were retained by the embankments. Under these flow conditions the usual hydraulic gradient from the floodplain to the river was reversed so that subsurface flows were directed towards the floodplain. Groundwater was depleted in dissolved oxygen (mean: 0.6 mg O2 L?1) and nitrate (mean: 0.5 mg NO3 ?-N L?1) relative to river water (mean: 10.8 mg O2 L?1 and 6.2 mg NO3 ?-N L?1, respectively). Removal of the embankments has reduced the channel capacity by an average of 60%. This has facilitated over-bank flow which is likely to favour conditions for improved flood storage and removal of river nutrients by floodplain sediments.

Editor Z.W. Kundzewicz; Associate editor K. Heal

Citation Clilverd, H.M., Thompson, J.R., Heppell, C.M., Sayer, C.D., and Axmacher, J.C., 2013. River–floodplain hydrology of an embanked lowland Chalk river and initial response to embankment removal. Hydrological Sciences Journal, 58 (3), 627–650.  相似文献   

17.
Soil detachment in concentrated flow is due to the dislodging of soil particles from the soil matrix by surface runoff. Both aggregate stability and shear strength of the topsoil reflect the erosion resistance of soil to concentrated runoff, and are important input parameters in predicting soil detachment models. This study was conducted to develop a formula to predict soil detachment rate in concentrated flow by using the aggregate stability index (As), root density (Rd) and saturated soil strength (σs) in the subtropical Ultisols region of China. The detachment rates of undisturbed topsoil samples collected from eight cultivated soil plots were measured in a 3.8 m long, 0.2 m wide hydraulic flume under five different flow shear stresses (τ = 4.54, 9.38, 15.01, 17.49 and 22.54 Pa). The results indicated that the stability index (As) was well related with soil detachment rate, particularly for results obtained with high flow shear stress (22.54 Pa), and the stability index (As) has a good linear relationship with concentrated flow erodibility factors (Kc). There was a positive linear relationship between saturated soil strength (σs) and critical flow shear stress (τc) for different soils. A significant negative exponential relationship between erodibility factors (Kc) and root density (Rd) was detected. This study yielded two prediction equations that allowed comparison of their efficiency in assessing soil detachment rate in concentrated flow. The equation including the root density (Rd) may have a better correlation coefficient (R2 = 0.95). It was concluded that the formula based on the stability index (As), saturated soil strength (σs) and root density (Rd) has the potential to improve methodology for assessing soil detachment rate in concentrated flow for the subtropical Chinese Ultisols. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

18.
A river at equilibrium is described by a statistically-stationary mean bed elevation profile that arises in response to steady supplies of relief, water and sediment. Outside of the profile shape, how is the equilibrium state of a river most reliably identified and rigorously defined? Motivated by a proposed link between equilibrium and physical processes, we use scaling theory to develop the dimensionless channel response number ξ=KUb/Up. ξ is a metric for the local disequilibrium state of gravel-bed mountain streams, which reflects a balance between the rate of topographic adjustment Ub, and the rate of bed sediment texture adjustment Up. The coefficient K can take one of two forms depending on choice of length scale for topographic adjustment. We hypothesize that equilibrium occurs where and when ξ≈O(1), and consequently, disequilibrium is the more general state captured by conditions of ξ≉O(1). The rates Ub and Up are controlled by the mechanics of sediment deposition and entrainment at the local scale of the channel width. The extent to which either process regulates disequilibrium depends on the bed strength, which is set by the time-varying grain size distribution and packing. We use flume experiments to understand ξ and find that in the limit ξ>>1, the time-varying response of an experimental channel depends sensitively on the spatially-averaged bed shear stress ratio τ/τref. When τ/τref≈1.5, Ub was the dominant control on disequilibrium. However, when τ/τref≈2.0, Up contributed more significantly to disequilibrium. These results suggest that after an upstream supply perturbation, the equilibrium timescale is governed by Up, which we show is consistent with expectations from linear damping theory. Our experimental test of ξ is promising, but inconclusive with respect to our hypothesis. This uncertainty can be readily addressed with numerical or additional physical experiments. © 2019 John Wiley & Sons, Ltd.  相似文献   

19.
Remote sensing of volumetric storage changes in lakes   总被引：1，自引：0，他引：1  
Three‐dimensional remote sensing promises a giant leap forward for surface‐water hydrology in much the same way that radar altimetry transformed physical oceanography. However, the complex geometries of small terrestrial water bodies introduce difficulties, particularly with respect to trade‐offs between changing water depth and inundation area. We use in situ measurements of water‐surface stage (ΔH/dt) and remotely‐sensed area (A) to compute time varying storage changes (ΔS) in nine lakes of the Peace‐Athabasca Delta, Canada. Despite their identical geomorphic setting, regression slopes between ΔH and A vary significantly between lakes, primarily from a predictable ‘area‐effect’ but also small bathymetric variations between basins. On average, lateral contraction/expansion (versus stage adjustment) contributes as little as 7% (versus 93%) to as much as 76% (versus 24%) of overall storage change ΔS. We conclude that both surface‐area and ΔH/dt, rather than just either alone, must be measured to confidently estimate ΔS from space. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

20.
A theoretical model is developed for predicting equilibrium alluvial channel form. The concept of greatest relative stability, achieved by maximizing resistance to ?ow in the ?uvial system, is presented as the basis for an optimization condition for alluvial systems. Discharge, sediment supply (quantity and calibre) and valley gradient are accepted as independent governing variates. The model is used to de?ne a dimensionless alluvial state space characterized by aspect ratio (W/d), relative roughness (D/d), and dimensionless shear stress (τ*) or, equivalently, channel slope (S). Each alluvial state exhibits unique values of Froude number and sediment concentration. The range of alluvial states for constant values of relative bank strength (parameterized by an apparent friction angle, ?′) forms a single plane in the state space (W/d, D/d, τ* or S). The scaling relations produced by the model are consistent with laboratory channels exhibiting a range of bank strengths, and with the behaviour of natural channels. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

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

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