CN104898164B - 一种基于地震相位微变分析的致密薄储层地震预测方法 - Google Patents
一种基于地震相位微变分析的致密薄储层地震预测方法 Download PDFInfo
- Publication number
- CN104898164B CN104898164B CN201510126506.1A CN201510126506A CN104898164B CN 104898164 B CN104898164 B CN 104898164B CN 201510126506 A CN201510126506 A CN 201510126506A CN 104898164 B CN104898164 B CN 104898164B
- Authority
- CN
- China
- Prior art keywords
- phase
- seismic
- change
- analysis
- earthquake
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000004576 sand Substances 0.000 claims description 69
- 238000003384 imaging method Methods 0.000 claims description 5
- 208000035126 Facies Diseases 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 abstract description 10
- 239000003208 petroleum Substances 0.000 abstract description 3
- 230000003595 spectral effect Effects 0.000 description 4
- 239000004575 stone Substances 0.000 description 4
- 239000002131 composite material Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 239000003245 coal Substances 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及石油地球物理勘探领域,特别是基于地震相位微变分析的致密薄储层地震预测方法,其特征是:所述方法包括4步骤实现基于岩性变化分析的地震层位解释追踪技术,本方法与普通层位解释的区别在于利用地震相位变化预测薄储层,判断是否是由区域构造引起,可以追踪同相轴的任意位置,依此来突出地震相位的微小变化,进而对致密薄储层进行预测。
Description
技术领域
本发明涉及石油地球物理勘探领域,特别是基于地震相位微变分析的致密薄储层地震预测方法。
背景技术
一般把厚度小于1/4波长的储层称为薄储层,目前薄储层地震预测方法有两类,一类是只用地震信息,如有分频成像(谱分解)、地层切片等,另一类就是井震联合反演,如井约束波阻抗反演、地质统计学反演等。
分频技术主要通过谱分解,将地震资料中的高频成份分离出来,通过不同频率的扫描,观察振幅能量变化,再依据不同频率对应的不同调谐厚度来识别薄储层(孙学栋,赵建儒,白军等.谱分解技术在营尔凹陷长沙岭地区薄储层预测中的应用.石油地球物理勘探,2011,46(增刊):72-75;徐胜峰,刘春园、季玉新等.分频技术在塔河油田石炭系薄储层预测中的应用.石油天然气学报,2011,33(2):65-69;徐丽英,徐鸣杰,陈振岩.利用谱分解技术进行薄储层预测.石油地球物理勘探,2006,41(3):299-303;邓和平.薄层预测技术概述.内蒙古石油化工,2009,17:82-83)。国内各主要盆地目的层的地震波速度为3500-5000m/s之间,地震波的频宽10-75Hz,因此该方法理论上能识别的最薄储层为11.7m,而在实际应用中,由于高频资料信噪比低、能量弱,不可能用到75Hz频率数据预测储层,因此,该方法适合于 15m以上薄储层预测;地层切片它以解释的两个等时界面为顶底,在地层的顶底界面间按照厚度等比例内插出一系列的层面,沿这些内插出的层面逐一生成切片,这些切片的属性就代表了同一沉积时期地质体的特征(王学习,毕建军,王利功等.地层切片穿时现象对地震属性的影响.物探与化探,2012,36(1):94-98)。一般认为地层切片技术可以突破地震1/4波长分辨率的极限,达到1/8波长(张军华,周振晓,谭明友等.地震切片解释中的几个理论问题.石油地球物理勘探,2007,42(3)348-352),即使这样,该方法分辨储层的能力也与分频成像相近。
第二类薄储层预测的方法就是井震联合反演,该方法借助于井上丰富的低频及高频信息,提高了地震的分辨率,它分辨薄储层的能力可以达到测井水平,但缺点也是显而易见的,模型占的比例大了,受井的影响较大,基本等同于井间内插,模型占的比例小了又达不到分辨薄层的目的(刘占族,张雷,霍丽娜等.地质统计学反演在煤层气薄储层识别中的应用.石油地球物理勘探,2012,47(1):30-34;张明振.对测井约束地震波阻抗反演的理解与应用.油气地球物理,2006,4(3):1-5.)。
鄂尔多斯盆地三叠系延长组是一套以湖相沉积为主的地层,自下而上分为十个油层组,分别为长10-长1油层组。岩性油藏主要分布在长4+5-长8油层组,为三角洲前缘沉积,沉积微相为水下分流河道,砂体薄、致密,为典型的致密储层。鄂尔多斯盆地油区地震资料主频30-40Hz,频宽为10-70Hz,大套油层组在地震剖面上为一个地 震反射同相轴,地震可识别,油层组中一般发育1-5套砂体,单砂体厚度5-20m,利用目前的薄储层预测技术难以预测油层组中的单砂体厚度。
发明内容
本发明的目的是提供一种基于地震相位微变分析的致密薄储层地震预测方法,解决储层厚度小于1/8波长致密薄储层的预测难题,为致密薄层勘探、油田开发单砂体预测、水平井轨迹设计及油藏精细描述提供依据,提高致密薄储层勘探开发效率。
本发明解决其技术问题采用的技术方案是:一种基于地震相位微变分析的致密薄储层地震预测方法,其特征是:所述方法包括以下步骤:
步骤1:正演分析,确定同一反射轴内砂体厚度变化引起地震反射轴相位变化幅度及同相轴在纵向上的时差变化大小,选择目的层不同砂体层厚度的井进行正演分析;
步骤2:基于岩性变化分析的层位解释追踪,对目的层位解释追踪;
利用地震相位变化预测薄储层,在地震相位稳定时,层位解释仍追踪波峰或波谷,当相位发生变化时,首先判断是否是由区域构造引起,如果相邻地震反射轴发生同样的变化,则是由区域构造引起,层位解释仍然与之前一致,如果相邻地震反射轴无变化,则是由砂体变化引起,层位解释不能再追踪波峰或是波谷,要按照原构造趋势追踪同相轴的任意位置;
步骤3:时窗选择,依据步骤2的层位解释结果,根据正演分析得到了不同砂体厚度变化同相轴的时差确定时窗;
步骤4:属性分析,在步骤3时窗内,利用地震均方根振幅、分频成像、地震相分析这些属性进行分析,并利用实际的井资料进行对比标定,来预测致密薄储层的砂体厚度。
所述的属性分析的时窗选择依据为做步骤1正演分析时将目的层放到同一时间位置,通过正演,比较不同砂体厚度地震反射同相轴相位的变化,找出同相轴在纵向上时差大小。
所述的正演分析选择目的层砂体较发育井,通过改变砂体的阻抗信息,使它变为目的层砂体厚度不同多口井,进行对比分析。
所述的时窗选择中,如果砂体厚度的变化引起相位增加,反射同相轴下移,沿层的时窗不能包含向下移动的同相轴的最大振幅,反之,如果砂体厚度的变化引起相位减小,反射同相轴上移,那边沿层的时窗不能包含向上移动的同相轴的最大振幅。
本发明的有益效果是:本方法基于岩性变化分析的地震层位解释追踪技术,以往对层位的解释是在井的标定下,多数是追踪波峰或是波谷,特别是某一套地层对应于一个反射同相轴的情况下,更是严格追踪波峰或是波谷,本方法与普通层位解释的区别在于利用地震相位变化预测薄储层,判断是否是由区域构造引起,可以追踪同相轴的任意位置,依此来突出地震相位的微小变化,进而对致密薄储层进行预测。
附图说明
下面结合实施例附图对本发明进一步说明。
图1是本发明实施例的流程图;
图2是Y1井长81段测井解释结果示意图;
图3是利用25Hz雷克子波制作的合成地震记录;
图4a是把长81段上部砂岩的阻抗变为与泥岩相当后的合成地震记录;
图4b是把长81下段部分砂岩的阻抗变为与泥岩相当后的合成地震记录;
图5a是x1井长31段测井解释成果图;
图5b是x2井长31段测井解释成果图;
图6是过x1、x2井地震剖面图。
具体实施方式
地震反射同相轴的相位变化多用来进行断层描述,当地层中存在断层时,地震反射同相轴的相位错断明显,依此来进行断层识别,但利用地震相位变化进行薄储层预测还未见报道。
地震的一个反射同相轴是同相轴所对应的一套地层内多个地震反射的叠合响应,当地层内砂泥岩结构稳定时,地震反射同相轴的相位比较稳定,当砂泥结构发生变化时,地震反射轴的相位就会发生微小变化,通过这一微小的变化,结合地质特点,就可以对单一砂体进行预测。
如图1的本发明基于地震相位微变分析的致密薄储层地震预测方法流程图所示,本发明通过以下步骤实现:
步骤1:正演分析,确定同一反射轴内砂体厚度变化引起地震反射轴相位变化幅度及同相轴在纵向上的时差变化大小,选择目的层不同砂体层厚度的井进行正演分析。
做正演时将目的层放到同一时间位置,通过正演,比较不同砂体厚度地震反射同相轴相位的变化,找出同相轴在纵向上时差大小,为后续属性分析的时窗选择提供依据;
在做正演分析时,也可以选择目的层砂体较发育井,通过改变砂体的阻抗信息,使它变为目的层砂体厚度不同多口井,这样更有利于对比分析。
通过地震正演并结合实际地震资料情况,建立不同厚度砂体变化与地震相位变化之间的关系。正演分析最方便快捷的方法是在strata中作合成记录,在正演过程中选择richer子波即可,子波的主频要依据实际工区地震资料的主频及频宽情况来确定。
步骤2:基于岩性变化分析的层位解释追踪,对目的层位解释追踪。利用地震相位变化预测薄储层,在地震相位稳定时,层位解释仍追踪波峰或波谷,当相位发生变化时,首先判断是否是由区域构造引起,如果相邻地震反射轴发生同样的变化,则是由区域构造引起,层位解释仍然与之前一致,如果相邻地震反射轴无变化,则是由砂体变化引起,层位解释不能再追踪波峰或是波谷,要按照原构造趋势追踪同相轴的意义位置。
本方法实现的关键就是基于岩性变化分析的地震层位解释追踪技术,以往对层位的解释是在井的标定下,多数是追踪波峰或是波谷, 特别是某一套地层对应于一个反射同相轴的情况下,更是严格追踪波峰或是波谷,本方法与普通层位解释的区别在于利用地震相位变化预测薄储层,判断是否是由区域构造引起,可以追踪同相轴的任意位置,依此来突出地震相位的微小变化,进而对致密薄储层进行预测。
步骤3:时窗选择,依据步骤2的层位解释结果,根据正演分析得到了不同砂体厚度变化同相轴的时差确定时窗。
此外在时窗选区时,要体现同相轴相位变化,如果砂体厚度的变化引起相位增加,也就是反射同相轴下移,沿层的时窗不能包含向下移动的同相轴的最大振幅,反之,如果砂体厚度的变化引起相位减小,也就是反射同相轴上移,那边沿层的时窗不能包含向上移动的同相轴的最大振幅。
步骤4:属性分析,在步骤3时窗内,利用地震均方根振幅、分频成像、地震相分析这些属性进行分析,并利用实际的井资料进行对比标定,来预测致密薄储层的砂体厚度。
该方法的薄储层分辨能力与地震资料品质、同一反射同相轴内砂体结构、砂泥岩的阻抗差异有关,当地震资料保真、保幅处理的较好、砂体结构相对简单(同一同相轴内发育两、三套砂体)、砂泥岩波阻抗差异明显的地区应用效果更好。
以鄂尔多斯盆地为例,鄂尔多斯盆地内部构造平缓,地层倾角约0.5°~0.7°,且同一套油层组厚度变化不大,因此在主砂带上,地震反射同相轴相位比较稳定。一套油层组发育1-5套小层砂体,开发的主力小层砂体一般占总砂厚的50%以上,主力小层砂体变薄后,会 引起地震相位角的增加或是减少。
如图2所示,Y1井长81段测井解释结果,其中Y1井长81段砂厚28.5m,发育两套砂体;由图3和图4a的Y1井长81段上部砂变薄后地震反射同相轴相位变化对比图可知,图3是利用25Hz雷克子波制作的合成地震记录,图4a是把长81段上部砂岩的阻抗变为与泥岩相当后的合成地震记录,即上部砂体变薄厚度后的合成记录,所用的子波与图3的一致,从图中可以看出,长81段的地震反射振幅强度基本没变(当砂体厚度大于15m,地震反射同相轴的振幅差异不大),而图4a的长81段反射同相轴的相位增加,同相轴明显向下移动。
由图3和图4b的Y1井长81段上部砂体变薄后地震反射同相轴相位变化对比图可知,图3是利用25Hz雷克子波制作的合成地震记录,图4b是把长81下段部分砂岩的阻抗变为与泥岩相当后的合成地震记录,即下段砂体变薄后的合成记录,所用的子波与图3的一致,从图中可以看出,长81段的地震反射振幅强度基本没变(当砂体厚度大于15m,地震反射同相轴的振幅差异不大),而图4b的长81段反射同相轴的相位减小,同相轴明显向上移动。
利用地震反射同相轴相位的变化,可以识别鄂尔多斯盆地5m以上的致密薄储层,远远突破的地震1/4或是1/8波长的限制。经过鄂尔多斯盆地现场试验论证5m以上单砂体的地震预测符合率达78.8%具体见图5、图6,特别是在该区水平井位优选中效果明显,水平井有效储层钻遇率达89.5%,较以往提高了4.3%。
其中图5a是x1井长31段测井解释成果图,图5b是x2井长31 段测井解释成果图,其中x1井长31 1砂厚9.5m,x2井长31 1砂厚3.9m,再根据图6的x1、x2井地震剖面图可知,长31段为一个地震反射同相轴,x1井处长31 1砂厚为9.5m,在图中地震相位变化点,长31的地震相位增加了90°,反射同相轴向下移动,而相邻的反射同相轴变化不大,因此长31 1相位的改变是由长31 1砂体厚度变化引起,变化点之后长31层位追踪不能沿着波峰,而是按照原来的趋势追踪零值点,这样沿长31层的地震属性就可以反映长31 1砂体的变化。根据地震解释结果,预测x1井长31 1砂厚为10m,实钻为9.5m;预测x2井长31 1砂体厚度为3.5m,实钻为3.9m,地震预测与实钻吻合较好。
本发明利用地震相位变化进行薄储层预测时,首先要根据不同地区具体地质特点,通过地震正演方法,确定砂体厚度变化量与地震相位的变化关系;其次,利用基于岩性变化分析技术对地震资料进行层位解释追踪,以往对目的层的解释多数是追踪波峰或是波谷,该方法与普通层位解释的区别在于层位追踪同相轴的任意位置,在地震相位稳定时,层位解释仍追踪波峰或波谷,当相位发生变化时,首先判断是否是由区域构造引起,如果相邻地震反射轴发生同样的变化,则是由区域构造引起,层位解释仍然与之前一致,如果相邻地震反射轴无变化,则是由砂体变化引起,层位解释不能跟着地震相位走,要按照原构造趋势走;再次进行沿层属性提取。在此过程中要注意时窗的选取,不能过大,也不能过小,要以突出目的层相位的细微变化为根本目标。这样在属性分析图上就可以识别薄储层的变化。
Claims (4)
1.一种基于地震相位微变分析的致密薄储层地震预测方法,其特征是:所述方法包括以下步骤:
步骤1:正演分析,确定同一反射轴内砂体厚度变化引起地震反射轴相位变化幅度及同相轴在纵向上的时差变化大小,选择目的层不同砂体层厚度的井进行正演分析;
步骤2:基于岩性变化分析的层位解释追踪,对目的层位解释追踪;
利用地震相位变化预测薄储层,在地震相位稳定时,层位解释仍追踪波峰或波谷,当相位发生变化时,首先判断是否是由区域构造引起,如果相邻地震反射轴发生同样的变化,则是由区域构造引起,层位解释仍然与之前一致,如果相邻地震反射轴无变化,则是由砂体变化引起,层位解释不能再追踪波峰或是波谷,要按照原构造趋势追踪同相轴的任意位置;
步骤3:时窗选择,依据步骤2的层位解释结果,根据正演分析得到了不同砂体厚度变化同相轴的时差确定时窗;
步骤4:属性分析,在步骤3时窗内,利用地震均方根振幅、分频成像、地震相分析这些属性进行分析,并利用实际的井资料进行对比标定,来预测致密薄储层的砂体厚度;
所述的步骤2中,利用地震相位变化预测薄储层,判断是否是由区域构造引起,追踪同相轴的任意位置,依此来突出地震相位的微小变化,进而对致密薄储层进行预测。
2.根据权利要求1所述的一种基于地震相位微变分析的致密薄储层地震预测方法,其特征是:所述的属性分析的时窗选择依据为做步骤1正演分析时将目的层放到同一时间位置,通过正演,比较不同砂体厚度地震反射轴相位的变化,找出同相轴在纵向上时差变化大小。
3.根据权利要求1所述的一种基于地震相位微变分析的致密薄储层地震预测方法,其特征是:所述的正演分析选择目的层砂体较发育井,通过改变砂体的阻抗信息,使它变为目的层砂体厚度不同多口井,进行对比分析。
4.根据权利要求1所述的一种基于地震相位微变分析的致密薄储层地震预测方法,其特征是:所述的时窗选择中,如果砂体厚度的变化引起相位增加,反射同相轴下移,沿层的时窗不能包含向下移动的同相轴的最大振幅,反之,如果砂体厚度的变化引起相位减小,反射同相轴上移,沿层的时窗不能包含向上移动的同相轴的最大振幅。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510126506.1A CN104898164B (zh) | 2015-03-23 | 2015-03-23 | 一种基于地震相位微变分析的致密薄储层地震预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510126506.1A CN104898164B (zh) | 2015-03-23 | 2015-03-23 | 一种基于地震相位微变分析的致密薄储层地震预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104898164A CN104898164A (zh) | 2015-09-09 |
CN104898164B true CN104898164B (zh) | 2017-10-17 |
Family
ID=54030924
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510126506.1A Active CN104898164B (zh) | 2015-03-23 | 2015-03-23 | 一种基于地震相位微变分析的致密薄储层地震预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104898164B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111158048B (zh) * | 2020-01-04 | 2022-10-04 | 杨林海 | 一种地震波形包络解释提高储层预测精度的分析方法 |
CN112099085B (zh) * | 2020-09-10 | 2022-01-18 | 成都理工大学 | 一种基于转换波的碳酸盐岩微古地貌雕刻及储层预测方法 |
CN114428322A (zh) * | 2020-10-29 | 2022-05-03 | 中国石油天然气股份有限公司 | 基于频率属性预测薄储层厚度的方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103777243A (zh) * | 2012-10-25 | 2014-05-07 | 中国石油化工股份有限公司 | 砂泥岩薄互层储层厚度预测方法 |
CN104280773B (zh) * | 2013-07-12 | 2017-04-05 | 中国石油天然气集团公司 | 利用随炮检距变化的时频谱交汇图预测薄层厚度的方法 |
-
2015
- 2015-03-23 CN CN201510126506.1A patent/CN104898164B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104898164A (zh) | 2015-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA3076280C (en) | Method and apparatus for determining thin interlayers | |
CN104502969A (zh) | 河道砂岩性油藏的识别方法 | |
Huang et al. | Use of nonlinear chaos inversion in predicting deep thin lithologic hydrocarbon reservoirs: A case study from the Tazhong oil field of the Tarim Basin, China | |
CN105527653B (zh) | 一种基于地质信息的虚拟井构建方法 | |
CN103527184A (zh) | 一种白云岩储层的预测方法和系统 | |
WO2017035104A1 (en) | Velocity model seismic static correction | |
Nascimento et al. | High-resolution acoustic impedance inversion to characterize turbidites at Marlim Field, Campos Basin, Brazil | |
CN104898164B (zh) | 一种基于地震相位微变分析的致密薄储层地震预测方法 | |
Nejadi et al. | Look ahead of the bit while drilling: Potential impacts and challenges of acoustic seismic while drilling in the mcmurray formation | |
Rodríguez-Pradilla | Microseismic monitoring of a hydraulic-fracturing operation in a CBM reservoir: Case study in the Cerrejón Formation, Cesar-Ranchería Basin, Colombia | |
Liu et al. | Seismic characterization of fault and fractures in deep buried carbonate reservoirs using CNN-LSTM based deep neural networks | |
Changzi et al. | Seismic prediction of sweet spots in the Da'anzhai shale play, Yuanba area, the Sichuan Basin | |
Aminzadeh et al. | Fundamentals of Petroleum Geophysics | |
Jafari et al. | Estimation of acoustic-impedance model by using model-based seismic inversion on the Ghar Member of Asmari Formation in an oil field in southwestern Iran | |
Cheng et al. | Reservoir prediction in a development area with a high-density well pattern using seismic sedimentology: An example from the BB2 block, Changyuan LMD oil field, Songliao Basin, China | |
CN104360386A (zh) | 一种针对花岗岩地层划分对比的面元法 | |
RU2145101C1 (ru) | Способ оценки эксплуатационных свойств нефтегазовой залежи | |
CN105989240A (zh) | 砂层预测方法及装置 | |
Torres et al. | Quantitative 3D seismic characterization of potential lacustrine-stratigraphic intervals in the Woodford Shale Oklahoma, using multiattribute seismic displays and applying a supervised neural network | |
Imam et al. | Extracting high‐resolution P‐wave reflectivity of the shallow subsurface by seismic interferometry based on autocorrelation of blast mining signals | |
Nakajima et al. | Characterization on reservoir complex and CO2 plume with Vp/Vs: Case study at Nagaoka site, Japan | |
Wood et al. | Seismic attribute and sequence stratigraphic integration methods for resolving reservoir geometry in San Jorge Basin, Argentina | |
Naseer | Quantitative appraisal of tectonically-influenced hydrocarbon-bearing Late-Cretaceous fluvial depositional system, Southwest Pakistan using spectral waveform-based instantaneous lateral thickness variability static simulations | |
House et al. | Understanding hydraulic fractures in tight-gas sands through the integration of borehole microseismic data, three-dimensional surface seismic data, and three-dimensional vertical seismic profile data: A Jonah field case study | |
Leiceaga et al. | A seismic reservoir characterization work flow for reducing risk utilizing AVO, simultaneous prestack inversion, and rock physics |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |