CO2驱油地震频散属性提取方法技术

本发明专利技术提供一种CO2驱油地震频散属性提取方法，该CO2驱油地震频散属性提取方法包括：步骤1，通过广义S变换谱分解提取注CO2前后两期资料的若干单频剖面；步骤2，注气前后地震道集时频分析，综合测井曲线和井旁道提取地震子波；步骤3，两期地震资料谱均衡处理，去除子波叠影的影响，恢复地层反射系数频散效应；步骤4，根据构建的频散属性提取方程，基于阻尼最小二乘提取纵波速度频散属性剖面。该CO2驱油地震频散属性提取方法能够从速度频散属性的角度检测注CO2气后的动态影响范围，利用速度频散的角度去解释CO2驱油的地震反射特征，进而落实注CO2气后的目标油藏影响范围。

【技术实现步骤摘要】

本专利技术涉及地震资料解释领域，特别是涉及到一种CO2驱油地震频散属性提取方法。
技术介绍
根据岩石物理分析,CO2注入地层驱油的过程中会发生速度频散，而速度频散是由于地震波传播过程中相速度随频率发生变化的现象，所以考虑含CO2气情况下速度频散的现象，进行速度频散属性的提取就很有必要。为此我们专利技术了一种新的CO2驱油地震频散属性提取方法，解决了以上技术问题。
技术实现思路
本专利技术的目的是提供一种利用速度频散的角度去解释CO2驱油的地震反射特征，进而落实注CO2气后的目标油藏影响范围的CO2驱油地震频散属性提取方法。本专利技术的目的可通过如下技术措施来实现：CO2驱油地震频散属性提取方法，该CO2驱油地震频散属性提取方法包括：步骤1，通过广义S变换谱分解提取注CO2前后两期资料的若干单频剖面；步骤2，注气前后地震道集时频分析，综合测井曲线和井旁道提取地震子波；步骤3，两期地震资料谱均衡处理，去除子波叠影的影响，恢复地层反射系数频散效应；步骤4，根据构建的频散属性提取方程，基于阻尼最小二乘提取纵波速度频散属性剖面。本专利技术的目的还可通过如下技术措施来实现：在步骤1中，首先对CO2驱油两期资料进行频谱分析确定其频率变化范围，进行广义S变换谱分解，分解为若干个若干单频剖面，为之后的谱均衡处理打下基础。在步骤2中，在注气后高频段信息相比低频段反射衰减明显，频带宽度变窄。在步骤3中，频变纵波反射系数公式为： r ( t , f ) = v p 2 ρ 2 - v p 1 ρ 1 v p 2 ρ 2 + v p 1 ρ 1 = 1 2 [ Δv p v p ( t , f ) + Δ ρ ρ ( t ) ] - - - ( 1 ) ]]>其中，r(t,f)为频变反射系数，t,f分别表示时间和频率，vp为纵波速度，vp2为上层纵波速度，vp1为下层纵波速度，ρ2为上层密度，ρ1为下层密度，Δvp为上下层速度差，Δρ为上下层密度差，vp＝(vp1+vp2)/2，ρ＝(ρ1+ρ2)/2，Δvp＝vp1-vp2，Δρ＝ρ1-ρ2；根据平稳地震记录褶积模型，即自激自收的地震记录等价于地震子波和反射系数的褶积，其在频率域相当于反射系数谱和子波谱的乘积，即：r(t,f)*w(t,f)＝s(t,f) (2)其中,r(t,f)为反射系数，w(t,f)为子波，s(t,f)为地震记录；将地震记录频谱展宽，即去除“子波叠影”的影响，进而恢复仅由复杂地下介质引起的声波反射信息，假设地层反射系数符合高斯白噪分布，则谱均衡具体过程为 r ( t , f ) = s ( t , f ) w ( t , f ) + ϵ - - - ( 3 ) ]]>其中，ε为正无穷小量，主要用来提高谱均衡过程的稳定性；通过子波的频域响应确定谱均衡系数，再分别对广义S变换谱分解后的剖面进行谱均衡处理。在步骤4中，将式(3)中r(f)在f0处泰勒展开，忽略掉二阶以及二阶以上高阶项，得到： r ( t , f ) ≈ r ( t , f 0 ) + ∂ [ r ( t , f ) ] ∂ f d f - - - ( 4 ) ]]>由于密度不会随着频率的变化而发生变化，因此密度对频率的偏导数忽略不计，将(1)式代入： r ( t , f ) - r ( t , f 0 ) = 1 2 ∂ Δv p v p 本文档来自技高网...

