地层应力状态下多孔岩石孔隙度校正

多孔岩石孔隙体积是随应力条件而变化的，由于实验室测定岩芯孔隙度时所用的压力条件与储层实际应力不心相同，故实际应用时需进行地层孔隙度校正。本文提出不同应力状态下孔隙度转换公式和方法。     

线性连续变化非弹性介质中震波射线轨迹与走时研究

地球内部非弹性介广泛存在，垂向线性连续变化非弹性介质模型具有重要理论与实际意义。本文利用复速度概念得到了复速度随深度线性连续变化的非弹性介质中地震波射线轨迹与等时线方程。理论结果表明：在线性连续变化非弹性介质中，射线轨迹呈圆形，等时面的形状为椭圆。最后对射线轨迹与等时线方程的适用性条件了讨论。     

关于波浪作用下粉土质海床中共振现象的解释

本文详细讨论了粉土质床在波浪作用下所呈现共振现象的力学机制。利用有限深海床下波浪载荷对土层响应的研究，系统地比较了理论计算结果和实验结果。提出在波浪作用下粉土质土层中存在着一硬土层，从而引发Ｂｒａｇｇ共振。     

3D多孔介质渗透率的格子Boltzmann模拟

多孔介质结构的复杂性使其内部的流动非常复杂,通常不加简化就难以求解。格子Boltzmann方法以其无需简化就能处理复杂孔隙结构的优势,成为研究多孔介质流动的一种有效的数值计算方法。从现场采集的彩色铸体剖面图中提取孔隙信息重构了一种接近真实岩石的3D孔隙介质数字岩心。在此基础上,使用格子Boltzmann方法进一步得到能准确求解多孔介质中流动问题的Navier-Stokes方程,用以研究孔隙岩石的渗流特性。该方法对简单模型渗透率的模拟结果与相应的解析解或实验测量结果吻合良好,对某油田6个样品的模拟结果与实验结果的误差不大,说明格子Boltzrnann方法可以用来计算实际问题的渗透率。     

无限弹性介质中隧道内爆炸动力响应的解析解

基于圆柱形薄壳结构理论,建立爆炸荷载作用下隧道结构的动力平衡方程.采用积分变换法求得Fourier变换域内爆炸荷栽作用下隧道位移场的积分表达式.结合具体算例,采用离散Fourier逆变换的数值计算方法,给出了爆心处隧道位移时程曲线、最大正应力时程曲线及某一时刻隧道位移幅值和最大正应力在x方向的变化曲线.其结果与数值计算结果相一致,验证了方法的可行性.     

气枪激发信号传播的谱元法数值模拟研究

为了研究气枪激发信号的波场,本文利用谱元法对双相介质中波的传播做了数值模拟,分析了波的传播特征.本文主要做了以下工作:(1)研究了使用谱元法(SEM)模拟孔隙弹性介质中波的传播,模拟结果表明,采用谱元法能有效解决双相介质的波场传播模拟问题.(2)验证了Biot理论中慢纵波的存在.双相介质中存在明显的慢纵波,流相波场的慢纵波比固相波场的明显,慢纵波具有很强的散射性质.(3)研究了地震波在孔隙度变化的沉积层中传播规律,验证了孔隙度对于相速度的影响.(4)针对有限的水库底部情况,研究了起伏界面对波传播的影响.     

黏弹性介质中瑞雷波有限差分数值模拟与波场分析

瑞雷面波经常被用来反演地表浅层横波速度,受到越来越广泛的关注。对瑞雷波的研究一般都基于完全弹性介质,而实际地层更接近黏弹性介质,对黏弹性介质中的瑞雷面波进行模拟更具实际意义。本文采用广义标准线性体模型来描述黏弹性介质,并采用交错网格有限差分法对考虑水平自由表面的黏弹介质进行正演模拟,再与弹性介质中的结果进行对比分析。首先采用非线性最优化算法根据期望常数品质因子直接求取松弛时间来拟合常Q模型,并给出广义标准线性固体的具体算例,实施自由表面条件时采用声学-弹性边界近似法,通过剪切模量不变来考虑自由表面上、下横向应力保持连续的条件。对于非自由表面,采用非分裂的多轴卷积完全匹配层来吸收波场。然后对几种典型的数值模型进行正演模拟计算,数值解与解析解的对比验证了本文方法的准确性与有效性,正演结果的对比表明波场尤其是面波频散会受黏弹性影响,因此有必要在面波勘探中考虑黏弹性因素。     

含泥质低孔渗各向异性黏弹性介质中的波频散和衰减研究

本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制（Biot流动和喷射流动机制）有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内（地震波勘探频率）黏弹性力学机制对波传播能量的衰减起主导作用.     

黏弹性与弹性介质中Rayleigh面波特性对比研究

Rayleigh面波的频散特性可以用来研究地表浅层结构. 本文使用时域有限差分法来模拟复杂黏弹性介质中的Rayleigh面波,研究了Q值对面波频散特性的影响.文中采用旋转交错网格有限差分,以非分裂卷积形式的完全匹配层为吸收边界,推出了求解二阶位移-应力各向同性黏弹性波动方程的数值方法.为了检验数值解的精度,首先将简单模型的正演结果与解析解对比,验证了方法的正确性;然后模拟了横向缓变层状介质和含有洞穴的介质中的面波,对弹性和黏弹性介质中的面波的频散特性进行对比分析.模拟结果表明浅层Q值对面波的频散特性有显著的影响;强吸收情况下,高阶面波的能量相对低阶面波能量显著增强.     

基于最大球算法定量表征土体空间孔隙网络

孔隙网络控制着土体渗流、排水固结与基质吸力等重要工程性质。本文介绍了多孔介质孔隙网络最大球建模基本原理与算法。以显微CT扫描振捣干法生成的南京粉砂试样为例,采用最大球算法建立了试样三维重构模型表征单元体(REV)的空间孔隙网络球棍模型,计算得到了样品REV尺度的孔隙网络参数,统计发现,孔隙半径、喉道半径、孔隙配位数、孔隙截面形状因子、喉道截面形状因子与喉道长度等孔隙参数均近似服从正态分布,孔隙体积近似服从衰减型指数分布。孔隙与喉道半径分别分布在100μm与65μm以内,两者数学期望分别为40. 0μm与18. 0μm;配位数分布在25以内,数学期望为5. 1;孔隙与喉道截面形状因子分别分布在0. 01~0. 04与0. 01~0. 05的区间内,两者数学期望分别为0. 019与0. 033;喉道长度分布在100~800μm以内,数学期望为292. 22μm。同时发现,样品中体积小于1. 5×10^7μm^3的小孔隙数量超过90%。本方法可应用于土体细观孔隙结构的定量表征。