首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 203 毫秒
1.
用旁侧声纳、浅地层剖面仪, 多普勒流速仪(ADP), 对长江中下游(武汉到河口段)河床进行了走航式的勘测工作, 获得了底床沙波形态的空间分布数据. 根据波高参数可将底床沙波分成4类(Ⅰ, Ⅱ, Ⅲ和Ⅳ), 并剖析了与水面坡降-过水面积的关系. 据沙波分布特征, 并结合100个断面沉积物采集数据, 以及水文实测和模拟参数, 可将河床分为5段(A, B, C, D和E段). A河段属中游蛇曲河型向下游分汊河型过渡, 受到下游B段—黄石基岩河床的影响, 堆积搬运明显, 大型沙波发育, 不对称性强; 基岩控制的B河段上方, 沙波由小变大, 但随着向下游能量释放, 迅速变小, 床底呈明显的冲刷搬运的特点. C河段是本区典型的分汊河型, 曲率高, 心滩发育. 随着上游能量的释放, 床底堆积搬运显著, 大型复合沙波发育, 波长可达150~260 m以上, 波高5~7 m, 上覆小型沙波, 反映季节性搬运的特点. D河段的沙波总体规模小于A~C, 对称性指数变小, 床底呈较弱堆积搬运特点, 并逐渐演变至河口区E段的、以小型沙波为主的河床; 受双向流的影响沙波对称性强. 研究还揭示了本区大型沙波多出现在水面坡降为0.2×10-5~2.0×10-5和过水面积为10×103~35×103 m2的河床区间, 底床沉积物多为中沙. 沿程模拟流速(流量为60000 m3·s-1下)表明形成A, B河段大沙波的流速需要2~3 m·s-1以上, 而下游C, D河段则只需要1.2~2.0 m·s-1. 研究成果为流域管理和环境工程研究提供了依据.  相似文献   

2.
链珠状沙波的发现及意义   总被引:1,自引:0,他引:1  
提出"沙波形态-伴生底形"的命名方法,将利用多波束测深系统对长江河口南北港分流口河槽床面形态测量时发现的由链状沙波和伴生底形-椭圆形凹坑组成的微地貌称为链珠状沙波.初步统计了链状沙波和椭圆形凹坑的几何参数,结果表明:链状沙波平均波高1.29 m,平均波长31.89 m,沙波指数(波长/波高)介于14~56.09,背流面常发育有次级沙波,每个沙波发育的椭圆形凹坑介于2~10个;椭圆形凹坑最大深度为1.98 m,最小深度为0.3 m,平均深度为0.98 m.对测区沙波发育环境进行了分析,主要得出:(1)在观测的一个涨潮和落潮过程中,落潮平均流速为0.78 m s~(-1),涨潮平均流速为0.27 m s~(-1),落潮历时大于涨潮历时;(2)底质粒径组成较复杂,粉砂(2~63mm)含量介于21.6%~23.4%,极细砂(63~125mm)含量介于28.2%~32.2%,细砂(125~250mm)含量在39.7%~41.6%;(3)沙波发育区水深在13~17 m.本研究丰富了沙波类型,为长江河口河槽动力地貌过程的研究提供了重要的参考资料,同时在河口工程应用中具有重要价值.  相似文献   

3.
三峡工程运用后,长江中游荆江河段持续冲刷,床沙与推移质、悬移质泥沙不断交换,从而造成该河段床沙发生不同程度的调整,对长江中下游河床演变及非平衡输沙机理的研究具有重要影响.在新水沙条件下,总结分析了沙波运动特性及床沙交换方式,引入Markov三态转移概率及非均匀沙隐暴系数,得到基于状态转移概率的沙质河段床沙级配调整的计算模型.研究结果表明:(1)20092014年,沙市站年内床沙中值粒径有先增大后减小的趋势,而监利站年内床沙中值粒径则先减小后增大,且荆江河段年际床沙中值粒径总体呈上升趋势,粗化程度约为6.9%~9.3%;(2)20092014年,沙市站床沙组成中粒径d<0.062 mm的泥沙所占比重不变,0.062 mm≤d<0.25 mm的泥沙所占比重逐年减少(累计减少11.4%),d≥0.25 mm的泥沙所占比重逐年增加(累计增加11.4%),而监利站床沙组成均存在波动性变化;(3)荆江河段床沙转换为推移质的概率随着泥沙粒径的增大而增大,床沙转换为悬移质的概率随着泥沙粒径的增大而减小,而推移质和悬移质转换为床沙的概率均随着泥沙粒径的增大而增大,河床发生冲刷粗化时泥沙输移的主要形式为悬移质(概率为81%~87%),而淤积细化时床沙补给主要来源于推移质(概率为8%~12%).通过验证,本文概率模型的计算结果与实测资料符合较好,能够应用于长江中游沙质河段年际床沙粗化及年内床沙级配调整过程预报,为进一步开展三峡工程下游非均匀悬移质泥沙沿程恢复机理的研究提供理论基础.  相似文献   

