首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Understanding the fundamental dynamics of interaction between multi-phase geophysical flows and engineering structures is crucial for mitigating geophysical hazards. Specifically, liquid phase between particles induces matric suction which could play a significant part in regulating flow dynamics and warrants further consideration. In this study, flume model tests were conducted to investigate the effects of water content (0–30%) on the impact behavior of granular flows. The particle image velocimetry technique was adopted to visualize the impact kinematics and the impact force was measured through a model barrier system. Results revealed that, besides geometric effects (kinetic sieving), mechanical effects (shearing and collision) are also vital for the mechanism of reverse segregation. At higher water contents, 20 and 30% in this study, discrete-surge impact, rather than a progressive impact process, was observed. The discrete surges induce impulses on the barrier. The discrete surges result from self-organization of unsaturated granular flows to overcome the enhanced shear strength induced by matric suction. Finally, a dimensionless index, namely the suction number, is used to quantify the effect of suction on the dynamic behavior of granular flows. Even for large-scale geophysical flows, if the content of fine particles is high, effect of suction should not be neglected.  相似文献   

2.
Flexible barriers undergo large deformation to extend the impact duration, and thereby reduce the impact load of geophysical flows. The performance of flexible barriers remains a crucial challenge because there currently lacks a comprehensive criterion for estimating impact load. In this study, a series of centrifuge tests were carried out to investigate different geophysical flow types impacting an instrumented flexible barrier. The geophysical flows modelled include covered in this study include flood, hyperconcentrated flow, debris flow, and dry debris avalanche. Results reveal that the relationship between the Froude number, Fr, and the pressure coefficient α strongly depends on the formation of static deposits called dead zones which induce static loads and whether a run-up or pile-up impact mechanism develops. Test results demonstrate that flexible barriers can attenuate peak impact loads of flood, hyperconcentrated flow, and debris flow by up to 50% compared to rigid barriers. Furthermore, flexible barriers attenuate the impact load of dry debris avalanche by enabling the dry debris to reach an active failure state through large deformation. Examination of the state of static debris deposits behind the barriers indicates that hyperconcentrated and debris flows are strongly influenced by whether excessive pore water pressures regulate the depositional process of particles during the impact process. This results in significant particle rearrangement and similar state of static debris behind rigid barrier and the deformed full-retention flexible barrier, and thus the static loads on both barriers converge.  相似文献   

3.
A quasi‐static homogeneous drained triaxial compression test on cohesionless sand under constant lateral pressure was simulated using a three‐dimensional discrete element method. Grains were modelled by means of particle clusters composed of rigid spheres or spheres with contact moments imitating irregular particle shapes. Attention was paid to the effect of initial void ratio and grain shape mixture on the shear strength, volume changes, force chains, kinetic, elastic and dissipated energies. In addition, the effect of the mean grain size, grain size distribution, grain size range, specimen size and roughness and stiffness of boundaries was numerically analysed in initially dense sand. Some numerical results were compared with available experimental results. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

4.
The onshore Georgina Basin in northern Australia is prospective for unconventional hydrocarbons; however, like many frontier basins, it is underexplored. A well-connected hydraulic fracture network has been shown to be essential for the extraction of resources from the tight reservoirs that categorise unconventional plays, as they allow for economic flows of fluid from the reservoir to the well. One of the fundamental scientific questions regarding hydraulic stimulation within the sub-surface of sedimentary basins is the degree to which local and regional tectonic stresses act as a primary control on fracture propagation. As such, an understanding of present-day stresses has become increasingly important to modern petroleum exploration and production, particularly when considering unconventional hydrocarbon reservoirs. This study characterises the regional stress regime in the Georgina Basin using existing well data. Wellbore geophysical logs, including electrical resistivity image logs, and well tests from 31 petroleum and stratigraphic wells have been used to derive stress magnitudes and constrain horizontal stress orientations. Borehole failure features interpreted from wellbore image and caliper logs yield a maximum horizontal stress orientation of 044°N. Integration of density log data results in a vertical stress gradient of 24.6 MPa km–1. Leak-off and mini-fracture tests suggest that this is the minimum principal stress, as leak-off values are generally shown to be at or above the magnitude of vertical stress. The maximum horizontal stress gradient is calculated to be in the range of 31.3–53.9 MPa km–1. As such, a compressional stress regime favouring reverse/reverse–strike-slip faulting is interpreted for the Georgina Basin.  相似文献   

