毕设答辩-粘声波正演模拟技术分析.ppt

* KAUST King Abdullah University of Science and Technology 粘声波正演模拟方法研究 姓 名：xxx 学 号： xxx 专 业： 地球科学与技术学院 导 师： xxx 地球科学与技术学院--本科毕设答辩 一、目的与意义 二、国内外研究现状 三、粘声波方程数值模拟过程 目录 四、模型试算 五、认识与讨论 一、目的与意义 二、国内外研究现状 三、粘声波方程数值模拟过程 目录 四、模型试算 五、认识与讨论 * 一、目的与意义 在地震勘探中，我们通常把地下介质为看成弹性各向同性介质。但是在实际中，地下介质为非完全弹性各向异性介质，存在吸收衰减现象。因此更趋向于粘弹性。粘弹介质正演模拟方法的研究对于还原地下介质的真实情况具有较强的理论意义和实际意义。由于粘弹性介质的正演模拟计算成本较高，计算繁琐，为了降低计算成本，本文进行了声学近似下的粘声波正演模拟方法研究。 * 一、目的与意义 正演模拟方法分类： 正演模拟 积分方程法 数值解法 射线追踪法 有限元法 伪谱法 有限差分法 通过有限差分算子将波动方程离散化，差分代替微分，求解代数方程组得到数值解 可以适应剧烈变化的地下介质，算法速度较快。 一、目的与意义 二、国内外研究现状 三、粘声波方程数值模拟过程 目录 四、模型试算 五、认识与讨论 二、国内外研究现状--粘声粘弹正演 1845年，Stocks在研究过程中，最早提出粘弹性理论，认为地层中非弹性机械能转化成热能的能量损失是因为地层的内摩擦； 2011年，廖建平等人进行了二维频率空间域粘声波正演模拟研究，采用了嵌套剖分网格排序法，以便于节省计算机内存； 2014年，郭成锋、杜启振进行了有关二维VTI介质粘声波正演模拟的研究 2014年，吴玉等人进行了基于分数阶拉普拉斯算子解耦的粘声介质地震正演模拟与逆时偏移的研究 20世纪40年代，N.H.Ricker提出粘弹介质更加接近实际地层，推动了粘弹理论的发展； 1988年，Carcione研究了粘滞声波在地层中传播的正演模拟方法； 1995年R·Knight等人提出了在二维随机介质中粘声波的传播以及吸收和波散衰减； 1995年，How-Wei Chen等人提出了一种并行算法，用来进行粘滞声波数值模拟，并在二维介质上得到了验证； 1996年，Jose M. Carcione等人进行了有关流体饱和粘声波孔隙介质中纵波的全频率瞬时解的研究 一、目的与意义 二、国内外研究现状 三、粘声波方程数值模拟过程 目录 四、模型试算 五、认识与讨论 * 三、粘声波方程数值模拟过程 本文主要进行了二维各向同性介质粘声波方程数值模拟，过程如下： 震源、边界条件、 稳定性、数值频散 编程、模型试算 结果分析 粘声波方程的推导 速度-应力方程的推导 计算差分系数 有限差分格式构建 * 三、粘声波方程数值模拟过程-方程推导 非均匀各向同性介质二维声波方程： 其中： 一阶速度-应力方程： * 三、粘声波方程数值模拟过程-网格选择 贴体网格等 网格 结构网格 非结构网格 规则网格 曲线网格 常规网格 交错网格 三角形网格等 网格分类： 计算量相同时，交错网格数值模比常规网格更精确，数值计算更稳定，收敛速度更快，压制数值频散更彻底，因此使用交错网格进行模拟。 * 三、粘声波方程数值模拟过程-差分格式推导 交错网格差分格式： 差分系数计算矩阵： * 二维粘声波方程正演模拟的交错网格有限差分格式： 三、粘声波方程数值模拟过程-差分格式推导 * 三、粘声波方程数值模拟过程-震源选择 震源的选择： 本文选择的是雷克子波作为加载的震源，来模拟激发地震波。表达式为：