首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 309 毫秒
1.
基于数据同化优化陆面过程模式参数,从物理机制上改进和完善陆面过程模式参数化方案,提高模式的区域适应性和模拟能力,对改善大气环流模式(GCM,General Circulation Model)和区域气候模式(RCM,Regional Climate Model)模拟和气候预测能力具有重要的现实意义.本文发展了一种基于数据同化优化陆面过程模式参数化方案的方法,并以优化非饱和土壤水模型为例,发展基于非饱和土壤水模型和扩展Kalman滤波(EKF,Extended Kalman Filter)算法并结合可变下渗能力模型(VIC,Variable Infiltration Capacity Model)的土壤湿度同化方案,以月为同化窗,利用复合形混合演化算法(SCE-UA,Shuffled Complex Evolution-University of Arizona)极小化目标函数,使模拟与同化的土壤湿度按给定的目标函数度量方式达到最佳拟合,优化非饱和土壤水模型参数:饱和水力传导率、饱和土壤湿度、饱和土壤水势、Clapp和Hornberger常数,得到同化时段1986年模型参数的最优值序列,然后用其改进非饱和土壤水模型,最后利用改进后模型进行1986~1993年的数值模拟试验,评估改进后模型的模拟能力,探讨基于数据同化优化陆面过程模式土壤水参数化方案的方法与实现途径.试验结果表明:基于数据同化优化模型参数,从物理机制上改进和完善模型,提高了模型的模拟能力;优化的模型参数在时间上具有较好的移植性,适应性更强,达到了改进模型的目的,因此,本方法是合理、可行的.本文为基于数据同化优化陆面过程模式参数化方案研究提供了范例,具有较好的借鉴作用.  相似文献   

2.
近地层能量闭合度对陆面过程模式影响   总被引:1,自引:0,他引:1       下载免费PDF全文
大量近地层观测试验表明,利用涡动相关法观测的湍流通量小于近地层可利用能量,即近地层能量是不闭合的,这种不闭合度一般为20%甚至更高.而陆面过程模式是基于地气间能量平衡建立,并且模式中的湍流边界层参数化方案通常根据实际观测的湍流通量来确定,因此能量不闭合必将对陆面过程模式造成一定的影响.本文利用2007年春季SACOL站的近地层观测资料,依据能量守恒将能量不闭合中的残余能量通过波文比分配到观测的湍流通量中,即修正涡动相关法观测的湍流通量使得近地层能量达到平衡;之后分别利用观测和修正的湍流通量,建立了能量不闭合和闭合情形下的湍流参数化方案,借助陆面过程模式SHAW,通过数值模拟和对比分析方法考察近地层能量闭合度对陆面过程模式的影响.研究结果表明近地层能量闭合对陆面过程模式有显著的影响:在陆面过程数值模拟中,当应用近地层能量不闭合的湍流通量形成的湍流参数化方案时,陆面过程模式会明显高估地表长波辐射及土壤温度;但当应用修正湍流通量使得近地层能量达到闭合形成的湍流参数化方案后,在不改变任何地表土壤物理生化属性的情况下,陆面过程模式能较好地模拟地表长波辐射和土壤温度.  相似文献   

3.
沙漠绿洲陆面物理过程和地气相互作用数值模拟   总被引:1,自引:0,他引:1  
利用一个已发展的陆面物理过程参数化方案与大气边界层数值模式耦合, 模拟了半干旱区沙漠绿洲非均匀下垫面的陆面物理过程及其与大气边界层的相互作用过程, 成功地模拟了局地气候效应和地表温度、净辐射、感热和潜热通量特征并与实测资料进行了比较. 给出“绿洲效应”这一自然现象的垂直剖面上更为清晰准确和细致的结构特征, 结果表明: “绿洲效应”具有明显的“冷岛效应”和“湿岛效应”; 它表现为在绿洲区域比戈壁沙漠区域环境温度低、湿度大、湍流动能输送弱, 具有下沉气流而导致与周围戈壁沙漠区域产生水平输送环流. 这些结果对于深入了解绿洲气候的形成和绿洲的维持机理具有重要的意义.  相似文献   

