CN1358289A - 按比例增大非结构化栅格(grid)的渗透率的方法 - Google Patents
按比例增大非结构化栅格(grid)的渗透率的方法 Download PDFInfo
- Publication number
- CN1358289A CN1358289A CN00809383.0A CN00809383A CN1358289A CN 1358289 A CN1358289 A CN 1358289A CN 00809383 A CN00809383 A CN 00809383A CN 1358289 A CN1358289 A CN 1358289A
- Authority
- CN
- China
- Prior art keywords
- grid
- node
- permeability
- line
- unit
- 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.)
- Pending
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 113
- 238000000034 method Methods 0.000 title claims abstract description 74
- 230000004907 flux Effects 0.000 claims abstract description 45
- 230000011218 segmentation Effects 0.000 claims description 46
- 238000013341 scale-up Methods 0.000 claims description 25
- 239000012530 fluid Substances 0.000 claims description 21
- 238000012546 transfer Methods 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 description 25
- 238000005516 engineering process Methods 0.000 description 9
- 230000010354 integration Effects 0.000 description 8
- 229910003460 diamond Inorganic materials 0.000 description 7
- 239000010432 diamond Substances 0.000 description 7
- 230000008569 process Effects 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 230000002596 correlated effect Effects 0.000 description 3
- 238000001764 infiltration Methods 0.000 description 3
- 230000008595 infiltration Effects 0.000 description 3
- 238000004215 lattice model Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000036962 time dependent Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 244000287680 Garcinia dulcis Species 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK 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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
-
- 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
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Medicines Containing Material From Animals Or Micro-Organisms (AREA)
- Image Generation (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
提供一种把与代表多孔介质的细分度栅格单元(图1)相关的渗透率按比例增大为与代表所述多孔介质的非结构化粗分度栅格单元(图3)相关的渗透率的方法。通过把粗分度栅格用作计算栅格的起源,产生面内非结构化的Voronoi计算栅格。随后利用与细分度栅格相关的渗透率填充计算栅格。导出计算栅格的流动方程,求解所述流动方程,随后计算所述计算栅格的节点间通量和压强梯度。这些节点间通量和压强梯度被用于计算与粗分度栅格相关的节点间平均通量和平均压强梯度。与粗栅格相关的节点间平均通量和平均压强梯度随后被用于计算与粗分度栅格相关的按比例增大的渗透率。
Description
技术领域
本发明涉及模拟多孔介质中的流体流动,更具体地说,涉及基于流动的,把与代表多孔介质的细格(fine-grid)***相关的渗透率按比例增大为与同样代表多孔介质的粗格(coarse-grid)***相关的渗透率的方法。
背景技术
在工业领域中,通过利用计算机,数值模拟被广泛用作物质***的模拟方法。在大多数情况下,需要模拟物质***中发生的迁移过程。迁移的通常是质量、能量、动量或它们的某些组合。通过利用数值模拟,能够在不进行实际的实验室实验或现场测试的情况下,再现和观察物理现象,并确定设计参数。于是,可显著减少设计时间和费用。
最感兴趣的一种模拟是根据含油气层的数值模型的性能,推断真实的含油气层的行为的方法。储油层模拟的目的是充分了解在储油层中发生的复杂化学、物理和流体流动过程,预测储油层的未来行为,从而使油气开采率达到最大。储油层模拟通常指的是储油层内流动的流体动力学,但是从广义来说,储油层模拟也可涉及整个石油***,包括储油层、注入井、采油井、地面出油管线及地面加工设施。
数值模拟的原理是利用计算机,数值求解描述物理现象的方程式。这种方程式通常是普通的微分方程和偏微分方程。这些方程式通常是利用诸如有限差分法、有限元法、有限体积法之类数值方法来求解的。在这些方法中的每一种中,要模拟的物质***被分成更小的单元(一组单元被称为栅格或网格),在每个单元中持续变化的状态变量由每个单元的多组数值表示。原始差分方程由一组代数方程代替,以便表述每个更小单元内的质量、能量和/或动量守恒的基本原理,以及表述单元之间的质量、能量和/或动量迁移的基本原理。这些代数方程的数目多达几百万。每个单元的有限数目的数值对连续变化数值的这种代替被称为“离散化”。为了分析随时间变化的现象,必须计算在称为时间步长(timestep)的离散时间间隔时的物理量,而不考虑随时间变化的连续变化条件。迁移过程的瞬态模拟按照时间步长的顺序进行下去。
在储油层的典型模拟中,在感兴趣范围中的特定点寻求主要未知量(通常是压强和相饱和度或相组成)的解答。这样的点被称为“栅格节点(gridnode)”或者更一般地被称为“节点”。围绕这样的节点构造单元,并且栅格被定义为这样的一组单元。诸如孔隙率和渗透率之类性质被假定为在单元内是恒定不变的。诸如压强和相饱和度之类的其它变量在节点处被确定。两个节点之间的链接被称为“连线”。两个节点之间的流体流动通常被模拟为沿着它们之间的连线的流动。
在常规的储油层模拟中,大多数栅格***是结构化的。即,单元具有相似的形状,以及相同数目的侧面或表面。最常用的结构化栅格是笛卡尔或放射状栅格,其中在二维空间中,每个单元具有四个侧面,在三维空间中,每个单元具有六个侧面。虽然结构化栅格易于使用,但是它们缺少适应储油层及油井几何形状方面的变化的灵活性,并且通常不能有效处理储油层中岩石和流体的物理性质的空间变化。已提出了供结构化栅格效率不高的这种情形中使用的灵活栅格。当栅格由具有可因位置而异的形状,大小及侧面或表面数目的多边形组成时,该栅格被称为是灵活的或非结构化栅格。和结构化栅格相比,非结构化栅格可更容易地与复杂的储油层特征相符,由于这个原因,已提出在储油层模拟中采用非结构化栅格。
储油层模拟中可使用的一种灵活栅格是Voronoi栅格。Voronoi单元被定义为与到其它任何节点的距离相比,更接近于其节点的空间区域,Voronoi栅格由这样的单元组成。每个单元与一个节点,以及一系列相邻单元相关。从几何学来说,Voronoi栅格是局部直交的;即,单元边界垂直于位于每条边界两侧的节点的连线。为此,Voronoi栅格也被称为垂直平分(PEBI)栅格。矩形栅格块(笛卡尔栅格)是Voronoi栅格的一种特殊情况。由于可自由选择节点的位置,因此PEBI栅格具有表现变化极大的储油层几何形状的灵活性。通过在指定范围内分配节点,并且以使各个单元含有与到其它任何节点位置的距离相比,更接近于其节点位置的所有点的方式,产生单元边界,产生PEBI栅格。由于PEBI栅格中节点间的连线被单元边界垂直平分,因此这显著简化了流动方程式的求解。关于PEBI栅格产生的更详细说明,请参见Palagi,C.L和Aziz,K.在66thAnnual Technical ConferencE and Exhibition,TX,Dallas(1991年10月6-9日)提交的论文SPE 22889:“Use of Voronoi Grid in ReservoirSimulation”。
连接PEBI单元的相邻节点形成的网格通常被称为Delaunay网格(如果仅由三角形构成)。在二维Delaunay网格中,储油层被分成三角形,栅格节点位于三角形的顶点,从而三角形填满储油层。当经过三角形顶点的圆(外心)在其内部并不含有其它任何节点时,则这种三角测量是Delaunay。在三维空间中,储油层区被分解成四面体,从而储油层体积被完全充满。当经过四面体顶点的球体(外接球)不含有其它任何节点时,这样的三角测量是Delaunay网格。Delaunay三角测量技术众所周知;例如参见Migdal等的美国专利5886702。
借助先进的储油层表征技术,利用几百万栅格单元模拟储油层的地质结构和志层学是常见的,每个栅格单元填充有一种储油层性质,所述储油层性质包括(但不局限于)岩石类型、孔隙率、渗透率、初始空隙流体饱和度及相对渗透率和毛细压强函数。但是,通常是利用远远少得多的栅格单元来进行储油层模拟的。直接把细格模型应用于储油层模拟通常是不可行的,因为细格模型的细微程度对计算资源提出了过高的要求。于是,需要一种把细格地质储油层模型转换或按比例增大为粗格模拟模型,同时尽可能地保持细格模型的流体流动特性。
储油层模拟的一个关键流体流动性质是渗透率。渗透率是岩石通过岩石中的互连小孔传输流体的能力。在含油气层中,渗透率可发生显著变化。一般,利用来自钻孔岩心的数据,产生用于细分度模型(地质模型)的渗透率。对于模拟单元来说,通过确定有效渗透率,解决地质模型的不均匀性。不均匀介质的有效渗透率一般被定义为等同的均匀介质的渗透率,对于相同的边界条件,所述等同的均匀介质将给出相同的通量(单位时间通过指定面积的流体流量)。确定有效渗透率(通常称为渗透率按比例增大)并不简单。主要的困难在于储油层中渗透率不均匀性和采用的边界条件的相互影响。
已提出了多种不同的按比例增大技术。这些技术的大多数可被表征为(1)直接法或(2)基于流动的方法。直接法的例子是各种类型的简单求平均(例如,算术求平均,几何求平均和调和求平均)及逐次重正化。基于流动的技术涉及求解流动方程,并且考虑渗透率的空间分布。一般地,基于流动的方法需要更大的计算量,但是精度高于直接法。
在下述论文:Wen,X.H.和Gomez-Hernandez,J.J.的“UpscalingHydraulic Conductivities in Heterogeneous Media:An Overview”,Journal of Hydrology,Vol.183(1996)9-32;Begg,S.H.;Carter,R.R.和Dranfield,P.的“Assigning Effective Values to Simulator GridblockParameters for Heterogeneous Reservoirs,”SPE Reservoir Engineering(1989年11月)455-465;Durlofsky,L.J.,Behrens,R.A.和Bernath,A.在Annual Technical Conference and Exhibition,Dallas,TX(1995年10月22-25日)提交的论文SPE 30709“Scale Up of Three DimensinalReservoir Descriptions”;和Li,D.,Cullick,A.,Lake,L.W.的“GlobalScale-up of Reservoir Model Permeability with Local GridRefinement”,Journal of Petroleum Science and Engineering,Vol.14(1995)1-13中提供了不同的按比例增大技术的综述。过去提出的按比例增大技术主要集中在结构化栅格上。需要一种把与细刻度地质模型相关的渗透率按比例增大为与非结构化的粗刻度储油层模拟模型相关的渗透率的方法。
发明内容
提供一种把与代表多孔介质的细分度栅格单元相关的渗透率按比例增大为与代表所述多孔介质的非结构化粗分度栅格单元相关的渗透率的方法。第一步是通过把粗分度栅格用作计算栅格的起源,产生面内非结构化的Voronoi计算栅格。计算栅格的单元小于粗分度栅格的单元,并且计算栅格及粗分度栅格的每个单元具有一个节点。随后利用与细分度栅格相关的渗透率填充计算栅格。导出计算栅格的流动方程,最好是单相、定态压强方程,求解所述流动方程,随后计算所述计算栅格的节点间通量和压强梯度。这些节点间通量和压强梯度被用于计算与粗分度栅格相关的节点间平均通量和平均压强梯度。与粗栅格相关的节点间平均通量和平均压强梯度随后被用于计算与粗分度栅格相关的按比例增大的渗透率。
在一个优选实施例中,由粗分度栅格构造计算栅格,产生平行于粗分度栅格的节点间连线的计算栅格的节点间连线。计算栅格的单元最好与细分度单元近似大小相同。最好通过把细分度栅格某一单元的预定渗透率分配给计算栅格的指定节点,利用渗透率填充计算栅格,如果使计算栅格叠加在细分度栅格上,则所述某一单元会包含指定节点的位置。为计算栅格导出的流动方程最好是单相、定态方程。最好通过只利用落入计算栅格预定子域内的计算栅格节点间连线,计算与粗分度栅格相关的节点间平均通量和平均压强梯度,并且最好只利用计算栅格的,平行于和所述子域相关的粗分度栅格的节点间连线的节点间连线,进行这样的计算。最好通过计算关于粗分度栅格计算的节点间平均通量与节点间平均压强梯度的比值,确定与粗分度栅格的节点间连线相关的渗透率。
附图说明
参考下述详细说明及附图,将更好地理解本发明及其优点,在附图中,相同的部分被赋予相同的数字和字母。
图1图解说明了通常用在多孔介质的地质模拟中的二维、结构化、细分度栅格。
图2图解说明了具有36个单元的二维、非结构化、粗分度栅格,所述36个单元中的四个单元被表示为具有标示为I,J,K和L的节点。
图3图解说明了由图2的粗分度栅格产生的二维、非结构化计算栅格;利用为2的栅格细化比产生计算栅格的单元。
图4图解说明了在本发明的实践中,由图2的具有节点I,J和K的三个非结构化单元产生的10个非结构化单元;通过利用为3的栅格细化比,由这三个单元产生所述10个单元。
图5图解说明了在本发明的实践中,由图2的具有节点I,J和K的三个非结构化单元产生的15个非结构化单元;通过利用为4的栅格细化比,由这三个单元产生所述15个单元。
图6图解说明了利用为4的栅格细化比,由图2的栅格产生的二维计算栅格,并显示了一个钻石形子域KILJ。
图7图解说明了图2的具有节点I,J,K和L的四个粗栅格单元,以及图6的钻石形子域R。
图8图解说明了为图6中的栅格的子域,并且只包括图6的积分域R的计算栅格。
图9图解说明了为图6中的栅格的子域的计算栅格,所述计算栅格包括图6的积分域R和围绕域R的缓冲区(复盖区)。
附图图解说明了实践本发明的方法的具体实施例。附图并不意味着把所述具体实施例的一般及预期修改得到的其它实施例排除在本发明的范围之外,
具体实施方式
本发明提供一种基于流动的,把与代表多孔介质的细刻度,地质栅格***相关的渗透率按比例增大为与代表相同多孔介质的粗刻度,非结构化栅格***相关的估计渗透率的方法。本发明特别适用于把通常由结构化单元代表的储油层地质模型相关的渗透率按比例增大为与粗刻度,非结构化Voronoi栅格(也称为垂直平分或PEBI栅格)相关的的渗透率。
作为执行按比例增大方法的第一步,利用粗栅格作为其起源,产生面内(areally)非结构化Voronoi栅格(这里通常称为“计算栅格”)。按照本发明的实践,粗栅格和计算栅格两者的节点间连线都相互排列成行(平行)。随后用来自细刻度的地质栅格的渗透率填充计算栅格。下一步骤是为非结构化的计算栅格建立流动方程,最好是单相、定态压强方程,求解该流动方程,并计算计算栅格的节点间通量和压强梯度。随后把这些通量和压强梯度用于计算与粗栅格相关的节点间平均通量和平均压强梯度。随后利用先前关于粗栅格计算的平均通量和平均压强梯度,计算与粗栅格连线相关的渗透。
下面将参考附图说明本发明的一个实施例。虽然附图只图解说明了二维(2-D)栅格***,但是本发明并不局限于2-D栅格。如同稍后所述那样,本发明也可应用于三维(3-D)栅格。
图1图解说明了供表现诸如蓄水层或含油气层之类的多孔介质之用的细刻度地质栅格的一个例子。虽然图1表示了64个单元的8×8笛卡尔栅格,但是应明白典型的地质栅格可含有几百万个单元。图1中的所有单元都是结构化单元,这是是一种典型的地质栅格,每个单元具有位于其中心的节点。图1中,以及在其它附图中,粗圆点表示单元节点;连续直线表示单元的边界;在图1-5和7中,节点间的虚线被称为节点间连线。由这些图中的虚线形成的三角形构成Delaunay网格。
图2图解说明了适用于储油层模拟的含有粗刻度单元的非结构化,PEBI栅格***的一个例子。图1和图2的栅格代表多孔介质的相同范围。在大多数储油层模拟应用中,图2的粗刻度单元的平均尺寸比图1的细刻度单元的平均尺寸大几倍。储油层模拟模型可涉及,例如每个储油层模拟单元约10到1000个地质单元,同时长度刻度从地质模型中的10米或更大变化为储油层模拟模型中的100米或更大。虽然图1中描绘的单元小于图2中描绘的单元(出于表示清楚的原因),在实际情况中,它们甚至可更小。用于产生细刻度地质栅格及粗刻度模拟栅格的技术众所周知。本发明发现一种基于流动的,把与细刻度栅格(例如图1的结构化栅格)相关的渗透率转换为与非结构化粗栅格***(例如图2的非结构化栅格)相关的渗透率的高效按比例增大方法。
图3图解说明了在本发明实践中产生的计算栅格的一个例子,本发明可用于把图1的地质模型的渗透率按比例增大为图2的非结构化储油层模拟模型的渗透率。本发明的方法中的第一步是为图1和2覆盖的相同范围产生计算栅格。
以这样的方式产生计算栅格,以便产生比粗栅格单元小的非结构化PEBI单元。计算栅格的各个节点间连线都平行于粗栅格的节点间连线。如同下面将参考优选实施例详细说明的那样,这种平行对准简化了压强梯度和通量的按比例增大,而压强梯度和通量的按比例增大是本发明方法中的一个重要步骤。
通过分配除粗栅格的节点之外的节点,并重新生成PEBI单元,产生计算栅格。所述另外的节点按比例分配在所有粗栅格连线上,并且位于各个Delauany三角形内。相同数目的新节点被增加到每个连线上,并且可增加任意数目的新节点。增大节点数目,导致增大可由粗栅格单元形成的较小PEBI单元的数目。这里把计算栅格的较小单元的数目的增大称为细化度的提高。通过计数每个粗栅格连线的额外节点数,可测量细化度。如果每个粗栅格连线增加一个节点,则细化比为2,因为它把粗栅格连线分成两个相等的小栅格连线;如果每个粗栅格连线增加2个节点,则细化比为3,如果添加3个节点,则细化比为4,依次类推。
在细化比为2的条件下,产生图3的PEBI栅格。参见图3,节点I,J和K表示图2中图解说明的具有节点I,J和K的三个单元的节点。为了产生图3的栅格,在节点K和I之间半路引入了另一个节点,在节点K和J之间半路引入了另一个节点,在节点I和J之间半路引入了另一个节点,对图2的所有其它节点间连线依次类推。随后基于图2的原始节点和新增加的节点,产生图3的PEBI栅格。
图4和5中分别表示了细化比为3和4的计算栅格的例子。图4表示由图2的具有节点I,J和K的三个单元产生的10个PEBI单元,表示为3的细化比。为了产生图4的PEBI单元,按比例在节点K和I之间增加两个节点10和11;按比例在节点K和J之间增加两个节点13和14;按比例在节点I和J之间增加两个节点15和16;并在由节点K,I和J形成的三角形内引入一个节点12。由图4的10个单元的节点间连线形成的Deaaunay网格(虚线)由具有相同角度的9个同样大小的的Delaunay三角形组成。图4的所有9个较小的Delaunay三角形类似于图2的Delaunay三角形IJK。
图5表示由图2的具有节点I,J和K的三个单元产生的15个PEBI单元,表示为4的细化比。为了产生图5中所示的15个PEBI单元,按比例在图2的节点K和I之间增加三个节点20、21和22;按比例在节点K和J之间增加三个节点23、24和25;按比例在节点I和J之间增加三个节点26、27和28;并在由节点K,I和J形成的三角形内相称地增加三个节点29、30和31。节点29、30和31位于通过节点20-28,平行于三角形KIJ的侧边画出的直线的交点。从而,图5中PEBI单元的所有节点间连线都平行于图2的节点K和I间的节点间连线,节点K和J间的节点间连线,以及节点I和J间的节点间连线之一。例如,参见图5,节点20和23间的节点间连线平行于图2的I和J间的节点间连线;节点23和24间的节点间连线平行于图2的节点K和J间的节点间连线;节点29和30间的节点间连线平行于图2的节点K和I间的节点间连线。
对于指定的细化比,可产生计算栅格的更小的相似Delaunay三角形的数目为n2与粗刻度栅格的Delaunay三角形数目的乘积,这时n是所需的整数细化比。例如,细化比为4的计算栅格将具有16个(42)类似于图2的Delaunay三角形KIJ,但是小于图2的Delaunay三角形KIJ的Delaunay三角形(图5的Delaunay栅格)。
细化比可以是大于1的任意整数。对于大多数储油层模拟应用,作为非限制性例子,细化比为2-10,更典型的是为2-5。所需的细化比将取决于地质栅格和粗栅格的相对大小,以及所需的计算栅格单元数目。细化比最好被选择为对于相同的储油层范围,产生大小基本和地质栅格的单元相同的计算栅格单元,最好细化比被选择为产生稍小于地质单元的计算栅格单元。
当每个Delaunay网格三角形的所有角度形成小于90°(锐角)的粗栅格连线时,或者当一个角度最多等于90°(直角)时,本发明的方法产生PEBI单元。但是,如果Delaunay三角形的任意角度大于90°(钝角)时,不存在真正的垂直平分。对于这种三角形,本发明的栅格产生技术将产生不平行于上层粗栅格连线的连线。如果产生了任意这种非平行连线,则在如下更详细说明的从计算栅格到粗栅格按比例增大渗透率中最好不使用这种非平行连线。最好这样产生粗栅格,使得在Delaunay网格中不存在任何钝角三角形。在大多数情况下,这可通过调节节点位置来实现。
本发明的小栅格生成技术并不局限于三角形连线网格(Delaunay网格);它还适用于连线网格包括矩形的结构化栅格。本计算栅格生成技术也适用于包括三角形和矩形的组合的连线网格。
渗透率填充
一旦产生了计算栅格,本发明的下一步骤是在计算栅格上填充渗透率。通过利用与细刻度地质模型(例如图1)相关的预定渗透率,实现该填充步骤。可使用各种填充方法。一种方法如下所述把来自地质栅格的渗透率分配给计算栅格的各个节点。为计算栅格上的指定节点,搜索该地质栅格,直到找到含有计算栅格的节点的地质栅格的单元为止,假定计算栅格叠加在地质栅格上。含有计算栅格的指定节点的地质单元的预定渗透率被指定为该指定节点的渗透率。随后通过求取形成节点间连线的两个节点的渗透率的平均值(最好采用调和求平均),计算计算栅格的每个节点间连线渗透率。利用渗透率填充计算栅格的另一种方法是把得自于地质单元的渗透率指定给计算栅格的节点间连线的中点。为计算栅格的节点间连线的指定中点,查找含有指定中点的地质栅格的单元(假定计算栅格叠加在地质栅格上),并把该地质栅格的渗透率分配给指定中点。如果地质模型具有对角线渗透张量,则在填充前,可根据连接方向组合这两个面渗透率。
流体流动方程
在用渗透率填充计算栅格之后,本发明实践中的下一步是导出计算栅格的流动方程,最好是计算栅格各个单元的单相、定态压强方程,并利用一组假定的边界条件,为各个计算栅格单元求解该方程。单相、定态压强方程为:
下述说明利用了数学符号,在说明书中内,多数数学符号当其出现时对其进行了定义。另外,为了完整性起见,在下述说明之后给出了含有这里使用的各种符号的定义的符号表。
假定(但不限定于)标量连线渗透率,对于PEBI栅格来说,方程(1)的离散形式为:
这里,T是传输率,下标j指的是所关心的节点,下标l指的是其所有邻近节点。本说明书中使用的术语“传输率”指的是指定粘度的流体在压差下,通过单元边界(或者节点间连线)的能力的量度。更具体地说,传输率被本领域技术人员认为是流体在多孔介质内两个相邻单元间的流动能力的量度。传输率被表述为
这里k是多孔介质的有效渗透率,A是相邻单元之间的边界面积,Δs是流体在这两个单元间移动时必须行进的平均或特征距离。
一组可接受的边界条件沿流动方向施加恒定的压强梯度,并假定在反方向上没有任何流动(无反向流动边界条件)。如图6中所示,流体流动被假定为只沿箭头方向产生(从左到右)。适于实践本发明的另一组边界条件是线性边界条件。在线性边界条件中,顺着所有边界,沿流动方向施加恒定的压强梯度。这可导致流过有助于按比例增大隔离区域中的渗透率的所有边界。上述边界条件不应被理解为对本发明的限制;也可使用其它适当的边界条件。
对于2-D应用来说,在两个主要方向的每个方向上求解压强方程(对于3-D应用来说,沿三个主要方向求解压强方程)。随后计算计算栅格的所有节点间连线的通量和压强梯度。储油层模拟领域中的技术人员通晓计算栅格各单元的恰当单相压强方程的导出,压强方程的求解,以及计算栅格的所有节点间连线的通量和压强梯度的计算。例如,参见Verma,S.在Reservoir Simulation Symposium,Dallas TX(1997年6月8-11日)提交的论文SPE 37999:“A Control Volume Scheme forFlexible Girds in Reservoir Simulation”。
按比例增大
一旦确定了计算栅格的通量和压强梯度,则计算与粗栅格的连线相关的平均通量和平均压强梯度。在与每个粗栅格单元相关的预定积分子域内求这些压强梯度和通量的平均值,最好在与粗栅格的各个节点间连线相关的子域内求这些压强梯度和能量的平均值。按比例增大后的通量与按比例增大后的压强梯度的比值给出按比例增大的渗透率。所述按比例增大的渗透率随后被用于计算传输率。
下面将提供从计算栅格到粗栅格按比例增大渗透率的数学基础。该说明利用了几个数学符号,其中一些符号当其在说明中出现时,进行了定义。另外,为了完整性起见,在详细说明之后给出了含有这里使用的各种符号的定义的符号表。
关于沿方向s的1-D单相流动的Darcy定律由下式给出:
相对于体积,对方程式(3)积分,计算子域R的体积内的平均通量,得到:
假定粘度恒定不变,则在粗刻度栅格,沿方向s,子域R的按比例增大的渗透率(
k)现在被定义为:
借助Rubin,Y.和Gomez-Hernandez,J.J.的论文“A StochasticApproach to the Problem of Upscaling ofConductivity in DisorderedMedia:Theory and Unconditional Numerical Simulations”(WaterResources Research,Vol.26,No.4,第691-701页,1990年4月)中描述的下述关系式,方程式(5)可表述为按比例增大的通量与按比例增大的压强梯度的比值:
在细分度下,近似方程式(6)中的积分,并替代连接体积Vf(Vf=AfΔsf),得到关于按比例增大的渗透率(kc)的下述表达式: 现在,连线传输率可被计算为:
通过利用根据压强方程解答计算的压强梯度和通量,可把方程式(7)应用于PEBI计算栅格。重要的一步是确定与将越过其按比例增大计算栅格渗透率的各个粗栅格连线相关的积分子域或区域。
参见图6和7,为了在PEBI栅格上按比例增大面连线渗透率,适当的子域是由图2的两个三角形(KIJ和IJL)组成的钻石形子域R。也可使用其它适宜的子域形状。现在,通过利用预定的选择标准,使用子域内的所有或部分的计算栅格节点间连线,方程式(7)可用于计算每个子域的按比例增大的渗透率。在一种简单而有效的方法中,按比例增大的渗透率只使用子域R内与代表子域R的粗栅格节点间连线平行的计算栅格节点间连线。在另一实施例中,按比例增大中可使用子域R中的所有计算栅格节点间连线。
参见图7,子域R中的计算栅格节点间连线被用于获得节点I和J间的节点间连线的按比例增大的渗透率。由于在面范围(areal domain)内,沿两个主要方向求解压强方程式,因此得到两个按比例增大的渗透率。如果使用无反向流动边界条件,则根据对于各个粗栅格面连线,具有最高平均压强梯度的面方向(areal direction),选择最终的按比例增大的渗透率。通过利用压强方程的这两个解答,其它技术也可用于按比例增大渗透率。例如,当使用线性边界条件时,可如下得到这两个压强梯度的线性组合:
常数a和b被确定成使最终得到的压强梯度与将关于其按比例增大渗透率的连线方向对齐。现在,按比例增大的渗透率可表述为:
本发明的按比例增大技术也可应用于3-D结构化及3-D分层PEBI栅格(也被称为2-D PEBI栅格)。也可扩展到3-D非结构化栅格。分层PEBI栅格在面内是非结构化的,在垂直方向上是结构化的(分层的)。产生这种栅格的一种方法是把2-D面PEBI栅格投射到地质顺序表面上。可逐层或多层按比例增大面连线渗透率。为了按比例增大垂直(层间)连线渗透率,可利用垂直方向上的压强梯度求解流动方程。(1)Heinemann,Z.E.等的“Modeling Reservoir Geometry With IrregularGrids”(SPE Reservoir Engineering,1991年5月),和(2)Verma,S.等的“A Control Volume Scheme for Flexible Grids in ReservoirSimulation”(SPE 37999,SPE Reservoir Simulation Symposium,Dalls,TX,1997年6月)描述了分层3-D栅格的构造。本发明的按比例增大方法也可扩展到完全的3-D PEBI栅格。
本发明非常有用的一方面是按比例增大方法可用于同时按比例增大任意数目的粗栅格连线的渗透率。可在对应于大的粗栅格域(按比例增大域)的计算栅格上,求解流体流动方程,从而允许同时按比例增大大量的渗透率。例如,可相对于图6的整个计算栅格求解流体流动方程,从而按比例增大图2中所示的所有粗栅格节点间连线的渗透率。另一方面,可在对应于一个或几个粗栅格连线的较小的计算栅格上,求解流体流动方程。图8中表示了用于按比例增大一个粗栅格连线渗透率的计算栅格的例子。这里,该计算栅格和按比例增大粗栅格连线IJ的渗透率中使用的钻石形积分域R相符。
也可在比按比例增大区域大的计算栅格上,求解流动方程,从而提供围绕按比例增大区域的缓冲区或复盖层,而不管所述按比例增大区域是由单个或多个连线组成。图9中表示了这样的一个例子,这里计算栅格包括连线IJ的钻石形积分域R(KILJ)和围绕域R的缓冲区或复盖区。这些选项提供了选择按比例增大域大小的灵活性,以及选择用于按比例增大的计算栅格大小的灵活性,在实际的场边界条件不同于按比例增大边界条件的某些应用中,对于按比例增大精度来说,这是重要的。
虽然已借助单相、无重力、定态、标量渗透率方程式给出了本发明的按比例增大方法的数学说明,但是该方法决不应局限于这些假定。利用本专利中公开的信息,本领域的技术人员可把本发明的按比例增大方法扩展到包括多相流动、非稳态流动、全渗透张量及重力。
虽然图中没有表示,但是,在实践本发明的过程中,利用适当的计算机化可视化设备,可在2-D和3-D非结构化栅格上随意显现被模拟区域的性质。性质的可视化有助于分析本发明的应用。在2-D中,可通过根据预定的色度,以颜色显示与各个连线相关的钻石形区域,呈现连线性质(例如通量、压强梯度或连线渗透率)。该钻石形区域的顶点包括形成所述连线的两个节点和所述连线的垂直平分线的两个端点。在3-D中,可用球或棒网表示Delaunay三角测量。可通过利用色度,对球着色,显示节点性质(例如压强、饱和度、组分或孔隙率)。类似地,通过对棒着色,可显示连线性质。色度可显示性质的量级。为了显示矢量性质(例如,压强梯度或通量),可使用箭头代表流动方向。
上面说明了本发明的原理及设想的应用该原理的最佳方式。对本领域的技术人员来说,在不脱离如下述权利要求中限定的本发明的精神和范围的情况下,显然可对上面说明的实施例做出各种变化。于是,本发明并不局限于所展示和说明的具体细节。
符 号
a,b=常数,无量纲
A=面积,平方英寸
K=渗透率,md
P=压强,psi
T=传输率,md-ft
u=Darcy通量,ft/天
V=体积,立方英寸
ΔP=压差,psi
Δs=连线长度,ft
μ=粘度,cp
下标c=粗栅格f=计算栅格i=连线下标j,l=节点下标R=积分子域
Claims (28)
1.一种把与代表多孔介质的细分度栅格的单元相关的渗透率按比例增大为与代表所述多孔介质的非结构化粗分度栅格的单元相关的渗透率的方法,包括下述步骤:
(a)通过把粗分度栅格用作计算栅格的起源,产生面内非结构化的Voronoi计算栅格,计算栅格的单元小于粗分度栅格的单元,并且所述计算栅格及所述粗分度栅格的每个单元具有一个节点;
(b)利用与细分度栅格相关的渗透率填充计算栅格;
(c)导出计算栅格的流动方程,求解所述流动方程,并计算所述计算栅格的节点间通量和压强梯度;
(d)利用在步骤(c)中计算的通量和压强梯度,计算与粗分度栅格相关的节点间平均通量和平均压强梯度;
(e)利用在步骤(d)中计算的平均通量和平均压强梯度,计算与粗分度栅格相关的按比例增大的渗透率。
2.按照权利要求1所述的方法,其中计算栅格单元的节点间连线平行于粗分度栅格的节点间连线。
3.按照权利要求1所述的方法,其中计算栅格的单元与细分度栅格的单元的大小近似相同。
4.按照权利要求1所述的方法,其中计算栅格的单元小于细分度栅格的单元。
5.按照权利要求1所述的方法,其中在步骤(b)中填充的渗透率被分配给计算栅格的节点,分配给计算栅格的指定节点的渗透率对应于细分度栅格的某一单元的预定渗透率,如果把计算栅格叠加在细分度栅格上,则所述某一单元会包含指定节点的位置。
6.按照权利要求5所述的方法,还包括通过求取形成指定节点间连线的两个节点处的渗透率的调和平均值,计算计算栅格的指定节点间连线的渗透率。
7.按照权利要求1所述的方法,其中在步骤(b)中填充的渗透率被分配给计算栅格的指定节点间连线的中点,分配给指定节点间连线的渗透率对应于细分度栅格的某一单元的预定渗透率,如果把计算栅格叠加在细分度栅格上,则所述某一单元会包含所述节点间连线的中点。
8.按照权利要求1所述的方法,其中通过只利用关于落入计算栅格预定子域内的计算栅格节点间连线计算的通量和压强梯度,确定步骤(d)中节点间平均通量和平均压强梯度的计算。
9.按照权利要求1所述的方法,其中通过只利用在步骤(c)中,关于计算栅格的平行于指定节点间连线的节点间连线计算的通量和压强梯度,确定步骤(d)中关于粗分度栅格的指定节点间连线的平均通量和平均压强梯度的计算。
10.按照权利要求1所述的方法,其中关于粗分度栅格的指定节点,确定在步骤(e)中计算的渗透率。
11.按照权利要求1所述的方法,其中通过计算在步骤(d)中,关于指定的节点间连线计算的平均通量与平均压强梯度的比值,关于粗分度栅格的指定节点间连线,确定在步骤(e)中计算的渗透率。
12.按照权利要求1所述的方法,其中细分度栅格的单元是结构化的。
13.按照权利要求1所述的方法,其中粗分度栅格是PEBI栅格。
14.按照权利要求1所述的方法,其中粗分度栅格和计算栅格都是PEBI栅格。
15.按照权利要求1所述的方法,其中粗分度栅格的节点间连线形成Delaunay三角形,并且在步骤(a)中产生的计算栅格含有相似的较小的Delaunay三角形,其数目等于n2与大分度栅格的Delaunay三角形数目的乘积,这里n是用于产生计算栅格的预定整数细化比。
16.按照权利要求1所述的方法,其中所有单元都是三维的。
17.按照权利要求16所述的方法,其中粗分度栅格和计算栅格都是面内非结构化,垂直方向结构化的栅格。
18.按照权利要求1所述的方法,还包括利用在步骤(e)中计算的渗透率,确定粗分度栅格的节点间连线传输率。
19.按照权利要求1所述的方法,其中步骤(c)的流动方程是单相、定态方程。
20.一种利用代表地下地质域的第二栅格的各个单元的预定渗透率,估计具有所述地下地质域的许多单元的第一栅格的各个单元的渗透率的方法,所述第二栅格含有数目远远多于第一栅格的大量单元,所述方法包括:
(a)构造代表所述地下地质域的非结构化第三栅格,所述第三栅格包含数目与第二栅格近似相同或更大的大量单元,第一、第二及第三栅格的每个单元都具有一个节点,并且相邻单元的两个节点间的各个链接是节点间连线,第三栅格的所有节点间连线平行于第一栅格的节点间连线;
(b)把渗透率分配给第三栅格的每个节点,所述渗透率对应于第二栅格的含有第三栅格所述节点的位置的单元的渗透率;
(c)为第三栅格***的各个单元,导出单相、定态压强方程;
(d)求解所述压强方程,并计算第三栅格的所有节点间连线的通量和压强梯度;
(e)利用第三栅格的节点连线,计算第一栅格的指定连线的估计渗透率;
(f)对第一栅格的所有连线,重复步骤(e)。
21.按照权利要求20所述的方法,还包括借助辅助步骤,在步骤(e)中计算渗透率,所述辅助步骤确定与第一栅格的指定栅格节点间连线相关的子域内的平均通量和平均压强梯度,并计算平均通量与平均压强梯度的比值,从而获得第一栅格的指定节点间连线的渗透率。
22.按照权利要求20所述的方法,其中步骤(e)的渗透率计算只使用第三栅格的与第一栅格的节点间连线平行的节点间连线。
23.按照权利要求20所述的方法,其中所有栅格单元都是三维的。
24.按照权利要求23所述的方法,其中第二栅格和第三栅格都是面内非结构化,垂直方向结构化的栅格。
25.一种利用与代表多孔介质中的流体流动的小分度栅格的单元相关的预定渗透率,估计与同样代表所述多孔介质中的流体流动的大分度栅格的单元相关的渗透率的方法,大分度栅格的每个单元具有一个节点,并且每个节点与相邻节点链接,从而形成节点间连线,这种连线构成Delaunay三角形,所述方法包括下述步骤:
(a)通过把大分度栅格的各个Delaunay三角形分成许多相似的较小Delaunay三角形,构造计算栅格,这种较小的Delaunay三角形的侧边为计算栅格的节点间连线,并且大分度栅格和计算栅格的节点间连线相互对齐;
(b)把和小分度栅格的预定渗透率相对应的渗透率分配给计算栅格;
(c)为计算栅格的各个单元,导出单相、定态压强方程,求解所述压强方程,并计算计算栅格的所有节点间连线的通量和压强梯度;
(d)利用在步骤(c)中计算的通量和压强梯度,计算大分度栅格的各个节点间连线的平均通量和平均压强梯度;和
(e)利用在步骤(d)中计算的平均通量和平均压强梯度,计算与大分度栅格的指定节点间连线相关的渗透率。
26.按照权利要求25所述的方法,其中计算栅格的节点数和小分度栅格的单元数近似相等。
27.按照权利要求25所述的方法,其中计算栅格的节点数小于小分度栅格的单元数。
28.按照权利要求25所述的方法,其中相似的较小Delaunay三角形的数目为n2与大分度栅格的Delaunay三角形数目的乘积,n是预定的整数细化比。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14070099P | 1999-06-24 | 1999-06-24 | |
US60/140,700 | 1999-06-24 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1358289A true CN1358289A (zh) | 2002-07-10 |
Family
ID=22492429
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN00809383.0A Pending CN1358289A (zh) | 1999-06-24 | 2000-06-21 | 按比例增大非结构化栅格(grid)的渗透率的方法 |
Country Status (10)
Country | Link |
---|---|
US (1) | US6826520B1 (zh) |
EP (1) | EP1212696A1 (zh) |
CN (1) | CN1358289A (zh) |
AU (1) | AU763537B2 (zh) |
BR (1) | BR0011886A (zh) |
CA (1) | CA2375506C (zh) |
EA (1) | EA003438B1 (zh) |
MX (1) | MXPA01012736A (zh) |
NO (1) | NO321064B1 (zh) |
WO (1) | WO2000079423A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101903803B (zh) * | 2007-12-14 | 2013-05-08 | 埃克森美孚上游研究公司 | 在非结构化栅格上对地下过程进行建模 |
CN105388524A (zh) * | 2014-09-04 | 2016-03-09 | 中国石油化工股份有限公司 | 一种以三角形表面网格为限定条件的复杂地质体精细描述方法 |
CN111622733A (zh) * | 2020-05-25 | 2020-09-04 | 中国石油天然气股份有限公司 | 一种致密气藏区域部署加密井的方法 |
CN112157911A (zh) * | 2020-09-11 | 2021-01-01 | 华侨大学 | 一种自支撑的3d打印梯度弹性多孔材料微结构设计方法 |
CN114371116A (zh) * | 2021-12-30 | 2022-04-19 | 北京红山信息科技研究院有限公司 | 一种基于渗透率的路测质量评估方法 |
Families Citing this family (84)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8428923B2 (en) | 1999-04-29 | 2013-04-23 | Schlumberger Technology Corporation | Method system and program storage device for simulating a multilayer reservoir and partially active elements in a hydraulic fracturing simulator |
US7509245B2 (en) * | 1999-04-29 | 2009-03-24 | Schlumberger Technology Corporation | Method system and program storage device for simulating a multilayer reservoir and partially active elements in a hydraulic fracturing simulator |
US6370491B1 (en) * | 2000-04-04 | 2002-04-09 | Conoco, Inc. | Method of modeling of faulting and fracturing in the earth |
US7369973B2 (en) * | 2000-06-29 | 2008-05-06 | Object Reservoir, Inc. | Method and system for representing reservoir systems |
US7761270B2 (en) * | 2000-12-29 | 2010-07-20 | Exxonmobil Upstream Research Co. | Computer system and method having a facility management logic architecture |
US7277836B2 (en) * | 2000-12-29 | 2007-10-02 | Exxonmobil Upstream Research Company | Computer system and method having a facility network architecture |
FR2842321B1 (fr) * | 2002-07-11 | 2008-12-05 | Inst Francais Du Petrole | Methode pour contraindre un champ de permeabilite heterogene representant un reservoir souterrain par des donnees dynamiques |
US7283941B2 (en) * | 2001-11-13 | 2007-10-16 | Swanson Consulting Services, Inc. | Computer system and method for modeling fluid depletion |
US7496488B2 (en) | 2003-03-06 | 2009-02-24 | Schlumberger Technology Company | Multi-scale finite-volume method for use in subsurface flow simulation |
US6823297B2 (en) | 2003-03-06 | 2004-11-23 | Chevron U.S.A. Inc. | Multi-scale finite-volume method for use in subsurface flow simulation |
US7096122B2 (en) * | 2003-07-22 | 2006-08-22 | Dianli Han | Method for producing full field radial grid for hydrocarbon reservoir simulation |
CN1898640A (zh) | 2004-01-30 | 2007-01-17 | 埃克森美孚上游研究公司 | 储层评价方法 |
CN1898675A (zh) | 2004-01-30 | 2007-01-17 | 埃克森美孚上游研究公司 | 储层模型建立方法 |
US7379852B2 (en) * | 2004-02-18 | 2008-05-27 | Chevron U.S.A. Inc. | N-phase interface tracking method utilizing unique enumeration of microgrid cells |
FR2869116B1 (fr) * | 2004-04-14 | 2006-06-09 | Inst Francais Du Petrole | Methode pour construire un modele geomecanique d'une zone souterraine destine a etre couple a un modele de reservoir |
CN100489558C (zh) | 2004-06-07 | 2009-05-20 | 埃克森美孚上游研究公司 | 用于求解隐式储层仿真矩阵的方法 |
US7526418B2 (en) * | 2004-08-12 | 2009-04-28 | Saudi Arabian Oil Company | Highly-parallel, implicit compositional reservoir simulator for multi-million-cell models |
US7596480B2 (en) * | 2005-04-14 | 2009-09-29 | Saudi Arabian Oil Company | Solution method and apparatus for large-scale simulation of layered formations |
FR2886742B1 (fr) * | 2005-06-02 | 2007-07-27 | Inst Francais Du Petrole | Methode de changement d'echelle des permeabilites absolues pour construire un modele de simulation d'ecoulement |
US7716028B2 (en) * | 2006-05-24 | 2010-05-11 | Schlumberger Technology Corporation | Method for modeling a reservoir using a 3D wettability map generated from a wettability logging tool |
WO2007149766A2 (en) | 2006-06-18 | 2007-12-27 | Chevron U.S.A. Inc. | Reservoir simulation using a multi-scale finite volume including black oil modeling |
US8078437B2 (en) * | 2006-07-07 | 2011-12-13 | Exxonmobil Upstream Research Company | Upscaling reservoir models by reusing flow solutions from geologic models |
WO2008020906A2 (en) * | 2006-08-14 | 2008-02-21 | Exxonmobil Upstream Research Company | Enriched multi-point flux approximation |
US7933758B2 (en) * | 2006-12-05 | 2011-04-26 | Conocophillips Company | Method and apparatus for geomodel uplayering |
CN101809627A (zh) * | 2007-06-01 | 2010-08-18 | 埃克森美孚上游研究公司 | 平面中约束沃罗努网格的生成 |
KR20090005638A (ko) * | 2007-07-09 | 2009-01-14 | 한국과학기술원 | 불일치 요소망을 해결하기 위한 다절점 천이 유한요소모델링 방법 및 기록매체 |
CN101896690B (zh) | 2007-12-13 | 2015-02-18 | 埃克森美孚上游研究公司 | 使用非结构化网格的储层模拟上的并行自适应数据分区 |
US8190414B2 (en) * | 2008-03-26 | 2012-05-29 | Exxonmobil Upstream Research Company | Modeling of hydrocarbon reservoirs containing subsurface features |
WO2009120409A1 (en) * | 2008-03-28 | 2009-10-01 | Exxonmobil Upstream Research Company | Computing a consistent velocity vector field from a set of fluxes |
RU2440591C2 (ru) | 2008-04-10 | 2012-01-20 | Шлюмбергер Текнолоджи Б.В. | Способ получения характеристик геологической формации, пересекаемой скважиной |
US8725477B2 (en) | 2008-04-10 | 2014-05-13 | Schlumberger Technology Corporation | Method to generate numerical pseudocores using borehole images, digital rock samples, and multi-point statistics |
CN102165449B (zh) * | 2008-05-16 | 2014-01-29 | 雪佛龙美国公司 | 用于多孔介质中的多相流的多尺度方法 |
US8095349B2 (en) * | 2008-05-30 | 2012-01-10 | Kelkar And Associates, Inc. | Dynamic updating of simulation models |
CA2724589C (en) * | 2008-06-09 | 2013-12-10 | Landmark Graphics Corporation | Distribution of properties in a 3d volumetric model using a maximum continuity field |
CA2720638C (en) * | 2008-06-16 | 2016-01-26 | Schlumberger Canada Limited | Streamline flow simulation of a model that provides a representation of fracture corridors |
WO2009155244A1 (en) * | 2008-06-16 | 2009-12-23 | Sunoco, Inc | Extensible spunbonded non-woven fabrics |
CA2729770C (en) * | 2008-07-03 | 2016-12-13 | Chevron U.S.A. Inc. | Multi-scale finite volume method for reservoir simulation |
EP2361414A2 (en) * | 2008-10-09 | 2011-08-31 | Chevron U.S.A. Inc. | Iterative multi-scale method for flow in porous media |
BRPI0922214A2 (pt) * | 2008-12-03 | 2015-12-29 | Chevron Usa Inc | método implementado por computador, e, sistema para gerar uma grade usada para construir um modelo de simulação de um reservatório de subsuperfície fraturado |
GB2477901A (en) * | 2008-12-03 | 2011-08-17 | Chevron Usa Inc | System and method for predicting fluid flow characteristics within fractured subsurface reservoirs |
US8311788B2 (en) | 2009-07-01 | 2012-11-13 | Schlumberger Technology Corporation | Method to quantify discrete pore shapes, volumes, and surface areas using confocal profilometry |
US9134457B2 (en) | 2009-04-08 | 2015-09-15 | Schlumberger Technology Corporation | Multiscale digital rock modeling for reservoir simulation |
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 |
BRPI0924623A2 (pt) * | 2009-06-30 | 2016-03-01 | Prad Res & Dev Ltd | método numérico para cálculo de transporte de calor de massa, químico e elétrico para um meio tridimensional poroso. |
CA2771865C (en) * | 2009-10-23 | 2016-04-05 | Exxonmobil Upstream Research Company | Method for optimization with gradient information |
WO2011056666A2 (en) * | 2009-10-28 | 2011-05-12 | Chevron U.S.A. Inc. | Multiscale finite volume method for reservoir simulation |
BR112012009045B1 (pt) * | 2009-11-12 | 2020-12-08 | Exxonmobil Upstream Research Company | métodos para gerar uma grade de simulação tridimensional para um modelo de reservatório, aparelho de gradeamento de simulação e simulador de reservatório |
CN102667804A (zh) * | 2009-11-23 | 2012-09-12 | 埃克森美孚上游研究公司 | 使用均匀化混合有限元建模地质性质的方法和*** |
CA2791642C (en) | 2010-04-06 | 2018-09-11 | Exxonmobil Upstream Research Company | Hierarchical modeling of physical systems and their uncertainties |
BR112012025995A2 (pt) | 2010-04-30 | 2016-06-28 | Exxonmobil Upstream Res Co | método e sistema para simulação de fluxo de volume finito |
US8775142B2 (en) | 2010-05-14 | 2014-07-08 | Conocophillips Company | Stochastic downscaling algorithm and applications to geological model downscaling |
US8463586B2 (en) | 2010-06-22 | 2013-06-11 | Saudi Arabian Oil Company | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5D unstructured grids |
US9187984B2 (en) * | 2010-07-29 | 2015-11-17 | Exxonmobil Upstream Research Company | Methods and systems for machine-learning based simulation of flow |
EP2599029A4 (en) | 2010-07-29 | 2014-01-08 | Exxonmobil Upstream Res Co | METHODS AND SYSTEMS FOR AUTOMATIC LEARNING FLOW SIMULATION |
WO2012015521A1 (en) * | 2010-07-29 | 2012-02-02 | Exxonmobil Upstream Research Company | Method and system for reservoir modeling |
US10198535B2 (en) | 2010-07-29 | 2019-02-05 | Exxonmobil Upstream Research Company | Methods and systems for machine-learning based simulation of flow |
EP2601642B1 (en) | 2010-08-04 | 2018-06-13 | Exxonmobil Upstream Research Company | System and method for summarizing data on an unstructured grid |
CN103210397B (zh) * | 2010-08-09 | 2015-12-09 | 科诺科菲利浦公司 | 具有保持传递率的油藏粗化方法 |
US8433551B2 (en) | 2010-11-29 | 2013-04-30 | Saudi Arabian Oil Company | Machine, computer program product and method to carry out parallel reservoir simulation |
US8386227B2 (en) | 2010-09-07 | 2013-02-26 | Saudi Arabian Oil Company | Machine, computer program product and method to generate unstructured grids and carry out parallel reservoir simulation |
BR112013002114A2 (pt) | 2010-09-20 | 2016-05-17 | Exxonmobil Upstream Res Co | formulações flexíveis e adaptáveis para simulações de reservatório complexas |
US10534871B2 (en) | 2011-03-09 | 2020-01-14 | Schlumberger Technology Corporation | Method and systems for reservoir modeling, evaluation and simulation |
FR2976099B1 (fr) * | 2011-06-01 | 2013-05-17 | IFP Energies Nouvelles | Methode pour construire un maillage d'un reseau de fractures a partir de diagramme de voronoi |
US9489176B2 (en) | 2011-09-15 | 2016-11-08 | Exxonmobil Upstream Research Company | Optimized matrix and vector operations in instruction limited algorithms that perform EOS calculations |
US9753180B2 (en) | 2012-03-28 | 2017-09-05 | Exxonmobil Upstream Research Company | Method for multiphase flow upscaling |
RU2582482C1 (ru) * | 2012-03-30 | 2016-04-27 | Лэндмарк Графикс Корпорейшн | Система и способ для автоматического локального измельчения сетки в системе моделирования резервуара |
WO2014051903A1 (en) | 2012-09-28 | 2014-04-03 | Exxonmobil Upstream Research Company | Fault removal in geological models |
US20150113379A1 (en) * | 2013-10-23 | 2015-04-23 | Schlumberger Technology Corporation | Representation of unstructured grids |
US10319143B2 (en) | 2014-07-30 | 2019-06-11 | Exxonmobil Upstream Research Company | Volumetric grid generation in a domain with heterogeneous material properties |
CN104331931B (zh) * | 2014-10-10 | 2017-05-17 | 北京航空航天大学 | 一种面向油藏数值模拟的地质断裂层面pebi网格处理方法 |
AU2015339884B2 (en) | 2014-10-31 | 2018-03-15 | Exxonmobil Upstream Research Company | Handling domain discontinuity in a subsurface grid model with the help of grid optimization techniques |
EP3213125A1 (en) | 2014-10-31 | 2017-09-06 | Exxonmobil Upstream Research Company Corp-urc-e2. 4A.296 | Methods to handle discontinuity in constructing design space for faulted subsurface model using moving least squares |
GB2566853B (en) * | 2016-06-28 | 2022-03-30 | Geoquest Systems Bv | Parallel multiscale reservoir simulation |
US11542784B2 (en) * | 2016-11-08 | 2023-01-03 | Landmark Graphics Corporation | Diffusion flux inclusion for a reservoir simulation for hydrocarbon recovery |
CA3043231C (en) | 2016-12-23 | 2022-06-14 | Exxonmobil Upstream Research Company | Method and system for stable and efficient reservoir simulation using stability proxies |
US11163923B2 (en) | 2017-02-14 | 2021-11-02 | Landmark Graphics Corporation | Automated upscaling of relative permeability and capillary pressure in multi-porosity systems |
CN107194070A (zh) * | 2017-05-23 | 2017-09-22 | 西南石油大学 | 一种基于pebi网格加密来确定水平井裂缝参数的方法 |
CN109632604B (zh) * | 2019-01-04 | 2021-06-15 | 中国海洋石油集团有限公司 | 一种孔隙尺度到岩心尺度聚合物驱相对渗透率粗化方法 |
EP3956699B1 (en) * | 2019-04-16 | 2023-03-29 | TotalEnergies OneTech | A method for upscaling of relative permeability of the phase of a fluid |
CN110006793B (zh) * | 2019-05-08 | 2020-02-07 | 中南大学 | 用于振动荷载下颗粒材料运动特性研究的试验装置和方法 |
US20220244424A1 (en) * | 2019-06-17 | 2022-08-04 | Totalenergies Se | Geological Grid Analysis |
CN110579433B (zh) * | 2019-09-30 | 2022-02-01 | 中国科学院力学研究所 | 一种颗粒样品两级渗透率的获取方法 |
US11402315B2 (en) | 2020-07-06 | 2022-08-02 | Landmark Graphics Corporation | Estimating relative permeability and capillary pressures of a geological formation based on multiphase upscaling |
CN112014291B (zh) * | 2020-08-19 | 2021-06-29 | 中国水利水电科学研究院 | 一种土样渗透变形特征判别方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4969130A (en) | 1989-09-29 | 1990-11-06 | Scientific Software Intercomp, Inc. | System for monitoring the changes in fluid content of a petroleum reservoir |
US5740342A (en) | 1995-04-05 | 1998-04-14 | Western Atlas International, Inc. | Method for generating a three-dimensional, locally-unstructured hybrid grid for sloping faults |
US5757663A (en) * | 1995-09-26 | 1998-05-26 | Atlantic Richfield Company | Hydrocarbon reservoir connectivity tool using cells and pay indicators |
US5710726A (en) * | 1995-10-10 | 1998-01-20 | Atlantic Richfield Company | Semi-compositional simulation of hydrocarbon reservoirs |
US5886702A (en) | 1996-10-16 | 1999-03-23 | Real-Time Geometry Corporation | System and method for computer modeling of 3D objects or surfaces by mesh constructions having optimal quality characteristics and dynamic resolution capabilities |
US6018497A (en) * | 1997-02-27 | 2000-01-25 | Geoquest | Method and apparatus for generating more accurate earth formation grid cell property information for use by a simulator to display more accurate simulation results of the formation near a wellbore |
US6106561A (en) * | 1997-06-23 | 2000-08-22 | Schlumberger Technology Corporation | Simulation gridding method and apparatus including a structured areal gridder adapted for use by a reservoir simulator |
FR2765708B1 (fr) * | 1997-07-04 | 1999-09-10 | Inst Francais Du Petrole | Methode pour determiner des parametres hydrauliques representatifs a grande echelle d'un milieu fissure |
US6052520A (en) * | 1998-02-10 | 2000-04-18 | Exxon Production Research Company | Process for predicting behavior of a subterranean formation |
US6230101B1 (en) * | 1999-06-03 | 2001-05-08 | Schlumberger Technology Corporation | Simulation method and apparatus |
-
2000
- 2000-06-13 US US09/593,121 patent/US6826520B1/en not_active Expired - Lifetime
- 2000-06-21 AU AU60534/00A patent/AU763537B2/en not_active Ceased
- 2000-06-21 CN CN00809383.0A patent/CN1358289A/zh active Pending
- 2000-06-21 MX MXPA01012736A patent/MXPA01012736A/es not_active Application Discontinuation
- 2000-06-21 EP EP00946837A patent/EP1212696A1/en not_active Withdrawn
- 2000-06-21 WO PCT/US2000/017101 patent/WO2000079423A1/en active IP Right Grant
- 2000-06-21 BR BR0011886-9A patent/BR0011886A/pt not_active Application Discontinuation
- 2000-06-21 CA CA2375506A patent/CA2375506C/en not_active Expired - Lifetime
- 2000-06-21 EA EA200200077A patent/EA003438B1/ru not_active IP Right Cessation
-
2001
- 2001-12-21 NO NO20016321A patent/NO321064B1/no not_active IP Right Cessation
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101903803B (zh) * | 2007-12-14 | 2013-05-08 | 埃克森美孚上游研究公司 | 在非结构化栅格上对地下过程进行建模 |
CN105388524A (zh) * | 2014-09-04 | 2016-03-09 | 中国石油化工股份有限公司 | 一种以三角形表面网格为限定条件的复杂地质体精细描述方法 |
CN111622733A (zh) * | 2020-05-25 | 2020-09-04 | 中国石油天然气股份有限公司 | 一种致密气藏区域部署加密井的方法 |
CN112157911A (zh) * | 2020-09-11 | 2021-01-01 | 华侨大学 | 一种自支撑的3d打印梯度弹性多孔材料微结构设计方法 |
CN114371116A (zh) * | 2021-12-30 | 2022-04-19 | 北京红山信息科技研究院有限公司 | 一种基于渗透率的路测质量评估方法 |
Also Published As
Publication number | Publication date |
---|---|
AU763537B2 (en) | 2003-07-24 |
CA2375506C (en) | 2010-08-03 |
NO20016321D0 (no) | 2001-12-21 |
EA200200077A1 (ru) | 2002-06-27 |
BR0011886A (pt) | 2002-03-05 |
EP1212696A1 (en) | 2002-06-12 |
AU6053400A (en) | 2001-01-09 |
NO321064B1 (no) | 2006-03-06 |
US6826520B1 (en) | 2004-11-30 |
WO2000079423A1 (en) | 2000-12-28 |
CA2375506A1 (en) | 2000-12-28 |
EA003438B1 (ru) | 2003-04-24 |
NO20016321L (no) | 2001-12-21 |
MXPA01012736A (es) | 2002-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1358289A (zh) | 按比例增大非结构化栅格(grid)的渗透率的方法 | |
Berre et al. | Flow in fractured porous media: A review of conceptual models and discretization approaches | |
Pollock | User guide for MODPATH Version 7—A particle-tracking model for MODFLOW | |
Karimi-Fard et al. | A general gridding, discretization, and coarsening methodology for modeling flow in porous formations with discrete geological features | |
Bogdanov et al. | Effective permeability of fractured porous media in steady state flow | |
Schaap et al. | Comparison of pressure‐saturation characteristics derived from computed tomography and lattice Boltzmann simulations | |
Lough et al. | An efficient boundary integral formulation for flow through fractured porous media | |
Hyman | Flow channeling in fracture networks: characterizing the effect of density on preferential flow path formation | |
Li et al. | Optimal uplayering for scaleup of multimillion-cell geologic models | |
Liu et al. | Multi-block lattice Boltzmann simulations of solute transport in shallow water flows | |
Bonduà et al. | 3D Voronoi grid dedicated software for modeling gas migration in deep layered sedimentary formations with TOUGH2-TMGAS | |
Bogdanov et al. | Pressure drawdown well tests in fractured porous media | |
Steward et al. | The Analytic Element Method for rectangular gridded domains, benchmark comparisons and application to the High Plains Aquifer | |
Khirevich et al. | Three-dimensional simulation of tracer transport dynamics in formations with high-permeability channels or fractures: Estimation of oil saturation | |
Pan | Wingridder-an interactive grid generator for TOUGH2 | |
CN107169227A (zh) | 一种分段压裂水平井的粗网格模拟方法及*** | |
LaForce et al. | Voronoi Meshing to Accurately Capture Geological Structure in Subsurface Simulations | |
Jamal et al. | Darcy’s model with optimized permeability distribution for the simulation of Stokes flow and contaminant transport in karst aquifers | |
Masciopinto | Particles’ transport in a single fracture under variable flow regimes | |
Mustapha | Finite element mesh for complex flow simulation | |
Haegland | Streamline methods with application to flow and transport in fractured media | |
Pierre et al. | Comparison of various discretization schemes for simulation of large field case reservoirs using unstructured grids | |
Bonduà et al. | 3D Voronoi pre-and post-processing tools for using the TOUGH2 family of numerical simulator for hydrocarbon gas migration | |
Sonier et al. | Full-Field Gas Storage Simulation Using a Control-Volume Finite-Element Model | |
Chattopadhyay et al. | Fractal characters of porous media and flow analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |