首页 | 官方网站   微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
以北京南部城市发展新区典型代表———房山、大兴、通州3处作为研究区,运用马尔科夫转移矩阵和景观格局指数的方法,进行土地利用变化和景观格局分析。通过对2001—2006年、2006—2011年两个时段的对比分析,发现2001—2006年研究区存在草地、林地、耕地等非建设用地类型面积下降,建设用地面积上升的趋势,且2006—2011年该趋势更加明显。结果表明研究区自然景观向人文景观变革迅速,并面临着后备土地资源不足等问题,亟需调整土地利用方式并优化景观系统空间格局。  相似文献   

2.
以具有重要生态地理位置的三峡库区(重庆段)为研究区,结合2000—2015年Landsat TM数据,运用统计和网络分析法,从时空角度研究景观格局和生态系统服务价值的变化及其两者间的相互关系。结果表明,整体而言,研究区耕地和林地所占比例较大,草地、水域、建设用地和其他用地所占比例较小,建设用地面积增加最快;从时间上看,2000—2015年斑块密度、最大斑块面积指数呈先减少后增加的趋势,聚集度最大值和香农多样性指数最小值出现在2010年;2000—2015年生态服务价值总量先增加后减少,2010年达最大;从空间上看,景观破碎度逐渐从库区东北、东南向各区县建成区及主城区转移,东北与东南地区的景观完整性逐渐提高;生态系统服务价值低值区集中于各区县建成区,主城区生态服务价值最低,长江流域西侧的生态服务价值低于东侧。合理的景观格局对提高该区域的生态系统服务价值具有重要意义。  相似文献   

3.
陕北黄土高原典型人类活动影响下景观生态风险评价   总被引:1,自引:0,他引:1  
植被恢复促进了景观格局和生态过程的改变,景观生态风险评价有助于揭示人类干预下生态系统结构、功能等在景观尺度产生的潜在不利生态效应。以生态脆弱、人与生态环境的相互作用敏感的陕北黄土高原为典型研究对象,根据2000、2005、2010和2015年夜间遥感地图、植被净初级生产力和土地利用数据,建立景观生态风险指数,对陕北黄土高原的综合风险状况及其变化进行定量评价,并对植被变化、城镇化分别与景观生态风险变化间的相关关系展开分析。结果表明,景观格局变化对生态稳定性有潜在影响,景观生态风险以增长趋势为主,2005—2010年景观生态风险增长的范围最大,涨幅为186.0%,其中,中风险和极高风险等级面积增长超2.7倍;2010—2015年景观生态风险增加有所缓解,增加12.9%;2000—2005年景观生态风险降低26.0%。景观生态风险低等级面积逐渐缩小,其他等级面积逐渐变大。其中,2005年后景观生态风险等级呈片状扩展,主要位于陕北西北部风沙区。植被变化、城镇化分别与景观生态风险值变化之间在西北风沙滩区均呈显著正相关。在景观生态风险增加区植被呈退化趋势;在景观生态风险低值区即退耕还林逐年拓展的区域,植被恢复促进景观生态风险的降低。植被恢复对景观生态风险的影响大于城镇化进程。研究成果对于区域发展和生态恢复的综合调控具有一定参考价值。  相似文献   

4.
丹东沿海地区土地利用变化对生境质量的影响研究   总被引:1,自引:0,他引:1  
土地利用变化与生态效应影响研究是全球环境变化的热点。以丹东沿海地区2000—2018年土地利用景观格局为基础,运用GIS技术和In VEST模型对丹东沿海地区生境质量进行定量评估,并结合土地利用变化规律分析了土地利用变化对生境质量的影响。丹东沿海地区作为连接海陆的重要过渡带,研究其土地利用变化对生境质量的影响,对于权衡生态保护与利用土地资源海陆统筹发展具有重要意义。结果表明,(1)耕地、草地及未利用地面积减少,林地、水域和建设用地面积增加,研究末期林地仍占研究区的62.895%,耕地比重由30.565%降至25.735%,是研究区第二大土地利用类型,其次是建筑用地、水域用地,草地与未利用地占比较小。(2)研究期间,生境质量低值区扩大,较低值区缩小,中等级、较高、高值区面积增加。生境质量呈现自东北向西南递减,自沿海向内陆递减的分布特征。(3)研究区生境质量受土地利用变化影响较大,2015—2018土地利用变化与生境质量变化都较为显著,以林地、河渠为主的土地利用类型受人类活动干扰较小,生境质量较好;沿海地区的建设用地对生态环境造成破坏,生境质量较差。(4)通过空间自相关以及热点分析发现,研究区生境质量分布状况存在空间聚集性,聚集效应近年来有分散的趋势,区域间的生境质量差异扩大。  相似文献   