5.
This paper reports an investigation on the performance of ethylene-vinyl acetate (EVA) foam when used as a cushion layer for rigid barriers used to resist debris flow. Large-scale pendulum impact tests were conducted to study the effects of layer thickness on cushion performance under six successive impacts. Results show that for the first impact at 70 kJ, the peak contact force with the EVA foam thickness of 500 mm is about twice larger compared to that of 1000 mm. Results also reveal that the cushion mechanism of elastic collapse of cell walls in the EVA foam results in peak contact forces and maximum transmitted loads that are up to 30 and 50% lower compared to gabions for the first impact at 70 kJ, respectively. The elastic behavior of EVA foam provides consistent cushioning efficiency. Furthermore, EVA foam is found to be susceptible to degradation by ultraviolet light so that a suitable coating layer is required for outdoor use. Polyurea was identified as a suitable coating material and a small-scale coating trial was performed to confirm this. Findings presented in this paper will have direct implications on the future design of cushion layers for rigid barriers used to intercept debris flows.  相似文献   

6.
The complexity of flow and wide variety of depositional processes operating in subaqueous density flows, combined with post‐depositional consolidation and soft‐sediment deformation, often make it difficult to interpret the characteristics of the original flow from the sedimentary record. This has led to considerable confusion of nomenclature in the literature. This paper attempts to clarify this situation by presenting a simple classification of sedimentary density flows, based on physical flow properties and grain‐support mechanisms, and briefly discusses the likely characteristics of the deposited sediments. Cohesive flows are commonly referred to as debris flows and mud flows and defined on the basis of sediment characteristics. The boundary between cohesive and non‐cohesive density flows (frictional flows) is poorly constrained, but dimensionless numbers may be of use to define flow thresholds. Frictional flows include a continuous series from sediment slides to turbidity currents. Subdivision of these flows is made on the basis of the dominant particle‐support mechanisms, which include matrix strength (in cohesive flows), buoyancy, pore pressure, grain‐to‐grain interaction (causing dispersive pressure), Reynolds stresses (turbulence) and bed support (particles moved on the stationary bed). The dominant particle‐support mechanism depends upon flow conditions, particle concentration, grain‐size distribution and particle type. In hyperconcentrated density flows, very high sediment concentrations (>25 volume%) make particle interactions of major importance. The difference between hyperconcentrated density flows and cohesive flows is that the former are friction dominated. With decreasing sediment concentration, vertical particle sorting can result from differential settling, and flows in which this can occur are termed concentrated density flows. The boundary between hyperconcentrated and concentrated density flows is defined by a change in particle behaviour, such that denser or larger grains are no longer fully supported by grain interaction, thus allowing coarse‐grain tail (or dense‐grain tail) normal grading. The concentration at which this change occurs depends on particle size, sorting, composition and relative density, so that a single threshold concentration cannot be defined. Concentrated density flows may be highly erosive and subsequently deposit complete or incomplete Lowe and Bouma sequences. Conversely, hydroplaning at the base of debris flows, and possibly also in some hyperconcentrated flows, may reduce the fluid drag, thus allowing high flow velocities while preventing large‐scale erosion. Flows with concentrations <9% by volume are true turbidity flows (sensu 4 ), in which fluid turbulence is the main particle‐support mechanism. Turbidity flows and concentrated density flows can be subdivided on the basis of flow duration into instantaneous surges, longer duration surge‐like flows and quasi‐steady currents. Flow duration is shown to control the nature of the resulting deposits. Surge‐like turbidity currents tend to produce classical Bouma sequences, whose nature at any one site depends on factors such as flow size, sediment type and proximity to source. In contrast, quasi‐steady turbidity currents, generated by hyperpycnal river effluent, can deposit coarsening‐up units capped by fining‐up units (because of waxing and waning conditions respectively) and may also include thick units of uniform character (resulting from prolonged periods of near‐steady conditions). Any flow type may progressively change character along the transport path, with transformation primarily resulting from reductions in sediment concentration through progressive entrainment of surrounding fluid and/or sediment deposition. The rate of fluid entrainment, and consequently flow transformation, is dependent on factors including slope gradient, lateral confinement, bed roughness, flow thickness and water depth. Flows with high and low sediment concentrations may co‐exist in one transport event because of downflow transformations, flow stratification or shear layer development of the mixing interface with the overlying water (mixing cloud formation). Deposits of an individual flow event at one site may therefore form from a succession of different flow types, and this introduces considerable complexity into classifying the flow event or component flow types from the deposits.  相似文献   

