学位论文 > 优秀研究生学位论文题录展示

三维合成震源记录交互剩余偏移速度分析

作　者: 杨其强
导　师: 冯恩民；张叔伦
学　校: 大连理工大学
专　业: 运筹学与控制论
关键词: 剩余偏移速度分析合成震源记录叠前深度偏移控制照明正演模拟傅里叶变换哈特莱变换
分类号: P631.4
类　型: 博士论文
年　份: 2009年
下　载: 133次
引　用: 0次
阅　读: 论文下载

内容摘要

本论文围绕三维合成震源记录交互剩余偏移速度分析课题开展研究并进行了拓展.内容包括地震波数值模拟,以提高计算效率为目的的合成震源记录叠前深度偏移算法和叠后偏移算法,以及利用控制照明的交互剩余偏移速度分析.地震数值模拟是在假定地下介质结构模型和相应物理参数己知的情况下,模拟研究地震波在地下各种介质中的传播规律,并计算在地面或地下各观测点所应观测到的数值地震记录中的一种地震模拟方法,它在地震勘探和天然地震领域中具有广泛的运用.在地震数据处理中,地震数值模拟计算出来的记录数据可以检验各种反演方法的正确性.波场数值模拟技术的发展主要体现在两个方面,即模拟精度和运算速度.傅里叶变换方法正演模拟由于具有较高的计算效率和精度而被广泛运用在地震波数值模拟中.本文对利用傅里叶变换的地震波数值模拟进行了改进和推广,并提出用哈特莱变换取代傅里叶变换的方法,以便利用当前的计算机资源有效快速地模拟地震波传播;尽管如此,用数值模拟方法提供三维合成数据以验证一些地震处理方法(如本文中的三维剩余偏移速度分析方法)仍因计算量庞大、计算时间长而有困难,为此,我们提出了合成震源正演模拟和线形震源正演模拟,合成震源正演模拟可以有效地提供合成震源记录以验证一些涉及合成震源记录相关的地震处理方法,而且不需要耗时较长的合成过程,利用线形震源正演模拟可以将合成数据用于需要利用控制照明的地震处理方法中.由于哈特莱变换具有比傅里叶变换效率高的特点,我们还把它引进到零炮检距正演模拟和偏移中.波动方程叠前深度偏移具有成像精确的特点,但速度较慢影响了其应用,特别是三维偏移.合成震源记录叠前深度偏移基于面炮合成技术减少了需要偏移的地震记录,提高了计算效率.在此基础上,我们利用大步长深度延拓和相位移波场插值,并与稀疏频率方法组合使用,大大减少了计算量,对二维盐丘模型、二维Marmousi模型、二维实际资料和三维盐丘模型试算表明,此方法计算时间可缩减一半以上而不影响成像质量.基于合成震源记录偏移的三维剩余偏移速度分析方法,建立在共成像点道集拉平的剩余曲率分析思想之上,并引进了控制照明技术.该方法使用大步长延拓与稀疏频率方法结合的合成震源记录叠前深度偏移方法,成像精确且具有较高的效率.偏移后可以从结果中直接抽取共成像点道集进行剩余曲率分析.每次迭代可对所有层位进行校正而不需要层剥离过程,这样可以减少迭代次数从而提高效率.控制照明技术的利用消除了平面波在非均匀介质中传播而发生的畸变,保证了剩余速度校正公式的平面波假设条件.三维合成数据和陆上实际资料的试算验证了该方法的有效性和实用性.并基于Qt编程技术和东方地球物理石油勘探公司开发的底层库编制了偏移速度分析软件包.

全文目录

摘要  4-6
Abstract  6-11
1 绪论  11-24
  1.1 地震波场数值模拟发展概况  11-15
    1.1.1 射线追踪技术  11-12
    1.1.2 波动方程数值模拟技术  12-15
  1.2 偏移发展概况  15-20
    1.2.1 偏移方法分类  16
    1.2.2 二维偏移和三维偏移  16
    1.2.3 叠前偏移和叠后偏移  16-17
    1.2.4 时间偏移和深度偏移  17
    1.2.5 积分类偏移方法和基于波场延拓的偏移方法  17-19
    1.2.6 波动方程地震偏移成像的发展趋势  19-20
  1.3 偏移速度分析发展概况  20-22
    1.3.1 深度聚焦分析方法  20-21
    1.3.2 共聚焦速度分析  21
    1.3.3 剩余曲率分析  21-22
  1.4 本论文的工作  22-24
2 地震波数值正演模拟  24-50
  2.1 一阶和二阶导数计算  24-30
    2.1.1 有限差分方法  24-25
    2.1.2 傅立叶变换方法  25-26
    2.1.3 哈特莱变换方法  26-27
    2.1.4 速度和精度比较  27-29
    2.1.5 利用常规傅里叶变换性质的算法  29-30
  2.2 声波方程正演模拟  30-34
    2.2.1 声波方程的数值解  30-31
    2.2.2 震源函数  31-32
    2.2.3 吸收边界  32
    2.2.4 精度  32-34
  2.3 合成震源声波方程正演模拟  34-36
    2.3.1 平面波震源正演模拟  34-35
    2.3.2 线形震源正演模拟与面炮记录合成  35-36
  2.4 零炮检距正演模拟  36-37
    2.4.1 相位移法  36
    2.4.2 上行波  36-37
    2.4.3 双程非反射方程  37
  2.5 数值试算  37-50
    2.5.1 叠前正演模拟  37-38
    2.5.2 零炮检距正演模拟  38
    2.5.3 时间比较  38-50
3 地震偏移  50-68
  3.1 叠前偏移  50-57
    3.1.1 大步长延拓与相位移波场插值  50-51
    3.1.2 大步长延拓与稀疏频率结合  51
    3.1.3 数值试算  51-57
  3.2 叠后偏移  57-68
    3.2.1 二维哈特莱变换  58
    3.2.2 哈特莱有限差分偏移  58-60
    3.2.3 逆时偏移  60-62
    3.2.4 最小二乘哈特莱有限差分偏移  62-63
    3.2.5 数值试算  63-68
4 三维交互剩余偏移速度分析  68-92
  4.1 深度-射线参数域的剩余校正公式  69-73
    4.1.1 二维剩余偏移速度分析  69-72
    4.1.2 深度-射线参数域的三维剩余校正公式  72-73
  4.2 利用控制照明的速度分析  73-75
  4.3 合成数据及野外实际资料试算  75-87
    4.3.1 合成数据  77-78
    4.3.2 陆上实际资料  78-87
  4.4 交互界面软件  87-92
结论  92-93
参考文献  93-104
攻读博士学位期间发表学术论文情况  104-106
致谢  106-107
作者简介  107-108

三维合成震源记录交互剩余偏移速度分析

内容摘要

全文目录

相似论文