4.
城市地表潜热通量数值模拟方法研究   总被引:2,自引:0,他引:2  
近年来,中尺度气象模式中城市建筑物影响的数值模拟方法日趋完善,城市地表潜热通量数值模拟方法是目前城市地气耦合研究面临的瓶颈问题之一.采用北京325 m气象塔140 m高度处的能量平衡观测资料,用陆面-单层城市冠层耦合模式进行了一年连续的离线模拟.根据模式对城市下垫面潜热通量模拟的系统性误差较大的特点,基于观测资料研究了城市下垫面潜热通量的数值模拟方法.采用观测分析与城市陆面模式相结合的方法,求出了城市中绿地的绿洲效应系数、城市不透水下垫面水份有效率的时间变化公式以及人为潜热排放的日变化曲线和四季最大值,可直接应用于城市陆面模式,有效提高了模式的整体模拟性能.该方法对其他城市地区有一定的参考意义.  相似文献   

5.
沙漠-绿洲陆-气相互作用和绿洲效应的数值模拟   总被引:17,自引:8,他引:9       下载免费PDF全文
本文利用一个已发展的陆面过程参数化方案与大气边界层模式耦合,模拟了半干旱区绿洲戈壁非均匀下垫面的陆面过程及其与大气边界层的相互作用,给出了“绿洲效应”这一自然现象垂直剖面上更为清晰、准确和细致的结构特征.数值模拟的结果与早前的许多观测实验结论相吻合,即“绿洲效应”具有明显的“冷岛效应”和“湿岛效应”;它表现为在绿洲区域比戈壁沙漠区域环境温度低、湿度大、湍流动能输送弱,具有下沉气流而导致与周围戈壁沙漠区域产生水平输送环流.而更加细致地研究这些现象对于深入了解绿洲气候的形成和绿洲的维持机理具有重要的意义.  相似文献   

6.
何文英  陈洪滨  李军 《地球物理学报》1954,63(10):3573-3584
复杂多变的陆地表微波比辐射率,造成陆面上星载微波观测反演大气参数较为困难,也使得许多卫星微波资料不易同化应用到数值模式,因此迫切需要提供准确可靠的陆面微波地表比辐射率信息.随着卫星观测技术的迅速发展,利用丰富的星载被动微波观测直接反演陆面微波比辐射率成为一种主要手段.国外针对星载微波成像仪和微波垂直探测器开展较为系统的陆面微波比辐射率研究,建立不同类型的地表比辐射率反演方法,开发地表比辐射率参数化方法并应用于辐射资料同化.对于卫星观测反演陆面微波比辐射率存在的问题,开展了评估分析和方法订正.国内利用卫星观测也开展了一些陆面微波比辐射率研究工作,尚需要系统、综合的提炼.对于地表特征复杂的中国地区,还需要评估认识不同陆面微波比辐射率反演方法在我国适用情况,需要增强陆面微波比辐射率数据质量的认识以及业务应用.  相似文献   

7.
陆面过程模式中有关土壤水热传输、植被冠层和空气动力学等过程的关键参数的不确定性严重制约着地表-大气相互作用模拟能力的提高.本文利用架设在我国吉林通榆和甘肃榆中典型半干旱区陆-气相互作用野外试验站的观测资料,结合大气边界层理论,利用多种方法系统估算了上述试验站地表空气动力学粗糙度(z0m)以及热传输附加阻尼(?B?1)的量值.结果表明,z0m在半干旱区具有明显的季节变化和年际变化特征,且在植被低矮的下垫面与现行通用的陆面模式中的默认值相差较大;而?B?1的日变化和季节变化特征明显.将修正后的z0m和?B?1参数化方案引入陆面过程模式,发现能够明显改善模式对于半干旱区地表感热通量的模拟能力.这些结果说明有必要进一步结合半干旱区的野外观测试验对陆面模式在该地区的缺省参数设置进行更广泛的评估,而基于外场观测试验和大气边界层理论估算的关键地表参数对于改进陆气相互作用的模拟体现出较大的应用潜力.  相似文献   

