CN103329225A - 地层中的不相融的两相流体的塑造 - Google Patents
地层中的不相融的两相流体的塑造 Download PDFInfo
- Publication number
- CN103329225A CN103329225A CN2011800619285A CN201180061928A CN103329225A CN 103329225 A CN103329225 A CN 103329225A CN 2011800619285 A CN2011800619285 A CN 2011800619285A CN 201180061928 A CN201180061928 A CN 201180061928A CN 103329225 A CN103329225 A CN 103329225A
- Authority
- CN
- China
- Prior art keywords
- increment
- fluid
- sample
- measure
- length
- 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
Images
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
Abstract
当将注水前缘注入如地下含有石油的岩层或岩心复合物的多孔媒介片段中时,其传播是在多个离散的时间步长中作为时间函数测量的。一种模型是通过沿着多孔媒介片段的长度以多个时间步长测量水饱和度曲线而形成的。该模型实际上将多孔媒介片段划分为多个距离相等的单独片段或分系统。每个分系统的饱和度是根据流体注入量、预定的分流、初始平均饱和度确定的。
Description
技术领域
本发明涉及地下储层的计算机化分析,特别是构造用于岩心样本可渗透性测试和油藏模拟的两相不相融流体的模型。
背景技术
已有的传统做法是在油气藏或储层开采期中的某些时间通过恢复技术来提高产量。这种技术中包括注水。水与油是不溶合的,因为它们不能彼此混合或彼此发生化学反应。存在于储层中的不同液体(石油、燃气、水)通过岩砂层的流速通常也会不同。
在储层的开采过程期中,通常的做法是对通过储层的液体流构造模型或模拟。这样做是为了通过储层准确地评估和分析潜在的或历史上的产量。
在构造储层液体流的模型或模拟时,必须要考虑液体不相融的特性。被称为Buckley Leverett模型的模型被广泛地使用了许多年。该技术最初是在《“Mechanism of Fluid Displacement in Sands”,S.E.Buckley and M.C.Leverett,Trans.AIME(1942),Vol.145,p.107-116》中。在后来的若干年中,使用该方法的作品中存在一些问题。一个具体的问题是,利用Buckley Leverett方法产生的地层流体饱和度数值对同一物理位置指示出多个流体饱和度数值,这通过定义是不会发生的。
发明内容
简单地说,本发明提供了一种新的和改进的计算机实施的用于获取地球岩层的多孔媒介片段对流体注入量的饱和度的测量方法。将多孔媒介片段的系统样本长度分割成多个样本长度增额,测量在选定时间增额中向样本长度增额中注入的流体注入量。随后,测量在选定的时间增额中在样本长度增额中由注入流体产生的液体分流量,并测量在选定的时间增额中注入流体在所述样本长度增额中的流体饱和度。对选定时间增额中注入流体在所述样本长度增额中的流体饱和度的测量进行记录,并测量选定时间增额中未在样本长度增额中达到饱和的剩余的流体容量。
本发明还提供了一种新的和改进的用于获取地球岩层的多孔媒介片段对流体注入量的饱和度的测量的数据处理系统。该数据处理系统包括数据存储器和处理器,该处理器执行步骤:将多孔媒介片段的系统样本长度分割成多个样本长度增额,测量选定时间增额中向样本长度增额中注入的流体的注入量。处理器还测量选定时间增额中在样本长度增额中由注入流体产生的分流曲线,并测量选定时间增额中注入流体在所述样本长度增额中的流体饱和度。处理器还对选定时间增额中注入流体在所述样本长度增额中的流体饱和度的测量进行记录,并测量选定时间增额中未在样本长度增额中达到饱和的剩余的流体容量。
本发明还提供了一种新的和改进的数据存储设备,其具有存储在计算机可读媒介中的计算机可操作指令,用于使数据处理系统对地球岩层的多孔媒介片段对流体注入量的饱和度进行测量,存储在数据存储设备中的指令使数据处理系统将多孔媒介片段的系统样本的长度分割成多个样本长度增额,并测量选定时间增额中向样本长度增额中注入的流体的注入量。存储在数据存储设备中的指令还包括用于测量选定时间增额中在样本长度增额中由注入流体产生的分流曲线,并测量选定时间增额中注入流体在所述样本长度增额中的流体饱和度的指令。存储在数据存储设备中的指令还包括用于对选定时间增额中注入流体在所述样本长度增额中的流体饱和度的测量进行记录,并测量选定时间增额中未在样本长度增额中达到饱和的剩余流体容量的指令。
附图说明
图1为作为水饱和度的函数的分流曲线的测量的图形显示。
图2为作为无因次距离的函数的注水饱和度的测量的图形显示,这是利用现有技术Buckley Leverett模型在不实施任何校正的情况下,通过使用用于显示图1的数据组而形成的。
图3为作为无因次距离的函数的激震前缘水饱和度曲线的测量的图形显示,这是使用利用平均水饱和度校正的现有技术BuckleyLeverett模型,通过使用用于显示图1的数据组形成的。
图4为根据本发明的用于为地表下结构的流体流动建模的计算机系统的示意图。
图5为根据本发明的在构造地表下结构流体流动模型时在图4中执行的数据组处理步骤的功能框图。
图6为作为水饱和度的函数的注入流体分流曲线的典型实例的合成图形显示。
图7为根据本发明对于若干孔隙容量(PV)比率的作为无因次距离函数的水饱和度曲线测量的图形显示,这是通过使用用于显示图6的数据组形成的。
图8为根据本发明实施平滑技术之前和之后的作为无因次距离函数的水饱和度曲线测量的图形显示,这是通过使用用于显示图6的数据组形成的。
图9是通过合成数据根据用于显示图6的数据组形成的饱和度曲线测量的绘图与根据现有技术Buckley Leverett方法的绘图对比的图形显示。
图10为作为饱和度的函数的合成分流曲线的图形显示。
图11为作为饱和度的函数的合成分流曲线的图形显示。
具体实施方式
首先,提供对两相流体流动的物理性质和关系的解释。已知为Buckley Leverett模型的模型是基于模型的某些物理条件的存在而获得的。流体位移是一维的,条件是平衡的。保持流体压力,流体是不相融的。认为重力和毛细压力被认为是可忽略的,流体是不能压缩的。图1为作为饱和度的函数的合成分流量曲线的图形显示,这通常是通过对岩层矿石的岩心样本进行实验室实验产生的。输入数据用于对现有技术Buckley Leverett模型方法给出一个理想的输出曲线。
对于水取代石油时的替代过程,在岩心栓或储层中任意点处的水分流量定义为:
其中, 以及
假定水和石油中的压力梯度相同并忽略毛细压力影响,上面的等式变为:
对于某段时间周期内的长度的控制量应用水质量守恒,可将该质量守恒写为:
假定水是不能压缩的,上面的等式变为:
如果Δx→0并且Δt→0并且将等式(1)中的分流项代入等式(3)中的流速,则传统的已知的Buckley Leverett等式模型为:
结合等式(4)得到下面的结构,其描述了注水前缘增量:
为了绘制注水前缘,可通过分流曲线计算f’w或该分流曲线是使用等式(2)的相对渗透率并随后反代入等式(6)中的数值生成的。图1是作为水饱和度Sw的函数的分流量曲线fw及其导数f’w的示例曲线。
但是,在原始的Buckley Leverett模型中,如图2所示,计算出的水饱和度曲线在任何距离具有3个饱和度值,即Sw1,、Sw2、Sw3。Buckley Leverett模型被改进,并且引入了激波前缘饱和度来增加图2中描绘的原始模型的现实意义。激波前沿前的原生的水饱和度曲线以及从Buckley Leverett等式得出的大部分饱和度曲线被消除并由激波前沿(图3)代替。前沿的数学解随后可以通过其它人利用平均水饱和度的概念获得。
如图2所显示,Buckley Leverett模型在沿着距离绘标的每个点上提供了多个饱和度,这在物理上是不可能的。其它人提出Buckley Leverett模型的这个问题存在于相对渗透率函数中。
Buckley Leverett模型是平衡条件下系统质量守恒的体现。该模型指示出特定时间间隔内驱替流体的累积,饱和度的变化等于进入系统的替换流体容量与离开系统的流体容量的差,如等式(4)中所示。这表示f’w可描述为:
其中qw为无因次形式
根据本发明,已经确定上述的注水前缘增加计算误差是由于模型没有被正确地实施。用于计算前缘的f’w(等式7)与从BuckleyLeverett模型获得的物理对象不同。图7的f’w是通过在实验室测试中的相对渗透率实验中测量到的数据生成的,其未考虑生成分流曲线(图1)中的入口注入量。在数学上,图7的f’w应表示为:
其中qw为无因次形式
于是,可以看到,等式(8)中的f’w与等式(7)中的f’w是不同的。前者说明系统出口处的速度的变化,而后者说明系统入口与出口之间的速度差。等式(7)中的f’w还违背了Buckley Leverett的平衡假设,因为入口与出口处的速度不应随时间变化。使用不正确的f’w时该方案的物理含义是系统中某段时间间隔内驱替流体的累积量等于所产生的流体体积上的变化,这在物理上不会发生。
于是可以证明,由于物理含义中的矛盾,f’w的数值不能通过来自相对渗透率实验的分离曲线直接采集并应用到Buckley Leverett模型。本发明提供了一种利用新的和改良的方法来塑造地球岩层中的流体前缘饱和度曲线的模型,其中f’w可以直接在该模型中使用,而不存在任何矛盾。
根据本发明,可将等式(6)的正确形式表示为:
对于具有注入点和产出点并且t0=0的驱水系统,等式(9)可重写为:
等式(10)
因为分子和分母代表相同的系统,因此因数应当代表注入系统中的水的无因次孔隙体积:
通过将等式(11)代入等式(10),等式变为:
为了追踪到作为已知分流曲线的注入系统的前缘的正向传播,可将系统划分为多个固定Δx的分系统。分流曲线会被获知并且将根据平均水饱和度描绘分流曲线。
等式中每个Δx的未知因数为(fw)p和【(S】w)Δt。注水比(fw)i和注入前的初始水饱和度〖(S〗wi)为便于测量的固定因数。注入的孔隙体积倍数〖(PV〗i)为可变因数,其为时间函数并且可利用等式(11)获得。只有两个未知数,(fw)p和【(S】w)Δt。未知数的值可通过利用分流曲线找到满足等式的适当值来获得。
同样的技术可用于注水前缘的反向追踪。这时,(fw)p和【(S】w)Δt为固定的已知参数,而(fw)i和〖(S〗wi)为未知数并需要使用分流曲线求解。本发明基于多孔媒介片段或样本,将上面的计算机化分析用于地下储层或岩层中的流体流动模型构造中。
如图4中所示,根据本发明的数据处理系统D包括计算机40,计算机40具有处理器42和耦接到处理器42的存储器44,存储器44用于存储操作指令、控制信息和数据库。如果需要,计算机40可以是便携式数字处理器,如膝上计算机、笔记本电脑形式的个人计算机,或者其它适当的已编程或可编程的数字数据处理装置,如台式电脑。同时,需要理解的是计算机40可以是具有多个结点的多芯片处理器,如来自Intel Corporation或Advanced Micro Devices(AMD)的处理器,或者是具有适当处理能力的任何普通类型的大型计算机,如来自Armonk,N.Y.的International Business Machines(IBM)或其它来源的可用计算机
计算机40具有用户接口46和用于显示输出数据或井测数据的处理记录的输出显示器48,该井测数据是根据本发明为了获得地下结构中流体的可传递性而进行测定的。输出显示器48包括如打印机及输出显示屏等设备,它们能够提供图形、数据表、图形图像、数据图等形式的打印的输出信息或可视化显示作为输出记录或图像。
计算机40的用户接口46还包括适当的用户输入设备或输入/输出控制单元50,用于为用户提供控制或访问信息和数据记录并操作计算机40的通道。数据处理系统D还包括存储于计算机存储器中的数据库52,它可以是内部的存储器44,或外部的、联网的或不联网的存储器,如位于相关联的数据库服务器56中以54标示的存储器。
数据处理系统D包括存储在计算机40的存储器44中的程序代码60。根据本发明,程序代码60为计算机可操作指令形式,用于使数据处理器42获取地下结构中流体的可传递性的测量,这将在下面进行详尽地解释。
应当注意到程序代码60可以是微代码、程序、例行程序的形式,或者是符号式的计算机可操作语言,该语言提供一套详细的排列好的控制数据处理系统D运行并引导其操作的操作指令。程序代码60可存储于计算机40的存储器44中,或者存储在计算机磁盘、磁带、传统的硬盘驱动器、电子只读存储器、光存储设备或其它的存储了计算机可用媒介的适当的数据存储设备上。如图所示,程序代码60还可包含于如以服务器64作为计算机可读媒介的数据存储设备上。
这里的图5的流程图F示出了通过计算机程序软件实现的本发明的逻辑结构。本领域内的技术人员应理解的是,该流程图示出了根据本发明运行的计算机程序代码元件的结构。本发明通过使用了程序代码指令的计算机组件实现其基本的实施例,这些程序代码指令指示数字数据处理系统D执行与流程图F中所示的相对应的处理步骤序列。
现参照图5,流程图F为高阶逻辑流程图,其示出了根据本发明的构造测量地下结构中流体的可传递性的方法。本发明的方法在计算机40中可利用存储在存储器44中的计算机程序步骤实施并可由计算机40的系统处理器42执行。处理系统D的输入数据为实验室或其它数据,包括初始水饱和度值、系统长度、孔隙率、注入量及比例数据、以及关于分流曲线(或岩层样本针对石油和水的相对渗透率)的数据。
如图5的流程图所示,示意性地示出了计算机执行方法或进程的步骤的优选顺序,该方法或进程用于获取地层岩石的多孔媒介片段对于流体注入量的饱和度的测量。
对于满足Buckley Leverett条件并且具有将要注入的用于替换另一种流体的流体(如,水)的多孔媒介片段或系统,如石油,可以下面的关系描述流体:
其中:
n:片段增额中的分系统或长度增额数量,其在注入点等于1
t:注入的时间步长,其在注入前等于0
Wn:在步骤t时注入分系统n的流体量
fi:注入流体的分流
【(f】p)n:产出流体的分流
(St)n:增额或分系统n的注入流体的饱和度
(St-1)n:前一时间步长中增额或分系统n的注入流体的饱和度
Q:注入流体的总量
通过应用下面所描述的并图5的进程序列中示意性示出的方法,水饱和度Sw可确定为时间函数以及片段中的一维空间,这在。在步骤100,多孔媒介片段或样本的长度在计算机数据中被划分为j个更小的长度相等的分系统,总注入量的计算机数据被分配为多个更小的容量。当满足下面的条件时,这些注入容量的离散化可代表某个时间步长中的注入量:
注水率(fi)、注入前的原始饱和度〖(S〗t=0)、注入量【(W】t)n=1为可通过实验测量的已知因数。如步骤102所示,这些初始参数是用于进一步处理的输入数据。
在步骤104,设置第一时间步长的初始计数,用于在注入点以第一长度增额执行进程,其中n=1,t=1。
在步骤106,对于长度增额为n时的产出流体(fp)、注入流体的饱和度〖(S〗t)n的创建,可通过利用预先确定的分流曲线(图6)来找到满足等式13的适当的(fp)和〖(S〗t)n的数值。这可通过多种方式来完成,如通过使用传统的计算机数值求解方法(如牛顿法)或通过其它计算机化的最优化或迭代尝试误差法。在步骤108,将确定的当前长度增额n下的分流量和注入流体的饱和度数值存入存储器。
在步骤110,当前长度增额n下的(fp)和〖(Si〗t)n的数值被用于材料平衡的计算,从而得到能够注入相邻系统的剩余容量,这是通过应用下面的等式进行的:Wn+1=Wn(fp)n。
在步骤112,确定在相邻时间步长(Wn+1)下注入的体积是否等于零。如果是,这表示仍然存在一些流体将要流入下一个相邻的长度增额n+1。在这种情况下,在步骤114,增加长度n、对相邻的长度增额建立(fp)和〖(S〗t)n、返回步骤106继续执行进程。
如果分系统数量n等于j,这表示已经在指定的时间步长对所有的分系统测量了饱和度。在步骤112,如果指定相邻时间步长(Wn+1)中的注入量等于0,则在指定时间步长t的总注入量已经包括之前的长度增额,不再留有可动流体进入下面的相邻长度增额。通过描绘作为距离的函数的长度增额1至n的〖(S〗t)n,可获得全部样本在时间步长t的注水前缘饱和度曲线。
在步骤116,确定长度增额中的累积注入量是否等于注入片段中的总注入量。如果是这种情况,则终止进一步的处理,饱和度曲线如步骤118中所描述的进行描绘。于是,描绘的输出显示代表作为时间和一维空间函数的注水前缘。如果在步骤116中的累积流体注入量尚未等于总注入量,则在步骤120对时间间隔计数器t进行增加。
为了确定下一个时间步长的饱和度曲线,〖(W〗t)等于本次时间步长中的注入量,[(S]t-1)n等于上次时间步长的[(S]t)n。
使用当前时间步长中的注水量针对下一时间步长执行进程,并且进程返回步骤106用于继续数据值确认。
图6示出了根据本发明的用于对注水前缘进行正向追踪的输入显示的实例。其通常是通过对岩层矿石的岩心样本进行实验而得到。在该实例中,数据源于定态岩心驱水实验。选择样本长度为Δx并具有无因次长度。执行进程来查看某个与整个样本具有相同岩石物理性质的片段的注水前缘性能。对不同数量的孔隙容量执行该进程,并追踪前缘增加,直至饱和度达到样本最初的原生水饱和度。与传统的Buckley Leverett前缘模型不同,图7中所绘制的每个前缘的解都是唯一的,并且不会产生多重值。显然前缘中出现了激震前缘现象,但是不必为了与实际匹配而强制在绘制中体现该现象。
图7的注水前缘曲线标示出,孔隙容量或PV比是根据最大孔隙容量决定的。PV比是基于具有1单位的无因次距离的片段的孔隙容量决定的。
图7的注水前缘曲线显示为每个连续长度增额的实际数值,但是也可以与实际计算增额值进行平滑的比较。原始的BuckleyLeverett模型似乎假设相关的分流为平均水饱和度,而不是任意点实际饱和度。
相比之下,本发明构造了一种基于非常小的Δx长度增额的样本的实际饱和度的水饱和度模型。在图7的平面图中,Δx增额可随意地选择为岩心样本的长度并且这些数值仅用于代表Δx的中间阶段以平滑曲线。图8上展示了原始曲线形状与平滑后的曲线形状的对比。应当注意的是,所标注的在数据图中每个连续长度下的确定模型饱和度曲线数值之间的差异可通过选择极小的Δx长度增额来避免。
图9展示了由Buckley Leverett与本发明的提出的饱和度曲线之间的另一个有益的对比。在两个前缘计算方案中使用了相同的注入量(无因次量=0.61)。这里的无因次距离代表样本的长度。
两种曲线之间的一个简单的可视化对比显露出某些情况。与用于前缘移动计算的原始量对比,图9的Buckley Leverett前缘模型未展示出相同的注水量。前缘曲线与上面的原始饱和度线条应当代表无因次注水量。图9的Buckley Leverett前缘模型下的区域展示出注入量为1.59,与用于向模型的注入量0.61不相等。这明确表示Buckley Leverett前缘模型违反了物质平衡法则。本发明的方法所描绘前缘展示出诸如量为0.61,这与用于前缘移动计算的量相同。
图9的Buckley Leverett前缘模型展示了在距离等于1的点的曲折。这是由于前缘对衍生出的分流曲线非常敏感,而根据本发明的前缘模型不存在该问题。
传统的Buckley Leverett前缘模型和根据本发明构造的前缘模型也对非常适合Buckley Leverett前缘模型的合成数据集进行检验。本文中的数据适合性涉及传统的单调的衍生形状,因为这会降低许多传统的Buckley Leverett前缘模型中的误差。图10展示了在样本中注入时的合成分流数据与饱和度前缘。图11中展示了两个模型之间的对比,其中本发明具有发展良好的更平滑的符合实际的激震前缘(曲线的右侧),而Buckley Leverett模型展示出了陡峭的激震前缘,其表现为直线,这是由于对Buckley Leverett模型的Weldge修改造成的产物。该产物在现有技术中是众所周知的,但是没有得到平滑地塑造,除了本发明。
图11输出的模型之间的另一个不太重要的差别在于曲线的左侧。与Buckley Leverett模型相比,本发明展示出了更符合实际的Sor的评估。这是由于Buckley Leverett模型直接将起始的一些点设置到Sor,而本发明在更慢和渐进的速度下分派高油饱和度值。新发明与实际更匹配,这是由于从Buckley Leverett模型的描述中达到Sor是不容易的过程。
已经对本发明进行了充分的描述,从而这方面具有普通知识的人能够再现和获得本发明所提到的结果。尽管如此,本领技术域内的技术人员可对本发明的主体实施一些这里未要求的修改,为了对这些修改提供确定的结构,或者相同的制造工艺,下面的权利要求中所要求的内容是需要,这些结构将覆盖在本发明的范围内。
应当注意和了解到,在不脱离如所附权利要求所述的本发明的精神和范围的情况下,可对上面所详细描述的本发明进行改进和修改。
Claims (24)
1.一种计算机执行的用于获取地球岩层的多孔媒介片段对流体注入量的饱和度测量的方法,其包括步骤:
将多孔媒介片段的系统样本长度分割成多个样本长度增额;
测量选定时间增额中向样本长度增额中注入的流体的注入量;
测量选定时间增额中在样本长度增额中由注入流体产生的分流曲线;
测量选定时间增额中注入流体在所述样本长度增额中的流体饱和度;
对选定时间增额中注入流体在所述样本长度增额中的流体饱和度的测量进行记录;
测量选定时间增额中未在样本长度增额中达到饱和的剩余流体容量。
2.如权利要求1所述的计算机执行的方法,还包括步骤:
确定所测量的剩余流体容量是否表示存在用于注入所述多孔媒介片段的相邻长度样本增额中的剩余流体容量;
如果是,重复如下步骤:测量注入流体的注入量,测量流体的分流曲线,测量流体饱和度,对测量的流体饱和度进行记录,测量选定时间增额中未在相邻样本长度增额中达到饱和的剩余流体容量;
如果不是,测量选定时间增额中通过注入流体带来的长度样本增额的饱和度曲线。
3.如权利要求2所述的计算机执行的方法,还包括步骤:在测量选定时间增额中通过注入流体带来的长度样本增额的饱和度曲线后,增加所选定的时间增额从而形成新的选定时间增额。
4.如权利要求3所述的计算机执行的方法,还包括步骤:
测量新的时间增额中向长度样本增额中注入的流体的注入量;
测量新的时间增额中由注入流体产生的长度样本增额中的分流曲线;
测量新的时间增额中注入流体在所述长度样本增额中的流体饱和度;
对新的时间增额中注入流体在所述长度样本增额中的流体饱和度的测量进行记录;
测量新的时间增额中未在长度样本增额中达到饱和的剩余流体容量。
5.如权利要求1所述的方法,其中注入的流体包括水。
6.如权利要求1所述的方法,其中多孔媒介片段包括岩心样本。
7.如权利要求1所述的方法,其中多孔媒介片段包括地下岩层片段。
8.如权利要求1所述的方法,还包括步骤:
对确定的注入流体的流体饱和度的测量进行输出显示。
9.一种用于测量地球岩层的多孔媒介片段对流体注入量的饱和度的数据处理系统,该数据处理系统包括:
数据存储器;
处理器,其用于执行如下步骤:
将多孔媒介片段的系统样本长度分割成多个样本长度增额;
测量选定时间增额中向样本长度增额中注入的流体的注入量;
测量选定时间增额中在样本长度增额中由注入流体产生的分流曲线;
测量选定时间增额中注入流体在所述样本长度增额中的流体饱和度;
在所述数据存储器中对选定时间增额中注入流体在所述样本长度增额中的流体饱和度的测量进行记录;
测量选定时间增额中未在样本长度增额中达到饱和的剩余流体容量。
10.如权利要求9所述的数据处理系统,其中处理器还执行如下步骤:
确定所测量的剩余流体容量是否表示存在用于注入多孔媒介片段的相邻长度样本增额中的剩余流体容量;
如果是,重复如下步骤:测量注入流体的注入量,测量流体的分流曲线,测量流体饱和度,对测量的流体饱和度进行记录,测量选定时间增额中未在相邻样本长度增额中达到饱和的剩余流体容量;
如果不是,测量选定时间增额中通过注入流体带来的长度样本增额的饱和度曲线。
11.如权利要求10所述的数据处理系统,其中处理器还执行如下步骤:
在测量选定时间增额中通过注入流体带来的长度样本增额的饱和度曲线后,增加所选定的时间增额从而形成新的选定时间增额。
12.如权利要求11所述的数据处理系统,其中处理器还执行如下步骤:
测量新的时间增额中向长度样本增额中注入的流体注入量;
测量新的时间增额中注入流体在长度样本增额中产生的分流曲线;
测量新的时间增额中注入流体在所述长度样本增额中的流体饱和度;
对新的时间增额中注入流体在所述长度样本增额中的流体饱和度的测量进行记录;
测量新的时间增额中未在长度样本增额中达到饱和的剩余流体容量。
13.如权利要求11所述的数据处理系统,还包括:
输出显示器,其用于对确定的注入流体的流体饱和度的测量进行输出记录。
14.如权利要求9所述的数据处理系统,其中注入的流体包括水。
15.如权利要求9所述的数据处理系统,其中多孔媒介片段包括岩心样本。
16.如权利要求9所述的数据处理系统,其中多孔媒介片段包括地下岩层片段。
17.一种数据存储设备,其具有存储在计算机可读媒介中的计算机可操作指令,该指令用于使数据处理系统对地球岩层的多孔媒介片段对流体注入量的饱和度进行测量,存储在数据存储设备中的指令使数据处理系统执行如下步骤:
将多孔媒介片段的系统样本长度分割成多个样本长度增额;
测量选定时间增额中向样本长度增额中注入的流体的注入量;
测量选定时间增额中在样本长度增额中由注入流体产生的分流曲线;
测量选定时间增额中注入流体在所述样本长度增额中的流体饱和度;
将选定时间增额中注入流体在所述样本长度增额中的流体饱和度的测量记录在数据处理系统中;
测量选定时间增额中未在样本长度增额中达到饱和的剩余流体容量。
18.如权利要求17所述的数据存储设备,还包括,所存储的指令包含用于使数据处理系统执行如下步骤的指令:
确定所测量的剩余流体容量是否表示存在用于注入多孔媒介片段的相邻长度样本增额中的剩余流体容量;
如果是,重复以下步骤:测量注入流体的注入量,测量流体的分流曲线,测量流体饱和度,对测量的流体饱和度进行记录,测量选定时间增额中未在相邻样本长度增额中达到饱和的剩余流体容量;
如果不是,测量选定时间增额中通过注入流体带来的长度样本增额的饱和度曲线。
19.如权利要求18所述的数据存储设备,还包括,所存储的指令包含用于使数据处理系统执行如下步骤的指令:
在测量选定时间增额中通过注入流体带来的长度样本增额的饱和度曲线后,增加所选定的时间增额从而形成新的选定时间增额。
20.如权利要求19所述的数据存储设备,还包括,所存储的指令包含用于使数据处理系统执行如下步骤的指令:
测量新的时间增额中向长度样本增额中注入的流体的注入量;
测量新的时间增额中在长度样本增额中由注入流体产生的分流曲线;
测量新的时间增额中注入流体在所述长度样本增额中的流体饱和度;
对新的时间增额中注入流体在所述长度样本增额中的流体饱和度的测量进行记录;
测量新的时间增额中未在长度样本增额中达到饱和的剩余流体容量。
21.如权利要求17所述的数据存储设备,其中注入的流体包括水。
22.如权利要求17所述的数据存储设备,其中多孔媒介片段包括岩心样本。
23.如权利要求17所述的数据存储设备,其中多孔媒介片段包括地下岩层片段。
24.如权利要求17所述的数据存储设备,还包括,所存储的指令包含用于使数据处理系统的输出显示器执行如下步骤的指令:
对确定的注入流体的流体饱和度的测量进行输出显示。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/974,434 US20120158309A1 (en) | 2010-12-21 | 2010-12-21 | Modeling Immiscible Two Phase Flow in a Subterranean Formation |
US12/974,434 | 2010-12-21 | ||
PCT/US2011/062015 WO2012087488A2 (en) | 2010-12-21 | 2011-11-23 | Modeling immiscible two phase flow in a subterranean formation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103329225A true CN103329225A (zh) | 2013-09-25 |
CN103329225B CN103329225B (zh) | 2016-11-16 |
Family
ID=45346557
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180061928.5A Active CN103329225B (zh) | 2010-12-21 | 2011-11-23 | 地层中的不相融的两相流体的塑造 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20120158309A1 (zh) |
EP (1) | EP2656360A2 (zh) |
CN (1) | CN103329225B (zh) |
CA (1) | CA2821004C (zh) |
WO (1) | WO2012087488A2 (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9823379B2 (en) | 2014-02-13 | 2017-11-21 | Groundmetrics, Inc. | System and method for mapping deep anomalous zones of electrical resistivity |
CA2940100C (en) | 2014-02-21 | 2020-05-26 | Groundmetrics, Inc. | Method for mapping the propagation of earth fractures |
CN113969768A (zh) * | 2020-07-23 | 2022-01-25 | 中国石油化工股份有限公司 | 一注多采井组定向赋能-差异释放式体积水驱方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1811413A (zh) * | 2004-12-20 | 2006-08-02 | 普拉德研究及开发股份有限公司 | 评估地质层中流体饱和度特性的方法 |
US20080046223A1 (en) * | 2006-08-16 | 2008-02-21 | Benoit Noetinger | Method of optimizing enhanced recovery of a fluid in place in a porous medium by front tracking |
US20100057413A1 (en) * | 2008-09-02 | 2010-03-04 | Chevron U.S.A. Inc. | Indirect-error-based, dynamic upscaling of multi-phase flow in porous media |
US20100312535A1 (en) * | 2009-06-08 | 2010-12-09 | Chevron U.S.A. Inc. | Upscaling of flow and transport parameters for simulation of fluid flow in subsurface reservoirs |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA1054911A (en) * | 1976-01-07 | 1979-05-22 | Exxon Production Research Company | Method for determining gas saturation in reservoirs |
US6052520A (en) * | 1998-02-10 | 2000-04-18 | Exxon Production Research Company | Process for predicting behavior of a subterranean formation |
-
2010
- 2010-12-21 US US12/974,434 patent/US20120158309A1/en not_active Abandoned
-
2011
- 2011-11-23 CN CN201180061928.5A patent/CN103329225B/zh active Active
- 2011-11-23 EP EP11794896.8A patent/EP2656360A2/en not_active Withdrawn
- 2011-11-23 CA CA2821004A patent/CA2821004C/en active Active
- 2011-11-23 WO PCT/US2011/062015 patent/WO2012087488A2/en active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1811413A (zh) * | 2004-12-20 | 2006-08-02 | 普拉德研究及开发股份有限公司 | 评估地质层中流体饱和度特性的方法 |
US20080046223A1 (en) * | 2006-08-16 | 2008-02-21 | Benoit Noetinger | Method of optimizing enhanced recovery of a fluid in place in a porous medium by front tracking |
US20100057413A1 (en) * | 2008-09-02 | 2010-03-04 | Chevron U.S.A. Inc. | Indirect-error-based, dynamic upscaling of multi-phase flow in porous media |
US20100312535A1 (en) * | 2009-06-08 | 2010-12-09 | Chevron U.S.A. Inc. | Upscaling of flow and transport parameters for simulation of fluid flow in subsurface reservoirs |
Also Published As
Publication number | Publication date |
---|---|
WO2012087488A3 (en) | 2012-10-18 |
US20120158309A1 (en) | 2012-06-21 |
EP2656360A2 (en) | 2013-10-30 |
CA2821004C (en) | 2020-01-07 |
WO2012087488A2 (en) | 2012-06-28 |
CA2821004A1 (en) | 2012-06-28 |
CN103329225B (zh) | 2016-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9187984B2 (en) | Methods and systems for machine-learning based simulation of flow | |
US10198535B2 (en) | Methods and systems for machine-learning based simulation of flow | |
US10087721B2 (en) | Methods and systems for machine—learning based simulation of flow | |
US10151855B2 (en) | Sequential fully implicit well modeling of transmissibility for reservoir simulation | |
EP3559401B1 (en) | Method and system for stable and efficient reservoir simulation using stability proxies | |
US20130096899A1 (en) | Methods And Systems For Machine - Learning Based Simulation of Flow | |
Jansen et al. | Closed-loop reservoir management | |
US11066911B2 (en) | Operating hydrocarbon wells using modeling of immiscible two phase flow in a subterranean formation | |
CN103329225A (zh) | 地层中的不相融的两相流体的塑造 | |
CA3013807A1 (en) | Sequential fully implicit well model with tridiagonal matrix structure for reservoir simulation | |
Natvig et al. | Multiscale Reservoir Simulation of High-Resolution Models | |
Møyner | Next generation multiscale methods for reservoir simulation | |
Trehan et al. | Machine-learning-based estimation of upscaling error, with application to uncertainty quantification |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |