首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 468 毫秒
1.
树干径向变化对茎液流以及气象因子的响应   总被引:3,自引:0,他引:3  
应用改进的SF-L型热扩散式树液流测定装置和DC型植物生长测量仪,结合北京2007年6-10月的气象数据,探讨了北京山区油松、侧柏树干茎周长动态因素与液流速率和气象因子之间的相关关系,并运用树木径向变化周期分段法将树木径向生长量从径向变化过程中分离开,结果表明:(1)油松、侧柏茎周长变化与液流速率的关系显著(P<0.01),油松为正相关,侧柏为负相关.(2)空气水汽压亏缺、日温差、空气湿度和日照时数与油松、侧柏茎周长变化关系显著,但在回归拟合过程中侧柏的解释因子只有空气水汽压亏缺和日平均风速,且解释程度比较低(0.356).(3)油松、侧柏茎周长生长量的拟合水平较高,且影响气象因子相同,侧柏的液流速率也成为解释因子之一.水汽压亏缺、空气温湿度、日照时数、土壤水势和风速对树木的生长起到决定性的作用.  相似文献   

2.
利用辽宁省阜新县1976-2012年赤松毛虫越冬死亡率数据及同期气象站的温度、降水量、风速、相对湿度、日照时数等资料,采用相关分析、逐步回归分析、主成分分析与逐步回归相结合等方法,确定影响赤松毛虫越冬死亡率的关键气象因子,并建立赤松毛虫越冬死亡率的预报模型,两种方法建立的模型都通过了0.05水平的显著性检验。结果表明,冬季温度是影响赤松毛虫越冬死亡率的主要因子,其次是日照时数和温湿度的配合情况。仅采用逐步回归分析方法(方法Ⅰ)所建模型得到的拟合值与观测值的相关系数为0.75,D2(越冬死亡率模拟值与实际值相对差值的平方)均值为0.58,而采用主成分与逐步回归相结合方法(方法Ⅱ)建立模型得到的拟合值与观测值的相关系数为0.76,D2均值为0.62。利用两种方法对2009-2012年松毛虫越冬死亡率进行模拟预测,其中2009,2010和2012年方法Ⅰ预测准确率在70%以上,方法Ⅱ的预测准确率在80%以上,表明方法Ⅱ所建模型的拟合及预报准确率均高于方法Ⅰ,但由于两种模型的适用条件不同,在各气象因子数据容易获取的情况下,应选用方法Ⅱ,反之则可选用方法Ⅰ进行预报。  相似文献   

3.
基于2011-2013年超级稻岳优6135分期播种田间试验数据,利用抽穗普期-成熟期每5d一次的千粒重动态观测资料并计算相应的灌浆速率,与同时段逐日气温和日照时数进行数理统计分析,建立基于不同温度要素、日照时数及温光组合的灌浆速率和千粒重的模拟模型,探究该水稻品种灌浆特性与各气象因子间的关联以及达到最优灌浆效果的气象条件。结果表明:灌浆速率、千粒重与各温度要素及日照时数间均呈一元二次函数关系,与各温光组合间呈二元二次函数关系。灌浆特性与各温度要素间的相关性优于其与日照时数间的相关性,不同温度要素及温光组合与千粒重的相关性优于与灌浆速率的相关性;日最低气温累积值与日照时数的组合是影响超级稻岳优6135灌浆特性的最重要温光因子。灌浆期积温达900℃·d,日照时数接近300h时,千粒重达到最大,灌浆过程趋于完成,籽粒趋于饱满,籽粒充实度最理想。  相似文献   

4.
石羊河流域近53 a参考作物蒸散量的敏感性分析   总被引:1,自引:0,他引:1  
利用国家气象信息中心提供的地面气候资料日值数据集,基于FAO Penman-Monteith公式计算了石羊河流域4个测站1959-2011年的逐日参考作物蒸散量(ET0)。利用敏感系数法计算了其对平均最高气温、平均最低气温、风速、平均相对湿度和日照时数的敏感系数,并分析了敏感系数的时空变化特征。结果表明,石羊河流域ET0对相对湿度最敏感,其次为风速和气温,而对日照时数的敏感性最低。由于气象要素分布不均,敏感系数的空间差异显著,相对湿度的敏感系数在上游祁连山区形成高值区,同时,气温在该区的敏感系数也相对较大,而风速的敏感系数在下游民勤盆地较大,日照时数的敏感系数在全区无明显差异。各气象因子的敏感系数均存在一定程度的波动,风速的敏感系数冬高夏低,气温和日照时数的敏感系数均为夏季最高,相对湿度敏感系数的绝对值持续上升在秋季达到最大。53 a来,相对湿度敏感系数波动变化,近20 a来其绝对值上升趋势显著,而风速、日照时数和气温的敏感系数无明显变化趋势。  相似文献   

