复杂介质地震波传播的褶积微分算子数值模拟h

将佛尔塞（Forsyte）广义正交多项式微分算子地震波场正演模拟算法运用于复杂非均匀介质模型的波场数值模拟中,并比较了该方法在计算效率和计算精度方面与有限差分方法和伪谱法的差异. 数值结果表明,这种广义正交多项式微分算子法计算速度快、精度高,对计算资源需求低,是一种颇具潜力的数值模拟方法.      

基于辛格式离散奇异褶积微分算子的弹性波场模拟

本文发展了基于辛格式离散奇异褶积微分算子(SDSCD)的保结构方法模拟弹性波场,求解弹性波动方程时,引入辛差分格式进行时间离散,采用离散奇异褶积微分算子进行空间离散.相比于传统的伪谱方法,该方法提高了计算精度和稳定性.数值结果表明SDSCD方法可以有效地抑制数值频散,为解决大尺度、长时程地震波场模拟问题提供了合适的数值方法.     

一次梅雨锋暴雨的边界层热通量输送与湍流动能收支分析

利用中尺度WRF模式对于2007年7月一次典型的梅雨锋暴雨过程进行了高分辨率数值模拟,对于边界层内的热通量输送和湍流动能的时空变化特征,以及湍流动能各收支项的分布及变化特征进行了分析。结果表明,降水发生时段内边界层热通量和湍流动能的时空分布特征与晴空日变化特征表现出显著不同,潜热通量随高度自下而上呈现"正—负—正"的分布,感热通量以负值为主,负值中心高度与潜热通量由负转正的高度相对应,湍流动能的发展高度与持续时间都有所增加,降水区近地面湍流动能弱于其他区域,但是在468 m以上高度则显著强于其他区域。降水区湍流动能的来源主要是平均风切变所产生的机械湍流,浮力作用与粘性耗散在降水期间消耗湍流动能,湍流输送作用将低层的湍流动能输送至较高的高度,使低层减小而高层增大,临界高度与湍流动能的大值中心高度对应。     

不同边界层参数化方案对一次梅雨锋暴雨过程湍流交换特征模拟的影响

利用WRF模式结合不同的边界层参数化方案,对2007年7月3—5日发生在江淮流域的一次梅雨锋暴雨过程进行多组数值模拟试验。结果发现,边界层方案的选取对于降水的落区和强度模拟会产生较显著的影响;在降水率及地面要素的模拟上,各方案在降水中后期的模拟差异明显大于降水发生阶段;不同边界层方案的选取对于降水时段内的水平风场、垂直运动和假相当位温的垂直分布都产生影响,直接影响降水时空分布的模拟;不同方案都模拟出了在降水发生之后不同于晴空日变化的湍流动能垂直分布,经分析发现与局地较强的垂直风切变和近地面强湍流气团被抬升有关,而浮力项起着耗散作用;各方案的湍流交换特征与湍流动能特征基本吻合,相比于其他方案,MYJ方案在降水区域的湍流动能及湍流交换强度明显偏弱,对热通量的输送也偏弱;GBM方案在边界层内的湍流混合偏弱而在边界层以上湍流混合显著偏强,热通量输送在边界层以上的高度上误差明显,影响了对降水区域气象要素的模拟能力,仍需要进一步改进。      

地震数据重建方法原理及运用

在数值分析领域中,许多方法都是借助于插值公式导出的.几乎所有的经典数值微分、数值求积和常微分方程的数值积分公式,都可以从插值公式中推导出来.地震勘探领域常见的插值重建方法按照数学原理大致可以分为三类:适用于均匀采样的预测误差滤波插值方法,基于算子的插值方法(如DMO和偏移算子)和在最小二乘意义下估算特定数学变换系数进行插值的方法.这里从数学原理上系统地阐述了这几种插值方法的基本原理,并利用数值试验和实际数据进行了检验.     

STUDY ON THE GENESIS OF ASYMMETRICAL DISTRIBUTION CHARACTERISTICS OF PRECIPITATION ASSOCIATED WITH THE TYPHOON HAITANG (2005) FROM THE VIEW OF ATMOSPHERIC FACTOR