4.
三峡工程运用后长江中游荆江河段河床持续冲刷,局部河段崩岸频发,影响河道内悬沙输移与河床形态调整.本研究采用实测长程河道地形及固定断面地形资料,确定了2002-2018年荆江河段的主要崩岸区域,估算了崩岸土体的泥沙总量,进而定量分析了河岸崩退对河床调整的影响.计算结果表明:荆江段累计河岸崩退体积约为2.0亿m3,约占该河...  相似文献   

5.
受上游水库运行的影响,自2003年后长江中游河段河床调整,引起河道自身过流能力的变化.本文选取长江中游城陵矶-汉口河段作为研究对象,根据实测资料计算了2003-2016年两个水文断面的水位-流量关系及特征流量变化,并采用一维水动力学模型计算了河段尺度的平滩流量.结果表明:(1)螺山站及汉口站2003-2016年的水位-流量关系呈现枯水流量下水位降低、洪水流量下水位抬升的特点;(2)两站年最大流量下对应的水位整体抬升,警戒水位下对应的流量分别减小9%和16%;(3)2003-2016年城汉河段平滩流量有增有减,无明显单向变化趋势,2014年后呈减小趋势.分析断面的过流能力变化,发现螺山站与汉口站警戒流量与动床阻力呈明显的反比关系,过流能力由于河道阻力的增大而减小;床沙粗化、河道相对水深的减小及洲滩植被覆盖度的增大是引起河道阻力增大、过流能力减小的原因.  相似文献   

6.
系统解析长江中游河湖交汇区域航道水深资源的变化,可支撑“黄金航道”的可持续发展。本研究以洞庭湖-荆江交汇段为对象,研究洞庭湖分汇区域水沙条件、荆江河段滩槽演变与河湖交汇段航道水深资源的关系。研究表明:三峡工程运行以来,洞庭湖三口分流量和分沙量延续了三峡工程运行前的减少态势,伴随洞庭湖分流量减少,相对增加的长江干流径流量增强了河床冲刷强度,反馈使得洞庭湖三口分流量和分沙量均处于较低值;荆江河段河床冲刷给航道尺度提升奠定了有利基础,河湖分汇区域不满足4.5 m×200 m(水深×宽度)的长度为12.6 km,占荆江河段碍航总长度的68.35%。碍航驱动机制上:松滋口分流区段(枝城-昌门溪)的4.5 m水槽贯通但宽度不足200 m,汊道分流关系不稳定及洲滩萎缩制约航道条件稳定,枯水位下降及“坡陡流急”现象仍然严峻,不利于航道条件保障及船舶安全航行;太平口分流区段水位下降、洲滩萎缩、汊道交替发展使得枯水航路不稳定或水深不足4.5 m;藕池口分流区域的洲滩崩退、洲滩冲刷引起的向下游泥沙输移,碍航表现为航宽不足200 m或出现水深不足4.5 m的浅滩;洞庭湖入汇影响区段(熊家洲-城陵矶)受弯道冲淤...  相似文献   