5.
以陕北黄土丘陵沟壑区延河、清涧河、汾川河3条流域作为研究区域,针对该区9个水文站点的6-9月径流量、输沙量的数据,探讨了时空尺度变化对径流与输沙相关性的影响.结果表明二者在不同的时间尺度和空间尺度上其相关程度差异明显.基于径流因子与输沙因子的关系所拟合的曲线模型中,7,8月份的模型拟合优度在整体上高于6月、9月,汛期的模型拟合优度一般低于7,8月份的模型拟合优度.汛期清涧河流域的模型拟合优度高于延河和汾川河流域,区域的模型拟合优度高于延河和汾川河流域,但略低于清涧河流域.  相似文献   

6.
华北地区杨树烂皮病发生发展气象适宜度预报模型   总被引:7,自引:0,他引:7  
利用河北省18个代表站点2002-2008年杨树烂皮病病情资料和相应气象资料,通过相关、回归分析等方法,分析河北省杨树烂皮病发生发展与气象条件的关系。结果表明:降水量、降水日数、相对湿度、温湿系数与杨树烂皮病发生率呈正相关,降水增加、雨日增多、相对湿度和温湿系数增大有利于杨树烂皮病发生发展;日照时数、夏季高温日数、极端最低气温与杨树烂皮病发生率呈负相关,天气晴朗、日照充足、极端气温升高不利于杨树烂皮病发生发展;影响河北省杨树烂皮病发生发展的关键气象因子分别是3月和4月温湿系数、3月下旬和4月下旬日照时数、4月上下旬气温、3月上旬和4月上旬以及上年12月上旬极端最低气温,利用归一化处理后的关键气象因子建立气象适宜度预报模型,通过回代检验,气象适宜度等级与杨树烂皮病实际发生等级相比,模型拟合准确率为91%;对河北省2009年5-6月杨树烂皮病盛发期气象适宜度等级进行预报,准确率为78%,效果较好;3-4月气象条件是影响华北杨树烂皮病发生发展的关键生态环境因子。模型可用于从气象角度对华北地区杨树烂皮病发生发展进行监测和预报。  相似文献   

7.
近40年来,山西省的气候条件发生重大变化。为了评估参考作物蒸散量(ET0)在气候变化条件下的时空变化趋势及原因,利用彭曼公式计算了晋北、晋中和晋南地区的ET0,用t检验法分析了不同区域气候因子及ET0演变趋势,利用偏相关分析法研究了各气候因子与ET0的相关关系及对ET0影响的贡献度。主要结论:晋北、晋中和晋南地区年均日照时数、年均最高气温、年均最低气温、年均风速和年均相对湿度等气象因子发生明显的上升或下降趋势,但ET0变化趋势不明显,仅呈缓慢下降趋势;ET0与年均日照时数、年均最高气温和年均风速均呈显著正相关关系,与年均相对湿度呈显著负相关关系。晋北地区年均最低气温和年均最高气温引起ET0显著上升,贡献度之和为31.36%;年均日照时数和年均风速导致ET0显著下降,贡献度之和为41%。晋中地区上升因子的贡献度之和为46.34%,下降因子为53.66%。晋南地区上升因子贡献度之和为27.78%,下降因子为50.15%。  相似文献   

