首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Experimental analysis on the impact force of viscous debris flow   总被引:1,自引:0,他引:1       下载免费PDF全文
A miniaturized flume experiment was carried out to measure impact forces of viscous debris flow. The flow depth (7.2–11.2 cm), velocity (2.4–5.2 m/s) and impact force were recorded during the experiment. The impact process of debris flow can be divided into three phases by analyzing the variation of impact signals and flow regime. The three phases are the sudden strong impact of the debris flow head, continuous dynamic pressure of the body and slight static pressure of the tail. The variation of impact process is consistent with the change in the flow regime. The head has strong–rapid impact pressure, which is shown as a turbulent‐type flow; the body approximates to steady laminar flow. Accordingly, the process of debris flows hitting structures was simplified to a triangle shape, ignoring the pressure of the tail. In order to study the distribution of the debris flow impact force at different depths and variation of the impact process over time, the impact signals of slurry and coarse particles were separated from the original signals using wavelet analysis. The slurry's dynamic pressure signal appears to be a smooth curve, and the peak pressure is 12–34 kPa when the debris flow head hits the sensors, which is about 1.54 ± 0.36 times the continuous dynamic pressure of the debris flow body. The limit of application of the empirical parameter α in the hydraulic formula was also noted. We introduced the power function relationship of α and the Froude number of debris flows, and proposed a universal model for calculating dynamic pressure. The impact pressure of large particles has the characteristic of randomness. The mean frequency of large particles impacting the sensor is 210 ± 50–287 ± 29 times per second, and it is 336 ± 114–490 ± 69 times per second for the debris flow head, which is greater than that in the debris flow body. Peak impact pressure of particles at different flow depths is 40–160 kPa, which is 3.2 ± 1.5 times the impact pressure of the slurry at the bottom of the flow, 3.1 ± 0.9 times the flow in the middle, and 3.3 ± 0.9 times the flow at the surface. The differences in impact frequency indicate that most of the large particles concentrate in the debris flow head, and the number of particles in the debris flow head increases with height. This research supports the study of solid–liquid two phase flow mechanisms, and helps engineering design and risk assessment in debris flow prone areas. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

2.
A new method to predict the runout of debris flows is presented. A data base of documented sediment‐transporting events in torrent catchments of Austria, Switzerland and northern Italy has been compiled, using common classification techniques. With this data we test an empirical approach between planimetric deposition area and event volume, and compare it with results from other studies. We introduce a new empirical relation to determine the mobility coefficient as a function of geomorphologic catchment parameters. The mobility coefficient is thought to reflect some of the flow properties during the depositional part of the debris‐flow event. The empirical equations are implemented in a geographical information system (GIS) based simulation program and combined with a simple flow routing algorithm, to determine the potential runout area covered by debris‐flow deposits. For a given volume and starting point of the deposits, a Monte‐Carlo technique is used to produce flow paths that simulate the spreading effect of a debris flow. The runout zone is delineated by confining the simulated potential spreading area in the down slope direction with the empirically determined planimetric deposition area. The debris‐flow volume is then distributed over the predicted area according to the calculated outflow probability of each cell. The simulation uses the ARC‐Objects environment of ESRI© and is adapted to run with high resolution (2·5 m × 2·5 m) digital elevation models, generated for example from LiDAR data. The simulation program called TopRunDF is tested with debris‐flow events of 1987 and 2005 in Switzerland. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

