无网格方法及其在地球物理勘探中的应用与展望

产生于20世纪70年代的无网格法是一类新兴的数值方法,现已被广泛运用于科学计算及工程问题并取得了显著成效.相较于有限元方法,无网格法具有网格依赖性弱、精度高、收敛速度快、方便构造高阶拟合函数等优点.经过几十年的发展,无网格法取得了长足的进步,诞生了多种方法.无网格方法之间的主要区别在于构造近似函数的不同和加权余量法的选择.本文概述了无网格法的发展历史和研究进展以及在地球物理学中的应用状况;详细阐述了多种无网格的近似方法,包括：移动最小二乘近似、径向基函数近似、核积分近似、再生核粒子近似等,以及常用的加权余量法—配点法、伽辽金法、局部伽辽金法.并对这些无网格方法的差异性作了系统对比分析,同时指出了在地球物理领域的适用范围.还例举了无网格法在地震勘探领域的波动方程、重磁勘探领域的位场延拓以及异常体正演的求解过程.最后讨论了无网格方法应用于地球物理勘探领域时存在的困难和挑战,并分析了未来无网格方法的应用前景和潜力方向.     

纬向切变基流中的非线性Rossby波

本文在半地转近似下讨论了纬向切变基流中的非线性正压Ｒｏｓｓｂｙ波，给出了存在波解的切变基流条件。     

计算最小走时和射线路径的界面网全局方法

用慢度分块均匀正方形模型将介质参数化，仅在正方形单元的边界上设置计算结点，这些结点构成界面网．根据Ｈｕｖｓｅｎｓ和Ｆｅｒｍａｔ原理，由不断扩张、收缩的波前点扫描代替波前面搜索，在波前点附近点的局部最小走时计算中对波前点之间的走时使用双曲线近似，通过比较确定最小走时和相应的次级源位置，记录在以界面网点位置为指针的３个一维数组中．借助这些数组通过向源搜索可计算任意点（包括界面网以外的点）上的全局最小走时和射线路径．这一方法不受介质慢度差异大小限制，占内存少，计算速度较快，适于走时反演和以Ｍａｓｌｏｖ射线理论为基础的波场计算．     

双曲线配合法在某路基沉降反分析中的应用

本文根据哈尔滨某公路路基实测沉降资料,应用双曲线配合法推算了地基土的最终沉降量,并用太沙基渗透固结理论反分析出地基中粘土层的平均固结系数和渗透系数。应用太沙基方法,将反算出的渗透系数计算路基沉降,结果比选用实验室测得的渗透系数计算的沉降更接近实际沉降值。     

边界层内冷锋流场的动力学特征

文中求解了锋面存在时地转动量近似下的大气边界层运动方程 ,得到了边界层内冷锋流场的一些特征 ,如冷锋坡度随地转涡度增加而增加 ,随地转风速时间倾向的增加而增加 ,随沿锋面传播方向的热成风分量的减少而增加。而边界层内冷锋面上下的流场与锋面坡度、地转风及其时空变化特征有关 ,共同特点是在冷锋面高度以下有下滑运动 ,而其上有一层上滑运动区。     

变化涡动沾性系数的Ekman动量近似模式风场结构

在边界层动力学中，涡动粘性系数是影响边界层风场结构的一个重要参数。本文利用边界层动力学中的Ekman动量近似理论，给出了涡动粘性系数随高度缓变条件下的Ekman动量近似边界层模式解，着重讨论了边界层的风场结构、水平散度、垂直涡度以及边界层顶部的垂直速度。结果分析表明：与常值涡动粘性系数情况相比，在边界层低层随高度增加的涡动粘性系数可以导致低层边界层风速随高度迅速增加，即风速垂直切变增加，同时风速矢与地转风之间的夹角减小。惯性项作用可以导致上述作用在气旋性区域减小、而在反气旋性区域增大。随高度增加的涡动粘性系数导致水平散度绝对值、垂直涡度绝对值以及边界层顶部的垂直速度绝对值在气旋性区域减小，而在反气性旋区域增大。涡动粘性系数与惯性之间的非线性相互作用是边界层动力学中重要过程。     

土坝下软基涵洞施工力学分析

通过逐级新增单元的有限元方法，模拟了土坝下软基涵洞周围土体分层填筑的施工过程。土体采用邓肯-张双曲线模型，涵洞与土体之间采用Goodman型摩擦接触单元，较好地反映了涵洞周围土体的位移场、应力场与涵洞结构内力，为实际工程设计提供了可靠依据，是涵洞结构计算的一种有效方法。     

合成平面波数据的Rytov近似叠前深度偏移

作者把Berkhout(1992)提出的合成平面波技术应用于线性Rytov近似叠前深度偏移，提高了波动方程叠前深度偏移的计算效率。通过在Marmousi模型数据试验，我们认为在达到基本相同成像效果的条件下，利用合成平面波技术的波动方程叠前深度偏移相对于基本共炮集数据的波动方程叠前深度偏移可提高计算效率约10倍。