CINRAD-SA双偏振雷达资料在降水估测中的应用初探

对基于水平反射率ZH和差分传播相移率K_(DP)的降水估测综合法R(C)进行了改进,并对广州S波段双偏振雷达2016年2次飑线和2次台风降水过程的Φ_(DP)使用小波分析进行滤波处理,在此基础上使用变距最小二乘法拟合得到K_(DP)的值。分别使用R(C)和R(Z_H)法对2次飑线和2次台风降水过程进行降水估算,将估算结果和雨量计小时雨量进行了对比,并将两种方法的评估结果进行了对比。结果表明:(1)对于飑线类型降水,R(C)法对5 mm·h~(-1)以上的降水估测精度要好于R(Z_H)法,且降水率越大,R(C)法优势越明显,当降水率≥20 mm·h~(-1)时,两次过程R(C)法比R(Z_H)法的平均相对误差(RE)降低了17. 2%,平均绝对误差(AE)减少了1.89 mm,平均均方根误差(RMSE)减少了1.66 mm;(2)对于台风类型降水,R(C)法对5 mm·h~(-1)以上的降水估测精度也好于R(Z_H)法,当降水率≥20 mm·h~(-1)时,两次过程R(C)法比R(Z_H)法的平均RE降低了33. 19%,平均AE减少了3. 95 mm,平均RMSE减少了4.05 mm;(3)对于飑线和台风两种类型降水R(C)法都明显改善了降水率较大时的R(Z_H)法低估问题,但R(C)法在降水率10 mm·h~(-1)时也存在低估,可能是由雨滴谱资料观测误差导致拟合的系数偏小或雷达硬件造成的观测偏差等造成的。     

飞机积冰预报算法对比及其集成预报模型研究

以人工增雨作业获取的飞机积冰实例资料为基础,利用WRF模式对51次飞机积冰过程进行数值模拟,对比分析了常用七种积冰预报算法对积冰潜势区和强度的预报效果,进而采用评分权重集成法建立了飞机积冰强度集成预报模型,并检验了其预报效果。结果表明:(1)假霜点温度经验法对2002年4月4日积冰个例的预报效果与实况一致,而其他积冰算法预报效果均与实况相差较大;(2)对51次飞机积冰预报效果进行统计检验发现,假霜点温度经验法的预报效果最好,积冰强度预报准确率为72.55%,其次是RAOB法,IC指数法和I积冰指数法次之,改进的IC指数法预报准确率最差,只有19.61%;(3)对比不同积冰算法建立的集成预报模型的预报效果发现,选用IC指数法、假霜点温度经验法、RAOB法进行集成预报时,预报准确率最高,且漏报率、偏弱率及偏强率均能控制在10%以内,比单一预报算法中的最高预报准确率提高了8%,且漏报率降低了4%,偏强率降低了8%。     

高分辨率模式雷达回波预报能力分析

利用2018年7—8月GRAPES_3 km、东北短临(WRFRUC)高分辨率模式综合雷达回波预报数据和辽宁省SWAN雷达组合反射率(MCR)实况,基于邻域法FSS评分指数,分析模式在台风北上和副热带高压边缘暴雨过程中的雷达回波预报能力。结果表明:两家模式在不同降水过程中对小阈值雷达回波有较好的预报技巧,随着回波量级增大,模式预报FSS逐渐减小,雷达回波55 dBz时,FSS甚至为0。当邻域半径是3时,35 dBz以下的回波预报中GRAPES模式在台风北上暴雨中的预报技巧低于副热带高压边缘,35 dBz则相反。WRFRUC模式始终表现为台风北上暴雨中预报较好。当邻域半径9时,WRFRUC模式在台风暴雨中的FSS评分高于GRAPES模式,GRAPES模式在副热带高压暴雨中的FSS评分始终高于WRFRUC模式。GRAPES和WRFRUC模式的最大FSS评分技巧均出现在邻域半径是11时,分别为0.239和0.195。GRAPES模式中FSS评分在12 h逐小时预报中前3个时次较强,WRFRUC模式则表现为中间时次强,两头弱。     

山东滨州市冰雹天气分型和预报方法研究

利用山东滨州市7个国家气象站与61个人影作业站点观测资料,结合高空观测及探空、ECMWF再分析等资料,对2001—2011年滨州降雹时空分布特征、天气系统和物理量特征、降雹形势分型和预报方法进行研究。得出:(1)降雹日数年均8.6次,总体呈现明显下降的年际变化特征;4—10月可降雹,6—7月降雹最多;降雹主要出现在14时—翌日02时;北部沿海相对较多。(2)降雹形势主要有5种类型:冷涡型降雹、低槽型降雹、横槽型降雹、西北气流型降雹、其他小范围降雹。根据冷涡中心位置冷涡型划分为两个关键区;低槽型可分为前倾槽、阶梯槽、较深低槽、与中低纬度共同作用的槽;横槽型降雹范围广、破坏性大;西北气流型存在连续性。(3)4类13种物理量具有不同分布特征和变异系数,均具有较好的代表性。不同月份、不同降雹影响程度和影响系统,物理量具有较明显差别。(4)0℃层高度在1370~5331m时,7种物理量可用于预报冰雹,K≥17℃、T850-500≥25℃、LI≤2℃、SRH≥0.1m^2·s^-2、SSI≥240、SWEAT≥100、Cape≥2J·kg^-1时可能降雹。6月、7—8月和其它月分别有3种、1种、3种物理量指标组合可用于预报冰雹,物理量的组合和数值有差异。     