3.
Tree‐ring records from conifers have been regularly used over the last few decades to date debris‐flow events. The reconstruction of past debris‐flow activity was, in contrast, only very rarely based on growth anomalies in broad‐leaved trees. Consequently, this study aimed at dating the occurrence of former debris flows from growth series of broad‐leaved trees and at determining their suitability for dendrogeomorphic research. Results were obtained from gray alder (Alnus incana (L.) Moench), silver birch and pubescent birch (Betula pendula Roth and Betula pubescens Ehrh.), aspen (Populus tremula L.), white poplar, black poplar and gray poplar (Populus alba L., Populus nigra L. and Populus x canescens (Ait.) Sm.), goat willow (Salix caprea L.) and black elder (Sambucus nigra L.) injured by debris‐flow activity at Illgraben (Valais, Swiss Alps). Tree‐ring analysis of 104 increment cores, 118 wedges and 93 cross‐sections from 154 injured broad‐leaved trees allowed the reconstruction of 14 debris‐flow events between AD 1965 and 2007. These events were compared with archival records on debris‐flow activity at Illgraben. It appears that debris flows are very common at Illgraben, but only very rarely left the channel over the period AD 1965–2007. Furthermore, analysis of the spatial distribution of disturbed trees contributed to the identification of six patterns of debris‐flow routing and led to the determination of preferential breakout locations of events. The results of this study demonstrate the high potential of broad‐leaved trees for dendrogeomorphic research and for the assessment of the travel distance and lateral spread of debris‐flow surges. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

4.
In many Alpine catchments the monitoring and hazard mitigation of debris‐flow events require enormous economic and social resources. To confront these problems, a subjective estimation of the most hazardous zones of the basin could be useful in the best, sustainable planning of protective measures. In this paper, a new methodology is proposed that develops a Management Priority Index (MPI) to rank sediment source areas by their quantitative capability to deliver debris‐flow volumes to a point of interest within the catchment. The MPI sets the intervention priority based on a combination of three sub‐indicators: a susceptibility indicator evaluating the overall catchment predisposition to generate debris flow, a triggering indicator and a volume budget indicator assessing the rate of deliverable volume to a selected outlet. MPI was applied to the basin of the Rio Gadria catchment (Venosta Valley, Bolzano, Italy), an alpine basin with an unlimited sediment supply that is characterized by multiple, very active, shallow landslides and bare soil zones. The proposed ranking method was successfully verified using post‐event surveys and through evidence from consolidation check dams built over many years in the basin. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

5.
Debris flows can grow greatly in size by entrainment of bed material, enhancing their runout and hazardous impact. Here, we experimentally investigate the effects of debris‐flow composition on the amount and spatial patterns of bed scour and erosion downstream of a fixed to erodible bed transition. The experimental debris flows were observed to entrain bed particles both grain by grain and en masse, and the majority of entrainment was observed to occur during passage of the flow front. The spatial bed scour patterns are highly variable, but large‐scale patterns are largely similar over 22.5–35° channel slopes for debris flows of similar composition. Scour depth is generally largest slightly downstream of the fixed to erodible bed transition, except for clay‐rich debris flows, which cause a relatively uniform scour pattern. The spatial variability in the scour depth decreases with increasing water, gravel (= grain size) and clay fraction. Basal scour depth increases with channel slope, flow velocity, flow depth, discharge and shear stress in our experiments, whereas there is no correlation with grain collisional stress. The strongest correlation is between basal scour and shear stress and discharge. There are substantial differences in the scour caused by different types of debris flows. In general, mean and maximum scour depths become larger with increasing water fraction and grain size, and decrease with increasing clay content. However, the erodibility of coarse‐grained experimental debris flows (gravel fraction = 0.64) is similar on a wide range of channel slopes, flow depths, flow velocities, discharges and shear stresses. This probably relates to the relatively large influence of grain‐collisional stress to the total bed stress in these flows (30–50%). The relative effect of grain‐collisional stress is low in the other experimental debris flows (<5%), causing erosion to be largely controlled by basal shear stress. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

6.
Servo‐hydraulic actuators have been widely used for experimental studies in engineering. They can be controlled in either displacement or force control mode depending on the purpose of a test. It is necessary to control the actuators in real time when the rate‐dependency effect of a test specimen needs to be accounted for under dynamic loads. Real‐time hybrid simulation (RTHS) and effective force testing (EFT) method, which can consider the rate‐dependency effect, have been known as viable alternatives to the shake table testing method. Due to the lack of knowledge in real‐time force control, however, the structures that can be tested with RTHS and EFT are fairly limited. For instance, satisfying the force boundary condition for axially stiff members is a challenging task in RTHS, while EFT has a difficulty to be implemented for nonlinear structures. In order to resolve these issues, this paper introduces new real‐time force control methods utilizing the adaptive time series (ATS) compensator and compliance springs. Unlike existing methods, the proposed force control methods do not require the structural modeling of a test structure, making it easy to be implemented especially for nonlinear structures. The force tracking performance of the proposed methods is evaluated for a small‐scale steel mass block system with a magneto‐rheological damper subjected to various target forces. Accuracy, time delay, and resonance response of these methods are discussed along with their force control performance for an axially stiff member. Overall, a satisfactory force tracking performance was observed by using the proposed force control methods.  相似文献   

