一种层析成像的正演和反演方法

作者在本文中提出了一种层析成像的正演和反演方法，首先对求解区作网络划分，由给定的正演速度模型在网格点上双线性插值，获得整个网格节点的速度分布，然后，对波阵面作一近似，在此基础上，用最小走时最短路径射线追踪法完成正演，反演则是通过正演的方法实现的，理论模型算例表明，本方法效果较好。     

多道瞬态瑞雷波法在探测煤矿采空区中的应用

采空区探测是工程建筑和煤矿生产中亟待解决的问题。常规物探方法在对落后的“房—柱”式开采技术造成的面积较小、埋深较浅的采空区进行探测时存在漏报或误报的情况。根据采空区的地球物理特征,探讨多道瞬态瑞雷波技术的探测机理,并通过勘探实例阐述了该方法在浅层煤矿采空区探测中的实用性和有效性。     

断层带爆炸裂隙区范围及裂纹扩展长度的研究

常用求算岩石内爆炸裂隙区半径的计算公式中忽略了冲击波的作用、粉碎区的存在以及爆生气体的准静态作用。醉对该公式进行修正后并引和断层带应力波衰减指数α，并对尖力波损伤后的裂隙区的扩展长度进行分析计算，从而得出：考虑“预报伤”作用的爆炸裂纹长度大于不考虑“预损伤”情况下的裂纹长度。     

略论珠穆朗玛峰重力值的推估

仅根据邻近点的重力与高程资料,采用了4种有关公式,有效地推估了珠穆朗玛峰顶上的重力值,该值为(976 970±7)×10-5m     

包络地形和重力波拖曳对区域气候模拟效果的影响

用NCAR区域气候模式(RegCM2)的数值试验，研究了包络地形和地形重力波拖曳作用对1991年5、6月份中国区域内气候模拟效果的影响。通过试验结果的对比分析发现，在模式中引入地形重力波拖曳作用可在一定程度上改善模式对气候系统和要素的模拟效果，使模拟结果更符合气候实况。包络地形对模拟结果也有一定的改善，但对于细网格的区域气候模式(RCM)来说，其效果不如重力波拖曳明显。不论是地形重力波拖曳还是包络地形，在模式积分的初期，它们的作用并不明显，随着积分时间的增长，它们对模拟结果的影响程度增大。     

濮阳市11月份寒潮天气预报

利用1965-1995年天气图和气象观测资料，分析了濮阳市11月份寒潮的时空分布特征、路径及天气形势，确定了预报指标。     

遗传算法中的光滑约束反演及其在青藏高原面波研究中的应用

光滑约束技术在线性反演中具有重要的作用，但在遗传算法的反演中则很难直接施加于模型参数，其原因是采用光滑处理后的模型参与迭代后，模型的多样性受到很强的压制，并在少量的迭代过程中使种群的各个模型趋向一致，从而得不到满足条件的最优解．本文给出了一种可用于遗传算法反演的间接光滑约束方法．该方法将遗传算法迭代过程中产生的模型经处理后得到的光滑模型，作为误差函数计算的输入模型．迭代过程仍采用原模型，避免了模型的多样性损失，在面波反演和接收函数反演的试验中取得了良好的效果．我们利用该方法对青藏高原地区的瑞利波相速度资料进行了反演，揭示了青藏高原中部地区S波速度结构的横向变化特征．结果表明，青藏高原北部地区地壳Ｓ波速度较南部地区低；大多数路径在15～40ｋｍ 深度范围内，存在12～25ｋｍ 厚的地壳低速层；上地幔低速层位于100ｋｍ 深度以下，厚度主要在40～80ｋｍ 范围内变化，个别路径可达100ｋｍ 以上．安多台以北、玛沁和玉树以西之间，在上地幔90～230ｋｍ 深度范围存在明显的低速层，最低速度约4.2～ 4.3km／s．根据不同路径的S波速度结构和前人的资料，我们认为印度板块的俯冲可能以雅鲁藏布缝合带附近为界．      

利用强震记录分析场地的地震反应

利用唐山近场强震记录对如何分析场地反应进行研究，分析比较了多种方法:线性反演法、传统的谱比法以及接收函数法．通过比较研究，发现各种方法都能给出土层场地的卓越周期，但每种方法给出的场地反应不完全相同．其中，传统的谱比法与线性反演法的结果较为接近，而接收函数法的结果则同其它两种方法的结果相差较大．在线性反演法中，给出了Ｓ波的品质因子QS，QS 在0.5～32Ｈz范围内与频率有关:QS（f）＝67 f 1.1．      

基于螺旋线上谱因式分解的地震波场隐式辛算法

均匀介质、复杂各向同性介质和各向异性介质中的地震波传播过程,可用统一形式的标量声波方程描述.考虑到在无损耗条件下,地震波方程描述了地震波场这一个无穷维的哈密顿体系随时间的演化过程,该过程为一个单参数连续辛变换,因而可以在其哈密顿形式表述下导出其辛格式.与显式辛算法相比,隐式辛格式对应的隐式辛几何算法具有无条件稳定的特点,可以允许较大的计算步长.但是由于隐式算法不可避免地面临高阶矩阵的求逆,其每一步的计算速度较慢.为实现矩阵快速求逆,文中采用了螺旋边界条件下谱因式分解的方法.在螺旋边界条件下,需要求逆的矩阵化为带状矩阵,而且其各列非零元素的位置和大小具有非常好的相似性,因而可以采用谱因式分解的方法实现快速LU分解.文中采用二阶精度的隐式蛙跳辛格式和谱因式分解方法,计算了常速度、层状介质和Marmousi模型中的波场.计算表明,隐式辛算法不失为波场计算的一种好方法.     

中层大气中重力波饱和机制的数值分析

通过对向上传播的重力波波包在中层大气中的非线性传播过程进行数值模拟,讨论中层大气中重力波的饱和机制.数值模拟结果表明,向上传播的重力波波包的扰动振幅在接近波包的本征水平相速度之前随高度单调增长,而当波振幅接近本征水平相速度时,在对流不稳定区域出现等位温面的翻转,同时波振幅的增长达到饱和（波振幅随高度不发生变化）.小尺度对流在等位温面的翻转和波饱和发生后产生,随后波包开始破碎,这些非线性过程的最终结果产生湍流.表明导致重力波饱和的关键因素是等位温面的翻转而不是诸如波破碎、湍流、波-流相互作用等其他的一些物理过程.