7.
三峡工程运用后城陵矶-武汉河段河床调整及崩岸特点   总被引:1,自引:1,他引:0  
为研究近期城汉河段河床调整及崩岸特点,利用实测水沙及地形等资料,采用河段平均的方法,计算了城汉河段断面形态的调整过程,主要包括平滩河槽形态调整及其与前期水沙条件之间的关系.计算结果表明:城汉河段平滩河宽由2003年的1710 m增加至2016年的1732 m,增幅为1.28%,平滩水深由2003年的16.47 m增加至2016年的17.95 m,增幅为9.0%;白螺矶、界牌、簰洲及武汉等河段河床调整以纵向冲深为主,但陆溪口河段河床调整横向展宽与纵向冲深同步发展;2006-2016年城汉河段多年平均崩退速率为5.5 m/a,崩岸总长19.6 km,占岸线总长的8.3%,右岸占55.3%.簰洲河段岸线崩长占城汉河段岸线崩退总长的75.9%.此外还分析了河床边界与水沙条件等因素对重点河段(簰洲河段)崩岸过程的影响,来水来沙条件占主导地位,局部区域崩岸的发生依赖于河床边界条件;建立了典型断面平滩河宽与前期水沙条件之间的经验关系,较好地反映了水沙条件变化对崩岸过程的影响.  相似文献   

8.
北冰洋海冰厚度穿透雷达探测与下表面形态特征分析   总被引:10,自引:2,他引:10  
通过对北冰洋海冰穿透雷达的探测与分析, 讨论了雷达电磁波探测海冰厚度和海冰下表面形态特征的能力. 结果表明, 雷达波可穿透超过6 m厚度的北极盛夏海冰, 雷达波在海冰内的传播速度是0.142~0.154 m/ns. 雷达图像表现出海冰下表面粗糙起伏的微形态特征, 这种微形态特征与海冰类型关系紧密, 反映出雷达探测有助于识别和划分海冰类型. 基于雷达图像同时显示海冰上、下界线的特点, 统计分析了二者界线实际长度的差异, 提出一个新的长度因子数值概念, 可用于海冰上、下表面面积差异的研究.  相似文献   

9.
闫峰  王艳姣 《湖泊科学》2008,20(5):655-661
针对悬浮泥沙影响水体遥感测深精度的问题,选择长江口南港至南槽为研究区,通过对遥感测深方法研究,结合悬浮泥沙光谱特性分析,把"泥沙因子"引入到水体遥感测深反演模型中,研究表明:1)单因子非线性模型中,指数模型对0-2m的水深反演效果较好,对数模型对2-7m的水深反演较好,二次回归模型对7-14m的水深反演效果较好:2)建立的BP人工神经网络水深反演模型综合了多个波段具有的水深信息,模型的反演效果好于单因子非线性模型;3)实验构建的泥沙遥感参数综合了不同波段具有的悬沙信息,削弱了叶绿素和外界环境条件对泥沙信息的干扰,可较好地反映悬沙浓度变化特征;4)建立的BP人工神经网络泥沙因子水深反演模型削弱了悬浮泥沙对遥感测深的影响,模型实际反演能力明显优于单因子非线性模型和多因子BP人工神经网络水深反演模型.  相似文献   

10.
郭羽羽  李思悦  刘睿  张静 《湖泊科学》2021,33(3):737-748
为探究不同时空尺度土地利用与水质的相互关系,以黄河流域为研究对象,基于2018年河流水质与土地利用数据,提取3种空间尺度(1000 m河段缓冲区、500 m河岸带缓冲区、子流域)的景观格局指数,结合4个季节河流水质指标,采用冗余分析方法定量探讨土地利用对河流水质的多时空尺度影响,获得水质资源管理最有效的时空尺度.结果表明:(1)在3种空间尺度下,耕地和草地占比最大.斑块密度(PD)在1000 m河段缓冲区尺度下更大,香农多样性指数(SHDI)在3种空间尺度下变化不大;(2)水质指标季节性和空间性存在显著差异,夏季水质优于冬季;(3)土地利用在夏季和1000 m河段缓冲区尺度下对水质影响最显著;(4)不同景观格局指数对河流水质产生的影响不同,耕地面积占比和林地面积占比、未利用地面积占比和SHDI对水质影响最为显著,除了耕地面积占比和PD值的增加会导致水质下降外,其余土地利用指数在大多数时空尺度下均可改善水质.该研究基于不同的时空尺度对土地利用与水质的关系进行定量研究,为合理规划土地资源提供参考,从而保护流域水质安全.  相似文献   

