Please wait a minute...
E-mail Alert Rss
 
物探与化探  2018, Vol. 42 Issue (5): 1013-1025    DOI: 10.11720/wtyht.2018.1559
     方法研究·信息处理·仪器研制 本期目录 | 过刊浏览 | 高级检索 |
三维黏声最小二乘逆时偏移方法模型研究
李金丽1,2, 曲英铭3,4, 刘建勋1,2, 李振春3, 王凯1,2, 于博5
1. 中国地质科学院 地球物理地球化学勘查研究所,河北 廊坊 065000
2. 国家现代地质勘查技术研究中心,河北 廊坊 065000
3. 中国石油大学(华东) 地球科学与技术学院,山东 青岛 266580
4. 中国石化地球物理重点实验室,江苏 南京 211030
5. 中国石油天然气集团公司 北京油气调控中心,北京 100007
A model study of three-dimensional viscoacoustic least-squares reverse time migration
Jin-Li LI1,2, Ying-Ming QU3,4, Jian-Xun LIU1,2, Zhen-Chun LI3, Kai WANG1,2, Bo YU5
1. Institute of Geophysical and Geochemical Exploration,Chinese Academy of Geological Sciences,Langfang 065000,China
2. National Center for Geological Exploration Technology,Langfang 065000,China
3. School of Geosciences,China University of Petroleum (East China),Qingdao 266580,China
4. SINOPEC Key Laboratory of Geophysics,Nanjing 211030,China
5. PetroChina Oil and Gas Pipeline Control Center,China National Petroleum Corporation,Beijing 100007,China;
全文: PDF(6805 KB)   HTML
输出: BibTeX | EndNote (RIS)      
摘要 

针对黏介质成像存在的问题及最小二乘逆时偏移方法的优势,通过广义标准线性固体(GSLS)的三维黏声波动方程,基于三维黏声波动方程的伴随算子及最小二乘逆时偏移框架,实现了三维黏声最小二乘逆时偏移方法。采用极性编码技术大幅度降低黏声最小二乘逆时偏移算法的计算量与内存,使三维黏声最小二乘逆时偏移算法的实用化成为可能。在模型试算部分,首先通过Marmousi模型验证了该方法对传统声波最小二乘逆时偏移方法的优势,最后对三维含Q平层模型及盐丘模型进行试算,证明了三维黏声最小二乘逆时偏移算法的正确性。

服务
把本文推荐给朋友
加入引用管理器
E-mail Alert
RSS
作者相关文章
李金丽
曲英铭
刘建勋
李振春
王凯
于博
关键词 三维黏声介质最小二乘逆时偏移伴随算子相位编码    
Abstract

In view of the problem of viscous imaging and the advantage of least-squares reverse time migration (LSRTM),a 3D viscoacoustic LSRTM method is proposed by using the viscoacoustic wave equation of the General Standard Linear Solid (GSLS),which is based on the adjoint operator of the 3D viscoacoustic wave equation and the frame of LSRTM.The polar encoding technology is used to dramatically decrease the computational cost and memory storage of the 3D viscoacoustic LSRTM,which can make it possible to be practical.In the numerical examples,the authors firstly verify the advantage of the proposed method over the traditional acoustic LSRTM using the Marmousi model.Finally,the test of the 3D flat model and salt dome model with Q proves the accuracy of the proposed 3D viscoacoustic LSRTM.

Key words3D    viscoacoustic media    least-squares reverse time migration    adjoint operator    encoding
收稿日期: 2017-12-12      出版日期: 2018-10-24
:  P631.4  
基金资助:中国地质调查局地质调查项目(DD20160224);中国地质调查局地质调查项目(DD20160164);中央级公益性科研院所基本科研业务费专项资金资助项目(AS2017J11)
作者简介: 李金丽(1989-),女,硕士,主要从事地震数据处理及偏移成像方面的研究工作
引用本文:   
李金丽, 曲英铭, 刘建勋, 李振春, 王凯, 于博. 三维黏声最小二乘逆时偏移方法模型研究[J]. 物探与化探, 2018, 42(5): 1013-1025.
Jin-Li LI, Ying-Ming QU, Jian-Xun LIU, Zhen-Chun LI, Kai WANG, Bo YU. A model study of three-dimensional viscoacoustic least-squares reverse time migration. Geophysical and Geochemical Exploration, 2018, 42(5): 1013-1025.
链接本文:  
https://www.wutanyuhuatan.com/CN/10.11720/wtyht.2018.1559      或      https://www.wutanyuhuatan.com/CN/Y2018/V42/I5/1013
  三维黏滞声介质交错网格示意
网格点 1 2 3 4
黏声波波场分
量与对应参数
p, M, τP u, ρ-1 l, ρ-1 w, ρ-1
  黏滞声波波场分量和黏弹性参数的空间位置
  黏声最小二乘逆时偏移成像技术路线流程
  Marmousi模型
a—速度模型;b—Q模型;c—偏移速度模型;d—慢度扰动模型
  偏移成像结果
a—声波LSRTM第1次迭代;b—黏声LSRTM第1次迭代;c—声波LSRTM第30次迭代;d—黏声波LSRTM第30次迭代
  x=2 760 m处反射系数与成像结果单道对比
  黏声最小二乘逆时偏移成像结果
a—本文方法;b—低通滤波方法
  不同误差的偏移速度模型
a—误差=3.84%;b—误差=5.63%;c—误差=7.52%;d—误差=10.13%
  不同误差的迭代收敛曲线
  平层模型
a—速度模型;b—慢度扰动模型
  平层模型不同迭代次数下黏声最小二乘逆时偏移成像结果
a—第1次迭代;b—第20次迭代
  盐丘模型偏移成像结果
a—盐丘模型(从三维盐丘模型中截取);b—反射系数;c—第1次迭代成像结果;d—第20次迭代成像结果
  盐丘等深度切片模型偏移成像结果
