首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 177 毫秒
1.
崔云扬  周毓荃  蔡淼 《大气科学》2019,43(3):618-633
本文利用河北邢台测站Ka波段微雨雷达(MRR)观测到的一次冷锋云系降水过程分析降水的垂直分布及演变特征。将MRR观测结果与天气雷达、地面雨滴谱仪、雨量计观测结果进行对比以检验MRR数据的可靠性。同时将MRR与雨滴谱仪和激光云高仪结合,研究了不同相对湿度阶段特征量、雨滴谱的平均垂直分布特征和降水特征量随时间、高度的演变特征。结果表明:MRR与雨量计及雨滴谱仪累计雨量结果较为接近,趋势一致。MRR 200 m雨强值与地面雨滴谱仪雨强值偏差最小,平均偏差为0.05 mm h?1,相关系数为0.93。相比雨滴谱仪,MRR观测到的小滴数浓度出现高估,大滴数浓度出现低估,中滴数浓度较为一致。降水在云内和云外受不同微物理过程影响,垂直变化特征不同。降水初期平均反射率和雨强在云底以下明显减小,小滴和中滴平均数浓度明显减小,蒸发作用影响较强。而在其余时间段在云内随高度降低平均反射率和雨强略有增加,小滴平均数浓度变化较小,中滴大滴平均数浓度增加,表明云内有云滴与雨滴间的碰并发生。而在云外低层,随高度降低平均有效直径明显增加,平均雨滴总数浓度明显减小,小滴平均数浓度显著减小,大滴平均数浓度显著增加,表明在云外低层雨滴间的碰并作用较强。  相似文献   

2.
黄山层状云和对流云降水不同高度的雨滴谱统计特征分析   总被引:2,自引:1,他引:1  
李慧  银燕  单云鹏  金祺 《大气科学》2018,42(2):268-280
根据2011年6~7月在黄山不同高度采用PARSIVEL雨滴谱仪测得的雨滴谱数据,对不同海拔高度上两类(层状云和对流云)降水粒子谱的微物理特征量、Gamma函数拟合以及雨滴的下落速度进行对比分析,结果表明:对流云降水的雨水含量和降水强度、雨滴的各类尺度参数和数浓度都比相同位置上层状云降水的大,同类降水中,山腰的雨滴尺度大于山顶和山底,这可能与各观测点和云底相对位置的不同有关;随降水强度增加,雨滴的质量加权平均直径Dm逐渐增大,广义截距参数(log10Nw)的标准差逐渐减小。拟合结果表明各高度的雨滴谱都比较符合Gamma分布,由拟合参数分析雨滴谱的演变,发现相对于对流云降水,层状云降水粒子谱随高度的变化较小,雨滴谱的演变较为稳定。此外,本文还对两类降水中雨滴的下落速度及影响落速的因素进行了分析。  相似文献   

3.
大气垂直运动对雷达估测降水的影响   总被引:1,自引:0,他引:1       下载免费PDF全文
运用雷达反射率因子(Z)与降水强度(R)之间的关系定量估测降水,降水云体中的大气垂直运动(wa)不可忽视。PARSIVEL激光雨滴谱仪(简称PARSIVEL)在获取雨滴粒径分布的同时可以从测量的雨滴下降速度分布中提取wa,用于分析PARSIVEL高度上的大气垂直运动对雷达Z估测降水强度影响。使用2014年5-6月华南季风降水观测试验期间广东阳江5次层状云、6次对流性降水过程中PARSIVEL数据分析大气垂直运动对定量估测降水影响,wa对层状云和对流云降水强度影响范围分别为-0.18~1.05 mm·h-1和-5.44~24.81 mm·h-1,相对影响值分别为-13.61%~13.99%,-38.59%~25.92%。静止大气条件下,雷达估算降水Z-R关系式中系数A,b引起的层状云和对流云降水估测偏差平均分别为10.9%和25.5%。真实大气中雷达估测降水的偏差平均情况是层状云降水由于wa的对消作用降低为9.2%,对流云降水则增加到51.2%。对流性降水中大气垂直运动对雷达估测降水的影响较大。  相似文献   