11.
南海北部东沙海域巨型水下沙丘的分布及特征   总被引:3,自引:2,他引:1       下载免费PDF全文
本文基于多波束测深和高分辨率多道反射地震数据研究了东沙海域深水巨型水下沙丘的特征.巨型水下沙丘发育在230~830m水深的上陆坡范围内,呈斑块状分布.NW-SE向的近海底流体运动不仅冲蚀地层,形成了三条与水下沙丘间隔分布的冲蚀带,为水下沙丘提供了沉积物来源,同时也为水下沙丘的形成提供了动力源.研究区水下沙丘波长(L)范围55~510m,波高(h)范围1.5~20m,二者呈指数关系分布.沙丘的波长随水深增大而增大,波高则在500~700m水深范围内最大.水下沙丘NE—SW向展布的脊线和几何参数关系是与现今水动力条件相平衡的结果.  相似文献   

12.
Nearshore regions act as an interface between the terrestrial environment and deeper waters. As such, they play important roles in the dispersal of fluvial sediment and the transport of sand to and from the shoreline. This study focused on the nearshore of Poverty Bay, New Zealand, and the processes controlling the dispersal of sediment from the main source, the Waipaoa River. Hydrodynamics and sediment-transport in water shallower than 15 m were observed from April through mid-September 2006. This deployment afforded observations during 3-4 periods of elevated river discharge and 5 dry storms.Similar wind, river discharge, wave, current, and turbidity patterns were characterized during three of the wet storms. At the beginning of each event, winds blew shoreward, increasing wave heights to 2-3 m within Poverty Bay. As the cyclonic storms moved through the system the winds reversed direction and became seaward, reducing the local wave height and orbital velocity while river discharge remained elevated. At these times, high river discharge and relatively small waves enabled fluvially derived suspended sediment to deposit in shallow water. Altimetry measurements indicated that at least 7 cm was deposited at a 15 m deep site during a single discharge event. Turbidity and seabed observations showed this deposition to be removed, however, as large swell waves from the Southern Ocean triggered resuspension of the material within three weeks of deposition. Consequently, two periods of dispersal were associated with each discharge pulse, one coinciding with fluvial delivery, and a second driven by wave resuspension a few weeks later. These observations of nearfield sediment deposition contradict current hypotheses of very limited sediment deposition in shallow water offshore of small mountainous rivers when floods and high-energy, large wave and fast current, oceanic conditions coincide.Consistently shoreward near-bed currents, observed along the 10 m isobath of Poverty Bay, were attributed to a combination of estuarine circulation, Stokes drift, and wind driven upwelling. Velocities measured at the 15 m isobath, however, were directed more alongshore and diverged from those at the 10 m isobath. The divergence in the currents observed at the 10 and 15 m locations seemed to facilitate segregation of coarse and fine sediment, with sand transported near-bed toward the beach, while suspended silts and clays were exported to deeper water.  相似文献   

13.
We investigate how waves are transformed across a shore platform as this is a central question in rock coast geomorphology. We present results from deployment of three pressure transducers over four days, across a sloping, wide (~200 m) cliff‐backed shore platform in a macrotidal setting, in South Wales, United Kingdom. Cross‐shore variations in wave heights were evident under the predominantly low to moderate (significant wave height < 1.4 m) energy conditions measured. At the outer transducer 50 m from the seaward edge of the platform (163 m from the cliff) high tide water depths were 8+ m meaning that waves crossed the shore platform without breaking. At the mid‐platform position water depth was 5 m. Water depth at the inner transducer (6 m from the cliff platform junction) at high tide was 1.4 m. This shallow water depth forced wave breaking, thereby limiting wave heights on the inner platform. Maximum wave height at the middle and inner transducers were 2.41 and 2.39 m, respectively, and significant wave height 1.35 m and 1.34 m, respectively. Inner platform high tide wave heights were generally larger where energy was up to 335% greater than near the seaward edge where waves were smaller. Infragravity energy was less than 13% of the total energy spectra with energy in the swell, wind and capillary frequencies accounting for 87% of the total energy. Wave transformation is thus spatially variable and is strongly modulated by platform elevation and the tidal range. While shore platforms in microtidal environments have been shown to be highly dissipative, in this macro‐tidal setting up to 90% of the offshore wave energy reached the landward cliff at high tide, so that the shore platform cliff is much more reflective. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

