黄建平1,2,,,
李振春1,2,
黄金强3,
李庆洋4,
朱峰1,2,
刘梦丽1,2
1. 中国石油大学(华东)地球科学与技术学院, 山东青岛 266580
2. 海洋国家实验室海洋矿产资源评价与探测技术功能实验室, 山东青岛 266071
3. 贵州大学资源与环境工程学院, 贵阳 550025
4. 中国石化中原物探研究院, 河南濮阳 457001
基金项目: 中国科学院战略性先导科技专项(A)(XDA14010303),国家自然基金重点项目(41720104006),国家油气重大专项(2016ZX05002-005-07HZ,2016ZX05014-001-008HZ),山东省创新工程(2017CXGC1602),山东省青岛市科技计划项目(16-5-1-40-jch)及自主创新科研计划项目(理工科)(17CX05011)联合资助
详细信息
作者简介: 郭旭, 男, 1994, 硕士在读, 主要从事各向异性偏移成像工作.E-mail:guozexu1028@163.com
通讯作者: 黄建平, 1982, 教授, 博士生导师, 主要从事地震波正演及反演成像工作.E-mail:jphuang@upc.edu.cn
中图分类号: P631收稿日期:2018-01-02
修回日期:2018-05-15
上线日期:2019-06-05
Least-squares reverse time migration based on first-order velocity-stress wave equation in VTI media
GUO Xu1,2,,HUANG JianPing1,2,,,
LI ZhenChun1,2,
HUANG JinQiang3,
LI QingYang4,
ZHU Feng1,2,
LIU MengLi1,2
1. School of Geosciences, China University of Petroleum(East China), Qingdao Shandong 266580, China
2. Laboratory for Marine Mineral Resources, Qingdao National Laboratory for Marine Science and Technology, Qingdao Shandong 266071, China
3. Guizhou University Resource and Environmental Engineering College, Guiyang Guizhou 550025, China
4. Sinopec Geophysical Prospecting Institute of Zhongyuan, Puyang Henan 457001, China
More Information
Corresponding author: HUANG JianPing,E-mail:jphuang@upc.edu.cn
MSC: P631--> Received Date: 02 January 2018
Revised Date: 15 May 2018
Available Online: 05 June 2019
摘要
摘要:地下地层普遍存在各向异性,忽略介质各向异性会导致速度估计不准确,成像精度下降.基于二阶声波方程的最小二乘逆时偏移忽略了介质各向异性及密度变化的影响,致使模拟地震数据与实际观测数据不匹配,影响收敛速度和反演成像质量.VTI介质一阶速度-应力方程能较好适应各向异性变密度情况,为此,本文首先从VTI介质一阶速度-应力方程出发,进行波动方程线性化;其次推导了相应的扰动方程和伴随方程,并通过伴随状态法得到梯度更新公式;最终形成基于一阶方程的LSRTM算法理论及实现流程.在实现算法的基础上,通过数值试算及成像结果对比,验证了本文算法在处理变密度和VTI介质时的有效性和优越性.偏移速度以及各向异性Thomsen参数误差的敏感性测试及误差收敛曲线对比结果进一步表明:速度及Thomsen参数对成像结果存在明显影响,其中速度敏感性最强,参数epsilon次之,参数delta的敏感性最弱.
关键词: 一阶速度-应力方程/
VTI介质/
最小二乘逆时偏移/
伴随状态法/
敏感性
Abstract:Subsurface formation is generally anisotropic, which can affect the kinematic characteristics of seismic waves. Ignoring this feature will result in inaccurate velocity estimation, thus lowering the imaging resolution of target regions. The conventional least-squares reverse time migration (LSRTM) algorithm based on the second-order equation ignores the influence of anisotropy and density, thus resulting in mismatch of the simulated and observed seismic data, affecting convergent speed and imaging quality. To make improvement on this issue, this work proposes a new LSRTM algorithm based on the first-order velocity-stress equation in VTI media. First we linearize the wave equations from velocity-stress equation. Then we derive the corresponding perturbation equations and adjoint equations and obtain the updated gradient formula using the adjoint-state method. Finally, the framework and implementation process of the LSRTM method based on first-order velocity-stress equations are formed. We demonstrate the validity this new algorithm by numerical tests. Amount of the sensitivity tests indicate that velocity and Thomsen parameters have obvious influence on the imaging results, in which the velocity sensitivity is highest, next the parameter epsilon, and delta sensitivity is the least.
Key words:First-order velocity-stress equation/
VTI media/
LSRTM/
Adjoint-state method/
Sensitivity
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_15016