简介:传统上,有限差分的差分系数一般可以通过泰勒级数展开法或优化方法来极小化频散误差得到。基于泰勒级数展开的差分法在有限的波数范围内精度较高,但在这个范围之外会产生较强的数值频散;基于最小二乘的优化有限差分法能在更大的波数范围内达到较高的精度,并可以在较小的计算需求内获得全局最优解。本文将基于最小二乘的优化有限差分法从二维正演模拟推广到三维,形成了计算效率高、高精度范围宽、适合并行计算的三维声波优化有限差分方法。频散分析及正演模拟表明本文发展的有限差分方法可以很好地压制数值频散。最后,将本文发展的有限差分方法应用到三维逆时偏移的震源波场延拓和检波点波场延拓中,并结合有效边界存储策略与checkpointing技术在GPU集群上实现三维逆时偏移以提高计算效率、减少存储量。三维逆时偏移试算结果表明本文三维优化有限差分方法与传统的有限差分法相比可以获得更高精度的偏移成像结果。
简介:利用ATOVS卫星辐射率资料和NCEP资料,采用中尺度数值模式WRF(V3.3)及其同化系统WRF-3DVar,设计了几组同化试验方案,对由0908号台风"莫拉克"导致的台湾南部2009年8月8—9日特大暴雨过程进行数值预报。结果表明,直接同化ATOVS辐射率资料对初始温度场、湿度场和风场均有较明显的改善,其中连续循环同化试验的改善效果最显著,同化后初始场增加了暴雨中心上空对流不稳定性,同时风场也表现出更强的低层辐合、高层辐散的中尺度特征。对比不同的辐射率资料发现,同化MHS资料对湿度场改进效果较好,同化AMSUA资料对温度场和风场改进效果较好,同时同化MHS和AMSUA资料改进效果整体上要好于前两者。在对降水预报定量检验中,同化试验的整体评分要高于控制试验,特别是降水临界值超过300mm后评分提高最为显著。连续循环同化试验对这次台风暴雨的预报与实况最为接近。
简介:对于中尺度数值天气预报来说,初始条件的准确与否已成为影响预报技巧的主要因素之一。现有的大气观测资料在时空分布上的不均匀,以及存在的观测误差,使得我们必须引进资料同化方法,为中尺度数值模式提供最优的初始场。由于传统的三维变分同化(3DVar)方法缺乏模式约束以及背景误差协方差矩阵(B矩阵)不具有流依赖性,因此本文提出一种基于历史样本投影的3DVar(HSP-3DVar)方法,它不仅具有流依赖的B矩阵,而且比传统的3DVar简单易行。为了评价HSP-3DVar的同化性能,我们基于区域暴雨预报模式AREM(AdvancedRegionalEtaModel)对其进行了观测系统模拟试验(OSSE),结果表明:HSP-3DVar能够有效融合观测信息,模式初值在各层的均方根误差都显著地降低。
简介:利用WRF以及WRFDA3.6.1的ETKF-3DVAR混合集合—变分同化方案,对2012年第15号台风"布拉万"进行了模拟,并分析协方差权重系数对"流依赖"属性和台风强度、路径的模拟效果的影响。通过单点实验发现:随着集合协方差权重系数增加,分析增量的"流依赖"属性就越发明显。当集合协方差权重系数大于0.5时,由于集合成员数量以及预报质量的限制,将会引入虚假的"流依赖"信息,使得分析增量的范围形状不再发生变化,强度有所减弱。在对台风"布拉万"路径和强度的模拟中发现:集合协方差权重系数取0.5以下的预报结果普遍好于系数大于0.5的预报结果,其中,当集合协方差权重系数取0.25时路径预报误差最小,当集合协方差权重系数取0.2时强度预报误差最小。
简介:介绍了如何从某种形式的度规出发,利用Mathematica软件计算联络,Ricci张量和总曲率,最终写出了真空Einstein场方程,并在附录中给出了相应的Mathematica程度,它可有效地协助用户进行与广义相对论有关的理论,观测研究。