首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 71 毫秒
1.
基于遥感的区域蒸散量监测方法——ETWatch   总被引:11,自引:1,他引:10       下载免费PDF全文
ETWatch是用于区域蒸散遥感监测的业务化运行系统,集成了具有不同应用优势的遥感蒸散模型,并以Pen-man-Monteith公式为基础建立下垫面表面阻抗估算方法,利用逐日气象数据与遥感反演参数,获得逐日连续的蒸散分布图。所生成的从流域级到地块级的数据产品能动态反映区域蒸散发的时空变化规律。基于ETWatch方法,对2002-2005年海河流域和河北省馆陶县的蒸散状况进行了连续监测。利用地面观测资料对蒸散遥感监测产品的验证表明,涡度相关观测数据能量闭合率在0.9以上时,1 km级的日蒸散结果的平均偏差约10%;对于地块级(30 m)的ET,受混合像元的影响,相对于土壤水分消耗法测量的作物蒸散发,遥感监测作物蒸散发的平均误差在12.7%左右。  相似文献   

2.
基于多尺度遥感数据估算地表通量的方法及其验证分析   总被引:2,自引:0,他引:2  
地表水热通量(显热通量、潜热通量)的遥感估算在全球气候变化、水资源、生态环境等研究领域具有重要的应用价值.MODIS数据的空间分辨率较低(热红外波段星下点为1 km),而地球表面的几何物理属性又具有高度非均匀性,因而在实际应用中面临较严重的尺度问题.探讨了多源卫星数据(中高分辨率Landsat TM与中低分辨率MODIS)相结合佑算像元通量的2种方法,分别利用高分辫率的地表分类及植被指数信息在混合像元内部进行亚像元处理,以提高非均匀地表混合像元的通量估算精度.研究数据来自于2008年黑河流域综合实验获取的遥感数据和辅助数据,验证数据来自于实验期间获取的不同下垫面的地表通量数据,包括涡度相关(EC)数据,以及大孔径闪烁仪(LAS)数据.计算结果表明,2种方法皆可在下垫面不均匀或者地表类型较复杂的情况下得到比较明显的纠正效果,纠正后的通量与观测更加接近.相比之下,利用植被指数分解温度的方法适用性更广,纠正效果更好.在地面验证中,对比分析了EC和LAS数据在TM尺度和MO-DIS尺度通量验证的适用性.LAS数据测量尺度与MODIS卫星像元尺度相匹配,可以直接验证MODIS通量计算结果,EC数据虽然可以直接验证TM计算的通量,但与MODIS数据对比,还需要进行尺度转换,即先用EC验证TM通量,然后将TM通量降尺度,与MODIS进行对比.最后对利用LAS验证通量的不确定性进行了分析,发现图像中LAS测点的几何定位误差以及LAS测量路径中像元的选取都对验证结果有一定影响.  相似文献   

3.
海底滑坡、浊流等深海底地质灾害严重威胁海洋工程安全,是国家深海开发亟待解决的风险问题。为避免深海海底地质灾害对海底工程造成危害,解决深海海底地质灾害监测预警的难题,我们研发了一套复杂深海工程地质原位长期监测系统。该系统通过声学、电阻率、超孔隙水压力等方法监测深海海底沉积物的物理力学性质变化,实现了对深海海底地质灾害的监测和预警。该系统主要包括海床基搭载平台、监测系统、通信控制系统、供电系统等。其中监测系统主要通过原位长期监测海底沉积物的电阻率、声学、超孔隙水压力等的变化来获取海底沉积物的物理力学性质变化;通信控制系统可以实现海底到海面,再到陆地的双向通信和数据传输。其中供电系统通过独特设计的海水电池工艺,可以满足该系统在海底长期工作一年的电量需求。复杂深海工程地质原位长期监测系统已完成了近海测试,并搭载“海洋地质六号”“东方红三号” “张謇号”等科考船在南海进行了多次远海海试,获取了丰富的实测数据。电阻率监测系统采用温纳法滚动测量,测得的水土界面位置平均电阻率为0.207 Ω·m。超孔隙水压力监测系统采用开放式结构的压差式光纤光栅孔压测量方法,监测到孔压观测的4个标志性阶段:(1)贯入过程引起的超孔隙水压力累计,峰值为34.942 kPa,历时0.182 h;(2)贯入完成后累积的超孔隙水压力衰减,衰减到9.973 kPa,历时为0.810 h;(3)环境应力引起的超孔隙水压力实时响应,超孔隙水压力的变化范围为8.327~14.384 kPa;(4)残余孔隙水压力平均值为11.150 kPa。声学监测系统采用两个一发三收模式,测量的海水平均声速为1 533 m/s,测量的海底沉积物自上而下的平均声速依次为1 586、1 587、1 784、1 735、1 831 m/s。复杂深海工程地质原位长期监测系统的成功研制将显著提升目前海洋工程地质原位长期观测的技术能力,解决复杂深海工程地质评价及地质灾害监测预警的技术难题。  相似文献   