【技术保护点】
CO2驱油地震频散属性提取方法，其特征在于，该CO2驱油地震频散属性提取方法包括：步骤1，通过广义S变换谱分解提取注CO2前后两期资料的若干单频剖面；步骤2，注气前后地震道集时频分析，综合测井曲线和井旁道提取地震子波；步骤3，两期地震资料谱均衡处理，去除子波叠影的影响，恢复地层反射系数频散效应；步骤4，根据构建的频散属性提取方程，基于阻尼最小二乘提取纵波速度频散属性剖面。

【技术特征摘要】
1.CO2驱油地震频散属性提取方法，其特征在于，该CO2驱油地震频散属性提取方法包括：步骤1，通过广义S变换谱分解提取注CO2前后两期资料的若干单频剖面；步骤2，注气前后地震道集时频分析，综合测井曲线和井旁道提取地震子波；步骤3，两期地震资料谱均衡处理，去除子波叠影的影响，恢复地层反射系数频散效应；步骤4，根据构建的频散属性提取方程，基于阻尼最小二乘提取纵波速度频散属性剖面。2.根据权利要求1所述的CO2驱油地震频散属性提取方法，其特征在于，在步骤1中，首先对CO2驱油两期资料进行频谱分析确定其频率变化范围，进行广义S变换谱分解，分解为若干个若干单频剖面，为之后的谱均衡处理打下基础。3.根据权利要求1所述的CO2驱油地震频散属性提取方法，其特征在于，在步骤2中，在注气后高频段信息相比低频段反射衰减明显，频带宽度变窄。4.根据权利要求1所述的CO2驱油地震频散属性提取方法，其特征在于，在步骤3中，频变纵波反射系数公式为： r ( t , f ) = v p 2 ρ 2 - v p 1 ρ 1 v p 2 ρ 2 + v p 1 ρ 1 = 1 2 [ Δv p v p ( t , f ) + Δ ρ ρ ( t ) ] - - - ( 1 ) ]]>其中，r(t,f)为频变反射系数，t,f分别表示时间和频率，vp为纵波速度，vp2为上层纵波速度，vp1为下层纵波速度，ρ2为上层密度，ρ1为下层密度，Δvp为上下层速度差，Δρ为上下层密度差，vp＝(vp1+vp2)/2，ρ＝(ρ1+ρ2)/2，Δvp＝vp1-vp2，Δρ＝ρ1-ρ2；根据平稳地震记录褶积模型，即自激自收的地震记录等价于地震子波和反射系数的褶积，其在频率域相当于反射系数谱和子波谱的乘积，即：r(t,f)*w(t,f)＝s(t,f) (2)其中,r(t,f)为反射系数，w(t,f)为子波，s(t,f)为地震记录；将地震记录频谱展宽，即去除“子波叠影”的影响，进而恢复仅由复杂地下介质引起的声波反射信息，假设地层反射系数符合高斯白噪分布，则谱均衡具体过程为 r ( t , f ) = s ( t , f ) w ( t , f ) + ϵ - - - ( 3 ) ]]>其中，ε为正无穷小量，主要用来提高谱均衡过程的稳定性；通过子波的频域响应确定谱均衡系数，再分别对广义S变换谱分解后的剖面进行谱均衡处理。5.根据权利要求4所述的CO2驱油地震频散属性提取方法，其特征在于，在步骤4中，将式(3)中r(f)在f0处泰勒展开，忽略掉二阶以及二阶以上高阶项，得到： r ( t , f ) ≈ r ( t , f 0 ) + ∂ [ r ( t , f ) ] ∂ f d f - - - ( 4 ) ]]>由于密度不会随着频率的变化而发生变化，因此密度对频率的偏导数忽略不计，将(1)式代入： r ( t , f ) - r ( t , f 0 ) = 1 2 ∂ Δv p ...

【专利技术属性】
技术研发人员：谭明友，张军华，张云银，崔世凌，曲志鹏，于景强，李军，亓亮，雷蕾，李晓晨，张伟忠，吴明荣，肖文，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司物探研究院，
类型：发明
国别省市：山东;37