首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 187 毫秒
1.
土壤水分的高时空分辨率和高精度评估对干旱监测具有重要意义。为探究我国内蒙古荒漠草原区土壤水分遥感反演最优模型,基于Landsat和MODIS数据进行改进型自适应反射率时空融合(Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model,ESTARFM),结合下垫面因子、地形因子、气象因子、植被因子等多要素环境因子,通过极限学习机(Extreme learning machine,ELM)和随机森林(Random forest,RF)两种方法构建土壤含水率反演模型,并与Landsat(未进行融合)构建的土壤含水率反演模型进行对比,最终筛选得到最优土壤含水率反演模型,并对研究区不同土地利用类型土壤含水率分布特征进行应用分析。结果表明:归一化植被指数是土壤含水率环境因子相关分析中最重要的预测因子(0~10、10~20、20~30 cm土壤深度处R2=0.85、0.82、0.79),其次为降水量(R2=0.73、0.68、0.71)、高程(R2=0.71、0.70、0.71)、水体指数(R2=0.69、0.69、0.68)、归一化盐分指数(R2=0.68、0.67、0.65)。与未进行时空融合所构建的模型相比,利用ESTARFM时空融合所构建的模型精度均有所提升,考虑ESTARFM时空融合时,ELM模型的R2RMSEMAE分别为0.89、6.58%、3.93%,RF模型的R2RMSEMAE分别为0.78、7.25%、4.95%;未考虑ESTARFM时空融合时,ELM模型的R2RMSEMAE分别为0.75、7.37%、5.24%,RF模型的R2RMSEMAE分别为0.71、7.48%、5.30%。表明ELM模型比RF模型的土壤含水率反演效果更好,且ELM-ESTARFM为土壤含水率反演最优模型。在此基础上,运用改进后的ELM-ESTARFM遥感反演模型监测了乌审旗全域土壤含水率,发现研究区北部和西北部的土壤含水率较高,南部地区的土壤含水率较低;对于不同土壤深度,土壤含水率由大到小依次为耕地、林地、草地、沙地,耕地区域0~10、10~20、20~30 cm土层含水率分别为18.92%、19.34%、21.84%,林地为11.80%、11.87%、12.40%,草地为10.97%、11.02%、12.22%,沙地为5.07%、5.35%、5.67%。  相似文献   

2.
为明确秸秆还田对东北黑土水分特征及物理性质的影响,设置秸秆覆盖还田(FG)、秸秆翻埋还田(FM)和秸秆不还田翻耕(FD)3个处理,测定土壤含水量、水分特征曲线、容重、硬度、土壤三相比及结构稳定性等参数。结果表明:(1)秸秆覆盖还田可显著提高春季耕层(0~30 cm)土壤含水量,较秸秆不还田翻耕处理增幅为11.17%~150.84%;不同处理耕层土壤在水吸力中吸力段土壤含水量变化曲线平滑,秸秆覆盖还田处理具有较高的土壤持水性。(2)秸秆还田能显著提高土壤水分有效性,与秸秆不还田翻耕处理相比,秸秆覆盖还田处理0~10 cm土层土壤田间持水量提高4.85%~11.03%,土壤凋萎系数提高10.85%~18.00%;秸秆翻埋还田处理0~10 cm土层土壤重力水增加9.65%~80.73%。秸秆翻埋还田提升了土壤供水能力,土壤比水容量较秸秆不还田翻耕处理增加4.8%~10.0%。(3)与秸秆不还田翻耕处理相比,秸秆还田降低了收获后土壤紧实度,降低幅度为0.18~0.31 MPa;秸秆覆盖还田增加表层土壤容重,降低土壤孔隙度,促进三相结构趋于合理,显著增加土壤结构稳定性。(4)皮尔森相关分析表明,三相比R值与结构距离(r=0.73*)、土壤容重(r=0.70*)相关性显著,在一定范围内三相比R值的增加有利于改善并促进土壤结构稳定。综上可知,东北黑土农田实施秸秆还田是提高春季土壤含水量、增强土壤持水性、提升土壤供水能力、调节土壤紧实性、调控土壤三相比、改善土壤结构和提高土壤宜耕性的有效措施。  相似文献   