7.
Consolidation tests have been performed to determine the mechanical properties of crushed salt as affected by applied stresses and consolidation period. The crushed salt with particle sizes ranging from 0.075 to 4.76 mm mixed with 5% saturated brine are consolidated under axial stresses ranging from 2.5 to 10 MPa. The consolidation curves show instantaneous and transient strains. The densities, uniaxial compressive strengths and elastic moduli measured after consolidation for 3–180 days increase with applied stresses and time. The crushed salt properties are determined as a function of mean strain energy required during consolidation. The relations are used to predict the crushed salt properties installed in exploratory boreholes under various external pressures. Two main mechanisms strengthen and stiffen the crushed salt mass: (1) consolidation by volumetric change due to particle rearrangement and cracking, and (2) recrystallization and healing processes. The borehole depth and installation time are significant factors controlling the long-term density, strength and stiffness of the crushed salt seal.  相似文献   

8.
Wang  Dong-Wei  Zhu  Cheng  Tang  Chao-Sheng  Li  Sheng-Jie  Cheng  Qing  Pan  Xiao-Hua  Shi  Bin 《Acta Geotechnica》2021,16(9):2759-2773

Deep geological repository is a favorable choice for the long-term disposal of nuclear wastes. Bentonite–sand mixtures have been proposed as the potential engineered barrier materials because of their suitable swelling properties and good ability to seal under hydrated repository conditions. To investigate the effects of sand grain size on the engineering performance of bentonite–sand mixtures, we prepare five types of bentonite–sand mixtures by mixing bentonite with sand of varying particle size ranges (0.075–0.25 mm, 0.25–0.5 mm, 0.5–1 mm, 1–2 mm and 2–5 mm, respectively). We carry out sequential oedometer tests under different simulated repository conditions, including constant vertical stress (CVS), constant stiffness (CS) and constant volume (CV) conditions. The microstructural heterogeneity and anisotropy of these soil mixtures are characterized through the quantitative analysis of micro-CT scanning results. Experimental results reveal that both sand grain size and boundary condition significantly influence the swelling of soil mixtures. Under three conditions, the temporal evolutions of swelling stress and strain follow similar trends that they increase faster at the beginning and gradually stabilize afterward. Comparing the ultimate values, swelling strains follow CVS?>?CS?>?CV, while swelling stresses follow CV?>?CS?>?CVS. Under CS boundary conditions, as the stiffness coefficient increases, the swelling pressure increases and the swelling strain decreases. CT results further indicate that mixtures with larger sand inclusions are more structurally heterogeneous and anisotropic, resulting in increased inter-particle friction and collision and a higher energy dissipation during the swelling process. Moreover, the non-uniform distribution of bentonite in local zones would be intensified, which plays an important role in compromising swelling behavior. Therefore, soil samples mixed with larger sand particles present a smaller swelling stress and strain values. This study may guide the choice of engineered barrier materials toward an improved design and assessment of geological repository facilities.

  相似文献   

9.
The aim of the study was to investigate the effect of rigid second phases on grain growth of a matrix phase. For this purpose, variable mixtures of norcamphor as the matrix phase, with glass beads (0.08–0.51 volume fraction) as second phase, were used to perform see-through rock-analogue experiments under static conditions at constant temperatures (50°C). Irrespective of the second-phase content, grain-size evolution of all mixtures can be subdivided into a stage of continuous grain growth, a transient stage and a stage of a finally stabilized grain size. On the grain-scale, the second phases affect the migrating grain boundaries either by pinning by single particles, by multiple particles or even by particle clusters. Summed up over the entire aggregate, these pinning regimes affect the average bulk grain size of the matrix grains, such that the changes in matrix grain size directly correlate with the amount of second phases, their dispersion and their degree of clustering. In this way, the matrix grain size decreases with increasing second-phase content, which can be expressed as a Zener relationship. Originating from the modification of an ordinary grain growth law, a new mathematical expression is defined, which allows the calculation of changes in the matrix grain size as a function of different second-phase volume fractions and particle sizes. Such models will be helpful in the future to predict microstructural changes in polymineralic rocks at depth.  相似文献   