4.
沙修竹  丁建芳  程博 《气象》2019,45(11):1569-1578
采用河南省2016—2017年100个日降水资料,对比分析雨滴谱反演回波与雷达回波的差异、雨滴谱反演降水强度与雨量计观测降水强度的差异;进行雨滴谱Gamma拟合,以探究河南省雨滴谱分布及降水云系类型;进行Z-I关系拟合,以探究河南省降水回波与降水强度的关系。结果表明:(1)雨滴谱反演回波、雷达观测回波的变化趋势具有较好一致性。而前者普遍小于后者,其可能原因:一是雷达通过最低仰角观测到的地面雨滴谱仪上方回波与地面雨滴谱仪之间存在一定高度差,二是雨滴下落时的蒸发、破碎过程,使到达地面的雨滴直径减小。(2)雨滴谱反演的降水强度与雨量计观测的降水强度相比,存在一定差异,但无显著偏大或偏小规律性特征。(3)对流云及层积混合云的雨滴谱宽大于层状云,中等尺度雨滴数密度较大。层状云的小水滴数密度较大。河南省大部分降水过程为雨滴谱较窄的层状云降水。(4)河南省降水回波与降水强度的拟合公式:Z=262I~(1.34),层状云拟合公式:Z=219I~(1.30),对流云拟合公式:Z=307I~(1.38)。(5)雨滴数浓度较高月份为6—7月(1500个·m~(-3)左右),降水强度较高月份为8—10月(60 mm·h~(-1)),雨滴最大直径较高月份为4—8月(4.3~4.8 mm),雨滴平均直径较高月份为3—4月(3 mm左右)。雨滴数浓度、降水强度、最大直径、平均直径的月份特征变化无一致性。  相似文献   

5.
激光雨滴谱仪测速误差对雨滴谱分布的影响   总被引:2,自引:2,他引:0       下载免费PDF全文
该文研究PARSIVEL激光雨滴谱仪的测量误差并提出订正方法。对2014年在广东阳江的PARSIVEL激光雨滴谱仪采集的两次降水过程数据进行分析发现,雨滴下落速度V随粒径D变化与静止大气中雨滴下落末速度随粒径变化的Atlas-Ulbrich曲线分布趋势一致,但D < 1 mm及D > 3 mm的速度偏差较大。其主要原因是大粒径雨滴形变造成速度偏离较大,仪器测量误差造成小粒径测速偏大,激光雨滴谱仪所在高度的大气垂直运动影响雨滴下落速度。根据PARSIVEL激光雨滴谱仪测量原理,基于雨滴形变与粒径关系,给出形变订正后的Atlas-Ulbrich修正曲线,并用于对小粒径测速订正。比较订正前后的雨滴谱分布,订正后的小雨滴浓度明显增加,大雨滴浓度略有减小,订正后浓度参数和斜率参数均增加,形状参数变化不明显。  相似文献   

6.
利用宁波微雨雷达(MRR)数据,进行仪器观测雨强评估,在此基础上,分析了宁波地区三次典型过程的雨滴谱和降水特征量分布情况,并对降水特征量在雨强预报中的应用做了初步分析。结果表明,MRR具备较可靠的降水观测能力,三次过程中雨滴粒子浓度均随雨滴直径增大而减小。西风槽和梅雨过程雨滴谱基本相似,液态水含量均较小且变化不明显;而台风"黑格比"雨滴粒子浓度更高,粒径范围更大,液态水含量显著提高且随高度降低而减小,粒子有效直径随高度降低明显增大,粒子数浓度明显减小,说明台风中雨滴碰并增长作用较强,大量小粒子相互碰并成大粒子。整层积分的液态水含量和整层平均的粒子有效直径对降水强度预报有较好的指示作用,二者变化均提前于降水强度变化约10min左右。  相似文献   