5.
为研究城市海岸带景观类型脆弱度变化特征,揭示海岸带生态脆弱性的时空分布规律,以青岛市2009和2018年遥感影像数据为基础,选取聚集度、分维数倒数、斑块密度、植被覆盖指数和生态适宜度等指标,计算景观类型脆弱度并构建生态脆弱性指数模型,进而对青岛市海岸带生态脆弱性进行评价。结果表明:(1)2009—2018年青岛市海岸带景观类型中,林地、水域和城乡工矿用地景观类型面积呈增加趋势,耕地、草地和未利用地景观类型面积趋于减少,其中,未利用地减少幅度最大(55.18%),城乡工矿用地大幅度增加(58.99%);(2)未利用地和草地的景观类型脆弱度较大,极易受到人为因素和外界环境的干扰,林地景观类型脆弱度次之,水域最小;(3)青岛市海岸带生态良好区和重度脆弱区面积减少,每年比上一年的年平均减少率分别为16.24%和16.83%,而轻度脆弱区、一般脆弱区、脆弱区都呈增加趋势,每年比上一年的年平均增长率分别为0.95%、10.74%和11.24%。生态脆弱性整体呈东部和西南部高、中部低的格局,轻度脆弱是青岛市海岸带的主要脆弱性等级。  相似文献   

6.
基于景观格局的草海流域生态风险评价   总被引:3,自引:0,他引:3  
喀斯特高原景观生态风险变化研究,对建立生态风险预警机制、降低生态风险概率和促进喀斯特高原地区景观格局优化具有重要意义.以2009、2014和2017年的3个时期高分辨率遥感影像为景观格局分析主要数据源,构建景观生态风险评价模型,研究了草海流域景观生态风险时空变化及其空间模式.结果表明,2009—2017年间,研究区景观结构变化较大,耕地和灌草地景观类型面积均呈减少趋势,林地、沼泽地、建设用地和水域景观类型面积呈增加趋势.2009—2017年间,高和较高生态风险主要分布在流域北部、东北部,中生态风险主要分布在流域东部、西南部,较低生态风险主要分布在流域南部-西部-西北部段,低生态风险主要集中在湖泊及其周围.2009—2017年间,研究区生态风险以低和较低生态风险为主,占流域面积的72.38%以上;2009—2014年间,高生态风险、较高生态风险和中生态风险的面积呈不同程度增加;2014—2017年间,高生态风险、较高生态风险和中生态风险的面积明显减少.2009—2017年期间,研究区景观生态风险存在显著空间正相关,"高-高"和"低-低"聚集是主要的生态风险集聚模式,"低-高"或"高-低"零散分布.  相似文献   

7.
探究景观格局变化及其驱动力是自然保护区景观格局优化、生态保护和建设管理的重要基础。基于广东省丰溪省级保护区2010、2015、2020年的遥感影像数据,借助Arc GIS、ENVI和FRAGSTATS等软件,通过分析动态度、转移矩阵、景观格局指数等探究景观格局的演变规律,并对其驱动力进行探讨。结果表明,1)保护区土地利用类型以林地为主,约占区域总面积的61.9%;2010-2020年间,林地和水域持续减少,耕地和建设用地持续增多,土地利用综合动态度由1.87%减少至1.75%,整体呈减少趋势,景观格局变化趋于平稳。保护区林地的转出量最大,转出面积有2.03×103 hm2;草地的转入量最大,转入面积有1.28×103 hm2;耕地的转入面积约为转出面积的3.3倍,整体增加约800 hm2。2)林地是保护区的主体景观,其斑块优势度最大,破碎度最小,景观连通性较好。2010-2020年间,耕地、建筑用地和水域的破碎度降低,连通度升高;草地的破碎度最大,连通性能减低。2020年,景观...  相似文献   