8.
考虑次网格变异性和土壤冻融过程的土壤湿度同化方案   总被引:3,自引:0,他引:3  
集合Kalman滤波以其简单有效的特点在陆面数据同化中广泛应用,通常作为预报模型的陆面过程模式往往要考虑模式次网格变异性和土壤冻融过程,若对此不加考虑而直接对土壤湿度进行同化可能会使得同化结果发生偏差.将双集合Kalman滤波应用于土壤湿度的同化,基于NCAR/CLM陆面过程模式建立了一个考虑次网格变异性和土壤冻融过程的土壤湿度同化方案:在同一个时间步内用状态滤波对模式网格内某片上液态水分含量进行优化,用参数滤波对该片上的固态水分含量和其他片上的液态/固态水分含量进行优化,由此考虑模式次网格变异性和土壤冻融过程的影响,从而实现对整个模式网格上土壤湿度的同化.初步的同化试验表明:其同化效果在有、无土壤冻融阶段都优于一般的不考虑次网格变异性和土壤冻融变化的同化方案;该同化方案不仅能够提高那些有直接观测信息的土壤层的土壤湿度模拟精度,还能在一定程度上改善那些没有任何观测信息的土壤层的模拟效果;另外,土壤湿度同化结果的改善还能在一定程度上提高陆面模式对于土壤温度的模拟精度.  相似文献   

9.
利用中尺度模式WRFV2.2及NCEP再分析资料,以24h短期天气模拟的方式,模拟了我国2003年7月下旬江南及华南地区的高温天气,检验了模拟高温对不同陆面方案的敏感性.结果表明,模式总体上能模拟出≥35℃的气温分布形势,但模拟高温对不同陆面方案比较敏感,在模拟强度、TS评分及误差上存在较大差异,如从评分来看,SLAB方案模拟最好,RUC次之,NOAH再次之;不同陆面方案能模拟出方案之间的系统性差异,如三个方案中SLAB方案模拟平均绝对误差最小,模拟低层西太平洋副热带高压脊线最为偏东;不同的陆面方案主要通过模拟地表感热通量的变化影响地面气温变化,由不同方案引起的环流变化与地面气温变化之间的正反馈机制不重要,甚至在主要高温区存在"负"反馈.改进陆面方案对提高地面气温的数值预报能力具有重要意义.  相似文献   

10.
动力学粗糙度长度是大气数值模式和微气象计算的重要物理参数,其准确性会影响大气数值模式的性能和微气象计算水平.然而影响动力学粗糙度长度的因素很多,通常给出的动力学粗糙度长度参数往往只考虑单个因素的作用,其适应能力具有明显的局限性,往往给陆面动量通量(摩擦速度)的估算带来很大误差.本文利用近几年对动力学粗糙度长度与各种不同影响因素的参数关系的研究成果,构建了一个包含近地层大气动力特征、大气热力特征、植被自然生长规律、降水年际波动的生态效应和植被类型等多个影响因子的动力学粗糙度长度参数化方案,该参数化方案几乎考虑了平坦低矮植被动力学粗糙度长度的全部影响因素.同时,还利用兰州大学半干旱气候与环境观测站多年观测资料对该动力学粗糙度长度参数化方案的应用效果与其他实验方案进行了系统比较分析.发现该参数化方案计算的摩擦速度不仅远比用常数动力学粗糙度长度计算的摩擦速度误差小,而且也明显比仅考虑部分因素变化的其他参数化方案计算的摩擦速度误差小.与以往用固定不变的动力学粗糙度长度计算的摩擦速度相比,其与观测值的相关系数由0.752提高到了0.937,标准差和偏差分别降低了近20%和80%;其与观测值的平均值仅相差0.004 m s?1,相对误差仅1.6%左右,大大减少了近地层动量通量的估算误差.检验表明,本文构建的针对平坦低矮植被的动力学粗糙度长度的多因子普适性参数化方案可以为大气数值模式提供更加科学的参数化方案.  相似文献   

