CN105093303A - 一种河流相砂体不连续界线精细表征的地震解释方法 - Google Patents
一种河流相砂体不连续界线精细表征的地震解释方法 Download PDFInfo
- Publication number
- CN105093303A CN105093303A CN201510490652.2A CN201510490652A CN105093303A CN 105093303 A CN105093303 A CN 105093303A CN 201510490652 A CN201510490652 A CN 201510490652A CN 105093303 A CN105093303 A CN 105093303A
- Authority
- CN
- China
- Prior art keywords
- sand body
- attribute
- plane
- discontinuous boundary
- boundary line
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
Abstract
本发明涉及一种河流相砂体不连续界线精细表征的地震解释方法,包括以下步骤:1)根据河流相砂体正韵律沉积特征,结合后续属性分析对属性提取时窗的需求,采用砂体包络法对砂体单元的构造进行解释;2)对目标储层段提取蚂蚁追踪属性和曲率属性;3)以解释的砂体包络顶底为时窗提取振幅平面属性、曲率平面属性和蚂蚁追踪平面属性;4)取地震剖面的曲率平面属性和蚂蚁追踪平面属性,进行对比分析,从而判定砂体不连续界线;5)振幅类属性纵向切片演绎,辅助判别砂体纵向叠置导致的不连续界线以及划分砂体沉积期次;6)利用步骤4)获得的不连续界线平面的判断结果和步骤5)获得的纵向砂体期次划分结论,共同辅助沉积相认识的建立。
Description
技术领域
本发明涉及一种河流相砂体不连续界线精细表征的地震解释方法,属于油田开发技术领域。
背景技术
河流相储层是我国已发现油田的主要油气资源储层类型,也是油气最为富集的储油层。据统计,河流相储层原油地质储量占我国已开发油田动用地质储量的40%以上,但这些油田经过几十年注水开发,水驱平均采收率只有30%,地质储量的采出程度则更低,为23%左右,这意味着储层内有四分之三的地质储量无法有效采出。这种现象除了开发因素导致,河流相储层内部的非均质性是主要原因。
河流相储层岩性的纵向和横向变化大,连续性差。纵向上,不同期次的河道沉积砂体相互叠置或切削;平面上,河道迁移摆动,形成复杂的交错和接触关系。砂体间的不连续界线成为隔断油水关系,阻碍水驱油的主要屏障。河流相砂体不连续界线的精细识别和表征有助于正确认识油藏的成藏规律,有效提高探井成功率及开发方案设计。
在很多地区,河流相砂体属于典型的薄互层砂泥岩沉积,在地震上通常表现为一定层段内薄互层的复合响应,加上噪音干扰,很难通过地震剖面或常规反演剖面预测砂体的厚度变化及砂体分布的界线所在。此外,短轴状或点状的地震反射特征也很难在地震剖面上实现连续追踪。针对河流相砂体研究的地震解释方法通常以开展高分辨率地震反演和属性分析为主,但目前为止尚未形成一套专门针对河流相砂体不连续界线识别的地震解释方法。
发明内容
针对上述问题,本发明的目的是提供一种能够实现砂体展布横向边界的刻画和纵向沉积期次的划分的河流相砂体不连续界线精细表征的地震解释方法。
为实现上述目的,本发明采用以下技术方案:一种河流相砂体不连续界线精细表征的地震解释方法,包括以下步骤:1)根据河流相砂体正韵律沉积特征,结合后续属性分析对属性提取时窗的需求,采用砂体包络法对砂体单元的构造进行解释;2)对目标储层段提取蚂蚁追踪属性和曲率属性;3)以解释的砂体包络顶底为时窗提取振幅平面属性、曲率平面属性和蚂蚁追踪平面属性;4)取地震剖面的曲率平面属性和蚂蚁追踪平面属性,进行对比分析,从而判定砂体不连续界线;5)振幅类属性纵向切片演绎,辅助判别砂体纵向叠置导致的不连续界线以及划分砂体沉积期次;6)利用步骤4)获得的不连续界线平面的判断结果和步骤5)获得的纵向砂体期次划分结论,共同辅助沉积相认识的建立。
所述步骤4)中,判定砂体不连续界线的过程如下:①以地震剖面-属性平面结合方式对提取的曲率属性和蚂蚁追踪属性进行对比,从中选取出能够表征大部分不连续界线的属性作为主控属性;②将主控属性以外的另一属性根据识别砂体不连续界线识别需求设置提取阈值,将大于或小于该阈值的信息提取并平面显示,作为补充属性;③将主控属性和补充属性进行平面叠合显示,得到最终表征砂体不连续界线的融合属性,融合属性能够全面反映地震响应同相轴的横向不连续特征;④以地震剖面-属性平面结合方式对融合属性进行精细解剖,剔除断层、岩性突变等因素对砂体不连续界线判断造成的干扰,筛除干扰之后,融合属性平面图上的信息作为判别砂体不连续界线的最终依据。
所述步骤5)中,振幅类属性纵向切片演绎的过程如下:①以砂体包络法解释的砂体单元顶底为界,利用层间等比例差值法对砂体单元进行劈分,劈分厚度不小于一个地震采样点;②对劈分后的每个层序单元提取振幅类平面属性,即切片属性,用以表征储层的纵向变化特征;③将地震属性切片特征相似的劈分单元合并后再次提取步骤②所提取的振幅类平面属性,得到的平面属性特征认定为同一时间段的储层沉积特征;④将不同历史时期的平面属性按照沉积时间先后进行对比分析,用于辅助砂体纵向沉积期次划分和砂体纵向叠置导致的不连续界线判别。
本发明由于采取以上技术方案,其具有以下优点:1、本发明技术针对河流相砂体的地震响应特征提出“砂体包络解释法”,相比现有构造解释方法更加针对砂体不连续界线的表征需求,具有更高的解释精度。2、本发明创新提出了“主控属性”与“补充属性”进行平面融合的方法来进行平面砂体不连续界线表征,相比常规技术具有更精确的识别能力。
附图说明
图1是本发明方法的流程示意图;
图2是本发明针对河流相正韵律沉积的地震反射特征分析图;
图3是本发明采用“砂体包络法”解释的目标砂体单元的顶底;
图4是本发明对目标储层段所提取的针对某个地震剖面的蚂蚁追踪属性和曲率属性对比图;其中,图(a)是目标储层段的地震剖面;图(b)是目标储层段的蚂蚁追踪属性;图(c)是目标储层段的曲率属性;
图5是本发明对目标储层段所提取的针对另一个地震剖面的蚂蚁追踪属性和曲率属性对比图蚂蚁追踪属性曲率属性;其中,图(a)是目标储层段的地震剖面;图(b)是目标储层段的蚂蚁追踪属性;图(c)是目标储层段的曲率属性;蚂蚁追踪属性曲率属性蚂蚁追踪属性
图6是本发明以解释的砂体包络顶底为时窗提取的蚂蚁追踪属性和曲率属性的平面对比图;其中,图(a)是提取的曲率平面属性;图(b)是提取的蚂蚁追踪平面属性;
图7是本发明蚂蚁追踪平面属性提取阈值后的平面显示;
图8是本发明将“主控属性”和“补充属性”进行平面叠合显示的示意图;
图9是本发明以地震剖面-属性平面结合方式对“融合属性”进行精细解剖,剔除断层、岩性突变等因素对砂体不连续界线判断造成的干扰,筛除干扰之后,融合属性平面图上的信息可以作为判别砂体不连续界线的最终依据;
图10是本发明利用层间等比例差值法对砂体单元进行劈分得到的小层单元;
图11是本发明对劈分后的每个层序单元提取振幅类敏感属性切片;其中,图(a)~图(g)分别是对第一~第七层序所提取的振幅类敏感属性切片(1)~(7);
图12是本发明将地震属性切片特征相似的劈分单元合并后再次提取属性所得到的属性切片;其中,图(a)是切片(6)(7)合并后再次提取的属性切片;图(b)是切片(3)(4)(5)合并后再次提取的属性切片;图(c)是切片(1)(2)合并后再次提取的属性切片。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提出了一种河流相砂体不连续界线精细表征的地震解释方法,包括以下步骤:
1)根据河流相砂体正韵律沉积特征,结合后续属性分析对属性提取时窗的需求,采用“砂体包络法”对砂体单元的构造进行解释。其中,后续属性具体包括步骤2)中所涉及到的蚂蚁追踪属性和曲率属性。
如图2所示,通过对目标储层地震反射特征连井对比发现,对于河流相正韵律沉积,砂体单元的泥顶和砂底的地震响应特征明显,而砂体单元顶由于处于沉积渐变位置,反而不易解释。如图3所示,对比还发现,砂体沉积的不连续性主要表现为地震波形的横向变化。因此,为满足后续属性分析需求,采用“砂体包络法”解释砂体单元的泥顶和砂底界面,实现砂体单元的构造解释。
2)对目标储层段提取蚂蚁追踪属性和曲率属性。
曲率是曲线的一种二维特征,用于描述曲线上某一点的弯曲程度,在地震属性分析研究方面,曲率属性主要用于断层、裂缝及河道等的刻画和表征。蚂蚁追踪属性是一种地震属性,主要用于断裂系统的自动识别和解释,通过调整追踪参数,也可有效识别砂体不连续接触导致的地震同相轴不连续特征。
图4给出了对目标储层段所提取的针对某一地震剖面的蚂蚁追踪属性和曲率属性对比结果,从图中可以看出,在该地震剖面上,蚂蚁追踪属性和曲率属性在地震剖面上表现出的波形不连续及能量强弱的表征效果上具有较好的一致性,并且,曲率属性对于波形的振幅、形态等各种细节变化的表征能力更强,尤其侧重同相轴扭曲的表征。
图5给出了对目标储层段所提取的针对另一个地震剖面的蚂蚁追踪属性和曲率属性对比结果,,从图中可以看出,在该地震剖面上,蚂蚁追踪属性对于表现为同相轴错断的波形不连续表征能力更强。
3)以解释的砂体包络顶底为时窗提取振幅平面属性、曲率平面属性和蚂蚁追踪平面属性。
图6给出了以解释的砂体包络顶底为时窗提取的蚂蚁追踪属性和曲率属性的平面图对比情况,其中,图(a)是提取的曲率平面属性;图(b)是提取的蚂蚁追踪平面属性,从图中可以看出,两种属性表征的信息有所重叠,曲率平面属性的信息更为丰富。
4)针对砂体不连续界线开展综合属性分析研究,其具体步骤为:
①以地震剖面-属性平面结合方式对提取的曲率属性和蚂蚁追踪属性进行对比(结合图4~6),从中选取出能够表征大部分不连续界线的属性作为主控属性;在本实施例中,从图4~6对比结果,选取曲率属性作为主控属性。
②将主控属性以外的另一属性根据识别砂体不连续界线识别需求设置提取阈值,将大于或小于该阈值的信息提取并平面显示,作为补充属性(如图7所示);
③将主控属性和补充属性进行平面叠合显示,得到最终表征砂体不连续界线的融合属性(如图8所示),融合属性可以全面反映地震响应同相轴的横向不连续特征;
④以地震剖面-属性平面结合方式对融合属性进行精细解剖,剔除断层、岩性突变等因素对砂体不连续界线判断造成的干扰(如图9所示),筛除干扰之后,融合属性平面图上的信息可以作为判别砂体不连续界线的最终依据。
5)振幅类属性纵向切片演绎,辅助砂体纵向叠置导致的不连续界线判别和砂体沉积期次划分,具体包括以下步骤:
①以“砂体包络法”解释的砂体单元顶底为界,利用层间等比例差值法对砂体单元进行劈分(如图10所示),劈分厚度不小于一个地震采样点;
②对劈分后的每个层序单元提取振幅类平面属性,用以表征储层的纵向变化特征,(如图11所示);振幅类地震属性是指根据地震振幅信息得到的属性,例如最大振幅、平均振幅、均方根振幅等。
③将地震属性切片特征相似的劈分单元合并后再次提取步骤②所提取的振幅类平面属性,得到的平面属性特征认定为同一时间段的储层沉积特征(如图12所示),从反映的储层沉积特征来看,可以将切片(1)(2)合并为一期砂体沉积,如图中(c)所示,切片(3)(4)(5)合并为一期,如图中(b)所示,切片(6)(7)合并为一期,如图中(a)所示;
④将不同历史时期的平面属性按照沉积时间先后进行对比分析,用于辅助砂体纵向沉积期次划分和砂体纵向叠置导致的不连续界线判别。
6)利用步骤4)获得的不连续界线平面的判别结果和步骤5)获得的纵向砂体期次划分结论,共同辅助沉积相认识的建立。
上述各实施例仅用于对本发明的目的、技术方案和有益效果进行了进一步详细说明,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种河流相砂体不连续界线精细表征的地震解释方法,包括以下步骤:
1)根据河流相砂体正韵律沉积特征,结合后续属性分析对属性提取时窗的需求,采用砂体包络法对砂体单元的构造进行解释;
2)对目标储层段提取蚂蚁追踪属性和曲率属性;
3)以解释的砂体包络顶底为时窗提取振幅平面属性、曲率平面属性和蚂蚁追踪平面属性;
4)取地震剖面的曲率平面属性和蚂蚁追踪平面属性,进行对比分析,从而判定砂体不连续界线;
5)振幅类属性纵向切片演绎,辅助判别砂体纵向叠置导致的不连续界线以及划分砂体沉积期次;
6)利用步骤4)获得的不连续界线平面的判断结果和步骤5)获得的纵向砂体期次划分结论,共同辅助沉积相认识的建立。
2.如权利要求1所述的一种河流相砂体不连续界线精细表征的地震解释方法,其特征在于:所述步骤4)中,判定砂体不连续界线的过程如下:
①以地震剖面-属性平面结合方式对提取的曲率属性和蚂蚁追踪属性进行对比,从中选取出能够表征大部分不连续界线的属性作为主控属性;
②将主控属性以外的另一属性根据识别砂体不连续界线识别需求设置提取阈值,将大于或小于该阈值的信息提取并平面显示,作为补充属性;
③将主控属性和补充属性进行平面叠合显示,得到最终表征砂体不连续界线的融合属性,融合属性能够全面反映地震响应同相轴的横向不连续特征;
④以地震剖面-属性平面结合方式对融合属性进行精细解剖,剔除断层、岩性突变等因素对砂体不连续界线判断造成的干扰,筛除干扰之后,融合属性平面图上的信息作为判别砂体不连续界线的最终依据。
3.如权利要求1所述的一种河流相砂体不连续界线精细表征的地震解释方法,其特征在于:所述步骤5)中,振幅类属性纵向切片演绎的过程如下:
①以砂体包络法解释的砂体单元顶底为界,利用层间等比例差值法对砂体单元进行劈分,劈分厚度不小于一个地震采样点;
②对劈分后的每个层序单元提取振幅类平面属性,即切片属性,用以表征储层的纵向变化特征;
③将地震属性切片特征相似的劈分单元合并后再次提取步骤②所提取的振幅类平面属性,得到的平面属性特征认定为同一时间段的储层沉积特征;
④将不同历史时期的平面属性按照沉积时间先后进行对比分析,用于辅助砂体纵向沉积期次划分和砂体纵向叠置导致的不连续界线判别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510490652.2A CN105093303B (zh) | 2015-08-11 | 2015-08-11 | 一种河流相砂体不连续界线精细表征的地震解释方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510490652.2A CN105093303B (zh) | 2015-08-11 | 2015-08-11 | 一种河流相砂体不连续界线精细表征的地震解释方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105093303A true CN105093303A (zh) | 2015-11-25 |
CN105093303B CN105093303B (zh) | 2017-10-03 |
Family
ID=54574157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510490652.2A Active CN105093303B (zh) | 2015-08-11 | 2015-08-11 | 一种河流相砂体不连续界线精细表征的地震解释方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105093303B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107247289A (zh) * | 2017-06-07 | 2017-10-13 | 中海石油(中国)有限公司 | 一种河流相地层的中期沉积旋回末期等时界面恢复方法 |
CN107966731A (zh) * | 2017-11-08 | 2018-04-27 | 西南石油大学 | 一种基于地震波形结构属性的河流相砂体叠置区识别方法 |
CN110441818A (zh) * | 2019-07-31 | 2019-11-12 | 中国石油天然气股份有限公司 | 单河道砂体的识别方法、装置和计算机存储介质 |
CN113255155A (zh) * | 2021-06-11 | 2021-08-13 | 中海石油(中国)有限公司 | 一种储层不连续界限的阻渗能力等级的评价方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101408624A (zh) * | 2007-10-08 | 2009-04-15 | 陶庆学 | 三维地震最佳时窗河道砂体储层预测评价技术 |
CN103675941A (zh) * | 2012-09-26 | 2014-03-26 | 中国石油天然气集团公司 | 一种获得复杂储层的砂体三维形态分布的方法 |
CN104502969A (zh) * | 2014-12-30 | 2015-04-08 | 中国石油化工股份有限公司 | 河道砂岩性油藏的识别方法 |
CN104636980A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 |
-
2015
- 2015-08-11 CN CN201510490652.2A patent/CN105093303B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101408624A (zh) * | 2007-10-08 | 2009-04-15 | 陶庆学 | 三维地震最佳时窗河道砂体储层预测评价技术 |
CN103675941A (zh) * | 2012-09-26 | 2014-03-26 | 中国石油天然气集团公司 | 一种获得复杂储层的砂体三维形态分布的方法 |
CN104636980A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 |
CN104502969A (zh) * | 2014-12-30 | 2015-04-08 | 中国石油化工股份有限公司 | 河道砂岩性油藏的识别方法 |
Non-Patent Citations (1)
Title |
---|
陈飞 等: "渤海海域W油田新近系明化镇组河流相砂体结构特征", 《地学前缘(中国地质大学(北京);北京大学)》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107247289A (zh) * | 2017-06-07 | 2017-10-13 | 中海石油(中国)有限公司 | 一种河流相地层的中期沉积旋回末期等时界面恢复方法 |
CN107247289B (zh) * | 2017-06-07 | 2019-01-01 | 中海石油(中国)有限公司 | 一种河流相地层的中期沉积旋回末期等时界面恢复方法 |
CN107966731A (zh) * | 2017-11-08 | 2018-04-27 | 西南石油大学 | 一种基于地震波形结构属性的河流相砂体叠置区识别方法 |
CN107966731B (zh) * | 2017-11-08 | 2019-04-30 | 西南石油大学 | 一种基于地震波形结构属性的河流相砂体叠置区识别方法 |
CN110441818A (zh) * | 2019-07-31 | 2019-11-12 | 中国石油天然气股份有限公司 | 单河道砂体的识别方法、装置和计算机存储介质 |
CN110441818B (zh) * | 2019-07-31 | 2021-03-30 | 中国石油天然气股份有限公司 | 单河道砂体的识别方法、装置和计算机存储介质 |
CN113255155A (zh) * | 2021-06-11 | 2021-08-13 | 中海石油(中国)有限公司 | 一种储层不连续界限的阻渗能力等级的评价方法 |
CN113255155B (zh) * | 2021-06-11 | 2024-05-28 | 中海石油(中国)有限公司 | 一种储层不连续界限的阻渗能力等级的评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105093303B (zh) | 2017-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104237964B (zh) | 古河道型砂岩铀矿定位方法 | |
CN102466815B (zh) | 三叠系碎屑岩油气藏识别方法 | |
CN105093303A (zh) | 一种河流相砂体不连续界线精细表征的地震解释方法 | |
CN104965979A (zh) | 一种致密砂岩有效储层识别方法 | |
CN104111483A (zh) | 一种断陷湖盆斜坡带油藏分布的勘探方法 | |
CN104020509B (zh) | 基于Bayes判别分析的砂质滩坝沉积微相测井识别方法 | |
CN105866832A (zh) | 深层致密砂岩气藏的低级序断层的识别方法和装置 | |
CN103744109A (zh) | 一种无井覆盖区碎屑岩风化壳结构的识别方法 | |
CN109633743A (zh) | 一种基于波形分类地震相技术预测煤层厚度的方法 | |
CN102454400A (zh) | 碳酸盐岩缝洞型储层识别方法 | |
CN110927817B (zh) | 一种薄层叠置河道水平井轨迹优化设计及控制的方法 | |
CN109239783A (zh) | 一种井震结合中取舍地震信息的方法 | |
CN109324171A (zh) | 一种基于岩性统计的沉积相定量识别方法 | |
RU2493362C1 (ru) | Способ разработки нефтяной залежи | |
CN111830558A (zh) | 一种断裂带雕刻方法 | |
CN109581489B (zh) | 嵌套式地震相的提取方法及系统 | |
CN112147714A (zh) | 河口坝内部结构分析方法 | |
CN109085646A (zh) | 一种基于eps相位体属性的三角洲相叠置砂体识别方法 | |
CN113805227B (zh) | 一种隐伏地层尖灭线识别方法 | |
CN115932967B (zh) | 一种基于蚂蚁体算法的裂缝识别方法 | |
Jiang et al. | The Application of Seismic Technology For Identification of The Small Faults And Analysis of The Remaining Oil Area in Mature Oilfield: A Case Study of DaqingChangyuan Oilfield | |
Liu | Discussion on reservoir description method of well vibration combined with tight well pattern area | |
Roberto et al. | New oil discoveries in the Ceuta area, SE Lake Maracaibo, Venezuela | |
Wang et al. | Control function of Sedimentary Facies on Tight Oil Distribution—Taking Gaotaizi Oil Layer in Qijia Area as An Example | |
Li et al. | Application of high-definition maximum positive and negative curvature in fault opening description of shale oil reservoir in Jimsar Sag, Junggar Basin, Xinjiang |
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 | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC Research Institute Patentee before: China National Offshore Oil Corporation |
|
CP01 | Change in the name or title of a patent holder |