10.
The 2008 Wenchuan earthquake induced a large number of landslides, and a vast amount of loose landslide materials deposited on steep hill slopes or in channels. Such loose materials can become sources of deadly debris flows once triggered by storms. On 13 August 2010, a storm swept Yingxiu and its vicinity, triggering a catastrophic debris flow with a volume of 1.17?million?m3 in Xiaojiagou Ravine. The debris flow buried 1,100?m of road, blocked a river and formed a debris flow barrier lake. A detailed field study was conducted to understand the initiation mechanisms and runout characteristics of this debris flow. Two types of debris flows are identified, namely hill-slope debris flow and channelized debris flow. The hill-slope debris flow was triggered in the forms of firehose effect, rilling and landsliding, whereas the channelized debris flow was triggered in the form of channel-bed failure. This debris flow was a water?Crock flow since most particles were gravel, cobble or larger rocks and the fraction of silt and clay was less than 2%. Grain contact friction, pore-pressure effects and inertial grain collision were the three most important physical interactions within the debris flow. Such interactions yielded a smaller runout distance (593?m) compared with those of mud?Crock flows of similar size.  相似文献   

11.
The elastic constants of a crystal under stress, defined as the second derivative of the crystal free energy with respect to strain, require a correction related to the static pressure at non-zero pressures. The corrections required for the elastic constants calculated by the free energy minimisation code PARAPOCS are described and tested by comparison with the elastic constants calculated numerically by applying small stresses in the appropriate orientations to simulated crystals of fluorite, forsterite, α-quartz and albite. The corrected elastic constants are then used to investigate the extrapolation of the bulk and shear moduli (and hence also the seismic wave velocities V p and V s) of β-spinel and forsterite to upper mantle pressures. A Murnaghan equation, thirdorder Eulerian finite strain equation, second order polynomial equation and a logistic equation were all fitted to the simulated bulk and shear moduli between 0 and 3 GPa pressure. The parameters derived for these equations are used to extrapolate the bulk and shear moduli to 14 GPa and the results are compared to the simulated high pressure moduli. Over this pressure range, the second order polynomial provides the best extrapolation of the bulk modulus, but the use of the logistic equation results in the best extrapolation of the shear modulus.  相似文献   

12.
The principal numerical approach to describing sedimentation from multi-component suspensions (Mirza and Richardson, 1979) has been applied only to systems containing two or three particle size populations, but is theoretically applicable to suspensions containing a wider range of particle sizes. In order to adapt this model to the simulation of sedimentation from natural high-density sediment-laden flows, we have conducted computational tests on the sedimentation of suspensions with up to ten particle size populations. The tests run smoothly with binary systems and sometimes ternary systems, but fail with systems containing more than three particle size populations. One cause of the discrepancy between theoretical predictions and computational tests arises from the method used to calculate the changing concentration of particles in each zone of a stratified settling suspension. In quantifying the changing sediment budget for each sedimentation zone, the equation set fails to include the apparent outflow of finer particle populations across the upper boundary of the zone. In the present study, we express the sediment budget of each particle population within each sedimentation zone as the net sediment flux, the algebraic sum of the apparent particle inflow from the zone's lower boundary and the apparent particle outflow across the zone's upper boundary. This revised model successfully predicts the evolution of multi-component suspensions containing up to ten particle size populations over the length of time required for complete sedimentation, up to 8×105 seconds. It can be used to predict the sediment sorting and the vertical textural variation of beds formed by the simulated sedimentation of a multi-component suspension. The model provides the basis for future computer simulation of sedimentation from highly concentrated sediment flows and for the prediction of downslope textural and structural variations of turbidites.  相似文献   

13.
Interaction between dry granular flow and deflectors   总被引：2，自引：0，他引：2  
  相似文献   

