CN111898091B - 油水界面的确定方法及装置 - Google Patents
油水界面的确定方法及装置 Download PDFInfo
- Publication number
- CN111898091B CN111898091B CN202010574623.5A CN202010574623A CN111898091B CN 111898091 B CN111898091 B CN 111898091B CN 202010574623 A CN202010574623 A CN 202010574623A CN 111898091 B CN111898091 B CN 111898091B
- Authority
- CN
- China
- Prior art keywords
- oil
- depth value
- layer
- water interface
- water
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Bioinformatics & Computational Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种油水界面的确定方法及装置,其中,该方法包括:根据饱和度平面分布图中油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值;以所述估读油水界面深度值作为初始值,循环确定最优油水界面深度值,每个循环周期均执行以下操作:确定每一小层中深度值大于当前周期油水界面深度值的油层数据点数目与油层数据点总数目的第一比值,及深度值大于当前周期油水界面深度值的油水同层数据点数目与油水同层数据点总数目的第二比值;将第一比值、第二比值的平均值反向代入各个小层数据中,得到各个小层的下一周期的油水界面深度值,直到得到满足预设条件的最优油水界面深度值。上述技术方案提高了确定油水界面的精度。
Description
技术领域
本发明涉及油气田开发技术领域,特别涉及一种油水界面的确定方法及装置。
背景技术
在油藏中,由于流体的分异调整作用,石油占据油藏的高部位,水体则位于油藏的底部或边部。石油与水体之间的接触面,即称为油水界面。确定油水界面对油水系统的判别、储量计算、井位部署和射孔方案设计等至关重要。在构造落实和储层分布基本清楚的条件下,依据对油藏类型的正确认识,按油水界面水平外推并在相应的油层顶面构造图上圈定油藏边界,能够确保含油面积的准确性。因此,落实油水界面至关重要。
目前,国内外确定油水界面的主要方法有压力测试法、测井解释法、试油法以及其他间接测试方法。然而,现有方法对于资料精度要求较高,需要在油田投产初期获取系统规范的测试结果,但此种情况面对较为复杂油水界面或多层多油水系统的状况时往往难以奏效。而且,在开发中后期,钻井数量逐渐增加,且井位可更好控制整个含油面积,甚至部分拓边井逐步进入油水界面之外的区域,这给综合判定油水界面提供了资料便利,但同时,油藏生产对现有测试方法用到的压力、原始含油饱和度数据影响较大,再加之饱和度解释本身的不确定性,导致后续RFT线性回归和饱和度判定油水界面的结果不能准确代表地层真实状况。
针对上述问题,目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种油水界面的确定方法,用以提高确定油水界面的精度,该方法包括:
根据饱和度平面分布图中油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值;
以所述估读油水界面深度值作为初始值,循环确定最优油水界面深度值,每个循环周期均执行以下操作:
根据各小层的当前周期油水界面深度值,以及每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值,确定每一小层中深度值大于对应小层的当前周期油水界面深度值的油层数据点数目与对应小层的油层数据点总数目的第一比值,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目与对应小层的油水同层数据点总数目的第二比值;
根据每一小层对应的第一比值,得到所有小层关于油层数据的第一比值的平均值;根据每一小层对应的第二比值,得到所有小层关于油水同层数据的第二比值的平均值;
将第一比值的平均值和第二比值的平均值反向代入各个小层数据中,得到各个小层的下一周期的油水界面深度值,直到得到满足预设条件的油水界面深度值作为最优油水界面深度值。
本发明实施例还提供了一种油水界面的确定装置,用以提高确定油水界面的精度,该装置包括:
获取单元,用于根据饱和度平面分布图中油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值;
确定单元,用于以所述估读油水界面深度值作为初始值,循环确定最优油水界面深度值,每个循环周期均执行以下操作:
根据各小层的当前周期油水界面深度值,以及每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值,确定每一小层中深度值大于对应小层的当前周期油水界面深度值的油层数据点数目与对应小层的油层数据点总数目的第一比值,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目与对应小层的油水同层数据点总数目的第二比值;
根据每一小层对应的第一比值,得到所有小层关于油层数据的第一比值的平均值;根据每一小层对应的第二比值,得到所有小层关于油水同层数据的第二比值的平均值;
将第一比值的平均值和第二比值的平均值反向代入各个小层数据中,得到各个小层的下一周期的油水界面深度值,直到得到满足预设条件的油水界面深度值作为最优油水界面深度值。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述油水界面的确定方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述油水界面的确定方法的计算机程序。
本发明实施例提供的技术方案通过:根据饱和度平面分布图中油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值;确定单元,用于以所述估读油水界面深度值作为初始值,循环确定最优油水界面深度值,每个循环周期均执行以下操作:根据各小层的当前周期油水界面深度值,以及每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值,确定每一小层中深度值大于对应小层的当前周期油水界面深度值的油层数据点数目与对应小层的油层数据点总数目的第一比值,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目与对应小层的油水同层数据点总数目的第二比值;根据每一小层对应的第一比值,得到所有小层关于油层数据的第一比值的平均值;根据每一小层对应的第二比值,得到所有小层关于油水同层数据的第二比值的平均值;将第一比值的平均值和第二比值的平均值反向代入各个小层数据中,得到各个小层的下一周期的油水界面深度值,直到得到满足预设条件的油水界面深度值作为最优油水界面深度值,实现了定量和定性确定油水界面,结合统计数据结果对估读油水界面深度值进行修正反演得到最优油水界面深度值,提高了确定油水界面的精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例中油水界面的确定方法的流程示意图;
图2是本发明实施例中小层砂体饱和度解释点与估读油水界面相对位置关系示意图;
图3是本发明实施例中小层内估读油水界面海拔与高于界面深度值数据百分比对比示意图;
图4是本发明实施例中各小层内百分比值与指导平均值对比示意图;
图5是本发明实施例中各层估读油水界面海拔高度示意图;
图6是本发明实施例中模型中各层油水界面示意图;
图7是本发明实施例中油水界面的确定装置的结构示意图;
图8是本发明实施例中计算机设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
由于发明人发现了背景技术提到的技术问题,综合利用油田开发中后期丰富的测井、测压和生产数据,提出了一种基于数理统计综合判断油藏油水界面的新方案,该方案依靠经验法在各小层平面构造图上依据层内饱和度解释结果(一般小层含水饱和度解释结果分为油层、水层和处于两者之间的油水同层)读取油水界面,结合统计数据结果,对估读结果进行修正反演,提高油水界面判定精度。具体地,现有方法中在后期压力梯度测试结果受干扰的情况下,仅凭经验估读油水界面误差较大,并且需要估读者具备丰富的油水界面判断经验和扎实的测井解释功底,无法系统推广。而本发明实施例提供的基于数据统计结果的研究方法更加具有系统性,它能够大幅度减小个别坏点数据对于整体趋势的影响,把定量和定性方法结合起来,利用数据点丰富的优势,应用数据统计的方法,修正研究者个人经验不足导致的误差,使得依据经验读图的结果获得更加合理的应用,更符合客观实际,更易于定量表示,同时可以把饱和度解释、测井曲线、压力数据和生产动态数据充分利用起来,得出更可靠的油水界面判别结果。
因此,综上所述,本发明实施例提供的油水界面确定方案为一种基于数据统计结果结合饱和度平面分布图估算油水界面的方案,即基于饱和度解释统计的油水界面确定方案。下面对该油水界面的确定方案进行详细介绍如下。
图1是本发明实施例中油水界面的确定方法的流程示意图,如图1所示,该方法包括如下步骤:
步骤101:根据饱和度平面分布图中油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值;
步骤102:以所述估读油水界面深度值作为初始值,循环确定最优油水界面深度值,每个循环周期均执行以下操作:
步骤1021:根据各小层的当前周期油水界面深度值,以及每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值,确定每一小层中深度值大于对应小层的当前周期油水界面深度值的油层数据点数目与对应小层的油层数据点总数目的第一比值,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目与对应小层的油水同层数据点总数目的第二比值;
步骤1022:根据每一小层对应的第一比值,得到所有小层关于油层数据的第一比值的平均值;根据每一小层对应的第二比值,得到所有小层关于油水同层数据的第二比值的平均值;
步骤1023:将第一比值的平均值和第二比值的平均值反向代入各个小层数据中,得到各个小层的下一周期的油水界面深度值,直到得到满足预设条件的油水界面深度值作为最优油水界面深度值。
本发明实施例提供的油水界面的确定方法实现了定量和定性确定油水界面,结合统计数据结果对估读油水界面深度值进行修正反演得到最优油水界面深度值,提高了确定油水界面的精度。
本发明实施例提供的油水界面的确定方法进一步可以包括:应用所述最优油水界面深度值进行石油天然气钻井,例如进行油水系统的判别、储量计算、井位部署和射孔方案设计等。
下面结合附图2至图4,对本发明实施涉及的各个步骤进行详细介绍如下。
一、首先介绍上述步骤101。
在一个实施例中,根据饱和度平面分布图中油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值,可以包括:基于饱和度平面分布图,去除干扰数据,根据去除干扰数据后的油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值。
具体实施时,去除干扰数据的原理在于每个小层依据一般的单一油藏油水分布理论,水样点(水层数据点)不会存在于油水界面以上远离油水界面的地方。去除干扰数据后,再获取油藏各小层的估读油水界面深度值,进一步提高了确定油水界面的精度。
在一个实施例中,基于饱和度平面分布图,去除干扰数据,根据去除干扰数据后的油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值,可以包括:
基于饱和度平面分布图,去除油藏内部远离油水界面的水层数据点,根据去除油藏内部远离油水界面的水层数据点后的油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值。
具体实施时,去除油藏内部远离油水界面的水层数据点后,再获取油藏各小层的估读油水界面深度值,进一步提高了确定油水界面的精度。
具体实施时,如图2所示,基于饱和度平面分布图(各小层数据点在等高线图上的分布),去除干扰数据,初步估读油藏各小层沿构造线的油水界面高度(深度)。干扰点去除时,考虑到后期注水开发过程中注入水会导致局部水淹,因而去除油藏内部远离油水界面的水层,考虑到油田开发过程中对压力的干扰可能导致边底水突入原始油水界面内部,因而估读油水界面时,尽可能以油层和油水同层解释点(数据点)分布为判断依据,在选择数据点时,选择油层和油水同层的数据点数目远大于选择水层数据点的数目,具体地,两者之差大于总数据点数目的百分之九十,其理由在于:油相解释结果总会出现在含油区域之内,这样就大大扩展了可选数据范围,打破了此前确定原始油水界面只能使用油田开发初期井点数据的约束。
具体实施时,估读油水界面人为判断一下,根据等高线图上的数据分布,读出来一条等高线使其高度以上范围能尽可能包含全部油点和油水同层数据点(因为解释,分层,地层物性等等原因,会有例外情况),并且让水层数据尽可能在这条线以下(水层点不一定,因为注水会导致边水侵入,因此发明人考虑更多以油和油水同层数据点为主要判断依据)各小层的平均比例值的做法(详见下面步骤102的介绍),则尽可能从统计角度尽可能减小不同小层油水界面高度读值的不确定性。
二、其次,介绍上述步骤102。
1、首先介绍步骤1021。
在一个实施例中,根据各小层的当前周期油水界面深度值,以及每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值,确定每一小层中深度值大于对应小层的当前周期油水界面深度值的油层数据点数目与对应小层的油层数据点总数目的第一比值,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目与对应小层的油水同层数据点总数目的第二比值,可以包括:
获取每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值;
将每一油层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目;将每一油水同层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目;
根据深度值大于对应小层的当前周期油水界面深度值的油层数据点数目,对应小层的油层数据点总数目,得到第一比值;根据深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目,对应小层的油水同层数据点总数目,得到第二比值。
具体实施时,上述确定关于油层数据的第一比值、关于油水同层数据的第二比值的实施方案,进一步提高了确定油水界面的精度。
在一个实施例中,如图3所示,将每一油层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目;将每一油水同层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目,可以包括:
绘制多个油层数据点的深度值和多个油水同层数据点的深度值,每一深度值出现的数目,以及当前周期油水界面深度值的关系图;
根据所述关系图,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目。
具体实施时,上述通过绘制如图3所示的关系图,来确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目、深度值大于对应小层的当前周期油水界面深度值的水层数据点数目,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目的实施方案,直观、清晰,提高了确定油水界面的效率,同时也进一步提高了确定油水界面的精度。
在一个实施例中,所述第一比值可以为第一百分比值,所述第二比值可以为第二百分比值,所述第三比值可以为第三百分比值。
具体实施时,各个比值为百分比值,直观、清晰,方便得到精确的油水界面。
具体实施时,在上述步骤1021中,对每个小层的油层、水层和油水同层数据点开展测点深度统计分析,基于统计分析的频率分布和估读油水界面深度相对关系(该关系可以是图3),得出各小层内各类饱和度解释数据(解释点、数据点)中高于估测界面深度值(估读油水界面深度值)的数据百分比,汇总各小层三个百分比数,作为后续综合判断的依据。
2、接着为了方便介绍,一同介绍上述步骤1022和步骤1023。
在一个实施例中,所述第一比值的平均值可以为第一百分比值的平均值,所述第二比值的平均值可以为第二百分比值的平均值,所述第三比值的平均值可以为第三百分比值的平均值。
具体实施时,各个平均值为百分比值的平均值,直观、清晰,方便得到精确的油水界面。
具体实施时,根据统计所得油层、水层和油水同层的百分比值,计算层系(1个层系包含了多个小层)内各油层整体各相高于油水界面海拔高度(深度)百分比平均值,得出如图4中三条平均值线,观察可知油层和油水同层高于油水界面的分布比例相对稳定,而水层由于后期开发可能导致边底水入侵,因而其结果相对比较分散,也可以理解为,使用饱和度解释中油层和油水同层的解释结果来约束油水界面更为合理。将计算出的各小层平均百分比值反向代入各个小层数据中,可算得基于平均百分比指导值的油水分界线高度(上述提到的下一周期的油水界面深度值,即更新的油水界面深度值),具体地,反向代入指此前根据估读油水界面深度值确定了百分比值,现在根据百分比值,按照数据在油藏中的位置从上到下数,数到达到该百分比值时对应的深度值即得所要求的高度(下一周期的油水界面深度值),此时将这个百分比值认为是百分位数。由于数据分布方式(每类数据点的深度频率分布图)接近正态分布,其平均值约等于众数和中位数。在标准差相对较小的情况下,可以认为平均值结果在置信区间以内。对比估读分界线深度(估读油水界面深度值),基于差值再次对油水界线(油水界面)进行调整。多次重复以上平均值计算和反向代入的步骤,直至油水界面达到合理值,即得到最优最优油水界面深度值。
具体实施时,上述步骤1023中的预设条件可以根据大量油田数据实验得到,例如两个周期之间数据变化值小于一阈值(例如0.5%)等。
四、接着,介绍后续进一步提高确定油水界面的精度的方案。
在一个实施例中,上述油水界面的确定方法还可以包括:根据可获得的RFT压力数据或测井数据,进行交叉对比处理,验证所述最优油水界面深度值的可信度。
具体实施时,根据可获得的RFT压力数据,或其他测井数据,交叉对比,进而验证所述最优油水界面深度值的可信度,进一步提高了确定油水界面的精度。
具体实施时,RFT是钻井时测得的地层原始压力数据,压力拐点就是油水界面,测井曲线也可以认识出油水界面。不同方法的油水界面深度往往略有差异,交叉验证的过程是看哪个方法更可靠。
下面结合附图5至图6,再举一例子便于理解本发明如何实施。
在实际油田应用中,应当根据实际生产资料,调整选取数据以及油水界面测定方式。采用本发明实施例提供的油水界面确定方法,对青海柴达木盆地某油田进行了油水界面测定。该区块开发时间较早,油田资料不全且地下油水状况分布复杂,不具备统一油水界面;开发进行到中后期导致对油藏原始含油饱和度以及压力等测试结果产生了较大影响。油藏区块内小层较多,并且由于地质情况复杂形成了各层不一致的油水界面,对于直接使用现有测井法或压力测试法测定油水界面有一定困难。本实施例中最终复算结果与估测油水界面深度如下表1所示。
小层号 | 复算水相 | 复算过度相 | 复算油相 | 油相过度区平均 | 估测值 |
1 | 1168 | 1178 | 1173 | 1175.5 | 1175 |
2 | 1169 | 1190 | 1170 | 1180 | 1175 |
3 | 1168 | 1175 | 1166 | 1170.5 | 1170 |
4 | 1153 | 1173 | 1159 | 1166 | 1160 |
5 | 1171 | 1168 | 1154 | 1161 | 1155 |
6 | 1157 | 1162 | 1151 | 1156.5 | 1160 |
7 | 1145 | 1175 | 1162 | 1168.5 | 1165 |
8 | 1133 | 1153 | 1162 | 1157.5 | 1165 |
9 | 1148 | 1146 | 1153 | 1149.5 | 1160 |
10 | 1133 | 1144 | 1143 | 1143.5 | 1150 |
表1
由于该地区存在两处不规则水体和后期注水开发的影响,导致水相油水界面误差较大,参考价值较小。统计结果表明,按照本发明实施例提供的油水界面确定方法:基于油相与油水同层数据点确定的油水界面较为一致且精度较高,可以作为该油藏油水界面认识的最终结果,该区域由于剧烈的构造运动和局部水体的影响,油水分布复杂,整体呈现随深度逐渐降低的规则变化规律。各层估读油水界面海拔高度图如图5所示,模型中各层油水界面图如图6所示。
综上,根据实际应用证明,利用本发明实施提供的油水界面的确定方法得到的油水界面精度高。
基于同一发明构思,本发明实施例中还提供了一种油水界面的确定装置,如下面的实施例所述。由于油水界面的确定装置解决问题的原理与油水界面的确定方法相似,因此油水界面的确定装置的实施可以参见油水界面的确定方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图7是本发明实施例中油水界面的确定装置的结构示意图,如图7所示,该装置包括:
获取单元01,用于根据饱和度平面分布图中油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值;
确定单元02,用于以所述估读油水界面深度值作为初始值,循环确定最优油水界面深度值,每个循环周期均执行以下操作:
根据各小层的当前周期油水界面深度值,以及每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值,确定每一小层中深度值大于对应小层的当前周期油水界面深度值的油层数据点数目与对应小层的油层数据点总数目的第一比值,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目与对应小层的油水同层数据点总数目的第二比值;
根据每一小层对应的第一比值,得到所有小层关于油层数据的第一比值的平均值;根据每一小层对应的第二比值,得到所有小层关于油水同层数据的第二比值的平均值;
将第一比值的平均值和第二比值的平均值反向代入各个小层数据中,得到各个小层的下一周期的油水界面深度值,直到得到满足预设条件的油水界面深度值作为最优油水界面深度值。
在一个实施例中,所述获取单元具体可以用于:基于饱和度平面分布图,去除干扰数据,根据去除干扰数据后的油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值。
在一个实施例中,所述获取单元具体可以用于:基于饱和度平面分布图,去除油藏内部远离油水界面的水层数据点,根据去除油藏内部远离油水界面的水层数据点后的油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值。
在一个实施例中,所述第一比值为第一百分比值,所述第二比值为第二百分比值;所述第一比值的平均值为第一百分比值的平均值,所述第二比值的平均值为第二百分比值的平均值。
在一个实施例中,上述油水界面的确定装置还可以包括:验证单元,用于根据可获得的RFT压力数据或测井数据,进行交叉对比处理,验证所述最优油水界面深度值的可信度。
在一个实施例中,根据各小层的当前周期油水界面深度值,以及每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值,确定每一小层中深度值大于对应小层的当前周期油水界面深度值的油层数据点数目与对应小层的油层数据点总数目的第一比值,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目与对应小层的油水同层数据点总数目的第二比值,包括:
获取每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值;
将每一油层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目;将每一油水同层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目;
根据深度值大于对应小层的当前周期油水界面深度值的油层数据点数目,对应小层的油层数据点总数目,得到第一比值;根据深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目,对应小层的油水同层数据点总数目,得到第二比值。
在一个实施例中,将每一油层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目;将每一油水同层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目,包括:
绘制多个油层数据点的深度值和多个油水同层数据点的深度值,每一深度值出现的数目,以及当前周期油水界面深度值的关系图;
根据所述关系图,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目。
本发明实施例还提供了一种计算机设备,如图8所示,包括存储器302、处理器304及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述油水界面的确定方法。
具体的,该计算机设备可以是计算机终端、服务器或者类似的运算装置。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述任意的油水界面的确定方法的计算机程序。
具体的,计算机可读存储介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机可读存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读存储介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
本发明实施例提供的油水界面的确定方案达到的有益技术效果是:实现了定量和定性确定油水界面,结合统计数据结果对估读油水界面深度值进行修正反演得到最优油水界面深度值,提高了确定油水界面的精度。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (14)
1.一种油水界面的确定方法,其特征在于,包括:
根据饱和度平面分布图中油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值;
以所述估读油水界面深度值作为初始值,循环确定最优油水界面深度值,每个循环周期均执行以下操作:
根据各小层的当前周期油水界面深度值,以及每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值,确定每一小层中深度值大于对应小层的当前周期油水界面深度值的油层数据点数目与对应小层的油层数据点总数目的第一比值,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目与对应小层的油水同层数据点总数目的第二比值;
根据每一小层对应的第一比值,得到所有小层关于油层数据的第一比值的平均值;根据每一小层对应的第二比值,得到所有小层关于油水同层数据的第二比值的平均值;
将第一比值的平均值和第二比值的平均值反向代入各个小层数据中,得到各个小层的下一周期的油水界面深度值,直到得到满足预设条件的油水界面深度值作为最优油水界面深度值。
2.如权利要求1所述的油水界面的确定方法,其特征在于,根据饱和度平面分布图中油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值,包括:基于饱和度平面分布图,去除干扰数据,根据去除干扰数据后的油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值。
3.如权利要求2所述的油水界面的确定方法,其特征在于,基于饱和度平面分布图,去除干扰数据,根据去除干扰数据后的油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值,包括:
基于饱和度平面分布图,去除油藏内部远离油水界面的水层数据点,根据去除油藏内部远离油水界面的水层数据点后的油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值。
4.如权利要求1所述的油水界面的确定方法,其特征在于,所述第一比值为第一百分比值,所述第二比值为第二百分比值;所述第一比值的平均值为第一百分比值的平均值,所述第二比值的平均值为第二百分比值的平均值。
5.如权利要求1所述的油水界面的确定方法,其特征在于,根据各小层的当前周期油水界面深度值,以及每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值,确定每一小层中深度值大于对应小层的当前周期油水界面深度值的油层数据点数目与对应小层的油层数据点总数目的第一比值,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目与对应小层的油水同层数据点总数目的第二比值,包括:
获取每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值;
将每一油层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目;将每一油水同层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目;
根据深度值大于对应小层的当前周期油水界面深度值的油层数据点数目,对应小层的油层数据点总数目,得到第一比值;根据深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目,对应小层的油水同层数据点总数目,得到第二比值。
6.如权利要求5所述的油水界面的确定方法,其特征在于,将每一油层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目;将每一油水同层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目,包括:
绘制多个油层数据点的深度值和多个油水同层数据点的深度值,每一深度值出现的数目,以及当前周期油水界面深度值的关系图;
根据所述关系图,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目。
7.一种油水界面的确定装置,其特征在于,包括:
获取单元,用于根据饱和度平面分布图中油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值;
确定单元,用于以所述估读油水界面深度值作为初始值,循环确定最优油水界面深度值,每个循环周期均执行以下操作:
根据各小层的当前周期油水界面深度值,以及每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值,确定每一小层中深度值大于对应小层的当前周期油水界面深度值的油层数据点数目与对应小层的油层数据点总数目的第一比值,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目与对应小层的油水同层数据点总数目的第二比值;
根据每一小层对应的第一比值,得到所有小层关于油层数据的第一比值的平均值;根据每一小层对应的第二比值,得到所有小层关于油水同层数据的第二比值的平均值;
将第一比值的平均值和第二比值的平均值反向代入各个小层数据中,得到各个小层的下一周期的油水界面深度值,直到得到满足预设条件的油水界面深度值作为最优油水界面深度值。
8.如权利要求7所述的油水界面的确定装置,其特征在于,所述获取单元具体用于:基于饱和度平面分布图,去除干扰数据,根据去除干扰数据后的油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值。
9.如权利要求8所述的油水界面的确定装置,其特征在于,所述获取单元具体用于:基于饱和度平面分布图,去除油藏内部远离油水界面的水层数据点,根据去除油藏内部远离油水界面的水层数据点后的油藏各小层的饱和度解释结果,获取油藏各小层的估读油水界面深度值。
10.如权利要求7所述的油水界面的确定装置,其特征在于,所述第一比值为第一百分比值,所述第二比值为第二百分比值;所述第一比值的平均值为第一百分比值的平均值,所述第二比值的平均值为第二百分比值的平均值。
11.如权利要求7所述的油水界面的确定装置,其特征在于,根据各小层的当前周期油水界面深度值,以及每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值,确定每一小层中深度值大于对应小层的当前周期油水界面深度值的油层数据点数目与对应小层的油层数据点总数目的第一比值,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目与对应小层的油水同层数据点总数目的第二比值,包括:
获取每一小层中的多个油层数据点的深度值和多个油水同层数据点的深度值;
将每一油层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目;将每一油水同层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目;
根据深度值大于对应小层的当前周期油水界面深度值的油层数据点数目,对应小层的油层数据点总数目,得到第一比值;根据深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目,对应小层的油水同层数据点总数目,得到第二比值。
12.如权利要求11所述的油水界面的确定装置,其特征在于,将每一油层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目;将每一油水同层数据点的深度值与对应小层的当前周期油水界面深度值进行比较,确定深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目,包括:
绘制多个油层数据点的深度值和多个油水同层数据点的深度值,每一深度值出现的数目,以及当前周期油水界面深度值的关系图;
根据所述关系图,确定深度值大于对应小层的当前周期油水界面深度值的油层数据点数目,以及深度值大于对应小层的当前周期油水界面深度值的油水同层数据点数目。
13.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至6任一所述方法。
14.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至6任一所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010574623.5A CN111898091B (zh) | 2020-06-22 | 2020-06-22 | 油水界面的确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010574623.5A CN111898091B (zh) | 2020-06-22 | 2020-06-22 | 油水界面的确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111898091A CN111898091A (zh) | 2020-11-06 |
CN111898091B true CN111898091B (zh) | 2023-09-26 |
Family
ID=73207792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010574623.5A Active CN111898091B (zh) | 2020-06-22 | 2020-06-22 | 油水界面的确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111898091B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104252551A (zh) * | 2013-06-28 | 2014-12-31 | 中国石油化工股份有限公司 | 不规则油水界面建模实现方法 |
CN106777524A (zh) * | 2016-11-24 | 2017-05-31 | 成都北方石油勘探开发技术有限公司 | 一种非统一油水界面油藏储量快速计算方法 |
CN107589469A (zh) * | 2017-08-04 | 2018-01-16 | 中国石油天然气股份有限公司 | 油水界面的确定方法和装置 |
CN110056346A (zh) * | 2019-04-17 | 2019-07-26 | 中国石油大港油田勘探开发研究院 | 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1729151B1 (en) * | 2005-05-27 | 2011-06-29 | Services Petroliers Schlumberger | Method and apparatus for measuring the wettability of geological formations |
US10557783B2 (en) * | 2015-03-27 | 2020-02-11 | Schlumberger Technology Corporation | Method for determining equilibrium wettability of an interface between a void space and a solid phase of a rock sample |
-
2020
- 2020-06-22 CN CN202010574623.5A patent/CN111898091B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104252551A (zh) * | 2013-06-28 | 2014-12-31 | 中国石油化工股份有限公司 | 不规则油水界面建模实现方法 |
CN106777524A (zh) * | 2016-11-24 | 2017-05-31 | 成都北方石油勘探开发技术有限公司 | 一种非统一油水界面油藏储量快速计算方法 |
CN107589469A (zh) * | 2017-08-04 | 2018-01-16 | 中国石油天然气股份有限公司 | 油水界面的确定方法和装置 |
CN110056346A (zh) * | 2019-04-17 | 2019-07-26 | 中国石油大港油田勘探开发研究院 | 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111898091A (zh) | 2020-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109034447B (zh) | 一种水驱油藏含水上升率变化的预测方法及装置 | |
CN107292074B (zh) | 一种井间连通性判断方法 | |
CN111665560B (zh) | 油气层识别方法、装置、计算机设备及可读存储介质 | |
CN107832574B (zh) | 一种基于随钻测井的水平井水淹层解释方法 | |
CN107956465B (zh) | 基于关联井的全区多井测井曲线标准化方法及装置 | |
CN111898091B (zh) | 油水界面的确定方法及装置 | |
CN110782187A (zh) | 基于含水率的枯竭油气藏型储气库动态评价方法及系统 | |
CN110632654B (zh) | 断块圈闭含油边界确定方法及装置 | |
CN113775326A (zh) | 可动水饱和度的评价方法、装置、电子设备及介质 | |
US20220342108A1 (en) | Method, device and medium for acquiring logging parameters | |
CN111241652A (zh) | 一种确定地层原油粘度的方法及装置 | |
CN114060015B (zh) | 一种致密砂岩含气性的评价方法及装置 | |
CN112711076B (zh) | 提取石油钻井中泥浆侵入地层深度的方法和装置 | |
CN111814015A (zh) | 煤层注水效果的检测方法、装置、电子设备及存储介质 | |
CN110984975A (zh) | 一种判断储层流体连通性的方法、设备及可读存储介质 | |
Tømmerås et al. | Prewell and postwell predictions of oil and gas columns using an iterative Monte Carlo technique with three-dimensional petroleum systems modeling | |
CN109033508A (zh) | 一种油藏含水率及含水上升率的确定方法和装置 | |
CN108920781A (zh) | 一种油藏含水率及含水上升率的确定方法和装置 | |
CN114139386B (zh) | 一种基于多层合采的连通性评估方法及装置 | |
CN109063228A (zh) | 一种水驱油藏含水上升率变化的确定方法及装置 | |
CN117991361A (zh) | 流体变化指示因子提取方法、装置、电子设备及存储介质 | |
CN112883530B (zh) | 理论水驱曲线确定方法和装置 | |
CN113740928B (zh) | 一种致密砂岩储层的流体识别方法 | |
CN114183130B (zh) | 基于电测井的碳同位素预测方法、装置、电子设备及介质 | |
CN113738319B (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 |