CN116953797A - 基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法 - Google Patents
基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法 Download PDFInfo
- Publication number
- CN116953797A CN116953797A CN202311019363.5A CN202311019363A CN116953797A CN 116953797 A CN116953797 A CN 116953797A CN 202311019363 A CN202311019363 A CN 202311019363A CN 116953797 A CN116953797 A CN 116953797A
- Authority
- CN
- China
- Prior art keywords
- frequency
- typical high
- gyratory
- thin
- sandstone reservoir
- 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
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000011161 development Methods 0.000 claims abstract description 20
- 238000005553 drilling Methods 0.000 claims abstract description 19
- 238000010586 diagram Methods 0.000 claims abstract description 17
- 238000009826 distribution Methods 0.000 claims abstract description 17
- 230000008021 deposition Effects 0.000 claims abstract description 16
- 239000011435 rock Substances 0.000 claims description 33
- 230000007480 spreading Effects 0.000 claims description 12
- 238000003892 spreading Methods 0.000 claims description 12
- 238000011160 research Methods 0.000 claims description 10
- 239000004576 sand Substances 0.000 claims description 9
- 238000010587 phase diagram Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 230000004044 response Effects 0.000 abstract description 2
- 230000009466 transformation Effects 0.000 abstract description 2
- 230000015572 biosynthetic process Effects 0.000 description 16
- 238000000151 deposition Methods 0.000 description 15
- 239000006071 cream Substances 0.000 description 7
- 239000010440 gypsum Substances 0.000 description 7
- 229910052602 gypsum Inorganic materials 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 6
- 235000002639 sodium chloride Nutrition 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 150000003839 salts Chemical class 0.000 description 3
- 239000006028 limestone Substances 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 239000011505 plaster Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005137 deposition process Methods 0.000 description 1
- 235000021185 dessert Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
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
- 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
-
- 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. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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
本发明涉及基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法,包括步骤:根据钻录井取心资料、测井曲线资料及二维地震剖面特征,建立等时地层格架;建立四级层序地层格架;识别典型高频旋回,选择薄层砂岩储层发育的典型高频旋回类型,建立薄层砂岩储层典型高频旋回识别模板;划分钻井中薄层砂岩储层典型高频旋回,统计典型高频旋回中砂岩储层厚度数据,建立典型高频旋回对比剖面和典型高频旋回控制下的沉积相平面图;利用Double FoxDraw软件开展湖相薄层砂岩储层预测。本发明利用高频旋回对沉积环境转化的响应关系,建立优势储层发育的典型高频旋回识别模板,追踪识别区域内优质储层分布特征,完成湖相薄层砂岩储层预测。
Description
技术领域
本发明涉及数据预测技术领域,特别涉及一种基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法。
背景技术
湖盆是地球上重要的油气资源聚集区之一,在油气勘探和开发过程中,湖泊三角洲一直是陆相油气储存的主要部位之一。受构造演化、气候条件和物源供给的影响,湖盆沉积储层大多为薄层的砂岩储层,由于其厚度薄、分布范围小、连通性差,在过去通常不被认为是优质的储层,勘探开发力度相对较小。但是随着油气成藏理论和勘探开发技术的发展,薄层湖相砂岩以其砂岩纯度高、粒度分选合适、上覆下伏泥岩封隔作用强、易形成“甜点”式油气藏等特征,开始成为今年来勘探开发的重点。
薄层湖相砂岩储层因其厚度薄、延伸范围短、沙体尖灭快等特征,在过去常用的地震反射剖面、地震反演等技术精度不足以完成相关的勘探预测要求,存在一些局限性,如数据获取成本高、时间和空间分辨较低、勘探精度低、预测结果不准确、反演结果多解性强等问题。随着高频旋回技术的发展,研究者开始探索更加高精度的方法去实现湖相沉积过程恢复,而基于此衍生出的等时对比地层层序格架,让湖相沉积的研究站在了等时界面上。高频旋回是湖泊等水体中由风力、水流和波浪等作用形成的小尺度地形特征。这些旋回特征的形成与湖盆沉积、水动力条件密切相关。湖盆中的湖泊地貌、沉积物输运和沉积作用等过程会留下高频旋回的记录。
因此,通过分析高频旋回特征,揭示湖盆演化过程以及沉积体系的动力学特征,从而达到高频旋回等时地层对比下的储能预测,可以为油气勘探和开发提供重要的地质依据。高频旋回的应用有助于揭示湖盆演化的复杂性,为油气资源的勘探和开发提供更准确的评价和预测。
发明内容
本发明的目的在于以层序地层学为理论依据,充分利用丰富的基础数据资料,提供一种基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法,追踪识别区域内优质储层分布特征,完成湖相薄层砂岩储层预测。
为了实现上述发明目的,本发明实施例提供了以下技术方案:
基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法,包括以下步骤:
步骤1,根据研究区钻录井取心资料、测井曲线资料及二维地震剖面特征,识别等时地层格架标志层、等时地层格架顶界面和等时地层格架底界面,从而建立等时地层格架;
步骤2,在步骤1建立的等时地层格架的约束下,建立四级层序地层格架;
步骤3,基于步骤2建立的四级层序地层格架,识别典型高频旋回,选择薄层砂岩储层发育的典型高频旋回类型,建立薄层砂岩储层典型高频旋回识别模板;
步骤4,根据步骤3建立的薄层砂岩储层典型高频旋回识别模板,划分钻井中薄层砂岩储层典型高频旋回,统计典型高频旋回中砂岩储层厚度数据,建立典型高频旋回对比剖面和典型高频旋回控制下的沉积相平面图;
步骤5,根据步骤4建立的典型高频旋回对比剖面和典型高频旋回控制下的沉积相平面图,利用Double FoxDraw软件开展湖相薄层砂岩储层预测。
所述步骤1具体包括以下步骤:
步骤1-1,基于伽马与电阻率、声波时差的测井曲线资料,结合钻井岩心资料,进行岩心归位与深度矫正,明确岩性与测井曲线的匹配关系;
步骤1-2,基于二维地震剖面特征,识别连续反射界面为等时地层格架顶界面、等时地层格架底界面,以及识别等时地层格架标志层;
步骤1-3,基于等时地层格架顶界面和等时地层格架底界面之间的地层为等时地层格架。
所述步骤2具体包括以下步骤:
步骤2-1,基于步骤1中等时地层格架顶界面、等时地层格架底界面以及等时地层格架标志层的约束,通过对岩石类型识别和测井曲线突变峰值界面,识别三级层序地层界面;
步骤2-2,基于识别的三级层序地层界面,建立三级层序地层格架;
步骤2-3,基于建立的三级层序地层格架,识别四级层序地层界面;
步骤2-4,基于识别的四级层序地层界面,建立四级层序地层格架。
所述步骤3具体包括以下步骤:
步骤3-1,在步骤2建立的四级层序地层格架中,根据岩性组合、岩石类型、测井曲线特征识别划分典型的高频旋回类型;
步骤3-2,在识别出的典型高频旋回类型中,选择具有薄层砂岩储层岩性组合、岩石类型和测井曲线特征的高频旋回类型,确定该类典型高频旋回发育的四级层序;
步骤3-3,根据选择出的薄层砂岩储层典型高频旋回,建立薄层砂岩储层典型高频旋回识别模板。
所述步骤4具体包括以下步骤:
步骤4-1,根据步骤3建立的典型高频旋回序列识别模板,划分研究区钻井中薄层砂岩储层典型高频旋回,并统计薄层砂岩储层典型高频旋回中砂岩厚度数据;
步骤4-2,在步骤2识别的四级层序地层格架顶界面拉平,在步骤3中识别的具有薄层砂岩储层典型高频旋回的四级层序地层格架建立薄层砂岩储层典型高频旋回对比剖面;
步骤4-3,根据步骤4-1统计的典型高频旋回中砂岩储层厚度数据,绘制研究区发育薄层砂岩储层典型高频旋回的四级层序地层格架沉积相平面图。
所述步骤5具体包括以下步骤:
步骤5-1,根据步骤4中统计的典型高频旋回中砂岩储层厚度数据和建立的薄层砂岩储层典型高频旋回对比剖面,圈定薄层砂岩储层分布范围;
步骤5-2,根据步骤4绘制的发育薄层砂岩储层典型高频旋回四级层序地层格架沉积相平面图,明确砂岩储层分布规律与展步方向;
步骤5-3,利用Double FoxDraw软件在圈定的薄层砂岩储层分布范围内,结合砂岩储层分布规律和展布方向,绘制由钻井控制的带有砂体展布规律的等值线图;
步骤5-4,根据绘制的带有砂体展布规律的等值线图导出带名称的数据体;
步骤5-5,利用Double FoxDraw软件的等值线勾画模块,采用100×50网格密度,加密次数为2,收敛精度为0.0015自动插值,处理并生成等值线间隔为1,步长为5的砂岩储层厚度等值线图,完成储层预测。
与现有技术相比,本发明的有益效果:
本发明建立等时地层格架,确保了高频旋回的可对比性,充分利用高频旋回对沉积环境转化的响应关系,从而建立优势储层发育的典型高频旋回识别模板,追踪识别区域内优质储层分布特征,完成湖相薄层砂岩储层预测。
本发明操作简单、步骤清晰、且识别模板具有区域可对比性,关键操作步骤做到定量化,提高湖相薄层砂岩储层预测尺度与精度,为深化地质认识、精确指导油气勘探开发提供充足的理论依据,降低勘探风险和成本投入。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明实施例等时地层格架顶、底界面和标志层图;
图2为本发明实施例三级、四级层序地层界面图;
图3为本发明实施例典型高频旋回序列类型图;
图4为本发明实施例钻井中薄层砂岩储层典型高频旋回识别图;
图5为本发明实施例薄层砂岩储层典型高频旋回对比剖面图;
图6为本发明实施例四级层序格架Esq1沉积相平面图;
图7为本发明实施例四级层序地层格架Esq1中薄层砂岩储层预测厚度等值线图。
具体实施方式
下面将结合本发明实施例中附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。同时,在本发明的描述中,术语“第一”、“第二”等仅用于区分描述,而不能理解为指示或暗示相对重要性,或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。另外,术语“相连”、“连接”等可以是元件之间直接相连,也可以是经由其他元件的间接相连。
实施例
本实施例以塔里木盆地新和地区古近系库姆格列木群地层为例,本发明通过下述技术方案实现,如图1所示,提出一种基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法,包括以下步骤:
步骤1,根据研究区钻录井取心资料、测井曲线资料及二维地震剖面特征,识别等时地层格架标志层、等时地层格架顶界面和等时地层格架底界面,从而建立等时地层格架。
首先,基于伽马与电阻率、声波时差等测井曲线资料,结合新和地区钻井岩心资料,测井曲线与岩性对应关系体现为砂岩对应低伽马、低声波时差曲线和高电阻率,泥岩对应高伽马、高声波时差曲线和低电阻率,膏岩对应低伽马、高声波时差和高电阻率,通过上述对应关系进行岩心归位与深度矫正,明确岩性与测井曲线的匹配关系。
接着,基于二维地震剖面特征,识别连续强反射界面为等时地层界面,在本实施例中,如图1所示,识别强反射轴T2 4作为等时地层格架顶界面,识别强反射轴T3 0作为等时地层格架底界面。
请继续参见图1,基于深度矫正后的岩心资料、测井曲线资料和二维地震剖面特征识别的等时地层格架顶、底界面,识别出等时地层格架标志层(以下简称标志层)。标志层以广泛分布、分布厚度均匀、特征明确的岩层为主,具备等时性,是等时地层格架建立的主要标志,本实施例识别出3套标志层,分别为标志层1、标志层2、标志层3,具体的识别特征如下:
①标志层1:以识别出来的等时地层格架底界面T3 0向上,位于库木格列姆群下部,岩性为碳酸盐岩,具有低伽马、高电阻率测井曲线特征的岩层识别为标志层1;
②标志层2:以识别出来的标志层1向上,位于库木格列姆群中部,岩性主要为膏盐岩,具有低伽马、高电阻率特征,二维地震剖面特征显示为强连续反射,同相轴交叉合并,具有膏岩的反射特征的岩层识别为标志层2;
③标志层3:以识别出来的标志层2向上,位于库木格列姆群上部,岩性主要为膏质泥岩、泥岩,具有高伽马、低电阻率特征的岩层识别为标志层3。
如图1所示,为本实施例识别出来的等时地层格架标志层、等时地层格架顶界面T2 4和等时地层格架底界面T3 0,基于T2 4和T3 0之间的地层为等时地层格架。
步骤2,在步骤1建立的等时地层格架的约束下,建立四级层序地层格架。
基于步骤1中等时地层格架顶界面T2 4、等时地层格架底界面T3 0以及标志层的约束,通过对岩石类型识别和测井曲线突变峰值界面,识别三级层序地层界面。具体来说,本实施例中三级层序地层界面主要有岩性转换界面、不整合面,共识别出3个三级层序地层界面,分别为Ⅰ号三级层序地层界面、Ⅱ号三级层序地层界面、Ⅲ号三级层序地层界面,具体的识别特征如下:
①Ⅰ号三级层序地层界面:将步骤1中识别的等时地层格架底界面T3 0识别为Ⅰ号三级层序地层界面;
②Ⅱ号三级层序地层界面:将步骤1中识别的标志层2的底部识别为Ⅱ号三级层序地层界面;
③Ⅲ号三级层序地层界面:将步骤1中识别的等时地层格架顶界面T2 4识别为Ⅲ号三级层序地层界面。
请参见图2,基于识别的Ⅰ号三级层序地层界面、Ⅱ号三级层序地层界面、Ⅲ号三级层序地层界面建立三级层序地层格架,在本实施例中,共识别出2套三级层序地层,分别为ESQ1和ESQ2,具体特征如下:
①ESQ1:将识别的Ⅰ号三级层序地层界面与Ⅱ号三级层序地层界面之间的岩层识别为ESQ1;
②ESQ2:将识别的Ⅱ号三级层序地层界面与Ⅲ号三级层序地层界面之间的岩层识别为ESQ2。
基于建立的三级层序地层格架,识别四级层序地层界面,在本实施例中,共识别出8个四级层序地层界面,分别为ESB1、ESB2、ESB3、ESB4、ESB5、ESB6、ESB7、ESB8,具体特征如下:
①ESB1:将识别的Ⅰ号三级层序地层界面识别为ESB1;
②ESB2:将伽马值由低至高、电阻率由低至高,岩性由膏盐岩转变为含膏砂岩的岩性转换面识别为ESB2;
③ESB3:将识别的Ⅱ号三级层序地层界面识别为ESB3;
④ESB4:将伽马值由低至高、电阻率由高至低,岩性由膏盐岩转变为泥岩的岩性转换面识别为ESB4;
⑤ESB5:将伽马值由低至高、电阻率由高至低,岩性由含膏泥岩转变为粉砂岩的岩性转换面识别为ESB5;
⑥ESB6:将伽马值由低至高、电阻率由高至低,岩性由含膏泥岩转变为泥岩的岩性转换面识别为ESB6;
⑦ESB7:将伽马值由低至高、电阻率由高至低,岩性由泥岩转变为膏盐的岩性转换面识别为ESB7;
⑧ESB8:将识别的Ⅲ号三级层序地层界面识别为ESB8。
基于识别的四级层序地层界面ESB1、ESB2、ESB3、ESB4、ESB5、ESB6、ESB7、ESB8,建立四级层序地层格架,请参见图2,共识别出7套四级层序地层格架,分别为Esq1、Esq2、Esq3、Esq4、Esq5、Esq6、Esq7,具体特征如下:
①Esq1:将识别的四级层序地层ESB1和四级层序地层ESB2之间的地层识别为Esq1;
②Esq2:将识别的四级层序地层ESB2和四级层序地层ESB3之间的地层识别为Esq2;
③Esq3:将识别的四级层序地层ESB3和四级层序地层ESB4之间的地层识别为Esq3;
④Esq4:将识别的四级层序地层ESB4和四级层序地层ESB5之间的地层识别为Esq4;
⑤Esq5:将识别的四级层序地层ESB5和四级层序地层ESB6之间的地层识别为Esq5;
⑥Esq6:将识别的四级层序地层ESB6和四级层序地层ESB7之间的地层识别为Esq6;
⑦Esq7:将识别的四级层序地层ESB7和四级层序地层ESB8之间的地层识别为Esq7。
步骤3,基于步骤2建立的四级层序地层格架,识别典型高频旋回,选择薄层砂岩储层发育的典型高频旋回类型,建立薄层砂岩储层典型高频旋回识别模板。
在步骤2建立的四级层序地层格架中,根据岩性组合、岩石类型、测井曲线特征识别划分典型的高频旋回类型,在本实施例中,共识别出4类典型高频旋回序列,请参见图3,为本实施例典型高频旋回序列类型图,包括典型高频旋回序列Ⅰ、典型高频旋回序列Ⅱ、典型高频旋回序列Ⅲ、典型高频旋回序列Ⅳ,具体特征如下:
①典型高频旋回序列Ⅰ:将自下而上岩性依次为粉砂岩-纹层状泥岩-灰质云岩-石膏岩的岩性组合识别为典型高频旋回序列Ⅰ;
②典型高频旋回序列Ⅱ:将自下而上岩性依次为粉砂岩-薄层状泥岩-膏质泥岩/膏质砂岩的岩性组合识别为典型高频旋回序列Ⅱ;
③典型高频旋回序列Ⅲ:将自下而上岩性依次为砂岩/泥岩-膏盐岩的岩性组合识别为典型高频旋回序列Ⅲ;
④典型高频旋回序列Ⅳ:将自下而上岩性依次为砂岩/泥岩-含膏质岩(包括石膏岩和膏质泥岩)的岩性组合识别为典型高频旋回序列Ⅳ。
在识别出的典型高频旋回类型中,选择具有薄层砂岩储层岩性组合、岩石类型和测井曲线特征的高频旋回类型,确定该类典型高频旋回发育的四级层序,在本实施例中,具有薄层砂岩储层典型高频旋回为典型高频旋回序列Ⅰ和典型高频旋回序列Ⅳ,主要分布在四级层序地层格架Esq1、Esq2、Esq5、Esq6中。
根据选择出的薄层砂岩储层典型高频旋回,建立薄层砂岩储层典型高频旋回识别模板,在本实施例中,建立两套具有薄层砂岩储层典型高频旋回识别模板,具体特征如下:
①典型高频旋回序列Ⅰ识别模板:将具备自下而上岩性依次为粉砂岩-纹层状泥岩-灰质云岩-石膏岩的岩性组合,岩石类型以粉砂岩、泥岩、云岩和石膏岩为主,且自然伽马和声波时差测井曲线以箱型、钟型为主特征的序列识别为典型高频旋回序列Ⅰ识别模板;
②典型高频旋回序列Ⅳ识别模板:将具备自下而上岩性依次为砂岩/泥岩-含膏质岩(包括石膏岩和膏质泥岩)的岩性组合,岩石类型以粉砂岩、膏盐岩为主,且自然伽马和声波时差测井曲线以钟型、齿化钟型为主特征的序列识别为典型高频旋回序列Ⅳ识别模板。
步骤4,根据步骤3建立的薄层砂岩储层典型高频旋回识别模板,划分钻井中薄层砂岩储层典型高频旋回,统计典型高频旋回中砂岩储层厚度数据,建立典型高频旋回对比剖面和典型高频旋回控制下的沉积相平面图。
根据步骤3建立的典型高频旋回序列识别模板,划分研究区钻井中薄层砂岩储层典型高频旋回,并统计薄层砂岩储层典型高频旋回中砂岩厚度数据,如图4为本实施例钻井中薄层砂岩储层典型高频旋回识别图。
在步骤2识别的四级层序地层格架顶界面拉平,在步骤3中识别的具有薄层砂岩储层典型高频旋回的四级层序地层格架Esq1、Esq2、Esq5、Esq6,建立薄层砂岩储层典型高频旋回对比剖面,如图5为本实施例中薄层砂岩储层典型高频旋回对比剖面图。
根据统计的典型高频旋回中砂岩储层厚度数据,绘制研究区发育薄层砂岩储层典型高频旋回的四级层序地层格架沉积相平面图,在本实施例中,采取传统沉积学绘制沉积相方法,以典型高频旋回发育的四级层序地层格架中优势相为指导,钻井沉积相为控制点,薄层砂岩储层典型高频旋回对比剖面为控制线,绘制研究区发育薄层砂岩储层典型高频旋回的四级层序地层格架沉积相平面图,如图6为本实施例中四级层序地层格架Esq1沉积相平面图。
步骤5,根据步骤4建立的典型高频旋回对比剖面和典型高频旋回控制下的沉积相平面图,统计砂岩储层厚度数据,利用Double FoxDraw软件开展湖相薄层砂岩储层预测。
根据步骤4中统计的典型高频旋回中砂岩储层厚度数据和建立的薄层砂岩储层典型高频旋回对比剖面,圈定薄层砂岩储层分布范围。
根据步骤4绘制的发育薄层砂岩储层典型高频旋回四级层序地层格架沉积相平面图,明确砂岩储层分布规律与展步方向。
利用Double FoxDraw软件在圈定的薄层砂岩储层分布范围内,结合砂岩储层分布规律和展布方向,绘制由钻井控制的带有砂体展布规律的等值线图。根据绘制的带有砂体展布规律的等值线图导出带名称的数据体。
利用Double FoxDraw软件的等值线勾画模块,采用100×50网格密度,加密次数为2,收敛精度为0.0015自动插值,处理并生成等值线间隔为1,步长为5的砂岩储层厚度等值线图,完成储层预测。如图7为本实施例四级层序地层格架Esq1中薄层砂岩储层预测厚度等值线图。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (6)
1.基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法,其特征在于:包括以下步骤:
步骤1,根据研究区钻录井取心资料、测井曲线资料及二维地震剖面特征,识别等时地层格架标志层、等时地层格架顶界面和等时地层格架底界面,从而建立等时地层格架;
步骤2,在步骤1建立的等时地层格架的约束下,建立四级层序地层格架;
步骤3,基于步骤2建立的四级层序地层格架,识别典型高频旋回,选择薄层砂岩储层发育的典型高频旋回类型,建立薄层砂岩储层典型高频旋回识别模板;
步骤4,根据步骤3建立的薄层砂岩储层典型高频旋回识别模板,划分钻井中薄层砂岩储层典型高频旋回,统计典型高频旋回中砂岩储层厚度数据,建立典型高频旋回对比剖面和典型高频旋回控制下的沉积相平面图;
步骤5,根据步骤4建立的典型高频旋回对比剖面和典型高频旋回控制下的沉积相平面图,利用Double FoxDraw软件开展湖相薄层砂岩储层预测。
2.根据权利要求1所述的基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法,其特征在于:所述步骤1具体包括以下步骤:
步骤1-1,基于伽马与电阻率、声波时差的测井曲线资料,结合钻井岩心资料,进行岩心归位与深度矫正,明确岩性与测井曲线的匹配关系;
步骤1-2,基于二维地震剖面特征,识别连续反射界面为等时地层格架顶界面、等时地层格架底界面,以及识别等时地层格架标志层;
步骤1-3,基于等时地层格架顶界面和等时地层格架底界面之间的地层为等时地层格架。
3.根据权利要求2所述的基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法,其特征在于:所述步骤2具体包括以下步骤:
步骤2-1,基于步骤1中等时地层格架顶界面、等时地层格架底界面以及等时地层格架标志层的约束,通过对岩石类型识别和测井曲线突变峰值界面,识别三级层序地层界面;
步骤2-2,基于识别的三级层序地层界面,建立三级层序地层格架;
步骤2-3,基于建立的三级层序地层格架,识别四级层序地层界面;
步骤2-4,基于识别的四级层序地层界面,建立四级层序地层格架。
4.根据权利要求3所述的基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法,其特征在于:所述步骤3具体包括以下步骤:
步骤3-1,在步骤2建立的四级层序地层格架中,根据岩性组合、岩石类型、测井曲线特征识别划分典型的高频旋回类型;
步骤3-2,在识别出的典型高频旋回类型中,选择具有薄层砂岩储层岩性组合、岩石类型和测井曲线特征的高频旋回类型,确定该类典型高频旋回发育的四级层序;
步骤3-3,根据选择出的薄层砂岩储层典型高频旋回,建立薄层砂岩储层典型高频旋回识别模板。
5.根据权利要求4所述的基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法,其特征在于:所述步骤4具体包括以下步骤:
步骤4-1,根据步骤3建立的典型高频旋回序列识别模板,划分研究区钻井中薄层砂岩储层典型高频旋回,并统计薄层砂岩储层典型高频旋回中砂岩厚度数据;
步骤4-2,在步骤2识别的四级层序地层格架顶界面拉平,在步骤3中识别的具有薄层砂岩储层典型高频旋回的四级层序地层格架建立薄层砂岩储层典型高频旋回对比剖面;
步骤4-3,根据步骤4-1统计的典型高频旋回中砂岩储层厚度数据,绘制研究区发育薄层砂岩储层典型高频旋回的四级层序地层格架沉积相平面图。
6.根据权利要求5所述的基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法,其特征在于:所述步骤5具体包括以下步骤:
步骤5-1,根据步骤4中统计的典型高频旋回中砂岩储层厚度数据和建立的薄层砂岩储层典型高频旋回对比剖面,圈定薄层砂岩储层分布范围;
步骤5-2,根据步骤4绘制的发育薄层砂岩储层典型高频旋回四级层序地层格架沉积相平面图,明确砂岩储层分布规律与展步方向;
步骤5-3,利用Double FoxDraw软件在圈定的薄层砂岩储层分布范围内,结合砂岩储层分布规律和展布方向,绘制由钻井控制的带有砂体展布规律的等值线图;
步骤5-4,根据绘制的带有砂体展布规律的等值线图导出带名称的数据体;
步骤5-5,利用Double FoxDraw软件的等值线勾画模块,采用100×50网格密度,加密次数为2,收敛精度为0.0015自动插值,处理并生成等值线间隔为1,步长为5的砂岩储层厚度等值线图,完成储层预测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311019363.5A CN116953797B (zh) | 2023-08-14 | 2023-08-14 | 基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311019363.5A CN116953797B (zh) | 2023-08-14 | 2023-08-14 | 基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116953797A true CN116953797A (zh) | 2023-10-27 |
CN116953797B CN116953797B (zh) | 2024-07-02 |
Family
ID=88449243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311019363.5A Active CN116953797B (zh) | 2023-08-14 | 2023-08-14 | 基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116953797B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120075951A1 (en) * | 2010-09-29 | 2012-03-29 | Chevron U.S.A., Inc. | Device and method for imaging of non-linear and linear properties of formations surrounding a borehole |
CN104141481A (zh) * | 2013-05-06 | 2014-11-12 | 中国石油天然气股份有限公司 | 一种超低渗透致密油藏水平井布井方法 |
CN105089662A (zh) * | 2015-07-24 | 2015-11-25 | 中国石油天然气股份有限公司 | 碳酸盐岩地层录井综合柱状图的校正方法及系统 |
CN108873070A (zh) * | 2017-05-12 | 2018-11-23 | 中国石油化工股份有限公司 | 恢复盆地内部地层古产状的方法 |
CN112832301A (zh) * | 2020-12-31 | 2021-05-25 | 成都理工大学 | 一种基于柱坐标系的电磁灌注桩、预制桩检测方法 |
-
2023
- 2023-08-14 CN CN202311019363.5A patent/CN116953797B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120075951A1 (en) * | 2010-09-29 | 2012-03-29 | Chevron U.S.A., Inc. | Device and method for imaging of non-linear and linear properties of formations surrounding a borehole |
CN104141481A (zh) * | 2013-05-06 | 2014-11-12 | 中国石油天然气股份有限公司 | 一种超低渗透致密油藏水平井布井方法 |
CN105089662A (zh) * | 2015-07-24 | 2015-11-25 | 中国石油天然气股份有限公司 | 碳酸盐岩地层录井综合柱状图的校正方法及系统 |
CN108873070A (zh) * | 2017-05-12 | 2018-11-23 | 中国石油化工股份有限公司 | 恢复盆地内部地层古产状的方法 |
CN112832301A (zh) * | 2020-12-31 | 2021-05-25 | 成都理工大学 | 一种基于柱坐标系的电磁灌注桩、预制桩检测方法 |
Non-Patent Citations (1)
Title |
---|
朱茂,等: "四川盆地潼南—合川区块中二叠统白云岩储层形成主控因素与勘探区带预测", 中国石油勘探, vol. 27, no. 4, 31 July 2022 (2022-07-31), pages 149 - 161 * |
Also Published As
Publication number | Publication date |
---|---|
CN116953797B (zh) | 2024-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111766630B (zh) | 一种盆地古地貌的恢复方法 | |
Parolai et al. | New relationships between Vs, thickness of sediments, and resonance frequency calculated by the H/V ratio of seismic noise for the Cologne area (Germany) | |
CN111257926B (zh) | 一种利用老旧地震资料进行古河谷铀储层预测的方法 | |
CN105899749A (zh) | 偏斜井眼和水平井眼的地层和结构解释 | |
CN103901478A (zh) | 一种井震信息联合确定储层沉积特征和分布的方法 | |
CN109541685B (zh) | 一种河道砂体识别方法 | |
CN110056346B (zh) | 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 | |
CN107133670A (zh) | 一种基于决策树数据挖掘算法的复杂岩性识别方法及系统 | |
CN106772587A (zh) | 基于同位多相协同克里金的地震弹性参数相控建模方法 | |
CN109164487A (zh) | 一种基于模型建立平均速度场的方法及精细成图的方法 | |
CN111983683B (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
CN105549084A (zh) | 一种三维高精度速度建模方法和系统 | |
CN108121009B (zh) | 开发后期复杂断块变速构造成图方法 | |
CN106526671A (zh) | 一种碳酸盐岩储层地震波响应特征分析方法 | |
WO2023124912A1 (zh) | 碳酸盐岩沉积相类型的预测方法及装置 | |
CN116953797B (zh) | 基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法 | |
CN111413737A (zh) | 一种密井网区地震层位解释到沉积单元的方法 | |
CN113885096A (zh) | 高频层序划分与小层对比方法、装置、电子设备及介质 | |
CN112987126A (zh) | 一种岩石物理相的地球物理预测方法 | |
CN113805227B (zh) | 一种隐伏地层尖灭线识别方法 | |
CN115857015B (zh) | 一种火山地层中定量预测沉凝灰岩分布的方法 | |
CN108732621A (zh) | 一种基于ffc-电阻率的随钻精细时深预测方法 | |
CN118068418A (zh) | 一种用于密井网区的地层对比方法 | |
CN116184525A (zh) | 一种地层格架下海陆过渡相页岩三维岩相建模方法 | |
CN117950017A (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 |