取振幅包络体属性，生成目的层段的振幅包络数据体R ;图2为时窗t即1500-2100毫秒之内的原始地震剖面，剖面中椭圆标识的位置为储层发育段，砂岩和火山岩不能有效区分；图3为对图2所示的原始地震数据提取的振幅包络体属性剖面R。
[0032]3)在步骤I)确定的时间窗口 t之内，对原始地震数据体提取瞬时相位余弦体属性，生成目的层段的瞬时相位余弦数据体I ;图4为在1500-2100毫秒时窗内对图2所示的原始地震数据提取的瞬时相位余弦体I的属性剖面。
[0033]4)应用傅里叶变换提取步骤2)生成的目的层段的振幅包络数据体R的振幅谱，振幅谱横坐标为频率，纵坐标为振幅，振幅谱呈随频率增加振幅下降的曲线，曲线由大斜率部分和小斜率部分组成，大斜率和小斜率部分的拐点为突变点，以突变点为界，右侧的小斜率部分为高频段，左侧的大斜率部分为低频段；突变点对应的横坐标为高频段和低频段的门槛频率f ;图5为本例中对图3所示的振幅包络数据体R提取的振幅谱，虚线左边为大斜率部分，虚线右边为小斜率部分，虚线对应的点为突变点，虚线对应的频率为门槛频率，本例中门槛频率为2HZ。
[0034]5)在步骤I)确定的时间窗口 t内，对步骤2)生成的目的层段的振幅包络数据体R，以步骤4)确定门槛频率f为界进行高通滤波，生成目的层段振幅包络数据体R的高通滤波数据体H ;所述的高通滤波是将高于f频率段的振幅保留，低于f频率段的振幅数据滤除。图6是对图3所示振幅包络数据体R根据图5确定的2HZ的门槛频率进行高通滤波的结果，图中大于2HZ的成分保留，小于2HZ的成分被滤掉。
[0035]6)对步骤5)生成的目的层段高通滤波数据体H进行负数赋值，将所有负数赋值为0，所有的正值保持不变，得到高通滤波正值数据体H’ ；图7是对图6所示的高通滤波数据体H进行负数赋值生成的高通滤波正值数据体H’，原有的负数全部被清零。
[0036]7)将步骤6)的到的高通滤波正值数据体H’与步骤3)生成的目的层段瞬时相位余弦数据体I进行乘法运算，生成能量差异属性数据体N ;所述的乘法运算是H’乘以I。图8是将图7所示的高通滤波正值数据体与图4所示的瞬时相位余弦数据体相乘的结果剖面，即能量差异属性数据体N的剖面，剖面中椭圆标识的火山岩和砂岩能够很好区分。
[0037]8)在原始地震数据体上根据步骤I)的井震时间深度匹配，完成勘探区域内目的层段储层或其他地质体的顶、底的层位追踪，生成储层顶界层位Cl和储层底界层位C2 ;所述的其他地质体是火山岩、礁。图9上所示的Cl、C2是根据图1所示的标定结果在图2所示的原始地震数据体上追踪的储层顶界和底界的层位。
[0038]9)将步骤7)生成的能量差异属性数据体N以步骤8)中追踪的层位Cl、C2为间隔时窗提取均方根振幅，生成预测储层或其他地质体平面分布图；图10为以C1、C2为间隔在图8所示的能量差异属性体上提取的最终储层分布图，本例中的储层为火山岩，图中强振幅为火山岩发育段，火山岩和砂岩能够很好区分。
【主权项】
1.一种利用地震能量差异属性预测储层的方法，特点是通过以下步骤实现: 1)激发合成地震并记录，通过勘探区域已知井的钻井分层、声波时差数据、密度数据、地震数据体进行井震时间深度匹配，用地震剖面中目的层段对应的时间范围作为地震体属性提取的时间窗口 t; 2)在时间窗口t之内，对原始地震数据体提取振幅包络体属性，生成目的层段的振幅包络数据体R ； 3)在时间窗口t之内，对原始地震数据体提取瞬时相位余弦体属性，生成目的层段的瞬时相位余弦数据体I ; 4)应用傅里叶变换提取目的层段的振幅包络数据体R的振幅谱，振幅谱横坐标为频率,纵坐标为振幅,振幅谱呈随频率增加振幅下降的曲线，曲线由大斜率部分和小斜率部分组成，大斜率和小斜率部分的拐点为突变点，以突变点为界，右侧的小斜率部分为高频段，左侧的大斜率部分为低频段；突变点对应的横坐标为高频段和低频段的门槛频率f ； 5)在时间窗口t内，对目的层段的振幅包络数据体R，以门槛频率f为界进行高通滤波，生成目的层段振幅包络数据体R的高通滤波数据体H ； 6)对目的层段高通滤波数据体H进行负数赋值，将所有负数赋值为O，所有的正值保持不变，得到高通滤波正值数据体H’ ； 7)将高通滤波正值数据体H’与目的层段瞬时相位余弦数据体I进行乘法运算，生成能量差异属性数据体N; 8)在原始地震数据体上根据步骤I)的井震时间深度匹配，完成勘探区域内目的层段储层或其他地质体的顶、底的层位追踪，生成； 9)将能量差异属性数据体N以储层顶界层位Cl和储层底界层位C2为间隔时窗提取均方根振幅，得到预测储层分布。2.根据权利要求1所述的方法，特点是步骤5)所述的高通滤波是将高于f频率段的振幅保留，低于f频率段的振幅数据滤除。3.根据权利要求1所述的方法，特点是步骤7)所述的乘法运算是高通滤波正值数据体H’乘以目的层段瞬时相位余弦数据体I。4.根据权利要求1所述的方法，特点是步骤8)所述的其他地质体是火山岩或火山岩礁。
【专利摘要】一种利用地震能量差异属性预测储层的方法，在时间窗口内，提取振幅包络体、瞬时相位余弦体属性，生成目的层段的振幅包络数据体R、余弦数据体I，包络数据体R的振幅谱横坐标为频率，纵坐标为振幅，呈随频率增加振幅下降的曲线，其中突变点对应的横坐标为高频段和低频段的门槛频率f，高通滤波生成包络数据体R的高通滤波数据体H，负数赋值，正值保持不变，得到高通滤波正值数据体H’，将其与余弦数据体I乘法运算，生成能量差异属性数据体N，根据井震时间深度匹配，完成勘探区域内目的层段储层或其他地质体的顶、底的层位追踪，将数据体N以储层顶界、底界层位为间隔时窗提取均方根振幅，得到预测储层分布。本发明可以剔除较弱地震振幅，保留地震强振幅。
【IPC分类】G01V1/30
【公开号】CN104977612
【申请号】CN201410142534
【发明人】陈志刚, 于京波, 马辉, 潘良云, 王恺
【申请人】中国石油集团东方地球物理勘探有限责任公司
【公开日】2015年10月14日
【申请日】2014年4月9日