CN103902751B - 油藏数值模拟的边界处理方法及装置 - Google Patents
油藏数值模拟的边界处理方法及装置 Download PDFInfo
- Publication number
- CN103902751B CN103902751B CN201210575346.5A CN201210575346A CN103902751B CN 103902751 B CN103902751 B CN 103902751B CN 201210575346 A CN201210575346 A CN 201210575346A CN 103902751 B CN103902751 B CN 103902751B
- Authority
- CN
- China
- Prior art keywords
- region
- reservoir
- interface
- oil
- dirty oil
- 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.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种油藏模拟方法,所述方法包括:区域划分装置,用于将油藏划分为多个区域,所述多个区域至少包括组分区域和黑油区域;区域划分装置,用于对每一个区域进行网格划分,并在区域之间形成界面;映射装置,用于将所述组分区域以及所述黑油区域分别映射至中间媒介物质以建立关联;界面处理装置,用于根据物质守恒定律,对所述界面进行界面处理;矩阵构建装置,根据所述多个区域的映射结果以及所述界面处理结果构建统一矩阵;矩阵求解装置,用于对所述统一矩阵进行求解。采用本申请的方法,可以在模拟时提供不同区块之间数据的实时交换,杜绝了区域之间流体的不连续性,还可提高了模拟的效率,增加了以后模拟器的扩展性。
Description
技术领域
本发明涉及一种黑油、组分一体化的油藏数值模拟技术,更具体地,本发明涉及一体化的油藏模拟中边界处理的方法及装置。
背景技术
油藏数值模拟技术在石油工业的应用起源于20世纪30年代,其基础理论是基于达西渗流定律。近年来随着计算机技术的飞速发展,现有的商用油藏数值模拟技术在功能配套和集成方面虽有了长足进步,但对多相流体流动现象的综合描述上,以及对油气从地下到地面的生产过程量化模型上,仍在沿用“黑油”,“组分”和“热采”的简化建模理念与定制方法。后黑油、组分、热采时代的数模研发工作,也即新一代油藏数值模拟器的研发起始于本世纪初,诸多的国外软件服务公司和大石油公司均投入了巨大的人力物力,取得了显著的进展。主要工作集中在以下三个领域:
数值模拟计算技术:计算技术的发展主要体现在基础模型的扩展、各类复杂网格生成、高性能解法、及物理模型集成等方面。新型、高性能解法器可提高非结构网格和节点异性体系的求解稳定性和速度,同时降低对平台和操作系统的要求。
多学科技术集成:多学科集成的发展已从勘探开发一体化平台的研发逐步向多种学科技术相辅相成方向延伸,跨越油藏的一体化耦合与集成模型已形成雏形。渐进、混合式更新流程,实时模拟、整体优化、自修正算法,以及辅助、自动历史拟合算法均获得了不同程度的成功运用。从地质到工程,从数据到决策的工作模式开始取代传统的单向模拟模式。
工作流程和使用方法:基于数据流的功能设计及面向目标(OOP)和面向服务(SOA)的体系结构已成为现代模拟器的主流软件架构方式。一键式的数据加载和智能使用导引降低了用户的使用门槛和难度。与传统模拟器相比,新架构在软件的可维护性、扩展性和再开发性上获得了跨越式的进步。
中国油气藏模拟软件的研发已有20年以上的历史,由于未能有效地建立油藏模拟持续发展的体系、系统地建立模拟算法的核心技术,目前国内三大石油公司主要依赖于引进的油气藏数值模拟软件。国外软件在特定功能、模拟流程、模型数据等方面,常常不符合国内复杂地质条件的特定需求,加之软件的使用难度和价格昂贵等因素,在很大程度上限制了数值模拟技术在国内的广泛应用,中国石油企业提高本国石油可采性和总产量的能力上受到了一定程度的限制。
而国外油藏数值模拟软件开发商均具有30年以上的技术积累,掌握并垄断着诸多的模拟计算核心技术。许多计算处理技巧是在长期的研发、服务、应用中逐步探索得到的,并未公开发表。但商业上流行的通常是对不同类型的油藏开发出相应的模拟器,例如黑油模拟器、组份模拟器、热采模拟器、化学驱模拟器、流线法模拟器等。
举例而言,斯伦贝谢ELIPSE油藏模拟软件中有ECLIPSE100黑油模拟器、ECLIPSE300组份模拟器、FRONTSIM流线法等。其中,E100黑油模拟器,适用于黑油、挥发性油、干气、湿气等油藏。E300组分模拟器,当相变化非常重要时使用组分模拟器适用于凝析气藏、挥发性油藏、注气开发等油藏。E500热采模拟器,适用于包含油、气、水三相的稠油热采。FrontSim流线模拟器,更快速的模拟流程、流线的直观显示、减少建模的不确定性、减少历史拟合时间。
但是无论是斯伦贝谢旗舰产品Eclipse模拟器系列下一代的Intersect,还是兰德马克的新一代油藏数值模拟器Nexus,虽然将常规的三相(油、气、水)流体流动扩展到了四相(固相),均只能进行同种类的油藏模拟,例如黑油、组分等,虽然也可以对油藏进行区域分解,但是只是针对同种类油藏,而且区域的边界基本都是连续的,无需特别处理。
综上所述,现有的油藏模拟技术均是依据单一的模型,而本申请就是基于黑油和组分油藏混合模拟油藏中,特别的是,需要给出不同区域边界处的处理方案,故现有技术仍然存在着改进的空间。
发明内容
考虑到上述存在的问题,本发明的目的之一在于提供在同时考虑黑油模型以及组分模型的油藏模拟技术方案时进行边界处理的技术方案。
根据本发明的一个方面,提供了一种油藏模拟方法,所述方法包括:将油藏划分为多个区域,所述多个区域至少包括组分区域和黑油区域;对每一个区域进行网格划分,并在区域之间形成界面;将所述组分区域以及所述黑油区域分别映射至中间媒介物质以建立关联;根据物质守恒定律,对所述界面进行界面处理;根据所述多个区域的映射结果以及所述界面处理结果构建统一矩阵;对所述统一矩阵进行求解。
根据本发明的另一方面,提供了一种油藏模拟系统,所述系统包括:区域划分装置,用于将油藏划分为多个区域,所述多个区域至少包括组分区域和黑油区域;区域划分装置,用于对每一个区域进行网格划分,并在区域之间形成界面;映射装置,用于将所述组分区域以及所述黑油区域分别映射至中间媒介物质以建立关联;界面处理装置,用于根据物质守恒定律,对所述界面进行界面处理;矩阵构建装置,根据所述多个区域的映射结果以及所述界面处理结果构建统一矩阵;矩阵求解装置,用于对所述统一矩阵进行求解。
根据本发明所提供的方法或系统,可以在模拟时提供不同区块之间数据的实时交换,杜绝了区域之间流体的不连续性,还可提高了模拟的效率,增加了以后模拟器的扩展性。
附图说明
通过对结合附图所示出的实施方式进行详细说明,本发明的上述以及其他特征将更加明显,本发明附图中相同的标号表示相同或相似的部件。在附图中,
图1示出了适于用来实现本发明实施方式的示例性计算系统100的框图。
图2示出了油藏的区域划分的一个示例。
图3示出了根据本发明一实施例的用于进行油藏模拟的方法流程图。
图3A示出了根据本发明另一实施例的区域划分与矩阵构建的示意图。
图3B示出了根据本发明一实施例的网格划分以及界面的示意图。
图3C示出了根据本发明一实施例的通过中间媒介物质进行映射的示意图。
图3D示出了根据本发明另一实施例的通过中间媒介物质进行映射的示意图。
图3E示出了根据本发明另一实施例的在雅各比矩阵中如何体现界面处物质的交换。
图4示出了根据本发明一实施例的用于油藏模拟的系统的框图。
具体实施方式
在下文中,将参考附图通过实施方式对本发明提供的用于软件快速配置的方法和系统进行详细地描述。
所属技术领域的技术人员知道,本发明的多个方面可以体现为系统、方法或计算机程序产品。因此,本发明的多个方面可以具体实现为以下形式,即,可以是完全的硬件、完全的软件(包括固件、驻留软件、微代码等)、或者本文一般称为“电路”、“模块”或“系统”的软件部分与硬件部分的组合。此外,本发明的多个方面还可以采取体现在一个或多个计算机可读介质中的计算机程序产品的形式,该计算机可读介质中包含计算机可用的程序码。
可以使用一个或多个计算机可读的介质的任何组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电的、磁的、光的、电磁的、红外线的、或半导体的系统、装置、器件或任何以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括以下:有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任何合适的组合。在本文件的语境中,计算机可读存储介质可以是任何包含或存储程序的有形的介质,该程序被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可包括在基带中或者作为载波一部分传播的、其中体现计算机可读的程序码的传播的数据信号。这种传播的信号可以采用多种形式,包括——但不限于——电磁信号、光信号或任何以上合适的组合。计算机可读的信号介质可以是并非为计算机可读存储介质、但是能发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序的任何计算机可读介质。计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括——但不限于——无线、电线、光缆、RF等等,或者任何合适的上述组合。
计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括——但不限于——无线、电线、光缆、RF等等,或者任何合适的上述组合。
用于执行本发明的操作的计算机程序码,可以以一种或多种程序设计语言的任何组合来编写,所述程序设计语言包括面向对象的程序设计语言-诸如Java、Smalltalk、C++之类,还包括常规的过程式程序设计语言-诸如“C”程序设计语言或类似的程序设计语言。程序码可以完全地在用户的计算上执行、部分地在用户的计算机上执行、作为一个独立的软件包执行、部分在用户的计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在后一种情形中,远程计算机可以通过任何种类的网络——包括局域网(LAN)或广域网(WAN)-连接到用户的计算机,或者,可以(例如利用因特网服务提供商来通过因特网)连接到外部计算机。
以下参照按照本发明实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述本发明的多个方面。要明白的是,流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合,都可以由计算机程序指令实现。这些计算机程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器,从而生产出一种机器,使得通过计算机或其它可编程数据处理装置执行的这些指令,产生实现流程图和/或框图中的方框中规定的功能/操作的装置。
也可以把这些计算机程序指令存储在能指令计算机或其它可编程数据处理装置以特定方式工作的计算机可读介质中,这样,存储在计算机可读介质中的指令产生一个包括实现流程图和/或框图中的方框中规定的功能/操作的指令装置(instruction means)的制造品。
也可以把计算机程序指令加载到计算机或其它可编程数据处理装置上,使得在计算机或其它可编程数据处理装置上执行一系列操作步骤,以产生计算机实现的过程,从而在计算机或其它可编程装置上执行的指令就提供实现流程图和/或框图中的方框中规定的功能/操作的过程。
下面参看图1。图1示出了适于用来实现本发明实施方式的示例性计算系统100的框图。如所示,计算机系统100可以包括:CPU(中央处理单元)101、RAM(随机存取存储器)102、ROM(只读存储器)103、系统总线104、硬盘控制器105、键盘控制器106、串行接口控制器107、并行接口控制器108、显示控制器109、硬盘110、键盘111、串行外部设备112、并行外部设备113和显示器114。在这些设备中,与系统总线104耦合的有CPU 101、RAM 102、ROM 103、硬盘控制器105、键盘控制器106、串行控制器107、并行控制器108和显示控制器109。硬盘110与硬盘控制器105耦合,键盘111与键盘控制器106耦合,串行外部设备112与串行接口控制器107耦合,并行外部设备113与并行接口控制器108耦合,以及显示器114与显示控制器109耦合。应当理解,图1所述的结构框图仅仅为了示例的目的而示出的,而不是对本发明范围的限制。在某些情况下,可以根据具体情况而增加或者减少某些设备。
在石油开采过程中,通常先开始利用地层压力进行衰竭式开采,等压力降到一定程度之后,油藏开采率就会下降,这时就需要采用一些其他技术来提高采收率。例如通常采用的注水开采就是为了提高地层压力。其中混相驱的机理主要是对重油起到两个作用:一个是实现油膨胀;另一个是降低油的粘度。混相驱可以注入CO2、N2、聚全物,也可以注入干气来实现。
然而,对于一个已经开采过一段时间的油藏,油藏内部属性会变得比较复杂,如果整个油藏区域都采用组分,虽然可以进行模拟,但是计算机耗时较大,而且对于以后的扩展也是不利的;即使采用将整个油藏分开模拟,挥发油藏采用组分,黑油区域采用黑油模拟器,但是当整个油藏的弛豫时间和每个区域的弛豫时间在一个量级时,这种做法也是毫无精度可言的。
图2示出的是包括挥发油藏和黑油油藏的油藏区域示意图。其中的黑油油藏即常规油藏,即含有非挥发组分的黑油和挥发性组分的原油溶解气两个系统。其中的挥发油藏,类似于凝析气,但含有较多重质化合物组分,在地层条件下类似液体,它的相包络线更宽,且此类油藏温度接近临界温度,所以挥发油藏也称为近临界油藏。
接下来介绍图3,图3所示的方法从步骤310开始,示出了根据本发明一实施例的用于本申请的油藏模拟方法的流程图,可以实现组分和黑油油藏混合同时模拟的技术方案,并且,申请人在界面的处理方面还作出了进一步的改进。
在步骤310中,将油藏划分为多个区域,所述多个区域至少包括组分区域和黑油区域。
在一个实施例中,通过步骤310,可以用于模拟以组分变化为主要特征的油气田开发过程,同时兼容扩展型黑油模拟器的功能。从而在对复杂油藏进行模拟时,可以同时应用到黑油和组分混合油藏的优势。
黑油模型是模拟含有非挥发组分的黑油和挥发性组分的原油溶解气两个系统在油藏中流动的数值模型,主要用于各种常规油藏,特别是不发生反凝析现象的油藏模拟,通常用于油、气田开发方案的优选、地下剩余油分布状况的分析和老油田的开发调整等等。在一个具体的实施例中,应用黑油模型的黑油区域包括下述种类中的至少一种:黑油、挥发性油、干气、湿气等油藏。
而组分模型可模拟含有多个高挥发性碳氢化合物在油藏中流动的模型,主要用于凝析气藏、轻质油和挥发油藏的开发设计。在一个具体的实施例中,应用组分模型的组分区域包括下述种类中的至少一种:凝析气田、轻质油油田,挥发油藏,CO2注入区域,氮气注入区域以及混相区等。
在油藏描述和油藏模拟时,通常将油藏中各种碳氢化合物进行简单分组,以提高模拟的效率。在数值模拟过程中对流体的表征也会随着条件的不同而发生变化。举例而言,在衰竭式开采时,用两个组分就能达到很好的模拟效果;而在强化开采,例如气体回注开采、混相驱开采,通常需要5个甚至10个以上的组分。
需要说明的是,组分区域和黑油区域的数目均可以是至少一个或多个,本领域技术人员可以根据具体的情况设定不同的数目。
在步骤320中,对所述多个区域中的每一个区域分别进行网格划分,并在所述多个区域之间形成界面。
图3B示出了根据本发明一实施例的网格划分以及界面的示意图,在这个实施例中,左边Bmt1,…,Bmt4…代表黑油区域界面处的网格,右边Cmt1,…,Cmt4,…代表组分区域界面处的网格。而阴影部分则示出了黑油区域和组分区域的界面,可见,界面可以是两个区域相邻的网格部分。
在步骤330中,将所述组分区域以及所述黑油区域分别映射至中间媒介物质以建立关联。
图3C示出了根据本发明一实施例的通过中间媒介物质进行映射的示意图。在这个实施例中,为便于扩展,提供一组标准物质组成,如3C所示的实施例中,作为中间媒介物质的“标准物质”包括:独立组分、独立气相和独立油相,所属领域技术人员在本申请公开的基础上当然也可以结合公知技术对其进行扩展。例如,由于标准物质是作为接口使用的,初步采用相态作为划分标准,这里只举例气相和液相,实际上液体也可以有不同的相,也可以加入固相(例如高压下水的固态冰就有很多不同的相),用于模拟木屑、水泥等的注入情况。当然,这些扩展的实施例也应该在本申请的保护范围之内。
在该实施例中,黑油区域的物质可以通过特定对应方式映射至“标准物质”相应的部分。具体到图3C而言,黑油区域中的“水聚合物”对应着标准物质中的“独立组分”,“气溶解气”部分对应着标准物质中的“独立气相”,黑油区域中的“油”则对应着标准物质中的“独立油相”。
同样组分区域的物质也可以通过特定对应方式映射至“标准物质”相应的部分。具体到图3C而言,组分物质中的“水”、“气”以及“轻组分,中间组分,重组分”则分别对应着标准物质中的“独立组分”、“独立气相”以及“独立油相”。
通过上述方式,“标准物质”作为中间媒介物质起到“接口”的作用,从而将黑油区域和组分区域的物质联系起来。
图3D示出了根据本发明另一实施例的通过中间媒介物质进行映射的示意图。在这个实施例中,可直接取黑油区域的物质作为中间媒介物质。如图3D所示的实施例中,黑鱼区域B所包括的“水”、“气”以及“油”,等物质作为中间媒介物质。进而,将组分区域C中的物质直接通过映射与黑油区域所包括的物质进行联系,具体而言,是将组分区域C中所包括的“水”映射到该中间媒介物质中的“水”、将组分区域C中所包括的“气”映射到该中间媒介物质中的“气”,以及图4所示组分区域C中的剩余组分映射到该中间媒介物质中的“油”,在图4所示的实施例中,剩余组分包括轻组分(例如C1以及N2),中间组分(例如,CO2,H2S,C2-C6),以及重组分(例如C7)。
在步骤340中,在所述油藏动态流动时,根据物质守恒定律对所述界面进行处理。
本申请所提供的技术不同于组分中“赝组分”的划分,这里所要处理的是在至少包括组分区域和黑油区域的混合油藏模拟时,可通过物质守恒定律对界面进行处理。
接下来,申请人将对本申请中在界面处如何体现了物质守恒定律进行详细的说明。在油藏动态流动时,如果有组分网格Cinti的物质减小(去掉源汇项的影响),则根据物质守恒,则必有物质进入黑油区域,通过映射关系黑油区域相应的物质增加;相反,如果有黑油网格Bmt1的物质减小(除去源汇项),则根据物质守恒,则必有物质进入组分区域,通过映射关系组分区域相应的物质增加。
在步骤350中,根据所述多个区域以及所述界面处理结果构建统一矩阵。如下图所示右边为生成的雅各比矩阵。
在本申请请的一个实施例中,在构建矩阵时,可将各区域网格中相应的物理量按同种类型排序。
图3A示出了根据本发明另一实施例的区域划分与矩阵构建的示意图,如图3A所示,将黑油区域的油藏放在矩阵的左上角,即将黑油区域的网格排序在前面;将其与组分类型区域如图中所示B组分、C和D区域的网格依次排列,进而生成该系数矩阵。
由于所生成的系数矩阵大多数为0元素,所以需要对该系数矩阵进行特殊的处理,例如D4排序、双循环排序等,以使的所生成的系数矩阵便于求解。
在油藏的模拟过程中,如果发生了油藏的动态流动,则在所生成的雅各比矩阵的非对角元中体现了界面处物质的交换。
图3E示出了根据本发明另一实施例的在雅各比矩阵中如何体现界面处物质的交换。将网格排序编号之后(例如整个油藏有N个网格),将边界处的网格划分给黑油区域,如图3E所示(设网格数为m,界面处为k个网格,图中k=4),所以雅各比矩阵元1到(m+k)为黑油部分的矩阵元素;同样对于组分部分在网格排列最后加k网格,即雅各比矩阵元N+1到N+k。矩阵中m+k和N+k的交叉项目,即非对角元素,体现了界面处物质的交换。
在步骤360中,对上述构建的统一矩阵进行求解。
然而,方程组求解技术并不是本申请发明的重点,所属领域技术人员具体可采用如下方法之一:
其一,线性方程组求解。包括高斯(Gauss)消元法、简单迭代(Jacobi)法、松弛迭代法、共轭梯度法。
其二,非线性方程组求解。包括梯度法、拟牛顿法等。
其三,多重网格算法。多重网格算法是一种迭代方法,其原理是:对于方程(组)求解时,误差的傅里叶变换之后各级的衰竭程度不同。高频项是局域网格之间的相互耦合,是局部行为;而低频项来源于边界条件,是全局行为。传统迭代方法,块松弛都是对局域的处理,对于全局误差衰减较慢,误差趋于光滑。将边界信息传递至计算点需要N量级的迭代,收敛速度非常慢。所以采用将网格方程的残值限定于粗网络上。在粗网络上精确求解之后,在将解延拓至细网格上,与原近似解组合得到网格新的近视解,这就是粗网格迭代矫正。粗网格格点少,边界条件能较快传递至计算点,但是同样纯在高频与低频的误差,可以通过进一步粗化网络,直到计算结果达到所需精度。
由于具体的矩阵求解方法是本领域技术人员熟知的技术,申请人对此不加赘述。
通过构建以及求解的统一的矩阵,能够实现组分和黑油的线性区域分解模拟,就能够得到模拟油藏方程中的变量结果,从而可以实现并测试能够求解黑油与组分的混合油藏模型。
基于以上公开,本领域技术人员明了,图3所示出的油藏数值模拟方法的应用范围至少包括如下各种具体的情形:
1.挥发油和凝析气自然衰竭开采;
2.凝析气藏的循环注干气计算;
3.挥发油和注气开发;
4.通过注贫气对残余油再蒸发开采;
5.凝析气田的注N2气开发;
6.油藏挥发油开发;
7.水驱油藏的混相驱开发;
8.黑油及天然气的开发模拟。
图4示出了根据本发明一实施例的用于油藏模拟的系统的框图。图4所示的系统包括:
划分装置410,用于将油藏划分为多个区域,所述多个区域至少包括组分区域和黑油区域。
区域划分装置420,用于对每一个区域进行网格划分,并在区域之间形成界面;
映射装置430,用于将所述组分区域以及所述黑油区域分别映射至中间媒介物质以建立关联;
界面处理装置440,用于根据物质守恒定律,对所述界面进行界面处理;
矩阵构建装置450,根据所述多个区域的映射结果以及所述界面处理结果构建统一矩阵;
求解装置460,用于对所述统一矩阵进行求解。本领域技术人员应当理解,图4所示的系统中的装置410-430分别对应于图3所示的方法中的步骤310-360,在此不再赘述。
采用本申请所提供的是一种可以实现组分和黑油油藏混合同时模拟的技术方法,增加了以后模拟器的扩展性。可以实现黑油模型和组分模型可以整合在一起同时模拟,即设计的算法能够实现组分和黑油的线性区域分解模拟,初步实现并测试能够求解黑油与组分的混合油藏模型。
在本申请中,特别提供的是针对实现组分和黑油油藏混合同时模拟的算法,所需要的界面处理技术,例如,提供了两种具体的实施例:标准物质接口技术和直接映射黑油方法。在一个改进的实施例中,如果中间媒介接口应用的是标准物质,就可以增加油藏模拟器的扩展性。
附图中的流程图和框图,图示了按照本发明各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分,所述模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个接连地表示的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或操作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
虽然以上结合具体实例,对本发明的利用远程应用处理本地文件的系统及方法进行了详细描述,但本发明并不限于此。本领域普通技术人员能够在说明书教导之下对本发明进行多种变换、替换和修改而不偏离本发明的精神和范围。应该理解,所有这样的变化、替换、修改仍然落入本发明的保护范围之内。本发明的保护范围由所附权利要求来限定。
Claims (10)
1.一种油藏模拟方法,其特征在于,所述方法包括:
将油藏划分为多个区域,所述多个区域至少包括组分区域和黑油区域;
对每一个区域进行网格划分,并在区域之间形成界面;
将所述组分区域以及所述黑油区域分别映射至中间媒介物质以建立关联;
在所述油藏动态流动时,根据物质守恒定律对所述界面进行处理;
根据所述多个区域以及所述界面的映射处理结果构建统一矩阵;
对所述统一矩阵进行求解。
2.根据权利要求1所述的油藏模拟方法,其特征在于,其中,所述中间媒介物质选择以下两种方式之一:
标准物质以及黑油区域所包括的物质。
3.根据权利要求2所述的油藏模拟方法,其特征在于,其中,
所述标准物质包括:独立组分、独立气相以及独立油相;
将所述组分区域以及所述黑油区域中所包括的物质分别映射至所述标准物质中的各个部分。
4.根据权利要求2所述的油藏模拟方法,其特征在于,其中,
将黑油区域所包括的多个物质作为所述中间媒介物质;
将组分区域的所包括的物质分别映射至所述黑油区域所包括的多个物质。
5.根据权利要求1所述的油藏模拟方法,其特征在于,其中,
在所述统一矩阵的非对角元中体现所述界面处物质的交换。
6.一种油藏模拟系统,其特征在于,所述系统包括:
区域划分装置,用于将油藏划分为多个区域,所述多个区域至少包括组分区域和黑油区域;
网格划分装置,用于对每一个区域进行网格划分,并在区域之间形成界面;
映射装置,用于将所述组分区域以及所述黑油区域分别映射至中间媒介物质以建立关联;
界面处理装置,用于根据物质守恒定律,对所述界面进行界面处理;
矩阵构建装置,根据所述多个区域的映射结果以及所述界面处理结果构建统一矩阵;
矩阵求解装置,用于对所述统一矩阵进行求解。
7.根据权利要求6所述的油藏模拟系统,其特征在于,其中,所述中间媒介物质选择以下两种方式之一:
标准物质以及黑油区域所包括的物质。
8.根据权利要求7所述的油藏模拟系统,其特征在于,其中,
所述标准物质包括:独立组分、独立气相以及独立油相;
将所述组分区域以及所述黑油区域中所包括的物质分别映射至所述标准物质中的各个部分。
9.根据权利要求7所述的油藏模拟系统,其特征在于,其中,
将黑油区域所包括的多个物质作为所述中间媒介物质;
将组分区域的所包括的物质分别映射至所述黑油区域所包括的多个物质。
10.根据权利要求6所述的油藏模拟系统,其特征在于,其中,
在所述统一矩阵的非对角元中体现所述界面处物质的交换。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210575346.5A CN103902751B (zh) | 2012-12-27 | 2012-12-27 | 油藏数值模拟的边界处理方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210575346.5A CN103902751B (zh) | 2012-12-27 | 2012-12-27 | 油藏数值模拟的边界处理方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103902751A CN103902751A (zh) | 2014-07-02 |
CN103902751B true CN103902751B (zh) | 2017-11-14 |
Family
ID=50994073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210575346.5A Expired - Fee Related CN103902751B (zh) | 2012-12-27 | 2012-12-27 | 油藏数值模拟的边界处理方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103902751B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104912538B (zh) * | 2015-04-16 | 2018-01-05 | 中国石油天然气股份有限公司 | 一种油田区块边界流量确定方法及装置 |
CN107463754B (zh) * | 2017-08-16 | 2020-06-02 | 中国海洋石油集团有限公司 | 用于flng装置液化工艺的海上边界条件仿真模拟方法 |
CN108959767B (zh) * | 2018-07-02 | 2020-05-22 | 中国地质大学(北京) | 一种窄河道型气藏不同井型凝析油伤害数值模拟方法 |
CN109882164B (zh) * | 2019-03-28 | 2022-11-25 | 中国石油大学(华东) | 一种裂缝性碳酸盐岩油藏的大尺度酸化模拟方法 |
CN110821453B (zh) * | 2019-11-07 | 2021-11-23 | 成都北方石油勘探开发技术有限公司 | 基于三维地质数值模型的注气油藏开发方案设计方法 |
US11754745B2 (en) | 2020-06-30 | 2023-09-12 | Saudi Arabian Oil Company | Methods and systems for flow-based coarsening of reservoir grid models |
CN117166996B (zh) * | 2023-07-27 | 2024-03-22 | 中国地质大学(北京) | 地质参数门槛值的确定方法、装置、设备及存储介质 |
-
2012
- 2012-12-27 CN CN201210575346.5A patent/CN103902751B/zh not_active Expired - Fee Related
Non-Patent Citations (5)
Title |
---|
基于预处理GMREs算法的油藏数值模拟研究;刘刚等;《天然气与石油》;20111031;第29卷(第5期);全文 * |
离散裂缝性油藏数值模拟方法;张允等;《大庆石油学院学报》;20100630;第34卷(第3期);全文 * |
缝洞型碳酸盐岩油藏三维三相数值模拟新方法;康志江等;《特种油气藏》;20100630;第17卷(第3期);全文 * |
缝洞型碳酸盐岩油藏数值模拟新方法;康志江;《大庆石油地质与开发》;20100228;第29卷(第1期);摘要和正文第1节 * |
非均质油藏多尺度混合有限元数值模拟方法;姚军等;《石油学报》;20120531;第33卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103902751A (zh) | 2014-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103902751B (zh) | 油藏数值模拟的边界处理方法及装置 | |
Flemisch et al. | Benchmarks for single-phase flow in fractured porous media | |
Du et al. | Field study: embedded discrete fracture modeling with artificial intelligence in Permian basin for shale formation | |
US10961826B2 (en) | Subsurface modeler workflow and tool | |
Moncorgé et al. | Modified sequential fully implicit scheme for compositional flow simulation | |
Beckner et al. | EMpower: new reservoir simulation system | |
Yang et al. | High resolution hydraulic fracture network modeling using flexible dual porosity dual permeability framework | |
EP2588894A2 (en) | Reservoir simulation of giant subsurface reservoirs with accelerated phase equilibrium determination | |
MX2011005108A (es) | Sistemas y metodos para desarrollar dinamicamente planos de pozo de perforacion con un simulador de deposito. | |
CN108153933A (zh) | 一种致密储层生产预测方法及装置 | |
Gange et al. | Optimal k-level planarization and crossing minimization | |
Sankaran et al. | A flow network model based on time of flight for reservoir management | |
CN108564264B (zh) | 一种确定注水开发效果的数据处理方法及装置 | |
CN103902752B (zh) | 一种新型黑油、组分一体化混合数值模拟方法及装置 | |
Edwards et al. | Representing hydraulic fractures using a multilateral, multisegment well in simulation models | |
Liao et al. | River environmental decision support system development for Suzhou Creek in Shanghai | |
Kang et al. | Development of improved 4D CAD system for horizontal works in civil engineering projects | |
Shahvali et al. | Efficient coupling for nonlinear multiphase flow with strong gravity | |
Ahmed et al. | Cracking the code of unconventional reservoir simulation | |
Anwar et al. | An effective streamline history matching technique | |
Cominelli et al. | Deployment of high-resolution reservoir simulator: Methodology & cases | |
Hammadi et al. | Unlocking the Production Potential of Brown Fields through Gas Lift Optimization GLO | |
Szyndel et al. | Implementing A Hardware Agnostic Commercial Black-oil Reservoir Simulator | |
Bahaa et al. | Soft computation application: utilizing artificial neural network to predict the fluid rate and bottom hole flowing pressure for gas-lifted oil wells | |
Ekkawong et al. | SMARTDRAIN: an intelligent workflow for well-placement optimisation by integrating streamlines time-of-flight and exhaustive search algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171114 Termination date: 20211227 |