The distribution of precipitation field from the typhoon Haitang (2005) during its landing on Fujian province shows obvious asymmetric feature. Based on the NCEP/NCAR FNL (Final Analysis) data, this study reveals the contributions of atmospheric factor to the asymmetrical distribution characteristics of precipitation associated with the typhoon, through the analysis of water vapor condition, vertical ascending motion condition, the calculation of the dry Q vector and its decomposition, and adiabatic heating in the air column of 1000hPa -600hPa (lower atmosphere) and 500hPa-100hPa (upper atmosphere). The results are as follows: (1) In the lower atmosphere, the humidity on both sides of typhoon path can be equivalent, while it is more wet on the right side than left in the upper atmosphere, which obviously presents asymmetric distribution characteristics. (2) Both range and intensity of the vertical motion on the right side are wider and stronger than counterparts on the left side no matter in the lower or upper atmosphere. (3) In the upper atmosphere, forcing role of atmosphere in vertical upward motion on the right side of typhoon path is the same as that on the left, while it is significantly different in the lower atmosphere, which is significantly broader in scope and stronger in the intensity, along with obvious asymmetric distribution characteristics. In addition, the further analysis of the Q vector decomposition indicates that the forcing effect of mesoscale weather systems on vertical upward motion is stronger than that of large scale weather systems in the lower atmosphere. (4) The adiabatic heating always exists on both lower and upper atmosphere, and the range and intensity of the adiabatic heating forcing showed asymmetric distribution on both lower and upper atmosphere. (5) In a summary, the upper atmosphere humidity conditions, the forcing role of lower atmosphere in vertical upward motion, especially, to mesoscale weather system, and adiabatic heating in the lower atmosphere, all show similar asymmetric distribution characteristics to that of precipitation field from the typhoon Haitang (2005), that is to say, the atmospheric factors as mentioned above are all contributed to genesis of the asymmetric distribution characteristics of precipitation.     

基于改进粒子群算法的地震标量波方程反演

针对标准粒子群优化(PSO)算法存在易出现早熟而陷入局部最优以及进化后期收敛速度慢等缺陷,通过考虑粒子所处位置间相互作用,提出了一种改进的并行粒子群优化算法.由于引入粒子位置间的相互影响,减少了粒子搜索过程盲目性,因此能有效提高算法的收敛速度.数值试验表明,这种改进的粒子群算法适用于二维标量波方程的速度反演,且算法具有...     

华南一次飑线升尺度增长过程的机制分析

利用中国自动站与CMORPH降水产品融合的逐时降水量0.1°×0.1°网格数据集、1°×1°NCEP再分析气象资料、常规观测资料和高分辨率的WRF数值模拟资料对2016年4月13日前后发生在华南的一次飑线升尺度增长过程进行了动力和热力环境条件分析。飑线形成前,对流层高层有急流,华南位于急流入口区右侧,有利于高空辐散。高空有弱槽存在,并且温度槽落后于高度槽,有利于槽发展。近地面上有低压系统控制,并且有风场的辐合,加强了抬升运动。模拟的结果表明:WRF模式成功模拟了此次飑线升尺度过程,飑线的降水量和落区较为接近实况。与以往的飑线过程相比,此次飑线过程中,升尺度阶段有效位能较初始对流阶段要高,但是整个华南地区对流有效位能始终处于中低水平。除了天气尺度系统的影响,多种中尺条件的配合对这次飑线升尺度发展很有利。低层的垂直风切变、低位涡和较高温度递减率是触发条件,提供了一定的不稳定条件和抬升条件。低层风场转向加强,垂直风切变的方向改变,充足的水汽供应和后部入流的加深加剧了不稳定和垂直抬升运动,是促进β中尺度飑线升尺度的重要条件。冷池在升尺度阶段不明显,但冷池对维持成熟阶段的α中尺度飑线结构很重要。     

各向异性介质地震波场的优化褶积微分算子法数值模拟

在前人工作基础上,通过对窗函数参数进行优化实现了对基于Shannon奇异核理论的交错网格褶积微分算子的优化过程.应用这种优化褶积微分算子方法对各向异性介质进行了数值模拟,讨论了优化褶积微分算子法模拟的PML吸收边界条件以及稳定性条件,分析了弹性波在此类介质中的传播特征,并与高阶交错网格有限差分方法进行了对比.数值实验结果表明,该方法适用于各向异性介质中弹性波场模拟,精度高,稳定性好,是一种研究复杂介质中地震波传播的有效数值方法.