本发明属于油气地球物理
技术领域：
，尤其涉及一种基于叠前地震数据的裂缝储层快速预测方法。
背景技术：
：根据采集的地震数据，现有的基于叠前地震数据裂缝预测技术，可分为针对宽方位角叠前地震数据的裂缝预测方法和针对窄方位角叠前地震数据的裂缝预测方法两类：(1)针对宽方位角叠前地震裂缝预测，国内外学者大多采用的是基于P波各向异性地震属性椭圆拟合技术，如Grechka和Tsvankin等(1999年)提出对水平层状介质中关于横向速度变化的动校正时差椭圆拟合，来预测随深度变化的裂缝倾向；曲寿利等(2001年)提出了根据阻抗随方位角变化来检测裂缝的方法；Gray(2004年)提出AVO梯度随全方位角变化的方法来预测裂缝(其过程与拟合椭圆相近)；孙炜等(2014年)提出了一种改进的方位各向异性裂缝预测方法。以上方法实现过程均较复杂，且在处理窄方位角地震数据的裂缝预测时，拟合椭圆受到限制；(2)针对窄方位角地震资料的裂缝预测，国内学者也提出了一些技术方法，如田立新等(2010年)通过远近偏移距属性差异地震正演和方位地震属性差异的地震正演，模拟得到针对窄方位地震采集区块的叠前裂缝储层预测方法；苏世龙等(2015年)通过采取限炮检距接收及一些特殊的数据规则化技术改善窄方位资料的性质从而进行叠前裂缝预测。这些针对窄方位角地震数据裂缝预测的技术方法，采取分方位角处理，过程复杂，理论不够完善。此外，多波技术也用于裂缝预测中，多波资料相对于裂缝的敏感程度更强，但是采集多波分量数据的施工成本太高，资料处理难度大，在实际生产中多波采集较少，因此多波技术受到很大的限制，不利于普遍适用。综上所述，现有的基于叠前地震数据裂缝预测方法存在实现过程复杂，拟合椭圆受到限制，施工成本太高，资料处理难度大。技术实现要素：本发明的目的在于提供一种基于叠前地震数据的裂缝储层快速预测方法，旨在解决现有的基于叠前地震数据裂缝预测方法存在实现过程复杂，拟合椭圆受到限制，施工成本太高，资料处理难度大的问题。本发明是这样实现的，一种基于叠前地震数据的裂缝储层快速预测方法，所述基于叠前地震数据的裂缝储层快速预测方法根据叠前地震资料的品质，采用去噪-拉平-切除的组合优化对地震资料进行处理，在优化处理后的道集上，拟合某一时间点各道的振幅值，将该时间点各道的拟合值减去看作均匀介质的各道近道集振幅值的平均值，得到各道的差值；求其标准差作为异常值，可用该异常值有效地表征裂缝储层发育强度；异常值看作是裂缝和流体的综合影响的结果；以采样率为步长，移动时间点，得到整个地震数据的异常值数据体。进一步，所述采用去噪-拉平-切除的组合优化对地震资料进行处理具体包括：步骤一，采用四阶多项式拟合进行去噪处理；步骤二，采用非地表一致性剩余静校正实现道集拉平处理；步骤三，叠前远道集在经过去噪、拉平处理后，若仍然存在杂乱，道集数据缺损、噪音干扰严重等现象，需要对数据采取切除30°以外的道集。进一步，所述采用四阶多项式拟合进行去噪处理包括：(1)基于AVO曲线拟合，采用四阶多项式拟合方法拟合叠前道集某一时刻的数据点振幅值，得到拟合曲线；(2)选择初始参数，将拟合曲线分别向上向下移动，对比同一偏移距的数据，去除高于向上移动值和低于向下移动值的采样点值，用相同偏移距的多项式拟合值代替；(3)检验拟合的AVO曲线，是否与正演道集AVO曲线趋势一致，一致则输出去噪道集，否则调整参数值大小，重复(1)和(2)，输出去噪后道集。进一步，所述采用非地表一致性剩余静校正实现道集拉平处理包括：(1)将一定偏移距内的道集进行叠加作为初始道S(t)；在一定时窗t内，分别计算各道与初始道的相关系数R(t,I)，采用如下公式求取，其中t1、t2为时窗的起止时间，X(t-τ,I)为输入道，I为道号，τ＝0，1，2……n,n为互相关函数的长度；(2)选择相关系数最大的一道作为模型道，从某道开始依次在指定可移动时窗T范围内，以某个时窗t滑动求取与相邻道的相关系数,所用公式与(1)中所述相同；(3)将相关系数最大时所滑动的时窗作为当前道的剩余时差，在指定时间段T内，各道减去看作剩余时差的时间之后，输出拉平道集。进一步，所述统计计算AVO异常值的方法包括：第一步，对优化处理之后的叠前道集一条测线同一时间点上的振幅值进行多项式拟合，采用最小二乘法原理来实现；第二步，根据均匀介质与各向异性介质在入射角0°-6°时其值基本相同，首先根据第一步得到某一时刻各道(这里各道指的是叠前道集包括CRP、CDP、CMP)的拟合值，同时根据深度和偏移距的关系计算近道集角度，并求取各道近道集(入射角0°-6°)振幅值的平均值，将道集各道的拟合值减去看作均匀介质的各道近道集振幅值的平均值，得到各道的差值；第三步，利用如下公式：其中bi为原始值，bii为拟合值，M为叠前道集某一道的总道数，i为道号，S为某道的标准差。统计各道的差值，计算标准差S作为各道此时间点的异常值，异常值可看作是裂缝和流体综合影响的结果，可用该异常值有效地表征裂缝储层发育强度；第四步，以采样率为步长，移动时间点，循环第一步-第三步，得到整条线的地震数据的异常值。第五步，改变测线，循环第一步-第四步，得到整个地震数据的异常数据体，通过剖面和平面成图，得到裂缝储层预测的图形效果，实现裂缝储层预测。进一步，多项式拟合，最小二乘法原理包括：1)在某时刻，若观测到的数据为一个入射角xi对应一个振幅值yi如下表所示：xix1x2……xnyiy1y2…….yn2)假设拟合函数不是线性函数，用如下的公式来进行拟合：y＝a0+a1x+a2x2+a3x3+.......+anxn,n＝0,1,2,3........；3)设拟合值和真实值之间的残差平方和为δi：4)当残差平方最小时，拟合曲线与真实值之间最为接近。所以，其关键问题就是求解多项式的系数。转化问题，关于多项式系数的函数，为得到最小值则求取其导数，令经过一系列的推导可以得到如下公式：求出ai，然后就可以进行拟合曲线。本发明提供的基于叠前地震数据的裂缝储层快速预测方法，适用范围较广、不受宽窄方位角地震资料的限制、实现过程快捷；基于各向异性理论，考虑在均匀介质(不含裂缝和流体的介质)中，根据Zoeppritz方程及其简化式可知随入射角的变化，振幅大小基本上不变；而在各向异性介质(含裂缝和流体)中，振幅随入射角的变化而变化；并验证分析两者的AVO曲线随入射角变化的特征，得知裂缝方向无论如何展布，AVO都会表现一定的差异，常规的裂缝预测方法基本忽视了这一点；本发明将这种差异定性的作为裂缝储层预测的基础，综合考虑了裂缝及流体的综合影响；快速有效地解决宽窄方位角地震数据裂缝储层预测的问题；是在优化处理后的道集上，拟合某一时间点各道的振幅值，将该时间点各道的拟合值减去看作均匀介质的各道近道集(约0°-6°之间)振幅值的平均值，得到各道的差值；求其标准差作为异常值，可用该异常值有效地表征裂缝储层发育强度。如先计算1000ms处的异常值之后，若采样率为2ms,则以采样率为步长，移动时间点到1002ms，计算得到异常值，依次循环到整个记录时长，得到整个时长的异常值。本发明对优化处理之后的叠前道集同一时间点上的振幅值进行多项式拟合；将道集各道的拟合值减去看作均匀介质的各道近道集(约0°-6°之间)振幅值的平均值，得到各道的差值；统计各道的差值计算标准差作为各道此时间点的异常值，此方法相对于常规的基于P波各向异性地震属性椭圆拟合技术的计算量少，计算效率高；且不依赖于各向异性椭圆拟合的裂缝储层预测分析技术，解决宽窄方位角地震数据裂缝储层预测的问题。假设优化处理步骤一样，基于椭圆拟合的方法，假设也用异常值来拟合椭圆，需要全方位角度的计算，还要进行拟合椭圆，来预测裂缝，而本发明的方法，求得异常值，即可作为裂缝预测的依据。本发明的方法对某工区进行实际资料(CRP道集)分析，在井的周围都有裂缝的发育，与实际钻井吻合(应用本发明裂缝储层预测的符合率达到83％)；有效地验证了本发明的实用性。附图说明图1是本发明实施例提供的基于叠前地震数据的裂缝储层快速预测方法流程图。具体实施方式为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。下面结合附图对本发明的应用原理作详细的描述。如图1所示，本发明实施例提供的基于叠前地震数据的裂缝储层快速预测方法包括以下步骤：S101：根据叠前地震资料的品质，首先采用“去噪-拉平-切除”的组合优化技术对地震资料进行处理；S102：对优化处理之后的叠前道集同一时间点上的振幅值进行二次以上的多项式拟合；S103：根据均匀介质与各向异性介质在入射角非常小(约0°-6°之间)时其值基本相同，将道集各道的拟合值减去看作均匀介质的各道近道集(入射角约0°-6°之间)振幅值的平均值，得到各道的差值；S104：统计各道拟合值与原始值的标准差作为各道此时间点的异常值，可用该异常值有效地表征裂缝储层发育强度；S105：以采样率为步长，移动时间点，循环步骤得到整个地震数据的异常值数据体。下面结合具体实施例对本发明的应用原理作进一步的描述。本发明实施例提供的基于叠前地震数据的裂缝储层快速预测方法，具体实施步骤如下：(1)叠前道集优化处理：根据叠前地震资料的品质，首先采用“去噪-拉平-切除”的一个组合优化技术对地震资料进行处理，其具体步骤如下：步骤1：采用四阶多项式拟合进行去噪处理：①基于AVO曲线拟合，通过四阶多项式拟合叠前道集某一时刻的数据点，得到拟合曲线；②选择初始参数，将拟合曲线上下移动，对比同一时刻的数据采样点值，将上下移动值之外的值去除，用拟合值代替；③检验拟合的AVO曲线，是否与正演道集AVO曲线趋势一致，一致则输出去噪道集，否则调整参数值大小，重复①和②，输出去噪后道集。步骤2：采用非地表一致性剩余静校正实现道集拉平处理：①将一定偏移距内的道集进行叠加作为初始道；在一定时窗内，分别计算各道与初始道的相关系数，选择相关系数最大的一道作为模型道；②从某道开始依次在指定可移动时窗范围内们，以时窗t滑动求取与相邻道的相关系数；③将相关系数最大的滑动数作为当前的剩余时差，在指定时间段T内，各道减去剩余时差之后，输出拉平道集。步骤3：叠前远道集在经过去噪、拉平处理后，若仍然存在杂乱，道集数据缺损、噪音干扰严重等现象，需要对数据采取切除30°以外的道集，这样能保证叠前道集品质及有效信息，且减少了后续分析的时间。若不存在上述问题，则省略此步骤。(2)统计计算AVO异常值步骤1：对优化处理之后的叠前道集同一时间点上的振幅值进行多项式拟合，采用最小二乘法原理来实现；步骤2：根据均匀介质与各向异性介质在入射角非常小(约0°-6°之间)时其值基本相同，将道集各道的拟合值减去看作均匀介质的各道近道集(入射角约0°-6°之间)振幅值的平均值，得到各道的差值；步骤3：利用如下公式(1)：其中M＝1,2,3,4.......(1)(bi为原始值，bii为拟合值)统计各道的差值，计算标准差作为各道此时间点的异常值，异常值可看作是裂缝和流体综合影响的结果，可用该异常值有效地表征裂缝储层发育强度；步骤4：以采样率为步长，移动时间点，循环步骤1、步骤2、步骤3，得到整个地震数据的异常值数据体，通过剖面和平面成图，得到裂缝储层预测的图形效果，从而实现裂缝储层预测。进一步，多项式拟合，最小二乘法原理包括：1)在某时刻，若观测到的数据为一个入射角xi对应一个振幅值yi如下表所示：xix1x2……xnyiy1y2…….yn2)假设拟合函数不是线性函数，用如下的公式来进行拟合：y＝a0+a1x+a2x2+a3x3+.......+anxn,n＝0,1,2,3........；3)设拟合值和真实值之间的残差平方和为δi：4)当残差平方最小时，拟合曲线与真实值之间最为接近。所以，其关键问题就是求解多项式的系数。转化问题，关于多项式系数的函数，为得到最小值则求取其导数，令经过一系列的推导可以得到如下公式：求出ai，然后就可以进行拟合曲线。本发明对优化处理之后的叠前道集同一时间点上的振幅值进行多项式拟合；将道集各道的拟合值减去看作均匀介质的各道近道集(约0°-6°之间)振幅值的平均值，得到各道的差值；统计各道的差值计算标准差作为各道此时间点的异常值，异常值可看作是裂缝和流体的综合影响的结果，此方法相对于常规的基于P波各向异性地震属性椭圆拟合技术的计算量少，计算效率高；且不依赖于各向异性椭圆拟合的裂缝储层预测分析技术，解决宽窄方位角地震数据裂缝储层预测的问题。假设优化处理步骤一样，基于椭圆拟合的方法，假设也用异常值来拟合椭圆，需要全方位角度的计算，还要进行拟合椭圆，来预测裂缝，而本发明的方法，求得异常值，即可作为裂缝预测的依据。本发明的方法对某工区进行实际资料(CRP道集)分析，在井的周围都有裂缝的发育，与实际钻井吻合(应用本发明裂缝储层预测的符合率达到83％)；有效地验证了本发明的实用性。以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。当前第1页1 2 3