CN111399056A - 一种基于划分方位角滤波预测裂缝强度的方法 - Google Patents
一种基于划分方位角滤波预测裂缝强度的方法 Download PDFInfo
- Publication number
- CN111399056A CN111399056A CN202010355818.0A CN202010355818A CN111399056A CN 111399056 A CN111399056 A CN 111399056A CN 202010355818 A CN202010355818 A CN 202010355818A CN 111399056 A CN111399056 A CN 111399056A
- Authority
- CN
- China
- Prior art keywords
- attribute
- frequency
- fracture
- stack
- coherent
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/512—Pre-stack
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于划分方位角滤波预测裂缝强度的方法,所属预测裂缝强度的方法至少包括如下步骤:S1:统计待测区域裂缝发育方向并确定中心角和方位角范围,根据方位角范围对输入道集数据提取中心角道集数据并进行叠前偏移处理,从而得到第一属性道集数据体;S2:对步骤S1中获得的第一属性道集数据体进行针对裂缝的优势频率分析,并设定频带范围后进行带通滤波处理,得到第二属性道集数据体;S3:对步骤S2获得的第二属性道集数据体进行叠前或叠后属性计算,获得若干用于预测不同规模裂缝发育强度的属性数据体。通过利用这些属性数据体分别对不同裂缝规模进行分析,从而实现准确预测地下裂缝储层的目的,更好为油气勘探提供服务。
Description
技术领域
本发明属于地震勘探数据处理领域,尤其涉及一种基于划分方位角滤波预测裂缝强度的方法。
背景技术
裂缝是地下油气聚集和运移的重要通道,在油气储层中发育裂缝时对于油气的沟通及采集具有积极意义。例如,四川盆地的须家河组致密砂岩及长兴组礁滩、海相页岩储层中均存在裂缝,钻遇裂缝储层的大部分钻井经测试均可获得高产工业气流,而没有钻遇裂缝储层的钻孔测试大部分获得低产工业气流。这样的勘探实践证明,裂缝预测在油气勘探中具有重要意义。当然,裂缝的发育强度与油气产能具有一定的联系,寻找储层中的强裂缝发育区域是油气勘探研究者的工作之一。
目前利用地震资料进行裂缝预测大致分为横波方法和纵波方法。其中,横波方法主要利用横波走时差异和垂直入射反射振幅来进行裂缝预测,理论上比纵波对裂缝产生的方位各向异性具有更强的敏感性,预测精度较高,但是由于多分量勘探和横波勘探的成本昂贵且高质量的横波资料获取难度大,推广应用受到限制;多分量转换横波裂缝检测在技术上比较复杂,应用还有待完善。
纵波方法又分为叠后与叠前方法,叠后方法主要通过检测地震波的不连续性来预测裂缝,如相干技术、曲率技术、蚂蚁体技术、相位属性等,简单实用,但是属于大尺度的定性预测技术。叠前方法基于方位各向异性理论,认为裂缝的各向异性导致地震波的各向异性,通过计算地震波不同方向的属性差异来预测裂缝的方向及相对密度,是一种半定量的裂缝预测技术,且成本较低,因此得到广泛推广和应用。
针对裂缝的预测,有相当多的专利技术文献,由此可见对裂缝的探测一直是研究及探索的热点。一些专利包括如发明《用于精细断层解释的优势频带相干处理方法》(专利号:200410058167.X)的专利公开了一种相干处理预测小断层及裂缝,经过优势频带相干技术处理后的地震剖面,显示出新的遗漏的小断层,小断层的断点位置更加准确,伸展方向更加连续;发明《裂缝预测方法和装置》(专利号:201010205983.4)的专利提出利用拾取的目标层时窗获取每个地震道的反射振幅,通过获取的方位角及反射振幅来进行椭圆拟合,以此来预测裂缝的方向及裂缝密度。但是,这些方法实现起来难度相当大,操作起来易受各种因素影响,主要存在如下几个问题:
(1)采用常规全方位的叠后地震资料进行裂缝预测,由于在叠加中存在平均效应,对预测具有方向性的裂缝型储层不利。
(2)叠前预测如果使用全频段的地震资料,往往会造成一些中、高频段受到压制,而这些受压制的频段往往与中、小规模裂缝的响应相关。
发明内容
本发明的目的在于,为了克服上述常规技术方法的不足,解决对中、小规模的裂缝储层探测问题,进而寻找到有利的储层位置。本发明公开了一种基于划分方位角滤波预测裂缝强度的方法,实现利用现有的三维地震资料采集方便、简单及准确地对地下深层中、小规模级别的裂缝储层进行有效预测。
本发明目的通过下述技术方案来实现:
一种基于划分方位角滤波预测裂缝强度的方法,所属预测裂缝强度的方法至少包括如下步骤:S1:统计待测区域裂缝发育方向并确定中心角和方位角范围,根据方位角范围对输入道集数据提取中心角道集数据并进行叠前偏移处理,从而得到第一属性道集数据体;S2:对步骤S1中获得的第一属性道集数据体进行针对裂缝的优势频率分析,并设定频带范围后进行带通滤波处理,得到第二属性道集数据体;S3:对步骤S2获得的第二属性道集数据体进行叠前或叠后属性计算,获得若干用于预测不同规模裂缝发育强度的属性数据体。
根据一个优选的实施方式,所述步骤S1中所述统计待测区域裂缝发育方向是指根据地质、测井、构造应力场分析所得到研究区目的层的裂缝发育方向。
根据一个优选的实施方式,所述步骤S1中输入道集数据是指经过保真去噪、保幅能量补偿、反褶积和静校正、动校正处理后的道集。
根据一个优选的实施方式,所述步骤S2中优势频率分析为根据实际研究区目的层的裂缝响应情况从而确定其相关的优势频率,通过分析并识别出不同的裂缝规模所对应的优势频率,建立并得到优势频率数据集。
根据一个优选的实施方式,所述不同的裂缝规模所对应的优势频率,能够通过如下方法识别确定:确定研究区内不同的裂缝发育情况并在沿层相干属性平面图上确定其代表性的样点;利用叠后地震数据体进行分频处理,得到若干单频地震数据体,并基于相干体计算,对各个相干体提取沿层属性的平均值,按所述代表性的样点的位置提取各单频相干属性值;对代表性的样点的各个单频相干属性值进行比较,确定出各样点中最小的相干值,该相干值所对应的单频体的频率值就是该样点的优势频率值。
根据一个优选的实施方式,所述不同的裂缝规模所对应的优势频率,能够通过如下方法识别确定:对井中的裂缝强度进行分级别,确定出不同规模的裂缝强度,并根据井震合成记录标定,将裂缝段由深度域转换到时间域;利用叠后地震数据进行分频处理后,再进行相干处理,得到一系列的单频相干数据体;利用转换到时间域的裂缝段对各个单频相干体提取井中裂缝段的图像,当裂缝图像中清晰度最高并且整体上该裂缝段相干属性值的平均值最小时,该单频相干体所对应的频率值即为优势频率值。
根据一个优选的实施方式,所述步骤S2中的带通滤波处理是指根据频带范围、相关滤波器进行频率域的滤波处理。
根据一个优选的实施方式,所述步骤S3中叠前计算为对第二属性道集数据体进行相关的叠前反演计算;所述叠前反演计算包括但不限于AVO、FVO计算。
根据一个优选的实施方式,所述步骤S3中叠后属性计算为对第二属性道集数据进行叠加处理后得到的叠后数据体,并利用叠后属性计算方法对叠后数据体进行属性数据体计算;所述叠后属性计算方法包括但不限于进行相干、曲率和相位计算。
前述本发明主方案及其各进一步选择方案可以自由组合以形成多个方案,均为本发明可采用并要求保护的方案;且本发明,(各非冲突选择)选择之间以及和其他选择之间也可以自由组合。本领域技术人员在了解本发明方案后根据现有技术和公知常识可明了有多种组合,均为本发明所要保护的技术方案,在此不做穷举。
本发明的有益效果:本发明采取分方位角范围及频率滤波技术,利用地震数据中特定方位角范围的中、高频段对中、小规模的裂缝储层响应比低频数据强,在处理过程中对某些频段进行相应的带通滤波处理,尽量突出对不同规模裂缝型储层的响应频带数据,通过所得到的第二属性道集数据体实施反演或属性计算,得到一系列的属性数据体。利用这些属性数据体分别对不同裂缝规模进行分析,从而达到准确预测地下裂缝储层的目的,更好为油气勘探提供服务。
附图说明
图1是本发明实施流程示意图;
图2是本次发明技术的裂缝走向方位角及中心角关系示意图;
图3是本次发明技术的中心角与最大方位角、最小方位角关系示意图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
需要说明的是,为使本发明实施例的目的、技术方案和优点更加清楚,下面对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。
因此,以下对本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参考图1所示,本发明公开了一种基于划分方位角滤波预测裂缝强度的方法。所属预测裂缝强度的方法至少包括如下步骤:
步骤S1:统计待测区域裂缝发育方向并确定中心角和方位角范围,根据方位角范围对输入道集数据提取中心角道集数据并进行叠前偏移处理,从而得到第一属性道集数据体。
步骤S2:对步骤S1中获得的第一属性道集数据体进行针对裂缝的优势频率分析,并设定频带范围后进行带通滤波处理,得到第二属性道集数据体。
步骤S3:对步骤S2获得的第二属性道集数据体进行叠前或叠后属性计算,获得若干用于预测不同规模裂缝发育强度的属性数据体。
实施例1
其中,针对步骤S1。统计裂缝发育方向并确定中心角、方位角范围,根据方位角范围对输入道集数据提取中心角道集数据并进行叠前偏移处理,从而得到第一属性道集数据体。
包括以下步骤:
步骤ⅰ.统计目的层段的裂缝发育方向,裂缝方向可以有一组或多组。该步骤主要是根据地质、测井、构造应力场分析所得到的研究区的主要裂缝发育方向。具体上可以根据叠后地震资料进行目的层的构造应力分析,进行目的层裂缝发育方向统计,也可以根据井中裂缝发育方向及相关地质认识来统计研究区的主要裂缝发育方向。主要裂缝发育方向可以为一组或多组,具体应视实际裂缝预测要求及精度确定。本发明中的方向均是指以正北方向为0°,顺时针旋转一周的方位角大小为360°;鉴于对称原理,可将所有方位角及中心角的数据值归于0-180°范围。
步骤ii.确定中心角值的大小,根据主要裂缝方向进行数学运算——加、减90°,即得到对应的中心角数据值。数学计算公式计算如下:
步骤ⅲ.确定提取道集的方位角范围。该步骤主要是指在确定中心角数据值后,以该中心角数据值进行加、减后所得到的方位角范围。值的大小应根据地震的方位角资料、覆盖次数及预测精度所确定,为设计的方位角范围值的一半。该方位角的数据范围的数学计算公式如下:
式(2)、(3)中,和分别表示第i个方位角范围内的方位角的最大值和最小值,表示第i个中心角数据值,为第i个方位角所设定方位角范围。输入道集数据是指经过保真去噪、保幅能量补偿、反褶积和静校正、动校正等处理后的道集数据,要求道集数据具有无明显噪音和静态时差。依次类推,完成各个中心角的方位角范围的确定。一般情况下,小于25°。
步骤iv、根据各个方位角范围进行道集数据的提取并进行叠前偏移处理,从而得到第一属性道集数据体。具体为对分方位角后的道集数据,利用相关的速度进行道集上的叠前时间偏移处理,得到叠前偏移处理后的道集数据。
其中,步骤S2:对步骤S1中获得的第一属性道集数据体进行针对裂缝的优势频率分析,并设定频带范围后进行带通滤波处理,得到第二属性道集数据体。包括以下步骤:
步骤i.对不同规模的裂缝进行优势频率分析,并设定相对应的频带范围。优势频率分析是指根据实际研究区目的层的裂缝级别情况确定其对应的优势频率,分析并识别出不同的裂缝规模所对应的一系列优势频率,建立并得到相关的优势频率数据集。通过确定不同裂缝级别的优势频率数据集及设计频率固定值(△d),经公式计算就能设定出不同裂缝发育规模的频带范围,而这些频带范围能对相对应的不同规模级别的裂缝具有相对敏感的响应。
其中优势频率的确定有如下几种方法,现分别阐述如下:
a、当井资料相对较少时,可以采用如下方法求取优势频率值:
(1)确定研究区内不同的裂缝发育情况并在沿层相干属性平面图上确定其代表性的样点。要求对研究区设定不同的裂缝发育规模及级别,裂缝级别可根据实际情况及井资料显示、专家经验对其进行分级,不同级别的裂缝可以在相干属性平面图上选择其有代表性的样点。要求不同裂缝级别的样点数要大于或等于3个。相干属性计算可以采用确定好的相干计算方法(如第二代相干算法)及其参数进行,沿目的层开固定时窗提取相干属性的平均值。原则上,裂缝强度大则相干值越小,反之,则相干值越大。一般情况下,不同级别的样点数越多,则对频带范围的确定相对有利;反之,则对频带范围的确定的精度降低。
具体地,本实施例中,相干计算方法优选地采用C2相干算法来对地震数据体进行相干体计算。在利用C2相干体算法来对叠后地震数据体进行相干体计算的过程中,首先对C2相干算法中的分析时窗内中心点在局部反射面x方向和y方向的视倾角p和q进行相应测试,该参数测试以达到能够清晰描述测试区域(测试区域只是研究区域中的一小块)内的较大规模的裂缝形态为最好。待测试完毕后,便利用测试确定出的视倾角p和q对所有中心角叠后数据体进行相干体计算。
C2相干体算法主要是首先定义一个以分析点为中心并且含有J道地震数据的矩形或椭圆形分析时窗,将当地的坐标轴作为分析点的中心,因此相似系数σ(t,p,q)便可以根据如下表达式计算得到:
其中,下标j表示落在分析时窗内的第j道,xj和yj分别表示第j道与分析时窗内的中心点t在x方向和y方向的距离,K表示垂直分析时窗,H表示希尔伯特变换,p和q分别表示分析时窗内中心点t所在局部反射界面x方向和y方向的视倾角,uj表示分析时窗内第j道的地震数据。
根据表达式(4)可以看出,通过改变视倾角p和q的取值可以调整局部平面的倾斜方向,从而对应得到一个相似系数。通过查找最大的相似系数,可以使得局部平面达到与实际反射界面的最佳拟合,因此此时的像素系数也即为分析点的相干性估计值(即相干数据值)。
需要指出的是,在本发明的其他实施例中,该方法还可以在步骤a(1)中采用其他合理的算法(例如C1相体干算法、C3相干体算法或新一代相干体算法)来对叠后地震数据体进行相干体计算,本发明不限于此。
(2)利用叠后地震数据体进行分频处理,得到一系列单频地震数据体,并利用上述(1)中的相干算法及参数进行单频地震数据体的相干体计算,并对各个相干体提取沿层属性的平均值,按上述(1)中的样点位置提取各个相干属性值。分频处理一般以频率的整数频率值开始、及某一整数频率值结束,以一定的整数间隔进行计算,得到一系列的分频地震数据体。如在实际进行分频处理,以20hz开始,以2hz作为增量,到70hz结束,得到一系列单频地震数据体。
(3)对样点的各个单频相干属性值进行比较,确定出各个样点的最小相干值,该相干值所对应的单频体的频率值就是该样点的优势频率值。依次类推,完成各个样点的优势频率值。
b、当井资料相对较多时,可以采用如下方法求取优势频率值:
(1)对井中的裂缝强度进行分级别,确定出不同规模的裂缝强度。并根据井震合成记录标定,将裂缝段由深度域转换到时间域。
(2)利用叠后地震数据进行分频处理后,再进行相干处理,得到一系列的单频相干数据体。利用转换到时间域的裂缝段对各个单频相干体提取井中裂缝段的图像,当裂缝的相干图像相对最为清晰且整体上该裂缝段相干属性值的平均值最小时,该单频相干体所对应的频率值即为优势频率值。
c、也可以综合采用步骤a、b中相关的计算方法,求取不同裂缝级别的优势频率值。
所述的带通滤波处理是指根据确定的频带范围、相关滤波器进行滤波处理,主要对地震信号的某些频段进行滤除,保留感兴趣的某些频带的滤波方法。
滤波的方法主要包括三角滤波或奥姆斯比带通滤波等频率域滤波方法。频带范围设定的具体操作为对工区内某一裂缝发育级别的第i个裂缝段所对应的优势频率进行合并后建立优势频率数据集{fi}(i=l,2,…,n),其中n为优势频率的个数;从优势频率数据集中确定出最小频率值,也就是取fmin=min{fi};从优势频率数据集中确定出最大频率值,也就是取fmax=min{fi}。根据fmin、fmax及设计频率固定值(△d)经计算得到该优势频带的频带范围。频带范围的起、止频率的具体计算公式如下:
f止=fmax+△d (5)
f起=fmin-△d (6)
式(5)、(6)中f止为滤波后频段的止频率,f起为滤波后频段的起频率,△d为设计频率固定值,fmax为优势频率数据集中的最大频率值,fmin为优势频率数据集中的最小频率值。△d一般设定为2hz,不同的裂缝发育规模可以根据实际情况及专家经验,对不同的裂缝发育强度进行划分,如可以根据裂缝发育密度划分出强裂缝发育、较强裂缝发育、中等裂缝发育、弱裂缝发育等四种。也可以根据井中的裂缝发育情况,将其划分为不同裂缝发育强度级别。
其中,步骤S3:对步骤S2获得的第二属性道集数据体进行叠前或叠后属性计算,获得若干用于预测不同规模裂缝发育强度的属性数据体,包括以下步骤:
步骤i.对第二属性道集数据体进行叠前或叠后属性计算。叠前或叠后属性计算是指对第二属性道集数据体进行相关的叠前反演(如AVO、FVO等);或者对第二属性道集数据进行叠加处理后得到的叠后数据体,并利用相关的属性计算方法(相干、曲率、相位等)对叠后地震数据体进行属性数据体计算。在实际计算中,可以根据井中资料的实际情况及专家经验,确定出一种适合于研究区的裂缝预测方法。一般情况下,可以选择相干算法或P波各向异性处理,但P波各向异性处理需要有三个及以上的中心角数据。也可以根据实际情况,当某一裂缝级别的属性数据体有三个及以上时,可以有选择的进行颜色融合计算,得到一个用于评价该裂缝级别的属性数据体。
步骤ii.利用确定的相关的反演或属性计算方法对第一属性数据体进行相关计算,得到一系列的属性数据体,从而根据这些数据体确定研究区的不同方向及发育强度的裂缝预测及评价。
本发明应用案例:
在川东北某研究区的油气勘探中发现,须四段致密砂岩储层具有裂缝强发育的特征,打在裂缝段中的钻井经测试获得高产工业气流。因此,在致密砂岩中寻找裂缝发育的区域并对其发育强度进行评价,是该区油气勘探相当重要的研究方向。
通过步骤S1,利用三维叠后地震资料、测井数据、层位数据、层速度数据及FRS软件进行构造应力场分析,从而基本上确定研究区目的层致密砂岩段的主要裂缝走向。该区的主要裂缝走向共有两组,分别为45°、135°,呈共轭状。因此,经相关计算可以得到的地震中心角共有两个,分别为45°、135°,如图2所示。中心角45°的地震数据主要负责对走向为135°的裂缝进行清晰成像,而中心角135°的地震数据主要负责对走向为45°的裂缝进行清晰成像。案例中设定均为15°,则两个方位角范围分别为30°~60°、110°~150°,如图3所示。根据相关参数提取这两个方位角范围的道集数据体,并对道集数据进行叠前偏移处理,得到两个第一属性道集数据体。
通过步骤S2。对研究区井中目的层的裂缝进行级别划分及优势频率分析,确定不同裂缝强度所对应的优势频率值,从而得到相关的不同裂缝强度级别的优势频率数据集。
根据优势频率数据集及所设计的频率固定值,经计算可得到不同裂缝级别频带范围。在实例中,将研究区分为三个裂缝发育强度级别,分别为裂缝强发育、较强裂缝发育、弱裂缝发育。
并根据步骤S2中,设计相关的单频体及相干体计算,实际中设定单频地震数据体的起算频率为20hz,到频率70hz结束,频率增量为2hz,得到一系列单频地震数据体;并对这些数据体进行相干体计算,采用C2相干计算方法。在实例中,对三种裂缝级别有代表性的设定样点,样点数共计为36个,分别为每种裂缝级别为12个,并提取其相关目的层段的单频相干属性值,利用各个样点的最小相干值以确定各个样点的优势频率。实例中经相关计算方法确定出研究区裂缝强发育的优势频带范围为32hz~42hz,并设定设计频率固定值为2hz,从而确定该裂缝级别的频带范围为30hz~44hz;中等裂缝发育的优势频带范围为44hz~52hz,并设定设计频率固定值为2hz,从而确定该裂缝级别的频带范围为42hz~54hz;;弱裂缝发育的优势频带范围为56hz~62hz,并设定设计频率固定值为2hz,从而确定频带范围为54hz~64hz。
利用计算得到不同裂缝发育强度级别的频带范围及采用带通滤波中的奥姆斯比带通滤波,在相关频率域中分别对两个对第一属性道集数据体进行带通滤波,得到六个第二属性道集数据体。
通过步骤S3。对六个第二属性道集数据体分别进行叠加,并确定采用C2相干计算方法进行计算。实例中对得到的六个叠加数据体分别采用同一相干计算方法及参数实施相干计算,所得到的六个相干数据体分别对研究区的两个主要裂缝方向(45°、135°)实施不同裂缝强度级别的预测及评价。经对比分析发现,利用本发明所得到的成果优于常规叠后相干技术所得到的成果,预测砂层中、小规模的裂缝精度及清晰度都比常规的相干技术高。并与现有钻井井中的裂缝进行验证,吻合率可达到81.7%。
利用本发明技术所预测的裂缝储层分布状态,由于实施方位角范围及带通滤波处理,使其更能突出不同方向及强度的裂缝储层的响应,利于对裂缝储层的分布状态进行分析、研究,所得成果优于常规裂缝预测技术所取得的预测成果。况且利用本发明技术更能清楚地实施对四川盆地海相礁滩储层及致密砂岩、页岩气储层的探测,均取得良好的预测效果,与相关地区的实钻钻井资料吻合度较高。
本说明书中公开的所有特征,或公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合。
本说明书中公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换。即,除非特别叙述,每个特征只是一系列等效或类似特征中的一个例子而已。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (9)
1.一种基于划分方位角滤波预测裂缝强度的方法,其特征在于,所属预测裂缝强度的方法至少包括如下步骤:
S1:统计待测区域裂缝发育方向并确定中心角和方位角范围,根据方位角范围对输入道集数据提取中心角道集数据并进行叠前偏移处理,从而得到第一属性道集数据体;
S2:对步骤S1中获得的第一属性道集数据体进行针对裂缝的优势频率分析,并设定频带范围后进行带通滤波处理,得到第二属性道集数据体;
S3:对步骤S2获得的第二属性道集数据体进行叠前或叠后属性计算,获得若干用于预测不同规模裂缝发育强度的属性数据体。
2.如权利要求1所述的一种基于划分方位角滤波预测裂缝强度的方法,其特征在于,所述步骤S1中所述统计待测区域裂缝发育方向是指根据地质、测井、构造应力场分析所得到研究区目的层的裂缝发育方向。
3.如权利要求1所述的一种基于划分方位角滤波预测裂缝强度的方法,其特征在于,所述步骤S1中输入道集数据是指经过保真去噪、保幅能量补偿、反褶积和静校正、动校正处理后的道集。
4.如权利要求1所述的一种基于划分方位角滤波预测裂缝强度的方法,其特征在于,所述步骤S2中优势频率分析为根据实际研究区目的层的裂缝响应情况从而确定其相关的优势频率,通过分析并识别出不同的裂缝规模所对应的优势频率,建立并得到优势频率数据集。
5.如权利要求4所述的一种基于划分方位角滤波预测裂缝强度的方法,其特征在于,所述不同的裂缝规模所对应的优势频率,能够通过如下方法识别确定:
确定研究区内不同的裂缝发育情况并在沿层相干属性平面图上确定其代表性的样点;
利用叠后地震数据体进行分频处理,得到若干单频地震数据体,并基于相干体计算,对各个相干体提取沿层属性的平均值,按所述代表性的样点的位置提取各单频相干属性值;
对代表性的样点的各个单频相干属性值进行比较,确定出各样点中最小的相干值,该相干值所对应的单频体的频率值就是该样点的优势频率值。
6.如权利要求4所述的一种基于划分方位角滤波预测裂缝强度的方法,其特征在于,所述不同的裂缝规模所对应的优势频率,能够通过如下方法识别确定:
对井中的裂缝强度进行分级别,确定出不同规模的裂缝强度,并根据井震合成记录标定,将裂缝段由深度域转换到时间域;
利用叠后地震数据进行分频处理后,再进行相干处理,得到一系列的单频相干数据体;利用转换到时间域的裂缝段对各个单频相干体提取井中裂缝段的图像,当裂缝图像中清晰度最高并且整体上该裂缝段相干属性值的平均值最小时,该单频相干体所对应的频率值即为优势频率值。
7.如权利要求1所述的一种基于划分方位角滤波预测裂缝强度的方法,其特征在于,所述步骤S2中的带通滤波处理是指根据频带范围、相关滤波器进行频率域的滤波处理。
8.如权利要求1所述的一种基于划分方位角滤波预测裂缝强度的方法,其特征在于,所述步骤S3中叠前计算为对第二属性道集数据体进行相关的叠前反演计算;所述叠前反演计算包括但不限于AVO、FVO计算。
9.如权利要求1所述的一种基于划分方位角滤波预测裂缝强度的方法,其特征在于,所述步骤S3中叠后属性计算为对第二属性道集数据进行叠加处理后得到的叠后数据体,并利用叠后属性计算方法对叠后数据体进行属性数据体计算;所述叠后属性计算方法包括但不限于进行相干、曲率和相位计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010355818.0A CN111399056B (zh) | 2020-04-29 | 2020-04-29 | 一种基于划分方位角滤波预测裂缝强度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010355818.0A CN111399056B (zh) | 2020-04-29 | 2020-04-29 | 一种基于划分方位角滤波预测裂缝强度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111399056A true CN111399056A (zh) | 2020-07-10 |
CN111399056B CN111399056B (zh) | 2020-12-08 |
Family
ID=71433840
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010355818.0A Active CN111399056B (zh) | 2020-04-29 | 2020-04-29 | 一种基于划分方位角滤波预测裂缝强度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111399056B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111897008A (zh) * | 2020-08-07 | 2020-11-06 | 西南石油大学 | 一种基于地震分频技术的断裂分级预测方法 |
CN113391350A (zh) * | 2021-04-30 | 2021-09-14 | 成都北方石油勘探开发技术有限公司 | 一种半定量叠后地震裂缝预测方法 |
CN113466936A (zh) * | 2021-07-01 | 2021-10-01 | 中海石油(中国)有限公司上海分公司 | 一种断层阴影区crp道集获取方法、装置、设备及存储介质 |
CN113671566A (zh) * | 2021-08-16 | 2021-11-19 | 西南石油大学 | 一种基于深度域地震数据计算裂缝参数的方法 |
CN113945964A (zh) * | 2020-07-16 | 2022-01-18 | 中国石油天然气股份有限公司 | 叠前地震数据多频段相干值获取方法及断裂预测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6055482A (en) * | 1998-10-09 | 2000-04-25 | Coherence Technology Company, Inc. | Method of seismic signal processing |
CN103076623A (zh) * | 2011-10-25 | 2013-05-01 | 中国石油化工股份有限公司 | 一种基于叠前相干的裂缝检测方法 |
CN104297784A (zh) * | 2014-08-12 | 2015-01-21 | 中国石油化工股份有限公司 | 一种基于地震纵波方位各向异性的裂缝预测方法 |
CN107728214A (zh) * | 2016-08-11 | 2018-02-23 | 中国石油化工股份有限公司 | 一种裂缝预测方法 |
CN108931814A (zh) * | 2017-05-26 | 2018-12-04 | 中国石油化工股份有限公司 | 一种基于多属性融合的基岩裂缝预测的方法 |
CN109799531A (zh) * | 2019-03-20 | 2019-05-24 | 成都理工大学 | 一种基于地震分频相干属性的裂缝储层预测方法 |
-
2020
- 2020-04-29 CN CN202010355818.0A patent/CN111399056B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6055482A (en) * | 1998-10-09 | 2000-04-25 | Coherence Technology Company, Inc. | Method of seismic signal processing |
CN103076623A (zh) * | 2011-10-25 | 2013-05-01 | 中国石油化工股份有限公司 | 一种基于叠前相干的裂缝检测方法 |
CN104297784A (zh) * | 2014-08-12 | 2015-01-21 | 中国石油化工股份有限公司 | 一种基于地震纵波方位各向异性的裂缝预测方法 |
CN107728214A (zh) * | 2016-08-11 | 2018-02-23 | 中国石油化工股份有限公司 | 一种裂缝预测方法 |
CN108931814A (zh) * | 2017-05-26 | 2018-12-04 | 中国石油化工股份有限公司 | 一种基于多属性融合的基岩裂缝预测的方法 |
CN109799531A (zh) * | 2019-03-20 | 2019-05-24 | 成都理工大学 | 一种基于地震分频相干属性的裂缝储层预测方法 |
Non-Patent Citations (1)
Title |
---|
仇鹏 等: "裂缝性储层的地震响应特征――以川中地区致密砂岩储层为例", 《天然气与石油》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113945964A (zh) * | 2020-07-16 | 2022-01-18 | 中国石油天然气股份有限公司 | 叠前地震数据多频段相干值获取方法及断裂预测方法 |
CN111897008A (zh) * | 2020-08-07 | 2020-11-06 | 西南石油大学 | 一种基于地震分频技术的断裂分级预测方法 |
CN111897008B (zh) * | 2020-08-07 | 2022-03-01 | 西南石油大学 | 一种基于地震分频技术的断裂分级预测方法 |
CN113391350A (zh) * | 2021-04-30 | 2021-09-14 | 成都北方石油勘探开发技术有限公司 | 一种半定量叠后地震裂缝预测方法 |
CN113391350B (zh) * | 2021-04-30 | 2023-10-13 | 成都北方石油勘探开发技术有限公司 | 一种半定量叠后地震裂缝预测方法 |
CN113466936A (zh) * | 2021-07-01 | 2021-10-01 | 中海石油(中国)有限公司上海分公司 | 一种断层阴影区crp道集获取方法、装置、设备及存储介质 |
CN113466936B (zh) * | 2021-07-01 | 2023-07-18 | 中海石油(中国)有限公司上海分公司 | 一种断层阴影区crp道集获取方法、装置、设备及存储介质 |
CN113671566A (zh) * | 2021-08-16 | 2021-11-19 | 西南石油大学 | 一种基于深度域地震数据计算裂缝参数的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111399056B (zh) | 2020-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111399056B (zh) | 一种基于划分方位角滤波预测裂缝强度的方法 | |
CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
CN101551463B (zh) | 三维观测系统噪声压制估算方法 | |
CN102109611B (zh) | 一种利用低频地震属性预测油藏优质储层的方法 | |
WO1997033183A1 (en) | Method for geophysical processing and interpretation using instantaneous phase and its derivatives and their derivatives | |
CN105445800A (zh) | 一种厚层砂体顶部分异岩性油藏的识别方法 | |
CN110231652B (zh) | 一种基于密度的含噪声应用空间聚类的地震相提取方法 | |
CN103777245B (zh) | 基于地震资料的油气成藏条件定量评价方法 | |
CN110703354B (zh) | 一种基于井震结合的岩性边界识别方法 | |
CN111506861B (zh) | 一种目的层有利区域裂缝强度计算方法 | |
CN111399048A (zh) | 一种对断溶体计算相关属性及数据加权重构的方法 | |
Yuan et al. | 6D phase-difference attributes for wide-azimuth seismic data interpretation | |
CN112114380A (zh) | 一种用于砂砾岩油藏储层的预测方法 | |
CN112946743B (zh) | 区分储层类型的方法 | |
CN112505754B (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
Nardone et al. | Shallow velocity structure and site effects at Mt. Vesuvius, Italy, from HVSR and array measurements of ambient vibrations | |
CN115857047B (zh) | 一种地震储层综合预测方法 | |
CN111399047A (zh) | 一种基于道集数据重构的增强河道地质异常体成像的方法 | |
CN110703329A (zh) | 基于弱振幅地震反射形成机制的岩性油藏边界确定方法 | |
Sanda et al. | The integrated approach to seismic attributes of lithological characterization of reservoirs: case of the F3 Block, North Sea-Dutch Sector | |
Deng et al. | A new index used to characterize the near-wellbore fracture network in naturally fractured gas reservoirs | |
CN113671566B (zh) | 一种基于深度域地震数据计算裂缝参数的方法 | |
CN110795513B (zh) | 河流相源储异位型致密油气甜点区分布的预测方法 | |
CN110297264B (zh) | 一种低渗气藏薄储层甜点地震预测方法 | |
CN113050168B (zh) | 基于阵列声波测井和声波远探测测井资料的裂缝有效性评价方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |