风廓线雷达通用数据格式
民用机场风廓线雷达系统技术规范
管理程序中国民用航空局空管行业管理办公室编号:AP-117-TM-2013-01下发日期:2013年8月6日民用机场风温廓线雷达系统技术规范目录第一章总则 (2)第二章系统构成及功能 (2)第一节一般规定 (2)第二节产品输出功能 (4)第三章系统性能 (6)第一节整体性能 (6)第二节各子系统性能 (7)第四章环境适应性 (12)附录一信号功率的谱矩及信噪比计算方法 (14)附录二风速、风向及Cn2计算方法 (15)附录三风温廓线雷达数据格式 (16)民用机场风温廓线雷达系统技术规范第一章总则第一条为规范民用机场风温廓线雷达系统(以下简称风温廓线雷达)的建设和运行,根据《中国民用航空气象工作规则》,制订本规范。
第二条本规范适用于中华人民共和国境内民用机场和军民合用机场民用部分(以下称民用机场)的风温廓线雷达系统的建设和运行。
第三条民用机场风温廓线雷达系统的构成、功能、性能和环境适应性等技术要求应当符合本规范。
第四条风温廓线雷达按照安装方式不同分为固定式和可移式两种。
可移式主要有车载可移式和方舱可移式两种。
第二章系统构成及功能第一节一般规定第五条风温廓线雷达主要由天线分系统、发射分系统、接收分系统、信号处理分系统、监控分系统、标定分系统、通讯分系统、数据处理及应用终端、配电分系统、RASS(Radio-Acoustic Sounding System,无线电-声探测系统)分系统构成。
第六条风温廓线雷达采用全相参脉冲多普勒体制。
第七条风温廓线雷达采用相控阵技术、全固态发射技术、数字接收机技术、脉冲压缩技术等。
第八条风温廓线雷达具有探测大气虚温的功能。
第九条风温廓线雷达应当具有实时采集功率谱数据,对功率谱数据进行信噪比计算、噪声电平计算、干扰信号剔除、杂波剔除、信号提取、谱矩估计、数据质量控制功能。
第十条风温廓线雷达应当具有生成风速、风向、垂直气流、Cn2(Refractive Index Structure Constant,大气折射率结构常数)、谱宽、信噪比等数据产品的功能。
风廓线雷达通用数据格式
附件:风廓线雷达通用数据格式(V1.2)2007年9月目录1. 文件名编码规则 (3)1.1 原始数据文件 (3)1.2 产品数据文件 (3)2.功率谱数据文件 (4)3. 径向数据文件 (4)3.1 文件组成单位 (4)3.2 文件框架 (4)3.3 文件结构 (6)4. 实时的采样高度上的产品数据文件 (10)4.1 文件组成单位 (10)4.2 文件框架 (10)4.3 文件结构 (10)5. 半小时平均的采样高度上的产品数据文件 (12)5.1 文件组成单位 (12)5.2 文件框架 (12)5.3 文件结构 (12)6. 一小时平均的采样高度数据文件 (13)6.1 文件组成单位 (13)6.2 文件框架 (13)6.3 文件结构 (13)附件一功率谱数据格式 (15)1. 文件名编码规则根据实际需求,建议使用长文件名命名法,对各类文件名进行约定。
文件名中的观测时间均为观测结束时间。
1.1 原始数据文件原始数据文件包括功率谱数据文件、瞬时径向谱数据文件,对于原始数据文件,建议每次观测生成一个文件,文件名具体命名方法如下:Z_RADR_I _IIiii_yyyyMMddhhmmss_O_WPRD_雷达型号_数据类型.TTT 其中:Z:国内交换文件;RADR:表示雷达资料;I:表示后面的IIiii为风廓线雷达站的区站号;IIiii:区站号(按地面气象站的区站号);yyyy:观测时间(年) (20**—);MM:观测时间(月) (01—12);dd:观测时间(日) (01—31);hh:观测时间(时) (00—23);mm:观测时间(分) (00—59);ss:观测时间(秒) (00—59);O:表示观测数据;WPRD:表示风廓线雷达资料;雷达型号:见表1;数据类型:功率谱数据文件用FFT表示;径向数据文件用RAD表示;TTT:当TTT = BIN时,表示二进制文件;当TTT = TXT时,表示文件格式为ASCII。
风廓线雷达C^(2)_(n)资料和L探空资料确定济南夏季边界层高度的对比研究
第44卷㊀第2期气象与环境科学Vol.44No.22021年3月Meteorological and Environmental SciencesMar.2021收稿日期:2019-02-19;修订日期:2019-07-03基金项目:山东省气象局气象科学技术重点研究项目(2014sdqxz05);山东省气象科学研究所数值天气预报应用技术开放基金(SDQXKF2015Z02)作者简介:王栋成(1969),男,山东莱州人,高级工程师,学士,从事大气环境与气候应用研究.E-mail:hjpj2008@通讯作者:汤子东(1969),男,山东蓬莱人,高级工程师,学士,从事气候预测和气候变化评估研究.E-mail:156****6207@王栋成,汤子东,邱粲,等.风廓线雷达C 2n 资料和L 探空资料确定济南夏季边界层高度的对比研究[J].气象与环境科学,2021,44(2):96-105.Wang Dongcheng,Tang Zidong,Qiu Can,et parison Study of Wind Profile Radar C 2n Data and L Sounding Data to Determine the Height of SummerBoundary Layer in Ji nan[J].Meteorological and Environmental Sciences,2021,44(2):96-105.doi:10.16765/ki.1673-7148.2021.02.012风廓线雷达C 2n资料和L 探空资料确定济南夏季边界层高度的对比研究王栋成1,汤子东1,邱㊀粲1,2,董旭光1,曹㊀洁1(1.山东省气候中心,济南250031;2.上海师范大学地理系,上海200234)㊀㊀摘㊀要:基于济南站边界层风廓线雷达(WPR )观测的大气折射率结构常数(C 2n ),采用偏离度法确定夏季白天的对流边界层高度(H CBL ),并与L 探空虚位温梯度法㊁L 探空湿度梯度法的结果进行对比㊂结果表明:(1)WPR 偏离度法判别的H CBL 成功率高,07㊁13㊁19时分别可达96.9%㊁100.0%㊁89.3%,尤其19时的判别成功率远高于L 探空两方法的46.4%㊁57.1%㊂(2)WPR 偏离度法与L 探空两方法判别的同时刻H CBL 结果一致性好,相关系数高㊂07时WPR 判别的边界层高度平均值分别较后两者的差+14.3㊁+9.3m ,13时的差+36.4㊁+21.6m ,19时的差-27.3㊁-25.4m ;与后两者的标准偏差,07时的分别为84.555㊁77.412m ,13时的分别为125.978㊁126.783m ,19时的分别为52.743㊁68.755m ;与后两者的相关系数,07时的分别为0.894㊁0.911,13时的均为0.980,19时的分别为0.980㊁0.969㊂WPR 能够将夏季H CBL 确定结果的绝对误差降至近似其垂直分辨率水平㊂(3)济南6㊁7月逐小时H CBL 最高可达3460.0㊁3220.0m ,小时平均最高分别为2205.0㊁2157.1m ;逐小时最高值在1316时均可出现,且以15时的最多;通常0715时的H CBL 呈波浪起伏式缓慢增高态势,1519时的H CBL 则是断崖式快速降低㊂关键词:对流边界层高度;风廓线雷达;L 探空;偏离度法;虚位温梯度法;湿度梯度法中图分类号:P49㊀㊀㊀㊀㊀文献标志码:A㊀㊀㊀㊀㊀文章编号:1673-7148(2021)02-0096-10引㊀言大气边界层的高度是污染物扩散㊁气候及大气等模式的一个重要输入参数,边界层高度的变化对数值预报中的物理过程㊁天气预报的诊断分析㊁环境空气质量的预报㊁空气污染物的监控与预测㊁大气环境容量的确定等也有相当重要的作用[1-3]㊂基于不同的物理量,判定边界层高度的依据不同[4]㊂目前,边界层高度的确定方法主要有两类:一是基于参数化和简单的模型计算,二是基于廓线测量[5]㊂参数化方法主要包括国标法[6,7]㊁罗氏法[7,8]㊁联合频率法-罗氏法B [7]㊁萨默斯(Summents)经验公式法[7]㊁最小二乘曲线拟合法[9]等㊂参数化方法实用性强,但由于边界层高度变化的复杂性及以地面常规观测资料判断边界层变化的局限性,各种参数化方案的准确性均有待进一步验证和改进[5]㊂廓线测量法主要包括干绝热法[7,10-11]㊁实测法[7,12]㊁风速极值法[13-15]㊁位温法[4,7](虚位温梯度法[15-17])㊁湿度梯度法[1,18]㊁湍流能量法[4,15]㊁综合评定法[7]等㊂廓线测量是确定大气边界层结构最常用的方法,可较准确地判断边界层高度,但大气边界层的特性和结构资料并不是常规气象观测的内容,因此对它的研究只能依赖于特定的试验项目或每日2~3次的L 探空,观测资料在时间和空间上都非常有限[19]㊂近年来,随着快速发展的风廓线雷达(简称WPR)[20,21]㊁激光雷达[22]㊁微波辐射计[23]等技术方㊀第2期王栋成等:风廓线雷达C2n资料和L探空资料确定济南夏季边界层高度的对比研究法的应用,可连续系统且高时空分辨率获取大气边界层各参数的分布数据,边界层的观测研究已进入了一个新时代[15]㊂其中,WPR资料能够快速㊁细致地反映边界层的结构㊁厚度㊁湍流演变过程特征[4,22-28]㊂WPR产品数据中的大气折射率结构常数(C2n)是表征湍流变化强弱的一个重要参数㊂C2n 与气温㊁气压㊁湿度密切相关,且随着高度的增加呈指数递减规律[20],并表现为在边界层顶出现极大值或者偏离正常值幅度大值[15-16,20],因此从C2n垂直廓线的时间序列中可以判定边界层高度[20]㊂现阶段基于WPR资料进行边界层高度确定的研究尚较少,且多属方法探索㊁个例或短期分析㊂如Bianco 等[17]提出了基于C2n等综合指标的模糊逻辑判别方法判别对流边界层高度(以下简称H CBL),并给出了3天的结果对比㊂Bianco等[29]对该方法进行了改进,并提供了个例验证㊂Allabakash等[30]采用改进的模糊逻辑方法给出了Gadanki测站两年各月07 20时H CBL日变化㊂国内曹晓彦①㊁王敏仲[31]㊁蒋德海[32]㊁张坚[15]㊁金莉莉[4]㊁王栋成[33]等,对C2n极大值法及求导后放大其跃变特征的偏离度法进行了方法研究或个例验证,均属单一指标判别㊂模糊逻辑法判别指标和规则较多,且相对复杂,不易执行,有的指标受杂波干扰较大,且对于WPR原始数据和产品数据质量要求高,适用性尚存局限㊂C2n极大值法采用单一指标判断,易受杂波干扰㊁信噪比异常㊁C2n值太小廓线突变特征不易找到等影响,以C2n廓线值直接判断边界层高度难度较大㊂偏离度法对C2n求导以放大其变量特征,根据C2n廓线总体判断其连续突变或跃变特征来判断边界层高度,方法简易可行,与实测对比准确度高,且可避免C2n极大值单一高度层突变有可能属异常值而简单地作为边界层顶的判断失误㊂在H CBL判别时,偏离度法总体优于极大值法[33]㊂可见,确定大气边界层高度的方法有很多,不同方法得出的边界层高度也存在差异,且不同探测手段得到的气象要素也有可能对应不同的层结[15]㊂利用长期实测数据对比验证,则是判断这些方法优劣㊁给出准确的边界层高度,以及分析其演变规律的最好途径㊂本文基于济南站2014年67月40天的0719时的WPR观测的C2n数据,采用偏离度法[33]确定H CBL,并与同时段的基于L探空资料虚位温梯度法[15-17]和湿度梯度法[1,18]确定的结果进行对比,进而探索以WPR高时空分辨率资料确定较长期的H CBL并分析其时空演变规律,以期为改进边界层高度确定参数化法的缺陷,研究边界层精细化结构与大气污染的关系,提高环境空气质量预报准确率等提供依据㊂1㊀资料与方法1.1㊀资㊀料研究采用的2014年6月130日和7月1524日逐小时整点后18min㊁高度层1003940m的C2n数据,来自济南WPR站实时每6min产品数据ROBS 文件㊂该站固定式边界层WPR型号为CLC-11-D 型,主要技术指标为:时间分辨率ɤ6min,有效数据最低探测高度100m㊁最高探测高度ȡ3km,雷达工作模式包括低模式(分辨率30㊁60㊁120㊁240m)㊁高模式(分辨率为60㊁120㊁240m),输出产品数据的高度分辨率为60m(820m以下)㊁120m(8202020 m)㊁240m(2020m以上)三种㊂采用的2014年6月130日和7月1524日每日07㊁13㊁19时整点后15min㊁高度层04000m的气温㊁相对湿度等数据,来自济南L波段探空雷达记录文件㊂L探空数据均经质量控制,符合中国气象局‘高空气象探测规范“㊂L探空秒级和分钟级数据的高度分辨率分别为6~7m㊁400m㊂考虑L探空整点后15min开始探测,气球按照400m/min的速度上升,10min的高度足以涵盖WPR探测每6min㊁3980m的高度范围,因此本研究采用的整点后18min的ROBS数据与整点后15min的L探空数据的时空匹配一致㊂每10m高度层间隔的气温㊁相对湿度资料,利用L波段高空气象探测系统数据处理软件内插获取[35]㊂采用的地面气象逐小时数据来自章丘气象站,数据均经自动㊁人工审核与质量控制,符合中国气象局‘地面气象观测规范“要求㊂章丘气象站与济南WPR观测站㊁济南L探空观测站属于相同站场,观测场经纬度㊁海拔高度基本一致,本文统称为济南站㊂1.2㊀C2n数据质量控制(1)济南WPR站2014年C2n产品数据非原始值,而是10ˑlog(C2n)且取整后的结果,C2n突变值被削平后无法再用以判断H CBL㊂本研究基于该站经质量控制后业务上传的原始产品数据和WPR79①曹晓彦,李炬,张京江,等.由风廓线数据分析北京地区夏季边界层高度的特征[C]//中国气象学会2007年年会气象综合探测技术分会场论文集.广州.2007:68-77.气象与环境科学第44卷设备参数等,采用中国气象局‘风廓线雷达通用数据格式(V1.2)“‘风廓线仪功能规格需求书(试行)“‘风廓线雷达观测规定(试行)“及QX/T78 2007‘风廓线雷达信号处理规范“等推荐的雷达气象散射方程[20],对该站2014年逐日逐6min的C2n 数据进行重新计算,并经一致性平均㊁时空连续性检验质量控制[20,21],按通用数据格式的要求统一定标输出㊂(2)实测C2n廓线数据常出现个别高度层值缺测㊁为0㊁连续2层以上为同一个数值㊁异常大或异常小值等现象㊂由于气象场的变化具有连续性[35],因此可根据上下层和前后3个时次的数据,进行奇异点判断,剔除异常数据,并用前后或上下点平均值来替代该点的数据,从而得到较为完整连续的随高度变化的C2n廓线㊂(3)因WPR高模式㊁低模式探测衔接问题, 940㊁2020m层C2n数据常为较大值,导致单一特征高度层的C2n偏大㊂对此采用上㊁下层值内插方法暂替代异常值㊂(4)除单时次C2n数据的质量控制外,该站2014年8月10日10月17日统计的1km以下各高度层风向㊁风速均为异常㊂分析发现,其对应的C2n数据也存在异常㊂因此,研究避开了该时段㊂1.3㊀对比方法济南站2014年6月130日和7月1524日,每日07:15㊁13:15和19:15均有L探空记录,共120组探空记录可进行H CBL的判别,并可与WPR观测的同时段整点后18min的C2n数据判别结果进行对比㊂Bianco[17]㊁Angevine[18]等以L探空虚位温廓线上第一个梯度值大于0.5ħ/(60m)的点的高度值定义为边界层高度㊂分析济南站实测资料发现,L 探空结果有许多是达不到该判别指标的,07时的30个有效样本中有21个达不到(占比70.0%),13时的32个有效样本中有18个达不到(占比56.3%), 19时的13个有效样本中有8个达不到(占比61.5%)㊂对此提出修正指标,以实现无遗漏地判别H CBL,即对于未出现大于0.5ħ/(60m)的情形,以虚位温梯度廓线上最大值对应的高度值定义为边界层高度㊂本文将该虚位温梯度法称为方法一㊂魏浩[1]㊁Angevine[18]等以L探空湿度梯度廓线上梯度最大值的点的高度值,定义为边界层高度㊂本文将湿度梯度法称为方法二㊂王栋成等[33]将WPR观测的10ˑlog(C2n)廓线的实测值与同时刻拟合曲线值的差定义为偏离度㊂实际偏离度最大值对应的高度层判断为对流边界层顶高度,即偏离度(D)=实测值-拟合值,H CBL=h{D max}㊂本文将偏离度法称为方法三㊂因WPR对雨滴很敏感[18,20],当降水出现时C2n 明显增大[36,37],所以当有雨时边界层高度无法用WPR数据测量[18]㊂本研究也证实,地面或空中降雨时段㊁各高度层相对湿度连续较大时段(多高度层连续大于85%),WPR观测的C2n数据均不适用于判定H CBL,同时L探空数据也不适用于判定H CBL㊂因此,首先剔除这两种因素影响的样本,得到同时适用于三种方法判定H CBL的样本数分别为07时32组㊁13时34组㊁19时28组(表1)㊂然后,分别用三种方法进行H CBL的识别判定,其中判断成功的称为有效样本(表1)㊂最后,再对有效样本中三种方法同时刻均有H CBL判定结果的样本进行对比,结果为07时30组㊁13时31组㊁19时11组㊂由表1样本统计结果可见,每日早中晚120组总样本,剔除降雨和湿度较大的样本后为94组㊂其中,方法三的H CBL 判别成功率最高,07㊁13㊁19时,以及合计判别成功率分别为96.9%㊁100.0%㊁89.3%㊁95.7%;方法一和方法二在07时㊁13时的成功率均略低于方法三的,但在19时则差异较大,分别只有46.4%㊁57.1%㊂早晨和中午时段,是对流边界层的发展和升高阶段,边界层顶的湍流㊁大气折射率㊁温湿结构等参数的变化存在显著极大值或突变值,易于被WPR和L探空雷达捕获到㊂在傍晚,对流边界层则处于崩溃阶段,边界层顶的大气结构与参数极大值或突变值多数情形下并不十分清晰稳定,由于L探空气球以400m/min上升,升速过快,捕获的难度很大,WPR的观测能力显然远高于L探空的㊂2㊀结果分析2.1㊀个例分析(1)三种方法确定的结果对比以济南WPR站2014年7月18日为例,对比三种方法确定的H CBL结果㊂该日日出㊁日落时间分别为05:05㊁19:25,10m高度0719时逐小时风速为0.5~4.1m/s,对流边界层气象条件以B㊁C类稳定度为主㊂根据三种方法确定的H CBL高度(图1㊁表2),方法三确定的07时H CBL分别较方法一㊁方法二的高20㊁20m,13时的H CBL分别低80㊁110m,19时的H CBL分别低40㊁20m㊂考虑WPR输出产品数据高度分辨率820m以下为60m㊁2020m以上为240m,L 探空的为10m级,可见,07㊁13㊁19时WPR的C2n偏离度法判别结果与L探空两方法的结果基本一致,差异均很小㊂89㊀第2期王栋成等:风廓线雷达C 2n 资料和L 探空资料确定济南夏季边界层高度的对比研究表1㊀2014年济南夏季三种方法判断H CBL 的能力对比时次方法总样本数/组剔除降雨和湿度较大数据后的样本数/组三种方法各自判断的有效样本数/组(成功率)三种方法同时刻的可对比样本数/组方法一403230(93.8%)07时方法二403230(93.8%)30方法三403231(96.9%)方法一403432(94.1%)13时方法二403431(91.2%)31方法三403434(100.0%)方法一402813(46.4%)19时方法二402816(57.1%)11方法三402825(89.3%)方法一1209475(79.8%)合计方法二1209477(81.9%)72方法三1209490(95.7%)㊀㊀另外,07时和19时在20003000m 高度层,即边界层和自由大气层之间均有一层覆盖逆温层(又称夹卷层)[38],逆温层厚度㊁虚位温梯度㊁湿度梯度㊁C 2n 偏离度均可能较大,但它们各自的突变或跃变层高度并不一致,在判断边界层高度时应注意识别与剔除㊂图1㊀L 探空虚位温梯度法(a )㊁湿度梯度法(b )㊁WPR 偏离度法(c )三种方法确定的济南站2014年7月18日07时(a1a6)㊁13时(b1b6)㊁19时(c1c6)H CBL 对比99气象与环境科学第44卷表2㊀三种方法确定的济南站2014年7月18日H CBL结果对比时次方法一H CBL/m L探空虚位温梯度/(ħ/(60m))方法二H CBL/m L探空湿度梯度/(%/(60m))方法三H CBL/m WPR偏离度/m-2/307时6200.362056407.87 13时21000.62130302020 5.92 19时8000.17806760 4.49㊀㊀(2)逐小时H CBL的日变化基于方法三确定的济南站2014年7月18日0719时逐小时H CBL变化见图2,并同时给出了该日地表温度㊁气温日变化和方法一㊁方法二确定的早中晚H CBL㊂由图2可见,该日H CBL自07时的640m逐渐波浪起伏式增高,至午后1516时达到最高值2740m,其后快速降低,至19时为760m㊂与通常的H CBL最高值出现在15时略有不同,该日最高值持续至16时㊂分析原因是该日在太阳辐射作用下,地表长波辐射较强,地表温度最高达64.3ħ(13时),近地面气温1316时滞后响应且缓慢持续增高达36.0ħ(16时),大气边界层内的湍流运动使近地层热量向上传递,上下层热交换充分且对流较强㊁热浮力和动力抬升作用滞后所致㊂同时也显示, H CBL最高值对近地面热通量作用的响应时间尺度为1h或更短,而对地表热通量作用的响应时间尺度则可滞后达3h,这与张坚[15]㊁卢萍[39]等的研究结果一致㊂图2㊀偏离度法确定的济南WPR站2014年7月18日白天的逐小时H CBL变化2.2㊀夏季67月三种方法可对比有效样本的统计分析基于济南WPR站和L探空站2014年6月1 30日和7月1524日共40天的每日07㊁13㊁19时同时段(WPR观测的C2n记录为整点后18min,L探空记录为整点后15min)的资料进行H CBL确定,三种方法同时刻的可对比有效样本数见表1㊁表3,三种方法确定的H CBL结果对比见表3㊁图3㊁图4㊂(1)三种方法确定的结果统计分析三种方法确定的H CBL均值㊁最小值㊁最大值虽个别有差异但差异较小,总体一致性很好,表明WPR实测和L探空实测结果均能很好地反映大气边界层高度变化的实际状况㊂方法三确定的济南2014年67月白天平均H CBL为1124.4m,其中07时的为510.0m,13时的为1884.5m,19时的为658.2m㊂07时,方法三的平均值较前两者的差值分别为+14.3㊁+9.3m,最小值的差值分别为+20㊁0m,最大值的差值分别为-20㊁-20m㊂30组有效样本中,方法三与前两者的差值范围在-180~+250m(76.7%样本的差值在-80~+100m)㊁-180~+200m (80.0%样本的差值在-100~+100m)㊂方法一和方法二的差值范围很小,为-50~+80m,且83.3%样本的差值在-20~+30m,平均值差仅5.0m㊂13时,方法三的平均值较前两者的差值分别为+36.4㊁+21.6m,最小值的差值分别为+20㊁+10m,最大值的差值分别为+50㊁+80m㊂31组有效样本中,方法三与前两者的差值范围在-170~+490m (70.0%样本的差值在-80~+100m)㊁-270~+310m (63.3%样本的差值在-100~+100m)㊂方法一和方法二的差值范围为-380~+290m,且80.6%样本的差值在-40~+80m,平均值差仅14.8m㊂19时,方法三的平均值则较前两者的差值分别为-27.3㊁-25.4m,最小值的差值分别为-70㊁-30m,最大值的差值分别为+40㊁+40m㊂11组有效样本中,方法三与前两者的差值范围在-130~ +40m(90.9%样本的差值在-70~+40m)㊁-150~ +50m(81.8%样本的差值在-60~+50m)㊂方法一001㊀第2期王栋成等:风廓线雷达C2n资料和L探空资料确定济南夏季边界层高度的对比研究和方法二的差值为-60~+110m,且90.9%样本的差值在-60~+10m,平均值差仅1.9m㊂表3㊀三种方法判定的2014年济南夏季同时刻H CBL统计结果对比时间方法三种方法可对比组数/组平均H CBL/m最小H CBL/m(样本序号,日/月)最大H CBL/m(样本序号,日/月)方法一495.7200(24,24/6)960(15,15/6) 07时方法二30500.7220(24,24/6)960(15,15/6)方法三510.0220(24,24/6)940(15,15/6)方法一1848.1680(36,20/7)2930(29,29/6) 13时方法二311862.9690(36,20/7)2900(29,29/6)方法三1884.5700(36,20/7)2980(29,29/6)方法一685.5230(35,19/7)1020(4,4/6) 19时方法二11683.6190(35,19/7)1020(4,4/6)方法三658.2160(35,19/7)1060(4,4/6)方法一1106.9200(24,24/6)2930(29,29/6)合计方法二721115.1190(35,19/7)2900(29,29/6)方法三1124.4160(35,19/7)2980(29,29/6)图3㊀三种方法判定的2014年济南夏季07时(a)㊁13时(b)㊁19时(c)H CBL结果逐样本对比㊀㊀(2)三种方法结果的相关性分析图4为三种方法结果的相关性对比㊂由图4可见,三种方法确定的H CBL结果的统计特征总体一致,且相关系数较高,均通过0.01的显著性检验,证实基于WPR的C2n实测数据和L探空温湿数据确定的H CBL结果一致㊂07时,方法三与方法一㊁方法二的相关系数R 均较高,为0.894㊁0.911,标准偏差σ分别为84.555㊁77.412m;方法一和方法二的R=0.989,σ=25.831m㊂13时,方法三与方法一㊁方法二的R均为0.980,σ分别为125.978㊁126.783m;方法一和方法二的R=0.983,σ=108.409m㊂19时,方法三与方法一㊁方法二的R分别为0.980㊁0.969,σ分别为52.743㊁68.755m;方法一和方法二的R=0.986,σ=42.619m㊂早中晚全部72组样本,方法三与方法一㊁方法二的R分别为0.992㊁0.992,σ分别为102.418㊁100.790m;方法一和方法二的R=0.995,σ=74.386m㊂以上得出的相关系数R,与曹晓彦等①的0.76㊁张坚[15]的季风间断期0.985和季风活跃期0.920㊁Bianco等[17]的0.91相比,属于比较理想的情况㊂分析原因是,济南站WPR观测的C2n数据和L探空数据质量均较好且观测样本足够多,两者观测时段同步匹配性好,且因属相同站场,探空测量气球漂移范围较小,L探空数据和WPR数据均能代表测站上空的实际情况㊂考虑到WPR数据输出60㊁120㊁240m的高度分辨率,以及WPR探测与L探空两种遥感手段探测原理的不同,方法三的偏离度法已能够将H CBL确定结果的绝对误差降至近似WPR垂直分辨率水平㊂101气象与环境科学第44卷图4㊀三种方法确定的2014年济南夏季07时(a1a3)㊁13时(b1b3)㊁19时(c1c3)㊁全部样本(d1d3)H CBL 的相关性对比单位:m2.3㊀济南夏季67月H CBL 日变化以WPR 观测数据C 2n 为基础,采用方法三确定了济南站2014年67月共40天的0719时逐小时H CBL ,统计结果见表4㊁图5,其中典型数据较全的逐小时变化见图5㊂(1)各典型日H CBL 逐小时值是波浪起伏式变化的,而非平均值那样的持续增高或降低;0715时的H CBL 总体呈现缓慢增高态势,而15时后有一些H CBL是断崖式降低的,如6月3㊁6㊁24㊁28日和7月19㊁22日等,而另一些则是延续滞后至1617时才迅速降低201㊀第2期王栋成等:风廓线雷达C2n资料和L探空资料确定济南夏季边界层高度的对比研究的,如6月17㊁23㊁27日和7月15㊁16日等㊂原因是对流边界层的形成阶段是由于热力抬升增强导致的热泡式起伏增长,而对流边界层的崩溃阶段则是热力抬升减弱导致坍塌式快速降低[16]㊂(2)济南站夏季6月小时H CBL最高为3460.0m,出现在6月6日15时(图5a);逐小时H CBL最高值在1316时均可出现,其中以15时的最多,为18天,其次是14时的11天㊁13时的8天㊁16时的3天㊂7月小时H CBL最高为3220.0m,出现在7月17日15时(图5b);逐小时H CBL最高值在1416时均可出现,其中以15时的最多,为4天,其次是16时的3天㊁14时的2天㊂6月㊁7月平均H CBL日变化规律基本一致,均是早晚低㊁中午高,最高值均出现在15时,0715时H CBL均为缓慢增高,1519时H CBL则快速降低㊂6月㊁7月的07时平均高度分别为498.2㊁526.7m,19时则分别为586.3㊁573.3m,15时最高分别为2205.0㊁2157.1m㊂7月的最高值略低于6月的,可能是时段短㊁样本数少所致㊂(3)WPR与L探空确定的H CBL平均值早㊁中㊁晚均较为一致㊂6月07㊁13㊁19时的差值分别为+16.2㊁+40.7㊁-82.3m,7月07㊁13㊁19时的差值分别为+3.7㊁0㊁-86.7m㊂表4㊀偏离度法确定的济南夏季6月㊁7月H CBL统计结果m 时次07080910111213141516171819 6月平均498.2/22632.2/23830.9/221117.1/211480.0/211692.7/221975.0/242197.1/212205.0/241637.1/211316.7/18816.7/18586.3/19 7月平均526.7/9722.5/8982.9/71060.0/61400.0/61551.4/71600.0/102088.6/72157.1/71900.0/71711.4/71210.0/8573.3/6㊀注:斜线后数据为样本数图5㊀偏离度法确定的2014年济南6月(a)和7月(b)H CBL月均值和典型日值的逐小时变化3㊀结论与讨论(1)基于WPR观测的C2n资料判别边界层高度的成功率高,尤其在19时,成功率远高于L探空方法的成功率;WPR观测的时间分辨率可达6min,也远高于L探空的㊂WPR与L探空资料判别的同时刻对流边界层高度结果一致性好,相关系数高,证实基于WPR判别边界层高度准确㊁可行㊂这些优势可使WPR实时㊁连续㊁高分辨率地判别边界层顶高度,精细化认知边界层的空间结构和演变规律,且对改进边界层高度的参数化方法和提高环境空气质量预报准确率均有很好的参考价值㊂(2)济南6月㊁7月对流边界层高度15时最高值可达3460.0㊁3220.0m,平均值可达2205.0㊁2157.1m㊂0715时的对流边界层高度为波浪起伏式缓慢增高态势,1519时的变化则是断崖式降低,这与热力抬升的增强和减弱过程直接相关㊂(3)虽然WPR观测数据具有高时间分辨率的优点,但其空间分辨率不及L探空数据的高㊂而且,WPR探测能力和产品数据的质量,仍是制约确定大气边界层高度的关键因素㊂当有降水时或各高度层的相对湿度连续较大时,或夜间边界层高度低于WPR最小探测高度100m时,或冬季探测能力低于边界层顶高度时,或低层信号受地面杂波干扰时,或数据缺测较多㊁异常情况大量出现时,不同观测模式㊁数据运算方法等自身存在的诸多问题,均可能会导致无法判定边界层高度的情形㊂当然,相比3km 探测高度的WPR,6km探测高度的WPR探测能力进一步提高,可一定程度地解决大气边界层高度有时无法判定的问题㊂(4)WPR产品数据文件中的C2n应严格统一为通用数据格式且统一标定方法,以使单站或多站边界层高度确定结果具有可比性㊂应当对高㊁低探测模式C2n数据衔接方法进行改进,使C2n廓线具有时301气象与环境科学第44卷空连续性和一致性㊂应尽可能地提高输出产品数据的空间分辨率,以适应精细化边界层结构应用的实际需求㊂基于C2n㊁信噪比(SNR)㊁速度谱宽㊁垂直径向速度等产品数据的综合指标判定边界层高度,应是未来研究发展的方向㊂参考文献[1]魏浩,胡明宝,艾未华.小波变换在WPR探测大气边界层高度中的应用研究[J].热带气象学报,2015,31(6):811-820. [2]滑申冰,师华定,王堃,等.20162017年冬季华北地区一次重污染过程的气象条件分析[J].气象与环境科学,2018,41(4):47-53. [3]姚琳,温新龙,沈竞.江西山地风电场风速数值模拟方法研究[J].气象与环境科学,2018,41(3):120-125.[4]金莉莉,李振杰,何清,等.乌鲁木齐市边界层日变化特征[J].气象科技,2017,45(4):686-697.[5]王耀庭,李威,张小玲,等.北京城区夏季静稳天气下大气边界层与大气污染的关系[J].环境科学研究,2012,25(10):1092-1098. [6]中华人民共和国标准.GB/T13201-91,制定地方大气污染物排放标准的技术方法[S].北京:标准出版社,1992:17-18. [7]廖国莲.大气混合层厚度的计算方法及影响因子[J].中山大学研究生学刊(自然科学医学版),2005,26(4):66-73. [8]马金,郑向东.混合层厚度的经验计算及与探空观测对比分析[J].应用气象学报,2011,2(5):567-576.[9]程水源,席德立,张宝宁,等.大气混合层高度的确定与计算方法研究[J].中国环境科学,1997,17(6):512-516.[10]吴祖常,董保群.我国陆域大气最大混合层厚度的地理分布与季节变化[J].科技通报,1998,14(3):158-163. [11]李二杰,刘晓慧,李洋,等.一次重污染过程及其边界层气象特征量分析[J].干旱气象,2015,33(5):856-860. [12]程水源,席德立.关于确定大气混合层高度的几种方法[J].环境科学进展,1997,5(4):63-67.[13]赵鸣,苗曼倩,王彦昌.边界层气象学教程[M].北京:气象出版社,1991:218-219.[14]洪钟祥,钱敏伟,胡非.由地基遥感资料确定大气边界层特征[J].大气科学,1998,22(4):613-624.[15]张坚.WPR在南海季风边界层结构及深圳海陆风研究中的应用[D].南京:南京信息工程大学,2013:1-35.[16]Hyun Y K,Kim K E,Ha K J.A comparison of methods to estimatethe height of stable boundary layer over a temperate grass-land[J].Agricultural and Forest Meteorology,2005,132(1/2):132-142.[17]Bianco L,Wilczak J M.Convective boundary layer depth:Improvedmeasurement by doppler radar wind profiler using fuzzy logic meth-ods[J].Journal of Atmospheric and Oceanic Technology,2002,19(11):1745-1758.[18]Angevine W M,White A B,Avery S K.Boundary-layer depth andentrainment zone characterization with a boundary layer profiler[J].Boundary-Layer Meteorology,1994,68(4):375-385. [19]Seibert P,Beyrich F,Gryning S E,et al.Review and intercompari-son of operational methods for the determination of the mixing height [J].Atmospheric Environment,2000,34(7):1001-1027.[20]何平.相控阵风廓线雷达[M].北京:气象出版社,2006:46-115.[21]胡明宝.WPR探测与应用[M].北京:气象出版社,2015:95-113.[22]李红,马媛媛,杨毅.基于激光雷达资料的小波变换法反演边界层高度的方法[J].干旱气象,2015,33(1):78-88. [23]党张利,张京朋,曲宗希,等.微波辐射计观测数据在降水预报中的应用[J].干旱气象,2015,33(2):340-343. [24]陈浩君,黄兴友,王亚东,等.上海TWP3型边界层WPR探测性能评估[J].气象科技,2015,43(3):355-360.[25]汪学渊,任雍,李栋.闽北地区边界层移动风廓线雷达对比试验评估[J].气象与环境科学,2014,37(3):108-113. [26]张丽,李磊.深圳市空气污染气象条件标准体系的研究[J].气象与环境科学,2016,39(3):112-116.[27]陆琛莉,李海军,张雪慧,等.2013年12月浙北北部两次重度霾过程的对比分析[J].气象与环境科学,2018,41(1):47-55. [28]戴竹君,高辉,李力,等.新型探测资料在2014年春节污染天气中的应用[J].气象与环境科学,2017,40(1):78-86. [29]Bianco L,Wilczak J M,White A B.Convective boundary layer depthestimation from wind profilers:Statistical comparison between an au-tomated algorithm and expert estimations[J].Journal of Atmospheric and Oceanic Technology,2008,25(8):1397-1413. [30]Allabakash S,Yasodha P,Bianco L,et al.Improved boundarylayerheight measurement using a fuzzy logic method:Diurnal and season-al variabilities of the convective boundary layer over a tropical sta-tion[J].J Geophys Res Atmos,2017,122:9211-9232. [31]王敏仲,魏文寿,何清,等.WPR对塔克拉玛干沙漠晴天边界层的探测分析[J].气象,2012,38(5):577-584.[32]蒋德海,王成刚,吴兑,等.利用WPR资料对广州地区边界层日变化特征的分析研究[J].热带气象学报,2013,29(1): 129-135.[33]王栋成,邱粲,董旭光,等.偏离度法确定济南夏季对流边界层高度研究[J].气象与环境学报,2019,35(3):29-36. [34]李伟,刘凤琴,徐伟,等.L波段高空气象探测系统软件[J].气象科技,2008,36(2):237-239.[35]周雯嫣,王凌震,胡帆,等.干扰对边界层风廓线的影响及对策[J].气象水文海洋仪器,2011(2):90-95.[36]史珺,赵玉洁,王庆元,等.风廓线雷达在一次短时暴雨过程中的应用[J].气象与环境科学,2017,40(4):83-89. [37]张彩英.基于风廓线雷达资料的暴雪天气过程分析[J].气象与环境科学,2016,39(4):80-85.[38]Stull R B.边界层气象学导论[M].杨长新,译.青岛:青岛海洋大学,1991:478-519.[39]卢萍,杨康权,李英.温江站夏季大气边界层垂直结构特征[J].沙漠与绿洲气象,2017,11(5):70-76.401。
MICAPS数据格式(全部)剖析
数据: 区站号 经度 纬度 跋海高度 单站内容长度 风速 第二层气压 除风向风速外缺值时整个层次取消掉,风向风速缺值用 注:单站内容长度为层数× 6
第一层气压 高度 9999 表示
温度
露点 风向
例子: diamond 5 98 年 08 月 21 日 08 时温度对数压力图 98 08 21 08 348 53068 112.00 43.65 966 96
一位小数,当降水量大于 1 时只填整数。
-2
表示填 24 小时降水量。 当降水量小于 1mm时不填,大于等于 1mm时只填整数。
-3
表示填温度。只填整数。
数据: 区站号(长整数) 符串)
经度
纬度 拔海高度(均为浮点数)
站点值 1 站点值 2 (均为字
注意按照 MICAPS3.2 扩展的数据格式定义,在 6 小时雨量中, 0.0 表示微量降水,而 不是无降水,上述类别数据填图属性中设置小数位数不起作用。考虑到实际业务中使用的 数据格式, 修改为 0.0 时表示无降水, 大于 0 并且小于 0.1 为微量降水。 任意使用负值或 9999 表示降水为 0,可能会导致数据分析中出现异常结果。
A1.3 第三类数据格式:通用填图和离散点等值线 (注意:数据中一定不能有经纬度相同的站点,否则生成三角网时将出错)
文件头: diamond 3 数据说明(字符串) 年 月 日 时次 等值线条数(均为整数) 等值线值 1 等值线值 2 数)
层次 平滑系数
加粗线值(均为浮点
剪切区域边缘线上的点数(整数) 边缘线上各点的经度值 1 纬度值 1 经度值 2 纬度值 2
A1.4 第四类数据格式 : 格点数据
文件头: diamond 4
数据说明(字符串)
新一代天气雷达产品数据格式
长度 2 2 4 4 2 2 2
名称 productCode productDate productTime fileLength radarCode receiveCode blockCount
含义 产品代号 产品生成日期,1970 年 1 月 1 日以来的天数 产品生成时间,00:00 以来的秒数 文件包含的字节数 雷达站代号 接收站代号 文件包含的数据块数目
扫描中心点 I 坐标
组成一条径向数据的库数
扫描中心点 I 坐标
比例因子
径向数据条双字节数
径向数据条数
扫描起始角
径向数据条 1 径向数据条 2
…… 径向数据条 n
扫描展角 径向小块 1、2 径向小块 3、4
…… 径向小块 n-1、n
包识别码 1 包识别码 2 包识别码 3 X 坐标起始位置
Y 坐标起始位置
3.1.1.3.1 首尾相连的数据包(Linked Vector Packet)
长度 2 2 2 2 2 2 2
名称 packetCode blockLength colorLevel
2 2
含义
包识别码, 6/9 数据块字节数
矢量颜色等级, 可选域,packetCode=9 时有此域 第 1 点 X 坐标 第 1 点 Y 坐标 第 2 点 X 坐标 第 2 点 Y 坐标 … 第 n 点 X 坐标 第 n 点 Y 坐标
径向小块 0、1,高字节为小块 0,低字节为小块 1
2 data[2]、data[3]
径向小块 2、3
…
…
2 data[n-2]、data[n-1] 径向小块 n-2、n-1
3.1.1.1.1.1 径向小块
长度
名称
大气中颗粒物检测搭配风廓线激光雷达技术要求
风廓线激光雷达1.1基本要求基于光学脉冲相干多普勒频移检测原理,实现中下层对流层(包括大气边界层)三维风场的精细化探测,包括但不限于风廓线、湍流强度及气溶胶后向散射系数、消光系数等的连续探测。
1.2技术参数要求1)★探测范围:45m~6000m ;2)★径向探测库:15m/30m可设置,探测层数不小于300层;3)★激光发射波长:1550nm,基于相干多普勒探测原理,激光人眼安全符合 EN 60825-1:2014 Class 1M标准;4)脉冲宽度:100ns~400ns可调;5)★单脉冲能量≥160µJ;6)★径向速度积分时间:1s~10s可调;7)水平风速测量范围:0-75m/s;水平风向测量范围:0~360°;8)★风速精度:≤0.1m/s,风向精度:<3°;9)扫描模式:定点/DBS/VAD/PPI/RHI/CAPPI脚本编程;10)伺服扫描范围:水平方向:0~360°,垂直扫描:-90~+270°;11)伺服精度:0.1°;12)伺服扫描速度:0~55°/s 可调;13)数据存储:≥1T固态硬盘,可存储12个月以上数据;14)平均功率:<200W(常温);15)供电:220V/50Hz;16)★电气安全性:所投同型号产品满足低电压指令(LVD)标准;17)数据产品:DBS/VAD风廓线、垂直气流、RHI/PPI/VOL径向速度场及反演风场、虚拟测风塔、气溶胶后向散射系数、消光系数、气溶胶光学厚度(AOD)、边界层高度、云底高度、3D流场等;18)数据格式:ASCII/二进制;19)★扫描摄像系统:伺服扫描头需加装光学摄像头,与扫描头一体化安装并可与伺服扫描头随动,具备天气实况记录、安防功能,并能保存图片及影像;20)★整机重量<90kg,整机体积(L*W*H):≤800*800*1200mm,环境防护:防水防尘:IP66。
风廓线雷达
频谱非相参累积后频谱(左)与相参累积后(右)比较
应用现状
水平风测定:通过测量水平风廓线,研究暴 雨、锋面、急流、重力波等天气现象。1989 年中国气象科学研究院研制成功了我国首部 UHF 多普勒风廓线仪,并用于北京中尺度灾 害性天气预报基地的业务试验,证明了其有 效性和可靠性。王欣等用风廓线仪资料分析 安徽梅雨期间的强降水过程,表明其对水平 风的垂直结构有较强的探测能力,能实时监 测中尺度降水期间风的垂直切变和对流特征, 提高临近天气预报的精度。
风廓线雷达信息提取
相参累积:风廓线雷达所处理的信号一般为 弱信号 ,需要进行回波信号积累,一般即如 何从噪声之中获取所需的信息 。对单个脉冲 进行适当相参积累之后,提取风廓线信息就成 为一种可能。信号相参积累可以较好的提高 其SNR ,直观的看,只要进行足够数量的相参 积累就可以获得比较理想的SNR ,继而就可以 比较容易的获得所需信息。
探测湍流、边界层高度、空气污染,进行大 气环境研究。洪钟祥等用低层大气廓线仪 LAP3000等设备获取的遥感资料研究了湍流 热通量的垂直分布、地面热通量以及边界层 高度等。 ZHONG Shiyuan 等用布于 California 中部山谷的22 部风廓线系统,对 该区域盛行夏季风环流变化特征作了观测分 析,研究了该山谷三维风温结构的日变化特 征。
湿度廓线的推算:风廓线仪能够探测三维风 速场,无线电探声系统(RASS) 又能够获得有 效温度廓线。通过RASS 雷达获得的温度廓 线可以连续地估算出湿度廓线,从而能够对对 流层中水汽在大气运动中的作用进行更深入 的研究。
选址和维护
地理环境的选择:风廓线雷达的理想站点是 周围没有树木、高压输电线和建筑物的浅谷 或者开阔的平坦区域,且远离地面和空中交通, 并将电磁干扰降到最小。 电磁净空环境的选择:由于风廓线雷达为高 频发射和接收设备,信号容易受到附近电磁波 的影响,频率选定前应尽早到当地无线电管理 部门办理频率适用许可并办理无线电电台执 照
气象micaps-数据格式
气象micaps-数据格式数据格式说明第一类数据格式:用于地面填图diamond 1 屏幕上需显示的内容年月日时次总站点数区站号经度纬度拔海高度站点级别总云量风向风速海平面气压(本站气压)3小时变压过去天气1 过去天气2 6小时降水低云状低云量低云高露点能见度现在天气温度中云状高云状船向船速第二类数据格式:用于高空填图diamond 2 屏幕上需显示的内容年月日时次层次总站点数区站号经度纬度拔海高度站点级别高度温度温度露点差风向风速第三类数据格式:用于通用填图和离散点等值线diamond 3 屏幕上需显示的内容年月日时次,层次,线条数等值线值1 等值线值2 ......平滑系数,加粗线值,裁剪框的点数,经纬度值1,经纬度值2, ...... 单站内容长度总站点数区站号经度纬度拔海高度量值第四类数据格式:用于格点等值线diamond 4 屏幕上需显示的内容年月日时次时效层次经度格距纬度格距起始经度终止经度起始纬度终止纬度X-DIM Y-DIM 等值线间隔等值线起始值终止值平滑系数加粗线值第五类数据格式:用于TLOGP和剖面图diamond 5 屏幕上需显示的内容年月日时次总站点数区站号经度纬度跋海高度单站内容长度百帕高度温度露点风向风速百帕第六类数据格式:用于传真图1728 X 2400的点阵文件文件名按国际电码规定命名第七类数据格式:用于台风路径diamond 7 屏幕上需显示的内容台风名称台风编号发报中心总项数年月日时次时效中心经度中心纬度中心最低气压最大风速七级风圈半径十级风圈半径移向移速第八类数据格式:用于城市站点预报diamond 8 屏幕上需显示的内容年月日时次时效总站点数区站号经度纬度拔海高度天气现象风向风速最低温度最高温度天气现象风向风速第九类数据格式:用于底图投影变换(不可定义为综合图)diamond 9投影方式标准经度标准纬度 X放大系数 Y放大系数 X偏移 Y偏移预先保留预先保留预先保留本轮廓线点数标识字符串颜色线宽线型{X Y}···或者:diamond 90 0 0 0 0 0 0 0 0 0本轮廓线点数标识字符串颜色线宽线型{经度纬度}···四种投影的缺省标准配置如下:投影种类投影编码标准经度标准纬度 X放大系数 Y放大系数 X 偏移 Y偏移Lambert 1 79 29.99999 1.89999 1.89999-581 -3411.5999Mecator 2 无影响 0 0.558 0.573 478.0 1276.0bbq 3 20 无影响 0.25 0.25 0 0nbq 4 20 无影响 0.25 0.25 0 0系统定义的地图文件和参数文件:1.LAMBERT lambertM.dat para1.dat2.MECATOR mecatorM.dat para2.dat3.BBQ bbqMapXY.dat para3.dat4.NBQ nbqMapXY.dat para4.dat云图左下角经纬度(86.4,1.3)第十类数据格式:用于综合图定义(不可再次定义为综合图)diamond 10 综合图中所含的图类数数据文件路径可带统配符的文件名数据类型代码···第十一类数据格式:用于流线图diamond 11 屏幕上需显示的内容年月日时次时效层次经度格距纬度格距起始经度终止经度起始纬度终止纬度X-DIM Y-DIM第十二类数据格式:用于单点雷达图象(待定义)第十三类数据格式:用于图象(卫星云图、雷达拼图、地形图等) diamond (8个字符) 13(3个字符)屏幕上需显示的内容(40个字符)年(5个字符)月(3个字符)日(3个字符)时次(3个字符)X方向图象大小(5个字符)Y方向图象大小(5个字符)图象左下角经度坐标(8个字符)图象左下角纬度坐标(8个字符)投影方式(2个字符):1-lambert;2-mecator;3-bbq;4-nbq;放缩系数(5个字符)图象种类(2个字符):1-红外;2-雷达拼图;3-地形图;4-可见光;5-水汽对照表文件名(12个字符)系统保留(16个字符)(文件头部分共128个字符)下面是256级灰度表示值···第十四类数据格式:用于记录修改后的等值线diamond 14 屏幕上要显示的内容年月日时次时效LINES:条数线宽点数X Y Z···标号个数X Y Z······LINES_SYMBOL:条数编码线宽点数 NoLabel 0···SYMBOLS:个数编码 X Y Z 风向角度或字符串···CLOSED_CONTOURS:个数线宽点数X Y Z···标号个数X Y Z······STATION_SITUATION站号属性···第十五类数据格式:用于调色板设置diamond 15 屏幕上显示内容序号红色成分值绿色成分值蓝色成分值···系统色彩配置说明:0-1 清屏色和传真图颜色2-255 云图256-272 线条273-289 底图290-321 填图322-338 符号339 TLOGP站点色彩340 3、16、17类站圈颜色341 17类中的汉字底图: 274轮廓线颜色,275经度纬度线的颜色,279河流地面填图:风290 温度291 现在天气292 能见度293 露点294 总云量295 高云状296 中云状297低云状298 低云高299 气压300 3小时变压301 过去天气一302 过去天气二303 低云量304 6小时雨量305 站点306 高度307 离散点填图值308(-) 309(0) 310(+)高空填图:风311 温度312 高度313 露点314符号:雨雪322 冻雨、冰雹、沙暴323 风雾324 阴、晴、多云325 注326 L/G327 N/D328槽线329 暖锋330 冷锋331 固球锋332 高温区333 霜冻334 划线中间结果335划线确认后336 霜冻点337第十六类数据格式:用于确定预报区域的站点diamond 16 stations 总站点数区站号纬度经度级别···第十七类数据格式:用于站点信息diamond 17 STATION 总站点数区站号纬度(度分) 经度(度分) 高度级别单站内容长度 {单站内容} ···第十八类数据格式:用于数值预报剖面图diamond 18 屏幕上需显示的内容年月日时次时效等分数等值线间隔等值线起始值终止值平滑系数加粗线值显示层格点场(全路径文件名)层数格点场(全路径文件名)···第十九类数据格式:用于初始化参数文件(不可定义为综合图)diamond 19显示中心位置的经度显示中心位置的纬度放大倍数地图数据文件名要素设置缺省值要素设置缺省值定义如下:#define OBSERVATION_NONE 0X00000000#define OBSERVATION_WIND 0X00000001#define OBSERVATION_TEMPTURE 0X00000002#define OBSERVATION_HEIGHT 0X00000004 #define OBSERVATION_DEW_POINT 0X00000008#define OBSERVATION_CLOUD_VOLUME 0X00000010#define OBSERVATION_PRESSURE 0X00000020#define OBSERVATION_CHANGED_PRESSURE_3H 0X00000040#define OBSERVATION_LAST_WEATHER_1 0X00000080#define OBSERVATION_LAST_WEATHER_2 0X00000100#define OBSERVATION_RAIN_FALL_6H 0X00000200 #define OBSERVATION_LOW_CLOUD_FORM 0X00000400#define OBSERVATION_LOW_CLOUD_VOLUME 0X00000800 #define OBSERVATION_LOW_CLOUD_HEIGHT 0X00001000 #define OBSERVATION_TRANSPARENCY 0X00002000#define OBSERVATION_CURRENT_WEATHER 0X00004000#define OBSERVATION_MIDDLE_CLOUD_FORM 0X00008000 #define OBSERVATION_HIGH_CLOUD_FORM 0X00010000#define OBSERVATION_STATION_ID 0X00020000 #define OBSERVATION_SHIP_DIRECTION 0X00040000 #define OBSERVATION_SHIP_SPEED 0X00080000 #define OBSERVATION_ALL 0XFFFFFFFFMICAP主要功能1.资料显示l 地面填图l 高空填图l 通用填图l 格点等值线l 流线l 离散点等值线l 台风路径l 传真图l 卫星云图l 雷达资料l TlnP图l 站点空间剖面图l 格点空间剖面图l 站点信息2.图形操作l 放大l 缩小l 动画l 漫游l 隐现l 翻页3.图形编辑l 修改等值线l 修改站点预报4.预报制作l 城市预报产品制作l 区域预报产品制作。
气象 micaps-数据格式
数据格式说明第一类数据格式:用于地面填图diamond 1 屏幕上需显示的内容年月日时次总站点数区站号经度纬度拔海高度站点级别总云量风向风速海平面气压(本站气压) 3小时变压过去天气1 过去天气2 6小时降水低云状低云量低云高露点能见度现在天气温度中云状高云状船向船速第二类数据格式:用于高空填图diamond 2 屏幕上需显示的内容年月日时次层次总站点数区站号经度纬度拔海高度站点级别高度温度温度露点差风向风速第三类数据格式:用于通用填图和离散点等值线diamond 3 屏幕上需显示的内容年月日时次,层次,线条数等值线值1 等值线值2 ......平滑系数,加粗线值,裁剪框的点数,经纬度值1,经纬度值2, ...... 单站内容长度总站点数区站号经度纬度拔海高度量值第四类数据格式:用于格点等值线diamond 4 屏幕上需显示的内容年月日时次时效层次经度格距纬度格距起始经度终止经度起始纬度终止纬度X-DIM Y-DIM 等值线间隔等值线起始值终止值平滑系数加粗线值第五类数据格式:用于TLOGP和剖面图diamond 5 屏幕上需显示的内容年月日时次总站点数区站号经度纬度跋海高度单站内容长度百帕高度温度露点风向风速百帕第六类数据格式:用于传真图1728 X 2400的点阵文件文件名按国际电码规定命名第七类数据格式:用于台风路径diamond 7 屏幕上需显示的内容台风名称台风编号发报中心总项数年月日时次时效中心经度中心纬度中心最低气压最大风速七级风圈半径十级风圈半径移向移速第八类数据格式:用于城市站点预报diamond 8 屏幕上需显示的内容年月日时次时效总站点数区站号经度纬度拔海高度天气现象风向风速最低温度最高温度天气现象风向风速第九类数据格式:用于底图投影变换(不可定义为综合图)diamond 9投影方式标准经度标准纬度 X放大系数 Y放大系数 X偏移 Y偏移预先保留预先保留预先保留本轮廓线点数标识字符串颜色线宽线型{X Y}···或者:diamond 90 0 0 0 0 0 0 0 0 0本轮廓线点数标识字符串颜色线宽线型{经度纬度}···四种投影的缺省标准配置如下:投影种类投影编码标准经度标准纬度 X放大系数 Y放大系数 X 偏移 Y偏移Lambert 1 79 29.99999 1.89999 1.89999-581 -3411.5999Mecator 2 无影响 0 0.558 0.573 478.0 1276.0bbq 3 20 无影响 0.25 0.25 0 0nbq 4 20 无影响 0.25 0.25 0 0系统定义的地图文件和参数文件:1.LAMBERT lambertM.dat para1.dat2.MECATOR mecatorM.dat para2.dat3.BBQ bbqMapXY.dat para3.dat4.NBQ nbqMapXY.dat para4.dat云图左下角经纬度(86.4,1.3)第十类数据格式:用于综合图定义(不可再次定义为综合图)diamond 10 综合图中所含的图类数数据文件路径可带统配符的文件名数据类型代码···第十一类数据格式:用于流线图diamond 11 屏幕上需显示的内容年月日时次时效层次经度格距纬度格距起始经度终止经度起始纬度终止纬度X-DIM Y-DIM第十二类数据格式:用于单点雷达图象(待定义)第十三类数据格式:用于图象(卫星云图、雷达拼图、地形图等) diamond (8个字符) 13(3个字符)屏幕上需显示的内容(40个字符)年(5个字符)月(3个字符)日(3个字符)时次(3个字符)X方向图象大小(5个字符)Y方向图象大小(5个字符)图象左下角经度坐标(8个字符)图象左下角纬度坐标(8个字符)投影方式(2个字符):1-lambert;2-mecator;3-bbq;4-nbq;放缩系数(5个字符)图象种类(2个字符):1-红外;2-雷达拼图;3-地形图;4-可见光;5-水汽对照表文件名(12个字符)系统保留(16个字符)(文件头部分共128个字符)下面是256级灰度表示值···第十四类数据格式:用于记录修改后的等值线diamond 14 屏幕上要显示的内容年月日时次时效LINES:条数线宽点数 X Y Z···标号个数 X Y Z······LINES_SYMBOL:条数编码线宽点数 NoLabel 0···SYMBOLS:个数编码 X Y Z 风向角度或字符串···CLOSED_CONTOURS:个数线宽点数 X Y Z···标号个数 X Y Z······STATION_SITUATION站号属性···第十五类数据格式:用于调色板设置diamond 15 屏幕上显示内容序号红色成分值绿色成分值蓝色成分值···系统色彩配置说明:0-1 清屏色和传真图颜色2-255 云图256-272 线条273-289 底图290-321 填图322-338 符号339 TLOGP站点色彩340 3、16、17类站圈颜色341 17类中的汉字底图: 274轮廓线颜色,275经度纬度线的颜色,279河流地面填图:风290 温度291 现在天气292 能见度293 露点294 总云量295 高云状296 中云状297低云状298 低云高299 气压300 3小时变压301 过去天气一302 过去天气二303 低云量304 6小时雨量305 站点306 高度307 离散点填图值308(-) 309(0) 310(+)高空填图:风311 温度312 高度313 露点314符号:雨雪322 冻雨、冰雹、沙暴323 风雾324 阴、晴、多云325 注326 L/G327 N/D328槽线329 暖锋330 冷锋331 固球锋332 高温区333 霜冻334 划线中间结果335划线确认后336 霜冻点337第十六类数据格式:用于确定预报区域的站点diamond 16 stations 总站点数区站号纬度经度级别···第十七类数据格式:用于站点信息diamond 17 STATION 总站点数区站号纬度(度分) 经度(度分) 高度级别单站内容长度 {单站内容}···第十八类数据格式:用于数值预报剖面图diamond 18 屏幕上需显示的内容年月日时次时效等分数等值线间隔等值线起始值终止值平滑系数加粗线值显示层格点场(全路径文件名)层数格点场(全路径文件名)···第十九类数据格式:用于初始化参数文件(不可定义为综合图)diamond 19显示中心位置的经度显示中心位置的纬度放大倍数地图数据文件名要素设置缺省值要素设置缺省值定义如下:#define OBSERVATION_NONE 0X00000000#define OBSERVATION_WIND 0X00000001#define OBSERVATION_TEMPTURE 0X00000002#define OBSERVATION_HEIGHT 0X00000004 #define OBSERVATION_DEW_POINT 0X00000008#define OBSERVATION_CLOUD_VOLUME 0X00000010#define OBSERVATION_PRESSURE 0X00000020#define OBSERVATION_CHANGED_PRESSURE_3H 0X00000040#define OBSERVATION_LAST_WEATHER_1 0X00000080#define OBSERVATION_LAST_WEATHER_2 0X00000100#define OBSERVATION_RAIN_FALL_6H 0X00000200 #define OBSERVATION_LOW_CLOUD_FORM 0X00000400#define OBSERVATION_LOW_CLOUD_VOLUME 0X00000800#define OBSERVATION_LOW_CLOUD_HEIGHT 0X00001000#define OBSERVATION_TRANSPARENCY 0X00002000#define OBSERVATION_CURRENT_WEATHER 0X00004000#define OBSERVATION_MIDDLE_CLOUD_FORM 0X00008000#define OBSERVATION_HIGH_CLOUD_FORM 0X00010000#define OBSERVATION_STATION_ID 0X00020000 #define OBSERVATION_SHIP_DIRECTION 0X00040000 #define OBSERVATION_SHIP_SPEED 0X00080000 #define OBSERVATION_ALL 0XFFFFFFFFMICAP主要功能1.资料显示l 地面填图l 高空填图l 通用填图l 格点等值线l 流线l 离散点等值线l 台风路径l 传真图l 卫星云图l 雷达资料l TlnP图l 站点空间剖面图l 格点空间剖面图l 站点信息2.图形操作l 放大l 缩小l 动画l 漫游l 隐现l 翻页3.图形编辑l 修改等值线l 修改站点预报4.预报制作l 城市预报产品制作l 区域预报产品制作。
《风廓线雷达单站数据 NetCDF格式》编制说明
气象行业标准《风廓线雷达单站数据 NetCDF格式》编制说明一、工作简况1.任务来源本标准由中国气象局提出,全国气象基本信息标准化技术委员会(SAC/TC 346)归口。
2019年由中国气象局下达中国气象局气象探测中心,项目编号QX/T-2019-91,计划项目名称为《风廓线雷达单站数据NetCDF格式》。
2.起草单位本标准起草单位为:中国气象局气象探测中心。
3.标准主要起草人及其工作分工本标准主要起草人为李瑞义、吴蕾、杨馨蕊、周薇、赵世颖、汪学渊、董德宝。
其分工如下:李瑞义,项目负责人,负责标准总体设计,标准编写审定,编制说明的编写审定;吴蕾、汪学渊、董德宝负责调研,标准制定;杨馨蕊、周薇、赵世颖,负责资料收集整理。
4.主要工作过程(1)、2017年9月,应气象信息化标准建设要求,开始编制风廓线雷达标准数据格式。
在官方收集、调研国外相关标准、规范等资料基础上,编制组完成《风廓线雷达标准数据格式说明》讨论稿(第一稿)。
(2)、2017年12月,根据中国气象局预报与网络司、中国气象局综合观测司、国家气象信息中心、国家气象中心、中国气象局公共气象服务中心、国家气候中心多部门专家的意见,形成讨论稿(第二稿),标准名称改为《风廓线雷达台站产品数据格式说明(NETCDF)》。
(3)、2018年3月,再次征求中国气象局预报与网络司、中国气象局综合观测司、国家气象信息中心、国家气象中心、中国气象局公共气象服务中心、国家气候中心多部门专家的意见,形成讨论稿(第三稿),标准名称改为《风廓线雷达单站数据NETCDF格式》。
(4)、2018年6月,编写标准格式转换软件,对NETCDF中组、维、变量和属性进行了微调,形成《风廓线雷达单站数据NETCDF格式》(第四稿)。
(5)、2018年10月起-12月底,标准格式转换软件在北京大兴南郊观象台试点。
(6)、 2019年7月起-10月底,优化的标准格式转换软件在北京七个风廓线雷达站试点。
风廓线雷达通用数据格式.doc
风廓线雷达通用数据格式(V1.2)2007年9月目录1. 文件名编码规则 (3)1.1 原始数据文件 (3)1.2 产品数据文件 (3)2.功率谱数据文件 (4)3. 径向数据文件 (4)3.1 文件组成单位 (4)3.2 文件框架 (4)3.3 文件结构 (6)4. 实时的采样高度上的产品数据文件 (10)4.1 文件组成单位 (10)4.2 文件框架 (10)4.3 文件结构 (10)5. 半小时平均的采样高度上的产品数据文件 (12)5.1 文件组成单位 (12)5.2 文件框架 (12)5.3 文件结构 (12)6. 一小时平均的采样高度数据文件 (13)6.1 文件组成单位 (13)6.2 文件框架 (13)6.3 文件结构 (13)附件一功率谱数据格式 (15)1. 文件名编码规则根据实际需求,建议使用长文件名命名法,对各类文件名进行约定。
文件名中的观测时间均为观测结束时间。
1.1 原始数据文件原始数据文件包括功率谱数据文件、瞬时径向谱数据文件,对于原始数据文件,建议每次观测生成一个文件,文件名具体命名方法如下:Z_RADR_I _IIiii_yyyyMMddhhmmss_O_WPRD_雷达型号_数据类型.TTT 其中:Z:国内交换文件;RADR:表示雷达资料;I:表示后面的IIiii为风廓线雷达站的区站号;IIiii:区站号(按地面气象站的区站号);yyyy:观测时间(年) (20**—);MM:观测时间(月) (01—12);dd:观测时间(日) (01—31);hh:观测时间(时) (00—23);mm:观测时间(分) (00—59);ss:观测时间(秒) (00—59);O:表示观测数据;WPRD:表示风廓线雷达资料;雷达型号:见表1;数据类型:功率谱数据文件用FFT表示;径向数据文件用RAD表示;TTT:当TTT = BIN时,表示二进制文件;当TTT = TXT时,表示文件格式为ASCII。
MICAPS数据格式(全部)
附录一MICAPS数据文件格式MICAPS系统得数据结构就是建立在文件系统基础上得。
其特点就是:利用目录来区分不同得数据来源、要素与层次,即不同得数据来源、要素与层次得数据要放在不同得目录中。
同一目录中得数据只能有时次或时效上得不同。
系统根据不同得数据格式来显示不同类型得图象。
除第六类数据(传真图)外,每个数据文件都有一个文件头,描述该数据文件属于哪一类数据格式、数据得日期、时次、时效及其它有关参数.除第6与13类数据(图象)外,数据文件均为文本文件。
MICAPS3、2定义或其它常用气象数据格式有:闪电定位数据(MICAPS3、2定义得第41类数据)GPS水汽数据(MICAPS3、2定义得第42类数据)地图信息数据(第9类数据扩展格式)自动站Z文件风廓线Z文件AWX格式卫星云图及产品HDF格式卫星云图标称图产品GPF格式卫星云图数据雷达基数据雷达PUP产品netCDF数据另外MICAPS3、2使用得通用数据格式:MIF格式地理信息数据、SHP格式地理信息数据。
A1、1第一类数据格式:地面全要素填图数据文件头:diamond 1 数据说明(字符串)年月日时次总站点数(均为整数)注:此类数据用于规范得地面填图数据:区站号(长整数) 经度纬度拔海高度(均为浮点数)站点级别(整数) 总云量风向风速海平面气压(或本站气压) 3小时变压过去天气1 过去天气2 6小时降水低云状低云量低云高露点能见度现在天气温度中云状高云状标志1 标志2(均为整数) 24小时变温24小时变压注:缺值时用9999表示,以后相同.站点级别表示站点得放大级别,即只有当图象放大到该级别时此站才被填图.以后相同。
当标志1为1,标志2为2时,说明后面有24小时变温变压.否则说明后面没有24小时变温变压。
例子:diamond 1 99年06月15日08时地面填图99 06 15 08 301650468127、45 50、25 166 16 7 340 6 975 48 0、1 38 7 600 9、1 25、0 0 14、7 9999 999912 1 -352533 98、48 39、77 1478 1 8 0 0 98 78 0、0130 8 2500 10、7 30、0 60 16、8 27 999912 2 352652 100、43 38、93 1483 4 8 270 3115 11 6 0、530 4 2500 12、6 15、061 16、0 24 17 12 1 2注意:地面自动站数据如果写为该类格式,可以在文件说明中加入可以识别得文字,默认使用“自动”作为识别文字,但自动站文件名定义可以使用8、3格式得“年月日时。
民用机场风廓线雷达系统技术规范
民用机场风廓线雷达系统技术规范管理程序中国民用航空局空管行业管理办公室编号:AP-117-TM-2013-01下发日期:2013年8月6日民用机场风温廓线雷达系统技术规范目录第一章总则 (2)第二章系统构成及功能 (2)第一节一般规定 (2)第二节产品输出功能 (4)第三章系统性能 (6)第一节整体性能 (6)第二节各子系统性能 (7)第四章环境适应性 (12)附录一信号功率的谱矩及信噪比计算方法 (14)附录二风速、风向及Cn2计算方法 (15)附录三风温廓线雷达数据格式 (16)民用机场风温廓线雷达系统技术规范第一章总则第一条为规范民用机场风温廓线雷达系统(以下简称风温廓线雷达)的建设和运行,根据《中国民用航空气象工作规则》,制订本规范。
第二条本规范适用于中华人民共和国境内民用机场和军民合用机场民用部分(以下称民用机场)的风温廓线雷达系统的建设和运行。
第三条民用机场风温廓线雷达系统的构成、功能、性能和环境适应性等技术要求应当符合本规范。
第四条风温廓线雷达按照安装方式不同分为固定式和可移式两种。
可移式主要有车载可移式和方舱可移式两种。
第二章系统构成及功能第一节一般规定第五条风温廓线雷达主要由天线分系统、发射分系统、接收分系统、信号处理分系统、监控分系统、标定分系统、通讯分系统、数据处理及应用终端、配电分系统、RASS(Radio-Acoustic Sounding System,无线电-声探测系统)分系统构成。
第六条风温廓线雷达采用全相参脉冲多普勒体制。
第七条风温廓线雷达采用相控阵技术、全固态发射技术、数字接收机技术、脉冲压缩技术等。
第八条风温廓线雷达具有探测大气虚温的功能。
第九条风温廓线雷达应当具有实时采集功率谱数据,对功率谱数据进行信噪比计算、噪声电平计算、干扰信号剔除、杂波剔除、信号提取、谱矩估计、数据质量控制功能。
第十条风温廓线雷达应当具有生成风速、风向、垂直气流、Cn2(Refractive Index Structure Constant,大气折射率结构常数)、谱宽、信噪比等数据产品的功能。
L波段边界层风廓线雷达(含RASS)
适应温度:室外为-40~50℃,室内为 0~30℃; 适应湿度:室外为≤100%,室内为≤95%; 抗风能力:能经受不大于 50m/s 的阵风风速(或不 大于 30m/s 的平稳风速),在上述情况下,天线和电磁 屏蔽网不产生永久性变形或破坏; 其它环境适应性:防盐雾、防霉、防沙尘和防雷击 能力;绝缘电阻≥10Ω,接地电阻≤4Ω; 电磁兼容性:具有静电屏蔽、电磁屏蔽设计,设备 地线中模拟地线(Ga)与数字地线(Gd)和安全地线(Gp) 要严格分开,以增强设备的抗干扰能力。 五波束,一个铅垂方向波束和四个方位相互正交、 具有相同仰角的倾斜波束。
支持远程监控及数据传输
1.具备 24 小时不间断探测的能力; 2.具备雷达数据自动处理功能: 根据工作模式和配置参数,生成 0~2 级雷达数据产 品;其中 0 级雷达数据产品是设备直接输出的数据,不 经过任何处理; 3.具备雷达数据产品综合显示功能: 可按照操作员的要求,实现 0~2 级雷达数据产品的 图形化展示,图形保存与输出等功能; 4.具备探测资料存储和网络传输能力: 实现雷达数据的分类存储、备份,自动/手动清理数 据缓存; 5.具备无人值守自动观测功能: 可根据工作模式和配置文件,实现长期连续自动观 测,外电中断时自动保存数据并关机,外电恢复后自动 启动并继续观测; 6.响应文件应包含以下方案: (1)技术方案(电讯总体和结构总体); (2)测试大纲和细则; (3)出厂验收细则; (4)系统现场验收细则。
倾斜波束倾角
15 5°
★天线分系统(包 括但不限于) ★接收分系统(包括 但不限于)
方舱要求
基建要求
天线增益: ≥30dB 波束宽度: ≤4.5° 动态范围: ≥90dB 噪声系数: ≤2 dB 2.5 米×3 米及以上。根据采购人要求,对方舱进行 相关补充设计,实现采购人在联合观测方面(如采购人 自己增加的探测设备)的功能需求。
风廓线雷达应用
改进数值预报初值场
(a)6月09日00时925百帕实测风,流线和位势高度客观分析 图,
6月09日00时925百帕的风场形势 (a) 控制实验,(b)风廓线试验,
风廓线雷达应用
1 风廓线多普勒相控阵雷达探测原理 2 风廓线的基本分析方法 3 应用实例
风廓线应用实例(锋面)
1990年8月7日风廓线图(0~3km) b、1990年8月7日11时、12时30分温度平流廓线图(0~3km)
低涡暖切变暴雨
1989年7月21日风廓线 图
1989年7月23日风廓线
强对流暴雨
上海风廓线雷达测得的水平风场:a 、b分别为99年9月6日用风矢量表示的水 平风随时间-高度分布图,横坐标表示时间,纵坐标为高度,风矢量上的颜 色等级表示垂直速度大小(向上为正)。
是水平 风, 不是垂 直环流 !!
时间: 从右 向左
锋面和切变线的监测
• 1989年8月9日~10日风廓线图(0~10km)
锋面和高空槽过境
• 1990年8月1日~3日风廓线图(0~10km)
低空急流和行星边界层风场
• 1989年8月9日~10日风廓线图(0~3km)
天气系统 短时变化 (中尺度雨
风廓线雷达应用
1 风廓线多普勒相控阵雷达探测原理 2 风廓线的基本分析方法 3 应用实例
Height(m) 1800
600
150 15:00
ห้องสมุดไป่ตู้
14:00 13:00
Height(km)
6 5.5
5 4.5
4 3.5
3 2.5
2 1.5
风廓线雷达观测规定(试行)
查看并保存运行状态日志
有()无()
7
前一日风羽图探测高度(km)
8
发射机峰值功率(kw)
9
发射机温度(℃)
10
系统有无故障
有()无()
11
故障报警时间(有故障时填写)
12
故障报警信息(有故障时填写)
附件2:
风廓线雷达每月标定记录表
站名:值班员:年 月 日
标定内容
标定结果
系统灵敏度
动态范围
系统相干性
第三章观测模式
第十一条风廓线雷达一般采用五波束观测,当受观测环境限制时,经中国气象局业务主管部门同意后,也可采用三波束观测。
第十二条风廓线雷达观测模式由若干高度上的观测模式组合而成,在设备安装调试完成后即确定组合。
第十三条风廓线雷达观测模式一般不得更改,如遇特殊情况需要更改,应当事先征得中国气象局业务主管部门的同意。
(试行)》的通知
各省、自治区、直辖市气象局:
《风廓线雷达观测规定(试行)》已经中国气象局批准,现予印发,请遵照执行。
二〇一一年十一月七日
风廓线雷达观测规定(试行)
第一章总则
第一条为加强对风廓线雷达观测业务的管理,根据《中华人民共和国气象法》,考虑到风廓线雷达的功能及特点,制定本规定。
第二条本规定适用于纳入中国气象局统一布局的用于气象业务的固定式边界层风廓线雷达和对流层风廓线雷达。移动风廓线雷达参照本规定执行。
第三十条本规定自发布之日起施行。
附件1:
风廓线雷达值班日记
站名:值班员:年 月 日
序号
值班时查看内容
检查情况
1
机房和天线阵地状况是否正常
正常()异常()
2
机房电源、UPS电源输出是否正常
《气象探测环境保护规范 风廓线雷达站》编制说明
《气象探测环境保护规范风廓线雷达站》编写说明一、工作简况1、任务来源本标准由全国气象仪器与观测方法标委会(SAC/TC 507)提出并归口。
2016年2月29日由中国气象局发布《中国气象局政策法规司关于下达2019 年气象观测装备相关标准制定计划的通知》(气法函〔2018〕62 号),项目编号:QX/T-2019-38。
2、编制单位本标准编制单位为中国气象局气象探测中心。
3、主要起草人及所做的工作4、主要工作过程1)2018年9月,编制单位召开本标准编制启动会,成立编制组,商讨工作计划;2)2018年12月,收集整理国内外有关风廓线雷达安装、运行等环境要求,综合考虑,保障风廓线雷达对环境的基本要求,包括场地、遮蔽等情况;3)2019年1月,开始本标准的编写;4)2019年12月,基本完成标准的编制工作,形成征求意见稿并在全国范围征求意见;二、标准的编制原则和主要内容1、标准的编制原则本标准在编制方面,坚持了以下几项原则:1) 科学性原则借鉴、参考和吸收国际先进经验和做法,分析、比较不同型号、频段雷达的情况,结合我国实际情况进行分析,使提出的《气象探测环境保护规范风廓线雷达站》标准具有科学性。
2) 适应性原则早前,中国气象局发布了中气函〔2012〕185号《风廓线雷达功能设计规范(L波段)》和《风廓线仪功能规格需求书》等一些指导性文件,用于指导现阶段风廓线雷达的设计与生产。
《气象探测环境保护规范风廓线雷达站》标准的编制,其总体性能指标与指导文件基本保持一致的前提下,对观测站点周边环境保护提出了更为具体的要求,保证风廓线雷达的探测环境。
3) 规范性原则按照GB/T 1.1-2009的要求进行编写。
2、标准主要内容本标准规定了风廓线雷达站探测环境保护的范围和要求。
本标准首先对风廓线雷达站探测环境保护涉及的一系列专有名词进行解释。
例如风廓线雷达、风廓线雷达站、障碍物、电磁干扰、斜波束、波束宽度等指标。
本标准其次对地理环境和电磁环境提出保护范围和要求。
SA SB CB雷达数据格式说明
CINRAD SA/SB 雷达基数据格式说明:1.数据的存储方式每个体扫存储为一个单独的文件2.数据的排列方式按照径向数据的方式顺序排列,对于CINRAD SA/SB雷达,体扫数据排列自低仰角开始到高仰角结束。
3.径向数据的长度径向数据的长度固定,为2432字节。
4.距离库长和库数反射率距离库长为1000米,最大距离库数为460;速度和谱宽距离库长为250米,最大距离库数为920。
CINRAD CB雷达基数据格式说明:5.数据的存储方式每个体扫存储为一个单独的文件6.数据的排列方式按照径向数据的方式顺序排列,对于CINRAD CB雷达,体扫数据排列自低仰角开始到高仰角结束。
7.径向数据的长度径向数据的长度固定,为4132字节。
8.距离库长和库数反射率距离库长为500米,最大距离库数为800;速度和谱宽距离库长为125米,最大距离库数为1600。
程序中的重要数据说明1.文件名Filename[],输入需要读取的基数据的文件名。
需将该文件放在执行程序所在的目录中才能读出其中的数据。
2.保存反射率、速度、谱宽,各层仰角的数组。
文件中读取的基数据存放在下列数组中:float V olRef[MaxCuts][MaxRads][RGates]; //反射率(浮点型,单位:DBZ)float V olVel[MaxCuts][MaxRads][VGates]; //速度(浮点型,单位:M/S)float V olSpw[MaxCuts][MaxRads][WGates]; //谱宽(浮点型,单位:M/S)float Elvation[MaxCuts]; //各层仰角(浮点型,单位:度)数组中无效数据标记为-999.0,距离折叠标记为999.0。
其中,1)MaxCuts=20,为最大层数;2)MaxRads为方位数,每度保存一个径向;3)Rgates为每个径向上反射率的距离库数,C波段为800,对应分辨率为0.5公里;S 波段为460,对应分辨率为1公里;4)Vgates为每个径向上径向速度的距离库数,C波段为1600,对应分辨率为0.125公里;S波段为920,对应分辨率为0.25公里;5)Wgates为每个径向上谱宽的距离库数,C波段为1600,对应分辨率为0.125公里;S波段为920,对应分辨率为0.25公里;3.读取不同波段的基数据文件的方法在头文件DataFormat.h中,对距离库数的定义为,用来读取S波段的基数据:const int RGates = 460; //反射率距离库数const int VGates = 920; //速度距离库数const int WGates = 920; //谱宽距离库数若要读取C波段的基数据时,只需将上述定义修改为:const int RGates = 800; //反射率距离库数const int VGates = 1600; //速度距离库数const int WGates = 1600; //谱宽距离库数注意:1)关于仰角层的说明:SA,SB,CB雷达在低层每个仰角上扫描两次,程序中,在保存基数据到数组中时,记为一个仰角层。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
其中:
Z:
国内交换文件;
RADR:
表示雷达资料;
I:
表示后面的 IIiii 为风廓线雷达站的区站号;
IIiii:
区站号(按地面气象站的区站号);
yyyy:
观测时间(年) (20**—);
MM:
观测时间(月) (01—12);
dd:
(00—23);
四位小数
测站的纬度,以度为单位,其中第一位为符
3 纬度
8 字节 号位,北纬取正,南纬取负,两位整数,四
位小数
4 观测场拔海高度
7 字节
观测场拔海高度,以米为单位,其中第一位 为符号位,四位整数,一位小数
5 风廓线雷达型号
2 字节 风廓线雷达型号,具体标识见表 2
6 回车换行
2 字节
▲ 第 3 段为低模式雷达性能参数,本段每个采集站点有且仅有一条记录,
名方法如下:
Z_RADR_I_IIiii_yyyyMMddhhmmss_P_WPRD_雷达型号_产品标识.TXT
其中:
Z:
国内交换文件;
RADR:
表示雷达资料;
I:
表示后面的 IIiii 为风廓线雷达站的区站号;
IIiii:
区站号(按地面气象站的区站号);
yyyy:
观测时间(年) (20**—);
MM:
附件:
风廓线雷达通用数据格式 (V1.2)
2007 年 9 月
目录
1. 文件名编码规则...................................................3 1.1 原始数据文件 ............................................................................................... 3 1.2 产品数据文件 ............................................................................................... 3
注:观测时间用世界时表示。
表 2 风廓线雷达产品标识
产品 实时的采样高度上的产品数据文件
半小时平均的采样高度上的产品数据文件 一小时平均的采样高度上的产品数据文件
产品标识 ROBS
HOBS OOBS
2.功率谱数据文件
功率谱数据文件由文件标识、测站基本参数、性能参数、观测参数及观测数 据组成,全部为二进制格式,功率谱数据文件根据需求实时动态生成。
注:观测时间用世界时表示。
表 1 风廓线雷达型号标识符
雷达种类
说明
标识符
风廓线雷达
P 波段,对流层 I 型风廓线雷达
PA
P 波段,对流层 II 型风廓线雷达
PB
L 波段,边界层风廓线雷达
LC
1.2 产品数据文件
产品数据文件包括实时的采样高度上的产品数据文件、半小时平均的采样高
度上的产品数据文件,一小时平均的采样高度上的产品数据文件,文件名具体命
波束 5 观测数据 NNNN RAD SIXTH 波束 6 观测数据 NNNN 中模式雷达性能参数 中模式观测参数 RAD FIRST 波束 1 观测数据 NNNN RAD SENCOND 波束 2 观测数据 NNNN RAD THIRD 波束 3 观测数据 NNNN RAD FOURTH 波束 4 观测数据 NNNN RAD FIFTH 波束 5 观测数据 NNNN RAD SIXTH 波束 6 观测数据 NNNN 高模式雷达性能参数 高模式观测参数 RAD FIRST 波束 1 观测数据 NNNN RAD SENCOND 波束 2 观测数据 NNNN RAD THIRD 波束 3 观测数据 NNNN RAD FOURTH 波束 4 观测数据 NNNN RAD FIFTH 波束 5 观测数据 NNNN RAD SIXTH 波束 6 观测数据
mm:
观测时间(分) (00—59);
ss:
观测时间(秒) (00—59);
O:
表示观测数据;
WPRD:
表示风廓线雷达资料;
雷达型号: 见表 1;
数据类型: 功率谱数据文件用 FFT 表示;
径向数据文件用 RAD 表示;
TTT:
当 TTT = BIN 时,表示二进制文件;
当 TTT = TXT 时,表示文件格式为 ASCII。
2 字节
▲ 第 2 段为测站基本参数,本段每个采集站点有且仅有一条记录,记录内
容参见表 4。
表 4 第 2 段记录格式说明表
序号
各组含义
额定长度
说明
1 区站号
5 字节
五位数字或第一位为字母,第二-五位为数 字
测站的经度,以度为单位,其中第一位为符
2 经度
9 字节 号位,东经取正,西经取负,三位整数,
观测时间(月) (01—12);
dd:
观测时间(日) (01—31);
hh: mm:
观测时间(时) 观测时间(分)
(00—23); (00—59);
ss:
观测时间(秒) (00—59);
P:
表示产品数据;
WPRD:
表示风廓线雷达资料;
雷达型号: 见表 1;
产品标识: 见表 2;
TXT:
表示文件格式为 ASCII。
2.功率谱数据文件................................................... 4 3. 径向数据文件.....................................................4
3.1 文件组成单位 ............................................................................................... 4 3.2 文件框架 ....................................................................................................... 4 3.3 文件结构 ....................................................................................................... 6 4. 实时的采样高度上的产品数据文件..................................10 4.1 文件组成单位 ............................................................................................. 10 4.2 文件框架 ..................................................................................................... 10 4.3 文件结构 ..................................................................................................... 10 5. 半小时平均的采样高度上的产品数据文件............................12 5.1 文件组成单位 ............................................................................................. 12 5.2 文件框架 .....................................................................................................12 5.3 文件结构 .....................................................................................................12 6. 一小时平均的采样高度数据文件....................................13 6.1 文件组成单位 ............................................................................................. 13 6.2 文件框架 ..................................................................................................... 13 6.3 文件结构 ..................................................................................................... 13 附件一 功率谱数据格式..............................................15
一位小数
5 南波束与铅垂线的夹角 4 字节 南波束与铅垂线的夹角(度),两位整数,
一位小数
6 北波束与铅垂线的夹角 4 字节 北波束与铅垂线的夹角(度),两位整数,
一位小数
7 中(行)波束与铅垂线 4 字节 中(行)波束与铅垂线的夹角(度),两位
的夹角(度)
整数,一位小数
8 中(列)波束与铅垂线 4 字节 中(列)波束与铅垂线的夹角(度),两位
格式说明见附件一。
3. 径向数据文件
3.1 文件组成单位 一次探测形成一个文件。