我国华南江南春季雷暴气候特征分析

基于我国华南江南地区274个基本地面气象观测站数据、全国闪电定位数据以及欧洲中心的全球大气再分析数据(ERA-Interim),对1981—2017年华南江南地区的春季雷暴日采用经验正交函数分解方法(EOF),并与气象要素场做回归分析。得出以下主要结论:(1)我国华南江南地区春季雷暴活动高发区主要在广西东部至广东西部;其高峰期在下午18:00和凌晨4:00左右,且大多数雷暴活动持续时间不超过3 h;山区雷暴活动主要在傍晚至夜间;平原雷暴活动主要在白天,高峰在17:00及06:00前后;(2)华南江南地区的雷暴活动存在着3~5年的短周期和16年左右的长周期变化;(3)雷暴日距平EOF分析的前3个主成分累计方差贡献达到72.3%。按其向量场的方差贡献分型,Ⅰ型表现为华南江南雷暴活跃特征呈现较统一的变化规律。深厚西南低涡槽前、上干下湿的水汽层结、上冷下暖的温度层结为华南江南地区发生大范围雷暴天气提供良好的动力、水汽和位势不稳定条件,是华南江南雷暴活跃异常的主要模态;Ⅱ型表现为从华南南部到江西与浙江南部有一条西南-东北向、下宽上窄的雷暴活跃正距平异常区,而两侧为负距平异常区。其环流特征表现为温度整层偏冷,水汽整层偏湿,而西南槽前动力抬升有利于水汽抬升凝结触发对流形成雷暴;Ⅲ型表现为华南和江南地区雷暴活跃特征呈南北反位相异常,其分界线在26 °N附近。其环流特征表现为较强的干冷空气南下与南方暖湿空气在南岭山区对峙形成异常的垂直环流圈。在其上升支,低层干冷空气被卷入中高层使得中高层暖湿空气凝结释放潜热形成对流,造成华南地区多雷暴发生,而江南地区处于垂直环流的下沉支,整层湿度偏干,造成江南地区雷暴相对偏少。      

摄动理论在腐蚀管线随机地震反应分析中的应用

埋地管线的腐蚀情况十分复杂,具有随机性的特点,从而导致其在地震激励下的反应必然具有随机性。本文在管线腐蚀离散状态模型和弹性地基梁原理的基础上,利用随机摄动理论推导给出了腐蚀管线在地震激励下位移和应力的解析表达式,并进行了均值和标准差的计算。实例分析中将该方法计算结果与Monte Carlo模拟方法计算结果进行了对比。结果表明,利用随机摄动方法可以快速、精确地求得腐蚀管线的随机地震反应。     

层状地基端承桩扭转动力阻抗研究

本文采用Novak薄层法推导粘弹性地基的扭转动力阻抗,并将其用于考虑桩土相互作用的单桩扭转动力阻抗,又通过传递矩阵法将此公式推广到求取层状地基单桩扭转动力阻抗。而且本文以工程中常用的端承桩为例,推导了层状地基中端承桩扭转动力阻抗的简化计算公式。根据此公式,分析了频率、上覆软土层厚度、上覆软土层刚度等因素对单桩扭转动力阻抗的影响,分析结果表明,随着振动频率提高,层状地基中端承单桩的扭转动力阻抗的实部有缓慢下降的趋势,而虚部则增大;随着上覆土层厚度的增加,层状地基中端承单桩的扭转动力阻抗的实部和虚部均减小;随着上覆土层刚度的减小,层状地基中端承单桩的扭转动力阻抗的实部减小,虚部在低频段减小,而在较高频率段则增大。     

某长江大桥深度超过15米的砂土层液化势的综合评价

不同抗震设计规范的砂土液化判别方法或国内外其他有代表性的液化判别方法所采用的地震动参数和土性指标及其埋藏条件是不同的,因而采用这些方法对同一工程场地进行液化势预测时其评价结果通常有一些差异,甚至会得到相反的结论。为了给重大工程建设提供较为合理、可信的地基液化势预测结果,采用多种液化判别方法进行场地液化势的综合评价是比较客观的,也是必要的。本文结合某长江大桥桥基工程,采用建筑抗震设计规范的砂土液化判别方法、国内外有代表性的液化判别方法、有限元数值分析法等多种方法逐一对该工程场地砂性土层进行液化判别,并结合室内动三轴液化试验结果,对主桥墩不考虑冲刷条件和考虑一般冲刷深度5m条件时的砂性土层进行了液化势的综合评价,并将各土层的液化势分为液化、可能液化和不液化3个等级,得到了较为合理可靠的判别结果。     

工程结构三维疲劳裂纹最大应力强度因子计算

介绍了裂纹的类型、裂纹尖端应力场的奇异性。以一维问题为例,推导论证了奇异单元能够很好的反映裂纹尖端应力场的奇异性。应力强度因子一般表达式表明应力强度因子与载荷呈线性关系,并依赖于物体和裂纹的几何形状和尺寸。本文借助大型通用有限元软件ANSYS,采用位移外插法计算了三维表面裂纹前沿不同位置处的应力强度因子,并与《应力强度因子手册》基于实验的理论公式计算结果相比较。结果表明:有限元结果与理论解误差较小,裂纹最深处应力强度因子最大。     

Time domain method for calculating free field motion of a layered half-space subjected to obliquely incident body waves

In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.