3.
基于Van Genuchten模型测定土壤水动力学参数,采用定点监测方法,在沟道中不同位置以及对照坡面进行土壤水分观测,分析了延安市典型治沟造地项目沟道造地土壤水分的时空变异特征,阐明治沟造地工程对沟道土壤水分的影响。结果表明:(1)土壤水力学参数在沟道土层深度为40 cm附近发生了显著改变,0~40 cm土层土壤容重1.12 ~1.25 g·cm-3,导水率达到40 mm·min-1以上,入渗能力强,同时饱和含水率较大,40 cm以下土层土壤容重在1.5 g·cm-3左右,导水率在1.25~1.41 mm·min-1之间,入渗速率明显减小;(2)沟道土壤水分显著大于对照坡面,其季节变化稍滞后于降水的季节变化,整个生长季在15.76%~21.91%间波动,高出对照坡面5%左右,垂直分布随土层深度的增加而增加,表层最低,为15.07%,160 cm土层最高,为22.84% ,深层土壤含水量优势更加显著;沟道土壤水分变异系数在0.131~0.234之间,相比较坡面,沟道表层以下土壤水分存在较强的空间异质性,100 cm以下土层变异系数在0.2左右,土壤水分变化活跃;(3)沟口土壤含水量显著高于沟头,土层深度40 cm以下的土壤含水量长期达到甚至超过田间持水量。通过治沟造地工程水分综合调控体系,沟道能够为作物生长提供水分充足的生境条件,提高水资源利用效率的同时改善农业生态环境。  相似文献   

4.
为构建适用于干旱区膜下滴灌条件的土壤水盐动态分布和棉花生长模型,基于2020—2021年的田间试验,经过对SWAP模型的土壤、土壤水力功能和作物生长等模块进行率定和验证,对灌溉水矿化度为1、2、3、4、5、6 g·L-1时的土壤水盐分布特征、作物生长过程和干物质累积分配进行数值模拟。结果表明:(1)土壤含水率与土壤含盐量的模拟精度以20~100 cm土层较好,0~20 cm土层模拟精度较差,其中土壤含水量的模拟效果优于土壤含盐量;随着灌溉水源矿化度的增加,土壤含水率和含盐量的模拟误差逐渐变小。(2)不同矿化度水源膜下滴灌棉花叶面积指数模拟效果较好(R2=90.72%,RMSE=0.35 cm2·cm-2NRMSE=8.73%,IOA=0.98)。(3)不同矿化度水源膜下滴灌棉花茎干物质累积量模拟效果较好(R2=89.08%,RMSE=6.12 g,NRMSE=23.16%,IOA=0.96)。研究结果表明,SWAP模型可以较好地对不同矿化度水源膜下滴灌的土壤水盐动态分布和棉花生长过程进行模拟。  相似文献   

5.
基于HYDRUS模型的暗管排水水盐运移模拟   总被引:4,自引:0,他引:4  
为了探索暗管排水条件下膜下滴灌农田的水盐运移规律,本文设计了埋深1 m,间距4 m的暗管排水模型试验,研究分析灌水淋洗过程中土壤水分和盐分的动态变化规律,并利用HYDRUS模型对暗管排水条件下的水盐运移规律进行数值模拟分析与验证。结果表明:经过3次灌水淋洗后,表层0~20 cm土层内盐分含量下降至2 g·kg-1,达到了非盐化土水平,20~60 cm土体内上层土壤脱盐效果高于下层,总盐分含量下降至8 g·kg-1以下;经过实测值与模拟值的验证,土壤盐分和水分的均方根误差RMSE最大分别为0.632和1.324,决定系数R2最小分别为0.992和0.906,说明模拟结果与实测结果吻合度较好,HYDRUS模型能够较好地模拟暗管排水过程中水盐运移规律。通过模型模拟6次灌水(共90 d)后暗管排水条件下不同土层深度的水盐动态变化特征,模拟结果表明,0~40 cm土层内盐分含量下降至2 g·kg-1,40~80 cm土层内盐分含量下降至4 g·kg-1左右,基本达到轻度盐化水平;距离暗管不同间距处的土壤剖面盐分含量呈波动变化,距离暗管越远,土壤剖面含盐量越大,盐分含量在0~8 g·kg-1范围内变化。  相似文献   