7.
基于伊犁河谷伊宁站内布设的微雨雷达(MRR)、地面OTT-PARSIVEL雨滴谱仪以及雨量筒观测资料,针对2019年9月30日伊宁地区一次大雨过程,对比检验MRR数据的可靠性,在此基础上探究降水不同时期MRR的微物理量垂直分布特征。结果表明:3种仪器观测的降水量较为接近、变化趋势一致,MRR观测的近地面35、70和105 m高度层雨强与OTT-PARSIVEL雨滴谱仪的观测值有较好的相关性,决定系数分别为0.9233、0.9289和0.9186,且当雨强低于3 mm·h~(-1)时,MRR与雨滴谱仪的雨强拟合收敛程度更高。不同雨强阶段,雨滴谱的垂直分布存在差异。低雨强阶段,降水初期环境湿度较低、蒸发较大,MRR的雷达反射率因子Z、液态水含量LWC及雨强R_i随高度降低而减小;降水中期,环境湿度较大、蒸发较小,Z、LWC、R_i随高度变化不大;降水末期,空中水汽和动力条件供应不足,靠近地面R_i明显减弱。中、高雨强阶段,粒子下落速度W较为稳定,而Z、LWC、R_i受雨滴间碰并作用影响较大,整体随高度降低而增大。此次大雨过程主要以小雨滴为主,其平均数浓度占总数浓度90%以上,且随高度降低而减少;中雨滴对雨强的贡献最大,不同雨强阶段的贡献率均在60%以上,且随高度降低而增大;大雨滴数浓度占比及对雨强的贡献均最小。  相似文献   

8.
两次降水过程的微降雨雷达探测精度分析   总被引:3,自引:3,他引:0  
温龙  刘溯  赵坤  李杨  李力 《气象》2015,41(5):577-587
垂直指向微降雨雷达(MRR)能够测量从近地面至高空的雷达反射率因子和雨滴谱分布特征,对认识降水微物理结构,改进雷达定量降水估计精度有重要作用。为评估MRR探测的雨滴谱分布、降水和雷达回波精度,利用南京地区夏季观测的两次降水过程,将MRR与业务S波段天气雷达、二维视频雨滴谱仪、常规雨量筒观测进行层状云降水和对流性降水下的定量对比分析。结果表明,MRR垂直探测的雷达反射率因子与S波段雷达观测在中低层(<4 km)平均差异<1 dB, 但高层(>4 km)出现显著低估,且该现象随降水强度增强更明显,这主要是雷达回波衰减导致。MRR在回波强度<35 dBz时对降水率的探测精度较高,但在>35 dBz时低估降水。其中,层状云降水的降水率比对流性降水更接近雨量筒观测。常规雨量筒对0.1 mm以下的降水无探测能力,而MRR探测敏感度较高,对于微弱降水率的估计效果也很好。由于MRR最大探测范围的限制,相对于2DVD而言,MRR探测的最大粒子直径低估、最小粒子浓度高估,但在中间段的探测效果和2DVD雨滴谱观测一致性较高。总体而言,MRR是一个有效的降水探测仪器,其探测结果在层状云降水过程中优于对流性降水过程。  相似文献   

9.
北京冬季降水粒子谱及其下落速度的分布特征   总被引:2,自引:0,他引:2  
为了深入探讨北京冬季云降水的微物理特征,提高雷达反演冬季固态降水的精度和冬季降水的预报水平,利用PARSIVEL(Particle Size and Velocity)降水粒子谱仪所观测的冬季降水粒子谱,结合地面显微镜粒子图像和云雷达数据,对比分析了北京海坨山地区冬季过冷雨滴、霰粒、雪花、混合态降水的粒子谱和下落速度特征,得到主要结论如下:(1)霰粒降水过程的云顶最高,整层的含水量最大,低层的退偏振比(LDR)最小,粒子更接近于球形;降雪过程的云顶最低,云中含水量最少,低层的退偏振比较大;混合态降水过程的雷达回波强度和高度特征介于两者之间,但低层的退偏振比最大;(2)在云中上升或下沉气流及湍流的影响下,过冷雨滴、霰粒和雪的下落速度均对称分布于各自理论下落末速度曲线的两侧。因此可根据粒子浓度相对于其直径和速度分布的中轴线位置,判断出该段降水过程中的主要粒子形态;(3)冬季雪花、霰粒和混合态降水粒子下落速度分布的散度较雨滴更大,其原因是由于冷云降水过程的粒子形态复杂,且固态粒子下落过程中更容易受破碎、聚并和凇附等微物理过程影响;(4)在4种降水类型中,雪的平均直径和离散度最大,雨滴最小;混合态降水粒子的总数浓度最大,雨滴的总数浓度最低,并且4种降水类型的粒子数浓度、平均直径和离散度均随降水强度的增大而增大。   相似文献   

10.
刘彦  苏德斌  杨宁  王亮 《高原气象》2023,(3):748-757
布设在内蒙古自治区巴彦淖尔市五原县气象局的二维雨滴谱仪(2DVD)在2020年7月27日观测到一次含雹混合降水事件,基于粒子直径与本地下落末速度对2DVD数据质控后,分析了此次过程不同直径-速度和不同直径-轴比下的粒子数量分布以及粒子谱、中值体积直径、质量加权平均直径、粒子数浓度和降水强度等参数随时间的演变。结果表明:(1)直径0.1~0.5 mm的雨滴实测下落末速度偏大于经验公式计算的下落末速度,直径5~10.4 mm的冰雹粒子下落速度为5.5~11.6 m·s-1,直径小于6 mm的冰雹粒子下落速度分布较广。直径0.1~2 mm的小雨滴轴比为0.9~1.1,直径2~5 mm的大雨滴轴比为0.7~1.0,直径5~10.4 mm的冰雹粒子轴比为0.5~1。(2)本次降雹类型为先雨后雹,雹雨混降,雹后持续降雨,冰雹谱谱宽为10.4 mm。逐分钟降水粒子谱存在3个直径极大值,在降水和降雹之间存在短暂的无降水,谱宽和粒子数浓度随时间同时增大和减小,并在降雹阶段,各参数陡增至峰值。(3)冰雹谱分布呈单调递减型,通过M-P分布函数分段模拟了降雹时段的雨滴谱和冰雹谱,模拟的粒...  相似文献   

11.
研究不同云系降水的微物理参数特征,对研究降水机制、人工影响天气、雷达定量测量降水、数值预报模式中微物理参数化方案的选择等都有一定意义。本文针对2015年济南地区的液态降水过程,基于微降水雷达(Micro Rain Radar,简称MRR)资料,研究不同云系降水的微物理参数。在400 m高度上,层状云降水0.02~0.2 mm h-1雨强样本数很大,但对累计降水量的贡献很小。混合云和对流云降水在大粒子端数浓度较高。在垂直方向上,层状云降水中的粒子的尺度较集中,中值体积直径D0平均在1 mm左右,随高度的变化不大。对流云降水在雨强大于20 mm h-1时,强垂直气流(包括上升气流和下沉气流)对粒子直径的影响较大,进而影响空中微降水雷达反演降水参数的数据质量。而垂直气流的影响对层状云降水影响较小,在层状云降水时,微降水雷达可以用来分析零度层亮带以下雨滴谱在垂直方向上的演变。  相似文献   

12.
The polarimetric radar network in Jiangsu Province has just been operationalized since 2020. The first intense precipitation event observed by this polarimetric radar network and disdrometer occurred during August 28-29, 2020 and caused severe flooding and serious damage in eastern Jiangsu Province. The microphysics and kinetics for this heavy precipitation convective storm is diagnosed in this study, in order to promote the application of this polarimetric radar network. Drop size distribution (DSD) of this event is estimated from measurements of a ground disdrometer, and the corresponding three-dimensional atmospheric microphysical features are obtained from the multiple polarimetric radars. According to features of updraft and lighting, the evolution of the convective storm is divided into four stages: developing, mature with lightning, mature without lightning and dissipating. The DSD of this event is featured by a large number of raindrops and a considerable number of large raindrops. The microphysical characteristics are similar to those of warm-rain process, and ice-phase microphysical processes are active in the mature stages. The composite vertical structure of the convective storm indicates that deep ZDR and KDP columns coincide with strong updrafts during both mature stages. The hierarchical microphysical structure retrieved by the Hydrometeor Identification Algorithm (HID) shows that depositional growth has occurred above the melting level, and aggregation is the most widespread ice-phase process at the -10℃ level or higher. During negative lightning activity, the presence of strongest updrafts and a large amount of ice-phase graupel by riming between the 0℃ and -35℃ layers generate strong negative electric fields within the cloud. These convective storms are typical warm clouds with very high precipitation efficiency, which cause high concentration of raindrops, especially the presence of large raindrops within a short period of time. The ice-phase microphysical processes above the melting layer also play an important role in the triggering and enhancing of precipitation.  相似文献   

