基于小波分析的西北区智能网格气温客观预报方法的检验评估

基于中国气象局(China Meterological Administration,CMA)高分辨率数值预报产品、欧洲中期数值预报中心(the European Center for Medium-range Weather Forecast,ECMWF)精细化数值预报产品和国家级地面观测站数据,采用小波分析方法及滑动训练、最优融合等技术对模式误差序列进行时频处理,实现了对模式系统误差和局地误差的订正,发展了西北区智能网格气温客观预报方法(northwest intelligent grid temperature objective prediction method,NWTM)。以2017年3月—2018年2月数据作为训练样本,对2018年3月—2019年1月西北区239个国家基本站进行检验。结果表明:1)NWTM对CMA和ECMWF两种模式产品的气温预报能力有显著的提升;随着预报时效增长,两种模式订正产品的误差增大。2)NWTM对ECMWF西北区最高气温的订正效果要明显优于CMA,但就最低气温而言,NWTM对CMA的订正效果更为显著。其中,就24 h最高气温而言,ECMWF在宁夏的订正效果最好,CMA在青海的订正效果最好;而对于24 h最低气温的预报,CMA在西北4省的订正效果相差不大,ECMWF在陕西的订正效果最好。3)空间误差检验表明:针对最高气温的预报,ECMWF订正产品的订正能力明显优于CMA,特别是在甘肃河西走廊和中东部、陕西北部和南部、宁夏中南部及青海大部。就最低气温的预报而言,ECMWF和CMA对甘肃河东和陕西南部的订正能力较好;ECMWF订正产品在宁夏中南部及青海南部的订正能力高于CMA,而CMA订正产品在陕西中部的订正能力更优。     

智能网格预报时空协调一致关键技术研发

自2016年起,我国大力推进气象业务向"智能化、无缝隙、精准化"发展。2017年宁夏积极推进智能网格预报业务体系建设,在网格预报智能订正、主客观融合订正、格点/站点/落区一体化、灾害天气智能识别等方面开展了一些尝试,取得了一定的成果,目前各项成果均引入到宁夏智能化天气预报业务平台中。其中一项关键技术,是智能网格预报的时空协调一致技术。该项技术在主客观融合订正、格点/站点/落区一体化中发挥了重要作用。本文重点针对这一关键技术的设计思路和实现方法进行介绍、讨论。     

自适应高斯频域滤波器在天气雷达中的应用

多普勒天气雷达的基数据(反射率因子、速度和速度谱宽)经常会被地物杂波所污染,这对基数据估计的准确性影响很大,进而又会影响导出产品。因此,抑制地物杂波是所有多普勒天气雷达获取高质量数据产品的必要步骤。作者首先介绍了自适应高斯频域滤波器的研究背景,其次分析了多普勒天气雷达地物杂波及天气回波信号的特征,再次给出了IIR椭圆地物杂波滤波器的原理以及自适应高斯频域滤波器的滤波算法,最后根据实际的雷达回波信号,分析了IIR椭圆地物杂波滤波器和自适应高斯频域滤波器的杂波抑制性能,并对结果进行了分析和比较。结果表明:不使用杂波图的自适应高斯频域滤波器的杂波抑制性能优于使用杂波图的IIR椭圆滤波器,且满足实时处理的要求。     

空间网格结构多维多点随机地震响应分析的高效算法

将林家浩教授提出的“虚拟激励法”进一步推广应用于空间网格结构多维多点非平稳随机地震响应分析,推导了多维虚拟激励随机振动分析方法的理论公式,给出了峰值响应估计方法,并讨论了多维地震动的随机模型及参数选取,通过编制的专用计算机程序分析了网壳结构的随机地震响应。本方法自动包含了参振振型间及各输人地震分量间的相关项,计算精确、快速,非常适合分析频率密集型空间网格结构的随机地震响应,是一种高效的随机振动分析算法。     

基于波前面三角形网格剖分的波前重建法三维射线追踪

在许多地震反演和偏移成像方法中,都要涉及到射线路径和旅行时的计算.本文将波前面三角形网格剖分和三维波前重建法射线追踪技术结合使用,实现了射线路径和旅行时的准确快速计算.三维波前重建法射线追踪过程中可以保证稳定合理的射线密度,克服了常规射线追踪方法存在阴影区的问题.波前面三角形网格剖分在描述和拆分波前面时更加准确有效,而且不需太多的网格数目,从而提高了射线追踪的精度和效率.该方法在三维复杂构造成像方面有独特的优势,目前在实际的Kirchhoff 偏移中的已经有相关应用.     