6.
通过蒸发装置模拟试验,将雨水引入到表土层以下不同深度,从而改变土壤接受雨水的深度,测定土壤蒸发量,分析土壤蒸发强度、蒸发量、蒸发率和蒸发量的动态变化,确定改变雨水接受深度对土壤蒸发规律的影响。结果表明:雨水接受深度增加,土壤蒸发强度大幅度减少,雨水接受深度(H)为30 cm时,降水量(P)为30、20、10 mm处理的第1天蒸发强度分别为1.6、1.3、0.9 mm·d-1,仅为地表接受雨水处理蒸发强度的21.02%、15.85%和12.98%。土壤蒸发量与雨水接受深度和降雨量大小有关,随时间的延长,相同接受深度条件下,蒸发量随雨量增大而增大。地面接受雨水的土壤蒸发量与时间呈对数相关(R2=0.9959),雨水接受深度为30 cm则呈线性相关(R2=0.9924)。与地面接受雨水相比,地表以下土层接受雨水可降低土壤蒸发损失,随深度增加,土壤蒸发率降低,P10H0、P10H10、P10H20和P10H30处理18 d蒸发率分别为141.14%、109.03%、49.16%和24.75%;土壤水分的蒸发率也与降雨量有关,雨水接受深度为20 cm,降雨量分别为10、20、30 mm时,6 d内土壤水分蒸发率分别为26.76%、13.57%和9.71%%。  相似文献   

7.
为研究不同苹果种植模式下土壤养分变化情况,在陕西省千阳县2种不同种植模式(矮化和乔化种植)的苹果园区采集土样,通过测定0~100 cm土层土壤中的全氮、全磷和有效磷含量,分析和比较其土壤氮磷养分含量差异及剖面分布特征。结果表明:在0~100 cm深度的土层中,除全磷外,矮化苹果园土壤各养分含量显著低于乔化苹果园(P<0.05),矮化富士、嘎啦苹果园和乔化富士苹果园全氮含量分别为0.60、0.63、0.76 g·kg-1,全磷含量分别为0.58、0.56、0.63 g·kg-1,有效磷含量分别为6.98、3.48、51.00 mg·kg-1;矮化苹果园全氮、全磷以及有效磷含量主要在0~10 cm的表层土壤中聚集显著(P<0.05);土壤养分含量与土层深度相关性显著(P<0.05),矮化苹果园和乔化苹果园土壤养分含量与土层深度的变化过程分别符合幂函数和对数函数的变化趋势。  相似文献   

8.
为改善煤矸石山土壤水分和养分状况,提高景观生态型草本植物的水分利用效率,以蜀葵和沙打旺为试验材料,采用三因素五水平二次回归通用旋转组合设计,通过盆栽模拟试验研究了土壤相对含水量、施氮量、施磷量三者耦合对两种景观生态型草本植物叶片水分利用效率的影响,并通过模型寻优得出了最优水肥方案。研究表明:三因素对蜀葵叶片水分利用效率的影响顺序为土壤相对含水量>施氮量>施磷量,三因素对沙打旺叶片水分利用效率的影响顺序为土壤相对含水量>施磷量>施氮量;二因素耦合对蜀葵水分利用效率的影响不显著(P>0.05),但对沙打旺水分利用效率有显著影响(P<0.05),其顺序为水氮>水磷>氮磷;土壤水分条件是影响蜀葵和沙打旺叶片水分利用效率的首要因素,施氮对蜀葵更重要,施磷对沙打旺叶片水分利用效率的影响更大,蜀葵的水肥耦合优化方案为土壤相对含水量37.6%~40.9%,施氮量(纯N)153.1~178.5 kg·hm-2,施磷量(P2O5)89.4~120.4 kg·hm-2;沙打旺的水肥耦合优化方案为土壤相对含水量69.0%~70.6%,施氮量(纯N)61.8~94.6 kg·hm-2,施磷量(P2O5)109.9~170.1 kg·hm-2。  相似文献   