4.
Heterogeneity effects often limit the accuracy of synchrotron X-ray fluorescence microprobe elemental analysis data to ± 30%. The difference in matrix mass absorption at Kα and Kβ fluorescence energies of a particular element can be exploited to yield information on the average depth-position of the element or account for heterogeneity effects. Using this technique, the heterogeneous distribution of Cu in a simple layered sample could be resolved to a 2 × 2 × 10 (x, y, z, where z is the depth coordinate) micrometer scale; a depth-resolution limit was determined for the first transition metal series and several other elements in calcite and iron oxide matrices. For complex heterogeneous systems, determination of average element depth may be computationally limited but the influence of heterogeneity on fluorescence data may still be assessed. We used this method to compare solid-state diffusion with sample heterogeneity across the Ni-serpentine/calcite boundary of a rock from Panoche Creek, California. We previously reported that Ni fluorescence data may indicate solid state diffusion; in fact, sample heterogeneity in the depth dimension can also explain the Ni fluorescence data. Depth heterogeneity in samples can lead to misinterpretation of synchrotron X-ray microprobe results unless care is taken to account for the influence of heterogeneity on fluorescence data.  相似文献   

5.
基于不同空间尺度、长时间序列的地表反照率产品探究黑河流域2000—2012年的反照率时空变化特征。首先基于角度格网化(AB)算法对黑河流域30 m环境卫星一号HJ-1/CCD大气层顶方向反射率进行了地表反照率估算,作为高空间分辨率的地表反照率产品;选择同种算法计算的1 km空间分辨率的Global LAnd Surface Satellite(GLASS)反照率产品作为低空间分辨率反照率产品。结果表明:1利用AB算法反演的HJ卫星反照率具有较高精度,满足流域尺度反照率时空特征分析的精度要求;2黑河流域地表反照率空间分布差异显著,流域上游植被覆盖区域反照率较低,中下游荒漠地表反照率较高;3流域地表反照率的年内变化与季节性降雪和作物物候周期性变化一致。从季节变化角度,黑河流域反照率月平均值的年变化呈"U"字型,其中,冬季反照率最高,春季和秋季次之,夏季最低。从年际变化角度,黑河下游反照率呈上升趋势,上游和整个流域呈下降趋势。  相似文献   

6.
黑河流域观测通量的空间代表性研究   总被引:7,自引:0,他引:7  
地表湍流通量包括显热通量和潜热通量,它的准确测定对天气气候预报、农业节水和水资源管理等意义重大。目前应用较为广泛的通量测量技术有涡动相关仪(EC)和大孔径闪烁仪(LAS),前者是单点观测,观测范围通常只有几百米;后者近年来得到长足的发展,可以观测从几百米到十几公里尺度上的区域湍流通量。建立了EC和LAS通量解析足迹模型,详细介绍了足迹模型的理论依据和建立方法,并对二者模型的异同进行了阐释。在黑河流域遥感—地面观测同步试验中,以高寒与干旱区伴生的黑河流域为试验区,包括寒区水文试验、森林水文试验和干旱区水文试验,使用足迹模型分别分析了临泽草地站、阿柔冻融观测站的EC、LAS观测通量的源区,结合风向变化规律,对二者通量观测差异进行分析解释;分析了盈科灌区绿洲站、大野口关滩森林站的EC观测数据,以及阿柔冻融观测站的EC、LAS观测数据,选择具有代表性的月份,进行观测站点的空间代表性分析。结果表明:对通量观测站点进行空间代表性分析是十分必要的,可得到观测站点通量源区的时空变化特征,同时足迹模型在通量观测数据的分析中有很大的实用价值,可为今后通量观测数据的应用提供参考。  相似文献   

