[0039] (Ⅲ)通过最小二乘法回归分析将地震属性与裂缝属性拟合后进行融合;
[0040] (Ⅳ)加入岩性属性,即裂缝影响因子数据,对融合体进行进一步约束融合,形成新的裂缝预测平面图。
[0041] 所述的方法,其中,步骤(Ⅲ)包括以下步骤:
[0042] ①加权值计算:获得n组样本数据(xi1,xi2,xi3;yi),i=1、2、……、n并根据公式:F’=a0+a1x1+a2x2+a3x3提取钻井所在位置裂缝密度的数值,将这些数值代入关系式中得a1,a2,a3,a0值,其中:F’为裂缝密度计算值,x1,x2,x3分别为AFE、分频相干和蚂蚁体属性值,a1,a2,a3为属性权重系数,a0为常数;
[0043] ②对a1,a2,和a3进行加权处理,使得a’1+a’2+a’3=1。
[0044] ③将步骤三至五得到的地震属性与述步骤得到的加权因子相乘并相加得到裂缝属性融合体,即a’1x1+a’2x2+a’3x3进行计算。
[0045] 所述的方法,其中,步骤(Ⅳ)包括:提取研究区目标层段的裂缝融合体属性,即计算目标时窗内属性体的RMS(均方根)值,形成一个二维面数据,将所述的二维面数据与步骤6计算得到的裂缝岩性影响因子ε数据相乘,最终得到研究区基岩裂缝预测分布图。
[0046] 通过研究确定基岩裂缝分布特点,落实裂缝型储层的储集性能,进一步揭示基岩油气藏的成藏规律,为基岩油气藏勘探潜力提供理论支撑,为继续寻找大中型油气田提供依据。
[0047] 通过本发明的研究,大大提高了基底裂缝预测的效率,能够较准确地确定裂缝分布特征,进一步明确基岩潜山勘探靶区。对加快油气产能建设的步伐具有重要的现实意义。该技术的研究成果将对中国东部其他盆地和区块基岩裂缝预测具有重要的借鉴作用。
附图说明
[0048] 在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中:
[0049] 图1是基于多属性融合的基底岩性裂缝预测流程图;
[0050] 图2是研究区基底岩性分布图;
[0051] 图3是研究区基底岩性裂缝预测分布图;
[0052] 在附图中,相同的部件使用相同的附图标记。附图并未按照实际的比例。
具体实施方式
[0053] 下面将结合附图对本发明作进一步说明。
[0054] 如图1所示,为基于多属性融合的基底岩性裂缝预测流程图,所述基于多属性融合的基岩裂缝预测方法,包括以下步骤:
[0055] 步骤(一):搜集研究区基岩资料,采集野外露头和钻井岩心样品,观察和统计野外露头和岩心裂缝信息,如密度、方向,充填程度等。
[0056] 其中,所述搜集的资料包括已有钻遇基岩的钻井和野外露头资料,尤其是钻井成像测井资料和裂缝密度曲线;对采集的钻井和露头样品开展镜下薄片鉴定,了解研究区的基岩裂缝微观特征和期次;研究总结目标区及其邻区大地构造背景,分析其应力环境和断裂分布特征,了解研究区基岩岩性分布情况。
[0057] 步骤(二):开展基底顶面层位标定,对研究区基底顶面界线进行地震解释工作,同时通过解释骨干剖面,搭建研究区的基岩地层格架。
[0058] 所述层位标定,是对研究区钻遇基底的井开展合成记录工作,同时确定基底顶面反射界面的地震反射特征。
[0059] 所述地震解释工作包括层位解释和断层解释,对基底顶面界线开展至少8*8 网格解释,然后对其网格化处理,为接下来的地震属性计算做好数据准备。
[0060] 步骤(三):计算Automatic fault Extraction(AFE)断裂预测属性。运用 AFE裂缝预测技术,对研究区尺度较大裂缝的展布进行预测。
[0061] 其中,步骤(三)所述AFE断裂预测技术是在三维地震不连续数据体上自动提取裂缝的技术。其技术实施步骤为:
[0062] 1)对目标区三维地震数据进行倾角,方位角计算,获得研究区基岩方位倾角信息;
[0063] 2)对原始地震数据进行倾角控制增强处理(Dip Steered Enhancement Attr),倾角控制增强模块能够改善地震资料品质以利于更好的层位、断层解释,在计算过程中充分考虑每一个采样点的倾角变化,并按计算样点为中心选取对比视窗,所参与计算的每一个采样点均定义相应的计算权度。
[0064] 3)不连续数据体解释,对2)得到的处理后的地震数据体,结合1)步得到的倾角,方位角信息,创建不连续数据体。
[0065] 4)对3)得到的不连续数据体做裂缝自动拾取计算,得到结果。
[0066] 步骤(四):计算分频相干属性。运用频谱分解方法,计算研究区高分辨率相干体,通过比较各个分频相干切片成像特征,结合实测数据,选择更加能突出小尺度裂缝的那个数据体。
[0067] 步骤(四)所述的频谱分解方法是薄层反射的调谐理论。频谱分解地震属性采用Gabor-Morley小波变换,将各地震道分解成不同的频带成分,有助于识别地震数据的局部岩性特征。
[0068] 步骤(四)所述的高分辨率相干体(Coherence Cube)是利用较最新的本征值算法,与过去的相关(Correlation)、相似(Semblance)相比,能更好的区分出断层细节特殊地质现象。
[0069] 步骤(四)所述的分频相干方法的原理是基于不同频带的地震数据可以检测不同尺度的断层及裂缝,其中的高频组分对微断层或微挠曲反映敏感,因此可以利用高频组分检测地层的横向不连续性(微断层或裂缝),帮助检测到利用全频段地震数据不易发现的小-微裂缝。
[0070] 所述的分频相干属性体主要实现流程为:
[0071] (1)通过分析找出有效地震频带范围,在该范围内,利用频率域滤波技术生成一系列单频三维地震数据体。
[0072] (2)对上述单频三维地震数据体分别进行计算,利用相同的本征值相干算法和处理参数生成对应的一系列相干数据体。再结合具体的地质情况,对这些频带的相干数据进行分析比较,选择更能表征小-微裂缝的分频属性体。
[0073] 步骤(五):计算蚂蚁追踪属性。对研究区开展蚂蚁追踪属性提取,高精度刻画研究区裂缝分布。
[0074] 所述蚂蚁体追踪技术是一种叠后微裂缝预测激素,能够有效识别小断层。所述的蚂蚁追踪属性体的生成流程为以下几个步骤:
[0075] 首先对地震数据体进行构造平滑处理,构造平滑的目的是采用中值滤波算法,消除或减弱地震数据的随机噪音干扰,即对地震数据进行“去噪”处理,使地震反射序列的连贯性得到增强。在对地震数据进行构造平滑处理时,通过设置主线方向、联络线方向和垂向的步长来确定平滑程度。另外,要根据实际情况决定是否选择倾角导向控制以及断层增强控制。在正常情况下,步长取1-2.5,不宜过大,而且一般勾选倾角导向控制和断层增强控制。
[0076] 然后生成方差体。通过对构造平滑后的地震体进行边缘增强处理,可以生成方差体,用来寻找边界异常。在追踪不连续点的过程中,利用方差体分析技术可以突出地震数据体中不连续的信用,并且方差体不连续性的增强效果明显优于其他地方。
[0077] 最后生成蚂蚁属性体。
[0078] 对地震资料进行预处理后,通过设定适当的追踪参数,运用蚂蚁算法对生成的方差体进行计算,提取蚂蚁追踪属性体。在这一个过程中涉及到6个关键参数,分别是初始蚂蚁分布边界、蚂蚁追踪背离、蚂蚁搜索步长、非法步长、合法步长和终止标准。
[0079] 对研究区的原始地震资料依次进行“去噪”处理,得到构造平滑属性体、方差属性体;在方差属性体基础上采用被动追踪的方式,调整6个参数,经过多次反复计算,得到蚂蚁追踪属性体。其中蚂蚁属性值的大小能够很好地表征裂缝的发育程度,通常属性值在-1至1之间。
[0080] 步骤(六):运用重磁震综合解释技术识别基底岩性,将获得的该区基底岩性分布图数据化。根据统计结果,计算不同岩性裂缝影响因子ε,生成影响因子ε平面数据体。
[0081] 所述的岩性识别方法具体步骤在此不再详述。
[0082] 步骤(六)所述的基底岩性分布图数据化是指将岩性图转化为离散化数据,即横坐标x、纵坐标y和岩性代码γ三列数据,例如:花岗岩γ赋值为1;千枚岩、板岩γ赋值为2;石英(片)岩γ赋值为3;火山岩γ赋值为4;碳酸盐岩γ赋值为5,;变余砂砾岩γ赋值为6等。具体赋值根据研究区实际情况和研究需要对γ进行赋值。
[0083] 所述裂缝影响因子ε是指不同岩性产生裂缝的难易程度,因子越大越容易产生裂缝,因子越小越不容易产生裂缝,其中
[0084] 在相同构造应力下,暗色矿物含量低的岩性(优势岩性)塑形较弱,容易产生裂缝成为储集层,具体分类见表1,序列越靠前越可能成为优势岩性。
[0085] 表1
[0086]
[0087] 本发明通过统计研究区及邻区岩性和裂缝对应关系,来获得该岩性的相对影响因子εi。具体方法如下:
[0088] Fi=εi(γi)+Ti,i=1,2,……,n;其中Fi为实测裂缝密度值,Ti为构造应力,εi(γi)为某岩性的裂缝影响因子。从表达式中可以看到,Fi与εi(γi) 成正相关,我们将实测裂缝密度与岩性代码带入公式,即可得到该研究区各个岩性的裂缝影响因子εi。
[0089] 将这个裂缝影响因子εi替换对应岩性γi,形成一个裂缝岩性影响因子εi平面数据体。
[0090] 步骤(七):基于地质属性的多属性融合WFGI计算。首先对上述步骤三至步骤五所述的属性进行去异常值和归一化处理,其次通过最小二乘法回归分析将地震属性与裂缝属性拟合后进行融合,然后加入岩性属性,即裂缝影响因子数据,对融合体进行进一步约束融合,形成新的裂缝预测平面图。
[0091] 步骤(七)所述的去除异常值处理是指去除一些属性运算中产生的异常值,具体方法为:当f’>fave+3σ或者f’
[0092] 步骤(七)中所述的归一化处理是一种无量纲处理手段,使物理系统数值的绝对值变成某种相对值关系。简化计算,缩小量值的有效办法。本发明采用的归一化方法是离差标准化和标准差标准化。
[0093] 离差标准化是对原始数据的线性变换,使结果值映射到[0,-1]之间,转换函数如下:y=(x-Min)/(Max-Min)。其中max为样本数据的最大值,min为样本数据的最小值。
[0094] 标准差标准化,经过处理的数据符合标准正态分布,即均值为0,标准差为 1,其转化函数为:y=(x-μ)/σ,其中μ是所有样本数据的均值,σ为所有样本数据的标准差。
[0095] 运用两种归一化处理手段主要是因为有时一种归一化方法会将异常平均化,因此需要对比这两种归一化方法获得的数据,选取异常区突出,且边界清晰的方法参与后续融合运算。
[0096] 步骤(七)所述的通过最小二乘法回归分析将地震属性与裂缝属性拟合后进行融合是指充分利用已有的测井和岩心裂缝资料,将地震属性与这些裂缝数据进行拟合确定每个属性的权值,然后将地震属性按照加权值融合在一起,形成一个新地震属性体的方法。具体步骤如下:
[0097] ①、加权值计算。假设有n次独立观测,得到n组样本数据(xi1,xi2,xi3;yi), i=1、2、……、n。设定一个关系式F’=a0+a1x1+a2x2+a3x3(F’为裂缝密度计算值,x1,x2,x3分别为AFE、分频相干和蚂蚁体属性值,a1,a2,a3为属性权重系数,a0为常数),提取钻井所在位置裂缝密度的数值,将这些数值代入关系式中,数据越多越好;
[0098] 应用最小二乘法原理,将实测值F与计算值F’的离差的平方和 最小优化判据。
[0099] 令
[0100]
[0101] 由于 是关于a1,a2,a3,a0的非负2次函数,比存在最小值,利用微积分的极值求法,得到:
[0102]
[0103] 这里 (i=0,1,2,……,p)是ai的最小二乘法估计,对上述求偏导,求得正规方程组的过程可用矩阵代数运算进行,得到正规方程组的矩阵表示:
[0104]
[0105] 移项得
[0106] 依据假定
[0107] R(X)=p+1=4,所以R(XTX)=R(X)=4,故(XTX)-1存在,解方程组得到:
[0108]
[0109] 最终求得a1,a2,a3,a0值。为了判断回归的效果好坏,还要计算相似系数R,剩余标准偏差S,R越趋于1越好,S越趋于0越好。
[0110] ②、对a1,a2,和a3进行加权处理,使得a’1+a’2+a’3=1。
[0111] ③、将步骤3-5得到的地震属性与述步骤得到的加权因子相乘并相加得到裂缝属性融合体,即a’1x1+a’2x2+a’3x3进行计算。
[0112] 步骤(七)所述加入岩性属性,对融合体进行进一步约束融合是指将新的裂缝属性融合体,与影响因子εi数据体进一步融合,形成最终的研究区基岩裂缝预测分布图,具体方法如下:
[0113] 首先提取研究区目标层段的裂缝融合体属性,即计算目标时窗内属性体的 RMS(均方根)值,形成一个二维面数据;然后将这个二维面数据与步骤(六) 计算得到的裂缝岩性影响因子ε数据相乘,最终得到研究区基岩裂缝预测分布图。
[0114] 实施例:
[0115] 采用东北LS断陷SN地区进行说明。如图1所示的震多属性融合的裂缝预测流程图,(一)搜集LS断陷SN地区基岩裂缝资料,采集野外露头和钻井岩心样品,观察和统计野外露头和岩心裂缝。根据实测结果编制了裂缝充填特征分布图、裂缝密度图等基础图件,分析了研究区裂缝的主要方向。(二)对LS断陷SN地区地震反射界面T5开展层位标定,开展地震解释工作,解释基底顶面界线T5,形成8×8网格,并对其网格化;(三)计算AFE断裂预测属性。运用AFE裂缝预测技术,对研究区尺度较大的裂缝展布进行预测。
[0116] 1)对目标区三维地震数据进行倾角,方位角计算,其中顶部和底部剖面的垂直采样率分别为10和20样本数,获得研究区基岩方位倾角信息;
[0117] 2)对原始地震数据进行倾角控制增强处理,倾角计算方式选择外飘型,y轴和x轴的生成平滑范围都为3。
[0118] 3)相干体解释,对上面步骤2得到的地震数据体,结合步骤1)得到的倾角,方位角信息,创建相干体,其中垂向视窗范围是59样本数。
[0119] 4)对得到的相干体做裂缝自动拾取(AFE)计算,参数去燥类型选择两个都要,得到结AFE属性体。
[0120] (四)计算分频相干体属性。运用频谱分解方法,计算研究区高分辨率相干体,通过比较各个分频相干切片成像特征,结合实测数据,选择更加能突出小尺度裂缝的那个数据体。
[0121] 分频相干属性体主要步骤:
[0122] (1)利用频率域滤波技术生成20Hz、30Hz、40Hz、50Hz、60Hz、70Hz和80Hz、一系列单频三维地震数据体。
[0123] (2)对上述单频三维地震数据体分别进行计算,利用相同的本征值相干算法和处理参数生成对应的一系列相干数据体。
[0124] (3)在地震解释的基础上,针对基岩,分别提取上述一系列相干数据体的沿层属性,再结合LS断陷SN地区的地质情况,对这些频带的相干数据进行分析,最终选择更能表征小-微裂缝的70Hz分频属性体。
[0125] (五)、计算蚂蚁追踪属性。对研究区开展蚂蚁追踪属性提取,精细刻画研究区裂缝分布。
[0126]
[0127] 表2
[0128] 蚂蚁追踪属性体的具体步骤是:
[0129] ①对地震数据体进行构造平滑处理,即对地震数据进行“去噪”处理,使地震反射序列的连贯性得到增强。在对地震数据进行构造平滑处理时,通过设置主线方向、联络线方向和垂向的步长来确定平滑程度。另外,要根据实际情况决定是否选择倾角导向控制以及断层增强控制。在正常情况下,步长取1-2.5,不宜过大,而且一般勾选倾角导向控制和断层增强控制。
[0130] ②生成方差体。通过对构造平滑后的地震体进行边缘增强处理,可以生成方差体,用来寻找边界异常。
[0131] ③生成蚂蚁属性体。在这一个过程中涉及到6个关键参数,分别是初始蚂蚁分布边界、蚂蚁追踪背离、蚂蚁搜索步长、非法步长、合法步长和终止标准。如表2所示:
[0132] 研究区裂缝走向和倾向的输入主要依据研究区内的成像测井资料以及岩心露头统计结果如图2所示。
[0133] 对梨树断陷的原始地震资料依次进行“去噪”处理,得到构造平滑属性体、方差属性体;在方差属性体基础上采用被动追踪的方式,调整6个参数,经过多次反复计算,得到蚂蚁追踪属性体。其中蚂蚁属性值的大小能够很好地表征LS断陷SN地区裂缝的发育程度。
[0134] (六)、如图2所示,运用重磁震综合解释技术识别出研究区岩性有花岗岩,火山岩,石英岩和千枚岩四种岩性,将岩性分布图数据化。计算不同岩性裂缝影响因子ε,生成影响因子ε平面数据体。
[0135] 将千枚岩γ赋值为1;石英(片)岩γ赋值为2;花岗岩γ赋值为3;火山岩γ赋值为4,将岩性分布图转化为横坐标x、纵坐标y和岩性代码γ三列数值的平面数据。根据前人研究和研究区和邻区钻井露头裂缝统计结果(下表),依据Fi=εi(γi)+Ti,i=1,2,……,n模型,计算得到这四种岩性的影响因子为千枚岩ε为0.22;石英片岩ε为0.16;花岗岩ε为0.32;
火山岩ε为0.3。将每种岩性ε值替换岩性符号γ,形成一个裂缝岩性影响因子εi平面数据体。
[0136] (七)、将步骤(三)-(五)所述的属性进行归一化,将三个属性体中当f’>fave+3 σ或者f’
[0137] 表3
[0138]
[0139] 提取JG1井裂缝曲线数值和JG101井岩心裂缝统计数值,运用最小二乘法对同一深度的三种地震属性与裂缝数值进行拟合,具体将这四组数据分别代入关系式F’=a0+a1x1+a2x2+a3x3,其中x1,x2,x3分别为AFE、分频相干和蚂蚁体属性值, F’为裂缝密度计算值。
[0140] 应用最小二乘法原理,将实测值F与计算值F’的离差的平方和 最小优化判据。通过一系列计算这三种属性的加权因子分别为a’1=0.23, a’2=0.46,a’3=
0.41。将三个属性体按照加权因子相乘后相加,得到新的裂缝融合体。对获得的融合体提取目标层位T5向下500ms时窗的RMS值,将获得平面数据与上述步骤(六)获得的影响因子数据相乘。最终形成入图3所示的基底岩性裂缝预测分布图
[0141] JG1井位于高异常区,所在高异常区范围较大,JG101井在中异常区边缘, JG2位于空白区,这都与已知的裂缝资料基本吻合。利用已知3口钻井进行检验,正确率在95%以上。
[0142] 虽然已经参考优选实施例对本发明进行了描述,但在不脱离本发明的范围的情况下,可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是,只要不存在结构冲突,各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文中公开的特定实施例,而是包括落入权利要求的范围内的所有技术方案。