9.
针对Hydrus-2D软件在红壤区涌泉根灌土壤水分运移模拟的适用性问题,依据非饱和土壤水动力学理论,并结合红壤区涌泉根灌土壤水分运动特征建立了涌泉根灌土壤水分的入渗模型,利用Hydrus-2D软件对模型进行求解,并对湿润锋运移距离以及土壤含水率的模拟值和实测值进行了对比验证。结果表明:在灌水结束时,Hydrus-2D软件对竖直向下方向湿润锋的模拟值和实测值之间相对误差为5.21%,水平方向湿润锋的模拟值和实测值之间相对误差为-7.28%,且湿润锋模拟值和实测值的相关系数(R2)均大于0.980,RMSE均在1.300 cm以内,F检验P值也均大于0.05;在灌水结束时,距离灌水器不同距离处土壤含水率剖面分布的模拟值和实测值基本一致,均表现为随着土层深度的增加而先增大后减小,在距离灌水器不同位置处,Hydrus-2D软件对剖面土壤含水率的模拟值和实测值之间的相对误差均在±10%以内,且土壤含水率的模拟值和实测值相关系数(R2)均大于0.990,RMSE在0.030 cm3·cm-3以内,F检验P值也均大于0.05。说明模拟值和实测值具有较好的一致性,模拟结果可为红壤区涌泉根灌系统的合理设计及运行提供依据。  相似文献   

10.
为了提高宁夏扬黄灌区水资源利用效率,促进玉米滴灌技术推广应用,在位于西北黄土高原的宁夏盐池县扬黄灌区开展了玉米滴灌制度研究。试验设计了单因素施肥水平2 131~3 847 m3·hm-2共9个梯度灌水总额,灌水次数10~11次,单次灌水定额150~450 m3·hm-2。通过对玉米生育期灌水量、土壤含水量、产量等指标的监测,分析了玉米生育期土壤含水量变化规律、各生育阶段耗水量强度、不同灌水定额下土壤水蓄存情况、年灌水量与产量的关系。结果显示:灌水前后0~40 cm土壤含水量有显著差异,40~100 cm差异不显著;单次灌水定额为225~375 m3·hm-2时灌溉水在土壤中得到有效蓄存,而单次灌水定额为450 m3·hm-2时,灌溉水有深层渗漏和侧渗现象。灌水量与玉米籽粒、青贮产量均呈正相关,R2分别为0.75和0.82,,在总灌水量为3 402、3 628 m3·hm-2和3 847 m3·hm-2时,产量分别为14 062、11 717、12 889 kg·hm-2,籽粒和青贮产量均随灌水量增加而减少,灌溉水生产效率分别为2.46、2.01、2.14 kg·m-3;灌浆期、抽雄期、拔节期是玉米生长的关键需水期,灌水定额为300~375 m3·hm-2;盐池县扬黄灌区适宜滴灌次数10~11次,灌水总额3 405~3 660 m3·hm-2。  相似文献   

11.
以6 a生灵武长枣为研究对象,设置微孔渗灌不同埋设深度D10(10 cm)、D20(20 cm)、D30(30 cm)和铺设半径R30(30 cm)、R40(40 cm)、R50(50 cm)的2因素3水平试验,分析了微孔渗灌不同环形布设方式对灵武长枣土壤水分分布、产量及水分利用效率的影响。结果表明:随着微孔渗灌埋深的增加,土壤含水率最大分布范围和湿润锋均向下移动,土壤含水率垂直方向最大分布范围可达到10~55 cm。随着微孔渗灌铺设半径的增大,土壤含水率最大分布范围和湿润锋均向外移动,并且各处理的最大湿润区域土壤含水率值均在20%以上。R40D20处理下枣树的枣吊长度、每吊开花数、每吊坐果数、坐果率、单果重、单株结果数、产量、水分利用效率均最大,分别为31.24 cm、66朵、39个、59.70%、19.20 g、444个、7 666.82 kg·hm~(-2)、2.83 kg·m~(-3),并且均与其他处理存在显著差异(p0.05)。因此,微孔渗灌埋深为20 cm、铺设半径为40 cm是灵武长枣最佳的铺设方式。  相似文献   