7.
土地利用/土地覆盖变化研究是近年来全球变化研究的焦点之一。全球和区域尺度的土地覆盖特征对全球环境状况的评估、模拟未来全球环境的情景有重要的作用。2000年在Internat ionalJournalofRemoteSensing杂志上出版了题为"GlobalandRegionalLandCoverCharacterizat ion from Remotely Sensed Data"的专辑。在此基础上,介绍、总结了国际上利用遥感影像进行全球和区域等大尺度土地覆盖研究的新进展。分别从数据源与制图的时空尺度、制图方法(数据预处理、分类、精度评估)等方面进行了介绍,并对现今的两个全球土地覆盖数据库进行了比较分析。  相似文献   

8.
区域蒸散发遥感估算方法及验证综述   总被引:7,自引:0,他引:7  
张荣华  杜君平  孙睿 《地球科学进展》2012,27(12):1295-1307
蒸散发是地表水热平衡的重要参量,也是农作物生长状况和产量的重要指标。与传统的蒸散发计算方法相比,遥感技术经济、适用、有效,在非均匀下垫面的区域蒸散发监测上具有明显的优越性。系统回顾了5种常用的区域蒸散发遥感估算方法,包括经验统计模型、与传统方法相结合的遥感模型、地表能量平衡模型、温度—植被指数特征空间模型以及陆面过程与数据同化等,分析了这些模型的最新研究进展及各自的优缺点,并对地表蒸散发的地面验证方法进行了概述。最后简要分析了区域蒸散发遥感估算存在的问题,并展望了其未来发展趋势。多源遥感数据协同反演与非遥感参数遥感化、蒸散发模型改进与多模型集成、陆面过程与遥感数据同化、遥感蒸散发估算及地面验证中的尺度问题与空间代表性问题研究等将会是未来区域蒸散发研究中的重点发展方向。  相似文献   

9.
黑河流域陆地—大气相互作用研究的几点思考   总被引:4,自引:2,他引:2  
陆地和大气之间是相互作用、相互影响的,不仅陆面分布与地表过程对大气变化有着响应过程,而且陆面物理过程、地表特征分布对大气过程也有着重要影响。我国西北地区地形分布极不均匀,是水资源变化和气候变化的敏感区与脆弱地区。这里,由于海拔高度起伏造就了以水为主线的上游山区冰雪、冻土—中游森林、绿洲—下游戈壁荒漠的多个自然景观带共存的内陆河流域,是开展陆—气相互作用研究的理想场所。文章集中黑河流域水循环过程、冰雪/冻土与大气相互作用、降水异质性、人类活动影响以及尺度等问题,阐述了流域研究现状以及研究热点,指出发挥各学科优势,将大尺度的大气过程与中间尺度和小尺度的陆面水文/生态过程结合,建立流域尺度大气—水文—生态过程真正的双向耦合模式,不仅研究气候变化对陆面水文、生态的影响,同时探讨陆面小尺度过程对大气的反馈机理,探索流域尺度的水循环过程和驱动机理,为水资源合理利用提供理论基础。  相似文献   

10.
准噶尔盆地东部中二叠统平地泉组(相当于芦草沟组)发育一套陆内裂谷背景下的湖泊相砂岩、泥岩、灰岩、白云岩、凝灰岩互层以及它们的过渡岩石,是区内最主要的烃源岩和储集层。近年来,在帐北断褶带、石树沟凹陷和吉木萨尔凹陷中二叠统平地泉组暗色泥岩、泥晶白云岩中发现了一类特殊的沉积岩,该类岩石发育类似火山岩的斑状结构。“斑晶”主要为粗晶白云石、方解石、黄铁矿及方沸石等,在岩心标本上呈“树枝状”或“雪花状”散布于基质中;基质则主要由泥晶白云石或泥质沉积物构成,富含有机质,发育水平层理及小型变形层理等。岩石学和矿物学分析表明,“斑晶”往往为粗晶方解石或白云石的集合体,也常见黄铁矿、方沸石与碳酸盐矿物共生,“斑晶”方解石发育环带,而白云石未见明显环带。电子探针分析表明“斑晶”方解石具有低镁、低铁及锶分布不均匀的特点,而“斑晶”白云石(FeO含量介于7.272%~11.086%之间)与“基质”泥晶白云石(FeO含量为1.027%)相比具有明显富铁的特点。流体包裹体分析表明“斑晶”方解石均一温度平均为180.68,℃,“斑晶”白云石均一温度平均为320.95,℃。这种特殊 “斑状”白云岩和“斑状”泥质岩很可能是湖底热液喷流作用的结果。当湖水沿深大断裂下渗至地下深处,与围岩发生物质交换并被加热后再沿断裂返回地表喷涌而出,热液流体携带的离子达到过饱和后就会析出方解石、白云石及黄铁矿等集合体,随热液的喷涌作用上升并破碎散落于湖底细粒沉积物内。“斑状”白云岩的发现对新疆北部中二叠统热水喷流沉积作用及该区油气的成因研究具有重要的意义。  相似文献   

