基于二维雨滴谱仪的巴彦淖尔地区降雹谱个例分析_刘彦.pdf
《基于二维雨滴谱仪的巴彦淖尔地区降雹谱个例分析_刘彦.pdf》由会员分享,可在线阅读,更多相关《基于二维雨滴谱仪的巴彦淖尔地区降雹谱个例分析_刘彦.pdf(10页珍藏版)》请在咨信网上搜索。
1、第 42 卷 第 3 期2023 年 6 月高原气象PLATEAU METEOROLOGYVol.42 No.3June,2023刘彦,苏德斌,杨宁,等,2023.基于二维雨滴谱仪的巴彦淖尔地区降雹谱个例分析 J.高原气象,42(3):748-757.LIU Yan,SU Debin,Yang Ning,et al,2023.Case Study of Hail Size Distribution in Bayannur Area based on Two-Dimensional Video Disdrometer J.Plateau Meteorology,42(3):748-757.DOI
2、:10.7522/j.issn.1000-0534.2022.00059.基于二维雨滴谱仪的巴彦淖尔地区降雹谱个例分析刘彦1,2,3,苏德斌1,2,3,杨宁1,3,王亮1,3(1.成都信息工程大学电子工程学院,四川 成都 610225;2.中国气象局云雾物理环境重点开放实验室,北京 100000;3.中国气象局大气探测重点开放实验室,四川 成都 610225)摘要:布设在内蒙古自治区巴彦淖尔市五原县气象局的二维雨滴谱仪(2DVD)在2020年7月27日观测到一次含雹混合降水事件,基于粒子直径与本地下落末速度对2DVD数据质控后,分析了此次过程不同直径-速度和不同直径-轴比下的粒子数量分布以及粒
3、子谱、中值体积直径、质量加权平均直径、粒子数浓度和降水强度等参数随时间的演变。结果表明:(1)直径0.10.5 mm的雨滴实测下落末速度偏大于经验公式计算的下落末速度,直径510.4 mm的冰雹粒子下落速度为5.511.6 m s-1,直径小于6 mm的冰雹粒子下落速度分布较广。直径0.12 mm的小雨滴轴比为0.91.1,直径25 mm的大雨滴轴比为0.71.0,直径510.4 mm的冰雹粒子轴比为0.51。(2)本次降雹类型为先雨后雹,雹雨混降,雹后持续降雨,冰雹谱谱宽为10.4 mm。逐分钟降水粒子谱存在3个直径极大值,在降水和降雹之间存在短暂的无降水,谱宽和粒子数浓度随时间同时增大和减
4、小,并在降雹阶段,各参数陡增至峰值。(3)冰雹谱分布呈单调递减型,通过M-P分布函数分段模拟了降雹时段的雨滴谱和冰雹谱,模拟的粒子谱大于2DVD实测的粒子谱,且模拟雨滴谱与实测谱的数浓度对数均方根误差为1.02,模拟冰雹谱与实测谱的数浓度均方根误差为2.21。关键词:二维雨滴谱仪;雨滴谱;冰雹谱;M-P分布文章编号:1000-0534(2023)03-0748-10 中图分类号:P426.64 文献标识码:ADOI:10.7522/j.issn.1000-0534.2022.000591 引言 强对流天气是我国主要的灾害性天气之一,冰雹灾害则是强对流天气下发生比较常见的气象灾害,它的发生会给农
5、业生产和人民的生命财产安全构成威胁(李红斌等,2010;万红莲等,2017;刘新伟等,2021;吴海英等,2021)。巴彦淖尔市北部为乌拉特草原,中部为阴山山地,南部为河套平原,阴山山脉狼山段位于巴彦淖尔市中部,最高海拔约2800 m,将巴彦淖尔市一分为二,山前为富饶的河套平原,山后为荒漠化草原。受此特殊地理环境的影响,冰雹成为巴彦淖尔市夏季主要气象灾害之一。对雨滴谱的观测方法有滤纸法,面粉法和使用一维雨滴谱仪与二维雨滴谱仪器的光学法,Kathiravelu et al(2016)对雨滴观测技术和原理做了总结和分类。2DVD 可以快速记录每个粒子的二维信息,国内外许多学者使用2DVD研究降水的
6、微物理过程(Luo et al,2021;王梦瑶,2019)。观测冰雹的仪器有测雹板、雹雨分测计、一维雨滴谱仪和二维雨滴谱仪等。对冰雹的研究,国外 Macklin and Ludlam(1961)直接测量了从2 km高空落下的人造收稿日期:20211009;定稿日期:20220531资助项目:西北区域人工影响天气能力建设项目(ZQC-R18217);国家自然科学基金项目(42075001)作者简介:刘彦(1997-),男,四川宜宾人,硕士研究生,主要从事气象探测领域研究.E-mail:通信作者:苏德斌(1966-),男,江苏扬州人,教授,主要从事气象探测数据分析及应用、云和降水物理等领域研究E
7、-mail:刘彦等:基于二维雨滴谱仪的巴彦淖尔地区降雹谱个例分析3 期冰雹的最终下落速度。Auer(1972)在北美高原2.1 km的海拔测量了球形和锥形软冰雹及冰雹的下落末速度。Thurai et al(2011)指出一维雨滴谱仪和2DVD对小冰雹进行观测时,一维雨滴谱仪观测的质量加权平均直径和降雨率值高于 2DVD。Bringi et al(2017)结合2DVD和S波段双偏振雷达对一次冬季霰雨进行研究发现霰颗粒主要呈块状,平均形状较长。Dieling et al(2020)在实验室内测试探究冰雹密度、拖拽系数、冰雹形状、空气密度和大风对冰雹下落速度的影响。国内牛生杰等(1999)在固原使
8、用测雹板组成测雹板网对冰雹谱分布和冰雹落地动能进行观测分析。石安英等(1989)在张家口使用测雹板观测到315份雹谱资料,拟合出6类雹谱分布类型。孙玉稳等(2012)使用雹雨分测计在张家口地区记录的资料总结出了 5种降雹类型及相关因子分布特征。岳治国和梁谷(2018)使用一维雨滴谱仪对陕西渭北一次降雹过程进行了粒子谱特征分析,拟合了冰雹直径与平均冰雹末速度的关系。陶然亭(2020)结合2DVD和称重式降水传感器,分析了南京冬季降雪的雪花下落速度、密度和粒子谱分布等特征。王俊等(2021)基于Parsivel激光雨滴谱仪和CINRADA/SA多普勒雷达观测资料对冰雹云降水个例进行雨滴谱特征分析。
9、由于2DVD布设相对较少,且降雹过程不易捕获,目前国内少有使用2DVD观测降雹的研究。2DVD相较于测雹板具有可以直接获得冰雹下落速度、轴比、记录时间和更高的测量精度等优势;相较于雹雨分测计具有直接获得冰雹下落速度、轴比和更高的测量精度等优势;相较于一维雨滴谱仪具有直接获得冰雹轴比、二维信息和更高的测量精度等优势。降水粒子尺度谱在不同时间,不同地点,不同降水过程的形态都不一样,但具有统计规律。为拟合实际观测的雨滴谱资料,常采用M-P分布函数和分布函数,其中M-P分布对于稳定降水拟合效果较好,分布对云中起伏较大的降水拟合效果较好(赵震等,2005;郑姣恒和陈宝君,2007;李山山等,2020;刘
10、胜男和王改利,2020)。针对冰雹谱拟合的研究,牛生杰等(1999)指出冰雹谱服从幂指数分布和指数分布,岳治国和梁谷(2018)使用 M-P分布分段拟合降雹过程的雨滴谱和冰雹谱,本文采用 M-P分布分段拟合降雹时段的雨滴谱和冰雹谱。冰雹谱是冰雹微结构的重要特征之一,也是进行人工防雹作业的重点,了解冰雹的微结构特征及其形成机理,对于防雹减灾具有重要意义。本文利用巴彦淖尔市五原县气象局 2DVD 在 2020 年 7 月27日观测到的1次含雹混合降水资料,对含雹混合降水事件进行不同直径-速度大小和不同直径-轴比的粒子数量分布以及粒子谱参数随时间演变的分析,同时结合地面测雹板对2DVD观测的冰雹直径
11、进行了验证。由于冰雹落点位置的不确定,提前布设雨滴谱仪来捕获冰雹是比较困难的,在此次西北人研究试验项目中,利用2DVD实测的冰雹谱和雷达结合起来研究冰雹的发生发展与消亡,发现降雹过程的微物理特征,提高对降雹的认知,另外对冰雹谱进行M-P分布拟合的结果也将有利于改进模式微物理参数化。2 数据来源与方法介绍 2.1数据来源中国气象科学研究院组织实施西北地区人工影响天气研究试验项目,在巴彦淖尔地区开展防雹技术研究试验,为加强外场试验的科学观测,增设了二维雨滴谱仪(由成都信息工程大学支持),考虑到设备运行及维护,遂将该设备布设在了五原县气象局(108.2936E,41.0583N,53337),2DV
12、D 位置如图1所示。2019年8月10日于巴彦淖尔五原气象局对2DVD进行了校正,2020年7月27日2DVD观测到2次降水过程,其中19:36(北京时,下同)-20:16为含雹混合降水过程。2.2仪器介绍奥地利 Joanneum 研究所研制的 2DVD是一个粒子成像系统,2DVD基于两个垂直摆放且具有高度差的快速扫描相机,可以记录每个粒子的尺寸、下落速度和形状信息。许多学者已经详细地介绍了 2DVD 的结构、设计和测量原理等(Kruger and Krajewski,2002;Schnhuber et al,2007,2008)。2DVD内部光学系统提供了两个光平面,两个光平面穿过大约10
13、cm10 cm的虚拟测量区域,并投影到两个快速线扫描相机上。如图2(a)所示,(左侧)t1时刻水凝物在光面上方,所有光电探测器都被照亮,t2,t3,t4时刻水凝物进入光面,遮挡产生阴影被光电传感器记录。(右侧)水凝物轮廓近似构建。水凝物穿过光平面产生的阴影被相机记录,2DVD将此阴影换算为水凝物的轮廓宽度。如图2(b)所示,两个光平面具有约7 mm的高度差,2DVD利用高度差和相机记录水凝物进入和离开第一,第二个平面的时间来计算出粒子的下落速度,即:v=2d(t2-t0)+(t3-t1)(1)749高原气象42 卷式中:v(单位:m s-1)为下落末速度;d(约7 mm)为相机A和B的垂直距离
14、;t0,t1(t2,t3)分别为水凝物第一次进入和完全离开相机 A(B)探测平面的时刻。线扫描相机约每18 us对每个平面进行一次采样,当水凝物落入测量区域时,2DVD从两侧记录水凝物经过测量区域的多个线扫描结果,获得粒子图1巴彦淖尔高程图Fig.1Bayannur elevation map图22DVD工作原理示意图Fig.22DVD working principle diagram750刘彦等:基于二维雨滴谱仪的巴彦淖尔地区降雹谱个例分析3 期两侧的轮廓,最后通过两侧轮廓计算出水凝物的等效体积直径和扁平率。2DVD自1991年开始开发以来,现在已经发布了三个版本:the Classic
15、Tall 2DVD,the Low_Profile 2DVD,the Compact 2DVD,本次试验使用的是第三代the Compact 2DVD。2DVD对粒子的水平探测分辨率为 0.2 mm,垂直分辨率为 0.10.2 mm。2DVD观测数据存在质量问题,尤其是对小粒子的观测。Kruger and Krajewski(2002)指出2DVD存在粒子误匹配问题,粒子速度异常偏大与偏小,并且在强降水时粒子错误匹配现象更严重。尽管2DVD已经通过设计入口(25 cm)比相机视野(位于入口中间,宽约为10 cm)宽来尽量减少雨滴从边缘外壳飞溅进入采样区域造成的不良影响(Schnhuber et
16、 al,2008),但研究发现2DVD仍存在对粒子的虚假检测,尤其是直径小于0.6 mm的粒子(Larsen and Schnhuber,2018)。Thurai et al(2017)通过比较气象粒子光谱仪和2DVD观测的雨滴谱,指出2DVD对直径大于0.7 mm的粒子提供更准确的测量。2.3数据质控本文参照Kruger and Krajewski(2002),Schnhuber et al(2008)和 Luo et al(2021)研究,雨滴从2DVD设备边缘处飞溅进入采样区,以及两个相机对小粒子图像的错误匹配会造成观测的粒子下落速度异常偏大与偏小,因此可以先通过经验公式计算出雨滴下落末
17、速度来确定下落速度可信区间,再将下落速度在可信区间外的粒子剔除,只保留下落速度在可信区间内的粒子,从而保证研究数据的质量。雨滴下落末速度经验公式是Atlas et al(1973)在标准大气压下总结出的,且雨滴下落速度会受到空气密度的影响。参照Niu et al(2010)本地化Atlas下落末速度经验公式的方法,应用到本次试验,得到 2DVD 观测站点的雨滴直径-下落末速度,再由理论下落末速度确定粒子下落速度可信区间对此次降水过程 2DVD 实测数据进行质控。Atlas雨滴直径与下落末速度公式:vA=|0,4.323()D-0.003,9.65-10.3e-6D,D 0.03 mm0.03
18、mm 0.6 mm(2)式中:vA为Atlas下落末速度;D为粒子直径。根据本地空气密度调整Atlas雨滴直径与下落末速度公式:vBA=vA(a0)b-1(3)式中:vBA为根据五原县空气密度调整后的Atlas下落末速度;a为巴彦淖尔五原空气密度;0为标准大气压下空气密度1.23 kg m-3;b为半理论修正系数,本文取0.4。基于vBA的质控公式:|vM-vBA|cvBA(4)式中:vM为2DVD观测的粒子下落速度;参数c为可信速度区间,本文取0.6。2.4粒子谱分布函数冰雹谱服从幂指数分布和 M-P 分布(牛生杰等,1999),M-P分布是Marshall and Palmer(1948)
19、首次提出:N(D)=N0e-D(5)式中:N(D)(单位:mm-1 m-3)为单位尺度间隔,单位体积内的粒子数;N0(单位:mm-1 m-3)为浓度参数;(单位:mm)为尺度参数。阶矩法估计粒子谱分布参数是数值模式中常采用的方法,定义i阶矩Mi为:Mi=0DiN(D)dD(6)计算 M-P 分布参数需要 2 个矩量,本文使用3、4阶矩量(M34)估计M-P分布参数N0和,其表达式为:N0=128M533M44,=4M3M4(7)2.5雨滴谱特征参数本文以 0.1 mm为间隔将粒子直径 Di(i=0,1,2,110)划分为 0.1 11 mm,二维雨滴谱仪 N(D)计算公式:N(Di)=1tDj
20、=1ni1Sjvj(8)式中:t采样时间,本文为60 s;D为粒子直径分辨率,本文为0.1 mm;ni为采样时间内直径为Di的粒子总数;Sj为粒子的有效采样面积,本文取固定值0.01 m2;vj为粒子的下落速度(单位:m s-1)。粒子总数密度Nt计算公式:Nt=j=1n1Sjtvj(9)式中:n为采样时间内的粒子总数。质量加权平均直径Dm计算公式:Dm=M4M3(10)中值体积直径 D0是指将 2DVD在采样时间内测得的降水量等分为一半,其中一半粒子直径小于D0,另一半直径大于D0,计算公式如下:751高原气象42 卷DminD0D3iN(Di)D=12D0DmaxD3iN(Di)D(11)
21、式中:Dmin,Dmax分别为采样时间内的粒子直径最小值和最大值。雨强R计算公式:R=36001t6i=1nd3iSi(12)式中:di为2DVD观测的粒子直径。3 含雹混合降水事件分析 巴彦淖尔五原县含雹混合降水过程发生在2020 年 7 月 27 日 19:36-20:16,当天平均气压为889.1 hPa,空气密度为1.04 km m-3。地面观测人员在19:38观测到降雹,地面测雹板测得冰雹最大直径为11 mm,降雹持续时间为4 min,2DVD在降雹时间段记录到31个直径大于5 mm的冰雹粒子。3.1降水粒子下落速度分布粒子下落末速度指的是粒子在下落过程中受到的重力、空气浮力和空气阻
22、力达到平衡时的速度。影响冰雹下落末速度的因素包括冰雹初始降落高度、形状、密度、大小以及环境气流等(岳治国和梁谷,2018),国内外学者对不同地区、不同降雹过程的冰雹下落末速度进行了研究(徐家骝,1978;岳治国和梁谷,2018;Dieling et al,2020),Heymsfield et al(2018)使用不同研究数据拟合了冰雹下落末速度公式:Vt=13.95(0.1D)0.51(1000Pa)0.55,拟合Laurie et al(1960)数据Vt=12.3(0.1D)0.57(1000Pa)0.55,拟合Heymsfield et al(2014)数据Vt=6.1(0.1D)0.
23、72(1000Pa)0.55,拟合Heymsfield et al(2018)数据(13)式中:Vt为冰雹下落末速度;Pa为本地气压。在降雹过程中一般将直径大于5 mm的粒子作为冰雹粒子(岳治国和梁谷,2018;Heymsfield et al,2018),本文将此次过程中直径小于5 mm视作雨滴,直径大于5 mm视作冰雹粒子。图3展示了2DVD观测的降水粒子在不同直径大小和速度区间下的数量分布,其中,紫色实线表示Atlas直径-下落末速度,黑色实线表示根据五原县空气密度调整后的Atlas直径-下落末速度,黑色虚线表示调整后的 Atlas 下落末速度上下 60%边界,蓝色、橙色和绿色实线分别表
24、示Laurie(1960),Heymsfield et al(2014),Heymsfield et al(2018)数据拟合的冰雹下落末速度。五原县此次含雹混合降水事件的空气密度小于 Atlas实验的空气密度,雨滴下落受到的空气阻力变小,从而可见根据五原空气密度调整后的Atlas雨滴下落末速度大于Atlas在标准大气压下的雨滴下落末速度。由于粒子的下落速度会受到空气运动(上升气流、下沉气流)、湍流以及粒子的碰并和破碎的影响(Niu et al,2010),因此粒子的下落速度往往不等于其下落末速度,由图3可见,雨滴的下落速度分布在根据五原空气密度调整后的Atlas下落末速度上下方,并且其下落速
25、度整体呈现为大于下落末速度。直径大于5 mm的冰雹粒子下落速度为5.511.6 m s-1,同一直径上的冰雹下落速度不同,直径56 mm冰雹下落速度分布最宽。不同降雹过程,冰雹下落速度差距较大,本次降雹过程冰雹下落速度分布在由Laurie(1960)和 Heymsfield et al(2018)研究数据拟合的下落末速度之间,与Heymsfield et al(2014)研究数据拟合的下落末速度有交叉,冰雹下落速度整体上随着冰雹直径增大而增加。2DVD 仪器对小雨滴存在虚假检测与错误匹配,目前许多国内外研究都基于直径与下落末速度的关系来进行数据质控(Kruger and Krajewski,图
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 二维 雨滴 巴彦淖尔 地区 降雹谱 个例 分析 刘彦
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。