14.
We present a two-dimensional vertical (2DV) flow and morphological numerical model describing the behaviour of offshore sand waves. The model contains the 2DV shallow water equations, with a free water surface and a general bed load formula. The water movement is coupled to the sediment transport equation by a seabed evolution equation. Using this model, we investigate the evolution of sand waves in a marine environment. As a result, we find sand wave saturation for heights of 10–30% of the average water depth on a timescale of decades. The stabilization mechanism, causing sand waves to saturate, is found to be based on the balance between the shear stress at the seabed and the principle that sediment is transported more easily downhill than uphill. The migration rate of the sand waves decreases slightly during their evolution. For a unidirectional steady flow the sand waves become asymmetrical in the horizontal direction and for a unidirectional block current asymmetrical in the vertical. A sensitivity analysis showed the slope effect of the sediment transport plays an important role herein. Furthermore, the magnitude of the resistance at the seabed and the eddy viscosity influence both the timescale and height of sand waves. The order of magnitudes found of the time and spatial scales coincide with observations made in the southern bight of the North Sea, Japan and Spain.  相似文献   

15.
地震波是一种随机的、不规则作用的动荷载脉冲,可分为振动型和冲击型。不同类型的地震波会对砂土液化和变形等产生重要影响,而传统的砂土震陷计算方法往往忽视这种因素,只考虑最大加速度幅值。通过编写UMAT子程序,在非线性有限元软件ABAQUS中开发亚塑性砂土边界面模型,对不同地震波类型下不同相对密度的砂土进行动单剪试验模拟,得到一系列砂土剪应变及竖向应变的时程曲线,并与室内试验结果进行对比分析。研究表明:在同一工况下,同类型的地震波引起的砂土竖向应变相近,不同类型引起的竖向应变差异明显;振动型地震波比冲击型引起的竖向应变更大。  相似文献   

16.
Tidal sand waves are dynamic bedforms found in coastal shelf seas. Moreover, these areas are inhabited by numerous benthic species, of which the spatial distribution is linked to the morphological structure of sand waves. In particular, the tube-building worm Lanice conchilegais of interest as this organism forms small mounds on the seabed, which provide shelter to other organisms. We investigate how the interactions between small-scale mounds (height ∼dm) and large-scale sand waves (height ∼m) shape the bed of the marine environment. To this end, we present a two-way coupled process-based model of sand waves and tube-building worm patches in Delft3D. The population density evolves according to a general law of logistic growth, with the bed shear stress controlling the carrying capacity. Worm patches are randomly seeded and the tubes are mimicked by small cylinders that influence flow and turbulence, thereby altering sediment dynamics. Model results relate the patches with the highest worm densities to the sand wave troughs, which qualitatively agrees with field observations. Furthermore, the L. conchilegatubes trigger the formation of sandy mounds on the seabed. Because of the population density distribution, the mounds in the troughs can be several centimetres higher than on the crests. Regarding sand wave morphology, the combination of patches and mounds are found to shorten the time-to-equilibrium. Also, if the initial bed comprised small sinusoidal sand waves, the equilibrium wave height decreased with a few decimetres compared to the situation without worm patches. As the timescale of mound formation (years) is shorter than that of sand wave evolution (decades), the mounds induce (and accelerate) sand wave growth on a similar spatial scale to the mounds. Initially, this leads to shorter sand waves than they would be in an abiotic environment. However, near equilibrium the wavelengths tend towards their abiotic counterparts again. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd  相似文献   