12.
为了解雨水集聚深层入渗(RWCI)系统土壤水分的入渗规律,设置不同灌水量(10 L、21 L和36 L)和RWCI设计坑深(40 cm和60 cm)的室内土箱试验,观测不同灌水量与不同水头变化情况土壤含水率变化和土壤湿润锋在径向和垂直方向上运移过程,依据非饱和土壤水动力学理论,建立HYDRUS-2D变水头边界条件土壤水分二维入渗模型。通过与实测数据对比,结果表明模型模拟值和实测值具有较好一致性:垂向湿润锋相对均方差(R_E)、平均绝对误差(MAE)和纳什系数(NE)分别为0.019、0.011 cm和0.994,径向湿润锋R_E、MAE和NE分别为0.018、0.851 cm和0.977,土壤含水率R_E、MAE和NE分别为0.188、0.016 cm~3·cm~(-3)和0.916。相比于设计深度为40 cm的RWCI系统, 60 cm RWCI系统在不同灌水量下能够更有效地增加果树根系分布层的土壤含水率,增加土壤水分入渗深度;相同灌水量下RWCI系统设计深度的径向湿润锋分布间无明显差异,而垂直方向的分布具有明显差异;RWCI系统在相同的设计深度下,随着灌水量增大湿润锋在垂向与径向的运移距离差异逐渐增大。  相似文献   

13.
为实现干旱区绿洲土壤含水量的快速、准确监测,利用采集自渭干河-库车河绿洲的84个表层(0~10cm)土壤样本,通过利用电磁感应仪(EM38)将所测解译后数据代替实测土壤含水量数据,将高光谱反射率重采样为Landsat8卫星遥感波段反射率,在选取光谱特征参数、提取敏感波段的基础上,利用偏最小二乘回归(PLSR)方法建立土壤含水量模型,将最优估算模型应用于遥感影像,实现研究区土壤含水量遥感反演。研究结果表明:(1)利用EM38所测水平模式土壤表观电导率与土壤含水量拟合效果最优,能够代替实测土壤含水量进行后续建模分析。(2)相比3种单一的光谱特征指数,利用多种光谱特征指数所建土壤含水量估算模型的建模效果更优,其干、湿各季建模集决定系数R~2大于0.7,均方根误差(RMSE)均小于0.5%,RPD均大于2,能够作为有效手段估算干旱区绿洲土壤含水量。(3)不同季节土壤含水量遥感反演值与实测值决定系数R~2均大于0.6,均方根误差(RMSE)均小于0.6%,显示了较高的预测精度,证明利用电磁感应技术与高光谱相结合能够实现对干旱区绿洲土壤含水量的精准、高效监测。  相似文献   

14.
关中平原干旱遥感监测指数对比和应用研究   总被引:1,自引:0,他引:1  
采用2000-2016年5月份MODIS数据,构建NIR-Red特征空间,对比分析基于该特征空间的垂直干旱指数(Perpendicular Drought Index,PDI)、改进型垂直干旱指数(Modified Perpendicular Drought Index,MPDI)、土壤湿度监测指数(Soil Moisture Monitoring Index,SMMI)及改进型土壤湿度监测指数(Modified Soil Moisture Monitoring Index,MSMMI)这四种干旱监测指数的有效性,并与实测土壤湿度进行相关性分析;最后采用精度最高的SMMI分析关中平原的旱情时空分布特征和规律。结果表明:(1) PDI、MPDI、SMMI及MSMMI均与10cm深土壤湿度存在负相关关系,可决系数R2分别为0.60、0.40、0.64、0.40,表明PDI、MPDI、SMMI及MSMMI均可作为旱情监测指标,且SMMI略优于其它三种监测指数;(2)关中平原东部、中部、西部部分地区旱情严重,西南部地区旱情较轻,且旱情呈年际波动显著的特征;(3) SMMI与月平均气温呈正相关关系区域占75.66%,与月降水量呈负相关关系区域占74.34%,其中通过90%显著性检验区域分别占总面积的27.36%、17.26%,说明降雨和温度不是导致旱情变化的主要影响因子。  相似文献   