a—盐丘等深度切片模型(从三维盐丘模型中截取);b—反射系数;c—第1次迭代成像结果;d—第20次迭代成像结果
[1] Carcione J M . Wave propagating in anisotropic linear viscoelastic media:theory and simulated wave-fields[J]. Geophys. J. In., 1990,101(3):739-742.
doi: 10.1111/gji.1990.101.issue-3
[2] 李振春, 王清振 . 地震波衰减机理及能量补偿研究综述[J]. 地球物理学进展, 2007,22(4):1147-1152.
[3] 李金丽, 李振春, 管路平 , 等. 地震波衰减及补偿方法[J]. 物探与化探, 2015,39(3):456-465.
doi: 10.11720/wtyht.2015.3.04
[4] 刘洋, 李炳秀, 刘财 , 等. 局部时频变换域地震波吸收衰减补偿方法[J]. 吉林大学学报:地球科学版, 2016,46(2):594-602.
[5] Bickel S H, Natarajan R R . Plane-wave Q deconvolution[J]. Geophysics, 1985,50(9):1426-1439.
doi: 10.1190/1.1442011
[6] Gelius L J . Inverse Q-filtering:a spectral balancing technique[J]. Geophysical Prospecting, 1987,35(6):656-667.
doi: 10.1111/gpr.1987.35.issue-6
[7] Varela C L, Rosa A L R,Ulrych T J .Modeling of attenuation and dispersion[J]. Geophysics, 1993,58(8):1167-1173.
doi: 10.1190/1.1443500
[8] Zhang X W, Han L G, Zhang F J , et al. An inverse Q-filter algorithm based on stable wavefield continuation[J]. Applied Geophysics, 2007,4(4):263-270.
doi: 10.1007/s11770-007-0040-9
[9] Hargreaves N D . Similarity and the inverse Q filter:some simple algorithms for inverse Q filtering[J]. Geophysics, 1992,57(7):944-947.
doi: 10.1190/1.1443307
[10] Wang Y H . A stable and efficient approach of inverse Q Filtering[J]. Geophysics, 2002,67(2):657-664.
doi: 10.1190/1.1468627
[11] Wang Y H . Inverse Q-filter migration[J]. Geophysics, 2008,73(1):S1-S6.
doi: 10.1190/gpysa7.73.1so_1
[12] 裴江云, 陈树民 . 近地表Q值求取及振幅补偿[J]. 地球物理学进展, 2001,16(4):18-22.
doi:
[13] 高军, 凌云 , 等. 时频域球面发散和吸收补偿[J]. 石油地球物理勘探, 1996,31(6):865-866.
[14] Bickel S H, Natarajan R R . Plane-wave Q deconvolution[J]. Geophysics, 1985,50(9):1426-1439.
doi: 10.1190/1.1442011
[15] 姚振兴, 高星, 李维新 . 用于深度域地震剖面衰减与频散补偿的反Q滤波方法[J]. 地球物理学报, 2003,46(2):229-233.
[16] Yan H Y, Liu Y . Estimation of inverse Q filtering for prestack reflected PP and converted PS waves[J]. Applied Geophysics, 2009,6(1):59-69.
doi: 10.1007/s11770-009-0009-y
[17] Dai N, West G F. Inverse Q Migration [C]//Expanded Abstracts of the 64 th Annual SEG Meeting,Society of Exploration Geophysicists , 1994: 1418-1421.
[18] Mittet R, Sollie R . Prestack depth migration with compensation for absorption and dispersion[J]. Geophysics, 1995,60(5):1485-1494.
doi: 10.1190/1.1443882
[19] Deng F, McMechan G A .True-amplitude prestack depth migration[J]. Geophysics, 2007,72(3):S155-S166.
doi: 10.1190/1.2714334
[20] Zhang Y, Zhang P, Zhang H Z. Compensating for visco-acoustic effects in reverse-time migration [C]//Expanded Abstracts of the 80 th Annual SEG Meeting,Society of Exploration Geophysicists , 2010: 3160-3164.
[21] Fletcher R P, Nichols D, Cavalca M. Wavepath-consistent effective Q estimation for Q compensated reserve-time migration [C]//Expanded Abstracts of 74 th Annual Internat EAGE Mtg , 2012.
[22] Qu Y, Li Z, Huang J W, et al. The least-squares reverse time migration for viscoacoustic medium based on a stable reverse-time propagator [C]//Expanded Abstracts of the 85 th Annual SEG Meeting,Society of Exploration Geophysicists , 2015a: 3977-3980.
[23] Qu Y, Li Z, Huang J, et al. The application of pseudospectral method and a stable reverse-time propagator for viscoacoustic RTM [C]//Expanded Abstracts of the 85 th Annual SEG Meeting,Society of Exploration Geophysicists , 2015b: 3977-3980.
[24] Qu Y, Huang J, Li Z, et al. A new reverse time migration method for viscoacoustic VTI medium. [C]//77 th EAGE Conference and Exhibition,Expanded Abstracts , 2015c.
[25] Qu Y, Li Z, Huang J , et al. Viscoacoustic anisotropic full waveform inversion[J]. Journal of Applied Geophysics, 2017,136:484-497.
doi: 10.1016/j.jappgeo.2016.12.001
[26] Qu Y, Huang J, Li Z , et al. Attenuation compensation in anisotropic least-squares reverse time migration[J]. Geophysics, 2017,82(6):S411-S423.
doi: 10.1190/geo2016-0677.1
[27] Zhu T Z, Harris J M, Biondi B . Q-compensated reverse-time migration[J]. Geophysics, 2014,73(9):S77-S87.
[28] Zhu T Z, Harris J M, Biondi B . Q-compensated reverse-time migration[J]. Geophysics, 2014,73(9):S77-S87.
[29] 李振春, 孙小东, 刘洪 . 复杂地表条件下共反射面元(CRS)叠加方法研究[J]. 地球物理学报, 2006,49(6):1794-1801.
[30] 李振春 . 地震叠前成像理论与方法[M]. 东营: 中国石油大学出版社, 2011.
[31] 李振春 . 地震偏移成像技术研究现状与发展趋势[J]. 石油地球物理勘探, 2014,49(1):1-21.
[32] 曲英铭, 黄建平, 李振春 , 等. 分层坐标变换法起伏自由地表弹性波叠前逆时偏移[J]. 地球物理学报, 2015,58(8):2896-2911.
[33] 李振春, 李庆洋, 黄建平 , 等. 裂缝型储层时空双变正演模拟研究[J]. 地球物理学进展, 2014,29(3):1183-1188.
doi: 10.6038/pg20140324
[34] 曲英铭, 黄建平, 李振春 , 等. 一种基于非规则网格的地震波射线追踪方法[J]. 石油物探, 2014,53(6):627-632.
[35] 曲英铭, 黄建平, 李振春 , 等. 基于单元交错网格的变坐标系正演模拟方法在声—弹介质中的应用[J]. 石油物探, 2015,54(5):582-591.
[36] 刘洋, 李向阳, 陈双全 . 高精度双吸收边界条件在地震波场模拟中的应用(英文)[J]. Applied Geophysics, 2015,12(1):111-119.
doi: 10.1007/s11770-014-0463-z
[37] 刘洋 . 波动方程时空域有限差分数值解及吸收边界条件研究进展[J]. 石油地球物理勘探, 2014,49(1):35-46.
[38] 刘洋, 李承楚, 牟永光 . 任意偶数阶精度有限差分法数值模拟[J]. 石油地球物理勘探, 1998,33(1):1-10.
[39] 赵岩, 刘洋, 任志明 . 基于优化时空域高阶有限差分方法的粘滞声波叠前逆时偏移(英文)[J]. Applied Geophysics, 2014,11(1):50-62.
doi: 10.1007/s11770-014-0414-8
[40] 陈可洋 . 地震波逆时偏移方法研究综述[J]. 勘探地球物理进展, 2010,33(3):153-159.
[41] 陈可洋 . 基于高阶有限差分的波动方程叠前逆时偏移方法[J]. 石油物探, 2009,48(5):475-478.
[42] 邓文志, 李振春, 王延光 , 等. 基于稳定逆时传播算子的黏声介质最小二乘逆时偏移[J]. 物探与化探, 2015,39(4):791-796.
doi: 10.11720/wtyht.2015.4.22
[43] 刘玉金, 李振春 . 扩展成像条件下的最小二乘逆时偏移[J]. 地球物理学报, 2015,58(10):3771-3782.
[44] Romero L A, Ghigliaz D C, Ober C C , et al. Phase encoding of shot records in prestack migration[J].Geophysics, 65(2):426-436.
[45] Schuster G T, Wang X, Huang Y , et al. Theory of multisource crosstalk reduction by phase-encoded statics[J]. Geophys. J. Int., 2011,184:1289-1303.
doi: 10.1111/gji.2011.184.issue-3
[46] Huang Y S, Schuster G T . Multisource least-squares migration of marine streamer and land data with frequency-division encoding[J]. Geophysical Prospecting, 2012,60:663-680.
doi: 10.1111/gpr.2012.60.issue-4
[47] 陈可洋, 吴清岭, 李来林 , 等. 松辽盆地三维地震资料连片处理关键技术及其应用效果分析[J]. 岩性油气藏, 2012,24(2):87-91.
[1] 郭培虹, 冯治汉, 王万银, 唐小平, 刘生荣. 北秦岭华阳川地区重磁三维反演及岩浆岩特征研究[J]. 物探与化探, 2021, 45(5): 1217-1225.
[2] 耿圣博, 闫红雨, 安战锋, 关海静, 王志博, 金久强, 徐明, 王鑫, 李冰, 郭琦. 中低山区直升机重磁缓起伏测量方法研究[J]. 物探与化探, 2021, 45(5): 1248-1255.
[3] 雍凡, 刘子龙, 蒋正中, 罗水余, 刘建生. 城市三维地震资料处理浅层成像关键技术[J]. 物探与化探, 2021, 45(5): 1266-1274.
[4] 刘海宁, 韩宏伟, 魏文, 张云银, 赵景蒲. 高密度三维区沙四段灰岩有利储层地震预测[J]. 物探与化探, 2021, 45(5): 1281-1287.
[5] 刘彦涛, 彭莉红, 孙栋华, 张伟盟, 王海红. 基于三维有限元的航空大地电磁倾子响应特征[J]. 物探与化探, 2021, 45(5): 1329-1337.
[6] 田郁, 乐彪. 复杂异常体模型下的三维MT倾子正演模拟[J]. 物探与化探, 2021, 45(4): 1021-1029.
[7] 智庆全, 武军杰, 王兴春, 孙怀凤, 杨毅, 张杰, 邓晓红, 陈晓东, 杜利明. 在瞬变电磁三维正演中的应用[J]. 物探与化探, 2021, 45(4): 1037-1042.
[8] 张小艳, 许慧, 姜水军. 基于UPSO-Kriging的综采工作面三维建模研究[J]. 物探与化探, 2021, 45(4): 1071-1076.
[9] 郭楚枫, 张世晖, 刘天佑. 三维磁场有限元—无限元耦合数值模拟[J]. 物探与化探, 2021, 45(3): 726-736.
[10] 李志强, 孙洋, 谭捍东, 张承客. 带地形的ZTEM倾子资料三维正反演研究[J]. 物探与化探, 2021, 45(3): 758-767.
[11] 刘昊娟. 地应力地震预测及其在南川页岩气开发中的应用[J]. 物探与化探, 2021, 45(3): 560-568.
[12] 周月, 官大维, 延海涛, 张小龙. 基于先验信息约束的重磁电联合三维交互反演技术实践——以彭山穹隆构造为例[J]. 物探与化探, 2021, 45(2): 308-315.
[13] 顾观文, 武晔, 石砚斌. 基于矢量有限元的大地电磁快速三维正演研究[J]. 物探与化探, 2020, 44(6): 1387-1398.
[14] 周俊杰, 陈聪, 喻翔, 高玲举, 陈涛. 矿区重力三维物性反演的参考模型构建方法及应用[J]. 物探与化探, 2020, 44(4): 878-885.
[15] 赵文举, 赵荔, 杨战军, 陶德强. 插值切割位场分离方法改进及其在资料处理中的应用[J]. 物探与化探, 2020, 44(4): 886-893.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址:北京市学院路29号航遥中心 邮编:100083
电话:010-62060192;62060193 E-mail:whtbjb@sina.com