8.
基于LUCC的密云水库上游流域人为干扰动态变化   总被引:2,自引:0,他引:2  
人类活动对地球表层的干扰所造成的影响和变化正日益增强,正确估计人为干扰对现在和未来环境变化影响的程度以及叠加在自然变化上的幅度具有重大意义。选择生态系统服务功能突出、人地关系矛盾剧烈的北京市密云水库上游流域为研究区,基于2005和2013年2期SPOT 5遥感影像数据,借助人为干扰度指数和GIS空间分析方法,从结构变化、类型变化、空间变化和人为干扰度4个角度,探讨密云水库上游流域人为干扰时空动态分异及变化特征。结果表明:(1)2005—2013年间,密云水库上游流域景观以林地-耕地-园地格局为主导,土地利用变化速度加快,人类活动对土地资源的集约利用强度增加。(2)农业结构调整、土地开发整理、经济发展、建设占用以及生态恢复工程是流域内土地利用与土地覆被变化的主要因素;旱地、园地、有林地、灌木林地面积减少,菜地、草地、疏林地、建设用地和滩地增加,森林生态系统质量呈退化趋势;工矿用地和未利用土地主要转为园地、林地和草地,生态恢复工程略显成效。(3)半干扰是流域内主导的干扰景观类型,全干扰和半干扰型景观面积增加,无干扰型景观面积减少,流域内景观的人为干扰强度呈现由深山区向浅山区逐渐增强的空间分异特征,交通道路沿线及人类活动活跃地区是人为干扰强度增加最显著的地区。  相似文献   

9.
根据开都孔雀河流域1990、2000和2010年Landsat-TM影像数据和中巴资源卫星数据,结合GIS技术,研究其近20 a土地利用/覆被及景观格局的动态变化特征及变化原因。结果表明,1990—2010年开都孔雀河流域土地利用类型变化最明显的是耕地、建设用地和水域湿地,年均变化率分别为4.76%、-1.61%和1.45%,而未利用地年均变化率最小,为-0.29%。草地面积的8.71%和未利用地面积的1.73%均转变为耕地,草地面积的2.41%转变为水域湿地,建设用地的54.74%转变为未利用地;林地面积在2000年前减少,主要转变为耕地,而在2000年后耕地和草地转变为林地,又使其面积增加。整个研究时段耕地、水域湿地分布趋于集中,逐渐成为研究区的优势景观类型;在2000年前林地分布破碎化明显,在2000年后趋于集中;1990—2010年,草地面积减少,景观优势度下降,景观形状趋于简单。认为人口的不断增加以及经济利益的驱动使得流域内耕地面积不断扩大,从而严重危及流域的生态环境安全。  相似文献   

10.
草海湿地区域景观生态脆弱性时空变化特征   总被引:3,自引:0,他引:3  
生态脆弱性作为生态学研究的热点和重要方向,对区域环境保护、土地利用规划及可持续发展具有重要指导意义。利用3S技术,基于景观生态学理论和方法构建景观脆弱度指数,从整体上对贵州草海湿地区域的景观类型及景观脆弱度变化特征进行分析。以研究区53个行政村作为评价单元,从较小尺度和可控单元出发,研究草海区域景观生态脆弱性变化及其时空分异特征。结果表明:(1)1995—2015年间草海湿地区域林地、草地面积减少较多,建设用地和水域面积持续增加;研究区景观脆弱度指数范围逐渐变大,但各景观类型脆弱度指数排序不变,由大到小依次为草地、林地、耕地、建设用地和水域。(2)1995—2015年间处于重度脆弱区、中度脆弱区的行政村主要集中在研究区四周,且数量有所减少;处于一般脆弱区和潜在脆弱区的行政村集中在研究区中部区域,数量有所增加。(3)各行政村的景观生态脆弱性水平呈现显著自相关,但景观生态脆弱性接近的行政村在空间上表现出的聚集性有所降低。(4)研究区行政村的景观生态脆弱性空间分布由低-低聚集类型转为高-高聚集类型,大部分行政村的空间自相关属于"不显著"类型,并没有出现低-高聚集区。  相似文献   