13.
刘黎平  谢蕾  崔哲虎 《大气科学》2014,38(2):223-236
本文首先利用数值模拟的方法,分析了利用毫米波云雷达功率谱密度反演雨滴谱时,降水粒子米散射效应、空气湍流、空气上升速度等对雨滴谱和液态水含量等参数反演的影响;建立了功率谱密度处理及其直接反演雨滴谱、液态水含量、降水强度和空气上升速度的方法;并利用2012年7月在云南腾冲观测的二次弱降水数据,采用毫米波雷达和Ku波段微降水雷达观测的回波强度、径向速度垂直廓线以及780 m高度上的功率谱密度对比的方法,以及毫米波云雷达观测的780 m高度上功率谱密度、回波强度与地面雨滴谱计算得到的这些量的对比方法,分析了毫米波雷达数据的可靠性;并将780 m高度上毫米波雷达反演的雨滴谱与地面雨滴谱数据进行了对比,分析了毫米波雷达反演的雨滴谱的准确性;分析了毫米波雷达回波强度偏弱的原因,讨论了该高度以下降水对毫米波雷达衰减的影响。结果表明:空气湍流对弱降水微物理参数反演影响不大,而空气上升速度和米散射效应均对反演结果有一定影响;毫米波雷达观测到的径向速度和功率谱密度与微降水雷达比较一致,回波强度的垂直廓线的形状与微降水雷达也比较一致,但毫米波雷达观测的回波强度偏弱;与雨滴谱计算值相比,毫米波雷达观测的低层的回波强度也偏弱,天线上的积水是造成毫米波雷达回波强度变弱的主要原因。毫米波雷达观测的低层的功率谱密度与地面雨滴谱观测的数据形状比较一致,但有一定的位移。毫米波雷达反演的雨滴谱与地面观测的谱型和粒子大小也比较一致。这些结果初步验证了毫米波雷达观测的功率谱密度及其反演方法的可靠性。  相似文献   

14.
为了更好地认识雹暴等强对流降水的微物理特征,利用降水现象仪观测资料和CINRAD/SA-D双偏振天气雷达的粒子相态识别和反射率因子等产品,分析2019年8月16日出现在山东北部的一次雹暴降水的雨、冰粒子谱(包括直径2~5 mm的霰粒子和直径大于5 mm的冰雹)的识别以及雨滴谱的演变特征,结果显示:济南双偏振天气雷达0.5°仰角PPI上在德州和陵县观测点附近识别的粒子主要为雨滴,临邑为冰雹加大雨;安装在3个观测点的降水现象仪均识别出少量冰粒子,3个观测点的冰粒子谱数密度低、分布不连续。雹暴降水开始阶段,出现少量大雨滴,这是风的筛选和蒸发作用导致的雨滴谱,具有低的小粒子数密度和较多大雨滴;雹暴强降水增大阶段,雨滴谱特征是小雨滴数密度偏低、大雨滴较多,即总的雨滴浓度低、雷达反射率因子高;雨强减弱阶段,小雨滴数密度显著增大、大雨滴数密度偏少,即总的雨滴浓度显著增大,但雷达反射率因子偏低。  相似文献   

15.
利用GPM卫星探测两个时次的资料,以1808号超强台风"玛利亚"为研究对象,分析了台风降水率、降水类型及台风高度水平分布,降水率垂直廓线变化特征,以及降水率三维结构分布特征。得出以下主要结论:两个时刻"玛利亚"均处在超强台风级,A时刻台风眼区为深厚对流区,B时刻眼区对流有所减弱,但是有强螺旋雨带出现。A、B时刻的降水率最大值与风暴顶高度并非一一对应,还与降水云系中微物理过程有关。GMI低频18.7 GHz探测的水粒子含量的大值区与强降水率对应较好,高频183.31±3 GHz探测的冰粒子信号与风暴顶高度分布一致。不同降水率对应的垂直廓线表明,降水率在5 km高度出现急剧变化,这是由于在该高度上雨滴碰并增长或者蒸发减小。从A时刻到B时刻,云墙区大于10 mm·h~(-1)的云墙半径内缩,B时刻眼壁与螺旋雨带之间存在着弱降水区及无降水区。  相似文献   