11.
基于卫星遥感资料的中国区域土壤湿度EnKF数据同化   总被引:6,自引:0,他引:6       下载免费PDF全文
土壤湿度在陆气相互作用过程中扮演着重要的角色,是气候、水文、农业、林业等研究中重要的地球物理参数之一.土壤湿度影响地面蒸散,径流、地表反射率、地表发射率以及地表感热和潜热通量,从而对气候有重要影响,它对大气的影响在全球尺度上仅次于海面温度,在陆地尺度其影响甚至超过海面温度.本文介绍了基于EnKF及陆面过程模型的中国区域陆面土壤湿度同化系统(CLSMDAS,China Land Soil Moisture Data Assimilation System),以及该系统应用于中国区域陆面土壤湿度同化试验的结果.CLSMDAS包括以下几个部分:1)陆面模式采用美国国家大气研究中心NCAR的陆面过程模型Community Land Model Version3.0(简写为CLM3.0);2)大气驱动场数据中的降水和地面入射太阳辐射数据来自FY2静止气象卫星每小时产品;3)陆面数据同化方法采用EnKF(Ensemble Kalman Filter)同化方法;4)观测数据包括AMSR-E卫星反演土壤湿度产品以及地面土壤湿度观测资料.利用CLSMDAS对2006年6~9月的土壤湿度同化试验结果的分析表明:陆面模式模拟和同化结果都能比较合理地反映出土壤湿度时空分布,同化的土壤湿度分布与2006年8月重庆、四川发生建国以来最严重的夏伏旱有非常好的对应关系,与发生在9月的湖北东部、广西南部等地的干旱区也有非常好的对应关系.  相似文献   

12.
沙漠陆面过程参数化与模拟   总被引:4,自引:2,他引:2       下载免费PDF全文
郑辉  刘树华 《地球物理学报》2013,56(7):2207-2217
沙漠地区植被稀疏、干旱少雨,其陆面物理过程具有与全球其它地区显著不同的特点.本文利用巴丹吉林沙漠观测资料,分析和计算了地表反照率、比辐射率、粗糙度和土壤热容量、热传导系数等关键陆面过程参数,建立了适合于沙漠地区的陆面过程模式DLSM (Desert Land Surface Model),并与NOAH陆面过程模式的模拟结果和观测资料进行了比较.结果表明:巴丹吉林沙漠地表反照率为0.273,比辐射率为0.950,地表粗糙度为1.55×10-3 m,土壤热容量和热扩散系数分别为1.08×106 J·m-3·K-1和3.34×10-7 m2·s.辐射传输、感热输送和土壤热传导过程是影响沙漠地区地表能量平衡的主要物理过程.通过对这三种过程的准确模拟检验,DLSM能够较准确地模拟巴丹吉林沙漠地气能量交换特征;短波辐射、长波辐射和感热通量的模拟结果与观测值间的标准差分别为7.98,6.14,33.9 W·m-2,与NOAH陆面过程模式的7.98,7.72,46.6 W·m-2的结果接近.地表反照率是沙漠地区最重要的陆面过程参数,地表反照率增大5%,向上短波辐射通量随之增加5%,感热通量则减小2.8%.本文研究结果对丰富陆面过程参数化方案,改进全球陆面过程模式、气候模式具有参考意义.  相似文献   

13.
基于微波亮温及集合Kalman滤波的土壤湿度同化方案   总被引:4,自引:0,他引:4       下载免费PDF全文
基于集合Kalman滤波及SCE-UA(shuffled complex evolution)算法发展了能够直接同化微波亮温的土壤湿度同化方案. 该方案以陆面过程模式CLM 3.0中的土壤水模型作为预报算子, 以辐射传输模型作为观测算子. 整个同化过程分为参数优化和土壤湿度同化两个阶段, 利用SCE-UA算法优化辐射传输模型中难以确定的植被光学厚度参数和地表粗糙度参数, 并利用优化参数作为观测算子的模型参数进行同化. 通过人工理想试验表明该同化方案可以明显改善表层土壤湿度的模拟精度, 并且对深层土壤湿度的模拟也有一定程度的改善; 利用AMSR-E亮温(10.65 GHz垂直极化)所进行的实际同化试验表明顶层(0~10 cm)土壤湿度同化结果与观测的均方根误差(RMSE)由模拟的0.05052减小到0.03355, 相对减小了33.6%, 而较深层(10~50 cm)平均减小了20.9%. 这些同化试验显示该同化方案的合理性.  相似文献   