8.
河南省参考作物蒸散量变化特征及其气候影响分析   总被引:1,自引:0,他引:1  
基于河南省111个气象站1971-2010年逐日平均气温、最高气温、最低气温、相对湿度、风速和日照时数等气候要素资料,应用Penman-Monteith模型计算各站点逐日参考作物蒸散量(ET0),结合数理统计方法,分析近40a来河南省年ET0的时空变化特征,并对其主要影响因子进行探讨.结果表明,Penman-Monteith模型对河南省ET0的模拟能力较强,模拟值与同期小型蒸发皿蒸发量的相关系数r=0.84(P <0.01).近40a,河南省年ET0平均值为796.1mm(±102.2mm,n=4169),在空间分布上,总体表现出北高南低的特征,并以24.7mm·10a-1(P <0.01)的线性倾向率减少,呈明显减少的站点主要分布在34°N以北地区.偏相关分析表明,全省各地(市)年ETo与各气象要素关系密切,除济源外,年ET0均表现出与风速呈负相关且相关系数最大.逐步回归分析显示,年ETo与平均气温、日照时数、风速和相对湿度的关系密切;风速、日照时数和平均气温对年ET0的贡献为正效应,而相对湿度为负效应.近40a,风速减小是导致河南省年ET0呈显著减小的主要原因;但从综合影响看,这是各气象因素综合作用的效果,且各因子的贡献存在区域差异.  相似文献   

9.
基于1956-2015年商丘市气象数据和霜冻害资料,研究黄淮冬麦区商丘市气象因子与霜冻害的关系,以及气温、地温、相对湿度、风速与草面温度的偏相关和多元线性回归关系。结果表明:随着全球气温变暖,商丘市冬春季积温以4.35℃·d·a~(-1)的速率逐年增加;晚霜冻害除受冷空气活动影响外,与前期积温和降水量密切相关,冬春季积温偏高或降水量偏少的情况下易发生晚霜冻害;随着小麦幼穗的发育,其对低温的敏感度增加,且可引起霜冻害的最低温度有逐渐升高的趋势;气温、地表温度、平均相对湿度、平均风速与草面温度呈极显著相关(P0.01),偏相关系数的大小表现为地表温度气温平均相对湿度平均风速,可见,除气温、地表温度外,平均相对湿度、平均风速也对霜冻害的发生及轻重程度起关键作用;各因子与草面温度可用模型表述为Y=0.558ST+0.482AT+0.087RH+1.304WS-12.704,经检验,线性回归方程成立,可通过该模型对草面温度进行监测,并为未安装草面温度传感器的地区提供可靠的冻害评估结果。  相似文献   

10.
基于1981-2010年江西省市(县)级稻瘟病监测资料和同期气象资料,利用符号相关法,分苗瘟、叶瘟、穗瘟分析各时期病害发生程度与气象因子的关系,并分别建立病害发生潜势气象预报模型。结果表明:与病害高峰同期的气象因子之间,早稻苗瘟与4月份的平均气温、雨量、日照时数呈正相关,与相对湿度呈负相关;早稻叶瘟与5月中旬-6月上旬的平均气温、雨量、日照时数呈负相关,与相对湿度呈正相关;早稻穗瘟与6月上旬-7月上旬的平均气温、日照时数呈负相关,与雨量、相对湿度呈正相关。晚稻苗瘟偏重发生的概率随着苗期高湿天数和雨日的增加而增大;晚稻叶瘟与8月中旬-9月上旬的平均气温、日照时数呈负相关,与雨量、相对湿度呈正相关;晚稻穗瘟与9月平均气温、相对湿度呈正相关,与雨量、日照时数呈负相关。通过相关关系和符号相同、相反的概率差值确定各气象因子距平值的量化值和权重系数,计算稻瘟病发生潜势气象指数,建立稻瘟病发生潜势气象预报模型。历史回代检验表明,模型准确率为60%~75%,且对偏重级别的准确率高于偏轻级别。  相似文献   