16.
利用设在伊宁的激光雨滴谱仪获取的2013年4月的降水资料,对层状云和混合云降水粒子谱的微物理参量平均值和Gamma函数拟合结果以及Z-I关系进行对比分析。计算结果表明,伊宁地区春季降水的微物理参量普遍偏小,小滴对降水浓度的贡献达到92%以上,即降水主要以小滴为主。层状云降水的雨强、雨滴数浓度、雨滴的各类微物理特征参量的平均值均大于混合云降水。函数拟合结果表明,混合云降水的雨滴谱宽大于层状云降水的雨滴谱宽,层状云和混合云降水的雨滴谱都比较符合Gamma分布,在小滴段Gamma分布对实际谱都有一定的低估,在大于1 mm的粒径段,拟合结果有一定的偏差。还讨论了雨滴大小因子Λ和形状因子μ之间的关系以及Z-I关系,Λ-μ关系与粒子尺度有关,根据拟合的二项式得到层状云降水粒子的平均直径大于混合云降水的平均直径。  相似文献   

17.
To analyze the effects of gas cannons on clouds and precipitation, multisource observational data, including those from National Centers for Environmental Prediction (NCEP) reanalysis, Hangzhou and Huzhou new-generation weather radars, laser disdrometer, ground-based automatic weather station, wind profiler radar, and Lin''an C-band dual polarization radar, were adopted in this study. Based on the variational dual-Doppler wind retrieval method and the polarimetric variables obtained by the dual-polarization radar, we analyzed the microphysical processes and the variations in the macro - and microphysical quantities in clouds from the perspective of the synoptic background before precipitation enhancement, the polarization echo characteristics before, during and after enhancement, and the evolution of the fine three-dimensional kinematic structure and the microphysical structure. The results show that the precipitation enhancement operation promoted the development of radar echoes and prolonged their duration, and both the horizontal and vertical wind speeds increased. The dual-polarization radar echo showed that the diameter of the precipitation particles increased, and the concentration of raindrops increased after precipitation enhancement. The raindrops were lifted to a height corresponding to 0 to -20 ℃ due to vertical updrafts. Based on the disdrometer data during precipitation enhancement, the concentration of small raindrops (lgNw) showed a significant increase, and the mass-weighted diameter Dm value decreased, indicating that the precipitation enhancement operation played a certain“lubricating”effect. After the precipitation enhancement, the concentration of raindrops did not change much compared with that during the enhancement process, while the Dm increased, corresponding to an increase in rain intensity. The results suggest the positive effect of gas cannons on precipitation enhancement.  相似文献   

18.
微降水雷达测量精度分析   总被引:2,自引:0,他引:2  
利用数值模拟的方法,讨论了利用微降水雷达MRR(Micro Rain Radar)雷达功率谱密度反演降水参数时,MIE散射(米散射)效应、垂直气流(包括上升气流、下沉气流)对数浓度N、雷达反射率Z、雨强I、液态含水量LWC等参数的影响。MIE散射主要影响直径为1.20~4.00 mm的粒子,MIE散射效应影响的N、Z、I、LWC偏差的平均值分别为2.74 m-3 mm-1、1.47 d BZ、0.0061 mm h~(-1)、0.0004 g m-3。下沉气流使反演液滴直径偏大,上升气流使得反演的液滴直径偏小,下沉气流的影响更大,尤其是对低层影响大于高层。例如,在300 m高度上,当液滴直径为2.67 mm时,下沉气流为2.00 m s-1时,理论上反演的直径为8.07 mm,超出了MRR探测的阈值,其相对误差值能接近200%。下沉气流使得反射率谱向大粒子方向平移,且谱型展宽;上升气流则相反。将MRR资料与同步观测的THIES雨滴谱仪数据进行比对,分析MRR资料的可靠性。选取2015年4月1日01~12时(协调世界时)山东济南的一次降水过程,将MRR在300 m高度上反演的雷达反射率因子、雨强、数浓度、中值体积直径与雨滴谱仪资料进行对比。结果表明:两种仪器探测的Z、I、N、中值体积直径D0在时间序列上都有较好的吻合度,变化趋势和幅度相近,Z、I、D0的平均偏差分别为1.19 d BZ、0.34 mm h~(-1)、0.36 mm。MRR反演的I值偏大,而粒子直径偏小,分析了产生偏差的主要原因,除了探测系统偏差、分析方法本身存在的偏差外,上升气流导致的偏差不容忽视。这些结果初步验证了微降水雷达观测的功率谱密度及其反演方法的可靠性。  相似文献   

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

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