14.
大气边界层研究进展   总被引:1,自引:0,他引:1  
大气边界层对云和对流的发展、演变有重要作用.本文回顾了在大气边界层高度计算方法,边界层的时空分布特征、结构和发展机理,以及边界层参数化方案等方面的主要研究进展.大气边界层高度计算方法主要分为基于大气廓线观测数据计算和基于模式参数化方案计算两大类;大气边界层高度频率分布形态具有明显的日变化特征,并且稳定、中性和对流边界层高度的频率分布呈现出不同的Gamma分布特征;地面湿度状况对边界层发展影响明显,对于不同的下垫面热力性质和地形状况,大气边界层高度呈现出明显的空间差异,青藏高原边界层高度明显高于一般平原地区;在强烈的地面加热驱动下,对流边界层与残余层通过正反馈机制循环增长可以形成4000 m以上的超高大气边界层;研制大气边界层、浅对流以及云物理方案的统一参数化框架是未来数值预报模式的发展趋势.  相似文献   

15.
土壤湿度不仅是地表水循环的重要组成部分,而且对天气和气候也有重要影响,它的模拟误差严重阻碍了人们对水循环的认知.本文首先评估了1°×1°水平分辨率的全球陆面数据同化产品(Global Land Data Assimilation System,GLDAS)对青藏高原中部那曲地区和东部玛曲地区土壤湿度的模拟性能;鉴于GLDAS较粗的分辨率无法精细描述分析区域土壤湿度空间分布特征,于是我们基于通用陆面过程模式(Community Land Surface Model,Version 4.5),开展了高分辨率0.1°×0.1°的模拟,并对高分辨率模拟土壤湿度误差的原因进行了深入分析.结果表明:(1)GLDAS陆面数据同化产品和高分辨率陆面模式模拟结果都可以反映出土壤湿度的季节变化特征,但在非冻结期均存在不同程度的干偏差,尤其是在玛曲地区;(2)对比观测和模拟的土壤湿度发现,观测数据表现出强烈的空间异质性,而模拟结果呈现出的是空间均一性.按照模拟误差进行归类分组,对比模拟性能优劣的两组站点发现,模式物理过程不是模拟性能差异的主要因素,而两组站点间地表特征参数中的土壤质地和地形参数,以及驱动数据均没有体现出空间异质性,这可能是土壤湿度模拟结果没有表现出空间异质性的原因.  相似文献   

16.
天气指某一个地区距离地表较近的大气层在短时间内的具体状态.大气中气象要素的空间分布可表现为各种瞬息万变的天气现象,这些天气的分布和变化是由不同时空尺度的天气系统引起的.天气与民生息息相关,其发展演变一直是大气科学研究和应用的重点领域.天气学的发展与观测系统、动力学理论和数值模式的发展密切相连.中国从20世纪50年代初开始建设观测网,到目前已建成门类齐全、布局合理的地基、空基和天基综合气象观测系统.特别是新一代稠密雷达网以及风云卫星系列的发展以及多次大型野外观测试验的实施使我们对天气的认识从宏观的天气形势深入到中小尺度天气系统精细热动力、云微物理结构和演变特征.观测系统的发展同时也促进了理论、数值模式和模拟的发展,中国已由初期主要以引进国外模式为主发展为目前主要发展具有中国自主知识产权的数值模式系统,基于高分辨数值模拟结果对不同尺度天气的发生发展机理和可预报性有了深入理解.此外,天气学已由初期的独立发展逐渐向多学科交叉方向转变,气候和环境的变化与天气演变之间的相互作用已成为大气科学的热点和前沿问题.文章重点回顾过去70年来中国在对天气演变起重要作用的天气现象及其短期变化过程的物理本质、演变规律和预报方法领域所取得的重大科学和技术成果,主要根据正式发表的文献从大气动力学、天气尺度天气特征、台风及热带天气、强对流天气特征、数值天气预报及资料同化,以及天气与气候、大气物理及大气环境等交叉领域六个方面分别加以综述.  相似文献   