11.
矿物-微生物交互作用广泛参与地球表层系统物质循环与能量流动过程,深刻地影响着一系列重要的地表生物地球化学进程。近年来地表半导体矿物的相关研究,为矿物-微生物交互作用提供了崭新研究方向,揭示地表“日光-半导体矿物-微生物”系统电子传递过程及其环境效应,是地质微生物学交叉领域研究的核心科学问题之一。本研究从地表不同生境“矿物膜”出发,以光电化学技术证实喀斯特、红壤、岩石漆“矿物膜”在1 000 min长时间循环实验中平均光电流值约为5.4、3.4、3.2 μA/cm2,证实“矿物膜”良好日光响应特性且铁锰氧化物矿物在其中发挥核心作用。基于笔者前期研究所发现的“矿物膜”电活性菌富集且与半导体矿物分布呈正相关性这一现象,本文进一步构建模拟光电子红壤细菌群落系统,20天后细菌群落α多样性显著提升,研究证实细菌群落具有模拟光电子响应活性,且电极与溶液群落均具有演化方向性;16S rRNA测序分析表明模拟光电子作用下ShewanellaPseudomonasStreptococcusLactobacillusAcinetobacter等电活性菌显著富集。综上,本文研究结果间接证实地表半导体矿物光电子可有效调控微生物群落结构并促进电活性菌在“矿物膜”中富集。  相似文献   

12.
国土整治生态修复中的农业景观生物多样性保护策略   总被引:1,自引:0,他引:1  
随着生态文明建设成为国家发展战略,国土整治修复需要更多生态科学和技术的支撑。农业景观生物多样性是重要的自然资源,为人类可持续发展提供了重要的生态系统服务,并提供了应对气候变化等风险的途径。本文分析了我国农业景观生物多样性保护及其对国土整治生态修复的重要性,国土整治生态修复中农业景观生物多样性保护面临的挑战,并总结了欧盟在农业景观生物多样性保护方面的经验和方法。最后,文章提出了将农业景观生物多样性保护整合到国土整治生态修复的四项策略,包括生态空间的“生态修复和管护-生态补偿-民生改善”策略、农业空间的“生态基础设施-生态系统服务提升-绿色发展”策略、城乡生活空间的“多功能景观-绿色宜居空间营造”策略、国土全域“生物多样性保护网络”策略,为生态文明背景下的国土整治生态修复提供参考。  相似文献   

13.
陆面过程研究的现状与发展趋势   总被引:25,自引:4,他引:25  
对目前陆面过程研究的现状及存在的问题进行了总结,并且提出了与边界层研究有关的几个前沿方向:①非均匀下垫面地表湍流通量的参数化及中尺度通量的参数化方法;②边界层理论在复杂地形条件下的应用,以及次网格地形对通量输送的影响;③进一步改进GCM中冠层地表水文参数化;④边界层的层云和层积云的参数化。  相似文献   

