首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 322 毫秒
1.
Rockfall, up to several hundreds of cubic meters, is a frequent and rapid landslide which menaces extensive areas in mountainous territories. Rockfall susceptibility zoning map at a large scale (1:5000–1:25 000) can be the first tool for land use planning in order to manage rockfall risk. A methodology allowing to analyze susceptibility in extensive areas with optimum cost/benefit relationship is needed. This work analyzes rockfall susceptibility in an extensive rocky mountain of the Principality of Andorra (Pyrenees Mountains), first on the rock slope and then on the exposed area located below. The rockfall record, obtained by means of geomorphological analysis, supplies the main data to analyze the susceptibility on the rock slope. An additional historical inventory verifies the accuracy of rockfall sizes recorded by means of the geomorphological analysis. According to the classification recommended by the Guidelines of Joint Technical Committee, the density of rockfall features on the rock slope assesses susceptibility in four levels. Subsequently, susceptibility on exposed areas has been analyzed by means of reach probability of rock blocks analysis using empirical models. Data acquired from thirteen recent events, from 1999 to 2004, have been used to verify the accuracy of the two empirical models mainly used (reach angle and shadow angle). Five reach probability limits (1, 0.5, 0.25, 0.01, and 0) establish boundaries between susceptibility levels. The resulting rockfall susceptibility zoning map allows: (a) to identify land areas and human elements exposed to rockfalls and, (b) to establish several exposition levels. This map can be a useful and cost-effective tool for administrations responsible to manage natural risk in order to guide urban grow in extensive areas or decide upon work programs based on in-depth analysis (hazard and risk).  相似文献   

2.
A semi-quantitative heuristic methodology is developed to map a rockfall detachment susceptibility zonation of El Hierro Island (Canary Archipelago). The rationalized procedure, which we called non-weighted bounded indicators, is based on overlapping thematic maps of conditioning factors to mass movement, which are appropriately and individually rescaled and then composed by addition to obtain a susceptibility numerical index through a GIS. As the consistency of the geomorphological analysis depends on the expert subjective criteria and the appropriate interpretation of the landscape, the use of this methodology reduces subjectivity and quantifies the degree of susceptibility. The main factors affecting the mass movement phenomena (rockfalls events), also recognized in the field and, therefore, considered in the presented GIS arrangement, are slope, profile curvature, lithology, vegetation cover and dykes density. To calculate the slope threshold or minimum angle characteristic of rockfall source areas, mixed Gaussian slope frequency decomposition is used. The curvature index reveals stepwise areas. Qualitative geomechanical characteristics are linked to a quantitative index according to a volcanic lithological-complexes classification. Both destabilization (root-wedging) and stabilization effects are considered into the vegetation cover index. The dyke density index incorporates the bearing rock capacity decrease produced in the halo around a dyke network intrusion. Slope, curvature and vegetation indexes thresholds have been fitted following field observations. A rockfall detachment susceptibility map is obtained and classified based on the histogram maxima. The rockfall inventory, based on rockfall events reported within the island, was used for the model validation. A 12?% of the whole island shows medium to very high susceptibility.  相似文献   

3.
本文以北京市怀柔区为例,通过现场调查,对688处崩塌灾害分别以面数据和点数据的形式获取了两套编目图。根据现场调查和资料分析,选取岩性、地形、断裂和道路建设作为该区崩塌灾害的主控因素,采用频率比(FR)模型对崩塌灾害的易发性进行了评价。为了对评价结果的预测性进行检验,采用随机分割法,选取了415处崩塌用于频率比模型的计算,剩余的273处崩塌用于评价结果预测性的验证。预测曲线表明,基于崩塌面数据的评价结果比基于点数据的评价结果具有明显的优越性。根据基于面数据的频率比模型评价结果,可以将研究区的崩塌灾害易发性划分为5个等级:较低易发(占全区14%)、低易发(占全区20%)、中等易发(占全区27%)、高易发(占全区22%)和极高易发(占全区17%)。相关工作和结论可以为区域地质灾害易发性评价中编目图的编制提供参考,并为怀柔区区域国土利用和防灾减灾提供指导。  相似文献   