11.
基于气候信息的喀斯特地区植被EVI模拟   总被引:1,自引:0,他引:1  
该研究以喀斯特地区植被为研究对象,分析各种气候因子与植被指数的相关性及作用机制,在此基础上建立基于气候因子的植被EVI拟合模型,为定量分析气候条件对植被的综合影响奠定基础。结果表明:气候因子对喀斯特地区植被EVI影响显著,植被EVI与水汽压、平均气温、露点温度、最低气温、最高气温的相关性均大于0.8且空间一致性好。除日照时数和风速外,该地区植被EVI对其他气候因子的响应均存在显著滞后性,滞后期约16 d。对植被EVI起直接作用的主要是温度类气候因子,水分类气候因子对植被EVI的直接作用不明显,但通过其他气象因子起了较强的间接作用。根据该地区植被与气候因子的关系建立了2个EVI拟合模型,其中基于同期气候因子的同期模型中入选的气候因子为水汽压(0期)、日照时数(0期)、露点温度(0期),基于同期、前期气候因子的混合模型入选气候因子为水汽压(-1期)、最高气温(-1期)、降水量(-1期)、露点温度(-1期)、日照时数(0期)。分别利用2001-2010年建模数据和2011年非建模数据对2个模型进行了单站点和片区两种尺度的精度验证。验证结果表明,2个模型对整个片区植被EVI的拟合精度高于单站点,且混合模型的拟合精度高于同期模型。2001-2010年同期模型和混合模型的片区拟合R2分别为0.843、0.892,站点拟合R2分别为0.765±0.033、0.801±0.021。2011年2个模型的片区拟合R2分别为0.797、0.873,站点拟合R2分别为0.716±0.073、0.746±0.064。对大多数站点而言,混合模型的拟合精度较高,但是由于2个模型的建模气候因子不同及各个站点植被的EVI与气候因子的综合响应也存在较大差异,同期模型对部分站点植被EVI拟合精度高于混合模型。  相似文献   

12.
基于北疆地区20个气象站点1966-2012年逐日气象数据,采用Penman-Monteith公式计算各站生长季逐日参考作物蒸散量(ET0),利用线性趋势分析、M-K突变检验、因子贡献率、敏感系数等方法,对该地区ET0的时空变化规律及其对气象要素的敏感性进行分析。结果表明,北疆地区生长季ET0多年平均值为939.43mm,47a间以16.30mm×10a-1的速率呈极显著下降趋势(P水汽压>风速>日照时数。水汽压升高、风速减小、日照时数减少对ET0的减小作用超过温度上升对ET0的增大作用,是导致近47a北疆地区ET0总体呈下降趋势的根本原因。  相似文献   

13.
欧阳习军      董晓华      魏榕      龚成麒      吴寒雨     《水土保持研究》2023,30(2):220-229
为了探究青藏高原植被覆盖时空演变特征及其驱动因子,对青藏高原的生态环境保护提供科学依据,基于1982—2015年青藏高原内部及其周边139个气象站点的气象数据和同期的GIMMS NDVI数据,研究了青藏高原生长季植被NDVI的时空变化特征及其与气候因子的响应关系。结果表明:(1)在研究期内,青藏高原生长季NDVI总体呈上升趋势,不同干湿地区生长季NDVI变化趋势有所差异,湿润地区植被退化面积占比相对较大,干旱地区植被改善面积占比相对较大。(2)研究区植被未来总体向改善方向发展,植被未来趋向改善面积占62.25%,趋向退化面积占37.58%。(3)研究区植被对各气候因子的响应存在一定的滞后性,草原、草甸、高山植被和灌丛4种主要植被对气温和相对湿度主要当月响应,对降水主要当月或滞后1个月响应,对日照时数主要滞后3个月响应。(4)气温、降水、相对湿度及日照时数4个气候因子对青藏高原植被NDVI变化的相对贡献率分别为37.19%,27.53%,20.30%和14.97%,其中,气温和降水是湿润/半湿润地区、半湿润地区、大部分半干旱地区及干旱地区植被NDVI变化的主要气候驱动因子,日照时数和相对...  相似文献   

14.
黑龙港流域参考作物蒸散量的时序变化和分形特征   总被引:1,自引:0,他引:1  
利用位于黑龙港流域的武强、深泽、饶阳、晋州、献县5站1957-2009年日最高气温、日平均气温、日最低气温、日平均相对湿度、日平均风速、日照时数资料,采用Penman-Monteith算法,计算各站不同时间尺度ET0,采用线性趋势分析法分析其趋势倾向,并应用滑动R/S分析方法研究该流域不同时间尺度ET0时间序列的分形特征。线性趋势分析显示,各站历史上自1957年以来的ET0年总值的气候倾向率在-33.81-10.79mm.10a^-1,即均呈下降趋势,但变化倾向率不同;各月ET0倾向率在-11.27-2.02mm.10a^-1,大多数为负值,其中5、6月份各站的下降趋势最大;春、夏、秋、冬季ET0倾向率为-15.87-1.30mm.10a^-1,且夏季各站之间的差异较大,尤以饶阳站与其他4站间的差异最大。气候要素倾向率的对比分析表明,5站参考作物蒸散量总体下降趋势的基本特征主要是由于风速下降、日照时数减少、日最高最低气温上升造成的;5、6月份下降趋势明显则主要是由于风速下降、日照时数减少、日平均气温上升减缓、相对湿度上升明显造成的;饶阳站与其他4站差异较大的主要原因在于2、3、11、12月风速下降和日照时数减少以及日最高最低气温上升趋势减缓、相对湿度和气压下降明显,以致这些月份的ET0倾向率大于0,形成秋季和冬季ET0倾向率大于0,年尺度ET0年际变化下降趋势不明显。R/S分析结果显示,5站全年和各季ET0时间序列的Hurst指数均大于0.5,相关系数均在0.98以上,分维数均小于1.5,说明各站全年和各季节ET0时间序列变化趋势在未来一段时间内具有持续性,即ET0在未来将呈较明显的下降趋势,这种趋势在除饶阳外的其他站月均有不同程度的表现。  相似文献   

15.
气候变化背景下内蒙古马铃薯关键生长期气候适宜性分析   总被引:5,自引:1,他引:4  
内蒙古是我国马铃薯主产区之一,在气候变化背景下探讨马铃薯的区域适应性有重要的现实意义。依据内蒙古地区119个站点逐日气象观测数据,采用气候适宜度评估模型,分1961—1990年、1991—2015年2个气候年代对内蒙古马铃薯开花—可收期气候适宜度的空间特征及变化趋势进行分析。结果表明:1)1961—2015年内蒙古马铃薯开花—可收期综合气候适宜度为0.62,气候条件总体优越,且日照适宜度温度适宜度降水适宜度。2)从空间分布上看,温度适宜度由西向东呈"低—高—中"的带状分布,降水适宜度由西向东、由北向南递增,日照适宜度西高东低;综合光温水三要素影响,阴山南北麓东段及大兴安岭南麓区为气候适宜度的高值区,具有适宜生产马铃薯的气候资源优势。3)从年代变化上看,1961—1990年,综合气候适宜度除东部偏北及西部偏北呈微弱的下降,总体呈增加的趋势;1991—2015年大部地区综合适宜度由上升转为下降趋势。4)气候变化背景下,1961—2015年内蒙古马铃薯开花—可收期气候适宜性以0.004 0×(10a)~(-1)的倾向率下降,其中日照适宜度基本持平,降水适宜度、气温适宜度呈不同程度的下降趋势,对马铃薯生长发育不利。因此,马铃薯开花—可收期应注意水热条件的管理,研究结果可为马铃薯气候资源利用与管理提供参考。  相似文献   

