CN109113732A - 储层非均质性的确定方法及装置 - Google Patents
储层非均质性的确定方法及装置 Download PDFInfo
- Publication number
- CN109113732A CN109113732A CN201810900739.6A CN201810900739A CN109113732A CN 109113732 A CN109113732 A CN 109113732A CN 201810900739 A CN201810900739 A CN 201810900739A CN 109113732 A CN109113732 A CN 109113732A
- Authority
- CN
- China
- Prior art keywords
- permeability
- reservoir
- median
- heterogeneity
- rectangular box
- 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
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供了一种储层非均质性的确定方法及装置,其中,该方法包括:在储层图上标识所有样品的位置;根据每一样品的渗透率,绘制所有样品渗透率的数据盒须图;数据盒须图用来表征渗透率的分布范围;数据盒须图的上须为所有样品渗透率中的最大值,下须为所有样品渗透率中的最小值;数据盒须图包括代表渗透率集中分布的矩形盒;将每一样品所代表的储层与数据盒须图关联,根据关联后的储层图与数据盒须图,确定储层非均质性。上述技术方案为储层非均质性的可视化确定方案,能直观、准确地反映储层的非均质性,可以预警油气田开发会出现的问题,为油田开发层系划分及开发方案设计提供依据,有利于高效开发油藏。
Description
技术领域
本发明涉及油气田开发技术领域,特别涉及一种储层非均质性的确定方法及装置。
背景技术
储层非均质性是指油气储层在漫长的地质历史中,经历了沉积、成岩及后期构造作用的综合影响,使储层的空间分布及内部的各种属性都存在极不均匀的变化。
目前,有一些学者利用成岩岩相分析法对轮南东斜坡东河砂岩非均质性进行了研究;有一些学者从宏观和微观两方面分析了东营凹陷古近系沙河街组碎屑岩储集层的非均质性特征;有一些学者以松辽盆地朝阳沟油田朝1-朝气3区块扶余油层为例,对高分辨率层序地层学,在储层宏观非均质性研究中的应用进行了研究;有一些学者对储层非均质性描述的地质统计学方法进行了研究,认为该方法具有一定的优越性和实用性;有一些学者利用变差函数方法,来分析储集层平面非均质性,研究中通过构建一种定量表征储集层平面非均质性的数学模型,来实现储集层平面非均质性表征,结果表明该方法优于传统的数理统计方法;有一些学者以江苏省扬州市真武油田戴南组二段碎屑岩储层为研究对象,利用数据包络分析法表征碎屑岩储层非均质性。
另外,中国专利申请CN201710256713.8提供了一种确定储层非均质性的方法,在确定非均质性的过程中,依据的是待测储层各个深度组中所有的地层内容,而不仅是岩心样品;中国专利申请CN201410675357.X提供了一种根据岩心刻度电成像测井、储层的岩性和层理构造以及核磁共振测井资料对储层的微观非均质性进行分析的方法。
然而,基于地质成因和地层学的方法不能定量表征储层非均质性,目前基于数学方法定量表征储层非均质性的方法又不够直观准确,往往与实际油藏开发相脱节。
发明内容
本发明实施例提供了一种储层非均质性的确定方法,用以直观、准确地反映储层的非均质性,该方法包括:
在储层图上标识所有样品的位置;
根据每一样品的渗透率,绘制所有样品渗透率的数据盒须图;所述数据盒须图用来表征渗透率的分布范围;所述数据盒须图的上须为所有样品渗透率中的最大值,下须为所有样品渗透率中的最小值;所述数据盒须图包括代表渗透率集中分布的矩形盒;
将每一样品所代表的储层与数据盒须图关联,根据关联后的储层图与数据盒须图,确定储层非均质性。
本发明实施例还提供了一种储层非均质性的确定装置,用以直观、准确地反映储层的非均质性,该装置包括:
标识单元,用于在储层图上标识所有样品的位置;
绘制单元,用于根据每一样品的渗透率,绘制所有样品渗透率的数据盒须图;所述数据盒须图用来表征渗透率的分布范围;所述数据盒须图的上须为所有样品渗透率中的最大值,下须为所有样品渗透率中的最小值;所述数据盒须图包括代表渗透率集中分布的矩形盒;
确定单元,用于将每一样品所代表的储层与数据盒须图关联,根据关联后的储层图与数据盒须图,确定储层非均质性。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述储层非均质性的确定方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行储层非均质性的确定方法的计算机程序。
本发明实施例提供的技术方案。上述技术方案实现了储层非均质性的数据可视化,能直观、准确地反映储层的非均质性,可以预警油气田开发会出现的问题,为油田开发层系划分及开发方案设计提供依据,有利于高效开发油藏。
本发明实施例提供的技术方案将多样的渗透率数据,通过绘图的方式展现,直观的反映储层渗透率在储层的展布,通过特殊的绘图方式能更高便捷的分析储层渗透率的非均质性,预警油气田开发会出现的问题,为油田开发方案设计提供依据:首先标识样品,其次绘制渗透率数据盒须图,通过盒须图中的矩形盒找到储层渗透率相对集中的范围,进而确定储层非均质程度。在本发明实施例中,油层的渗透率范围通过盒须图的上须、下须一目了然,矩形盒可以直观地展现储层渗透率的集中分布范围,通过矩形盒所展现的渗透率集中范围最终确定的非均质性,去除了反映储层非均质性次要矛盾的数据点,使结果更加准确,同时,该数据可视化确定储层非均质性的方法也为油田开发层系划分及开发方案设计提供依据,有利于高效开发油藏。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是本发明实施例中储层非均质性的确定方法的流程示意图;
图2是本发明实施例中储层非均质性的确定装置的结构示意图;
图3是本发明实施例中目标试验区纵向非均质性数据可视化的评价示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
发明人发现:现有储层非均质性研究的不足:(1)大部分基于地质因素,不能定量表征非均质性,很难对比评价储层;(2)很多应用数学计算方法确定储层非均质性的方案并不能准确表明储层非均质性;(3)应用数学方法定量表征非均质性的方案考虑了所有获取数据,并没有抓住主要矛盾,不能为油藏开发提供实用性依据。
因此,考虑到上述技术问题,发明人提出了一种储层非均质性的确定方案,即储层非均质性数据可视化的评价方案,该方案通过统计学及特殊绘图方法,解决现有技术储层非均质性评价中存在的问题,结果更加直观、准确,该方案主要包括:首先标识样品,其次绘制渗透率数据盒须图,通过盒须图中的矩形盒找到储层渗透率相对集中的范围,计算两个矩形盒内的加权平均渗透率,最终得到非均质性准数,确定储层非均质程度。在本发明实施例中,油层的渗透率范围通过盒须图的上须、下须一目了然,中位数(第一中位数)两侧的矩形盒可以直观的展现储层渗透率的集中分布范围,通过矩形盒所展现的渗透率集中范围最终确定的非均质性,去除了反映储层非均质性次要矛盾的数据点,使结果更加准确,同时,该数据可视化方法也为油田开发层系划分及开发方案设计提供依据,有利于高效开发油藏。下面对该储层非均质性数据可视化的确定方案进行详细介绍如下。
图1是本发明实施例中储层非均质性的确定方法的流程示意图,如图1所示,该方法包括如下步骤:
步骤101:在储层图上标识所有样品的位置;
步骤102:根据每一样品的渗透率,绘制所有样品渗透率的数据盒须图;所述数据盒须图用来表征渗透率的分布范围;所述数据盒须图的上须为所有样品渗透率中的最大值,下须为所有样品渗透率中的最小值;所述数据盒须图包括代表渗透率集中分布的矩形盒;
步骤103:将每一样品所代表的储层与数据盒须图关联,根据关联后的储层图与数据盒须图,确定储层非均质性。
本发明实施例提供的技术方案将多样的渗透率数据,通过绘图的方式展现,直观的反映储层渗透率在储层的展布,通过特殊的绘图方式能更高便捷的分析储层渗透率的非均质性,预警油气田开发会出现的问题,为油田开发方案设计提供依据:首先标识样品,其次绘制渗透率数据盒须图,通过盒须图中的矩形盒找到储层渗透率相对集中的范围,进而确定储层非均质程度。在本发明实施例中,油层的渗透率范围通过盒须图的上须、下须一目了然,矩形盒可以直观地展现储层渗透率的集中分布范围,通过矩形盒所展现的渗透率集中范围最终确定的非均质性,去除了反映储层非均质性次要矛盾的数据点,使结果更加准确,同时,该数据可视化确定储层非均质性的方法也为油田开发层系划分及开发方案设计提供依据,有利于高效开发油藏。
具体实施时,首先,可以采集代表每一储层的样品,然后,预先测量代表每一储层的样品的渗透率,该样品的渗透率表征其所代表储层的渗透率,储层渗透率可以用来研究储层的非均质性;再次,将样品与储层图上对应储层位置(标识)相关联,即在储层图上标识所有样品的位置。
首先,介绍上述步骤101。
在一个实施例中,在储层图上标识所有样品的位置,可以包括:
若是分析平面渗透率非均质性,在储层平面图上标识所有样品的位置;
若是分析纵向渗透率非均质性,在储层纵向图上标识所有样品的位置。
具体实施时,上述步骤101在储层图上标识所有样品的位置的过程就是找到每一样品所代表的储层的过程,这样就知道了每一储层的渗透率。若是分析平面渗透率非均质性,则将样品点在储层平面图上标识;若是分析纵向渗透率非均质性,则将样品点在储层纵向图上标识。
其次,介绍上述步骤102的绘制渗透率数据盒须图的过程。
在一个实施例中,根据每一样品的渗透率,绘制所有样品渗透率的数据盒须图,可以包括:
将所有样品渗透率按升序排序;
将样品渗透率最大值作为数据盒须图的上须,将样品渗透率最小值作为数据盒须图的下须;绘制所述上须和下须之间的对数刻度线,将所有样品渗透率对应在所述对数刻度线上;
计算所有样品渗透率的第一中位数;
计算所述下须与所述第一中位数之间渗透率的第二中位数;
根据下须、第一中位数以及所述第二中位数之间的关系,绘制代表渗透率集中分布的第一矩形盒;
计算所述上须与所述第一中位数之间渗透率的第三中位数;
根据上须、第一中位数以及所述第三中位数之间的关系,绘制代表渗透率集中分布的第二矩形盒。
具体实施时,本发明实施例中数据盒须图的上须代表渗透率最大阈值,下须代表渗透率最小阈值。通过上、下须可以直观的展现储层渗透率的分布范围。
在一个实施例中,根据下须、第一中位数以及所述第二中位数之间的关系,绘制代表渗透率集中分布的第一矩形盒,可以包括:
若下须与第二中位数相逼近,将下须与第二中位数作为两端边绘制一个矩形盒作为第一矩形盒;
若第一中位数与第二中位数相逼近,将第一中位数与第二中位数作为两端边绘制一个矩形盒作为第一矩形盒;
根据上须、第一中位数以及所述第三中位数之间的关系,绘制代表渗透率集中分布的第二矩形盒,可以包括:
若上须与第三中位数相逼近,则将上须与第三中位数作为两端边绘制一个矩形盒作为第二矩形盒;
若第一中位数与第三中位数相逼近,则将第一中位数与第三中位数作为两端边绘制一个矩形盒作为第二矩形盒。
具体实施时,绘制渗透率数据盒须图的过程可以包括:
①将所有样品渗透率由小到大排序,上须XN为渗透率最大值,下须X1为渗透率最小值,上须和下须之间用对数刻度线,用来表征储层渗透率的分布范围:
X1,……,XN。
②计算所有渗透率的中位数A(第一中位数):
当N为奇数时,A=X(N+1)/2;
当N为偶数时,A=(XN/2+XN/2+1)/2。
③计算下须X1与A之间渗透率的中位数a1(第二中位数)作为上枢纽,若下须X1与a1相逼近,则将下须X1与a1作为两端边绘制一个矩形盒作为第一矩形盒;若A与a1相逼近,则将A与a1作为两端边绘制一个矩形盒作为第一矩形盒,计算第二中位数a1的过程可以为:
当N和(N+1)/2都为奇数时,a1=X(N+1)/4;
当N为偶数,N/2为奇数时,a1=X(N+2)/4;
当N为偶数,N/2为偶数时,a1=(XN/4+XN/4+1)/2。
④计算上须XN与A之间渗透率的中位数a2(第三中位数)作为下枢纽,若上须XN与a2相逼近,则将上须XN与a2作为两端边绘制一个矩形盒作为第二矩形盒;若A与a2相逼近,则将A与a2作为两端边绘制一个矩形盒作为第二矩形盒,第三中位数a2的计算过程可以为:
当N和(N+1)/2都为奇数时,a2=X3(N+1)/4;
当N为偶数,N/2为奇数时,a2=X(3N+2)/4;
当N为偶数,N/2为偶数时,a2=(X3N/4+X3N/4+1)/2。
矩形盒内的渗透率即储层渗透率的相对集中范围。具体实施时,通过第一中位数A两侧的矩形盒(第一矩形盒和第二矩形盒)可以直观地展现储层渗透率的集中分布范围。
在一个实施例中,上述储层非均质性的确定方法,还可以包括:
根据第一矩形盒内的样品渗透率和第二矩形盒内的样品渗透率,以及第一矩形盒内样品对应储层的有效厚度和第二矩形盒内样品对应储层的有效厚度,计算储层的非均质性准数;将每一样品所代表的储层与数据盒须图关联,根据关联后的储层图与数据盒须图,确定储层非均质性,包括:
根据所述非均质性准数,关联后的储层图与数据盒须图,确定储层非均质性。
具体实施时,通过盒须图中的矩形盒找到储层渗透率相对集中的范围,计算两个矩形盒内的加权平均渗透率,最终得到非均质性准数,确定储层非均质程度,通过非均质性准数能更准确的定量评价储层渗透率,结合非均质性准数确定储层非均质性,进一步提高了储层非均质性的确定的准确性。计算非均质性准数的过程可以为:
假设两个矩形盒的渗透率为Ki,Kj,其所对应油层的有效厚度hi,hj,计算非均质性准数C并将其标注:
C=[∑(Kihi)/∑hi]/[∑(Kjhj)/∑hj]。
具体实施时,本发明实施例中,“相逼近”的含义是两数值接近,例如:下须X1与a1相逼近的含义即为X1的值接近a1的值,详见下述具体实施例的介绍。
最后介绍上述步骤103:将样品标识点与盒须图刻度相关联,完成储层非均质性数据可视化评价。
基于同一发明构思,本发明实施例中还提供了一种储层非均质性的确定装置,如下面的实施例。由于储层非均质性的确定装置解决问题的原理与储层非均质性的确定方法相似,因此储层非均质性的确定装置的实施可以参考上述储层非均质性的确定方法的实施,重复之处不再赘述。以下所使用的,术语“模块”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图2是本发明实施例中储层非均质性的确定装置的结构示意图,如图2所示,该装置包括:
标识单元02,用于在储层图上标识所有样品的位置;
绘制单元04,用于根据每一样品的渗透率,绘制所有样品渗透率的数据盒须图;所述数据盒须图用来表征渗透率的分布范围;所述数据盒须图的上须为所有样品渗透率中的最大值,下须为所有样品渗透率中的最小值;所述数据盒须图包括代表渗透率集中分布的矩形盒;
确定单元06,用于将每一样品所代表的储层与数据盒须图关联,根据关联后的储层图与数据盒须图,确定储层非均质性。
在一个实施例中,所述绘制单元可以包括:
排序单元,用于将所有样品渗透率按升序排序;
盒须图绘制单元,用于将样品渗透率最大值作为数据盒须图的上须,将样品渗透率最小值作为数据盒须图的下须;绘制所述上须和下须之间的对数刻度线,将所有样品渗透率对应在所述对数刻度线上;
第一计算单元,用于计算所有样品渗透率的第一中位数;
第二计算单元,用于计算所述下须与所述第一中位数之间渗透率的第二中位数;
第一矩形盒绘制单元,用于根据下须、第一中位数以及所述第二中位数之间的关系,绘制代表渗透率集中分布的第一矩形盒;
第三计算单元,用于计算所述上须与所述第一中位数之间渗透率的第三中位数;
第二矩形盒绘制单元,用于根据上须、第一中位数以及所述第三中位数之间的关系,绘制代表渗透率集中分布的第二矩形盒。
在一个实施例中,所述第一矩形盒绘制单元具体可以用于:
若下须与第二中位数相逼近,将下须与第二中位数作为两端边绘制一个矩形盒作为第一矩形盒;
若第一中位数与第二中位数相逼近,将第一中位数与第二中位数作为两端边绘制一个矩形盒作为第一矩形盒;
所述第二矩形盒绘制单元具体可以用于:
若上须与第三中位数相逼近,则将上须与第三中位数作为两端边绘制一个矩形盒作为第二矩形盒;
若第一中位数与第三中位数相逼近,则将第一中位数与第三中位数作为两端边绘制一个矩形盒作为第二矩形盒。
在一个实施例中,上述储层非均质性的确定装置还可以包括:
非均质性准数计算单元,用于根据第一矩形盒内的样品渗透率和第二矩形盒内的样品渗透率,以及第一矩形盒内样品对应储层的有效厚度和第二矩形盒内样品对应储层的有效厚度,计算储层的非均质性准数;
所述确定单元具体用于:根据所述非均质性准数,关联后的储层图与数据盒须图,确定储层非均质性。
在一个实施例中,所述标识单元具体可以用于:
若是分析平面渗透率非均质性,在储层平面图上标识所有样品的位置;
若是分析纵向渗透率非均质性,在储层纵向图上标识所有样品的位置。
下面再结合附图3以实例来进行说明,以便于理解如何实施本发明。
本实施例为研究储层非均质性数据可视化的评价,下面结合附图3进行说明。新疆油田某储层纵向油层渗透率及有效厚度如下表1所示。
油层 | S<sub>1</sub><sup>1-2</sup> | S<sub>1</sub><sup>2-1</sup> | S<sub>1</sub><sup>2-2</sup> | S<sub>2</sub><sup>2</sup> | S<sub>3</sub><sup>1-1</sup> | S<sub>3</sub><sup>1-2</sup> | S<sub>3</sub><sup>2-1</sup> | S<sub>4</sub><sup>1-1</sup> | S<sub>4</sub><sup>2-1</sup> | S<sub>4</sub><sup>2-2</sup> | S<sub>5</sub><sup>1-1</sup> | S<sub>5</sub><sup>1-2</sup> | S<sub>5</sub><sup>2-1</sup> | S<sub>5</sub><sup>2-2</sup> |
有效厚度(<sub>m</sub>) | 2.07 | 1.77 | 2.51 | 1.31 | 1.96 | 1.96 | 2.09 | 2.56 | 1.65 | 1.79 | 4.22 | 3.59 | 2.59 | 3.23 |
渗透率(<sub>m</sub>D) | 27 | 105.1 | 35.9 | 27.1 | 312.7 | 25.9 | 8.4 | 650 | 16.6 | 1.8 | 10.2 | 7.9 | 292 | 8.6 |
表1
(1)将渗透率由小到大排序如下表2,上须为650mD,下须为1.8mD,如图3。
油层 | S<sub>4</sub><sup>2-2</sup> | S<sub>5</sub><sup>1-2</sup> | S<sub>3</sub><sup>2-1</sup> | S<sub>5</sub><sup>2-2</sup> | S<sub>5</sub><sup>1-1</sup> | S<sub>4</sub><sup>2-1</sup> | S<sub>3</sub><sup>1-2</sup> | S<sub>1</sub><sup>1-2</sup> | S<sub>2</sub><sup>2</sup> | S<sub>1</sub><sup>2-2</sup> | S<sub>1</sub><sup>2-1</sup> | S<sub>5</sub><sup>2-1</sup> | S<sub>3</sub><sup>1-1</sup> | S<sub>4</sub><sup>1-1</sup> |
有效厚度(<sub>m</sub>) | 1.79 | 3.59 | 2.09 | 3.23 | 4.22 | 1.65 | 1.96 | 2.07 | 1.31 | 2.51 | 1.77 | 2.59 | 1.96 | 2.56 |
渗透率(<sub>m</sub>D) | 1.8 | 7.9 | 8.4 | 8.6 | 10.2 | 16.6 | 25.9 | 27 | 27.1 | 35.9 | 105.1 | 292 | 312.7 | 650 |
表2
(2)渗透率中位数A=(XN/2+XN/2+1)/2=(X7+X8)/2=26.45,如图3;
(3)a1=X(N+2)/4=X4=8.6,由于a1与下须相逼近,将下须X1与a1作为两端边绘制一个矩形盒,如图3;
(4)a2=X(3N+2)/4=X11=105.1,由于a2与A相逼近,将A与a2作为两端边绘制一个矩形盒,如图3;
(5)计算非均质性准数C=[∑(Kihi)/∑hi]/[∑(Kjhj)/∑hj]=6.7,并标注在图上,如图3;
(6)将油层与盒须图相关联,如图3,根据关联后的储层图与数据盒须图,确定储层非均质性。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述储层非均质性的确定方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行储层非均质性的确定方法的计算机程序。
本发明实施提供的技术方案的有益技术效果为:
1:将多样的渗透率数据,通过绘图的方式展现,直观的反映储层渗透率在储层的展布;
2:通过上、下须可以直观的展现储层渗透率的分布范围;
3:通过中位数两侧的矩形盒可以直观地展现储层渗透率的集中分布范围;
4:通过特殊的绘图方式能更高便捷地分析储层渗透率的非均质性,预警油气田开发会出现的问题,为油田开发方案设计提供依据;
5:通过非均质性准数能更准确地定量评价储层渗透率,进而更准确地确定储层的非均质性。
综上,本发明实施例提供的技术方案中,油层的渗透率范围通过盒须图的上、下须一目了然,第一中位数两侧的矩形盒可以直观地展现储层渗透率的集中分布范围,通过矩形盒所展现的渗透率集中范围最终确定的非均质性,去除了反映储层非均质性次要矛盾的数据点,使结果更加准确,同时,该数据可视化方法也为油田开发层系划分及开发方案设计提供依据,有利于高效开发油藏。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (12)
1.一种储层非均质性的确定方法,其特征在于,包括:
在储层图上标识所有样品的位置;
根据每一样品的渗透率,绘制所有样品渗透率的数据盒须图;所述数据盒须图用来表征渗透率的分布范围;所述数据盒须图的上须为所有样品渗透率中的最大值,下须为所有样品渗透率中的最小值;所述数据盒须图包括代表渗透率集中分布的矩形盒;
将每一样品所代表的储层与数据盒须图关联,根据关联后的储层图与数据盒须图,确定储层非均质性。
2.如权利要求1所述的储层非均质性的确定方法,其特征在于,根据每一样品的渗透率,绘制所有样品渗透率的数据盒须图,包括:
将所有样品渗透率按升序排序;
将样品渗透率最大值作为数据盒须图的上须,将样品渗透率最小值作为数据盒须图的下须;绘制所述上须和下须之间的对数刻度线,将所有样品渗透率对应在所述对数刻度线上;
计算所有样品渗透率的第一中位数;
计算所述下须与所述第一中位数之间渗透率的第二中位数;
根据下须、第一中位数以及所述第二中位数之间的关系,绘制代表渗透率集中分布的第一矩形盒;
计算所述上须与所述第一中位数之间渗透率的第三中位数;
根据上须、第一中位数以及所述第三中位数之间的关系,绘制代表渗透率集中分布的第二矩形盒。
3.如权利要求2所述的储层非均质性的确定方法,其特征在于,根据下须、第一中位数以及所述第二中位数之间的关系,绘制代表渗透率集中分布的第一矩形盒,包括:
若下须与第二中位数相逼近,将下须与第二中位数作为两端边绘制一个矩形盒作为第一矩形盒;
若第一中位数与第二中位数相逼近,将第一中位数与第二中位数作为两端边绘制一个矩形盒作为第一矩形盒;
根据上须、第一中位数以及所述第三中位数之间的关系,绘制代表渗透率集中分布的第二矩形盒,包括:
若上须与第三中位数相逼近,则将上须与第三中位数作为两端边绘制一个矩形盒作为第二矩形盒;
若第一中位数与第三中位数相逼近,则将第一中位数与第三中位数作为两端边绘制一个矩形盒作为第二矩形盒。
4.如权利要求2所述的储层非均质性的确定方法,其特征在于,还包括:
根据第一矩形盒内的样品渗透率和第二矩形盒内的样品渗透率,以及第一矩形盒内样品对应储层的有效厚度和第二矩形盒内样品对应储层的有效厚度,计算储层的非均质性准数;
将每一样品所代表的储层与数据盒须图关联,根据关联后的储层图与数据盒须图,确定储层非均质性,包括:
根据所述非均质性准数,关联后的储层图与数据盒须图,确定储层非均质性。
5.如权利要求1所述的储层非均质性的确定方法,其特征在于,在储层图上标识所有样品的位置,包括:
若是分析平面渗透率非均质性,在储层平面图上标识所有样品的位置;
若是分析纵向渗透率非均质性,在储层纵向图上标识所有样品的位置。
6.一种储层非均质性的确定装置,其特征在于,包括:
标识单元,用于在储层图上标识所有样品的位置;
绘制单元,用于根据每一样品的渗透率,绘制所有样品渗透率的数据盒须图;所述数据盒须图用来表征渗透率的分布范围;所述数据盒须图的上须为所有样品渗透率中的最大值,下须为所有样品渗透率中的最小值;所述数据盒须图包括代表渗透率集中分布的矩形盒;
确定单元,用于将每一样品所代表的储层与数据盒须图关联,根据关联后的储层图与数据盒须图,确定储层非均质性。
7.如权利要求6所述的储层非均质性的确定装置,其特征在于,所述绘制单元包括:
排序单元,用于将所有样品渗透率按升序排序;
盒须图绘制单元,用于将样品渗透率最大值作为数据盒须图的上须,将样品渗透率最小值作为数据盒须图的下须;绘制所述上须和下须之间的对数刻度线,将所有样品渗透率对应在所述对数刻度线上;
第一计算单元,用于计算所有样品渗透率的第一中位数;
第二计算单元,用于计算所述下须与所述第一中位数之间渗透率的第二中位数;
第一矩形盒绘制单元,用于根据下须、第一中位数以及所述第二中位数之间的关系,绘制代表渗透率集中分布的第一矩形盒;
第三计算单元,用于计算所述上须与所述第一中位数之间渗透率的第三中位数;
第二矩形盒绘制单元,用于根据上须、第一中位数以及所述第三中位数之间的关系,绘制代表渗透率集中分布的第二矩形盒。
8.如权利要求7所述的储层非均质性的确定装置,其特征在于,所述第一矩形盒绘制单元具体用于:
若下须与第二中位数相逼近,将下须与第二中位数作为两端边绘制一个矩形盒作为第一矩形盒;
若第一中位数与第二中位数相逼近,将第一中位数与第二中位数作为两端边绘制一个矩形盒作为第一矩形盒;
所述第二矩形盒绘制单元具体用于:
若上须与第三中位数相逼近,则将上须与第三中位数作为两端边绘制一个矩形盒作为第二矩形盒;
若第一中位数与第三中位数相逼近,则将第一中位数与第三中位数作为两端边绘制一个矩形盒作为第二矩形盒。
9.如权利要求7所述的储层非均质性的确定装置,其特征在于,还包括:
非均质性准数计算单元,用于根据第一矩形盒内的样品渗透率和第二矩形盒内的样品渗透率,以及第一矩形盒内样品对应储层的有效厚度和第二矩形盒内样品对应储层的有效厚度,计算储层的非均质性准数;
所述确定单元具体用于:根据所述非均质性准数,关联后的储层图与数据盒须图,确定储层非均质性。
10.如权利要求6所述的储层非均质性的确定装置,其特征在于,所述标识单元具体用于:
若是分析平面渗透率非均质性,在储层平面图上标识所有样品的位置;
若是分析纵向渗透率非均质性,在储层纵向图上标识所有样品的位置。
11.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至5任一所述方法。
12.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至5任一所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810900739.6A CN109113732B (zh) | 2018-08-09 | 2018-08-09 | 储层非均质性的确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810900739.6A CN109113732B (zh) | 2018-08-09 | 2018-08-09 | 储层非均质性的确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109113732A true CN109113732A (zh) | 2019-01-01 |
CN109113732B CN109113732B (zh) | 2022-03-29 |
Family
ID=64851625
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810900739.6A Active CN109113732B (zh) | 2018-08-09 | 2018-08-09 | 储层非均质性的确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109113732B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110500069A (zh) * | 2019-08-02 | 2019-11-26 | 中国矿业大学(北京) | 各向异性非均质储层煤层气开采完井技术选择方法 |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090260883A1 (en) * | 2008-04-16 | 2009-10-22 | Terratek Inc. | Continuous measurement of heterogeneity of geomaterials |
US20090260415A1 (en) * | 2008-04-16 | 2009-10-22 | Schlumberger Technology Corporation | Apparatus for continuous measurement of heterogeneity of geomaterials |
US20090319243A1 (en) * | 2008-06-18 | 2009-12-24 | Terratek, Inc. | Heterogeneous earth models for a reservoir field |
CN101864939A (zh) * | 2009-04-15 | 2010-10-20 | 中国石油天然气股份有限公司 | 一种裂缝向下延伸酸压方法 |
US20120165215A1 (en) * | 2009-06-26 | 2012-06-28 | The Regents Of The University Of California | Methods and systems for phylogenetic analysis |
WO2012125357A2 (en) * | 2011-03-16 | 2012-09-20 | Saudi Arabian Oil Company | Subsurface directional equalization analysis of rock bodies |
US20130025939A1 (en) * | 2011-07-26 | 2013-01-31 | Denis Heliot | System and method for controlling and monitoring a drilling operation using refined solutions from a panistic inversion |
CN103244087A (zh) * | 2013-05-21 | 2013-08-14 | 中国石油大学(北京) | 一种低渗透油藏调剖堵水选井决策方法 |
US20130231350A1 (en) * | 2011-03-08 | 2013-09-05 | Zalicus Pharmaceuticals Ltd. | Solid dispersion formulations and methods of use thereof |
CA2870735A1 (en) * | 2012-04-20 | 2013-10-24 | Chevron U.S.A. Inc. | System and method for calibrating permeability for use in reservoir modeling |
US20150186570A1 (en) * | 2013-12-31 | 2015-07-02 | Smith International, Inc. | Computing systems, tools, and methods for simulating wellbore departure |
US20150182679A1 (en) * | 2012-07-20 | 2015-07-02 | President And Fellows Of Harvard College | Tissue-engineered pumps and valves and uses thereof |
WO2015127349A1 (en) * | 2014-02-24 | 2015-08-27 | Saudi Arabian Oil Company | Systems, methods, and computer medium to produce efficient, consistent, and high-confidence image-based electrofacies analysis in stratigraphic interpretations across multiple wells |
US20160202390A1 (en) * | 2013-10-01 | 2016-07-14 | Landmark Graphics Corporation | In-Situ Wellbore, Core and Cuttings Information System |
CN106062310A (zh) * | 2014-02-28 | 2016-10-26 | 界标制图有限公司 | 使用无监督分类程序的相定义 |
CN107218024A (zh) * | 2017-05-02 | 2017-09-29 | 中国石油天然气股份有限公司 | 多层低渗透砂岩油藏高含水期菱形反九点井网的调整方法 |
CN107219159A (zh) * | 2017-04-19 | 2017-09-29 | 中国石油天然气股份有限公司 | 一种确定储层非均质性的方法和装置 |
CN107656903A (zh) * | 2017-08-23 | 2018-02-02 | 中国石油天然气股份有限公司 | 数据体的平均指标的确定方法和装置 |
-
2018
- 2018-08-09 CN CN201810900739.6A patent/CN109113732B/zh active Active
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090260883A1 (en) * | 2008-04-16 | 2009-10-22 | Terratek Inc. | Continuous measurement of heterogeneity of geomaterials |
US20090260415A1 (en) * | 2008-04-16 | 2009-10-22 | Schlumberger Technology Corporation | Apparatus for continuous measurement of heterogeneity of geomaterials |
US20090319243A1 (en) * | 2008-06-18 | 2009-12-24 | Terratek, Inc. | Heterogeneous earth models for a reservoir field |
CN101864939A (zh) * | 2009-04-15 | 2010-10-20 | 中国石油天然气股份有限公司 | 一种裂缝向下延伸酸压方法 |
US20120165215A1 (en) * | 2009-06-26 | 2012-06-28 | The Regents Of The University Of California | Methods and systems for phylogenetic analysis |
US20130231350A1 (en) * | 2011-03-08 | 2013-09-05 | Zalicus Pharmaceuticals Ltd. | Solid dispersion formulations and methods of use thereof |
WO2012125357A2 (en) * | 2011-03-16 | 2012-09-20 | Saudi Arabian Oil Company | Subsurface directional equalization analysis of rock bodies |
US20130025939A1 (en) * | 2011-07-26 | 2013-01-31 | Denis Heliot | System and method for controlling and monitoring a drilling operation using refined solutions from a panistic inversion |
CA2870735A1 (en) * | 2012-04-20 | 2013-10-24 | Chevron U.S.A. Inc. | System and method for calibrating permeability for use in reservoir modeling |
US20150182679A1 (en) * | 2012-07-20 | 2015-07-02 | President And Fellows Of Harvard College | Tissue-engineered pumps and valves and uses thereof |
CN103244087A (zh) * | 2013-05-21 | 2013-08-14 | 中国石油大学(北京) | 一种低渗透油藏调剖堵水选井决策方法 |
US20160202390A1 (en) * | 2013-10-01 | 2016-07-14 | Landmark Graphics Corporation | In-Situ Wellbore, Core and Cuttings Information System |
US20150186570A1 (en) * | 2013-12-31 | 2015-07-02 | Smith International, Inc. | Computing systems, tools, and methods for simulating wellbore departure |
WO2015127349A1 (en) * | 2014-02-24 | 2015-08-27 | Saudi Arabian Oil Company | Systems, methods, and computer medium to produce efficient, consistent, and high-confidence image-based electrofacies analysis in stratigraphic interpretations across multiple wells |
CN106062310A (zh) * | 2014-02-28 | 2016-10-26 | 界标制图有限公司 | 使用无监督分类程序的相定义 |
CN107219159A (zh) * | 2017-04-19 | 2017-09-29 | 中国石油天然气股份有限公司 | 一种确定储层非均质性的方法和装置 |
CN107218024A (zh) * | 2017-05-02 | 2017-09-29 | 中国石油天然气股份有限公司 | 多层低渗透砂岩油藏高含水期菱形反九点井网的调整方法 |
CN107656903A (zh) * | 2017-08-23 | 2018-02-02 | 中国石油天然气股份有限公司 | 数据体的平均指标的确定方法和装置 |
Non-Patent Citations (3)
Title |
---|
HAIJIAN YANG: "Nonlinearly preconditioned semismooth Newton methods for variational inequality solution of two-phase flow in porous media", 《JOURNAL OF COMPUTATIONAL PHYSICS》 * |
王子胜: "用动态数据确定缝洞型油藏的有效厚度和渗透率", 《新疆石油地质》 * |
陶诏灵: "概率论与数理统计", 《西北工业大学出版社》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110500069A (zh) * | 2019-08-02 | 2019-11-26 | 中国矿业大学(北京) | 各向异性非均质储层煤层气开采完井技术选择方法 |
CN110500069B (zh) * | 2019-08-02 | 2022-05-17 | 中国矿业大学(北京) | 各向异性非均质储层煤层气开采完井技术选择方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109113732B (zh) | 2022-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zeeb et al. | Evaluation of sampling methods for fracture network characterization using outcrops | |
Wilson et al. | From outcrop to flow simulation: Constructing discrete fracture models from a LIDAR survey | |
Bisdom et al. | Calibrating discrete fracture-network models with a carbonate three-dimensional outcrop fracture network: Implications for naturally fractured reservoir modeling | |
CN110847901B (zh) | 一种变矿化度地层水下的致密砂岩储层流体识别方法 | |
Lochbühler et al. | Structure-coupled joint inversion of geophysical and hydrological data | |
CN104155701B (zh) | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 | |
CN104573333B (zh) | 利用聚类分析进行模型优选的方法 | |
Qu et al. | Incorporation of deformation band fault damage zones in reservoir models | |
CN103064124B (zh) | 一种校正电磁勘探中地形影响的比值方法 | |
CN105593703A (zh) | 根据基于核心的岩石定型构建并通过岩石置换建模增强的模拟到地震工作流程 | |
CN104632202B (zh) | 确定干粘土三孔隙度测井参数值的方法及装置 | |
CN105425292A (zh) | 一种油气预测方法及装置 | |
CN107966732A (zh) | 基于空间结构导向的地震属性变化率求取方法 | |
CN107807396B (zh) | 地层基质电阻率的确定方法和装置 | |
CN105388525A (zh) | 储层预测方法及装置 | |
Chesnaux et al. | Regional fracture network permeability using outcrop scale measurements | |
CN108663714A (zh) | 一种沉积微相刻画方法 | |
Jonah et al. | Validity of an empirical rule for delineating aquifer prospects at the Gidan Kwano Campus Development Phase II, Federal University of Technology, Minna, Northcentral Nigeria | |
Xu et al. | Petrophysical rock classification in the Cotton Valley tight-gas sandstone reservoir with a clustering pore-system orthogonality matrix | |
Colombera et al. | Models for guiding and ranking well-to-well correlations of channel bodies in fluvial reservoirs | |
CN102928871B (zh) | 一种基于旋转菱形的属性提取与断层描述方法 | |
Jamshidi Gohari et al. | Improving permeability estimation of carbonate rocks using extracted pore network parameters: a gas field case study | |
CN109113732A (zh) | 储层非均质性的确定方法及装置 | |
Phelps et al. | Exploring viable geologic interpretations of gravity models using distance-based global sensitivity analysis and kernel methods | |
CN107991705B (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 |