中央高校基本科研业务费专项资金(14CX06072A) 作品数:9 被引量:79 H指数:7 相关作者: 黄建平 李振春 李娜 李娜 张艳艳 更多>> 相关机构: 中国石油大学(华东) 中国石油化工集团公司 华北油田 更多>> 发文基金: 中央高校基本科研业务费专项资金 国家自然科学基金 国家重点基础研究发展计划 更多>> 相关领域: 天文地球 电子电信 更多>>
基于一阶速度-应力方程的多震源最小二乘逆时偏移 被引量:9 2016年 最小二乘逆时偏移(Least-Square Reverse Time Migration,LSRTM)相比于常规偏移具有更高的成像分辨率、振幅保幅性及均衡性等优势,是当前研究的热点之一.然而,目前LSRTM算法大多是基于二阶常密度标量声波方程建立的,忽略了密度变化对振幅的影响,因而基于振幅匹配策略的常规LSRTM很难在变密度介质下取得保真的成像结果.一阶速度-应力方程能够很好地处理变密度介质,但简单地将一阶速度-应力方程应用到LSRTM中缺乏理论基础.为此,本文从LSRTM的正问题入手,提出了基于交错网格的一阶速度-应力方程LSRTM理论方法.首先将一阶波动方程线性化,建立了一阶方程LSRTM的目标泛函,随后推导其伴随方程,并借助伴随状态法给出了迭代更新流程,最终建立了基于一阶速度-应力方程LSRTM的理论框架.进一步,通过在相位编码LSRTM中引入随机最优化思想,极大地减小了计算量、提高了计算效率.最后,通过模型试算验证了本算法的正确性和有效性. 李庆洋 黄建平 李振春 雍鹏 李娜关键词:相位编码 伪深度域声波数值模拟方法及应用 被引量:8 2015年 地震勘探所面对的地下介质中常包括低速体、起伏地表、高陡构造、孔缝单元等,对这些区域进行地震波场模拟时,为保证模拟精度和计算稳定性,网格间距需要取得很小,从而导致基于笛卡尔坐标系的常规有限差分正演模拟方法的局部过采样问题。为此,本文在前人研究成果的基础上,将伪深度思想引入到正演模拟中,利用曲线坐标系下的梯度、散度公式推导出伪深度域一阶速度—应力方程,实现了伪深度域地震正演模拟。同时,鉴于横向过采样问题没有得到有效解决的实际情况,采用自适应变差分算子计算横向空间导数,提出了伪深度域自适应变差分算子算法,综合考虑了纵、横向的空间采样问题,进一步提高了计算效率。最后,将算子扩展到逆时偏移中,实现了伪深度域自适应逆时偏移算法。模型试算结果表明:相比常规正演模拟算法,伪深度域自适应变差分算子算法在保证精度的前提下,计算效率提高了约25%,存储占用降低了近30%。 李庆洋 黄建平 李振春 李娜 曲英铭 苏云关键词:曲线坐标系 正演模拟 逆时偏移 裂缝型储层时空双变正演模拟研究 被引量:3 2014年 本文将深层裂缝型储层视为各向同性背景下的小尺度非均质体,采用时空双变正演模拟算法模拟其地震响应,并将Lanczos滤波算子推广到该算法中,在保证最优的计算效率及内存的前提下,解决了常规变网格算法长时间采样下的不稳定问题.测试结果验证了本文算法对深层微构造的模拟精度及效率.最后将本算法应用到裂缝储层正演模拟中,定量分析了不同裂缝开度、密度和倾角下的地震振幅响应规律,为裂缝型储层的检测和解释提供了可靠的理论指导. 李振春 李庆洋 黄建平 李娜 田坤关键词:正演模拟 裂缝密度 稳定性 起伏地表贴体全交错网格仿真型有限差分正演模拟 被引量:13 2015年 为了精确模拟起伏地表下的波场特征,本文发展了曲线坐标系下基于仿真型有限差分的贴体全交错网格正演模拟算法。本算法主要有三个优势:1采用贴体网格离散起伏地表,不仅消除了阶梯离散问题,还可适应剧烈起伏情况;2将全交错网格引入到曲线坐标系中,避免了标准交错网格的插值误差和同位网格中奇偶失联引起的高频振荡现象,提高了模拟精度,且减小了算法实现的复杂度;3引入仿真型差分算子实施曲线坐标系下的自由边界条件,使边界和内部具有相同的差分精度,从而进一步提高了模拟精度,且与牵引力镜像法相比具备了更可靠的理论基础。模型试算结果表明:本文算法结果和解析解几乎完全吻合,具有较高的模拟精度,能够适应各种复杂模型;具有宽松的频散条件,在相同精度下本文算法所允许的网格间距比同位网格的网格间距更大。 李庆洋 黄建平 李振春 李娜 王成泉 张艳艳关键词:起伏地表 正演模拟 基于Student′s t分布的不依赖子波最小二乘逆时偏移 被引量:13 2017年 最小二乘逆时偏移(Least-Squares Reverse Time Migration,LSRTM)与常规偏移相比具有更高的成像分辨率、振幅保真性及均衡性等优势,是当前研究的热点之一.震源子波的估计直接影响LSRTM结果的好坏,在实际情况下考虑到震源子波的空变特性,其估计十分困难.为了消除子波对LSRTM结果的影响,本文发展了基于卷积目标泛函的不依赖子波LSRTM算法.目标泛函由观测记录卷积模拟记录的参考道以及模拟记录卷积观测记录的参考道组成,由于观测子波和模拟子波在目标泛函的两项中同时存在,从而消除了子波的影响.此外,常用的基于L2范数拟合的LSRTM算法对噪声非常敏感,尤其是当地震数据中含有异常值时,常规LSRTM无法得到满意的结果.Student′s t分布相比L2范数具有更好的稳健性,本文将其推广到不依赖子波LSRTM中,提升了算法的稳健性,最后通过理论模型及实际资料试算验证了算法的有效性和对复杂模型的适应性. 李庆洋 黄建平 李振春关键词:T分布 卷积 去均值归一化互相关最小二乘逆时偏移及其应用 被引量:12 2016年 最小二乘逆时偏移相对于常规逆时偏移具有更高的成像分辨率、振幅保幅性及均衡性等优势,在一定程度上可满足岩性油气藏勘探的需求,是目前研究的热点之一.然而由于实际地下介质的黏滞性和变密度,以及无法准确地估计震源子波等,基于振幅匹配的常规最小二乘逆时偏移算法很难在实际资料处理中取得好的效果.此外,实际数据常包含大量噪声,进一步限制了常规算法的应用.为此,本文通过修改目标泛函,提出了去均值归一化的互相关最小二乘逆时偏移算法,并给出了陆上资料的应用实例.研究表明,归一化策略减弱了算法对子波能量的苛求;互相关算法强调相位匹配,进一步弱化了子波的影响,提升了算法的稳定性和可靠性;去均值思想在处理低信噪比资料方面有较大优势.理论模型和陆上实际资料处理都验证了本算法的有效性和对复杂模型的适应性. 李庆洋 黄建平 李振春 李娜 李闯关键词:相位匹配 子波估计 基于贴体全交错网格的起伏地表正演模拟影响因素 被引量:5 2016年 贴体网格有限差分正演模拟算法不仅能够精确模拟任意起伏地形下的波场特征,且计算效率较高,是一种很有应用前景的处理西部复杂地表问题的方法;然而,目前求解波动方程时常用的同位网格和标准交错网格,在处理贴体网格起伏地表正演模拟时存在诸多问题。为此,将全交错网格引入到曲线坐标系下,避免了标准交错网格的插值误差和同位网格中奇偶失联引起的高频振荡现象,提高了模拟精度,减小了算法实现的复杂度。在自由边界条件实施时,采用牵引力镜像法计算速度分量,速度自由边界条件配合紧致交错差分格式更新应力分量,得到了较好的效果。随后,重点研究了贴体全交错网格正演模拟算法的影响因素,考虑了网格正交性、网格间距和网格拼接等的影响,并取得了如下认识:算法对网格的正交性没有过分要求;网格间距的突变会引起虚假反射的产生;不同类型的网格拼接对模拟结果不会造成明显的影响。 李庆洋 李振春 黄建平 李娜 苏在荣关键词:起伏地表 正演模拟 影响因素 一种稳定的高精度双变网格正演模拟与逆时偏移方法 被引量:20 2014年 双变网格正演模拟方法因其低存储、高效率等优势,在刻画地下小尺度异常体,尤其是超深部微结构方面发挥了越来越重要的作用。变网格正演模拟方法在实现过程中需要考虑两个问题:网格步长变化处的虚假反射和长时间采样下的不稳定性。从理论上推导了虚假反射误差的数学表达式,分析了虚假反射影响因素,并通过引入Lanczos滤波算子,实现了一种基于交错网格的稳定且高精度的双变网格正演模拟算法。模型正演试算结果表明,在超长采样时间迭代测试下本算法依然稳定,且虚假反射也得到明显压制。随后,将双变网格波场延拓算子应用于偏移成像过程,实现了双变网格逆时偏移算法。裂缝储层模型的成像结果表明,双变网格逆时偏移方法能够以较小的内存,较为准确、高效地成像地下微小构造。 李振春 李庆洋 黄建平 李娜 田坤优化的多震源最小二乘逆时偏移 被引量:16 2016年 最小二乘逆时偏移相比于常规偏移具有较大优势,但计算量过大限制了其推广、应用;相位编码技术通过将多炮数据组合成一个超道集,可有效提高计算效率、压制串扰噪声,编码策略种类繁多,但缺少对比分析。为此,在实现多种编码策略最小二乘逆时偏移算法的基础上,通过模型试算,详细对比、分析了四种较为常用的编码策略(振幅编码、极性编码、随机时延编码、平面波编码),并将随机最优化思想推广到相位编码最小二乘逆时偏移中,提出了优化的相位编码最小二乘逆时偏移算法。该算法考虑了梯度的随机特性,即每次迭代的梯度都是之前迭代梯度的指数衰减加权平均,通过合并之前的信息,减小了梯度的随机波动。模型试算证实,优化的相位编码最小二乘逆时偏移算法比传统相位编码最小二乘逆时偏移方法收敛更快,可节省计算成本、提高计算效率。 李庆洋 黄建平 李振春 李娜 李闯关键词:相位编码