16.
黄河上游参考作物蒸散量变化特征及其对气候变化的响应   总被引:11,自引:9,他引:2  
研究参考作物蒸散量对气候变化的响应对于辨析气候要素对蒸散发的影响具有重要意义。该文在验证FAO推荐的Penman-Monteith(P-M)方法在黄河上游地区适用性的基础上,分析了10个气象站点近50a来参考作物蒸散量的变化特征,计算了参考作物蒸散量对4种气候要素的敏感系数及其对气候变化的响应。结果表明:FAOP-M方法在研究区域具有较强的适用性;参考作物蒸散量随海拔升高而减少,主要集中在生长季的3-10月。高海拔站点参考作物蒸散量年值多呈显著增加趋势,低海拔则明显减少,且变化过程不同。气温和风速敏感系数的年内变化分别呈显著的波峰型和波谷型,日照时数的变化不明显,相对湿度的敏感性呈生长季略有增加趋势;年际变化方面,气温的敏感系数呈显著增加趋势,低海拔地区相对湿度的敏感系数呈显著增加趋势。气温、日照时数的增加和相对湿度的降低导致了高海拔站点参考作物蒸散量的增加,高海拔地区的蒸散发主要受气温、日照时数等能量制约;低海拔站点参考作物蒸散量的减少主要受日照时数和风速减小、相对湿度提高的影响,水分条件的限制更显著。该研究对于结合未来气候变化趋势开展黄河上游地区生态管理,促进生态系统的良性发展,具有重要的科学意义和应用价值。  相似文献   

