Walkaway井地联合地震资料Q值波形反演方法.docVIP

Wal ka wa y 井地联合地震资料 Q 值波形反演方法杨宇山 , 李媛媛 , 刘天佑 Wal ka wa y 井地联合地震资料 Q 值波形反演方法 杨宇山 , 李媛媛 , 刘天佑 (中国地质大学 (武汉) 地球物理与空间信息学院 , 湖北 武汉 430074) [ 摘要 ] 传统的品质因子 Q 值估算方法通常是基于 V SP 资料或者叠加的地面地震资料 , 在薄互层存在的情 况下 , 它们难以提取地震子波及有效地分离固有衰减和散射衰减 。为了获得更可靠的地震波衰减图像 , 研究 采用了一种基于 Wal kaway 井地联合地震资料的波形反演来求取 Q 值的方法 。通过理论模型的数值计算分 析 , 探讨了初始模型厚度 、Q 值及噪声对反演结果的影响 。合成的 Wal kaway 井地联合地震资料的衰减估计 结果表明 , 波形反演的 Q 值估算方法充分利用了各种波的运动学和动力学信息 , 避免了初至拾取和波场分 离等处理困难 , 并能减小可能的人为误差 , 提高计算的精度 。将该方法用于某油田的高分辨 Wal kaway 井地 联合地震资料 Q 值估计 , 结果表明 , 介质品质因子 Q 与介质性质和孔隙流体有良好的对应关系 。 [ 关键词 ] Wal kaway 井地联合地震勘探 ; 品质因子 ; 波形反演 [ 中图分类号 ] P6311 44 [ 文献标识码 ] A [ 文章编号 ] 100029752 (2009) 0220053206 当地震波在粘弹性介质中传播时 , 介质粘滞性引起的吸收作用是导致地震波能量衰减的主要因素之 一 , 这种介质所固有的衰减特性通常用品质因子 Q 值来描述 。地下介质的 Q 值是研究地震波能量衰减与 地下流体性质之间关系的重要参数[ 1 ] 。然而在实际中 , 品质因子 Q 值并没有得到充分的利用 , 其主要原因 是缺少一种精确可靠的 Q 值估算方法 。目前常用的 Q 值估算方法有谱比法 、升时法 、质心频移法 、互相关法 等[ 2～7 ] , 然而 , 这些方法在实际应用过程中常常会遇到以下两方面的困难[ 8 ] : ①一般所说的且感兴趣的地 震衰减是指由孔隙流体和岩石结构产生的固有 ( 内在) 的衰减 , 它必须与散射衰减区分开来 。然而在薄互 层情况下 , 所观测到的地震衰减通常是固有衰减和散射衰减的总和 , 由于这两种衰减对地震波的作用非常 相似 , 很难将它们完全分离开来 。②传统的 Q 值估算方法都要对子波进行提取 , 而对于实际资料而言 , 当 有薄层的情况下 , 子波的提取是非常困难的 , 并且由于传统谱分析的局限性 , 很难通过它来提取地震记录 中子波的频谱 。综上所述 , 研究新的 Q 值估算方法是十分必要的 。 笔者利用 Wal kaway 井地联合地震资料 , 采用波形反演的方法对介质 Q 值进行估算 。由于波形反演能 够最大限度地利用各种波的运动学和动力学信息 , 减轻初至拾取和波场分离等处理困难 , 减小可能的人为 误差 , 提高计算的精度 , 因此 , 在理论上该方法优于传统的 Q 值估算方法 。Wal kaway 井地联合地震勘探是 一种在同样的激发条件下 ,同时同步采集地面地震资料和井中 V SP 资料的新技术 , 通过井 、地资料的联 合 , 可以充分利用两种资料的优点 , 相互约束 , 相互补充 。 基本原理 在一维粘弹性介质中 , 上 、下行波之间有如下递推关系式[ 9 ] : 1 1/ w k t k c k w k / tk U k +1 D k +1 U k ( 1) = ck / w k t k w k / tk D k βh k - αh k 式中 , U k 、D k 、ck 、tk 分别表示在第 k 层界面上的上行波 、下行波 、反射系数和透射系数 ; w k = ei e / hk ( 其中β = ω ,α = ω ,ω为角频率 , v k 为相速度 , hk 为层厚度) [ 10 ] 。如果已知介质的层数 N 以及各层反 v k 2 Qv k [ 收稿日期 ] 2009202203 [ 作者简介 ] 杨宇山 ( 19772) , 男 , 2000 年大学毕业 , 博士 , 讲师 , 现主要从事地球物理反演理论及应用研究工作 。 射系数和厚度 ( 如图 1 所示) , 则有 :01/ w k t k c k w k / tkF11F21F12- R1 + 射系数和厚度 ( 如图 1 所示) , 则有 : 0 1/ w k t k c k w k / tk F11 F21 F12 - R 1 + R - R ∏ = = ( 2) E N