14.
Strength and deformation behaviors of rockfill materials,key factors for determining the stability of dams,pertain strongly to the grain crushing characteristics.In this study,single-particle crushing tests were carried out on rockfill materials with nominal particle diameters of 2.5 mm,5 mm and 10 mm to investigate the particle size effect on the single-particle strength and the relationship between the characteristic stress and probability of non-failure.Test data were found to be described by the Weibull distribution with the Weibull modulus of 3.24.Assemblies with uniform nominal grains were then subjected to one-dimensional compression tests at eight levels of vertical stress with a maximum of 100 MPa.The yield stress in one-dimensional compression tests increased with decreasing the particle size,which could be estimated from the single-particle crushing tests.The void ratio-vertical stress curve could be predicted by an exponential function.The particle size distribution curve increased obviously with applied stresses less than 16 MPa and gradually reached the ultimate fractal grading.The relative breakage index became constant with stress up to 64 MPa and was obtained from the ultimate grading at the fractal dimension(a?2:7).A hyperbolical function was also found useful for describing the relationship between the relative breakage index and input work during one-dimensional compression tests.  相似文献   

15.
In this paper the effects of maximum particle size, particle gradation/sorting and fabric on bulk mechanical behaviour of granular materials such as coarse grained soils and rockfills are investigated" from micromechanical considerations starting from the grain scale level, using numerical" simulations based on Discrete Element Modelling (DEM). Hydrostatic compaction and biaxial tests on 2-dimensional assemblies of discs with varying particle sizes and gradations were modelled using DEM. An examination of the constitutive behaviour of granular media considering" the particulate nature of the medium has been attempted to explain the effect of particle size and gradation. Simulation results on perfectly parallel graded assemblies indicate that with increase in the size of the particles, a marginal increase (or no increase) in the angle of internal friction is observed during biaxial loading conditions. A change to a wider gradation (keeping the minimum grain size the same) results in a decrease in the angle of internal friction and an increase in volumetric strain to a considerable extent. Based on micromechanical force and fabric parameters, the basis for the physical behaviour was established. This helps in understanding the physics of parallel gradation techniques.  相似文献   

16.
流态化运动是高速远程滑坡的主要运动形式,是揭示高速远程滑坡运动机理的重要基础。基于粒子图像测速（PIV）分析方法,采用物理模型试验对不同粒径组成条件下的颗粒流内部的速度分布、剪切变形及流态特征进行了研究,并对高速远程滑坡流态化运动特征进行了讨论分析。结果表明:碎屑流流态化运动特征与颗粒粒径呈显著的相关性,随着粒径的减小或细颗粒含量的增加,颗粒流底部相对于边界的滑动速度以及整体的运动速度均呈逐渐减小的趋势,颗粒流内部剪切变形程度增加,颗粒的运动形式由“滑动”向“流动”转变;当颗粒粒径较小或细颗粒含量较高时,颗粒流内部剪切速率增大的趋势在颗粒流底部更加显著,反映了粒径减小有助于促进颗粒流内部剪切向底部的集中;在同一颗粒流的不同运动阶段及不同纵向深度,其流态特征具有显著差别,颗粒流前缘及尾部主要呈惯性态,颗粒间以碰撞作用为主,而主体部分则主要呈密集态,颗粒间以摩擦接触作用为主;在颗粒流表面及底部,颗粒间相互作用方式主要是碰撞作用,中间部分则以摩擦作用为主;对于不同粒径的颗粒流,随着粒径的增大或粗颗粒含量的增加,颗粒流内部颗粒的碰撞作用加强,颗粒流整体趋于向惯性态转变。  相似文献   

17.
We established an equation of state for nanocrystalline forsterite using multi-anvil press and diamond anvil cell. Comparative high-pressure and high-temperature experiments have been performed up to 9.6 GPa and 1,300°C. We found that nanocrystalline forsterite is more compressible than macro-powder forsterite. The bulk modulus of nanocrystalline forsterite is equal to 123.3 (±3.4) GPa whereas the bulk modulus of macro-powder forsterite is equal to 129.6 (±3.2) GPa. This difference is attributed to a weakening of the elastic properties of grain boundary and triple junction and their significant contribution in nanocrystalline sample compare to the bulk counterpart. The bulk modulus at zero pressure of forsterite grain boundary was determined to be 83.5 GPa.  相似文献   

