首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 129 毫秒
1.
近30年来南四湖流域城市化进程中的水系变化分析   总被引:1,自引:1,他引:0       下载免费PDF全文
城市化对水系演化影响的研究国内目前多集中在快速、高速城市化地区,而对大流域、城市化发展较缓地区的研究比较薄弱.以南四湖流域为研究区,基于1987、2000和2014年3期遥感影像,分析了流域城市化进程中的下垫面变化特征;选取流域1980s、2003和2014年的地形图进行水系提取,从数量参数、结构参数和连通性参数3个角度分析近30年城市化进程中水系结构的时空变化特征.结果表明:(1)近30年来流域建设用地增加了1568.06 km~2,2000年以后城市建设用地扩张显著,2012年流域人口城市化率为32%;(2)1980s—2010s流域总河流长度、面积和河网密度均呈现出持续减少趋势,分别减少了135.46 km、2.75 km~2和0.49 km/km~2,各级河流表现出不同的变化特点,较低等级河流受到的影响较大;而流域水面率持续增加,近30年共增加了59.79%;(3)流域水系总体上还保持着自然状态下的空间格局,但结构特征发生了较大改变,河网结构稳定度减少了4.30%,连接率和实际结合度分别减少了21.82%和21.62%;子流域内部距湖区越远的空间城市扩展强度指数值越大,城市化对水系的影响越显著.该研究将补充对不同空间尺度、不同城市化水平地区河网水系演化影响的案例,并为研究区河网水系的保护提供支持与参考.  相似文献   

2.
丽江-小金河断裂位于青藏高原东南缘,是川滇菱形块体内重要的次级边界断裂。构造活动对区域河流水系的发育有重要影响,因此,定量研究水系地貌特征可在一定程度上反映区域内构造活动信息。本文采用地形坡度、地形起伏度、面积-高程积分值(HI)等地貌参数研究沿该断裂区域地貌对构造活动的响应。研究发现,断裂北侧坡度和地形起伏度较断裂南侧高,跨断裂4级流域盆地内亚流域盆地面积-高程积分值变化也指示断裂北侧流域成熟度低于断裂南侧,这反映了丽江-小金河断裂存在一定逆冲分量。  相似文献   

3.
为分析海拔高度大于1 300m的中高山地区地形地貌因子对滑坡发育的敏感性程度,运用基于因子综合贡献率的敏感性分析方法,研究黔东北中高山地区坡面形态、地形坡度、相对高差三类影响因子对滑坡发育的影响程度.在各类影响因子对滑坡发育频数、面积及体积分量贡献率的基础上,利用顺序赋值法量化各分量贡献率的贡献指数,并以各分量贡献指数的几何平均数构建敏感性系数,定量分析三类地形地貌影响因子对中高山地区滑坡发育的敏感性.研究结果表明,阶梯形和凸形坡面形态的斜坡中后部为地貌演化堆积区,能提供较大重力势能,属滑坡发育的敏感性坡面形态因子;地形坡度在18°~22°之间的较陡坡和相对高差在40~80m之间的矮坡常分布于受地表水强烈冲刷的地带,在一定重力势能和强烈地下水共同作用下,极易诱发浅层小规模滑坡;地形坡度大于27°的极陡坡和相对高差大于120m的高坡由于具备较大重力势能,常诱发深层大规模滑坡.  相似文献   