17.
Ebb-tidal deltas are highly dynamic environments affected by both waves and currents that approach the coast under various angles. Among other bedforms of various scales, these hydrodynamics create small-scale bedforms (ripples), which increase the bed roughness and will therefore affect hydrodynamics and sediment transport. In morphodynamic models, sediment transport predictions depend on the roughness height, but the accuracy of these predictors has not been tested for field conditions with strongly mixed (wave–current dominated) forcing. In this study, small-scale bedforms were observed in the field with a 3D Profiling Sonar at five locations on the Ameland ebb-tidal delta, the Netherlands. Hydrodynamic conditions ranged from wave dominated to current dominated, but were mixed most of the time. Small-scale ripples were found on all studied parts of the delta, superimposed on megaripples. Even though a large range of hydrodynamic conditions was encountered, the spatio-temporal variations in small-scale ripple dimensions were relatively small (height 0.015 m, length 0.11 m). Also, the ripples were always highly three-dimensional. These small dimensions are probably caused by the fact that the bed consists of relatively fine sediment. Five bedform height predictors were tested, but they all overestimated the ripple heights, partly because they were not created for small grain sizes. Furthermore, the predictors all have a strong dependence on wave- and current-related velocities, whereas the ripple heights measured here were only related to the near-bed orbital velocity. Therefore, ripple heights and lengths in wave–current-dominated, fine-grained coastal areas ( mm) may be best estimated by constant values rather than values dependent on the hydrodynamics. In the case of the Ameland ebb-tidal delta, these values were found to be m and m. ©2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.  相似文献   

18.

Saw-tooth bars are shore-oblique sand bars that are found along most ebb-tidal deltas of the Frisian Wadden Islands. Although they might significantly affect sediment transport pathways and volumes on the deltas, their general characteristics and dynamics are largely unknown. The main aim of this paper is to determine the typical saw-tooth bar heights, wave lengths, widths, orientations, migration rates and depths of occurrence. To this end, we analysed bathymetries from the Dutch and German Frisian Wadden Islands between 1970 and 2015. Bar heights range between 0.5 and 2 m, and wave lengths range between 300 and 900 m, with an average of 670 m, and the bar crests have a down contour orientation of approximately 25°. The bars are between 800 and 2200 m wide. Saw-tooth bars are located at depths from 3 to 12 m, depending on the slope of the area. Migration speeds of up to 112 m/year were found, using a spatial correlation method. Bar height and migration speed are positively linearly correlated, as well as bar wave length and bar width, bar height and the orientation of the depth contours and migration speed and the orientation of the area. The derived characteristics are compared to those of other bar types to evaluate hypotheses regarding the formation mechanism of the bars.

  相似文献   

19.
A comparison has been made between the hydraulic geometry of sand‐ and gravel‐bed rivers, based on data from alluvial rivers around the world. The results indicate a signi?cant difference in hydraulic geometry among sand‐ and gravel‐bed rivers with different channel patterns. On this basis, some diagrams for discrimination of meandering and braided channel patterns have been established. The relationships between channel width and water discharge, between channel depth and water discharge, between width–depth ratio and water discharge and between channel slope and water discharge can all be used for channel pattern discrimination. The relationship between channel width and channel depth can also be used for channel pattern discrimination. However, the accuracy of these relationships for channel pattern discrimination varies, and the depth–discharge relationship is a better discriminator of pattern type than the classic slope–discharge function. The cause for this difference has been explained qualitatively. To predict the development of channel patterns under different natural conditions, the pattern discriminator should be searched on the basis of independent or at least semi‐independent variables. The relationship between stream power and bed material grain size can be used to discriminate channel patterns, which shows a better result than the discriminator using the slope–discharge relationship. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

20.
基于多波束测深和高分辨率多道反射地震数据,以及遥感资料研究南海北部东沙海域深水沙丘,总结底形形态分布特征,分析其水动力成因.按照空间分布形态将沙丘分为平行沙丘、波状沙丘和网状沙丘,本文首次发现水下网状沙丘,重点讨论网状沙丘的特征分布和形成机制.网状沙丘发育在东沙环礁东北部水深约352~420 m之间,研究发现其延伸方向共有三个走向:NS、NE-SW和EW,前两个方向为西向传播的内波所致,而最后一个则很有可能是由洋盆西向传播的内波遇到东沙环礁后形成衍射,向北传播的内波引起了向南的近海底底流,冲刷海底底层,剥蚀出足够的沉积物颗粒,从而形成EW走向的水下沙丘.值得一提的是,在多边形网状沙丘的上方,反射地震剖面揭示水层中发育与下方沙丘相当波长的高频震荡内波,并且在遥感资料上也发现相关的衍射和内波干涉区域,而内波干涉很有可能导致高频震荡,从而形成网状沙丘.本研究丰富了水下沙丘形态类型,探讨了新型水下沙丘的形成机制,在陆坡底沙运动及工程应用中具有重要的价值.  相似文献   

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

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