11.
提升区域的碳汇能力是中国生态文明建设的重点战略方向,是促进经济社会发展绿色转型的重要举措。太行山区是中国华北地区重要的生态屏障,其生态系统拥有良好的碳汇能力。研究太行山区生态系统碳储量时空分异特征及其影响驱动机制,对华北地区落实国家“双碳”工程建设,提升区域释氧固碳能力,乃至全面提升区域生态环境质量具有重要的意义。以太行山区为例,基于2005、2010、2015、2020年太行山区四期土地覆盖及碳密度数据,使用InVEST模型估算研究区碳储量,使用地理探测器探索驱动碳储量空间分异的主要因子,分析驱动机制。研究结果表明,(1)在2005—2020年期间,太行山区的土地利用类型发生明显变化。林地、建设用地土地利用面积增加,耕地、草地土地利用面积减少。耕地和草地主要转化为建设用地,同时也有一部分耕地转化为林地。(2)太行山区碳储总量在1.48×109—1.50×109t之间,整体逐渐增加。从土地类型来看,碳储量占比由大到小依次为:林地、耕地、草地、建设用地、水域、未利用地。林地增加是太行山区碳储量增加的主要原因;(3)太行山区碳储量空间分异主要受地...  相似文献   

12.
洞庭湖区景观生态风险评价及其时空演化   总被引:1,自引:0,他引:1  
洞庭湖是长江流域重要的调蓄湖泊,该区生态系统服务极为重要,也是中国重要的农产品提供区域。科学评估洞庭湖区景观生态风险,揭示其时空演化特征,对探究其生态格局形成过程与机理,降低生态风险胁迫具有重要意义。利用1990、2000、2010、2015年4期LandsatTM多光谱遥感影像,分析了1990—2015年洞庭湖区土地利用演变特征;基于景观生态学理论引入景观生态风险评价模型,将研究区划分为1 165个生态风险评价单元,基于GIS和地统计学,构建景观生态风险指数模型,对洞庭湖区土地利用动态变化下的景观生态风险时空分布特征及空间关联格局进行了研究。结果表明:(1)1990—2015年,洞庭湖区各土地利用类型均发生了显著变化,其中建设用地面积增加较大,年均增幅达3.61%,主要由林地和耕地的转入;(2)1990—2015年,洞庭湖区的景观生态风险状况以较低生态风险、中等生态风险为主,其面积占景观总面积的70%以上。但高生态风险区域呈增加趋势,面积占比由1990年的4.95%,增加到2015年的8.70%;与2010年相比,中等风险区面积由46.09%降到2015年的34%,这期间低生态风险区和较低生态风险区的面积明显增加;(3)研究区各时期景观生态风险呈现正的空间相关性,且在空间上趋于集聚。景观生态风险分布与区域地理环境具有较高的耦合性,人类活动的集聚程度与风险的空间集聚程度是对应的。以湖滨地带为核心区和远离湖滨区为两大群体,高值-高值(H-H)风险聚集区主要分布湖滨周边,低值-低值(L-L)风险聚集区主要分布于远离湿地核心区的县市区域。研究结论以期为洞庭湖区生态保护与环境整治,调控生态风险提供参考。  相似文献   

13.
景观格局是由自然或人为因素形成的一系列大小、形状各异、排列不同的景观要素共同作用的结果。本文利用2005年3月份MODIS卫星数据,经过一系列处理后得到江苏省景观图,利用决策数分类法将研究区分为5种景观类型:耕地、城镇用地、水体、自然植被、其它用地;选取斑块数、斑块密度、破碎度、斑块分维数、多样性指数、优势度、均匀度、聚集度、分离度等景观指数,从斑块、斑块类型和景观3个层次上对研究区景观格局进行定量分析。结果表明:江苏省景观斑块共57030块,总面积为1018,4275万km2;在5种景观类型中,耕地和城镇用地是主要景观类型,占全省总面积的68.6451%;江苏省景观多样性程度不高,且斑块空间形状较简单。  相似文献   