基于渗流和管流耦合的管涌数值模拟

堤坝地基的渗透变形过程实际上是“土中水”转变为“水中土”的过程。在渗透变形发生的集中管涌通道区域,采用常规渗流分析理论,单纯增大管涌通道渗透系数的方法是不太合适的。在未发生渗透变形的区域,用常规渗流理论计算;在管涌通道区域,用管流理论,公共边界上两者之间水头相等、流量大小相等且方向相反,能够较好符合渗透变形的发展规律。为了适应计算过程中内部边界条件不断变化的特点,采用无网格法伽辽金法（element free Galerkin method,EFG）对渗流场进行计算。算例计算表明,这种渗流-管流耦合的方法能够模拟管涌通道绕过防渗墙等复杂的发展过程。     

网格状带齿加筋体界面特性的宏细观力学机制研究

网格状带齿加筋是一种新型的土工加筋技术,汲取了传统网格状与条带式带齿加筋的优点,其筋土界面特性是加筋结构设计的基础。利用室内拉拔试验从宏观角度分析筋材与土的界面特性,结合颗粒流程序,筋材用平行黏结模拟,从细观角度分析筋土接触界面的产生和发展演化规律。数值结果表明：齿筋的存在对界面处颗粒的移动及剪切带的宽度影响显著,其破坏模式呈现波浪型,在齿筋附近出现应力集中现象。筋土界面处筋材与土颗粒摩擦、咬合和阻抗作用提供“似黏聚力”是改变网格状带齿加筋体系作用机制的内在原因。研究成果对于明确网格状带齿加筋拉拔的机制及模型的建立都具有意义。     

油藏精细地质模型网格粗化算法及其效果

在前人研究基础上, 根据DP(Dykstra-Parsons)系数能定量评价储层非均质性, 微网格块的渗透率值粗化后, 其等效渗透率的上、下限(C_min、C_max)能反映渗透率的各向异性的特点, 提出了一种运算速度快和相对有效的网格粗化算法。该算法能考虑到储层非均质性对不同方向渗透率值的影响, 且求解过程相对简单。应用该方法对鄂尔多斯盆地中部某油藏的陆相储层的精细地质模型进行了网格粗化计算, 然后在粗化后的模型上进行油藏数值模拟研究, 同时针对研究区地质背景和产出流体微可压缩的物性特征, 首次利用流线模拟器对精细地质模型进行了油藏数值模拟研究, 并以此结果为标准, 对该网格粗化算法时效性进行了系统评价。分析表明, 该算法具有较快的计算速度和较高的可靠性, 是解决储层非均质强、物性差的陆相成因油藏精细油藏数值模拟的一种行之有效的手段。      

曲线坐标系下的完全匹配层吸收边界条件

在地震波数值模拟中,需要采用吸收边界条件以吸收人为边界反射。本文针对曲线坐标系下的二阶弹性波方程提出了一种完全匹配层(PML)吸收边界条件。与直角坐标系下的PML吸收边界条件类似,曲线坐标系下的PML吸收边界条件是一种在频率域中给出的人工边界条件,由相应的复坐标变换得到。在变换到时间域后,完全匹配层中将出现复杂的卷积运算。为了避免这些卷积运算,引入了4个中间变量。为了简化自由边界条件,采用正交贴体网格对起伏地表模型进行网格剖分。数值算例表明,该方法可以有效消除人为边界反射。     

基于移动网格技术的热融湖动态演化过程数值模拟

基于移动网格技术建立了热融湖动态演化有限元数值模型, 研究了青藏高原多年冻土区典型热融湖动态演化过程, 分析了热融湖半径、深度的变化过程及其对湖底及周围多年冻土温度状况的影响. 结果表明:在移动边界热融湖模型中, 热融湖半径以0.7 m·a^-1的速度近线性地增大; 随着下伏高含冰量冻土的融化, 热融湖深度增加先慢后快, 最后逐渐减小趋于稳定. 热融湖深度和半径从5月末至翌年1月末增加显著, 在2-5月间基本保持稳定. 伴随着热融湖的扩展, 地表边界逐渐演变为湖底边界, 热融湖的热影响范围逐渐增大. 在固定边界热融湖模型中, 其热影响会逐渐趋于稳定, 由于初始尺寸大, 其湖底多年冻土退化速率大于移动边界模型, 而远离湖边的多年冻土退化速率要小于移动边界模型. 如果不考虑热融湖边界随时间的变化, 可能会高估热融湖对湖底多年冻土的热影响, 而低估其对附近多年冻土的热影响.