15.
在干旱、半干旱地区,地下水是陆生植被生存的重要水源,而传统的地下水位监测方法费时费力,及时获取大尺度高精度的地下水位埋深显得十分重要。在温度植被干旱指数(TVDI)的基础上,分别提取三期Landsat 8遥感数据的归一化植被指数(NDVI)、修正土壤调节植被指数(MSAVI)、比值植被指数(RVI)、差值植被指数(DVI)和增强型植被指数(EVI)和地表温度(Ts),并引入DEM数据对Ts进行地形校正,减少地形起伏对能量二次分配的影响,并选择最佳TVDI反演地下水位埋深状况,结果表明:(1)在传统TVDI的基础上引入DEM进行地形校正,分析校正前后TVDI与地下水位埋深的决定系数,R~2从0.4381提高到0.5053,这说明地形校正能够有效地提高地下水位埋深的反演精度;(2)通过三期影像对比分析了五种不同的TVDI值分别与其对应栅格点的地表土壤湿度和地下水位埋深的决定系数,总体上都是Ts-MSAVI较好,最高R~2分别为0.5547、0.5202;(3)在实测土壤含水量缺失的情况下,可以根据反映土壤湿度高低的因子(TVDI)间接地反演地下水埋深分布状况。  相似文献   

16.
河套灌区农田地下水埋深普遍较浅且年内波动较大,明确不同膜下滴灌条件下深层土壤水分对根区的补给作用及作物根系吸水的响应差异有利于膜下滴灌技术的完善和推广。本研究开展了连续2 a(2017—2018年)的春玉米田间试验,设置3个膜下滴灌灌溉水平,分别控制土壤基质势下限为-10 kPa(S1)、-30 kPa(S3)和-50 kPa(S5)。利用HYDRUS-2D模型模拟0~120 cm深度土壤含水量、根层下边界(100 cm深度处)水分通量和作物根系吸水速率。结果表明,经过率定后的HYDRUS-2D模型对0~120 cm深度土壤含水量模拟结果的根均方差(RMSE)和决定系数(R~2)分别为0.039~0.042 cm~3·cm~(-3)和0.78~0.73,模拟结果可靠。膜下滴灌农田100 cm和120 cm深度处土壤含水量较高且处理间差异不大,说明不同滴灌条件对于100 cm以下深层土壤含水量影响较小;但不同处理显著影响根区下边界的水分通量和根系吸水速率。基质势下限控制水平越低,深层土壤水分对于根区的补给量(毛管上升)越大,S1、S3、S5生育期内累积补给量在31.9~49.6 mm之间。S5处理根系吸水速率较低,根系吸水受到显著抑制,从而造成作物生长指标和产量显著低于S1和S3处理(P0.05);而S1和S3之间籽粒产量差异不显著。综上,在本研究所设置的3个滴灌处理中,S3生育期内灌溉定额为240~300 mm,既较S1显著减少灌水量、提高水分利用效率,又具有较好的根系活力,有效利用深层土壤水分,因此建议该地区春玉米膜下滴灌的灌水下限为-30 kPa。  相似文献   

17.
3种保水措施影响下的风沙土水盐及地温特征   总被引:3,自引:0,他引:3  
对塔克拉玛干沙漠腹地不同保水措施(覆膜、覆沙、喷施Guilspare)和裸地(对照)土壤水分、盐分分布及耕作层温度的变化特征进行了分析.结果表明:夏季采取保水措施能够抑制土壤蒸发,提高土壤的蓄水保墒能力.同时,还能起到抑制土壤盐分表聚的作用,明显地改变了风沙土水盐空间分布,降低表层盐分被淋溶至林木根系层的概率.保水措施...  相似文献   

18.
探究施用生物炭对冬麦田土壤水热因子及土壤呼吸的影响,对生物炭在麦田应用及农田固碳减排有重要的实践意义。2018年10月至2021年6月在关中灌区连续3 a进行了麦田生物炭施用试验,试验设置生物炭施用量水平分别为:0 t·hm-2·a-1(C0)、10 t·hm-2·a-1(C10)、20 t·hm-2·a-1(C20),通过测定小麦生长季的土壤温度、土壤水分、土壤呼吸速率及产量,明确不同施炭量对冬小麦田土壤水热因子及土壤呼吸的影响。各处理生育期内土壤呼吸速率及全生育期CO2排放量存在显著性差异(P<0.05),均表现为C0>C20>C10。生物炭施入增加了生育期内的平均土壤温度,同时显著提高了0~20 cm土壤含水量(P<0.05),并减弱了土壤水分在生育期内的变化幅度。C10、C20处理3 a平均土壤含水量较C0分别增加了17.0%、29.0%。5 cm及10 cm土壤温度能分别解释土壤呼吸变化的54.7%~...  相似文献   

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

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