地震波逆时偏移方法研究综述_图文

1、第33卷第3期2010年6月文章编号 :1671-8585(2010 03-0153-07勘探地球物理进展P rogr ess in Explor ation Geophysics Vol. 33, No. 3Jun. , 2010地震波逆时偏移方法研究综述陈可洋(中国石油大庆油田有限责任公司勘探开发研究院, 黑龙江大庆163712摘要:地震波逆时偏移是目前精度较高的深度域成像方法, 不同学者从不同的角度对逆时偏移方法进行了研究。因此, 在回顾地震波逆时偏移方法发展历程的基础上, 概述了当前流行的几种叠前、叠后逆时成像方法及其适用性, 通过模型实例分析了几种叠后逆时成像方法存在的一些问题, 同

2、时对比了几种叠前逆时成像方法对Mar m 2ousi 模型的逆时成像效果, 最后指出了逆时偏移方法的发展方向。关键词:逆时偏移; 波动方程; 爆炸反射面; 半速替换; 初至旅行时; 相关法; 逆时成像中图分类号:P63114文献标识码:A地震波逆时偏移法是现行偏移方法中最精确的一种成像方法, 是一种基于波动理论的深度域偏移方法, 它与地震波正演数值模拟问题刚好相反, 将检波器接收到的波场进行逆时延拓, 利用逆时成像条件实现对地下各点的成像。由于逆时偏移方法采用全波场波动方程, 并对波动方程中的微分项进行差分离散(如有限差分法、有限元法、伪谱法等 , 对波动方程的近似较少(这取决于差分方程与微分

3、方程的逼近程度 , 因此不受倾角和偏移孔径的限制, 可以有效地处理纵横向存在剧烈变化的地球介质物性特征(如速度、密度等 。Whitmore 1(1983 最早在第53届SEG 年会上提出逆时偏移方法, Baysal 等(1983 提出了不同的逆时偏移概念, Loewenthal 和Mufti 3(1983 将逆时偏移方法应用在空间-频率域。Levin 等(1984 概括了逆时偏移的基本原理和实现方法。H ildebrand 等5(1987 将逆时偏移方法应用于波阻抗成像, 取得了很好的效果。Levy 和Esmersoy642三维逆时成像。Yoon 等15(2004 指出逆时偏移方法需要解决的问

4、题, 如逆时偏移方法可以消除多次波效应, 但是无法压制逆时偏移过程中的次生干扰。Zhang 等16(2007 提出真振幅逆时偏移的计算策略, Sava 等17(2005 提出了时移成像方法, 薛东川等(2008 采用有限元法实现了弹性波动方程叠前逆时偏移, Guan 等19(2008 采用多步法高效实现了逆时偏移, Liu 等20(2008 提出了反假频波动方程正演和逆时偏移方法, Wards 等2118(2008 提出了相移时间步逆时偏移方法, Soubaras 等22(2008 提出了两步法显式匹配逆时偏移方法, Jones 23(2008 提出了逆时偏移的预处理方法, Guo 等(200

5、8 采用三角网格有限差分法实现了叠后逆时偏移。Chattopadhyay 等25(2008 和Costa 等26(2009 提出了时空移相关法叠前逆时成像方法, 实现了复杂介质叠前逆时成像。陈可洋(2009 提出基于旅行时成像条件的高阶交错网格有限差分叠前逆时偏移方法, Zhang 等31(2009 采用一步内插法实现逆时深度偏移, Dus 2saud 等(2008 研究了逆时偏移的计算策略问题等。本文对逆时偏移方法进行了数值计算和归纳总结。32273024(1988 研究了逆时波场外推、成像和反演问题。Zhu 和Lines 7(1994 比较了逆时偏移与克希霍夫积分偏移方法, 得出前者对Ma

6、r mousi 模型成像精度更高, 但耗时多的结论。Wu 等(1996 研究了三维高阶有限差分法逆时偏移技术, Causse 和U rsin 9(2000 进行了粘弹性波动方程试算, 证明逆时偏移方法对粘弹性波一样适用。Sun 和McMechan (2001 进行了标量波动方程的逆时深度偏移, 对纵波和横波的成像结果表明, 弹性波逆时偏移比单纯声波效果更好。张美根等(2001 实现了各向异性弹性波叠前、叠后逆时深度偏移, 张会星等12(2002 实现了弹性波动方程叠前逆时偏移, Mul 2der 等(2003 比较了单程和双程波动方程偏移, 1413111081 逆时成像方法文中采用的地震波正

7、演递推算子和逆时延拓算子的具体计算公式可参见文献27, 在此不再赘述。下面分述4种逆时成像条件及其适用性。收稿日期:2009-12-28; 改回日期:2010-02-04。作者简介:陈可洋(1983 , 男, 助理工程师, 硕士, 主要从事地震资154勘探地球物理进展第33卷1. 1 基于爆炸反射界面原理和零时间的成像方法爆炸反射界面成像方法是最常用、最简单的一种成像方法, 适合于单程旅行时零偏移距剖面的叠后偏移成像, 不适合于叠前逆时成像。它把地下反射界面假设为爆炸源, 爆炸源的形状、位置与反射界面的形状和位置完全一致, 它所产生的脉冲波的强度和极性与界面反射系数的大小和正负一致, 所以称地

8、下反射界面为爆炸反射界面。在零时刻所有爆炸反射界面同时起爆, 在地表接收地震波场, 接着将地表接收到的波场以给定的模型速度沿逆时间方向延拓至零时刻, 此时的波场值就能够正确描述地下反射界面的位置, 自动实现叠后逆时深度成像。1. 2 半速替换原理和零时间的成像方法该方法同样基于爆炸反射界面的基本假设, 不同的是, 它考虑的是双程旅行时零偏移距剖面偏移成像, 不考虑界面反射系数的大小与极性, 同样不适用于叠前逆时成像。该方法将震源置于整个地表, 在零时刻所有震源同时起爆, 沿着速度模型向地下传播, 遇地质界面后反射至地表, 并被所有地表检波器所接收。接着, 检波器接收到的波场以给定介质速度的一半

9、向地下传播(这里采用半速替换的目的是使双程旅行时等效为单程旅行时, 因为这里只考虑反射界面至地表检波器的距离, 它是地表 地下界面 地表总距离的一半(基于水平介质假设 。速度减半后, 双程旅行时就自动等效于单程旅行时 , 并沿着逆时间方向延拓至零时刻, 此时的波场已经能够正确地反映出地下反射界面的位置, 从而实现了叠后逆时成像。1. 3 初至旅行时成像该方法(可参见陈可洋提出的基于旅行时成像条件的叠前逆时成像方法 不同于以上两种方法, 它通常适用于叠前声波和弹性波逆时偏移成像, 本文将其推广应用于声波叠后双程旅行时零偏移距逆时偏移成像。它采用波场传播的时间一致性成像原理, 即地震波到达地下某一

10、反射界面的时间与地表接收波场逆时延拓至该反射界面位置的旅行时一致, 否则该位置不在反射界面上, 由此筛选出符合条件的成像值, 组合后得到逆时偏移剖面。其具体过程可描述为:将震源置于地表(对叠后逆时偏移情况, 将震源置于整个地表, 而对于叠前逆时偏移情况, 则只置于点震源位置 , 在零时刻所有震源同时起爆, 沿着速度模型向地下传播, 遇地质界面后反射至地表被所有地表检波器接收, 接着27件, 并将检波器接收到的波场以给定介质速度沿逆时间方向向地下延拓, 在延拓的每一个时间步中, 根据已知的初至旅行时成像条件, 筛选出具有相同旅行时的波场值作为其成像值, 从而实现逆时成像。1. 4 相干成像该方法

11、不同于以上3种方法, 它适合于处理具有双程旅行时特征的叠前、叠后声波和弹性波逆时偏移成像。它采用波场传播的最大相干性成像原理, 即地震波到达地下某一反射界面的波场与地表接收到的波场逆时延拓至该反射界面的波场具有最大相干性。如果两者的相干性不是最大, 那么该位置不在反射界面上, 由此记录最大相干值并累加作为该位置处的成像值, 从而实现逆时偏移。其具体过程与初至旅行时成像方法相似, 差别在于不必求取初至旅行时成像条件, 而是对地表接收波场与震源波场在同一时刻作相干成像处理, 筛选出某一空间位置处的最大相干累加值作为其成像值, 从而实现逆时成像。目前典型的几种相干成像方法如下(以二维为例 。1. 4

12、. 1 非时空移相关法这是最为常用的成像方法, 它是将震源波场(S 与接收波场(R 在同一空间位置(x, z 作互相关运算再累加来估算该位置处的成像值I (x, z :I (x, z =E E Sst sXs(x, z , t *R s (x, z, t (1其频率域表达式为:I (x, z =S (x, z, X R (x, z, X E Ess(2式中:x 和z 分别表示水平和深度轴, t 为时间, 下标s 代表炮集个数。该成像方法对小阻抗差结构的成像是有效的, 如图1所示。式(2 与常规单程波偏移成像方法相比, 两者的计算公式相同, 这很好地说明了单程波偏移方法和双程波逆时偏移方法具有相

13、似的理论基础。1. 4. 2 非时空移反射系数相关法该方法借鉴了照明分析的基本思想, 并可进一步压制次生干扰, 其计算公式如下:震源照明I 1(x, z =E S (x, z, t *R (x, z, t EE S (x, z, tssts2s t(3接收照明I 2(x, z =E S (x, z, t *R (x, z, t EE R (x, z, tssts2s t(4第3期陈可洋1地震波逆时偏移方法研究综述155式(3 和式(4 在频域也与常规单程波偏移算法计算公式相同。为了避免分母自相关值接近于零, 导致数值溢出(发散 , 对式(3 和式(4 的分母进行修正处理, 增加一修正因子E (

14、也可称为稳定性因子 , 于是得到修正后的计算公式如下:震源照明I 1(x, z =接收照明I 2(x, z =E S (x, z, t *R (x, z, t EE R (x, z, t +Essts2s t(6式(5 和式(6 即为平滑逆时成像方程。Guit 2ton 等33(2007 指出, 在理论上, 修正因子E 取决于地震资料的信噪比, 而在实际计算中, 可通过选(5取不同的修正因子E 值作交互试验, 分析其剩余谱 特征来实现参数的合理选取。EsE S s (x, z , t *R s (x, z, ttE St2s(x , z, t +E图1 震源波场(a 、接收波场(b 和震源波场

15、与接收波场的相关处理(c1. 4. 3 时空移相关法时空移成像方法是通过震源波场与接收波场在空间方向和时间方向均作互相关运算来估算最佳成像值的一种方法。它由Sava 和Fomel 17(2006 提出, 是1. 4. 1方法的一般情况, 其计算公式如下:I (m =子波作为震源, 声波速度从上到下依次为2000m/s 、2500m/s 、3000m/s 、3500m/s, 时间步长为015ms, 纵横向空间步长均为5m 。采用一阶双曲型声波方程, 应用高阶交错网格有限差分法进行离散, 差分精度为时间2阶和空间10阶, 在边界处采用完全匹配层吸收边界条件, 制作出不同零偏移距剖面, 并实现相应的

16、叠后逆时深度偏移。(7EsE S s (m +h, t +S *tR s (m-h, t -S 在频域可表示为 I (m =Es2i X SS s (m+h, X R s (m-h, X e E X(8式中:i=-1, m 代表空间坐标位置(x, z , h 代表空间方向相关计算位移量(h x , h z , S 代表时间方向相关计算位移量。由该方程同样可以推导出时空移反射系数相关法逆时成像方程。Sava 等17(2006 还指出时移比空移更为有效, 可以大大减小计算量和存储量。2 实例分析2. 1 叠后逆时深度偏移以二维带凹陷的层状介质模型为例, 如图2所 10001000图2 层状介质速度

17、模型图3为基于爆炸反射界面原理和零时间成像156勘探地球物理进展第33卷深度偏移剖面, 采用的主频为30H z, 总旅行时间为1s; 图4和图5均为基于半速替换和零时间成像方法得到的零偏移距剖面和叠后逆时深度偏移剖面, 总旅行时间为2s, 采用的主频分别为40H z 和10H z; 图6是采用主频40H z 、基于自激自收方式合成双程旅行时零偏移距剖面, 再根据旅行时成像方法得到的叠后逆时深度偏移剖面。对比图3a 、图4a 、图5a 可知, 尽管3种方法采用的主频不同, 但剖面中均无明显的数值频散和边界反射现象, 说明合成零偏移距剖面的计算精度较高。对比图3b 、图4b 、图5b 和图6可知,

18、 采用爆炸反射界面原理和零时间成像方法的叠后逆时偏移剖面能够准确地恢复地质界面, 且信噪比较高。采用半速替换和零时间成像方法的叠后逆时偏移剖面也可以准确地恢复水平层位, 但剖面中存在较强的频散现象(即伪波动 , 其随着主频的降低而逐步减弱, 并且绕射波没有收敛。究其原因, 首先是由于检波器接收到的反射波是双程旅行时, 而绕射波传播的路径是从绕射点至地表检波器, 因此是单程旅行时, 采用半速替换后, 绕射波变成了/半单程旅行时0, 无法实现收敛。其次, 根据K =和N =(这里f $hK 为波长, v 为波速, f 为频率, $h 为空间步长, N 为每波长所占有的空间节点数 , 在给定频率f

19、的情况下进行半速替换, 波长K 必然减小, 而此时空间步长$h 大小不变, 那么每个波长所占空间样点数必然减少一半。通常每个波长占8个采样点时才能保证整个计算过程具有较弱的频散现象(只有精细网格剖分才能最大程度压制频散, 而这必然带来巨大的计算量 , 而这里半速替换后, 每个波长所占采样点数减小为4个, 从而导致了频散现象显著。由此可见, 该方法只适合于低频和水平介质的叠后逆时偏移。这也从侧面说明, 该方法适合于CDP 域的叠前逆时成像, 因为 CDP 第3期陈可洋1地震波逆时偏移方法研究综述157图5 主频10H z 时, 基于自激自收方式合成的零偏移距剖面(a 及其半速替换后的叠后逆时深度

20、偏移剖面(b域波场的偏移处理是基于水平层状介质假设的。为了解决反射波和绕射波所具有的单、双程旅行时问题, 本文采用垂直旅行时成像方法, 这样可以保证绕射波准确地收敛于绕射点, 反射波准确归位(图6 。由于不需要进行半速替换, 因此偏移剖面无频散现象。对比图4b 、图5b 和图6也不难发现, 剖面中均存在一些次生干扰, 这是基于双程波 动方程方法进行逆时偏移所难免的。另外还需指出, 基于旅行时成像条件的叠后逆时偏移剖面中, 断层会削弱部分深层同相轴的能量。2. 2 叠前逆时深度偏移采用3种逆时成像方法分别对Marmousi 速度模型(图7 进行叠前逆时成像试验。模型总大小为1915m 605m,

21、 空间纵横向网格间距均为5m, 最小速度为1500m/s, 最大速度为5500m/s, 密度为1. 0g/cm 3, 时间步长为014ms, 满足稳定性条件。采用最大频率为80H z 的零相位雷克子波作为震源, 采用右侧放炮多次覆盖的观测系统, 震源和检波器均置于地表, 第1炮位于模型的420m 处, 炮间距为10m, 道间距为5m, 排列长度为420m, 共激发了151炮。图8a图8c 分别为应用非时空移相关成像方法、非时空移反射系数相关成像方法和基于旅行时成像方法的叠前逆时深度偏移剖面。对比图8与图7可以看出, 各种地质构造形态成像清晰准确, 3条右倾主干大断层、向斜和背斜构造及不整合界面

22、清晰可见, 局部的小断块也清晰可辨, 说明这几种成像方法的逆时偏移噪声能量较小, 计算精 度较高。图6 根据零偏移距剖面(图4a 并基于旅行时成像条件生成的叠后逆时深度剖面158 勘探地球物理进展 第 33 卷 图8 叠前逆时深度偏移剖面 a 非时空移相关成像方法; b 非时空移反射系数相关成像方法; c 基于旅行时成像条件 和存储量均是非常庞大的, 特别是三维逆时偏移情 3 结束语 地震波逆时偏移方法是当前勘探或计算地球 况, 这主要是由相关成像运算造成的; 它的主频 有限, 在保证逆时偏移无数值频散的前提下, 要提 高逆时偏移剖面的主频必然要减小空间网格步长, 这势必增加更多额外的计算量和

23、存储量。 综上分析可知, 逆时偏移的工业化应用必需建 立在快速、 高效的计算机水平的基础上。 参 1 物理学界的研究热点, 该方法采用精确差分离散计 算方法来求解双程波动方程, 不需要进行上下行波 分解, 能有效处理多值走时问题, 并可以通过简单 的直达波切除和窗口相关成像处理, 或通过应用旅 行时成像方法拾取成像值来实现复杂介质高精度 逆时偏移成像。 目前, 逆时偏移仍存在一些有待继续研究的问 题, 如: 它是基于双程地震波动方程, 不可避免地 在逆时偏移过程中会产生次生干扰; 无论采用哪 种成像方法和数值离散计算方法, 它所需的计算量 考 文 献 Whitmore N D. Iterati

24、ve depth imaging by backwar d time pr opagation J . Expanded Abstracts of 53 rd An2 nual Inter national SEG Meeting, 1983, 382 384 2 3 Baysa l E, Kosloff D D, Sherwood J W C. Rever se time migration J . Geophysics, 1983, 48( 11 : 1 514 1 524 Loewenthal D, Mufti I R. Rever se 2time m igr ation in 第3期

25、 陈可洋1 地震波逆时偏移方法研究综述 159 spatial frequency domain J . Geophysics, 1983, 48( 5 : 627 635 4 5 Levin S A. Pr inciples of r everse t ime migr at ion J . Ge2 ophysics, 1984, 49( 5 : 581 583 Hildebr and S T. R everse 2time dept h migration: Imped2 ance imaging condition J . Geophysics, 1987, 52 ( 8 : 21 20

26、 of 78t h Annual Internationa l SEG Meet ing, 2008, 2 341 2 345 Liu F Q, Zhang G Q, Mor ton S A, et a l. An ant i dis2 2 per sion wave equation for modeling and r everse 2time migration J . Expanded Abst racts of 78 t h Annual In2 ternational SEG Meeting, 2008, 2 277 2 281 War ds B D, Mar gr ave G F

27、 . Phase 2shift time 2stepping for rever se ime migration J . Expanded Abstracts of 2t 78 th Annual International SEG Meeting, 2008, 2 262 2 266 22 Soubar as R, Zhang Y. Two 2step explicit marching method for r everse time migration J . Expanded Ab2 st racts of 78 t h Annual Inter nat ional SEG Meet

28、ing, 2008, 2 272 2 276 Jones I F . P re 2processing consider at ions for r ever se time migrat ion J . Expanded Abstracts of 78 t h Annual International SEG Meeting, 2008, 2 297 2 301 24 Guo S J, Li Z C, Sun X D, et a l. Post2stack r everse 2time migration using a finite difference method based on t

29、ri2 angular gr ids J . Applied Geophysics, 2008, 5 ( 2 : 115 120 25 Chattopadhyay S, Mcmechan G A. Imaging conditions for prestack r everse 2time migrat ion J , Geophysics, 2008, 73( 3 : S81 S89 Costa J C, Silva Neto F A, Alcantara M R M, et a l. Ob2 liquity 2correction imaging condition for rever s

30、e time migration J . Geophysics, 2009, 74( 3 : S57 S66 27 28 陈可洋. 基于高阶 有限差 分的波 动方程 叠前逆 时偏移 方法 J . 石油物探, 2009, 48( 5 : 475 478 陈可洋, 杨微, 吴清岭, 等. 几种地震波叠后深度偏移方 法的比较 J . 勘 探地球 物理 进展, 2009, 32( 4 : 257 260 陈可洋. 标量声波 波动方 程高阶 交错网 格有限 差分法 J . 中国海上油气, 2009, 21( 4 : 232 236 陈可洋. 高阶弹性 波动方 程正演 模拟及 其逆时 偏移成 像研究 D

31、. 黑龙江大庆 : 大庆石油学院, 2009 Zhang Y, Zang G. One 2step extr apolation method for reverse time migr ation J . Geophysics, 2009, 74 ( 4 : A29 A33 32 Dussaud E, Symes W W, Williamson P, et al. Computa 2 tional strategies for r everse 2time migration J . Ex2 panded Abstracts of 78 t h Annual International SE

32、G Meeting, 2008, 2 267 2 271 33 Guitton A, Valenciano A, Bevc D, et al. Smoothing im 2 aging condition for shot2pr ole migr ation J . Geo2 physics, 2007, 72( 3 : S149 S154 1 060 1 064 6 Levy B C, Esmer soy C. Var iable background Bor n in2 version by wavefield back pr opagation J . Jour nal on Appli

33、ed Mathematics, 1988, 48( 4 : 30 32 7 8 Zhu J, Lines L. Imaging of complex subsur face str uc 2 tur es by VSP migration J . CJEG, 1994, 30( 1 : 73 83 Wu W, Lines L R, Lu H . Analysis of higher 2order, fi2 nite 2differ ence schemes in 32D reverse 2time migrat ion J . Geophysics, 1996, 61( 3 : 845 856

34、 Causse E, U rsin B. Viscoacoustic rever se ime migra2 2t tion J . Journal of Seismic Exploration, 2000, 9 ( 1 : 23 9 165 184 10 Sun R, M cMechan G A. Scalar rever se 2time depth mi2 gration of elast ic seismic data J . Geophysics, 2001, 66 ( 5 : 1 515 1 518 11 张美根, 王 妙月. 各 向异 性弹 性波 有限 元叠 前 逆时 偏 移 J

35、 . 地球物理学报, 2001, 44( 5 : 711 719 12 张会星 , 宁书年. 弹性波动方程叠前逆时偏移 J . 中国 矿业大学学报, 2002, 31( 5 : 372 375 13 Mulder W, Plessix R. One 2way and two 2way wave 2e2 rd quation migration J . Expanded Abstr act s of 73 An2 nual International SEG Meeting, 2003, 1 292 1 295 Yoon K, Shin C, Suh S, et al. 3D reverse 2

36、time migra2 tion using acoustic wave equation: An experience with the SEG/ EAGE data set J . The Leading Edge, 2003, 15 22( 1 : 38 41 Yoon K, Marfurt K J, Houston U, et al. Challenges in rever se ime migration J . Expanded Abstr act s of 74t h 2t Annual Internationa l SEG Meeting, 2004, 1 057 1 060

37、16 Zhang Y, Sun J, Gray S. Reverse 2time migration: Am2 plitude and implementation issues J . Expanded Ab2 str act s of 77 th Annual Inter national SEG Meeting, 17 2007, 2 145 2 149 Sava P, Fomel S. T ime 2shift imaging condition J . Ex2 panded Abstracts of 75 t h Annual Inter nat ional SEG Meet ing

38、, 2005, 1 850 1 853 18 薛东川 , 王尚旭. 波动方程有限元叠前逆时偏移 J . 石 油地球物理勘探, 2008, 43( 1 : 17 21 19 Guan H, Li Z, Wang B, et al. A multi step approach for 2 efficient reverse 2time migration J . Expanded Abstracts 26 14 29 30 31 ( 编辑: 戴春秋 第 33 卷第 3 期 2010 年 6 月 勘探地球物理进展 P rogr ess in Explor ation Geophysics Vol.

39、33, No. 3 Jun. , 2010 ABSTRACT Review of seismic reverse time migr ation methods. Chen Keyang. PEG, 2010, 33( 3 : 153 159 Seismic rever se time migration is a method of high imaging precision in depth domain. Different realizations of reverse mi2 gration have been put forward. This paper review ed t

40、he history of seismic rever se time migration, and summarized several prev2 alent prestack and poststack reverse time migration schemes and their adaptability. We demonstrated different poststack re2 ver se time imaging conditions with model data, and compar ed the results fr om different reverse ti

41、me migration methods on Marmousi model. The outlook of r everse time migr ation was al2 so forecast. Keywords: reverse time migration; wave equation; exploding re2 flector inter face; half 2speed substitution; first arrival; correla2 tion method; reverse time imaging condition Chen K eyang, Explorat

42、ion and Development Research Institute of Daqing O ilfield Company Limited, Daqing 163712, China Review of quality control of seismic data processing in the course of integrated study of seismology & geology. Li Dawei, Yang Kai, Yang Shunjun. PEG, 2010, 33( 3 : 160 167 W ith more and more seismology

43、 & geology integrated work being implemented in geophysical exploration industry, the quality control of seismic data processing has been received great attention. In this paper, three typical QC methods used in industry: well2control Q C method based on VSP data, inver sion2 oriented QC which is ba

44、sed on rock physics modeling, and con2 sistency iented Q C based on a delicate processing to keep true 2or amplitude, fr equency, phase and waveform relatively, w er e re2 view ed in detail. Discussions on the future development of QC methods of seismic data processing were suggested. Keywords: inte

45、gr ated study of seismology & geology; quality control; w ell control seismic data processing Li Dawei, School of Ocean and Ear th Science, Tongji U niver si2 ty, Shanghai 200092, China Realization of time 2domain residual curvature migr ation velocity analysis in iCluster software platform. Zhang B

46、ing, Wang Hua2 zhong, Sun Chenglong, Liu Shaoyong. PEG, 2010, 33 ( 3 : 168 173 Residual cur vature analysis ( RCA based on common ima2 ging point gathers is one of the effective methods in impr oving the accuracy of migration velocity. By analyzing the relationship between the residual curvature of

47、Kirchhoff PSTM common imaging point gathers and migration velocity, we developed a PSTM velocity analysis module on iCluster seismic data pro2 cessing platform. We tested the module with model data and re2 al data. A ll the tests show that this module can adjust the RMS velocity field locally and qu

48、antificationally, and can wor k well with Kirchhoff prestack time migration. Keywords: r esidual curvatur e analysis; migr ation velocity analy2 sis; PSTM; iCluster Zhang Bing, School of Ocean and Ear th Science, Tongji Univer2 sity, Shanghai 200092, China Forward modeling based on anisotropic gas 2

49、bear ing sandstone model and characteristics analysis of AVO response. Zhou Hua2 ilai, Li Luming, Luo Shengxian, Wang Mingchun. PEG, 2010, 33 ( 3 : 174 178 In order to investigate the multi wave char acteristics of e2 2 lastic waves in anisotropic gas 2bearing sandstone r eser voirs to facilitate th

50、e inter pretation of multi wave seismic data, this pa 2 2 per adopts high2order staggered id finite difference method to 2gr simulate propagation of seismic waves in anisotropic gas 2bear 2 ing sandstone media based on anisotropic elastic wave equation. The effect of anisotropic coefficients on mult

51、i wavefield was an 2 2 alyzed, and the responses of the wavefield in isotropic and ani 2 sotropic gas2bear ing sandstone reservoir s were compared. Re 2 sults show that anisotropy significantly affects wavefield char 2 acteristics and AVO response, and the impact manifests different behaviors on P 2

52、waves and converted waves. The amplitude charac 2 teristics in seismic record is consistent with the AVO reponse, indi 2 cating that the combination of forward modeling and AV O anal 2 ysis can provide a solid base for identifying multi wavefield. 2 Keywords: multi wave and multi component; wavefiel

53、d char ac 2 2 2 teristics; AVO analysis; per fectly matched layer; forward modeling Zhou Huailai, State Key Laboratory of Oil and Gas Reser voir Geology and Exploitation and Key Lab of Ear th Explor ation & Information Techniques of Ministry of Education, Chengdu U2 niversity of Technology, Chengdu

54、610059, China Exper iment on ground resistance in high 2density electrical meas 2 ur ement. Zhang Lingyun, Liu Hongfu, Li Chengyou. PEG, 2010, 33( 3 : 179 183 A n experimental study of electrical prospecting with uni 2 form and non2uniform gr ound resistance on the same testing field w as introduced. The variability of appar ent resistance and the effect of ground resistance on sensiti