14.
河岸带是水陆生态系统之间的过渡带,对保护水体和改善生态环境质量具有重要作用。研究河岸带的土地利用和景观格局演变规律,对于维持城市河流生态系统健康、改善城市河流污染具有重要意义。以云南省高原湖泊城市昆明为研究对象,基于1996年、2007年和2015年3期Landsat TM遥感影像,在ArcGIS软件支持下,以城市河流两侧500 m范围作为研究区域,提取土地利用信息,借助土地利用类型特征参数、景观格局指数等对河岸带的土地利用和景观格局演变进行分析,结果表明:(1)从时间变化来看,20年来河岸带的土地利用类型发生了明显的改变,耕地的减少和建设用地的增加是主要的变化趋势,尤其自2007年以来,建设用地的动态度由前一阶段的4.66增长为5.53,表明建设用地呈现越来越快的增长趋势;(2)从空间变化来看,1996年耕地是河岸带各缓冲区的主要用地类型,至2007年,建设用地在各缓冲区的占比已跃居首位,到2015年,建设用地持续快速增长,占比达到各缓冲区的60%左右,成为各缓冲区的景观基质,河岸带的景观结构发生明显改变;(3)从景观类型水平的指数分析来看,20年间,耕地的最大斑块指数(LPI)快速减少且分离度指数(SPLIT)急剧增长,表明其分布呈现破碎化发展;绿地斑块的完整性有所提高,建设用地面积快速增加并呈大面积连片分布的趋势,这与土地利用变化分析结果一致。基于对不同缓冲区景观格局指数的分析,得出昆明市河流生态廊道功能的影响范围为距离河道200 m宽的缓冲区范围,200 m以外的区域景观连接度降低,破碎化程度加剧。  相似文献   

15.
珠江三角洲地区高速公路沿线景观格局变化研究   总被引:2,自引:0,他引:2  
根据1995年和2005年土地利用数据,借助地理信息系统软件和景观生态学方法研究珠江三角洲区域发展轴线——广深高速公路和广珠高速公路运营前后沿线区域的土地利用和景观格局变化,并进一步对比分析相同功能等级道路的生态效应差异。结果表明:(1)高速公路运营后,沿线区域建设用地增加迅速,耕地大量流失且迅速向建设用地转化;(2)公路的建设运营使得沿线区域人工干扰强度加大,景观斑块密度增加、面积减少、斑块形状趋向复杂,景观破碎化加剧;(3)道路途径区域的社会经济条件以及原有的景观格局特征是导致道路生态破碎效应差异的重要原因;(4)公路对自然景观的影响效应远大于其对人工景观的影响效应。  相似文献   

16.
河南省黄泛区农业景观格局动态及驱动因素   总被引:2,自引:0,他引:2  
以1987年和2002年TM影像及1997年河南省土地利用现状图为主要数据源,采用斑块密度、景观多样性指数、破碎化指数、优势度和分布质心5个景观格局指数定量分析近20 a来河南省黄泛区农业景观格局动态变化及其驱动因子.结果表明,1987-2002年河南省黄泛区农业景观动态变化显著,各景观类型面积、斑块数量和斑块密度都发生了较大变化.水浇地、旱地和建设用地景观类型面积增加,是该地区占优势的景观类型,2002年占全区总面积比例为83.51%.该区人类活动强度较大,景观多样性和破碎化指数增大,优势度减小,说明各景观类型所占面积比例差异正在缩小,区域整体景观破碎化程度仍在增加.农业发展、人口增长、经济发展、人民生活水平提高和城镇化进程加速等是该区农业景观格局变化的主要驱动因素.  相似文献   

17.
以五台山地区为研究区,基于生态敏感性-生态恢复力-生态压力度模型,较为全面地选取11个指标构建五台山地区生态脆弱性评价指标体系,以GIS技术为支撑,应用空间主成分分析、重心迁移、差异性分析等方法,对五台山地区2000、2005、2010、2015年4期的生态脆弱性进行评价,并进一步分析其时空变化、重心迁移及差异性特征。结果表明:(1)五台山地区2000—2015年生态脆弱性指数的均值为4.83,整体处于中度脆弱状态,各等级所占的面积虽处于波动状态,但仍以轻度和中度脆弱性为主,重度脆弱区域分布在北部和东北部地区;(2)2000—2015年生态脆弱指数的重心均落在繁峙县境内,且变化幅度较小;(3)2000—2015年3个县级行政区均属于中度脆弱性状态,从高到低依次为繁峙县、五台县和代县;(4)不同土地利用类型的生态脆弱性差异显著。耕地、未利用地和城乡建设用地的生态脆弱性总体水平偏高,处于中度脆弱状态;而林地和草地的生态脆弱性较低,处于轻度脆弱状态。对该地区的生态脆弱性进行定量研究,不仅能正确认识该地区生态脆弱的程度,而且能掌握其空间分布差异,从而进一步对生态脆弱区进行治理,促进资源的有效利用和生态系统的可持续发展。  相似文献   