7.
Debris flows generated from landslides are common processes and represent a severe hazard in mountain regions due to their high mobility and impact energy. We investigate the dynamics and the rheological properties of a 90 000 m3 debris‐flow event triggered by a rapid regressive landslide with high water content. Field evidence revealed a maximum flow depth of 7–8 m, with an estimated peak discharge of 350–400 m3 s?1. Depositional evidence and grain‐size distribution of the debris pose the debris flow in an intermediate condition between the fluid‐mud and grain‐flow behaviour. The debris‐flow material has silt–clay content up to 15 per cent. The rheological behaviour of the finer matrix was directly assessed with the ball measuring system. The measurements, performed on two samples at 45–63 per cent in sediment concentration by volume, gave values of 3·5–577 Pa for the yield strength, and 0·6–27·9 Pa s for the viscosity. Based on field evidence, we have empirically estimated the yield strength and viscosity ranging between 4000 ± 200 Pa, and 108–134 Pa s, respectively. We used the Flo‐2D code to replicate the debris‐flow event. We applied the model with rheological properties estimated by means of direct measurements and back‐analyses. The results of these models show that the rheological behaviour of a debris‐flow mass containing coarse clasts can not be assessed solely on the contribution of the finer matrix and thus neglecting the effects of direct grain contacts. For debris flows composed by a significant number of coarse clasts a back‐analysis estimation of the rheological parameters is necessary to replicate satisfactorily the depositional extent of debris flows. In these cases, the back‐estimated coefficients do not adequately describe the rheological properties of the debris flow. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

8.
Entrainment of underlying bed sediment by a debris flow can significantly increase the debris‐flow magnitude. To study this phenomenon, a theoretical approach to assessing bed‐sediment entrainment is presented. The approach is based on a static approximation that bed‐sediment entrainment occurs when the shearing stress of the flow is sufficiently high to overcome the basal resistance of the bed sediment. In order to delineate erodible zones in a channel, we analyze the critical condition of this static equilibrium model, and subsequently propose a new concept of a critical line to detect the entrainment reaches in a channel. Considering the spatial and temporal uncertainties of the input parameter, the approach is further incorporated within a Monte Carlo method, and the distribution of entrainment zones and post‐entrainment volumes can be analyzed. This approach is illustrated by back‐analysis of the 2010 Yohutagawa debris‐flow event, Japan. Results from 10 000 trials of Monte Carlo simulation are compared with the in situ surveys. It is shown that the present approach can be satisfactorily used to delineate erodible zones and estimate possible entrainment volume of the event. Discussion regarding the sensitivities and limitations of the approach concludes the paper. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

9.
Dendrogeomorphology was used to investigate past events on a cone affected by both debris flows and snow avalanches. We report on results of 520 cores from 251 injured Larix decidua Mill. and Picea abies (L.) Karst. trees sampled on the Birchbach cone (Swiss Alps). Detailed analysis of tree‐ring sequences allowed dating of 561 growth disturbances in individual trees for a 252 yr period, extending from 1750 to 2002, which could be attributed to 30 different event years. We then localized the position of rows of traumatic resin ducts (TRDs) within the tree ring so as to assess the intra‐seasonal position of damage. In agreement with data on the local growth period, TRDs located at the beginning of the new growth ring were considered the result of avalanche impacts that occurred during the dormant season or in earliest earlywood between late October and early May. In contrast, TRD found in late earlywood or within latewood were considered the result of periglacial debris‐flow activity, as these layers of the tree ring are locally formed between July and early October. For nine out of the 30 reconstructed event years, the intra‐seasonal timing of TRDs indicated that reactions must be the result of past snow avalanche activity. In 19 other event years, TRDs showed that damage has been caused between July and early October and, thus, through debris flows in the Birchbach torrent. Finally, the spatial patterns of trees showing reactions as a result of particular events were assessed so as to approximate the extent of past debris flows and snow avalanches. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

