多分量感应测井响应的交错网格有限差分法模拟

从电磁响应方程出发,应用交错网格有限差分法推导了三维各向异性地层多分量感应测井响应的差分计算格式,并利用所形成的算法考察了多种环境条件下多分量感应测井的响应规律。计算结果表明,多分量感应测井各磁场分量不仅与地层水平电阻率和垂直电阻率有关,还与仪器方位角、井斜角、井眼轴与地层走向的相对位置有关。根据仪器方位角与井斜角间的关系和各测量量,可反演出地层水平电阻率和垂直电阻率,进而改善电性各向异性储层的含油饱和度评价。计算结果表明该算法是准确可靠的。     

基于三维有限差分方法的三分量感应测井正演模拟

在三分量感应测井中,正演模拟是研究测井响应特征和资料处理的必要基础,其能够为快速有效分析地层物理特性提供依据,从而帮助找到最佳钻井区域。为系统考察多种复杂地层模型下的三分量感应测井的响应特征,从Maxwell电磁方程组出发,利用交错网格有限差分法推导三维频域电磁响应的差分计算格式。针对水平地层和倾斜地层等复杂模型中的三分量感应测井响应进行数值计算,总结发现三维感应测井工具对复杂地层具有很强的探测能力,且纵向磁场分量比水平磁场分量具有更高的纵向分辨率,对于地层变化更加敏感。正演模拟及分析结果可为三分量感应测井仪器参数的优化设计、资料处理等工作提供重要的理论基础及参考依据。     

非均匀介质地震波传播交错网格高阶有限差分法模拟

采用常规的二阶声波方程有限差分方法对于非均匀介质进行了数值模拟时，其数值模拟精度较低。而采用一阶双曲型标量波动方程，则无须对介质的弹性常数进行空间求导。根据Taylor级数展开式，推导出了交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式以及一阶双曲型标量波动方程交错网格任意偶数阶精度差分格式，并给出了该差分算法的稳定性条件。用该差分算法对均匀介质模型、非均匀介质模型和Mannousi模型进行了数值模拟试验，并与伪谱法进行了对比。结果表明，一阶双曲型标量波动方程交错网格高阶差分法的模拟精度与伪谱法的精度非常接近，计算效率高，且适合于模拟非均匀介质、复杂构造和复杂地质体的地震波场。     

三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟

三维数值模拟是研究各向异性介质中复杂弹性波波场特征和传播规律的重要手段。根据Taylor展开式,推导出了交错网格一阶空间导数的任意偶数阶精度有限差分近似式和相应的差分系数计算式,给出了三维各向异性介质中弹性波一阶双曲型应力速度方程交错网格任意偶数阶精度差分格式和稳定性条件,并推导出了三维各向异性介质PML法吸收层边界条件公式和相应的交错网格差分格式。对方位各向异性介质模型和正交各向异性介质模型中弹性波的传播进行的三维数值模拟结果表明,弹性波在三维各向异性介质中传播时存在拟P波、拟SV波和拟SH波,并出现横波分裂、横波分裂盲区、波面三分叉等特殊现象,另外,弹性波场在空间是变化的,其拟P波、拟SV波和拟SH波的耦合关系比较复杂。     

非均匀介质地震波传播交错网格高阶有限差分法模拟

采用常规的二阶声波方程有限差分方法对于非均匀介质进行了数值模拟时 ,其数值模拟精度较低。而采用一阶双曲型标量波动方程 ,则无须对介质的弹性常数进行空间求导。根据Taylor级数展开式 ,推导出了交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式以及一阶双曲型标量波动方程交错网格任意偶数阶精度差分格式 ,并给出了该差分算法的稳定性条件。用该差分算法对均匀介质模型、非均匀介质模型和Marmousi模型进行了数值模拟试验 ,并与伪谱法进行了对比。结果表明 ,一阶双曲型标量波动方程交错网格高阶差分法的模拟精度与伪谱法的精度非常接近 ,计算效率高 ,且适合于模拟非均匀介质、复杂构造和复杂地质体的地震波场     

三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟

基于Biot理论，给出了三维双相各向异性介质应力一速度弹性波方程交错网格任意偶阶精度有限差分解法，并对三维双相各向异性（TIH）介质中弹性波场进行了模拟。结果表明，弹性波在三维双相各向异性介质中传播时存在快纵波qP1、快横波qS1、慢横波qS2和慢纵波qP2，并清楚地观测到。T横波分裂、横波分裂盲点、波面三分叉等特殊现象。快纵波在固相和流相中的相位相同；而慢纵波在固相和流相中的相位相反，即慢纵波在流相中振幅大，而在固相中的振幅较小。另外，三维双相各向异性介质中弹性波场的快纵波和快横波的耦合关系、波的类型、能量分布和相位等都是在三维空间中变化的。     