18.
生态安全格局构建是快速城市扩张背景下缓解生态保护和土地开发矛盾、保障区域生态安全的有效途径。生态源地识别与空间阻力面构建一直是生态安全格局构建中的技术难点。以粤港澳大湾区为研究对象,利用1980、1995和2010年3期遥感影像数据,分析土地利用景观格局变化,提取生态源地,采用最小累积阻力模型判定生态廊道,从而综合构建生态功能网络。结果表明:(1)研究区土地利用景观格局变化剧烈,林地始终保持较大优势度,耕地、林地的斑块结合度和聚集度不断下降,破碎化程度有所增加,景观整体趋于不稳定。(2)湾区生态源地总面积约为2 638. 3 km~2,潜在生态廊道长约705 km,整体生态环境有所改善,但是珠江口东部地区生态连通性差。(3)基于保障自然生态系统的需求,生态源地建设应以北部连绵山体生态屏障和南部海岸等为基底,建设高质量的连片森林,统筹协调跨区域的山林屏障带建设,改善生态源地连通度并提升区域生态系统服务,建设世界级森林城市群。  相似文献   

19.
以海南岛为研究区域,基于2010年和2020年两期土地利用数据,利用CA-Markov模型预测了3种模拟情景下2030年和2040年土地利用空间分布格局,定量分析了海南岛2010—2040年土地利用变化及其生态环境效应。结果表明,(1)2010—2020年间,海南岛土地利用类型以林地为主,面积占比达62%以上,其次耕地。该时期土地利用类型转移的主要特征是林地、耕地、水域转为建设用地,建设用地面积增加483.55 km2,增长率高达53.27%。(2)利用CA-Markov模型模拟3种情景下海南岛2030年和2040年的建设用地均呈现较大幅度的增长,其中耕地保护情景模拟下的耕地面积最多,生态保护情景模拟下的林地和水域面积有所增加。(3)海南岛生态环境质量呈现出中部高、四周低的阶梯式空间分布格局,2010—2040年间海南岛生态环境质量逐期降低,由2010年的0.606 4下降至2040年的0.587 8。在3种模拟情景中,生态保护情景下的海南岛生态环境质量最高,其次是耕地保护情景,基准情景最差。(4)2010—2040年间,林地转为建设用地是导致海南岛生态环境质量恶化的主要类型和主导因素,耕地转为林地则是致使未来生态环境质量改善的主要土地利用变化类型,生态贡献率超过70%。总体而言,海南岛生态环境质量恶化的趋势大于生态环境质量改善的趋势。  相似文献   

20.
基于CA-Markov模型土地利用对景观格局影响辨识   总被引:1,自引:0,他引:1  
采用CA-Markov模型方法,研究挠力河流域不同时段土地利用对湿地景观格局干扰强度的差异性.结果表明:(1)流域湿地景观格局变化显著.1967-2000年,挠力河流域湿地类型景观最大斑块指数逐渐减小,耕地和水域斑块的面积周长分维数逐渐减小,而沼泽和居民地的散布与并列指数变大,草甸和耕地聚集度指数增大;湿地景观破碎化程度增加;湿地景观空间分布逐渐由流域周围高地势区向中心河流廊道退缩.(2)流域土地利用变化明显.耕地和居民用地数量迅速增加,耕地趋于集中连片分布,土地利用呈由自然湿地逐渐转化为耕地的变化过程,耕地成为主导的土地利用类型.(3)流域土地利用强度存在阶段差异.1967-2000年,流域内人类活动对湿地景观的干扰强度逐渐加剧,人为景观面积比例由1967年的26.5%上升至2000年的67.8%.土地利用影响强度指数由1967年的1.690上升至2000年的2.394.2000年模拟的流域土地利用影响强度比2000年实际土地利用影响强度小,由此可以识别1983-2000年人类对湿地景观的干扰强度比1967-1983年有所增大;利用居民地和沼泽面积的变化可以识别流域人口增长和排水活动对土地利用影响的时段差异性.  相似文献   

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

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

京公网安备 11010802026262号