SeisWay模型正演.ppt

模 型 正 演 汇报人：侯建全 实现功能 地震模型正演方法介绍 二维建模及软件运行流程 三维建模软件运行流程 核心技术 正演模拟方法 波动方程有限差分 较适用于强非均匀介质模型，任意复杂的地质模型，考虑波传播的动力学特性对地震响应的影响，模拟结果较精确，计算速度要求不高的情况。 射线追踪方法 较适用于近似水平层状介质模型，主要考虑反射界面的地震响应形态，对计算速度要求较高的情况。 射线追踪法原理 费马原理 求解特征方程 高斯射线束方法 高斯射线束 点源可以表示为高斯射线束的叠加 波动方程有限差分法 1）声波方程地震模型正演 2）弹性波方程地震模型正演 1）高阶差分声波模拟中数值频散问题理论分析 高阶差分法声波模拟 2）弹性波波动方程数值模拟 矢量波传播 弹性波方程较声波方程更准确地描述了地震波在固体介质中的传播 岩性勘探的要求 多波多分量勘探的发展 交错网格高阶差分法弹性波模拟 弹性波交错网格高阶差分频散理论分析 实现功能 地震模型正演方法介绍 二维建模及软件运行流程 三维建模软件运行流程 核心技术 模型正演流程 参考SEGY剖面建立地震模型 模型示例二(逆断层) 实现功能 地震模型正演方法介绍 二维建模及软件运行流程 三维建模软件运行流程 核心技术 绿山建模过程 1、设置模型范围 根据工区的实际,可以建立深度模型和时间模型.如下图,可以选择合适的深度单位。 通过手工输入模型范围，也可以通过导入该工区所设计的观测系统确定模型范围；模型范围确定后可以生成初始的三维模型。 2、创建剖面 3、层位编辑 介质定义 根据工区实际地质情况,可一估算出各层横波和纵波速度以及各层的密度.然后,可以对这一属性进行定义。 5、模型显示　 通过OPENGL库,实现对三维模型的放大 、缩小、旋转等操作，可以对地层进行材质和纹理处理，以达到真实模拟地层的目的。 三维模型观测系统加载 显示射线路径 开发三维模型正演模块 在Linux系统下开发并行的三维射线追踪方法、波动方程方法（声波、弹性波）的正演模拟功能； 进一步完善三维建模功能，并实现与多种软件系统输出数据的接口，丰富模型功能； 完善照明度分析方法。通过照明度分析、CRP面元分析，指导观测系统设计，提高软件系统技术含量。 Linux系统二维模型正演 实现功能 地震模型正演方法介绍 二维建模及软件运行流程 三维建模软件运行流程 核心技术 LandMark模型显示 关键技术 实现了二维、三维的声波、弹性波传播波动方程数值模拟，采用了交错网格的高阶差分方法（2－10阶）。 开发重点 射线追踪技术的发展 波动方程法模型示例一 水平层状模型 正演模拟炮记录 模型示例二 逆断层模型 正演模拟炮记录 射线追踪模拟记录 射线追踪路径 射线追踪法模型一 输入测线 编辑地层线 网格化地层生成2D 剖面 输入层的属性 根据 2D 剖面生成网格化界面 参考 SEG---Y 地震数据来拾取定义层位 根据网格化后的层位文本文件直接输入 可从几个业界标准(Geoquest, LandMark…)输入界面模型 生层三维模型 同类软件三维建模介绍 三维地震建模流程 显示初始模型 4、模型网格化 创建三维地震模型是根据有限的控制点数据通过一定的网格化算法,计算出符合模型要求的地层数据。经过对多种网格化算法优缺点的比较，选择了四种网格化算法，根根据地层特征的不同，选择不同网格化算法。 克里格 反距离 多元二次函数 Delaunay三角剖分 覆 盖 次 数 对 比 CMP CRP CRP面元信息分析 偏移距－覆盖次数 方位角 设计横纵方向的均匀网格，在深度方向上计算模型节点上相应的地震属性值，便于同处理软件的结合，运算速度快 多数软件采用的建模方法 绿山软件早期采用的建模方法 适合任意模型 * 汇报人：侯建全 2006－09 = 稳定值 其中 有限差分模拟遇到的主要问题： 模拟精度以及数值频散问题 人为边界反射问题 稳定性 计算效率问题 将平面谐波代入二维声波差分方程，可以得到相速度空间频散曲线（设 ）： 影响因素： 一个波长内离散点数、差分精度、传播方向 = 高阶差分系数的确定： 传统2阶 2M阶 数值频散试验dx=dz=10m, dt=1ms，高阶差分效果好 不同差分精度时间频散曲线 例子：V=3000m/s，f=50hz，则波长为60m，要保证空间频散较低，dx取10m。但稳定性要求 ，即 ，若取dt=1.8ms，则dt/T=0.09。实际模拟中的频散主要由空间离散引起。 不同差分精度空间频散曲线 算法稳定性问题 稳定性问题是数值求解波动