本发明涉及地震储层预测领域,更具体的说,它涉及一种基于分频相控的薄储层预测方法。
背景技术:
1、目前,关于薄储层识别预测,国内外学者进行了大量研究。总体上可分为两大类:基于地球物理属性的定性或半定量薄储层预测技术及地质综合描述单一岩溶改造作用薄储层预测技术。
2、基于地球物理属性的定性或半定量薄储层预测技术是根据原始地震数据隐含的大量地质信息,通过提取地震属性分析其反映的储层岩性、物性特征,进行储层的沉积、含油气性预测。由于地震属性种类繁多,提取的地震属性天差万别,如何提取、优化特色属性成为实现薄储层识别与预测的难题。
3、单一岩溶改造作用薄储层预测技术,以高分辨率地震资料处理及反演为主。由于薄储层研究需要在较高的分辨率的基础上开展工作,而现有地震资料往往频带窄、主频低,无法满足储层精细描述的精度。为了解决该类难题一般通过两方面进行,一方面是通过对原始资料进行拓频处理,但是该方法保幅保真存在不确定性,影响储层预测的可靠性;另一方面从地震资料数据采集阶段,通过高密度地震资料采集与处理提高地震资料成像精度,但周期较长,无法满足油藏开发快节奏的需求。
4、总而言之,上述每种薄储层预测方法都具有一定的可行性和局限性。因此,随着油田勘探开发的不断深入,急需形成一套适用于油藏开发快节奏的薄储层和非常规储层反演技术。
技术实现思路
1、针对现有技术存在的不足,本发明的目的在于提供一种基于分频相控的薄储层预测方法,其针对薄储层建立了在大数据波形驱动低频模型基础上开展变γ叠前确定性反演,进行相控叠前地质统计学反演组合,准确刻画薄砂层内部细节,解决了薄层油藏精细刻画识别难题。
2、为实现上述目的,本发明提供了如下技术方案:一种基于分频相控的薄储层预测方法,其包括如下步骤:
3、s1、对测井资料进行质控及一致性处理。对叠前地震资料进行质控和优化处理,提高叠前地震资料的精度。
4、s2、根据s1处理的资料,建立低频模型。
5、s3、以低频模型为基础,结合s1处理的资料,开展改进的变γ叠前联合确定性反演,获得低频、中频模型,公式如下:
6、vp(t)=vp(t0)exp[2∫0trp(τ)dτ]
7、vs(t)=vs(t0)exp[2∫0trs(τ)dτ]
8、vρ(t)=vρ(t0)exp[2∫0trρ(τ)dτ]
9、式中,t是时间采样点,t0是初始时间,τ=τ0(s-1),s为介质不均匀程度的参数,vp为纵波速度,vs为横波速度,vρ为密度,rp为纵波速度反射系数、rs为横波速度反射系数、rρ为密度反射系数。
10、rp、rs、rρ,通过如下公式获得:
11、rpp(θ)=(1+tan2θ)rp-8γsin2θrs+(1-4γsin2θ)rρ
12、式中,θ为入射角,γ为纵横波速度比,rpp(θ)为随着角度变化的纵波反射系数,θ、γ和rpp(θ)均根据s1处理的资料获得。
13、s4、根据反演结果开展砂体厚度刻画,制作平面砂地比图。
14、s5、根据s3反演结果,以平面砂地比图为约束进行地质统计学反演,将反演结果拓展到高频,获得全频带地质模型。
15、本发明进一步设置为:s2具体步骤为,首先,对s1结果中的地震波形进行大数据分析,获得低频数据。然后根据s1结果中的井点处纵波阻抗、纵横波速度比、密度,刻画地下介质的物理特性和结构特征,对低频数据进行限定,获得低频模型。
16、本发明进一步设置为:s2低频模型完成后需要分析低频模型是否符合地质认识,如果符合进行s3,如果不符合则返回s1。对于低频模型是否符合地质认识的判断,包括依据s1结果进行的核实和模型各部分一致性的检查。
17、本发明进一步设置为:s2建立低频模型时,需要根据已知的地震属性推算出未知的地震属性,不同井的同类型地震属性,通过空间转换公式获得,公式如下:
18、
19、式中,sn和si为不同井的地震属性,si根据s1结果获得。w为井道间关系函数,w通过s1结果推断获得,n为井的数量。
20、本发明进一步设置为:s2建立低频模型时,需要根据已知的地震属性推算出未知的地震属性,同一井不同类型地震物性属性,通过物性转换公式获得,公式如下:
21、
22、式中,lj和li为不同类型的地震物性属性,li根据s1结果获得。λij为属性相关匹配函数,λij通过s1结果推断获得,n为井的数量。
23、本发明进一步设置为:s2建立低频模型时,需要根据已知的地震属性推算出未知的地震属性,不同井不同类型的地震物性属性,通过全区物性参数转换公式获得,公式如下:
24、
25、式中,slj和sli为不同井不同类型的地震物性属性,wij为slj和sli对应的井道间关系函数,wij通过s1结果推断获得。λij为属性相关匹配函数,λij通过s1结果推断获得,n为井的数量。
26、本发明进一步设置为:s5中地质统计学反演采用变差函数,变差函数公式如下:
27、
28、式中,h为给定步长,n(h)为距离等于h的点对数,z(xi)为处于点xi处变量实测值,z(xi+h)为与点xi偏离h处变量的实测值,r(h)为反演结果。
29、本发明进一步设置为:s1中对测井资料的处理包括测井曲线误差校正和多井一致性处理。
30、本发明进一步设置为:测井曲线误差校正的步骤为:
31、对井筒具有明显变形的层段,采用castagna公式或smiths公式进行校正。
32、castagna公式:式中,p为地层密度。vp为纵波速度。a、b、c均为常数,依据井筒未变形层段的测井资料分析获得。
33、smiths公式:s(v)=dre,式中,s(v)为声波速度。r为电阻率。d、e均为常数,依据井筒未变形层段的测井资料分析获得。
34、对井筒具有明显垮塌的层段,采用测井响应的体积百分比模型进行校正,公式为:
35、式中,ci为第i种岩石成分未垮塌井段的测井响应值。vi为第i种岩石成分体积百分比。c为正演的测井响应值。
36、本发明进一步设置为:多井一致性处理公式为:式中,z(x,y)为所研究的地震属性,(x,y)是平面上某个点的坐标,为趋势值,g为校正值。
37、对测井资料中已知的z(x,y)进行回归分析,获得回归方程f(x,y),依据f(x,y)求出趋势值将z(x,y)与相减,求出g。
38、本发明进一步设置为:s1中对叠前地震资料的处理包括叠前随机噪声压制和时差校正。
39、本发明进一步设置为:叠前随机噪声压制采用geoeast系统进行三维叠前随机噪声衰减。
40、本发明进一步设置为:s1中对叠前地震资料采用非刚性匹配法进行时差校正。
41、本发明进一步设置为:平面砂地比图采用网格化插值法绘制。首先将平面砂地比图划分为一系列小网格,并将每个样点的砂土含量值分配到对应的网格节点上。然后根据已知数据点的空间位置和砂土含量值,采用插值方法计算每个网格节点的砂土含量值,绘制出砂土含量分布的等值线图。
42、本发明进一步设置为:所述网格化插值法为基于三角网格的等值线成图线性插值方法。
43、综上所述,本发明相比于现有技术具有以下有益效果:本发明针对薄储层建立了在大数据波形驱动低频模型基础上开展变γ叠前确定性反演,进行相控叠前地质统计学反演组合,准确刻画薄砂层内部细节,解决了薄层油藏精细刻画识别难题,并结合振幅属性、含油气检测等多方法综合定性定量相互验证,解决增产提速的地质问题,降低了勘探风险,具有普遍适用性和广阔的应用前景。
1.一种基于分频相控的薄储层预测方法,其特征在于,包括如下步骤:
2.根据权利要求1所述的一种基于分频相控的薄储层预测方法,其特征在于:s2具体步骤为,首先,对s1结果中的地震波形进行大数据分析,获得低频数据;然后根据s1结果中的井点处纵波阻抗、纵横波速度比和密度,刻画地下介质的物理特性和结构特征,对低频数据进行限定,获得低频模型。
3.根据权利要求2所述的一种基于分频相控的薄储层预测方法,其特征在于:s2低频模型完成后分析低频模型是否符合地质认识,如果符合,进行s3,如果不符合,则返回s1;对于低频模型是否符合地质认识的判断,包括依据s1结果进行的核实和模型各部分一致性的检查。
4.根据权利要求2所述的一种基于分频相控的薄储层预测方法,其特征在于:s2建立低频模型时,根据已知的地震属性推算出未知的地震属性,不同井的同类型地震属性通过空间转换公式获得,公式如下:
5.根据权利要求2所述的一种基于分频相控的薄储层预测方法,其特征在于:s2建立低频模型时,根据已知的地震属性推算出未知的地震属性,同一井不同类型地震物性属性通过物性转换公式获得,公式如下:
6.根据权利要求2所述的一种基于分频相控的薄储层预测方法,其特征在于:s2建立低频模型时,根据已知的地震属性推算出未知的地震属性,不同井不同类型的地震物性属性通过全区物性参数转换公式获得,公式如下:
7.根据权利要求1所述的一种基于分频相控的薄储层预测方法,其特征在于:s5中地质统计学反演采用变差函数,变差函数公式如下:
8.根据权利要求1所述的一种基于分频相控的薄储层预测方法,其特征在于:s1中对测井资料的处理包括测井曲线误差校正和多井一致性处理。
9.根据权利要求8所述的一种基于分频相控的薄储层预测方法,其特征在于:测井曲线误差校正的步骤为:
10.根据权利要求8所述的一种基于分频相控的薄储层预测方法,其特征在于:多井一致性处理公式为:式中,z(x,y)为所研究的地震属性,(x,y)是平面上某个点的坐标,为趋势值,g为校正值;
11.根据权利要求1所述的一种基于分频相控的薄储层预测方法,其特征在于:s1中对叠前地震资料的处理包括叠前随机噪声压制和时差校正。
12.根据权利要求11所述的一种基于分频相控的薄储层预测方法,其特征在于:叠前随机噪声压制采用geoeast系统进行三维叠前随机噪声衰减。
13.根据权利要求11所述的一种基于分频相控的薄储层预测方法,其特征在于:s1中对叠前地震资料采用非刚性匹配法进行时差校正。
14.根据权利要求1所述的一种基于分频相控的薄储层预测方法,其特征在于:平面砂地比图采用网格化插值法绘制;首先将平面砂地比图划分为一系列小网格,并将每个样点的砂土含量值分配到对应的网格节点上;然后根据已知数据点的空间位置和砂土含量值,采用插值方法计算每个网格节点的砂土含量值,绘制出砂土含量分布的等值线图。
15.根据权利要求14所述的一种基于分频相控的薄储层预测方法,其特征在于:所述网格化插值法为基于三角网格的等值线成图线性插值方法。