CN108254786B - 一种储层层序识别方法及装置 - Google Patents
一种储层层序识别方法及装置 Download PDFInfo
- Publication number
- CN108254786B CN108254786B CN201810001664.8A CN201810001664A CN108254786B CN 108254786 B CN108254786 B CN 108254786B CN 201810001664 A CN201810001664 A CN 201810001664A CN 108254786 B CN108254786 B CN 108254786B
- Authority
- CN
- China
- Prior art keywords
- unit
- sandy ground
- curve
- identifier
- ratio
- 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 49
- 238000012545 processing Methods 0.000 claims abstract description 62
- 239000004576 sand Substances 0.000 claims description 31
- 230000015572 biosynthetic process Effects 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 17
- 238000004458 analytical method Methods 0.000 claims description 12
- 239000011435 rock Substances 0.000 claims description 10
- 108010001267 Protein Subunits Proteins 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 11
- 238000009826 distribution Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 7
- 230000006872 improvement Effects 0.000 description 7
- 238000003860 storage Methods 0.000 description 7
- 238000012300 Sequence Analysis Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000009434 installation Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000001615 p wave Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000010348 incorporation Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000004062 sedimentation Methods 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 239000004317 sodium nitrate Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 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/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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
-
- 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/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6226—Impedance
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
本申请实施例公开了一种储层层序识别方法及装置,所述方法包括:获取岩性信息;对所述岩性信息进行数字化处理;根据所述数字化处理的岩性信息,确定单位砂地比曲线;对所述单位砂地比曲线进行第一处理,并根据所述第一处理结果进行层序识别。本申请实施例提供的储层层序识别的方法及装置,可以提高层序识别的准确性。
Description
技术领域
本申请涉及地震综合解释技术领域,特别涉及一种储层层序识别方法及装置。
背景技术
层序是由不整合面与不整合面相对应的整合面作为边界的、一个相对整合的、有内在联系的地层序列。一个层序可以是一级层序、二级层序和三级层序。准层序和准层序组(也称四级层序和五级层序)是层序的地层构成单元。准层序由一个湖(海)泛面或与之相对应的界面为边界的、相对整合的、有内在联系的岩层或岩层序列所组成。湖(海)泛面是一个将新老地层分开,其上下水深明显地急剧变化的一个界面。每一个层序、准层序都对应着相应的水体变化特征,这种水体变化特征体现在水平面上升或者下降,该时期某个位置所沉积的砂岩厚度与地层总厚度的相对关系也具有一定特点。
层序一般分为高水位体系域(高位域)、水进体系域、低水位体系域(低位域)。低水位体系域从不整合面开始沉积,水平面下降,沉积中心向盆地中心转移,同一地点随时间变化表现为砂岩逐渐变厚、变粗,占地层厚度比例增大。水进体系域从初始湖(海)泛面开始沉积,水平面上升,沉积中心向物源方向转移,同一地点随时间变化表现为砂岩逐渐变薄、变细,占地层厚度比例减小。高水位体系域从最大湖(海)泛面开始沉积,水平面下降,沉积中心向盆地中心转移,同一地点随时间变化表现为砂岩逐渐变厚、变粗,占地层厚度比例增大。
在地震资料综合解释领域的过程中,层序及准层序的识别是一项重要的工作。
现有的层序识别方法主要是通过分析岩性资料、测井曲线以及地震剖面中的响应特征来实现的,划分的准确性不高。例如,地层砂地比的对比分析方法中,因为砂地比是一个反映一段地层总体特征的结果,较难体现更精细的变化,在利用砂地比进行层序划分的过程中,一旦涉及到更加精细的的层序划分时往往划分结果不够准确。图1为某地区流一2亚段砂地比的平面分布图,图2为某地区流一3亚段砂地比的平面分布图,图中相同灰度所代表的砂地比的值是相同的。那么根据图1和图2,利用砂地比进行层序划分时,由于东南部从流一3亚段到流一2亚段地层砂地比增高,砂体远离物源方向沉积,可以判断这是一个水退的过程,故做出该地区流一2亚段和流一3亚段整体属于高位域的判断。但是流一2亚段、流一3亚段是两个相对较厚的地层,地层内部砂地比是如何变化的,从流一3亚段到流一2亚段的砂地比是如何过渡的,这些问题都无法得到精确的结果。
因此,发明人发现现有的通过分析岩性资料、测井曲线以及地震剖面中的响应特征来实现的层序识别方法,由于生产中层段特征不明显,往往划分结果不够准确。
发明内容
本申请实施例的目的是提供一种储层层序识别方法及装置,以提高层序识别的准确性。
为解决上述技术问题,本申请实施例提供一种储层层序识别的方法及装置是这样实现的:
一种储层层序识别方法,包括:
获取岩性信息;
对所述岩性信息进行数字化处理;
根据所述数字化处理的岩性信息,确定单位砂地比曲线;
对所述单位砂地比曲线进行第一处理,并根据所述第一处理结果进行层序识别。
优选方案中,所述岩性信息包括:单井砂岩分布深度。
优选方案中,所述获取岩性信息包括:从岩性录井资料获取所述单井砂岩分布深度,或者,根据现有的测井曲线获取所述单井砂岩分布深度。
优选方案中,所述对岩性信息进行数字化处理包括:
对所述岩性信息按照预设间隔进行采样,获取采样数据;
判断所述采样数据是否符合预设条件;
根据所述判断结果对所述采样数据进行标识。
优选方案中,所述预设条件为:该深度处的岩石为砂岩或砾岩。
优选方案中,所述对所述采样数据进行标识包括:将所述符合预设条件的采样数据标识为第一标识,将不符合预设条件的采样数据标识为第二标识。
优选方案中,所述第一标识为1;所述第二标识为0。
优选方案中,所述根据数字化处理的岩性信息确定单位砂地比曲线,包括:
以一个采样点为中心,统计所述地层厚度范围内第一标识的个数和第二标识的个数;
计算第一标识和第二标识的总个数;
将所述第一标识的个数与所述总个数的比值作为该采样点的砂地比;
将所述数字处理后的岩性信息对应的采样点的砂地比连接为曲线,即为单位砂地比曲线。
优选方案中,所述第一处理包括下述中的至少一项:单井分析、连井对比、反演。
一种储层层序识别装置,包括:信息获取单元、数字化处理单元、砂地比曲线确定单元
和识别单元;
所述信息获取单元,用于获取岩性信息;所述岩性信息包括:单井砂岩分布深度;
所述数字化处理单元,用于对所述岩性信息进行数字化处理;
所述砂地比曲线确定单元,用于根据所述数字化处理的岩性信息,确定单位砂地比曲线;
所述识别单元,用于对所述单位砂地比曲线进行第一处理,并根据所述第一处理结果进行层序识别。
优选方案中,所述数字化处理单元包括:采样子单元、判断子单元和标识子单元;
所述采样子单元,用于对所述岩性信息按照预设间隔进行采样,获取采样数据;
所述判断子单元,用于判断所述采样数据是否符合预设条件;所述预设条件为:该深度处的岩石为砂岩或砾岩;
所述标识子单元,用于根据所述判断结果对所述采样数据进行标识。
优选方案中,所述砂地比曲线确定单元包括:统计子单元、计算子单元、比值子单元和连接子单元;
所述统计子单元,用于以一个采样点为中心,统计所述地层厚度范围内第一标识的个数和第二标识的个数;
所述计算子单元,用于计算第一标识和第二标识的总个数;
所述比值子单元,用于将所述第一标识的个数与所述总个数的比值作为该采样点的砂地比;
所述连接子单元,用于将所述数字处理后的岩性信息对应的采样点的砂地比连接为曲线。
由以上本申请实施例提供的技术方案可见,本申请实施例提供的层序识别方法和装置,对于小尺度的单位砂地比曲线能够反映砂岩与非砂岩情况,储层反演中同相轴清晰,能够区分较薄的砂层;中小尺度的单位砂地比曲线,能够明显反映出砂层组合特征,便于连井层序对比;中大尺度的单位砂地比曲线,能够直观地体现水进水退特征,便于进行精细的层序对比;大尺度的单位砂地比曲线忽略了细节的变化,更加宏观的展现水体进退趋势,能够直观的观察层序变化。因此,利用本申请实施例提供的储层层序识别的方法,可以提高层序识别的准确性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为某地区流一2亚段砂地比的平面分布图;
图2为某地区流一3亚段砂地比的平面分布图;
图3是本申请储层层序识别方法一个实施例的流程图;
图4是本申请实施方式中单井层序分析时y3-1x井单井岩性、砂地比曲线和层序图;
图5是本申请实施方式中连井层序分析时y4x、y3-1x、y16x和y1x四井的单井岩性、砂地比曲线和层序图;
图6是本申请实施方式中纵波阻抗反演剖面的一个示意图
图7是本申请实施方式中1米单位砂地比曲线拟声波反演剖面示意图;
图8是本申请实施方式中10米单位砂地比曲线拟声波反演剖面示意图;
图9是本申请实施方式中25米单位砂地比曲线拟声波反演剖面示意图;
图10是本申请储层层序识别装置一个实施例的模块图;
图11是本申请装置实施例中数字化处理单元的子模块图;
图12是本申请装置实施例中砂地比曲线确定单元的子模块图。
具体实施方式
本申请实施例提供一种储层层序识别方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图3是本申请储层层序识别的方法一个实施例的流程图。参照图3,所述方法可以包括以下步骤。
S101:获取岩性信息。
所述岩性信息可以包括:单井砂岩分布深度。
在一个实施方式中,可以从岩性录井资料获取所述单井砂岩分布深度。
在另一实施方式中,可以根据现有的测井曲线获取所述单井砂岩分布深度。例如,可以利用砂岩特征识别伽马测井曲线中的砂岩部分,并获取所述伽马测井曲线中砂岩部分对应的深度。
S102:对所述岩性信息进行数字化处理。
可以对所述获取的岩性信息进行数字化处理。
在一个实施方式中,所述对岩性信息进行数字化处理可以包括:对所述岩性信息按照预设间隔进行采样,获取采样数据;判断所述采样数据是否符合预设条件;根据所述判断结果对所述采样数据进行标识。所述预设条件可以为:该深度处的岩石为砂岩或砾岩。
具体地,可以将所述符合预设条件的采样数据标识为第一标识,可以将不符合预设条件的采样数据标识为第二标识。
进一步地,所述第一标识可以为1。所述第二标识可以为0。
在一个实施方式中,所述预设间隔可以为0.05米-0.5米。例如,可以取值为0.1米。通过调整预设间隔,可以得到合适的预设间隔,以保证采样得到的采样数据可以完整反映岩石类别。
在一个应用场景中,可以利用excel实现对所述岩性信息的数字化处理。
具体地,可以将所述岩性信息输入到名称为“Sheet1”的表格中,所述表格中可以按列展现下述数据:井名、开始深度、结束深度和岩性信息描述。可以将井名、开始深度、结束深度和岩性信息描述分别作为A、B、C、D四列。
可以将标题作为第一行,并在E列输入如下代码:“=IF(COUNTIF(D2,"*砂岩*")>0,1,IF(COUNTIF(D2,"*砾岩*")>0,1,0))”,利用上述代码,可以在E列得到数字化处理的结果,即,一行中岩性描述若为砂岩或者砾岩,则该行E列的值为1,若为其他岩性则该行E列的值为0。进一步地,还可以删除E列中为0的数据。
进一步地,可以在空白的工作表(例如“Sheet2”)的A1中输入最浅砂岩的起始深度,根据预设间隔填充A列中的其他单元格,直至所述岩性信息的最大深度。
在E2输入代码“=INT(B2*10)/10”,并在E列按照E2的公式按列填充。可以将开始深度保留到小数点后第一位。
在F5输入公式“=INT(C2*10)/10”,并在F列按照F2的公式按列填充。可以将结束深度保留到小数点后第一位。
可以在G5、G6分别输入代码“=E2”、“=F2”,并将G5、G6按公式向下填充,这样在G列中G2为第一层开始深度,G3为第一层结束深度,依次类推。
可以在单元格H2中输入1,单元格H3中输入代码“=IF(F3=F2,H2+1,1)”,表示如果F列上下行数据相同,说明该行G列数据为本层岩性数据结束深度,其值可以为2,否则为起始深度,其值可以为1,G列的“1”值称为起止标识。
可以将Sheet1中G、H列中数值复制到Sheet2中。
可以Sheet2中将A列数值升序排列,并确定扩展选定区域。这样A列中可以包含所有砂岩起止深度中间的所有采样点。每层砂岩开始、终止深度都有两个采样点,这两个采样点所在行B列中数据分别为1、0或者2、0。
可以在Sheet2中的第一列前插入一列,插入后,可以在A2中输入代码:“=IF(ISBLANK(C2),IF(C1>0,1,0),IF(B1=B2,IF(C2=2,1,0),0))”,上述代码可以用于标识重复行。
可以将Sheet2中A列数据为1的删除。
可以将Sheet2中A列数据清空,并在第1行输入1,第二行输入代码“=IF(C2=1,1,IF(C1=2,0,A1))”,上述代码可以使A列砂岩或砾岩层开始和结束深度间的单元格值为1,其余单元格值为0。该列数据可以作为为砂岩标识,表示该采样点如果在砂岩或砾岩层内值为1,否则值为0。
S103:根据所述数字化处理的岩性信息,确定单位砂地比曲线。
根据所述数字化处理的岩性信息,可以确定单位砂地比曲线。所述单位砂地比可以表示以采样点为中心的地层厚度范围内砂层与地层厚度的比值。通常,地层厚度为1米时,可以将计算得到的单位砂地比称为1米单位砂地比;地层厚度为10米时,可以将计算得到的单位砂地比称为10米单位砂地比;地层厚度为25米时,可以将计算得到的单位砂地比称为25米单位砂地比。
在一个实施方式中,所述根据数字化处理的岩性信息,确定单位砂地比曲线,具体可以包括:以一个采样点为中心,统计所述地层厚度范围内第一标识的个数和第二标识的个数;计算第一标识和第二标识的总个数;将所述第一标识的个数与所述总个数的比值作为该采样点的砂地比;将所述数字处理后的岩性信息对应的采样点的砂地比连接为曲线,即为单位砂地比曲线。
在上述利用excel的应用场景中,可以在Sheet2中将砂岩标识所在列剪切到F列,并删除B列中的起止标识。进一步地,可以插入250行空白行。
可以在单元格B251中输入代码“=SUM(F246:F256)/11”,上述代码可以用于计算11个采样点中砂岩或者砾岩采样点的个数与总采样点个数的比值,该列数据即为1米单位砂地比曲线(例如,附图4中的1Sratio)。
可以在单元格C251中输入代码“=SUM(F201:F301)/101”,上述代码可以用于计算101个采样点中砂岩或者砾岩采样点的个数与总采样点个数的比值,该列数据即为10米单位砂地比曲线(例如,附图4中的10Sratio)。
可以在单元格D251中输入代码“=SUM(F126:F376)/251”,上述代码可以用于计算251个采样点中砂岩或者砾岩采样点的个数与总采样点个数的比值,该列数据即为25米单位砂地比曲线(例如,附图4中的25Sratio)。
可以在单元格E251中输入代码“=SUM(F1:F501)/501”,上述代码可以用于计算501个采样点中砂岩或者砾岩采样点的个数与总采样点个数的比值,该列数据即为50米单位砂地比曲线(例如,附图4中的50Sratio)。
S104:对所述单位砂地比曲线进行第一处理,并根据所述第一处理结果进行层序识别。
可以对所述单位砂地比曲线进行第一处理,并根据所述第一处理结果进行层序识别。
在一个实施方式中,所述第一处理可以包括下述中的至少一项:单井分析、连井对比、反演。
在一个实施方式中,在进行单井分析时,单位砂地比曲线却对其内部砂体富集程度的变化有直观的显示。具体地,1米单位砂地比曲线能够精细表述当前深度是否为砂层;10米单位砂地比曲线能够较好的反映砂层组信息;25米单位砂地比曲线能够较好体现砂体变化趋势;50米单位砂地比曲线忽略了微观的砂体变化,能更直观观察到砂体富集趋势的整体变化。
在一个实施方式中,在进行连井对比时,可以分析不同采样间隔对应的砂地比曲线。通常,相邻井位在同一时期沉积的地磁往往具有相似曲线特征,可以根据相似曲线特征进行层序识别。
在一个实施方式中,在进行反演时,可以将所述单位砂地比曲线作为纵波阻抗进行输入进行反演。
在一个实施方式中,进行层序分析时,可以用向上的三角表示地层沉积的过程中水面上升、砂岩粒度逐渐变细、砂岩减少、砂层变薄、砂体向近物源方向移动等属于水进体系域的过程;可以用向下的三角表示地层沉积的过程中水面下降、砂岩粒度逐渐变粗、砂岩增多、砂层变厚、砂体向远物源方向盆地中心移动等属于高位域或低位域的过程。
通常,单井层序分析时:1米单位砂地比曲线的大部分深度处曲线值为0或1,可以直接观察砂岩的具体分布;10米单位砂地比曲线中,一些砂岩比较密集的层段曲线值接近1,中间的薄泥岩层也不会突然降到0,可以从曲线中划分出高值区域和低值区域,并根据曲线峰值的变化,可以判断在某一层段内的水进水退趋势,划分四、五级层序;25米单位砂地比曲线中,曲线的最小值一般不为0,最高值一般不为1,根据砂地比曲线能够分辨出增大或者减小的趋势,可以根据曲线峰值和趋势的变化,判断在某一层段内的水进水退趋势,划分三、四级层序;50米单位砂地比曲线中,砂地比曲线值的增大或减小的趋势更加明显,根据所述砂地比曲线趋势的变化,可以判断在较长时期某一层段内的水进水退趋势,划分三四级层序。参照图4,图4为单井层序分析时y3-1x井单井岩性、砂地比曲线和层序图。
连井层序分析时:可以根据多口井的单位砂地比曲线所反映的砂体分布特征和砂地比变化趋势划分层序。参照图5,图5为连井层序分析时,y4x、y3-1x、y16x和y1x四井的单井岩性、砂地比曲线和层序图。
分析反演结果时,可以看到常规反演剖面同相轴往往体现的是砂泥岩的相对变化,图6是纵波阻抗反演剖面的一个示意图。图6中,井轨迹左侧为伽马曲线,右侧为电阻率曲线。从图6可见,中下部的砂体并不比上部更富集,强轴同相轴峰值幅度却比上部高。单位砂地比反演剖面同相轴反映的是砂体的富集程度,同相轴的强弱代表的是与砂岩在某一位置的富集程度是一一对应的。图7是1米单位砂地比曲线拟声波反演剖面示意图。参照图7,1米单位砂地比反演结果能分辨的砂体更薄、砂层更多,可利用该反演结果能进行精细的层序分析。图8是10米单位砂地比曲线拟声波反演剖面示意图。参照图8,10米单位砂地比反演剖面中,一些离的较近的砂层连成一片,但是内部仍可区分相对的强弱关系,在整个剖面上各区域的相对强弱关系明显,上下砂体富集变化趋势明显,可用于四级层序地层分析。图9是25米单位砂地比曲线拟声波反演剖面示意图。参照图9,25米单位砂地比反演剖面中,密集分布的砂体合并成幅度强、频率低的同相轴,其间分布较多泥岩的砂地层显示为强相位、高频率的较细的同相轴,泥岩密集分布的地层显示为负相位同相轴。由图9可以看出,25米单位砂地比反演剖面相对10米单位砂地比反演剖面,原来较难对比的两套强轴之间也有了相对的强弱变化,这样就能在更大的尺度上进行层序分析,可用于三、四级层序地层分析。
利用本申请实施例提供的层序识别方法,对于小尺度(如1米)的单位砂地比曲线能够反映砂岩与非砂岩情况,储层反演中同相轴清晰,能够区分较薄的砂层;中小尺度(如10米)的单位砂地比曲线,能够明显反映出砂层组合特征,便于连井层序对比;中大尺度(如25米)的单位砂地比曲线,能够直观地体现水进水退特征,便于进行精细的层序对比;大尺度(如50米)的单位砂地比曲线忽略了细节的变化,更加宏观的展现水体进退趋势,能够直观的观察层序变化。因此,利用本申请实施例提供的储层层序识别的方法,可以提高层序识别的准确性。
本申请还提供一种储层层序识别装置的实施例。
图10是本申请储层层序识别装置一个实施例的模块图。参照图10,所述储层层序识别装置可以包括:信息获取单元100、数字化处理单元200、砂地比曲线确定单元300和识别单元400。
所述信息获取单元100,可以用于获取岩性信息.所述岩性信息可以包括:单井砂岩分布深度。
所述数字化处理单元200,可以用于对所述岩性信息进行数字化处理。
所述砂地比曲线确定单元300,可以用于根据所述数字化处理的岩性信息,确定单位砂地比曲线。
所述识别单元400,可以用于对所述单位砂地比曲线进行第一处理,并根据所述第一处理结果进行层序识别。
图11是本申请装置实施例中数字化处理单元的子模块图。参照图11,所述数字化处理单元200可以包括:采样子单元201、判断子单元202和标识子单元203。
所述采样子单元201,可以用于对所述岩性信息按照预设间隔进行采样,获取采样数据。
所述判断子单元202,可以用于判断所述采样数据是否符合预设条件。所述预设条件为:该深度处的岩石为砂岩或砾岩。
所述标识子单元203,可以用于根据所述判断结果对所述采样数据进行标识。
图12是本申请装置实施例中砂地比曲线确定单元的子模块图。参照图12,所述砂地比曲线确定单元300可以包括:统计子单元301、计算子单元302、比值子单元303和连接子单元304。
所述统计子单元301,可以用于以一个采样点为中心,统计所述地层厚度范围内第一标识的个数和第二标识的个数。
所述计算子单元302,可以用于计算第一标识和第二标识的总个数。
所述比值子单元303,可以用于将所述第一标识的个数与所述总个数的比值作为该采样点的砂地比。
所述连接子单元304,可以用于将所述数字处理后的岩性信息对应的采样点的砂地比连接为曲线。
上述实施例公开的储层层序识别装置实施例与本申请实施例公开的储层层序识别方法相对应,可以实现本申请方法实施例,并取得方法实施例的技术效果。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的装置、模块,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。
为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中,内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括短暂电脑可读媒体(transitory media),如调制的数据信号和载波。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (9)
1.一种储层层序识别方法,其特征在于,包括:
获取岩性信息;
对所述岩性信息进行数字化处理;
根据所述数字化处理的岩性信息,确定单位砂地比曲线,所述单位砂地比表示以采样点为中心的地层厚度范围内砂层与地层厚度的比值;
对所述单位砂地比曲线进行第一处理,并根据所述第一处理结果进行层序识别,所述第一处理包括下述中的至少一项:单井分析、连井对比、反演;
其中,所述根据所述数字化处理的岩性信息,确定单位砂地比曲线,包括:
以一个采样点为中心,统计地层厚度范围内第一标识的个数和第二标识的个数;计算第一标识和第二标识的总个数;将所述第一标识的个数与所述总个数的比值作为该采样点的砂地比;将所述数字处理后的岩性信息对应的采样点的砂地比连接为曲线,获得单位砂地比曲线。
2.根据权利要求1所述的方法,其特征在于,所述岩性信息包括:单井砂岩分布深度。
3.根据权利要求2所述的方法,其特征在于,所述获取岩性信息包括:从岩性录井资料获取所述单井砂岩分布深度,或者,根据现有的测井曲线获取所述单井砂岩分布深度。
4.根据权利要求1所述的方法,其特征在于,所述对岩性信息进行数字化处理包括:
对所述岩性信息按照预设间隔进行采样,获取采样数据;
判断所述采样数据是否符合预设条件;
根据所述判断结果对所述采样数据进行标识。
5.根据权利要求4所述的方法,其特征在于,所述预设条件为:该采样数据对应深度处的岩石为砂岩或砾岩。
6.根据权利要求4所述的方法,其特征在于,所述对所述采样数据进行标识包括:将所述符合预设条件的采样数据标识为第一标识,将不符合预设条件的采样数据标识为第二标识。
7.根据权利要求6所述的方法,其特征在于,所述第一标识为1;所述第二标识为0。
8.一种储层层序识别装置,其特征在于,包括:信息获取单元、数字化处理单元、砂地比曲线确定单元和识别单元;
所述信息获取单元,用于获取岩性信息;所述岩性信息包括:单井砂岩分布深度;
所述数字化处理单元,用于对所述岩性信息进行数字化处理;
所述砂地比曲线确定单元,用于根据所述数字化处理的岩性信息,确定单位砂地比曲线,所述单位砂地比表示以采样点为中心的地层厚度范围内砂层与地层厚度的比值;
所述识别单元,用于对所述单位砂地比曲线进行第一处理,并根据所述第一处理结果进行层序识别,所述第一处理包括下述中的至少一项:单井分析、连井对比、反演;
其中,所述砂地比曲线确定单元包括:统计子单元、计算子单元、比值子单元和连接子单元;
所述统计子单元,用于以一个采样点为中心,统计地层厚度范围内第一标识的个数和第二标识的个数;
所述计算子单元,用于计算第一标识和第二标识的总个数;
所述比值子单元,用于将所述第一标识的个数与所述总个数的比值作为该采样点的砂地比;
所述连接子单元,用于将所述数字处理后的岩性信息对应的采样点的砂地比连接为曲线,获得单位砂地比曲线。
9.根据权利要求8所述的装置,其特征在于,所述数字化处理单元包括:采样子单元、判断子单元和标识子单元;
所述采样子单元,用于对所述岩性信息按照预设间隔进行采样,获取采样数据;
所述判断子单元,用于判断所述采样数据是否符合预设条件;所述预设条件为:该深度处的岩石为砂岩或砾岩;
所述标识子单元,用于根据所述判断结果对所述采样数据进行标识。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810001664.8A CN108254786B (zh) | 2018-01-02 | 2018-01-02 | 一种储层层序识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810001664.8A CN108254786B (zh) | 2018-01-02 | 2018-01-02 | 一种储层层序识别方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108254786A CN108254786A (zh) | 2018-07-06 |
CN108254786B true CN108254786B (zh) | 2019-09-10 |
Family
ID=62725433
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810001664.8A Active CN108254786B (zh) | 2018-01-02 | 2018-01-02 | 一种储层层序识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108254786B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109298465B (zh) * | 2018-10-23 | 2020-04-10 | 中国石油天然气股份有限公司 | 相控砂体结构主体块状砂岩定量区分方法及装置 |
CN112649887B (zh) * | 2019-10-11 | 2024-03-01 | 中国石油天然气股份有限公司 | 基于钻井资料定量划分层序地层的方法及装置 |
CN111058829B (zh) * | 2019-12-05 | 2021-06-25 | 中国矿业大学 | 基于图像处理的岩层分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105844082A (zh) * | 2016-03-11 | 2016-08-10 | 西安石油大学 | 一种沉积相图中主导相的判断方法及系统 |
CN107272081A (zh) * | 2017-07-18 | 2017-10-20 | 中国石油大学(华东) | 一种山前地区沉积相带展布范围预测方法 |
-
2018
- 2018-01-02 CN CN201810001664.8A patent/CN108254786B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105844082A (zh) * | 2016-03-11 | 2016-08-10 | 西安石油大学 | 一种沉积相图中主导相的判断方法及系统 |
CN107272081A (zh) * | 2017-07-18 | 2017-10-20 | 中国石油大学(华东) | 一种山前地区沉积相带展布范围预测方法 |
Non-Patent Citations (2)
Title |
---|
Sql Server 数据库技术实现海量录井地质数据的处理及应用;李斌 等;《石油工业计算机应用》;20160131(第1期);第7-12页 |
海量录井层序"砂地比"的自动计算;李斌 等;《高校地质学报》;20120331;第18卷(第1期);第189-192页 |
Also Published As
Publication number | Publication date |
---|---|
CN108254786A (zh) | 2018-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Neal et al. | Accommodation succession (δ A/δ S) sequence stratigraphy: observational method, utility and insights into sequence boundary formation | |
Steel et al. | Clinoforms, clinoform trajectories and deepwater sands | |
CN109061765B (zh) | 非均质薄砂岩互层油藏的圈闭评价方法 | |
CN108254786B (zh) | 一种储层层序识别方法及装置 | |
International Subcommission on Stratigraphic Classification (Amos SALVADOR, Chairman) | Unconformity-bounded stratigraphic units: INTERNATIONAL SUBCOMMISSION ON STRATIGRAPHIC CLASSIFICATION (AMOS SALVADOR, Chairman) | |
Compston et al. | Further SHRIMP geochronology on the early Cambrian of South China | |
CN111273348B (zh) | 基于更新概率比率恒定理论的多点地质统计叠前反演方法 | |
Lia et al. | Uncertainties in reservoir production forecasts | |
CN107765298B (zh) | 一种确定近地表速度模型的方法及装置 | |
Gunnell et al. | Shoulder uplift of the Western Ghats passive margin, India: a denudational model | |
Cooper et al. | Quantitative biostratigraphy of the Taranaki Basin, New Zealand: A deterministic and probabilistic approach | |
CN107092031A (zh) | 一种声波测井数据校正方法及装置 | |
CN109389154B (zh) | 洪积扇砂砾岩储层岩相识别方法及装置 | |
Mitchell et al. | Erosion of the submarine flanks of the Canary Islands | |
CN106526671B (zh) | 一种碳酸盐岩储层地震波响应特征分析方法 | |
CN108490491A (zh) | 一种基于波形指示反演的滩体预测方法 | |
CN105240006A (zh) | 一种适用于火山岩储层的油水层识别方法 | |
CN108386185A (zh) | 一种确定砂体连通性的方法及装置 | |
CN112505754B (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
Castillo et al. | Prehistoric earthquakes on the banning strand of the san andreas fault, north palm springs, California | |
Popova | Analogy in the World of Geological Uncertainties, or How Reservoir Analogs May Refine Your Probabilistic Geomodel | |
CN105089659A (zh) | 一种砾岩油藏渗流单元识别方法 | |
Clark et al. | Eocene turbidite-population statistics from shelf edge to basin floor, Spitsbergen, Svalbard | |
CN106772596A (zh) | 一种确定叠前时间偏移速度场的方法及装置 | |
CN107765318B (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 |