4.
地震台网的震中定位精度判据   总被引:2,自引:0,他引:2  
用已出版的地震报告中的数据对震中精度进行了可靠且保守的估计。因为大量的地震研究依赖于地震定位报告,所以要估算地震定位误差。地震定位和参数估算大多数都是通过传统的线性反演方法用一维地球模型计算走时来获得。假设误差是高斯分布、零均值和不相关的,报告给出的定位结果是不精确的。遗憾的是,这种假设通常情况下是不成立的,特别对高置信度而言,不能真实估算定位误差。研究发现,地震定位精度与地震台网的几何形状紧密相关。我们利用两个确定震中位置的爆破形成近震台网(0°~2·5°)的定位判据。利用蒙特卡罗模拟台网几何,我们发现当近震台网符合以下判据时近震台网定位能精确到5km之内,具有95%的置信度:(1)250km以内有10个或更多的台站;(2)方位角小于110°;(3)次级方位角小于160°;(4)30km以内至少有一个台站。为得到近区域台网(2·5°~10°)、区域台网(2·5°~20°)和远震台网(28°~91°)的地震定位精度判据,我们使用了大量定位很好的地震和核爆炸数据。当超过局部范围时,我们发现次级方位角足以控制定位精度。当次级方位角超过120°时,定位误差将增加。当台站覆盖的次级方位角标准小于120°时,近区域台网的定位精度为20km,具有90%的置信度,区域台网和远震台网的定位精度为25km,具有90%的置信度。  相似文献   

