用交错网格有限差分法计算三维频率域电磁响应

用交错网格有限差分法(SFD)，实现了三维频率域电磁场响应 的数值模拟. 该方法适用于任何方向的磁偶极子源. 经与解析方法、积分方程等 其他方法的计算结果对比表明，交错网格有限差分法结合散度校正和不完全乔累斯基分解预 处理的双共轭梯度迭代方法进行正演计算，速度快、精度高、结果稳定，能适应三维复杂介 质的数值模拟，为三维电磁反演奠定了基础.     

利用高阶交错网格有限差分法模拟地震波在非均匀孔隙介质中的传播

烃类储集层是一种复合多相介质，在固体颗粒的空隙中含有气体或液体. 研究弹性波在该类地层中的传播规律对于油气勘探开发，特别对于全波列声波测井有重要意义. 为了提高孔隙弹性介质数值模拟的计算效率，本文采用改进显式交错网格有限差分算法取代常用的空间域四阶和时间域二阶的速度 - 应力有限差分算法，算法的空间域为八阶、时间域为二阶. 虽然计算的时间步长略小于空间域四阶的情形，但高阶有限差分算法可以选择较粗糙的网格，因此补偿了计算的低效；同时高阶交错网格有限差分算法的空间频散性比低阶算法小. 利用该算法计算了一个两层模型的波场，同时还模拟了等效弹性和孔隙弹性模型中波的传播. 结果表明慢波及其影响明显，尽管慢波衰减很快，但被某一界面反射后，转换形成的P波和S波仍以正常的方式传播，且比慢波衰减小.     

地铁诱发地面运动的衰减规律的研究分析

以地铁车站间的隧道为研究对象，在研究和分析有关资料的基础上，用快速拉格朗日有限差分法软件FLAC建立了二维数值分析模型，并以列车振动加速度时程作为模型输入，进行了隧道衬砌－地层的整体动力学分析，分析了列车运行产生的振动在地面的衰减规律，得到了一些结论，可为地面建筑的防振减振控制提供参考。     

求解波动方程的任意差分精细积分法

在地震偏移成像技术中，常常要反向外推波场，因此涉及到波动方程的数值求解问题。本文提出了一种求解偏微分方程的新的半解析方法——任意差分精细积分(ADPI)法。其大体思路是：空间域上作坐标离散，但不采用古典差分法的等分离散方式，而是一种相对自由的、任意的离散法；时间域上则采用子域精细积分的方法，既保留了精细积分法的高精度，又克服了工作量大、占用存储大等缺点。该方法具有精度高、带宽小、稳定性好等多项突出优点，并且可以灵括处理各类边界条件。本文从简单一维、较普遍一维、以及二维3种情形讨论该方法对波动方程的具体应用，结合实际的算例，详细分析了各种算法的可行性和精度特点。     

中国煤和烃源岩镜质组的激光诱导荧光显微特征及其应用

应用自行组装的激光诱导荧光显微光度计(LIFM)系统地研究了中国煤和烃源岩中镜质组的荧光强度及荧光变化，建立了镜质组的荧光变化与成熟度之间的相关模式图(R-I-VRo)，该图可解决镜质组反射率抑制问题，并将烃源岩成熟度评价范围扩展到VRo达2．0％左右，该方法对于正确评价烃源岩成熟度具有广泛的应用前景。     

高精度差分GPS技术在海洋资源调查中的应用

本文以差分GPS技术多年来在海洋资源调查中使用为例，对差分技术在海洋资源调查中的应用和定位质量评定方法作探讨，并对两种综合导航定位系统作了介绍。     

应用激光喇曼光谱特征参数反映有机碳质的成熟度

激光喇曼光谱技术具有测试简便和准确等优点，有机碳质物的激光喇曼光谱参数可以作为成熟度指标。并计算出镜质体反射率与激光喇曼光谱参数的回归关系式，喇曼光谱成熟度指标对应的成熟度范围为2.1%-15.0%(Ro)。     

压力对镜质组反射率与烃类生成的影响

镜质组反射率是广泛使用的成熟度参数，长期以来人们普遍认为温度和时间是控制镜质组反射率的主要因素，而压力对反射率的影响人们知之甚少。本文综述了模拟实验地质实例中压力对镜质组反射率和烃类生成影响的不同观点，介绍了过压造成镜质组反射率延尺的两个动力学模型。     

TEM中心回线测深解释新参数△n的数学模型研究

层状大地上中心回线与重叠回线的近区视角电阻率公式中，延时t的幂指数n以往都当作常数，实际是一变数，笔者推导证明它的差值△n与视纵向电导S的对数差分相当，因此可以用于直接划分地层的解释，为此做了相当数量的数学模型计算，结果表明在一定条件下能较好地划分地层，但观测精度要求很高。     

预应力锚索抗滑桩内力计算的有限差分“m-k”法

在于地基系数“m-k”法的原理，提出进行锚索抗滑桩全桩内力计算的有限差分法，同时，编写了该法的计算和图形处理程序。实例计算结果表明：只要等量差分段取得足够小，虽然所耗机时稍长，但可获得高精度数值解，其精度高于采用传统手算法所得的精度，而且图形处理结果可优化抗滑桩的结构设计，并考虑了各种桩底支承条件。本方法基本原理同悬臂桩法，滑动面以下锚索桩的内力计算与普通桩完全相同，滑动面以上不同的只是此时的桩顶条件有所不同，必须考虑锚索拉力的作用，地基系数“m-k”法适用于滑动面以上地层为土层或严重风化岩层、滑动面以下为岩层的情况。