17.
使用地球观测系统的中分辨率成像光谱仪(EOS.MODIS)提供的归一化植被指数(NDVI)产品估算植被覆盖度和航天飞机雷达地形测绘任务(SRTM)制成的数字高程模型(DEM)数据遥感产品替换WRF模式默认的植被覆盖度和地形高度,并且利用WRF模式及其先进的三维变分同化系统(WRF.3DVar)循环同化东北半干旱区自动气象站近地面气象要素,对东北半干旱区的温度场、湿度场、风场和能量场的结构及其日变化特征进行了较为细致的模拟研究.通过4组数值模拟试验分别探讨了同化气象要素与改变模式地表参数引起的不同下垫面潜热、感热的分配关系和降水、土壤湿度变化弓f起的地表能量通量模拟效果,并利用通榆站、奈曼站、锦州站、和密云站2009年6-8月的通量观测资料与模拟结果对比检验.结果表明,WRF模式能够较好地模拟出东北半干旱区夏季的近地面温度、风向、净辐射、感热和潜热等要素的变化特征及日变化规律.同化试验(Case2)模拟的近地面气温、相对湿度、风速相比控制性试验(Case1)有所改善;陆面参数试验(Case3)和集合试验(Case4)改善了感热和地表热通量的模拟.WRF模式能较好地模拟出下垫面土壤湿度随时间变化的规律,集合试验(Case4)土壤湿度模拟结果与4个通量站观测值相比无太大差别,但降水的模拟有待改善.本研究利用卫星遥感资料改善模式下垫面陆面参数,利用气象资料同化改善近地面大气要素模拟精度,这是将各种不同空间和时间尺度的多源数据与数值模拟融合的有益尝试.此研究生成的东北地区资料同化数据集可用于气候变化、干旱监测等方面,对深入了解半干旱区气候的形成和维持机理具有重要的意义.  相似文献   

18.
积云对流参数化方案对热带大气季节内振荡模拟的影响   总被引:2,自引:0,他引:2  
本文利用中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室(LASG)的大气环流谱模式的最新版本SAMIL2.2.3系统分析了积云对流参数化方案对热带大气季节内振荡(ISO)模拟的影响.结果表明该模式的三种积云对流参数化方案都能模拟出热带地区气候平均季节内振荡(CISO)时空演变的基本特征;其中Tiedtke方案相比于Manabe方案和Zhang-McFarlane方案,模拟出的CISO与观测更为接近.而SAMIL2.2.3模式对ISO年际变化部分(TISO)的模拟结果显示尽管SAMIL2.2.3相对于之前的版本对Tiedtke方案和Zhang-McFarlane方案针对不同问题做了必要调整,但三种积云参数化方案中只有Zhang-McFarlane方案能够模拟出27~50天TISO模态的东传特征.另外,SAMIL2.2.3与其较早版本SAMIL2.08中的Manabe方案完全相同,物理过程的其他部分如辐射方案和气溶胶条件有所调整,但仅后者能模拟出100E以西热带27~50天TISO模态的东传特征.因此,模式发展过程中物理参数化方案的改进并非在每一方面上都产生正面效果.动力框架和物理过程,以及不同物理过程之间的协调发展是改进模式的重要途径.  相似文献   

19.
本文针对模式发展的需要, 在Farouki土壤热传导率参数化方案的基础上, 综合Johansen和Côté的参数化方案, 发展了一个用于青藏高原中部的土壤热传导率参数化方案, 用“全球协调加强观测计划之亚澳季风青藏高原试验计划(CAMP/Tibet)”中那曲布交(BJ)站实际资料对该方案进行了检验, 并将它用于公共陆面模式(CoLM)中, 对青藏高原那曲地区进行了单点数值模拟试验. 结果表明: 在未冻结及冻结土壤中, 新方案比Farouki方案计算的土壤热传导率小, 更接近实测值. 加入新方案的CoLM模式对土壤温度模拟的准确性比原模式有一定的提高.  相似文献   

20.
中尺度涡在大洋环流中有着重要的作用. 为检验高度计资料同化应用于中尺度涡模拟的效果, 提高环流模式对中尺度涡模拟的精度, 使用三维变分的OVALS同化系统结合POM模式进行了中尺度涡同化模拟实验. 该系统将卫星高度计资料同化反演成为温盐伪观测数据, 然后再次进行常规温盐同化, 得到温盐分析场. 根据中尺度涡现象的物理特性, 经过对比检验, 对系统中的同化参数作了合理设置. 使用T/P和ERS1/2的10 a卫星高度计的海表高度异常资料, 分别用同化与非同化两种方案对西北太平洋的中尺度涡进行对比模拟实验. 将实验结果与观测数据的比较表明, 加入高度计资料同化的模拟结果远远好于未使用同化的模拟结果. 说明高度计反演温盐场的同化方法用于对于模拟中尺度涡现象是非常有效的.  相似文献   

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

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