CN107870366B - 开发动态硬约束的砂砾岩连通体地震预测方法 - Google Patents
开发动态硬约束的砂砾岩连通体地震预测方法 Download PDFInfo
- Publication number
- CN107870366B CN107870366B CN201610843608.XA CN201610843608A CN107870366B CN 107870366 B CN107870366 B CN 107870366B CN 201610843608 A CN201610843608 A CN 201610843608A CN 107870366 B CN107870366 B CN 107870366B
- Authority
- CN
- China
- Prior art keywords
- connected component
- glutenite
- well
- prediction method
- hard constraint
- 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
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Abstract
本发明提供一种开发动态硬约束的砂砾岩连通体地震预测方法,该开发动态硬约束的砂砾岩连通体地震预测方法包括:步骤1,搜集前期基础资料;步骤2,进行连通体边界判识;步骤3,进行连通体波阻抗特征判别;步骤4,进行开发动态分析;步骤5,建立初始连通体模型;步骤6,进行井约束的地震反演;步骤7,利用反演波阻抗体进行连通体刻画。该开发动态硬约束的砂砾岩连通体地震预测方法通过分析已有井的开发注采对应情况,寻找可以确定的井间连通关系,以这些连通关系为硬约束条件进行地震反演,得到连通体的空间分布,有效的利用了开发动态信息,预测结果能够与开发井注采关系吻合。
Description
技术领域
本发明涉及油田开发地震研究领域,特别是涉及到一种开发动态硬约束的砂砾岩连通体地震预测方法。
背景技术
沉积凹陷北部陡坡带一般发育有大套的砂砾岩沉积体,由于砂砾岩沉积离物源较近,储量丰都大,往往具有很大的开发潜力。但是砂砾岩沉积体不仅纵横向变化十分复杂,同时由于砂砾岩扇体相变快、岩性多样,从而导致岩电关系复杂、储层非均质性极强、储层连通关系呈空间非均质变化,造成含油差异及产能变化较大,严重影响开发效果。如何有效将叠后地震资料与开发生产数据相结合,进行砂砾岩连通体的预测一直是困扰地球物理学家的重要问题。为此我们发明了一种新的开发动态硬约束的砂砾岩连通体地震预测方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种以开发动态为硬约束条件的砂砾岩连通体地震预测方法。
本发明的目的可通过如下技术措施来实现:开发动态硬约束的砂砾岩连通体地震预测方法,该开发动态硬约束的砂砾岩连通体地震预测方法包括:步骤1,搜集前期基础资料;步骤2,进行连通体边界判识;步骤3,进行连通体波阻抗特征判别;步骤4,进行开发动态分析;步骤5,建立初始连通体模型;步骤6,进行井约束的地震反演;步骤7,利用反演波阻抗体进行连通体刻画。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,需要搜集的资料包括:录井数据、测井数据、砂砾岩期次层位、生产井与注水井射开层位、注水井注水层段、生产井采油层段、地震数据。
在步骤2中,根据开发层段的分布与地质认识确定连通体发育的纵向深度范围、根据录井结果并结合前期地质认识确定储层发育的横向范围,在这个数据体内进行地震反演。
在步骤3中,根据录井岩性选出不同深度段不同岩性对应的声波、密度曲线,计算对应的波阻抗值,绘制不同岩性的概率直方图,其中构成连通体的为物性较好、孔隙度与渗透率较高的岩性,通过直方图确定该岩性的波阻抗分布区间。
在步骤4中,挑选出相邻的一口注水井、一口生产井,寻找它们在同一期次内两个射开段,分析在相邻时期注水井注水速度与生产井液量或者含水的变化,如果在开始注水一段时间之后生产井液量或者含水上升则确定为两井连通,如果液量或含水无变化甚至降低,而周边其他井相同层段却受这口注水井注水影响,认为这两口井在这个深度段不连通。
在步骤5中,首先确定两口井上同一套期次内对应的个砂体储层,在确定注水连通的情况下,设定连通值为1;在两套砂体确定不连通的情况下,在井旁道连通值设置为1,随两口井距离变小连通值逐渐降低0,其他区域连通值为0-1之间随机数值,对区块内所有能够确定连通关系的井都进行相同操作,形成空间连通体概率模型。
在步骤6中,首先利用测井曲线插值生成初始波阻抗模型,采用蒙特卡洛随机模拟算法进行地震反演,反演计算过程中,随机参数选用初始连通体模型中的参数,实现连通关系对反演过程的约束作用,反演过程中在完成一定次数迭代之后根据波阻抗计算相邻井之间砂体的展布范围,判断反演数据体的连通性,如果与初始给定的连通关系不一致,则重新进行计算。
在步骤7中,已知能够成为连通体岩性的波阻抗范围,在反演得到的波阻抗体上镂空对应岩性波阻抗范围之外的空间,只保留对应波阻抗范围内的区域,所保留的即为预测的井间连通体。
本发明中的开发动态硬约束的砂砾岩连通体地震预测方法,包括前期基础资料搜集、连通体边界判识、连通体波阻抗特征判别、开发动态分析、初始连通体模型建立、井约束的地震反演、模型判别修正、连通体刻画7个步骤。通过分析已有井的开发注采对应情况,寻找可以确定的井间连通关系,以这些连通关系为硬约束条件进行地震反演,得到连通体的空间分布。这种约束反演方法有效的利用了开发动态信息,预测结果能够与开发井注采关系吻合。
附图说明
图1为本发明的开发动态硬约束的砂砾岩连通体地震预测方法的一具体实施例的流程图;
图2为本发明的一具体实施例中岩性-波阻抗直方图;
图3为本发明的一具体实施例中连通体模型示意图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的开发动态硬约束的砂砾岩连通体地震预测方法的流程图。
步骤一,前期基础资料搜集:
需要搜集的资料包括:录井数据、测井数据、砂砾岩期次层位、生产井与注水井射开层位、注水井注水层段、生产井采油层段、地震数据等。
步骤二,连通体边界判识:
根据开发层段的分布与地质认识确定连通体发育的纵向深度范围、根据录井结果并结合前期地质认识确定储层发育的横向范围,在这个数据体内进行地震反演。
步骤三,连通体波阻抗特征判别:
根据录井岩性选出不同深度段不同岩性对应的声波、密度曲线,计算对应的波阻抗值,绘制不同岩性的概率直方图,其中构成连通体的一般为物性较好、孔隙度与渗透率较高的岩性,通过直方图确定该岩性的波阻抗分布区间,如图2所示,虚线表砂砾、砂质砾岩等有效储层出现的频数,实线表无效储层如泥岩,砾岩出现的频数;波阻抗在10500到12700之间时有效储层多与无效储层,认为在这个区间内波阻抗反映有效储层。
步骤四,开发动态分析:
挑选出相邻的一口注水井、一口生产井,寻找它们在同一期次内两个射开段,分析在相邻时期注水井注水速度与生产井液量或者含水的变化,如果在开始注水一段时间之后生产井液量或者含水上升则可以确定为两井连通,如果液量或含水无变化甚至降低,而周边其他井相同层段却受这口注水井注水影响,可以认为这两口井在这个深度段不连通。
步骤五,初始连通体模型建立:
首先确定两口井上一套期次内对应的个砂体储层,在确定注水连通的情况下,设定连通值为1;在两套砂体确定不连通的情况下,在井旁连通值设置为1,随两口井距离变小连通值逐渐降低0,,其他区域连通值为0-1之间随机数值。对区块内所有能够确定连通关系的井都进行相同操作,形成空间连通体概率模型,如图3所示,连通砂体内部赋值为1;不连通砂砾在井点处赋值为1,向周围逐渐降低至0,变化趋势如图中颜色深度变化;除连通砂体与非连通砂体外其他区域网格点值在0到1之间随机分布。
步骤六,井约束的地震反演:
首先利用测井曲线插值生成有初始波阻抗模型,采用蒙特卡洛随机模拟算法进行地震反演,反演计算过程中,随机参数选用初始连通体模型中的参数,实现连通关系对反演过程的约束作用,反演过程中在完成一定次数迭代之后根据波阻抗计算相邻井之间砂体的展布范围,判断反演数据体的连通性,如果与初始给定的连通关系不一致,则重新进行计算。
步骤七,利用反演波阻抗体进行连通体刻画:
已知能够成为连通体岩性的波阻抗范围,在反演得到的波阻抗体上镂空对应岩性波阻抗范围之外的空间,只保留对应波阻抗范围内的区域,所保留的即为预测的井间连通体。
Claims (6)
1.开发动态硬约束的砂砾岩连通体地震预测方法,其特征在于,该开发动态硬约束的砂砾岩连通体地震预测方法包括:
步骤1,搜集前期基础资料;
步骤2,进行连通体边界判识;
步骤3,进行连通体波阻抗特征判别;
步骤4,进行开发动态分析;
步骤5,建立初始连通体模型;
步骤6,进行井约束的地震反演;
步骤7,利用反演波阻抗体进行连通体刻画;
在步骤5中,首先确定两口井上同一套期次内对应的各砂体储层,在确定注水连通的情况下,设定连通值为1;在两套砂体确定不连通的情况下,在井旁道连通值设置为1,随两口井距离变小连通值逐渐降低至 0,其他区域连通值为0-1之间随机数值,对区块内所有能够确定连通关系的井都进行相同操作,形成空间连通体概率模型;
在步骤6中,首先利用测井曲线插值生成初始波阻抗模型,采用蒙特卡洛随机模拟算法进行地震反演,反演计算过程中,随机参数选用初始连通体模型中的参数,实现连通关系对反演过程的约束作用,反演过程中在完成一定次数迭代之后根据波阻抗计算相邻井之间砂体的展布范围,判断反演数据体的连通性,如果与初始给定的连通关系不一致,则重新进行计算。
2.根据权利要求1所述的开发动态硬约束的砂砾岩连通体地震预测方法,其特征在于,在步骤1中,需要搜集的资料包括:录井数据、测井数据、砂砾岩期次层位、生产井与注水井射开层位、注水井注水层段、生产井采油层段、地震数据。
3.根据权利要求1所述的开发动态硬约束的砂砾岩连通体地震预测方法,其特征在于,在步骤2中,根据开发层段的分布与地质认识确定连通体发育的纵向深度范围、根据录井结果并结合前期地质认识确定储层发育的横向范围,在这个数据体内进行地震反演。
4.根据权利要求1所述的开发动态硬约束的砂砾岩连通体地震预测方法,其特征在于,在步骤3中,根据录井岩性选出不同深度段不同岩性对应的声波、密度曲线,计算对应的波阻抗值,绘制不同岩性的概率直方图,其中构成连通体的为物性较好、孔隙度与渗透率较高的岩性,通过直方图确定该岩性的波阻抗分布区间。
5.根据权利要求1所述的开发动态硬约束的砂砾岩连通体地震预测方法,其特征在于,在步骤4中,挑选出相邻的一口注水井、一口生产井,寻找它们在同一期次内两个射开段,分析在相邻时期注水井注水速度与生产井液量或者含水的变化,如果在开始注水一段时间之后生产井液量或者含水上升则确定为两井连通,如果液量或含水无变化甚至降低,而周边其他井相同层段却受这口注水井注水影响,认为这两口井在这个深度段不连通。
6.根据权利要求1所述的开发动态硬约束的砂砾岩连通体地震预测方法,其特征在于,在步骤7中,已知能够成为连通体岩性的波阻抗范围,在反演得到的波阻抗体上镂空对应岩性波阻抗范围之外的空间,只保留对应波阻抗范围内的区域,所保留的即为预测的井间连通体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610843608.XA CN107870366B (zh) | 2016-09-23 | 2016-09-23 | 开发动态硬约束的砂砾岩连通体地震预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610843608.XA CN107870366B (zh) | 2016-09-23 | 2016-09-23 | 开发动态硬约束的砂砾岩连通体地震预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107870366A CN107870366A (zh) | 2018-04-03 |
CN107870366B true CN107870366B (zh) | 2019-07-12 |
Family
ID=61751080
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610843608.XA Active CN107870366B (zh) | 2016-09-23 | 2016-09-23 | 开发动态硬约束的砂砾岩连通体地震预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107870366B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104915530A (zh) * | 2014-03-10 | 2015-09-16 | 中国石油化工股份有限公司 | 油藏井间连通关系的建立方法 |
CN105019894A (zh) * | 2015-07-29 | 2015-11-04 | 长江大学 | 一种多层油藏井间连通性模型建立方法及系统 |
CN105467464A (zh) * | 2015-08-14 | 2016-04-06 | 中国石油大学(华东) | 近岸水下扇砂砾岩有效连通体划分和对比方法 |
CN105717541A (zh) * | 2016-03-17 | 2016-06-29 | 成都创源油气技术开发有限公司 | 储层流动单元划分方法 |
-
2016
- 2016-09-23 CN CN201610843608.XA patent/CN107870366B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104915530A (zh) * | 2014-03-10 | 2015-09-16 | 中国石油化工股份有限公司 | 油藏井间连通关系的建立方法 |
CN105019894A (zh) * | 2015-07-29 | 2015-11-04 | 长江大学 | 一种多层油藏井间连通性模型建立方法及系统 |
CN105467464A (zh) * | 2015-08-14 | 2016-04-06 | 中国石油大学(华东) | 近岸水下扇砂砾岩有效连通体划分和对比方法 |
CN105717541A (zh) * | 2016-03-17 | 2016-06-29 | 成都创源油气技术开发有限公司 | 储层流动单元划分方法 |
Non-Patent Citations (1)
Title |
---|
塔中缝洞型碳酸盐岩储层连通性分析及应用;韩剑发等;《科学技术与工程》;20160229;第16卷(第5期);第148页 |
Also Published As
Publication number | Publication date |
---|---|
CN107870366A (zh) | 2018-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105041307B (zh) | 一种碎屑岩油气储层优势渗流通道识别的工业化流程 | |
CN105095986B (zh) | 多层油藏整体产量预测的方法 | |
CN104879103B (zh) | 一种分层注水效果分析方法 | |
RU2761946C1 (ru) | Способ и устройство для определения подхода к комплексной разработке сланца и соседних нефтяных коллекторов | |
WO2016070587A1 (zh) | 基于trap-3d软件寻找油气藏的方法 | |
CA2890817C (en) | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties | |
CN110130882A (zh) | 一种基于测井测试资料的油藏区域评价方法 | |
CN104100263B (zh) | 一种非均匀水淹层剩余油饱和度的确定方法 | |
CN105445800A (zh) | 一种厚层砂体顶部分异岩性油藏的识别方法 | |
Haddad et al. | Modelling of Cold Heavy Oil Production with Sand (CHOPS) using a fluidized sand algorithm | |
CN109212627B (zh) | 一种基于细分层系的薄储层描述及增储规模确定方法 | |
Dapeng et al. | An independent fracturing water-flooding development method for shallow low-permeability thin oil layers in multi-layer sandstone reservoirs | |
CN105404972A (zh) | 一种油藏开发不确定性研究及风险控制的方法 | |
Sadeghnejad et al. | Field scale characterization of geological formations using percolation theory | |
Suarez et al. | Fracturing-to-Production Simulation Approach for Completion Optimization in the Vaca Muerta Shale | |
Ghods et al. | Automatic estimation of fracture properties in multi-stage fractured shale gas horizontal wells for reservoir modeling | |
CN107766689B (zh) | 开发动态约束的储层渗透率时变模型的建立方法 | |
Liu et al. | The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs | |
CN107870366B (zh) | 开发动态硬约束的砂砾岩连通体地震预测方法 | |
Uddin et al. | Seismic correlated Mallik 3D gas hydrate distribution: effect of geomechanics in non-homogeneous hydrate dissociation by depressurization | |
Liu et al. | Effect of sandstone and mudstone thickness on artificial fracturing for hydrocarbon extraction from low-permeability reservoirs | |
CN105239998A (zh) | 一种基于海上油田储层细分的部井方法 | |
CN110569512B (zh) | 一种泥页岩中渗透性夹层发育程度定量表征方法 | |
Temizel et al. | Optimization of production in slanted horizontal wells in tight reservoirs | |
Lamberghini et al. | Performance Based Reservoir Characterization in a Tight Gas Reservoir System-Case Study from Lajas and Punta Rosada Formations in the Neuquen Basin, Argentina |
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 |