4.
Dramatic effects resulting from landslides on human life and economy of many nations are observed sometimes throughout the world. Landslide inventory and susceptibility mapping studies are accepted as the first stage of landslide hazard mitigation efforts. Generally, these landslide inventory studies include identification and location of landslides. The main benefit is to provide a basis for statistical susceptibility zoning studies. In the present study, a landslide susceptibility zoning near Yenice (NW Turkey) is carried out using the factor analysis approach. The study area is approximately 64 km2 and 57 landslides were identified in this area. The area is covered completely by Ulus Formation that has a flysh-like character. Slope angle, elevation, slope aspect, land-use, weathering depth and water conditions were considered as the main conditioning factors while the heavy precipitation is the main trigger for landsliding. According to the results of factor analysis, the importance weights for slope angle, land-use, elevation, dip direction, water conditions and weathering depth were determined as 45.2%, 22.4%, 12.5%, 8.8%, 8.1% and 3.0% respectively. Also, using these weights and the membership values of each conditioning factor, the membership value for landslide susceptibility was introduced. In the study area, the lowest membership value for landslide susceptibility was calculated as 0.20. Consequently, combining all results, a landslide susceptibility map was obtained. Compared with the obtained map, a great majority of the landslides (86 %) identified in the field were found to be located in susceptible and highly susceptible zones.  相似文献   