10.
The evolution of a debris‐flow cone depends on a multitude of factors in the hydrogeomorphic system. Investigations of debris‐flow history and cone dynamics in highly active catchments therefore require an integrative approach with a temporal and spatial resolution appropriate for the goals of the study. We present the use of an orthophoto time series to augment standard dendrogeomorphic techniques to describe the spatio‐temporal dynamics of debris flows on a highly active cone in the western Austrian Alps. Analysis of seven orthophotos since 1951 revealed a migration of active deposition areas with a resulting severe loss of forest cover (> 80%) and a mean tree loss per year of 10·4 (range 1·3–16·6 trees per year). Analysis of 193 Pinus mugo ssp. uncinata trees allowed the identification of 161 growth disturbances corresponding to 16 debris flows since 1839 and an average decadal frequency of 0·9 events. As a result of the severe loss of forest cover, we speculate that < 20% of the more recent events were actually captured in the tree‐ring record, giving a decadal return interval of ~7·5 events for a period of 60 years. Based on three annual field observations, it is evident that this catchment (the Bärenrüfe) produces very frequent (< 1 yr), small (in the order of a few 10 to 100 m3) debris flows with minor material relocation. The specific challenges of tree‐ring analysis in this tree species and in highly active environments are explicitly addressed in the discussion and underline the necessity of employing complementary methods of analysis in an integrative manner. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

11.
12.
The quantification of debris‐flow hazard requires estimates of debris‐flow frequency and magnitude. Several methods have been proposed to determine the probable volume of future debris flows from a given basin, but most have neglected to account for debris recharge rates over time, which may lead to underestimation of debris‐flow volumes in basins with rare debris flows. This paper deals with the determination of debris recharge rates in debris‐flow channels based on knowledge of debris storage and the elapsed time since the last debris flow. Data are obtained from coastal British Columbia and a relation is obtained across a sample of basins with similar terrain and climatic conditions. For Rennell Sound on the west coast of the Queen Charlotte Islands, the power‐law relation for area‐normalized recharge rate, Rt, versus elapsed time, te was Rt = 0·23te?0·58 with an explained variance of 75 per cent. A difference in recharge rates may exist between creeks in logged and unlogged forested terrain. The power function for undisturbed terrain was Rt = 0·20te?0·49, while the function for logged areas was Rt = 0·30te?0·77. This result suggests that for the same elapsed time since the last debris flow, clearcut gullies tend to recharge at a slower rate than creeks in old growth forest. This finding requires verification, particularly for longer elapsed times since debris flow, but would have important implications for forest resource management in steep coastal terrain. This study demonstrates that commonly used encounter probability equations are inappropriate for recharge‐limited debris flow channels. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

13.
Detection of surface change is a fundamental task in geomorphology. Terrestrial laser scanners are increasingly used for monitoring surface change resulting from a variety of geomorphic processes, as they allow the rapid generation of high‐resolution digital elevation models. Irrespective of instrument specifics, survey design or data processing, such data are subject to a finite level of ambiguity in position measurement, a consideration of which must be taken into account when deriving change. The propagation of errors is crucial in change detection because even very small uncertainties in elevation can produce large uncertainties in volume when extrapolated over an area of interest. In this study we propose a methodology to detect surface change and to quantify the resultant volumetric errors in areas of complex topography such as channels, where data from multiple scan stations must be combined. We find that a commonly proposed source of error – laser point elongation at low incidence angles – has a negligible effect on the quality of the final registered point cloud. Instead, ambiguities in elevation inherent to registered datasets have a strong effect on our ability to detect and measure surface change. Similarly, we find that changes in surface roughness between surveys also reduce our ability to detect change. Explicit consideration of these ambiguities, when propagated through to volume calculations, allows us to detect volume change of 87 ± 5 m3, over an area of ~ ?4900 m2, due to passage of a debris flow down a 300 m reach of the Illgraben channel in Switzerland. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