14.
As we know, desert pavement is a superficial feature from wind erosion widespread throughout arid lands and plays a dynamic role in geomorphic, hydrologic, and ecologic processes. In these regions, the plant cover is low or inexistent. Desert pavement is usually covered with large, small, and angular stones. Spatial distribution of stones at the surface is reflecting rainfall distribution at the land surface. Our objective in this study was to determine the relationship between the surface physical characteristics of desert pavement (meter scale), water movement, subsoil, and plant cover. The density and type of cover in land surface are also investigated. For this purpose, two different areas (Hajaligholi playa, Damghan) with different geology were selected. In each region, profiles were dogged. In site one, six different profiles were dogged, and each profile samples in different horizons (0–10, 10–20, 20–30, 30–40, and 40–50 cm) were being taken. In site two, three profiles were dogged, and the same horizon samples were being taken. The type and density of vegetation were determined. Finally, two sites were being compared. The results showed that strong relationship is between the stone cover, soil genesis, and plant distribution. Physicochemical properties (particle sorting, stones rates, electrical conductivity, pH, and gypsum) highly depend on the upland geological characteristics of the region. Also, soil structure and salt concentration have strong relationship with stone component. The distribution and density of vegetation show strong linkage with these properties. Also, in sites with desert pavements and dense plant cover, soils are protected and, thus, reduce desertification (loss fertile soils).  相似文献   

15.
通过对喀斯特地区乔木林、灌木林和草丛不同植被类型覆盖下的土面、石面、石缝、石沟、石洞、石槽和石坑7类小生境土壤样品分析测定,探讨了小生境土壤矿物组成及有机碳的空间异质性特征。结果表明:(1)不同植被类型下土壤矿物组成存在明显差异,即从草丛→灌木林→乔木林,小生境土壤中的石英、高岭石基本保持不变,但草丛土壤中的白云石矿物全部消失,灌木林、乔木林土壤中逐渐出现绿泥石;(2)随着植被正向演替,土壤有机碳含量显著增加,小生境土壤有机碳含量在不同植被类型条件下存在明显的异质性,表现为草丛>乔木林>灌木林。与此同时,以小生境土壤有机碳含量及变异系数为评价因子,分别对乔木林、灌木林及草丛中小生境有机碳进行聚类分析,结果显示,小生境微地貌特征是土壤有机碳含量及空间变异性的重要影响因子。对于土壤有机碳保护,草丛及灌木林区应以提高植被覆盖率为主;乔木林区则主要是保持现有良好的植被状况,减少人为干扰。   相似文献   

16.
长江源区小冬克玛底冰川区积雪消融特征及对气候的响应   总被引:6,自引:3,他引:3  
2005年使用花杆法在小冬克玛底冰川上进行了两轮积雪融化观测,结果显示:6月中、下旬,冰川区积雪消融基本与气温同步变化,但到7月上旬积雪消融发生了变化:其融化时间提前了2 h,在12:00融化量就达到一定高度;融化量大,日平均融化量较6月中下旬大0.71 mm水当量;最大积雪融化量出现时间滞后于气温最高时间约2 h,于18:00达到最大值.雪坑雪层剖面真实记录了积雪积累与融化过程中的雪层变化.积雪融化除受温度影响之外,还与风速等其它因素有密切关系.温度越高,融化量越大;风速与积雪融化量的关系相对复杂.积雪融化量只与1.5 m高度风速有关.在2h尺度上,近地面风速大,积雪融化量也大;在日尺度上,风速与积雪融化量的关系较好,但在两个观测时段表现出不同的相关关系.研究区积雪量很少,受昼夜气温变化与下垫面状况的影响,积雪融水几乎不产生径流,对河川径流的影响程度极小.  相似文献   

17.
新时期,地质分析理论在国土空间规划领域的应用日益广泛。为进一步优化关中盆地国土空间格局,从抓住关键因素、科学定量评价、阈值标准有据、结果可信适用的思路出发,该研究提出了基于木桶理论、风险理论和边际理论的资源环境承载能力和国土空间开发适宜性评价(简称“双评价”)理论框架与技术方法,并应用于关中盆地的双评价实践工作中。基于木桶理论,识别关中盆地活动断裂、地质灾害、富硒土地、引汉济渭水资源等4项关键因素;基于边际理论,开展资源承载能力评价;基于风险理论,开展地质环境承载能力、国土空间适宜性和地质安全评价;在单因素评价的基础上开展综合评价;将综合评价结果与区域国土空间开发现状或规划结果做叠加分析,调整和优化区域发展“三区三线”。研究结果表明基于木桶-风理-边际理论的双评价方法可行,评价结果可信:关中盆地87%的土地属于适宜农业生产和城镇化建设;尤其是引汉济渭工程使关中盆地成为承接东部人口和GDP发展的最佳地区;国土空间开发保护中要珍惜集中连片分布的富硒土地,作为永久农田予以保护;禁止在活动断裂分布范围和地质灾害高危险区进行工程建设,以国土空间规划和用途管制为抓手,从源头防控地震与地质灾害风险。  相似文献   