18.
The Udo tuff cone of Cheju Island, South Korea, is a middle Pleistocene basalt tuff cone that has formed by early Surtseyan-type eruptions and later drier hydroclastic eruptions. The tuff cone comprises steep (20–30°) and planar beds of lapillistone, lapilli tuff and tuff that can be grouped into seven sedimentary facies (A-G). Facies A and B comprise continuous to lenticular layers of grain-supported and openwork lapillistone that are inversely graded and coarsen downslope. They suggest emplacement by grain flows that are maintained by gravity-induced stress and grain collisions. Facies C includes poorly sorted, crudely bedded and locally inversely graded lapilli tuff, also suggestive of rapid deposition from highly concentrated grain flows. Facies D includes thinly stratified and mantle-bedded tuff that was probably deposited by fallout of wind-borne ash. Other facies include massive lapilli tuff (Facies E), chaotic lapilli tuff (Facies F) and cross-bedded tuffaceous sandstone (Facies G) that were deposited by resedimentation processes such as debris flow, slide/slump and stream flow, respectively. The grain flows that produced Facies A, B and C are interpreted to have originated from falling pyroclasts, which initially generated highly dispersed, saltating avalanches, in which momentum was transferred by the particles themselves. This transport mechanism is similar to that of debris fall. As the slope gradient was too low to maintain a highly dispersed flow, the debris fall decelerated and contracted due to a decrease in dispersive pressure. The mode of momentum transfer changed to one of collision because contraction of the debris fall resulted in an increase in particle concentration. This transport mechanism is similar to that of common grain flows. Grain segregation occurred in several ways. Initial segregation of ash from lapilli occurred due to their differing terminal fall velocities, and their contrasting degrees of sliding friction with the bed. Percolation of ash into interstices of lapilli during flow (kinematic sieving) augmented further segregation of ash from lapilli. The latter process, along with a dispersive pressure effect, gave rise to vertical inverse size grading. Downdip inverse grading was produced by particle overpassing.  相似文献   

19.
岩石作为矿物颗粒的集合体,矿物粒径非均质性对其宏观力学特性影响比较明显。基于颗粒流程序PFC2D,通过设置不同种类粒径组合及粒径比来体现粒径非均质性,研究了粒径非均质性对岩石材料宏观力学特性（弹性模量、峰值强度、泊松比）的影响。研究中设计了6种粒径组合方案,粒径种类数分别为:连续粒径、10种、8种、5种、3种、2种,每种方案下设置5种平均粒径及4种粒径比,进行单轴压缩试验。结果表明,岩石内部存在颗粒尺寸效应和粒径非均质效应,岩石弹性模量和峰值强度随粒径增大均呈减小的趋势,随粒径非均质性的提高整体上也呈减小的趋势,但局部变化阶段受模型中细颗粒含量及数量的影响会呈增大的趋势。粒径对弹性模量的作用机制主要是通过影响模型孔隙率实现的。研究结果揭示了岩石宏观特性的变化是模型内颗粒尺寸效应和粒径非均质性效应共同作用的结果,为掌握矿物粒径对岩石强度及变形特性的影响提供了一定依据。  相似文献   

20.
Collision data are presented from coloured high-speed films of three size fractions of sand grains saltating over a bed of the total grain population. Each fraction was colour tagged and the proportion of each size ejected by grains colliding with the surface was recorded on a number of films taken as the bed was progressively eroded. The results confirm earlier findings that V3/V1?0.5–0.6, Vn/V1?.08 and the rebound angle increases with decreasing grain size. Ejected grains are examined in relation to their size, the impactor size, ejection speed and angle and the number of ejecta per collision. In addition, changes in grain parameters are observed with time. For fine impactors, ejection speeds generally increase with a decrease in ejecta size, but the fine fraction does not follow this trend for the coarse and medium impactors. Ejection angles are typically between 40° and 60°, with coarse grains having shallower mean angles than fine ejecta. The number of ejections per collision increases with a decrease in particle size for each impactor size. The general tendency for coarse particles to be ejected at lower speeds and shallower angles than fine particles will lead to sorting of the grain sizes. There is poor correlation between the forward momentum loss of the saltating grams at collision and both the forward momentum of the ejected grains and the number of ejected grains. Much of the forward momentum of the saltating grains is transfered to creeping grains. The composition and geometry of the bed are considered to be important factors in the evolution of the saltation cloud.  相似文献   

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

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