14.
The effect of the small fraction of clays on the rheological behaviour of alpine debris flow is poorly understood. This is partly due to the complexity of the debris flow mineralogy and the broad particle size distribution. This study has investigated this issue by simulating an alpine debris flow with a mixture of well characterized fractions and then varying the clay fraction composition. Four samples were tested, ranging from a clay fraction made up of only kaolinite (1:1 type clay) to samples where 80 per cent of the kaolinite is replaced by bedeillite (a 2:1 type clay similar to smectite). Changing the content of 2:1 type clay has a strong influence on the behaviour of the whole material, despite its low weight fraction of around 2 per cent. The tests carried out on these reconstituted materials were compared with the results obtained for natural debris flow materials and showed some common trends: in particular, the rheological parameters for materials with and without 2:1 clays with respect to yield stress as a function of solid content. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

15.
The rheology of debris flows is difficult to characterize owing to the varied composition and to the uneven distribution of the components that may range from clay to large boulders, in addition to water. Few studies have addressed debris flow rheology from observational, experimental, and theoretical viewpoints in conjunction. We present a coupled rheological‐numerical model to characterize the debris flows in which cohesive and frictional materials are both present. As a first step, we consider small‐scale artificial debris flows in a flume with variable percentages of clay versus sand, and measure separately the rheological properties of sand–clay mixtures. A comparison with the predictions of a modified version of the numerical model BING shows a reasonable agreement between measurements and simulations. As application to a field case, we analyse a recent debris flow that occurred in Fjærland (Western Norway) for which much information is now available. The event was caused by a glacial lake outburst flood (GLOF) originating from the failure of a moraine ridge. In a previous contribution (Breien et al., Landslides, 2008 , 5: 271–280) we focused on the hydrological and geomorphological aspects. In particular we documented the marked erosion and reported the change in sediment transport during the event. In contrast to the laboratory debris flows, the presence of large boulders and the higher normal pressure inside the natural debris flow requires the introduction of a novel rheological model that distinguishes between mud‐to–clast supported material. We present simulations with a modified BING model with the new cohesive‐frictional rheology. To account for the severe erosion operated by the debris flow on the colluvial deposits of Fjærland, we also suggest a simple model for erosion and bulking along the slope path. Numerical simulations suggest that a self‐sustaining mechanism could partly explain the extreme growth of debris flows running on a soft terrain.  相似文献   

16.
The dynamics and the surface evolution of a post‐LGM debris‐flow‐dominated alluvial fan (Tartano alluvial fan), which lies on the floor of an alpine valley (Valtellina, Northern Italy), have been investigated by means of an integrated study comprising geomorphological field work, a sedimentological study, photointerpretation, quantitative geomorphology, analysis of ancient to modern cartography and consultation of historical documents and records. The fan catchment meteoclimatic, geological and geomorphological characteristics result in fast rates of geomorphic reorganization of the fan surface (2 km2). The dynamics of the fan are determined by the alternation of low‐return period catastrophic alluvial events dominated by non‐cohesive debris flows triggered by extreme rainstorms which caused aggradation and steepening of the fan and avulsion of its main channel, with periods of low to moderate streamflow discharge punctuated by low‐ to intermediate‐magnitude flood events, causing slower but steady topographic reworking. The most ancient parts of the fan surface date back at least to the first half of the 19th century, but most of the fan surface has been restructured after 1911, mainly during the debris‐flow‐dominated events of 1911 and 1987. Phases of rapid fan toe incision and fan degradation have been recognized; since the 1930s or 1940s, the Tartano fan has been subjected to a state of deep entrenchment and narrowing of the main trunk channel and distributary area. Post‐Little Ice Age climate change and present‐day surface uplift rates have been considered as possible explanations for the observed geomorphic evolution, but tectonic or climatic controls cannot account for the order of magnitude of the erosional pace. Anthropogenic controls plausibly override the natural ones: in particular, the building of a dam in the late 1920s, about 2 km upstream of the fan, seems to have triggered fan dissection, having altered the sediment discharge through sediment retention. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