18.
Land use/land cover change is a global phenomenon which reflects natural resources degradation and/or utilization. Remote sensing and GIS have been widely used to monitor such changes at watershed level. The present study evaluates the LU/LC change during 1989 - 2001 in a semi-arid watershed of central India. Geocoded satellite data of 1989 and 2001 on 1:50,000 scale, were visually interpreted to prepare thematic maps which were later digitized using ArcGIS softwares. The analysis shows that vast tracts of cultivated land have become uncultivated and at some places even converted to wasteland. However, the land under dense forest and open forest has decreased due to expansion of built-up land and other anthropogenic activities. Increase in area of uncultivated land, wasteland and decrease in cultivated land and open scrub is also supported by rainfall analysis, which shows a declining trend and a fall of 186.93 mm in average annual rainfall for 1986-2003 period. The change detection map prepared using land use/land cover of 1989 and 2001 as inputs shows that out of the total geographical area of the watershed, 25.78% of the watershed area has seen a change from one land use category to another, however rest 74.22% has remained unchanged.  相似文献   

19.
One of the major uncertainties in the remote sensing estimates of regional evapotranspiration (ET) over heterogeneous landscapes is that significant estimation differences between different data scales always exist. It is necessary to predetermine the optimal remote sensing scale under such conditions. In this study, the effect of pixel scale on ET estimation over a typical oasis in north-western China was investigated. By the area-averaging and wavelet multi-resolution analysis aggregation techniques, four ??input?? up-scaling and ??output?? up-scaling procedures were implemented, respectively, to determine the optimal data scale. Validation suggested that high-resolution Landsat-based ET estimates could be used as the ground reference data when in situ measurements were not available. The differences between Landsat-based and MODIS-based ET estimates (i.e. the so-called Regional Error Distribution of the latter) showed that 1-km resolution MODIS data overestimated ET over the landscapes where NDVI was less than or equal to ~0.10; the data underestimated ET over the landscapes where NDVI was greater than or equal to ~0.40; ET could only be relatively accurately predicted by the data on the surfaces where NDVI ranged from 0.10 to 0.40. However, MODIS-based 1-km ET estimates could effectively reveal the predominant spatial distribution trend of regional ET. The results confirmed that the optimal remote sensing scale for modeling ET was approximately 480?m over the study area; land surface heterogeneity could cause significant errors in ET estimation once data scales exceeded this threshold; in addition, MODIS data at 250-m resolution were adequate for correcting the 1-km resolution ET estimates over heterogeneous landscapes.  相似文献   

20.
基于2001—2015年MOD10A1/MYD10A1、MOD13Q1以及相关气象数据,采用积雪持续时间比率法,监测了天山山区的季节雪线高程,分析了其变化特征及影响因子。结果表明:①近15年天山山区雪线整体呈显著上升趋势,平均高程3 680 m左右,其中,北坡、伊犁河谷、南坡季节雪线的稳定性依次减弱,平均高程分别为3 620 m、3 390 m及3 820 m;空间上雪线高程呈现南高北低、东高西低的纬度地带性分布特点。②年际尺度上,气温是影响天山山区雪线高程的主控因素,呈显著正相关,南北坡与之相同,但伊犁河谷则降水是影响其变化的主控因素,呈显著负相关;季节尺度上,夏季气温、冬季降水是影响雪线高程的主控因素,降水与其呈负相关,但气温较高的地区,夏秋季降水会促进积雪融化,使雪线高程上升;月尺度上,7月气温、1月降水对其影响最明显,且存在一定的滞后反应。③天山山区雪线高程比零度层低800 m左右,两者呈较好正相关;雪线高程与NDVI(Normalized Difference Vegetation Index)呈负相关,植被覆盖较好区域,同年NDVI与雪线高程相关性较好,植被覆盖较差区域,前一年NDVI与其相关性较好。  相似文献   

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

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