17.
松嫩平原西部生长季参考作物蒸散发的敏感性分析   总被引:26,自引:13,他引:13  
为了预测气候变量扰动引起的区域参考作物蒸散发的变化,该文以松嫩平原西部为研究区,运用FAO Penman-Montieh方程计算了34个站点1951~2001年的生长季参考作物蒸散发,并计算其对气温、风速、日照时数和相对湿度的敏感系数,最后分析了敏感系数的时空变化特征.结果表明:参考作物蒸散发对相对湿度最为敏感,其次是气温,最后是风速和日照时数;5~9月,气温和日照时数的敏感系数呈单峰型分布,在7月份达到最高值,风速敏感系数呈单谷型分布,在7月份达到最低值,相对湿度敏感系数的绝对值表现为持续上升趋势;生长季气候变量敏感系数的空间差异性较大.日照时数、气温、风速等3个气候变量的敏感系数都在研究区西南部形成高值区,而相对湿度则在研究区西南部形成低值区,在研究区东北部形成高值区.  相似文献   

18.
京津冀地区潜在蒸散量时空演变特征及归因分析   总被引:5,自引:5,他引:0  
为了深入认识京津冀地区潜在蒸散量的时空变化特征及其对气候变化的响应,该研究基于京津冀地区23个气象站57 a逐日气象观测资料,应用Penman-Monteith公式计算各站点日潜在蒸散量(ET0),剖析ET0的时空变化特征,运用敏感性分析法定量研究ET0对各气象要素的敏感性及其时空变化特征,定量识别各气象要素变化对ET0变化的贡献。研究结果表明:1)京津冀地区ET0空间分布整体呈由南向北递减趋势(除中部地区的塘沽站、黄烨站与保定站点ET0较高外)。ET0整体呈下降趋势,线性趋势率为-0.92 mm/a。ET0变化趋势空间分布由西北向东南递减,以春季减幅最为明显。2)京津冀地区ET0对相对湿度的最为敏感(-0.44),其次为风速(0.31)、日照时数(0.28)与平均气温(0.26)。随时间推移,ET0对平均风速与相对湿度敏感性整体呈下降趋势,而ET0对平均气温与日照时数的敏感性逐渐增强。敏感性系数空间分布从西北到东南:风速与平均气温敏感性系数逐渐递增,而日照时数与相对湿度敏感性系数逐渐递减。3)风速变化对京津冀地区ET0变化的贡献最大,平均气温次之。风速为主导因素的站点个数随时间呈下降趋势,平均气温与日照时数为主导的站点个数随时间呈上升趋势,说明近年来平均气温与日照时数对潜在蒸散量变化的影响愈加明显,这可能是由于近年来京津冀地区雾霾尤其是冬季雾霾对日照时数、气温与风速的产生一定影响,进而影响ET0。  相似文献   

19.
东北地区参考作物蒸散量对主要气象要素的敏感性分析   总被引:13,自引:1,他引:12  
利用国家气象局提供的地面气候资料日值数据集,通过FAO推荐的Penman-Monteith公式计算了东北地区1961-2008年生长季(5-9月)逐日的参考作物蒸散量(ET0),分析了ET0及主要气象要素的变化趋势,并通过响应曲线、敏感矩阵、敏感系数等方法分析了ET0对气温、日照时数、平均风速、平均相对湿度的敏感性。结果表明:(1)近50a来,东北地区的气温呈极显著上升趋势(P0.01),日照时数、平均风速、平均相对湿度呈极显著下降趋势(P0.01);东北地区生长季平均日ET0在以3.60mm.d-1为平均值、±0.3mm.d-1的范围内波动,总体上比较稳定,最大值出现在2001年(3.87mm.d-1),最小值出现在1990年(3.28mm.d-1);(2)当气温、日照时数、平均风速的变化量从-20%变化到20%时,ET0表现为逐渐增加的趋势,当平均相对湿度的变化量从-20%增加到20%时,ET0则逐渐减小;(3)气温、日照时数、平均风速、平均相对湿度的生长季平均日敏感系数均具有较强的空间分异特性,其中气温变化对ET0的影响最为明显,其次是平均相对湿度,日照时数、平均风速对ET0的影响较小。  相似文献   

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

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