基于有限差分的化学剂平面运移规律数值模拟

在化学驱中,各种化学剂在地层中的浓度对于驱油是极为关键的,理想的情况是合理设计化学剂注入参数（注 入浓度、注入量和注入速度）,使得地层各处化学剂的浓度接近实验室设计的最佳浓度,从而驱油效率最高。化学剂在 地层中除正常的对流外,还会发生扩散和吸咐,而平面渗流中,各处渗流速度的不同,使得化学剂在平面的浓度分布比 一维情况复杂得多。依据有限差分法和雅可比迭代方法,根据达西渗流定律确定平面渗流速度场和压力场分布,再根 据已求得的速度场结合化学剂运移的对流弥散方程,求解得到了化学剂的浓度场。数值模拟结果给出了不同时刻化 学剂在地层中的浓度分布,并针对不同的运移滞后系数、弥散系数和化学剂注入量,分析了这些参数对于化学剂分布 的影响。     

高放核废物处置库T-M耦合数值模拟

考虑到地下工程受地应力、温度、水的相互影响,导致岩体力学性质的改变,采用ANSYS有限元分析软件,基于Mohr-Coulomb破坏准则,建立高放核废物处置库的数值模型,在此基础上进行热-力（TM）耦合二维数值模拟．研究了高放核废物处置库9种方案运行一百年的温度、应力、位移变化,通过对比分析,在满足安全、经济的基础上,选择最佳的处置库间距和深度,为核废物处置库的设计提供参考依据．     

基于有限差分方法的水中爆炸三维数值模拟

针对水中爆炸的三维数值模拟问题,采用基于有限差分的拉格朗日与欧拉相结合的方法,在矩形网格上离散差分基本方程组;在欧拉步中引入模糊方法处理界面,计算各输运量;编写了数值模拟程序,进行了炸药在水中爆炸的三维数值模拟.爆炸场中考虑了有刚性障碍物及有反射有限边界计算域.结果表明,模拟结果与实际情况基本吻合,特别是冲击波在水中传播的模拟效果较好,证明所采用的三维数值模拟方法是有效和可行的.     

二维浅水数值模拟的有限差分MADI法

在贴体平面二维正交曲线网格基础上,采用对交替方向隐式ADI法变量位置分布进行了改进的MADI法对二维水深平均的浅水方程进行差分离散求解,并将这一算法应用于长江南通河段的水流模拟.数值模拟计算结果和实测资料比较结果表明,该算法能较好地模拟复杂条件下天然河道涨潮、落潮水流归槽和流速分布过程.     

用时域有限差分方法计算电磁散射

首先介绍了时域有限差分方法的基本原理,然后用该方法计算了球形目标的散射场,并与由ＭＩＥ理论所得到的精确解作了比较,两种结果的很好符合说明了时域有限差分方法的有效性,最后指出了使用时域有限差分方法应该注意的问题     

Euler方程有限差分方法数值模拟

研究了定义在二维水槽上带非线性自由面边界条件的Euler方程的数值解.通过合适的σ坐标变换对不规则的水槽液体区域变换为一个规则的正方形区域,建立流场变量的差分耦合迭代的算法,运用交错网格求解了无粘不可压缩的Euler方程的数值解,数值结果表明,与之前结果和解析解比较数值解较好,对水平激励和垂直激励下非线性的效果和波拍的现象非常明显.     

基于有限差分的隧道建模与分析技术

针对FLAC3D软件在多导洞分步开挖隧道建模方面存在的不足,首次将任意平面前沿推进法三角网生成技术应用于FLAC3D模型建立过程中,重点研究了适合软件内置Fish语言编程的数据结构和算法流程,引入邻边位置判据,提高算法效率;生成正三角形单元或等腰三角形单元,提高单元质量;以数学函数进行疏密控制,有效实现了自适应划分。成功实现了独立的、无需依赖其他建模软件的FLAC3D前处理程序;将网格生成技术应用于实际隧道工程的计算建模中,取得了很好的效果,充分验证了该方法的可行性。     

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

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

新疆花岗岩物理力学性质的矿物组成和微观结构分析

为分析不同类型花岗岩物理力学性质差异,通过选取不同粒径的新疆花岗岩试样进行电镜扫描和矿物成分测试试验,研究岩石矿物组成和微观结构与其物理力学性质之间关系。结果表明：花岗岩中石英和长石矿物胶结很好而具有孔隙率低和强度高特点;随着粒径增大,花岗岩中云母含量增加,导致粗粒花岗岩的孔隙率较高和强度偏低;通过电镜扫描得出花岗岩微观结构质密,矿物胶结很好,颗粒分布均匀,从而具有孔隙率低和强度高特点,同时由于细粒花岗岩微观结构更质密且晶粒镶嵌紧密而孔隙率最低但强度最高。可见天湖细粒花岗岩体由于质密、低孔隙率和高力学强度的特点更适合作为高放废物地质处置岩体。     