三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟

基于Biot理论,给出了三维双相各向异性介质应力-速度弹性波方程交错网格任意偶阶精度有限差分解法,并对三维双相各向异性(TIH)介质中弹性波场进行了模拟.结果表明,弹性波在三维双相各向异性介质中传播时存在快纵波qP1、快横波qS1、慢横波qS2和慢纵波qP2,并清楚地观测到了横波分裂、横波分裂盲点、波面三分叉等特殊现象.快纵波在固相和流相中的相位相同;而慢纵波在固相和流相中的相位相反,即慢纵波在流相中振幅大,而在固相中的振幅较小.另外,三维双相各向异性介质中弹性波场的快纵波和快横波的耦合关系、波的类型、能量分布和相位等都是在三维空间中变化的.     

多媒质区域介质子网格时域有限差分法

介绍了应用于多媒质区域电磁问题的介质子网格时域有限差分法．该方法在保证精度的前提下，节省了大量计算时间及计算机内存，尤其在处理多媒质界面方面具有很大的优势．     

用柱状成层各向异性介质的并矢Green函数模拟多分量感应测井仪器的响应

采用递推方法得到柱状成层各向异性介质(横向各向同性)中并矢Green函数的解析表达式。该表达式可用于模拟柱状成层各向异性地层中任意点源(包括电流源和磁流源)的响应,地层数目可以任意,源点和场点的位置可以在任意地层中。利用上述表达式模拟含金属心轴和绝缘保护层多分量感应测井仪器在有井眼和侵入带各向异性地层中的响应。为提高模拟精度,考虑各分量线圈系的具体形状。将金属心轴作为一层介质处理,既可以考虑其电导率有限,也可以考虑其电导率为无穷大的情况。数值模拟结果表明,共面线圈系具有与共轴线圈系完全不同的响应特性。共面线圈系的响应特性更为复杂,对钻井液电导率、侵入带电导率、地层电导率、地层各向异性的变化更为敏感,且在很多情况下其响应会随这些参数的变化出现符号改变。此外,由于共面线圈系的同一响应可对应各向同性地层或各向异性地层的不同电导率,从而使得对测量数据的解释处理变得复杂。     

三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟

三维数值模拟是研究各向异性介质中复杂弹性波波场特征和传播规律的重要手段。根据Taylor展开式，推导出了交错网格一阶空间导数的任意偶数阶精度有限差分近似式和相应的差分系数计算式，给出了三维各向异性介质中弹性波一阶双曲型应力-速度方程交错网格任意偶数阶精度差分格式和稳定性条件，并推导出了三维各向异性介质PML法吸收层边界条件公式和相应的交错网格差分格式。对方位各向异性介质模型和正交各向异性介质模型中弹性波的传播进行的三维数值模拟结果表明，弹性波在三维各向异性介质中传播时存在拟P波、拟SV波和拟SH波，并出现横波分裂、横波分裂盲区、波面三分叉等特殊现象，另外，弹性波场在空间是变化的，其拟P波、拟SV波和拟SH波的耦合关系比较复杂。     

变步长的有限差分线法

在等步长有限差分线法的理论基础上,构造出变步长的有限差分线法格式,并对其进行了误差分析。     

用有限差分法分析列车对隧道中电磁波传播特性的影响

文章推导出复杂边界条件下二维Helmholtz方程的差分公式，实例分析了列车对矩形隧道和半圆形隧道中电磁波截止频率的影响、列车-矩形隧道中主模的衰减常数。实例计算结果表明：有限差分法的计算结果与有限元法计算结果完全一致，有限差分法是分析列车-隧道中电磁波传播问题的一种有效方法。     

对流扩散方程的非一致网格有限差分方法

构造了一种二阶非等距网格差分格式,给出了截断误差及其稳定性.数值算例给出了几种不同网格处理情形下的计算结果,和已有的差分格式进行比较,表明新格式具有较好的平均误差分布.     

波导截止特性和衰减常数的有限差分分析

文章提出用有限差分法计算波导截止频率、截止波数、衰减常数,推导出用有限差分法计算相关参数的计算公式,编程计算了矩行波导、L型波导和单脊型波导的TE、TM波的截止频率、截止波数及衰减常数,计算结果与解析解非常吻合,证明文章提出的计算方法是正确和可行的。     

粗网有限差分法及其燃耗计算

本文介绍的粗网有限差分法求解二维双群中子扩散方程，可用于反应堆燃料管理中的燃耗计算，计算结果与通用的细网有限差分法、节块格林函数方法作了比较，其有效增殖系数的相对误差小于０．０８％，该方法具有计算时间短、满足一定精度要求等特点，适用于工程初步设计的燃料管理分析。