最优化广义离散Shannon奇异核交错网格褶积微分算子地震波场模拟

笔者基于离散Shannon奇异核褶积微分算子(GDSCD)计算地震波速度应力方程的空间导数,推导了一阶GDSCD的具体形式,并提出了优化方法,即在频率域逼近平面波的真实导数,得到了不同半径和采样下限的最优窗函数系数,通过滤波响应分析算子精度,与多种数值方法对比以及模型测试表明,笔者构造的最优化GDSCD模拟地震波具有较高的计算效率和精度.     

Numerical simulation of seismic wave propagation in complex media by convolutional differentiator

We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.     

地震波传播的褶积微分算子法数值模拟

为了解决地震波场数值模拟中的速度与精度的匹配及局部信息与全局信息的耦合等问题, 该文借助新推出的基于Forsyte广义正交多项式的褶积微分算子, 将计算数学中的Forsyte多项式褶积微分算子应用到地震波传播的数值模拟中.复杂非均匀介质模型数值模拟结果说明了该方法的可行性和优越性.该方法同时具有广义正交多项式褶积微分算子的高精度和有限差分短算子算法的高速度.通过对算子长度的调节及算子系数的优化, 可同时兼顾波场解的全局信息与局部信息.      

非均匀介质中地震波传播的数值模拟

地震波场数值模拟方法多种多样,各种方法都有各自的特点.这里推出一种全新地震波场模拟方法--基于Forsyte广义正交多项式的迭积微分算子法,该方法将计算数学中的Forsyte多项式,应用到地震波传播的数值模拟中,它同时具有广义正交多项式迭积微分算子的高精度和有限差分短算子算法的高速度.通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息.复杂非均匀介质模型数值模拟结果说明了该方法的可行性.将该方法的计算结果与傅氏变换伪谱法、错格高阶有限差分法相比较,结果说明,该方法在波场模拟方面具有较好的发展潜力,并具有自身独到的优越性.     

华南飑线升尺度增长过程中的多尺度能量相互作用分析

采用WRF模式对华南飑线的升尺度增长过程进行模拟,利用Barnes滤波将模式数据分解为三个尺度,分别代入相应的能量方程中进行计算,从能量角度研究飑线升尺度增长过程中动能和位能的变化,以及三个尺度系统能量的相互转化.研究表明:动能的变化与飑线过程中各尺度系统的演变有较好的对应,β中小尺度对流的发展对应β中小尺度系统动能的...     

基于群体搜索的串行蒙特卡罗反演方法的并行算法(英文)

随着并行计算技术的发展,非线性反演计算效率在不断提高,但对于基于单点搜索的非线性反演方法,其并行算法的实现则是一个难题。本文将群体搜索的思想引入到基于单点搜索的非线性反演方法,构建了并行算法,以量子蒙特卡罗方法为例进行了二维地震波速度反演及实际资料波阻抗反演,并测试了使用不同节点数进行计算的效率。计算结果表明:该并行算法在理论和实际资料反演中是可行的和有效的,具有很好的通用性;算法计算效率随着使用节点数的增加而提高,但算法计算效率的提高幅度随着使用节点数的增加逐渐减小。     

PS转换波界面二次源法射线追踪

界面二次源法是最近提出的一种最小走时射线追踪方法,尤其适合层状介质中走时和射线路径的计算.该方法相对于传统的最小走时树方法（如Moser法）,仅在物性界面上设置二次源,射线路径的方向只在层界面处发生改变,该方法最大程度地消除了射线路径的锯齿状现象,同时也避免了低变速区的射线路径多值现象,因此,它具有更高的追踪精度和效率.本文采用界面二次源法在各向同性介质中实现了PS转换波射线追踪,理论模型的计算证实了界面二次源法追踪PS转换波的准确性和高效性,同时该方法在各向异性介质中也很好地追踪出分离的PSV波和PSH波,因此该方法有利于横波分裂在地震勘探中的研究和应用     

一种飓风数值预报的初值化方法

本文使用具有Gallelean变换和详细动量、热量垂直输送通量描述的飓风动力方程组,由切向风速得到垂直运动场、加热场等,结果表明了垂直运动场和加热场具有明显的螺旋雨带特征,温度场有暖心结构及其位势高度场具有深低压特征,这些特征和实际观测结果相吻合。这个方法为飓风数值预报初始场提供了可行的方法。 用数值方法预报飓风的发生发展和路经及降水分布是飓风研究的一种主要方法。但是,对于数值预报飓风的初始场的确定仍是一个没有彻底解决的问题。在过去的初始场处理中,大致分几种,一种是初始场全部采用理想化的低压扰动,另一种则是用实际资料得到环境流场然后再加上一个理想化的     

热岛效应强迫下的中尺度环流的动力特征及极限风速的一种解析表达

本文采用考虑热岛造成的局地热力扭动的Boussinesq二维垂直剖面的动力方程组,利用高截断谱方法,得到一类稳定的中尺度非线性平衡态和极限风速的解析表达,证明了局地中尺度环流是热岛造成的局地热力扰动和平流非线性共同作用的结果;并进一步讨论了热岛强迫下的这类局地环流的某些动力特征。     

凝结潜热释放和地表热通量对一次飑线过程的影响

针对2013年7月4—5日一次产生大风、短时强降水强对流天气的飑线过程进行数值试验,研究了凝结潜热释放和地表热通量在中尺度对流系统的发生发展过程中的影响。结果表明:(1)凝结潜热释放对飑线系统有重要的作用,对飑线发展维持、移动及成熟阶段的垂直结构都有着一定影响。(2)当飑线系统进入成熟阶段后,小尺度的积云对流中的凝结加热作用于中高层大气,加强了高层辐散低层辐合的配置,而对流系统内垂直上升运动的加强又导致凝结加热作用更强。(3)凝结加热作用极大地促进了飑线的低层入流、高层出流的垂直结构,以及系统前方中层入流和高层出流之间形成的间接垂直反环流,促进高空动能下传的同时,使系统前方对流不稳定性增加,新的对流单体易于触发。(4)潜热释放间接增强了对流内部的冷性下沉气流,进而加强了低层的阵风锋,使得新的对流单体能在飑锋处触发,飑线以新老对流单体交替的方式向前移动。(5)在高低空急流的垂直耦合促进飑线发展的同时,成熟阶段飑线中较强的凝结加热对高低层急流有反馈作用,一定程度上增强了这种垂直耦合。(6)地表感热、潜热通量等边界层非绝热过程在对流系统的触发和发展中起到了较为重要的作用。地表热通量在白天加强了对流边界层的湍流混合作用,由此产生的特征维持到了夜间,形成了适合飑线触发的层结条件。另外,在对流形成之前,与地表潜热通量相关的边界层加湿作用为对流的爆发贮存了丰富的水汽和不稳定能量。