基于非稳态对流扩散方程高阶紧致差分的煤自燃数值模拟

为提高求解煤自燃非稳态对流扩散方程的有限差分法精度,构造高阶紧致差分格式,采用Richardson外推法进行求解;并针对煤自燃多场耦合模型的求解复杂度,提出一种具有承袭性的解耦策略。利用本文的解法和策略,对宁夏枣泉矿煤进行煤自燃多场耦合数值模拟,模拟计算的结果与相同实验条件下的实验结果相符。     

混合边界油藏压力分布的有限元法模拟

建立了稳定渗条件下的有限元方程，针对定边界和封闭边界合作用下的渗流问题进行了模拟计算，给出了圆形地层和矩形地层在以上条件下的压力分布，结果表明，混合油藏外边界将产生复杂的地层压力分布，压力分布随地层方位的不同而不同，同时说明，有限元法用于求解复杂边界条件的渗流问题是方便有效的。     

定常对流扩散反应方程的指数型高阶差分格式

提出了一种数值求解一维定常对流扩散反应方程的指数型四阶差分格式．首先,由常数变易法求得模型方程的通解并应用到xi-1,xi,xi＋13点上,得到模型方程的指数型差分格式．然后,利用源项f（x）在点xi处的二阶泰勒展开,得到定常对流扩散反应方程的指数型四阶差分格式．最后,用数值算例验证了该格式的高阶精度和可靠性．     

广义Korteweg-de Vries方程的高精度差分格式

对广义Korteweg-de Vries(generalized Korteweg-de Vries,GKdV)方程的初边值问题进行数值研究,提出一个2层非线性守恒差分格式,该格式的收敛阶为O(τ2+h4)。证明该格式在离散意义下保持原问题质量守恒和能量守恒,分别运用离散能量法和Von Neumann分析法证明该格式的可解性和绝对稳定性。数值实验结果表明,本文格式在时间和空间方向上分别具有2阶和4阶精度,且是质量和能量守恒的。     

非均匀孔隙压力下水压致裂的数值试验

运用数值试验方法研究了非对称分布的孔隙压力梯度下水力压裂裂缝的扩展行为 ,数值试验再现了孔隙水压作用下岩石微破裂诱致宏观破坏的演化过程·数值试验结果表明 ,水力压裂裂纹的扩展不仅受到裂纹尖端局部孔隙水压力大小的影响 ,同时也受控于宏观上孔隙水压梯度分布 ,在扩展方向上致裂裂纹往往偏向于高孔隙压力的区域 ,并且其压裂压力随着局部孔隙压力的增大而降低 ,数值试验研究结果与物理实验结果有较好的一致性     

裂缝多孔介质中达西流动的有限差分方法

用差分方法来模拟二维裂缝多孔介质中的单相达西流动问题.采用降维模型对二维区域内的裂缝进行建模,相比整个区域而言裂缝的宽度很小并且把裂缝视为一个一维的界面,流体通过裂缝会和周围多孔介质发生耦合现象.在裂缝和周围多孔介质中,流体流动均遵循达西定律和守恒定律.采用差分方法来求解降维模型中推导的流体流动方程.通过数值实验验证了该方法的有效性,并证明了裂缝是快速通道还是地质屏障取决于裂缝处渗透率张量的大小.     

随钻地震波场高阶交错网格有限差分数值模拟

 随钻VSP是利用钻头振动噪声作为震源,在地面布设检波器进行信号接收的一种新的物探技术。基于随钻地震波场的空间传播特征,从一阶速度-应力波动方程出发,利用高阶交错网格有限差分,模拟出不同炮点深度和不同偏移距情形的随钻地震波场,研究随钻地震波场的空间传播规律和特征。结果表明,数值模拟可以有效的解决现场试验周期长,费用高,地下结构未知因素多等问题,为随钻地震技术研究提供了一种高效的手段。     

油藏数值模拟技术的误差分析

为了使数值模拟技术成为认识油藏地质的重要手段 ,并成为方案实施的理论基础 ,提高模拟的可信度和真实性 ,必须将误差降低到最小程度 .以吐哈油田陵 2中东区为例 ,根据现场和理论相结合 ,归纳出了数值模拟过程中存在的几种主要误差 ,系统地剖析了产生误差的原因 ,并相应地提出了消除误差的办法 ,进一步完善了数值模拟工作体系 ,提高了模拟精度 ,为定量描述和评价油藏、预测开发趋势和优化技术对策及措施提供了有力的技术保障 ,也对认识油藏和油田开发有着十分重要的意义 .