5.
曾冰茹  李云良  谭志强 《湖泊科学》2023,35(5):1796-1807
由于气候变化和人类活动等多重影响,流域河湖水系格局与连通程度发生了显著变化,进而引发洪涝灾害等一系列水资源问题。本文以鄱阳湖流域为研究区,基于Google Earth Engine(GEE)提取1989—2020年5期水系数据,采用图论方法构建水系评价体系,定量分析该地区近30年来水系格局和结构连通性的时空演变特征,并结合该时期地形、土地利用和归一化植被指数(NDVI)等数据,利用连通性指数(index of connectivity,IC)评估功能连通性的动态变化,进而探讨水文连通与径流量和输沙量的联系。结果表明,近30年来鄱阳湖流域水系结构趋于复杂化,主要体现在流域北部。除干流外,其他等级河流的数量和长度均有所增加,其中Ⅲ级河流最为明显。河网密度、水面率、河网复杂度和发育系数均呈增加趋势,2000年后的变化率约为2000年前的两倍。水系连通环度、节点连接率和水系连通度总体增加,结构连通性呈好转趋势且变化幅度较小。功能连通分析表明,近30年来大部分流域IC减少,流域下游靠近主河道的平坦地区IC较高,上游远离河道的植被密集区域IC较低。此外,IC与年径流量和输沙量表现为显著的正相关性(...  相似文献   

6.
新疆天山地区壳幔S波速度结构特征及变形分析   总被引:1,自引:0,他引:1       下载免费PDF全文
天山地区地质构造复杂,地震活动频繁,其壳幔变形和深部结构一直受到学者们的高度关注.然而,由于天山地区地震台站资料较少,致使壳幔变形研究结果与解释存在诸多争议.本研究利用在天山地区(40°N-46°N,78°E-92°E)新布设的11个流动宽频带地震台站和该地区39个固定台站的观测资料,采用接收函数与面波联合反演方法,获得了研究区地壳厚度及壳幔S波速度结构.反演结果显示天山地区(41.5°N-44°N,78°E-88°E)平均地壳厚度为56 km,塔里木盆地(40°N-41.5°N,79°E-90°E)、准噶尔盆地(44°N-46°N,82°E-90°E)和吐鲁番盆地(42°N-43°N,88°E-90°E)具有较厚的沉积层,地壳平均厚度为43 km、53 km和46 km,整体表现为天山厚、盆地相对较薄的特征;在研究区南天山的最高峰(42°N,80.5°E)及北天山的最高峰(43.5°N,86°E)附近,中下地壳存在较厚的低速层,我们认为在强烈挤压作用下低速、低强度的中下地壳强烈变形可能是导致该区域快速隆升的主要原因.在研究区中部,位于塔里木盆地与准噶尔盆地之间的天山地区,中下地壳及上地幔均存在低速层,且盆地莫霍面向天山倾斜明显.结合前人的研究成果推测,在南北向构造挤压应力作用下,塔里木盆地与准噶尔盆地发生了向天山造山带方向的双向壳幔层间插入俯冲.在研究区东部,塔里木盆地东北缘与天山东部接触带的地壳内没有明显的低速层,推测应处在早期挤压变形状态,该区域的壳幔边界为缓变的速度梯度带,可能与上地幔热物质侵入或渗透有关.  相似文献   

7.
太湖流域城市化对平原河网水系结构与连通性影响   总被引:1,自引:1,他引:0  
太湖流域的快速城市化使河流形态发生了改变,影响了平原河网水系结构与连通性本文以近几年城市化发展较为迅速的张家港地区为例,基于景观生态学和图论的相关理论,利用2002与2015年的水系数据,对全区以及各行政分区水面率、河网密度、河频率等指标进行研究,旨在揭示该地区河网水系结构与连通性参数变化及城市化对其影响结果表明:研究区2015年与2002年相比全区各水系结构指标都呈减小趋势,水面率、河网密度、河频率的衰减率分别为15.7%、12.6%、33.3%,河网有主干化趋势,干流面积与长度发育的同步性较好;不同行政区的水系变化有明显的空间分异,张家港南部水系较北部水面率、河频率衰减更剧烈;东北部地区河网密度与河网复杂度变化明显;张家港全区河网连通性略有提高,但地区分异明显,沿江的西北部地区河网连通能力有所加强,而东南部内河区域河网畅通程度减弱,被圩区、水闸与泵站切断的水系需要更合理的调度方案才能增强并维持水体的连通度;水面率、河网密度与人口、GDP之间相关水平较高,各地区水系结构参数与其变化率的空间分异和该地区人口、经济发展水平密切相关;河网河链数、节点数受城市化影响较大,线点率、实际连通度的变化与城市化指标相关性较弱本研究可以为后期张家港地区优化水系结构及防洪排涝建设提供基础,并帮助基于各片区的特点开展因地制宜的管理.  相似文献   

8.
在收集了沈阳及其附近地区 (12 1°30′~ 12 4°30′E ,4 1°~ 4 3°N)已有地震震源机制解成果的基础上 ,利用MAPSIS系统提供的P波初动求解方法 ,求解了 1974年 12月~ 2 0 0 1年 12月期间 ,该区域小震复合震源机制解 (平均解 )。结果表明 ,研究区主压应力轴P轴集中分布在NEE SWW方向 ,主张应力轴T轴集中分布在NNW SSE方向 ,这与本地区区域应力场的方向一致。  相似文献   

9.
东北地区重力均衡异常特征的初步研究   总被引:3,自引:0,他引:3       下载免费PDF全文
本文利用地面实测重力资料和地形高程资料,采用普拉特-海福特(Pratt-ttayford)重力均衡理论模型,取1°×1°方格网,通过使用现成改正表格查取改正值与个别计算点用理论公式计算作校核的方法,计算了我国东北地区75个计算点的均衡重力异常值;并对局部第四系覆盖较厚地区作了第四系密度改正;在此基础上,构制了我国东北N39°—49°,E121°—131°大部分地区的均衡重力异常图;结合区域布格重力异常和区域空间重力异常特征以及莫霍界面的起伏特点作了对比分析和讨论  相似文献   

10.
杨素  万荣荣  李冰 《湖泊科学》2022,34(4):1055-1074
水文连通性作为连接河湖水体间物质、能量及信息传递与交换的关键纽带,对水环境、水生态和生境状况具有联动与触发反馈作用,已成为水文学、水利学和生态学等诸多领域的研究热点.太湖流域作为全国典型的流域性密集平原河网区,在快速城镇化背景下,河湖水文连通关系变化剧烈并引起了一系列生态环境效应.本文以水文连通性概念与内涵为背景,对太湖流域水文连通性研究进行了综述.太湖流域水文连通性评价方法以图论法、景观格局指数及水文连通性函数等方法联合使用为主,且聚焦于区域尺度研究;演变过程按人类活动影响强度大致划分为自然水系形成阶段、古代人类活动影响阶段和现代人类活动影响阶段;气候变化和人类活动共同影响着太湖流域水文连通性变化,近百年来水利工程建设和城市化进程等人类活动的影响尤为剧烈;良好的水文连通性有利于汛期减轻流域洪涝灾害及非汛期保障水资源供给,但水文连通性的提高对水环境和水生态的效应由于涉及因素众多尚存在争议.针对当前研究现状和存在问题,提出(1)平原河网区水文连通性的定量表征与评估是水文连通性研究的前提;(2)定量解析流域水文连通性的驱动机制是水文连通性研究的重点;(3)深入跟踪大型引水调水工程对流域水文连通影响及其效应是一项长期任务;(4)兼顾环境生态效应的水利工程生态化改造研究,开展工程控制背景下的流域水文连通多目标优化调控,是实现太湖流域洪水调蓄、水资源供给、水环境净化、生物多样性维持等生态系统服务协同提升的重要途径.  相似文献   

11.
Stream burning is a common flow enforcement technique used to correct surface drainage patterns derived from digital elevation models (DEM). The technique involves adjusting the elevations of grid cells that are coincident with the features of a vector hydrography layer. This paper focuses on the problematic issues with common stream burning practices, particularly the topological errors resulting from the mismatched scales of the hydrography and DEM data sets. A novel alternative stream burning method is described and tested using five DEMs of varying resolutions (1 to 30 arc‐seconds) for an extensive area of southwestern Ontario, Canada. This TopologicalBreachBurn method uses total upstream channel length (TUCL) to prune the vector hydrography layer to a level of detail that matches the raster DEM grid resolution. Network pruning reduces the occurrence of erroneous stream piracy caused by the rasterization of multiple stream links to the same DEM grid cell. The algorithm also restricts flow within individual stream reaches, further reducing erroneous stream piracy. In situations where two vector stream features occupy the same grid cell, the new tool ensures that the larger stream, designated by higher TUCL, is given priority. TUCL‐based priority minimizes the impact of the topological errors that occur during the stream rasterization process on modeled regional drainage patterns. The test data demonstrated that TopologicalBreachBurn produces highly accurate and scale‐insensitive drainage patterns and watershed boundaries. The drainage divides of four large watersheds within the study region that were delineated from the TopologicalBreachBurn‐processed DEMs were found to be highly accurate when compared with the official watershed boundaries, even at the coarsest grid resolutions, with Kappa index of agreement values ranging from 0.952 to 0.921. The corresponding Kappa coefficient values for a traditional stream burning method (FillBurn) ranged from 0.953 to 0.490, demonstrating a significant decrease in mapping accuracy at coarser DEM grid resolutions. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

12.
Distributed hydrological models require a detailed definition of a watershed's internal drainage structure. The conventional approach to obtain this drainage structure is to use an eight flow direction matrix (D8) which is derived from a raster digital elevation model (DEM). However, this approach leads to a rather coarse drainage structure when monitoring or gauging stations need to be accurately located within a watershed. This is largely due to limitations of the D8 approach and the lack of information over flat areas and pits. The D8 approach alone is also unable to differentiate lakes from plain areas.

To avoid these problems a new approach, using a digital river and lake network (DRLN) as input in addition to the DEM, has been developed. This new approach allows for an accurate fit between the DRLN and the modelled drainage structure, which is represented by a flow direction matrix and a modelled watercourse network. More importantly, the identification of lakes within the modelled network is now possible. The proposed approach, which is largely rooted in the D8 approach, uses the DRLN to correct modelled flow directions and network calculations. For DEM cells overlapped by the DRLN, flow directions are determined using DRLN connections only. The flow directions of the other DEM cells are evaluated with the D8 approach which uses a DEM that has been modified as a function of distance to the DRLN.

The proposed approach has been tested on the Chaudière River watershed in southern Québec, Canada. The modelled watershed drainage structure showed a high level of coherence with the DRLN. A comparison between the results obtained with the D8 approach and those obtained by the proposed approach clearly demonstrated an improvement over the conventionally modelled drainage structure. The proposed approach will benefit hydrological models which require data such as a flow direction matrix, a river and lake network and sub-watersheds for drainage structure information.  相似文献   


13.
Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.  相似文献   

14.
Watershed delineation is a required step when conducting any spatially distributed hydrological modelling. Automated approaches are often proposed to delineate a watershed based on a river network extracted from the digital elevation model (DEM) using the deterministic eight‐neighbour (D8) method. However, a realistic river network cannot be derived from conventional DEM processing methods for a large flat area with a complex network of rivers, lakes, reservoirs, and polders, referred to as a plain river network region (PRNR). In this study, a new approach, which uses both hydrographic features and DEM, has been developed to address the problems of watershed delineation in PRNR. It extracts the river nodes and determines the flow directions of the river network based on a vector‐based hydrographic feature data model. The river network, lakes, reservoirs, and polders are then used to modify the flow directions of grid cells determined by D8 approach. The watershed is eventually delineated into four types of catchments including lakes, reservoirs, polders, and overland catchments based on the flow direction matrix and the location of river nodes. Multiple flow directions of grid cells are represented using a multi‐direction encoding method, and multiple outflows of catchments are also reflected in the topology of catchments. The proposed approach is applied to the western Taihu watershed in China. Comparisons between the results obtained from the D8 approach, the ‘stream burning’ approach, and those from the proposed approach clearly demonstrate an improvement of the new approach over the conventional approaches. This approach will benefit the development of distributed hydrological models in PRNR for the consideration of different types and multiple inlets and outlets of catchments. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

15.
Drainage networks are the basis for segmentation of watersheds, an essential component in hydrological modelling, biogeochemical applications, and resource management plans. With the rapidly increasing availability of topographic information as digital elevation models (DEMs), there have been many studies on DEM‐based drainage network extraction algorithms. Most of traditional drainage network extraction methods require preprocessing of the DEM in order to remove “spurious” sink, which can cause unrealistic results due to removal of real sinks as well. The least cost path (LCP) algorithm can deal with flow routing over sinks without altering data. However, the existing LCP implementations can only simulate either single flow direction or multiple flow direction over terrain surfaces. Nevertheless, terrain surfaces in the real world are usually very complicated including both convergent and divergent flow patterns. The triangular form‐based multiple flow (TFM) algorithm, one of the traditional drainage network extraction methods, can estimate both single flow and multiple flow patterns. Thus, in this paper, it is proposed to combine the advantages of the LCP algorithm and the TFM algorithm in order to improve the accuracy of drainage network extraction from the DEM. The proposed algorithm is evaluated by implementing a data‐independent assessment method based on four mathematical surfaces and validated against “true” stream networks from aerial photograph, respectively. The results show that when compared with other commonly used algorithms, the new algorithm provides better flow estimation and is able to estimate both convergent and divergent flow patterns well regarding the mathematical surfaces and the real‐world DEM.  相似文献   

16.
湖北牛山湖小型鱼类的群落结构及多样性   总被引:1,自引:1,他引:1  
流域数据模型是流域特征的语义、行为和规则的表达,是进行集成流域模拟和管理的空间数据组织的重要内容.本文以长江三角洲太湖流域上游的西苕溪流域为例,将汇流单元分为自然流域、“大包围”、圩区三种形式,分别建立了各自内部的河湖网络关系,即山区由自然流域单元,内包含水库(或湖泊)、河流、水工点的树状河湖网络关系;平原区由人工汇流单元,内包含大包围、圩区、湖泊、湿地、河道、水工点的网状河湖网络关系,为建立适合我国的流域数据模型做了有益的探索.  相似文献   

17.
Drainage networks delineated from Digital Elevation Models (DEMs), are the basis for the modelling of geomorphological and hydrological processes, biogeochemical cycling, and water resources management. Besides providing effective models of water flows, automatically extracted drainage networks based on topography can diverge from reality to varying degrees. The variability of such disagreement within catchments has rarely been examined as a function of the heterogeneity of land cover, soil type, and slope in the catchment of interest. This research gap might not only substantially limit our knowledge of the uncertainty of hydrological prediction, but can also cause problems for users attempting to use the data at a local scale. Using 1:100000 scale land cover maps, Quaternary deposits maps, and 2 m resolution DEMs, it is found that the accuracy of delineated drainage networks tends to be lower in areas with denser vegetation, lower hydraulic conductivity, and higher erodibility. The findings of this study could serve as a guide for the more thoughtful usage of delineated drainage networks in environmental planning, and in the uncertainty analysis of hydrological and biochemical predictions. Therefore, this study makes a first attempt at filling the knowledge gap described above.  相似文献   

18.
The digital elevation model (DEM) has become an essential tool for an increasing array of mountain runoff analyses, particularly the derivation and mapping of stream channel networks. This study examines how well commonly applied DEM‐based channel derivation methods at different spatial resolutions can represent the channel network for a glaciated Rocky Mountain headwater catchment. The specific objectives are to (1) examine how differences in gridded DEM resolution affect spatially distributed values of local slope, specific contributing area, and topographic wetness index derived from both eight and infinite directional flow algorithms, (2) map the actual stream channel network to examine the influence of surface variables on channel initiation, and (3) assess accuracy of DEM‐derived networks compared with the field surveyed network. Results show that for the same contributing area threshold, increasing grid cell size leads to increased channelization of modeled networks. A plot of local slope versus contributing area reveals a negative relationship similar to that of prior studies in un‐glaciated areas but with breaks in slope at contributing areas that are too small to represent thresholds for channelization. Field survey results and evaluation of DEM‐derived channel networks suggest that channel network formation is not clearly related to surface topographic variables at Loch Vale. Digitally derived channel networks do not accurately predict low order channel locations, but approximations of the channel network with drainage density and headward extent of channelization similar to the observed network can be derived with both a 1 m and 10 m DEM using a contributing area threshold of approximately 4x104 m2. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

19.
Predicting runoff and erosion from watersheds burned by wildfires requires an understanding of the three-dimensional structure of both hillslope and channel drainage networks. We investigate the small- and large-scale structures of drainage networks using field studies and computer analysis of 30- m digital elevation model. Topologic variables were derived from a composite 30-m DEM, which included 14 order 6 watersheds within the Pikes Peak batholith. Both topologic and hydraulic variables were measured in the field in two smaller burned watersheds (3.7 and 7.0 hectares) located within one of the order 6 watersheds burned by the 1996 Buffalo Creek Fire in Central Colorado. Horton ratios of topologic variables (stream number, drainage area, stream length, and stream slope) for small-scale and large-scale watersheds are shown to scale geometrically with stream order (i.e., to be scale invariant). However, the ratios derived for the large-scale drainage networks could not be used to predict the rill and gully drainage network structure. Hydraulic variables (width, depth, cross- sectional area, and bed roughness) for small-scale drainage networks were found to be scale invariant across 3 to 4 stream orders. The relation between hydraulic radius and cross-sectional area is similar for fills and gullies, suggesting that their geometry can be treated similarly in hydraulic modeling. Additionally, the rills and gullies have relatively small width-to-depth ratios, implying sidewall friction may be important to the erosion and evolutionary process relative to main stem channels.  相似文献   

20.
Automatically extracting drainage networks from digital elevation models coupled with the constant stream threshold value is a regular method. These extracted networks can be verified by comparing the channel initiation points with those from real networks. From the results analysed, the differences in channel initiation points will affect the network geometries, geomorphological indices and hydrological responses. This paper develops two automatic algorithms, the headwater‐tracing method and the fitness index, to trace the flow paths from headwaters to the outlet and to calculate the reasonable stream threshold. Instead of the method determined by trial and error or field survey, the accurate channel initiation points can be obtained from airborne photographs coupled with high‐resolution SPOT images for suitable drainage network extraction. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

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

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