17.
Glacier recession and landform development in a debris‐charged glacial landsystem characterized by an overdeepening is quantified using digital photogrammetry, digital elevation model (DEM) construction and mapping of the Icelandic glacier Kvíárjökull for the period 1945–2003. Melting of ice‐cores is recorded by surface lowering rates of 0·8 m yr–1 (1945–1964), 0·3 m yr–1 (1964–1980), 0·015 m yr–1 (1980–1998) and 0·044 m yr–1 (1998–2003). The distribution/preservation of pushed and stacked ice‐cored moraine complexes are determined by the location of the long‐term glacial drainage network in combination with retreat from the overdeepening, into which glacifluvial sediment is being directed and where debris‐rich ice masses are being reworked and replaced by esker networks produced in englacial meltwater pathways that bypassed the overdeepening and connected to outwash fans prograding over the snout. Recent accelerated retreat of Kvíárjökull, potentially due to increased mass balance sensitivity, has made the snout highly unstable, especially now that the overdeepening is being uncovered and the snout flooded by an expanding pro‐glacial, and partially supraglacial, lake. This case study indicates that thick sequences of debris‐charged basal ice/controlled moraine have a very low preservation potential but ice‐cored moraine complexes can develop into hummocky moraine belts in de‐glaciated terrains because they are related to the process of incremental stagnation, which at Kvíárjökull has involved periodic switches from transport‐dominant to ablation‐dominant conditions. Glacier recession is therefore recorded temporally and spatially by two suites of landforms relating to two phases of landform production which are likely typical for glaciers occupying overdeepenings: an early phase of active, temperate recession recorded by push moraines and lateral moraines and unconfined pro‐glacial meltwater drainage; and a later phase of incremental stagnation and pitted outwash head development initiated by the increasing topographic constraints of the latero‐frontal moraine arc and the increasing importance of the overdeepening as a depo‐centre. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

18.
This paper focuses on a topographic methodology to characterize the amount of sediment stored in channels and the use of historical photographs for aerial survey by stereophotogrammetry, as part of wider research on debris‐flow magnitude prediction. The topographic methodology uses equidistant four‐point cross‐sections along the long profile of the channel. Each cross‐section is representative of a 50‐m reach of the channel. To calculate the volume of each reach, the difference is calculated between a reference level and the topographic surface. The reference level is the lowest level where the debris flow can erode, and in the current method this level is estimated from fixed points along the long profile of the channel. The accuracy of the method has been estimated by comparing results of a detailed topographic survey, with a standard deviation corresponding to about 6 per cent of the total calculated sediment volume. This topographic methodology has been used on aerial photographs by photogrammetry. This tool was applied to photographs taken on 12 past dates. The scales of the archive photographs used range from 1:3000 to 1:30 000, but results are consistent and permit us to calculate sediment states of the channel for different past dates with an uncertainty of about 6 per cent of the total volume. The application of the technique to the Manival debris‐flow torrent has permitted us to propose some partial sediment budgets and erosion‐rate estimates. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

19.
20.
Coupling morphological, sedimentological, and rheological studies to numerical simulations is of primary interest in defining debris‐flow hazard on alluvial fans. In particular, numerical runout models must be carefully calibrated by morphological observations. This is particularly true in clay‐shale basins where hillslopes can provide a large quantity of poorly sorted solid materials to the torrent, and thus change both the mechanics of the debris flow and its runout distance. In this context, a study has been completed on the Faucon stream (southeastern French Alps), with the objectives of (1) defining morphological and sedimentological characteristics of torrential watersheds located in clay‐shales, and (2) evaluating through a case study the scouring potential of debris flows affecting a clay‐shale basin. Morphological surveys, grain‐size distributions and petrographic analyses of the debris‐flow deposits demonstrate the granular character of the flow during the first hectometre, and its muddy character from there to its terminus on the debris fan. These observations and laboratory tests suggest that the contributing areas along the channel have supplied the bulk of the flow material. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

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

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