5.
The main objectives of this paper are to design and evaluate a hybrid approach based on Gaussian mixture model (GMM) and random forest (RF) for detecting rockfall source areas using airborne laser scanning data. The former model was used to calculate automatically slope angle thresholds for different type of landslides such as shallow, translational, rotational, rotational-translational, complex, debris flow, and rockfalls. After calculating the slope angle thresholds, a homogenous morphometric land use area (HMLA) was constructed to improve the performance of the model computations and reduce the sensitivity of the model to the variations in different conditioning factors. After that, the support vector machine (SVM) was applied in addition to backward elimination (BE) to select and rank the conditioning factors considering the type of landslides. Then, different machine learning methods [artificial neural network (ANN), logistic regression (LR), and random forest (RF) were trained with the selected best factors and previously prepared inventory datasets. The best fit method (RF) was then used to generate the probability maps and then the source areas were detected by combining the slope raster (reclassified according to the thresholds found by the GMM model) and the probability maps. The accuracy assessment shows that the proposed hybrid model could detect the potential rockfalls with an accuracy of 0.92 based on training data and 0.96 on validation data. Overall, the proposed model is an efficient model for identifying rockfall source areas in the presence of other types of landslides with an accepted generalization performance.  相似文献   

6.
地震型砾岩崩塌体运动学特征分析   总被引:1,自引:0,他引:1       下载免费PDF全文
汶川地震在工程区诱发了大量的崩塌灾害,提供了大量的"天然崩塌试验"现象。本文以汶川地震诱发的都江堰青城后山红岩度假区的大型砾岩崩塌体作为研究对象,主要分析了它的运动学特征。通过测量获得崩塌体最大位移块石运移边界角(reach of angle)。结合现场调查和数值模拟,得出坡体植被及人工构筑物的遮挡作用是导致研究对象运移边界角数值偏大的主要因素。统计该崩塌堆积体沿沟谷延伸方向的块石分布情况,并分析沟谷拐点处块石大小比例及各大小块石数占统计区相应大小块石总数的比例。得出大块石(>1m3)在滚落过程中,平均约有总数量的1/4被截留在了拐点处。作者认为,充分分析自然力作用下的典型崩塌体的分布及运移等特征对于崩塌的运动学研究及人类的防灾减灾工程具有重要的参考和指导。  相似文献   

7.
Rockfalls strongly influence the evolution of steep rocky landscapes and represent a significant hazard in mountainous areas. Defining the most probable future rockfall source areas is of primary importance for both geomorphological investigations and hazard assessment. Thus, a need exists to understand which areas of a steep cliff are more likely to be affected by a rockfall. An important analytical gap exists between regional rockfall susceptibility studies and block-specific geomechanical calculations. Here we present methods for quantifying rockfall susceptibility at the cliff scale, which is suitable for sub-regional hazard assessment (hundreds to thousands of square meters). Our methods use three-dimensional point clouds acquired by terrestrial laser scanning to quantify the fracture patterns and compute failure mechanisms for planar, wedge, and toppling failures on vertical and overhanging rock walls. As a part of this work, we developed a rockfall susceptibility index for each type of failure mechanism according to the interaction between the discontinuities and the local cliff orientation. The susceptibility for slope parallel exfoliation-type failures, which are generally hard to identify, is partly captured by planar and toppling susceptibility indexes. We tested the methods for detecting the most susceptible rockfall source areas on two famously steep landscapes, Yosemite Valley (California, USA) and the Drus in the Mont-Blanc massif (France). Our rockfall susceptibility models show good correspondence with active rockfall sources. The methods offer new tools for investigating rockfall hazard and improving our understanding of rockfall processes.  相似文献   

8.
This study shows a rockfall susceptibility assessment at local scale in north Beijing of China, including the identification of rockfall sources onsite by terrain and rock discontinuities analysis and run-out distance prediction by Rocfall? simulation. Two types of rockfall were defined including one type on the cliffs with long inclined slopes and another type on the road slopes with low height. Two historical rockfall events were used to back-calibrating the parameters used for run-out distance simulation. Based on the work, rockfall susceptibility map at local scale was created in GIS, which was compared with the map obtained at regional scale (entire Huairou district scale). Due to the difference of approaches applied, procedure of assessment and types of source data acquired, the two resulting rockfall susceptibility maps are proved to be different. Still, both of them are useful and could be used at different level’s decision for rockfall prevention and mitigation. Different types of uncertainties exist in the study of rockfall susceptibility assessment. To reduce the uncertainties, studies on both approaches and techniques are suggested.  相似文献   

9.
Rockfall hazard and risk analysis for Monte da Lua,Sintra, Portugal   总被引:1,自引:0,他引:1  
The prediction of rockfall trajectories below a rock cliff is essential in susceptibility, hazard and risk maps, particularity close to populated areas. The Monte da Lua hill area in Portugal, a tourist destination close to the historic city of Sintra (UNESCO World Heritage), is a typical granite boulder chaos landscape where from time to time rockfalls occur, the last such event having occurred on 29 January 2002. This area is therefore suitable to develop a rockfall study in order to provide hazard and risk maps a basis for mitigation measures. A preliminary investigation of the area leads to the identification of 188 potentially dangerous boulders. Detailed locations and geotechnical characteristics in terms of geometry, strength and context were sampled for each boulder. Digital elevations at 1 × 1 m resolution, known rockfall trajectory and building locations are provided in a GIS project for the study together with the spatial database of boulder characteristics. The modelling approach was conducted in two steps: (1) discrimination of the boulders in terms of static and dynamic mobility behaviour with multivariate analysis; (2) stochastic simulation of rockfall trajectories. The rockfall trajectory algorithm proposed is straightforward and is only dependent on elevation data, initial location of boulders and a friction angle. Due to the slope of the area, it assumes that rockfall is always of the rolling or sliding type. The friction angle was calibrated on the basis of the rockfall travel distance recorded on 29 January 2002 and generates simulated “realistic” trajectories. A smaller friction angle increases all simulated trajectories, leading to more “pessimistic” scenarios. The combined analysis of trajectories and potential damage to buildings and discrimination in terms of static and dynamic behaviour provides a final table in which all 188 sampled boulders are classified in one of the five risk grades.  相似文献   

10.
山区输变电工程崩塌随着国家“一带一路”和“全球能源互联网”战略的构建和实施,山区复杂地质环境条件下的输变电工程建设越来越多。然而,现阶段对于输变电工程崩塌(滚石)灾害识别、预测和评价方法等方面的研究相对较少,难以有效解决山区崩塌(滚石)灾害频发与输变电工程安全建设、运营之间不断增长的矛盾。鉴于此,作者借助航空遥感、无人机航测、岩体结构分析工具、崩塌(滚石)运动数值模拟软件等新技术方法等,对山区输变电工程中崩塌(滚石)灾害识别与预测方法等开展了相关研究。提出:(1)基于区域工程地质分区-区段遥感分析-山体无人机航拍-岩体结构特征分析的思路,实现对长距离、大区域山区输变电工程的崩塌(滚石)危岩体快速识别和分析;(2)利用三维崩塌(运动)模拟方法快速识别崩塌(滚石)在杆塔(变电站)范围的到达情况,然后再利用多种方法对比确定能到达杆塔(变电站)范围内崩塌(滚石)的冲击特征;(3)山区输变电工程崩塌(滚石)灾害风险计算方法。  相似文献   

11.
In this paper a tool for semi-quantitative susceptibility assessment at a regional scale is presented which is applicable at areas with complex geological setting. At a study area within the Northern Calcareous Alps geotechnical mappings were implemented into a Geographical Information System and analysed as grid data with a cell size of 25 m. The susceptibility to sliding and falling processes was considered according to five classes (very low, low, medium, high, very high). Susceptibility to sliding was analysed using an index method. The layers of lithology, bedding conditions, tectonic faults, slope angle, slope aspect, vegetation and erosion were combined iteratively. Dropout zones of rockfall material were determined with help of a Digital Elevation Model. The movement of rolling rock samples was modelled by a cost analysis of all potential rockfall trajectories. These trajectories were also divided into five susceptibility classes. The susceptibility maps are presented in a general way to be used by communities and spatial planners. Conflict areas of susceptibility and landuse were located and can be presented destinctively.  相似文献   

12.
山西壶关太行山大峡谷景区为中国最美十大峡谷之一,但景区落石灾害频发,严重威胁景区安全运营。本文基于高精度地形信息与岩土体强度特性,采用坡度角分布方法开展区域尺度潜在落石源区识别,并引入岩体破坏敏感性指标定量描述潜在落石源区失稳概率。然后,利用经验模型Flow-R模拟落石运动扩散过程,获取落石的传播概率与能量分布情况。最后,提出落石危险性双因子评价模型实现落石危险性定量评估。获得主要结论如下:(1)研究区内潜在落石源区面积为25.7 km2(35.7%),主要以条带状分布于峡谷两侧陡壁。其中岩体破坏高敏感性区为3.3 km2。(2)研究区落石高危险区面积达3.22 km2,主要威胁景区内游客集散地与交通线路,尤其在S327荫林线红豆峡入口处落石危险性最高。(3)野外调查验证结果表明了应用坡度角分布方法识别潜在落石源区的高效性与准确性,提出的双因子评价模型可为峡谷区落石危险性评估提供快速解决方案。本文提出的“区域落石源区识别-源区失稳概率分析-落石危险性评估”的一整套技术方案能够为类似的高山峡谷区落石灾害早期识别及风险防控提供技术参考。  相似文献   

13.
In this study we show the application of a long-range Terrestrial Laser Scanner (TLS) to a detailed rockfall study in a test zone at Vall de Núria, located in the Eastern Pyrenees. Data acquisition was carried out using TLS-Ilris3D, the new generation of reflector-less laser scanners with a high range, accuracy and velocity of measurements. Eight scans were performed at 3 stations to acquire coordinates of almost 4 million points. The results from the acquired data are a high accuracy Digital Elevation Model (DEM) and the reconstruction of the joint geometry. The former is used for inventory of rockfalls and for more accurate rockfall simulation (trajectories and velocities). The latter allows us to model the geometry and volume of the source area in recent rockfalls. Our findings suggest that TLS technology could be a tool of reference in rockfall studies in the near future.  相似文献   

14.
Predictive GIS-Based Model of Rockfall Activity in Mountain Cliffs   总被引:6,自引:1,他引:6  
Rockfall susceptibility has been analysed in mountain cliffs of the Cantabrian Range, North Spain. The main aim of this analysis has been to build a predictive model of rockfall activity from a low number of environmental and geological variables. The rockfall activity has been quantified in a GIS. The cartographic information used shows the spatial distribution of all the recent talus screes as well as their associated source areas in the rock-slopes. The area relation At/Ar (recent talus scree polygon/source basins) in the rock slopes has been used as the rockfall activity indicator. This relation has been validated in 50 pilot rock-slopes and compared with the relation number of recent rock fragments/source basin, obtained from field work. The environmental factors causing rockfall depend on the rock slope situation, and these are: altitude and sun radiation on the rock cliff. The geological factors considered are: lithology, relative position of the main discontinuities with respect to the topographic surface and two morphologic parameters: the roughness and slope gradient. A logistic regression analysis has been applied to a population of 442 limestone and quartzite rock cliffs. The dependent variable is the rockfall activity indicator, which allows the definition of two classes of rock cliff units: low and high activity. The independent variables are altitude, sun radiation (equinox radiation, summer solstice radiation, winter solstice radiation), slope roughness, slope gradient,anisotropy and lithology. Results suggest that it is possible tobuild a valid cartographic predictive model for rockfall activity in mountain rock cliffs from a limited number of easily obtainable variables. The method is especially applicable in massive rock slopes or in regions with uniform rock mass characteristics.  相似文献   

15.
滚石灾害是山区常见的地质灾害类型,研究滚石的运动特征对地质灾害调查及危险性评估有着重要意义。通过对尼泊尔某项目滚石灾害后现场进行工程地质调查,分析其灾害成因机制,查明事故原因。调查结果表明:9 ·15灾害非人类活动的影响,属自然地质灾害,造成事故的主要原因为超高位岩体崩塌,而滚石的范围又超过前期预测的危险区。通过现场痕迹分析,推测出滚石的运动路径。根据调查出的撞击点位置、物质组成及几何特征,作者提出运用运动学原理还原滚石运动轨迹,并利用rocfall软件对超高位危岩体崩落后的运动轨迹进行随机模拟分析,推算出超高位危岩体崩落后能量大小的变化,为防护措施方案提供可靠的依据。张口式帘式网韧性强,防护能级高,对高陡边坡滚石灾害能起到很好的效果。本文可为类似高陡边坡的危岩治理防护提供参考。  相似文献   

16.
竹山下铀矿床是粤北下庄铀矿田内大型铀矿床之一,铀矿化类型为"交点"型和硅化带型。在详细的野外地质调查基础上,对竹山下铀矿床4种不同类型黄铁矿进行元素含量分析及硫同位素测试,结果表明:"交点"型矿石中黄铁矿相对富集Pb、 Cu、 Co、 As、 Ni、 Se、 Bi、 U、 Sb、 Zn等微量元素;"交点"型铀矿化形成于中深部高温环境,成矿热液具有地幔流体特征,成矿过程硫来源与该区花岗岩中黄铁矿的硫来源一致或者相似,花岗岩中的黄铁矿可能为该期成矿事件的产物;竹山下矿床在垂向上表现出越往深部硫逸度越低的特征;"交点"型和硅化带型中黄铁矿具有相似的微量元素配分曲线,表明二者具有相同的成矿热液来源,且与辉绿岩中黄铁矿配分曲线相似,表明该区成矿热液具有深源性。  相似文献   

17.
浙西梅雨滑坡易发性评价模型对比   总被引:1,自引:0,他引:1       下载免费PDF全文
我国目前滑坡易发性评价研究主要集中在西南地区,对东南部降雨引发特别是梅雨引发的滑坡研究较少.选取浙江省西北部梅雨控制区淳安县为研究区,通过遥感解译结合野外详细调查,共确定滑坡596处,并建立滑坡编录数据库.选取高程、坡向、坡度、曲率、工程岩组、断层、道路、建设用地、植被等9个滑坡影响因子,基于GIS栅格分析方法,采用人工神经网络(ANN)、logistic回归和信息量3种评价模型,分别对32种不同影响因子组合进行滑坡易发性对比评价,得到滑坡易发性指数图.应用评价曲线下面积AUC(area under curve)对评价结果进行检验,ANN、logistic回归和信息量3种模型的正确率分别是93.75%、89.76%和90.06%;采用淳安县2014年梅汛期发生的13处滑坡作为预测样本,3种模型预测率分别是94.75%、94.33%和77.21%.上述分析结果表明:ANN模型优于其他两者.以ANN模型评价结果指数图为基础进行易发性分区,采用滑坡强度指标进行分区结果检验,滑坡强度值由易发性低、较低、中和高依次递增,说明分区结果合理.研究成果可以为浙西降雨型滑坡特别是由梅雨引发滑坡的易发性评价提供参考.   相似文献   

18.
The Tramuntana range, in the northwestern sector of the island of Mallorca (Balearic Islands, Spain), is frequently affected by rockfalls which have caused significant damage, mainly along the road network. In this work, we present the procedure we have applied to calibrate and validate rockfall modelling in this region, using 103 cases of the available detailed rockfall inventory (630 rockfall events collected since the eighteenth century). We have exploited STONE (Guzzetti et al. 2002), a GIS-based rockfall simulation software which computes 2D and 3D rockfall trajectories starting from a DTM and maps of the dynamic rolling friction coefficient and of the normal and tangential energy restitution coefficients. The appropriate identification of these parameters determines the accuracy of the simulation. To calibrate them, we have selected 40 rockfalls along the range which include a wide variety of outcropping lithologies. Coefficients values have been changed in numerous attempts in order to select those where the extent and shape of the simulation matched the field mapping. Best results were summarized with the average statistical values for each parameter and for each geotechnical unit, determining that mode values represent more precisely the data. Initially, for the validation stage, 10 well-known rockfalls exploited in the calibration phase have been selected. Confidence tests have been applied to their modelling results taking into account not only the success but also the mistakes. The best accuracy is obtained when the rockfall has a preferential trajectory and worse results when the rockfall follows two or more trajectories. Additionally, the greater the rockfall runout length, the less precise the simulation is. We have further validated the calibrated parameters along the Ma-road (111 km), the main transportation corridor in the range, using 63 rockfall events that occurred during the past 18 years along the road. Of the rockfalls where source areas were properly identified, 81.5 % are well represented by STONE modelling, as the travel paths and the depositional areas are successfully ascertained. Results of the analysis have been used by the Road Maintenance Service of Mallorca to assess hazard and risk posed by rockfall at regional scale to design the road management plan.  相似文献   

19.
本文系《七五》国家重点攻关项目中三级专题湘南锡铅锌多金属隐伏矿预测研究的部分成果。在系统研究多金属成矿地质环境后,从解剖重点矿田、矿床入手,研究区内成矿规律,分析两类(同熔、重熔)岩浆的同熔、重熔、分异演化机制和沉积成岩机理与成矿的内在联系,探讨物质来源,成矿溶液性质和矿床成因,进而划分出同熔岩浆、重熔岩浆和层控等三大成矿系列,并依据成矿方式、围岩蚀变、矿物组合、矿体产状等细分出22个矿化类型(亚类),建立了成矿模式。依据成矿与地质环境的关系及矿床分带特征,总结了寻找隐伏矿床的找矿模式(标志)和找矿系统,并进行了成矿预测。  相似文献   

20.
Multi-scenario Rockfall Hazard Assessment Using LiDAR Data and GIS   总被引:1,自引:0,他引:1  
Transportation corridors that pass through mountainous or hilly areas are prone to rockfall hazard. Rockfall incidents in such areas can cause human fatalities and damage to properties in addition to transportation interruptions. In Malaysia, the North–South Expressway is the most significant expressway that operates as the backbone of the peninsula. A portion of this expressway in Jelapang was chosen as the site of rockfall hazard assessment in multiple scenarios. Light detection and ranging techniques are indispensable in capturing high-resolution digital elevation models related to geohazard studies. An airborne laser scanner was used to create a high-density point cloud of the study area. The use of 3D rockfall process modeling in combination with geographic information system (GIS) is a beneficial tool in rockfall hazard studies. In this study, a 3D rockfall model integrated into GIS was used to derive rockfall trajectories and velocity associated with them in multiple scenarios based on a range of mechanical parameter values (coefficients of restitution and friction angle). Rockfall characteristics in terms of frequency, height, and energy were determined through raster modeling. Analytic hierarchy process (AHP) was used to compute the weight of each rockfall characteristic raster that affects rockfall hazard. A spatial model that considers rockfall characteristics was conducted to produce a rockfall hazard map. Moreover, a barrier location was proposed to eliminate rockfall hazard. As a result, rockfall trajectories and their characteristics were derived. The result of AHP shows that rockfall hazard was significantly influenced by rockfall energy and then by frequency and height. The areas at risk were delineated and the hazard percentage along the expressway was observed and demonstrated. The result also shows that with increasing mechanical parameter values, the rockfall trajectories and their characteristics, and consequently rockfall hazard, were increased. In addition, the suggested barrier effectively restrained most of the rockfall trajectories and eliminated the hazard along the expressway. This study can serve not only as a guide for a comprehensive investigation of rockfall hazard but also as a reference that decision makers can use in designing a risk mitigation method. Furthermore, this study is applicable in any rockfall study, especially in situations where mechanical parameters have no specific values.  相似文献   

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

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