CN111936841B - 确定多孔介质的流体流特点 - Google Patents
确定多孔介质的流体流特点 Download PDFInfo
- Publication number
- CN111936841B CN111936841B CN201980012848.7A CN201980012848A CN111936841B CN 111936841 B CN111936841 B CN 111936841B CN 201980012848 A CN201980012848 A CN 201980012848A CN 111936841 B CN111936841 B CN 111936841B
- Authority
- CN
- China
- Prior art keywords
- wettability
- local
- fluid
- representation
- rock sample
- 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
- 239000012530 fluid Substances 0.000 title claims abstract description 114
- 239000011435 rock Substances 0.000 claims abstract description 89
- 238000000034 method Methods 0.000 claims abstract description 85
- 238000012986 modification Methods 0.000 claims abstract description 60
- 230000004048 modification Effects 0.000 claims abstract description 60
- 230000032683 aging Effects 0.000 claims abstract description 33
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 30
- 239000011148 porous material Substances 0.000 claims abstract description 25
- 238000003860 storage Methods 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 238000004088 simulation Methods 0.000 claims description 50
- 230000008859 change Effects 0.000 claims description 21
- 238000000926 separation method Methods 0.000 claims description 21
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 19
- 239000011707 mineral Substances 0.000 claims description 19
- 238000009826 distribution Methods 0.000 claims description 17
- 239000000203 mixture Substances 0.000 claims description 13
- 230000006399 behavior Effects 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 5
- 238000005094 computer simulation Methods 0.000 claims description 4
- 239000011800 void material Substances 0.000 claims description 4
- 239000003129 oil well Substances 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 44
- 239000002245 particle Substances 0.000 abstract description 12
- 238000004590 computer program Methods 0.000 abstract description 9
- 238000009736 wetting Methods 0.000 description 11
- 230000035699 permeability Effects 0.000 description 9
- 239000000126 substance Substances 0.000 description 8
- 238000006073 displacement reaction Methods 0.000 description 7
- 239000012528 membrane Substances 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 239000004215 Carbon black (E152) Substances 0.000 description 4
- 239000012267 brine Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 229930195733 hydrocarbon Natural products 0.000 description 4
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 238000009533 lab test Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012797 qualification Methods 0.000 description 2
- 238000013515 script Methods 0.000 description 2
- 230000002269 spontaneous effect Effects 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000000356 contaminant Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 125000001183 hydrocarbyl group Chemical group 0.000 description 1
- 238000005213 imbibition Methods 0.000 description 1
- 238000011005 laboratory method Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
-
- 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
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/241—Earth materials for hydrocarbon content
-
- 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
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Fluid Mechanics (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Medicinal Chemistry (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Dispersion Chemistry (AREA)
- Radiology & Medical Imaging (AREA)
- Pulmonology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
公开了包括方法、装置和计算机程序产品的技术。这些技术包括在计算机存储介质上编码的用于确定润湿性的计算机指令。这些技术使用数字老化计算过程来提供在至少两个流体存在的情况下物理岩石样本的润湿性更改的表示。这些技术包括检索物理岩石样本的表示,该表示包括与物理岩石样本对应的孔隙空间和颗粒空间数据,计算孔隙空间中表面位置的局部曲率,根据计算出的局部曲率确定水膜破裂是否会发生,以及基于对于水膜破裂的确定来对物理岩石的润湿性进行分类。
Description
相关申请的交叉引用
本申请要求于2018年1月10日提交的标题为“Digital Rock Procedures forDetermining Wettability”的美国临时申请No.62/615,701的权益,该申请通过引用并入本文。
背景技术
两相相对渗透率(two-phase relative permeability)描述了一个流体在存在另一个流体的情况下移动通过多孔介质移动的难易程度。它是烃储层岩石(hydrocarbonreservoir rocks)的重要特点,并且是油气储层建模和模拟活动的关键输入。
术语“润湿性(wettability)”被用于测量:与相同孔隙空间中的第二流体相比,流体粘附于(或润湿)表面的趋势,这是一个流体在存在其它不混溶的流体的情况下在固体表面上扩散或粘附的趋势。在石油环境中,润湿性是储层岩石表面优先接触多相或两相流体系统中特定流体的趋势。
相对渗透率和通过多孔介质的多相流一般取决于流体-流体-岩石系统的各种特点,包括岩石表面特性、每个流体的物理特性和流动条件。一个流特点是无量纲的“毛细管数”(Ca,Capillary number),它表示粘性力(viscous force)与毛细管力(capillaryforce)之比。“毛细管数”(Ca)通常被定义为表观流体速度(superficial fluid velocity)乘以参考流体(reference fluid)的粘度(viscosity)除以流体-流体界面张力(fluid-fluid interfacial tension)。另一个特点是“润湿性”,它表示岩石表面比一个流体更喜欢另一个流体,并表现为可测量的特性,称为接触角(contact angle)。润湿性通常分为水润湿(平均接触角0-90度)、中性润湿(~90度)、油润湿(90-180度)或其某种变型(诸如强水润湿或弱油润湿)。
润湿性对多孔材料中的多相流动行为(诸如流动通过储层岩石的油和水、油和气、或水和气)有强烈影响。储层岩石(不同于人造材料)倾向于具有在空间上变化的润湿性,即,接触角在岩石的孔隙空间内的表面上随位置的不同而变化。接触角的分布是含烃岩石的矿物成分和地质历史以及表面纹理、与岩石接触的流体(例如,水、油)的化学成分等的结果。在与油接触之前,岩石几乎总是自然是水润湿的;随着油在地质时期内侵入孔隙空间,由于某些物质在油相中的粘附,初始的水润湿特点会在某些位置发生更改。润湿性更改的机制取决于各种局部系统特性,诸如压力、温度、矿物质类型和流体成分。
在钻井过程中从油/气储层中提取岩石样本时,它们常常被泥浆、其它流体和化学物质污染。在对这些岩石样本进行实验室实验以获得两相相对渗透率曲线或毛细管压(capillary pressure)曲线时,常见的实验室测试规程是首先清洁岩石,然后尝试恢复岩石的自然(提取前的)润湿性。尝试重新建立代表地下储层状况的流体相(fluid phase)分布和润湿性分布(例如,接触角)的过程称为“老化(aging)”。典型的老化规程将涉及清洁岩石样本中的污染物并用诱发水润湿条件的化学物质处理岩石样本。接下来,将样本浸入盐水(brine)中,并在称为“排水(drainage)”的规程中,使用估计的温度和压力的储层条件将油推入(因此将盐水推出)。然后将这个系统“老化”一段时间,例如4周,从而允许润湿性发生更改,大概与实际地下岩石形成中发生的变化足够相似,从而恢复真实的润湿条件。
发明内容
在一个方面,一种用于在至少两个流体存在的情况下确定物理岩石样本的润湿性更改的计算机实现的方法,包括由计算系统检索物理岩石样本的表示,该表示包括与物理岩石样本对应的孔隙(pore)空间和颗粒(grain)空间数据;由计算系统计算孔隙空间中的表面位置的局部曲率;以及基于对局部曲率的确定,由计算机系统对与物理岩石对应的表面位置的润湿性进行分类。
以下是这个方面内的一些附加特征。计算局部曲率包括由计算系统计算孔隙空间中的表面位置的局部毛细管压和局部分离压(disjoining pressure),并且该方法还包括根据计算出的局部曲率确定是否会发生水膜破裂(water-film breakage);对于表面位置中的一些,通过将毛细管压与基于局部曲率的压力和局部分离压之和进行比较,确定是否会发生水膜破裂;以及对于表面位置中的具有水膜破裂的至少一些表面位置,基于局部表面特性、局部流体化学成分和老化时间来确定润湿性更改的程度。
生产还包括将润湿性模拟所需的输入参数指派给计算机模拟环境。岩石样本的表示被建模为三角网格(triangulated mesh),并且流体由体元(voxel)表示,体元被指派有属性。可以将体元切割以使流体与三角网格相符。该方法还包括根据用户指定的条件执行排水模拟以建立流体相分布。用户指定的条件选自包括实验室条件和储层条件的代表的组。执行排水模拟还包括使用流体计算技术来模拟排水,以针对用户指定的条件建立通过岩石样本表示的流体相分布。执行排水模拟还包括在岩石样本表示中的每个表面处执行局部曲率计算,以及应用局部曲率计算的结果来计算每个表面位置处的局部毛细管压。
在执行确定润湿性更改的数值表示的迭代之后,该方法还包括允许流体相由于经修改的表面润湿性而重新分布,并且重复确定润湿性更改的数值表示。这以每个表面位置处计算出的接触角的形式提供了表面润湿性更改的数值预测以及所得到的润湿性空间分布。离散的表面元(surface element)和一个或多个流体元(fluid element)连通,并进行取决于矿物类型的局部特性、流体特性、流体成分、毛细管压和分离压的局部润湿性确定。润湿性更改被迭代地执行,然后进行流体相重新分布,润湿性更改被重复进行直到计算系统检测到收敛,所述收敛是通过后续重复没有显著改变所确定的接触角的值而测量出的。
其它方面包括有形地存储在非暂态计算机可读介质上的计算机程序产品和计算系统(诸如计算机系统、计算机服务器)。
上述方面中的一个或多个可以提供以下优点中的一个或多个。润湿性对多孔材料中的多相流动行为(诸如流过储层岩石的油和水、油和气、或水和气)有强烈影响,在钻井过程期间可以从油/气储层中提取储层岩石,被清洗以测量建模的物理、化学和几何特性。模拟的岩石样本经过模拟的老化和润湿性更改,而不是使样本经受常规的实验室规程处理,以获得两相相对渗透率曲线或毛细管压曲线。相较于适于上面提到的物理实验室方法的例如4周的“老化”时间,所公开的计算的“老化”和“润湿性更改”的过程可以在短得多的时间段内执行。
通过以下描述和权利要求,本发明的其它特征和优点将变得清楚。
附图说明
图1描绘了用于润湿性更改的模拟表示的系统。
图2描绘了示出数字老化的操作的流程图。
图3描绘了示出用于表示润湿性更改模拟的操作的流程图。
图4描绘了示出验证/校准润湿性更改模拟的操作的流程图。
图5、5A描绘了岩石样本中的孔的现有技术图,其示出了固定角度。
具体实施方式
参考图1,示出了用于模拟的系统10,以数值模拟代表地下储层条件的润湿性恢复或“老化”过程(即,“数字老化”)。在这个实施方式中的系统10基于客户端-服务器或基于云的体系架构,并且包括被实现为大规模并行计算系统12(独立的或基于云的)的服务器系统12和客户端系统14。服务器系统12包括存储器18、总线系统11、接口20(例如,用户接口/网络接口/显示或监视器接口等)和处理设备24。在存储器18中的是数字老化引擎32,其对物理岩石样本的数字表示进行操作,该数字岩石表示以数字方式表示物理材料(例如,物理岩石样本)的数字表示(数字岩石样本)的孔隙空间和颗粒空间。还在存储器中的是模拟引擎34,其模拟润湿性更改。
在一些实施例中,通过邻近气井或油井(例如,钻机37)的储层岩石的模拟多相流动行为发生。确定多相流动行为包括确定物理岩石样本的润湿性的更改。
物理岩石样本的数字表示可以是第三方应用,该第三方应用在与服务器12不同的系统上执行。系统10仅需要物理岩石样本32’的数字表示,以使数字老化引擎数字地准备物理岩石样本的数字表示。提供岩石样本的数字表示32’的一种方法是例如从岩石样本的微CT扫描生成的3D图像中获得该表示32’。
存储器18还存储由引擎32使用的参数,诸如通过将矿物类型33a指派给颗粒以确定每种矿物类型的表面特性而获得的颗粒表面特性,以及表面纹理和粗糙度特性。存储器18还存储参数,诸如流体特性33b,例如,每种预期流体(例如,水、气、油中的两种或更多种)的流体密度和粘度以及流体-流体界面张力特性。存储器18还存储用于特定矿物类型的参数,诸如流体的化学成分数据33c和流体成分的亲和度(affinity)数据33d。存储器18还存储与流体结合的每种矿物类型的分离压33e,以及由老化引擎32使用的选择的老化时间33f。此外,还存储储层压力和温度数据。被评估的矿物类型可以是在储层实际位置发现或预期的那些。
模拟引擎34包括用于建立岩石样本模拟环境的模块34a、用于执行排水模拟的模块34b以及用于计算孔隙空间中的表面的局部曲率的模块34c。系统10访问存储2D和/或3D网格、坐标系和库的数据储存库38,其可以被用于使用诸如计算流体动力学或所谓的晶格玻尔兹曼(Lattice Boltzmann)方法之类的任何众所周知的计算技术进行排水模拟。
现在参考图2,用于数字老化的过程40检索41岩石样本的孔隙空间和颗粒空间的数字表示,建立输入参数和条件(例如,初始和边界条件),通过向颗粒指派矿物类型并确定矿物类型的表面特性来确定42颗粒表面特性,确定44每个流体的流体密度和粘度以及流体-流体界面张力,确定46与流体结合的每个矿物类型的分离压,并通过执行润湿性更改60(图3)而在一定时间间隔内老化岩石的数字表示。在润湿性更改过程60的迭代之后,过程40相对于润湿性更改过程60的先前迭代来确定例如局部曲率、接触角等的改变量48,并将所述量与阈值改变量进行比较。颗粒表面特性包括表面纹理,通常称为表面“粗糙度”特性,其将对接触角产生影响。
如果改变是显著的,那么过程60在指定条件(图3)下使用后续的排水模拟来进行重复49a。如果改变不显著,那么过程40可以退出49b或建立52岩石样本模拟环境的后续迭代,从而指派后续的输入参数集。岩石样本模拟的后续迭代52可以被用于改变图2的任何元素,如由括号54所示的。
术语“润湿性”用于测量:与第二流体相比,一个流体粘附于(或润湿)表面的趋势。即,“润湿性”是一个流体在存在其它不混溶的液体的情况下在固体表面上扩散或粘附的趋势。在石油环境中,润湿性是储层岩石表面优先接触多相或两相流体系统中特定流体的趋势。润湿性可以受一个流体渗透到另一个流体中的量的影响,从而造成膜破裂(即,第一流体的膜保护表面不受第二流体侵害的能力的破坏)。下面讨论的润湿性更改过程被用于模拟一个流体随时间的递增以及表面和/或孔隙表面的递增而对另一个流体的影响。
现在参考图3,示出了用于模拟润湿性更改的过程60。在这个示例中,过程60建立62岩石样本模拟环境并指派如上面所指定的输入参数。岩石样本模拟环境包括在存储器中表示为节点集合32”的岩石样本的孔隙空间和颗粒空间的数字表示32’,节点集合32”与润湿性模拟所需的已指派的输入参数相关联。
例如,可以将岩石样本的微CT扫描生成的3D图像数字化并表示为例如存储器中的图结构,其中该图具有节点和叶子,例如,以逐像素(或像素的组)为基础。孔隙可以通过相对对比度水平来识别。该过程可以执行最佳路径算法(例如,使用Dijkstra方法)以找到孔隙长度。最佳路径可以被视为该孔隙的中心线。该过程识别中心线上的每个节点/像素的对比度水平,并且可以过滤对比度水平的高频值,该过程可以使用其它过滤和/或空间平滑技术来过滤高频值。节点可以用润湿性模拟所需的输入参数进行注释。
其它表示可以是体元。例如,在模拟之前,岩石模拟环境被建模为体元的集合。三角网格图像过程(计算机图形学中使用的一类多边形网格,其包括通过公共边或角连接的三角形的集合)和分割网格构造可以被用于在模拟环境中表示岩石样本的表面,并且流体由允许被“切割”的体元表示。体元被切割以使体元与网格相符。每个体元单元被指派有属性。
过程60执行排水模拟64以针对用户指定的条件(可以代表实验室和/或储层条件)建立流体相分布,这些条件包括储层的压力和温度。排水模拟是对岩石样本的部分迭代执行的。最常见的是,岩石样本最初是均匀水润湿的,但其它初始润湿性定义也是可能的。该过程可以访问例如3D网格和坐标系,它们与库一起被用于使用任何众所周知的计算技术(诸如计算流体动力学或所谓的晶格玻尔兹曼方法)来模拟排水。该模拟将通过岩石样本表示32’建立针对用户指定条件的流体相分布。
过程60在岩石样本表示32’中的每个表面(例如,外表面和内部孔隙表面)处执行局部曲率计算66。通过局部曲率计算,过程60计算68每个表面位置处的局部毛细管压,并通过将局部毛细管压与几何推导的、基于曲率的压力和分离压之和进行比较来确定70在该表面处是否可能发生水膜破裂。对于非常小的孔隙,由于非常大的曲率而引起的分离压效应将超过基于材料的分离压的效应。在这种情况下,只有孔隙的曲率才有意义。这种情况可能发生在粘土或页岩型地层中。
岩石表面的曲率是使用微分几何来计算的。一个特定的工具是“VTK”(可视化工具包,一种用于3D计算机图形的开源、可免费取得的软件系统)。可以使用其它工具等,也可以使用其它技术。Nira Dyn等人在“Optimizing 3D Triangulations Using DiscreteCurvature Analysis”中也讨论了微分几何(2001)(从https://www.inf.usi.ch/hormann/papers/Dyn.2001.OTU.pdf获得),其全部内容通过引用并入本文。过程60通过在排水模拟之后取得与非润湿相接触的位置以及适合于水膜破裂的表面位置的子集来计算71润湿性的空间分布,并基于局部流体的化学成分和老化时间确定72润湿性更改的程度,与较短的老化时间相比,较长的老化时间可能会导致更多的油润湿(更大的)接触角。然后,该过程可以将润湿性分类为水润湿(平均接触角0-90度)、中性润湿(~90度)、油润湿(90-180度)。
润湿性更改模拟可以是缓慢的动态过程,其发生在小的离散时间间隔内。在驱替(displacement)模拟期间的任何时间,每个表面元都被测试,以确定是否发生了例如接触角的更改并确定更改量。因此也捕捉了润湿性更改的空间膨胀的历史。在驱替模拟期间,在任何时间增量处,该过程测量所施加的局部压力并将其与基于曲率的局部压力和分离压之和进行比较。如果发生膜破裂,那么该过程将修改局部润湿性分类。在驱替模拟期间,在任何时间增量处,该过程都测量第二流体(例如,油)何时与表面元接触,并相应地修改润湿性分类。该过程通过注入流体或在排水期间毛细管压的增加来继续进行驱替模拟,例如相对渗透率。如上所述,该过程在最终压力表示储层压力的情况下执行,或者可以被细分为若干次排水压力迭代,其中在每个步骤之后执行润湿性更改,然后进行压力提升。之后,过程60可以允许流体相由于经修改的表面润湿性而重新分布,并重复过程60直到润湿性没有改变或存在不显著的(经验的或预定的)改变量。
通过过程60,该过程以每个表面位置处计算出的接触角的形式提供74表面润湿性更改的数值预测以及所得到的润湿性空间分布。在过程60中,离散的表面元和一个或多个流体元连通,从而使得局部润湿性确定能够取决于矿物类型的局部特性、流体特性、流体成分(例如、沥青质的存在)、毛细管压和分离压。缓慢的动态润湿性更改允许非润湿流体侵入非常小的孔隙中,这是因为润湿性可以以无需克服毛细管压的方式来被改变。
扩展地,有可能迭代地使用该方法,其中在执行润湿性更改之后进行流体相重新分布,重复过程60直到系统10收敛(进一步的重复没有显著改变结果(例如,接触角)而测量得出)为止。因此,过程60可以参考建立的阈值,该阈值提供了为了停止润湿性更改模拟的进一步迭代而需要改变接触角的量。
还可以在相对渗透率模拟中使用过程60来预测例如在物理相对渗透率实验室实验(这些实验可能花费几个月的时间,因此即使执行了润湿性恢复流程,润湿性的进一步更改也是可能的)期间发生的润湿性更改。
实现力平衡是限定润湿性更改的合格规则的基本机制。润湿性更改的合格规则涉及三个力:
(1)毛细管压(PC),即,非润湿(油)和润湿(水)流体的压力差。它可以直接来自边界条件,例如,油的入口压力–水的出口压力(P入口油-P出口水),但是其它测试方法也会造成得出毛细管压的油和水的分布,例如相对渗透率测试。
(2)基于曲率的压力(PV),即,由固体几何形状迫使油/水界面的形状变为该几何形状的形状(即,岩石表面/孔隙的形状等)从而引起的压力。曲率是使用几何算法进行计算的,几何算法通常可在各种软件库(诸如上面提到的VTK(可视化工具包))中获得。
(3)分离压(PD),即,防止膜破裂的压力,即,当施加到膜上的压力超过最大分离压时,膜将破裂。
如上文(3)所解释的,膜破裂的力平衡要求可以表示为:
(PP>PV+PD),或毛细管压大于基于曲率的压力加上分离压
或者
(PP-PV>PD),或毛细管压减去基于曲率的压力大于分离压。
润湿性指数模拟:
参考图4,示出了验证过程80。该验证过程通过使用82润湿性指数(诸如Amott指数和USBM方法(美国矿业局)指数,其是通过实验确定的We)形式的岩石样本的润湿性测量来提供对图2的润湿性计算的验证。这些值被用于验证和/或校准由图2的数字老化过程40产生的预测的润湿性分布。Amott测试是石油工程中储层岩心的众所周知的实验润湿性测量,它结合了两个自发吸收(spontaneous imbibition)测量和两个强制驱替(forceddisplacement)测量,并定义了两个不同的指数,即,从岩心分析实验中获得的Amott水指数和Amott油指数。
润湿性指数模拟过程80基于估计的输入参数(参见上文)执行84数字老化模拟以提供模拟结果。该过程使用该模拟结果确定86预测的润湿性指数的值Wp。润湿性指数模拟过程80将预测值与实验值We(来自82)进行比较88,并且当这些值与90不匹配时,Wp≠We±T充分在用户定义的公差T内。关于匹配的用户定义的该公差可以是基于归因于润湿性指数的实验确定的公差和/或用于确定预测的特性的公差等以经验地或数字方式确定的。
此后,过程80迭代修改92到数字老化模拟过程40的输入参数值中的一个或多个,并重复94数字老化模拟过程40直到来自数字老化模拟40的值与来自实验值的值相符从而指示验证98。
如以上针对图1所讨论的,参数值是由引擎33a获得的一个或多个颗粒表面特性33a’,该引擎将矿物类型指派给颗粒并确定这些矿物类型中的每一个的表面特性,诸如流体特性33b之类的参数,例如,每种预期流体的流体密度和粘度以及流体-流体界面张力特性、流体的化学成分数据和特定矿物表面的流体成分的亲和度数据以及每种矿物类型的分离压。通过系统10将参数值中的一个或多个增加既定的增量以完成迭代修改。递增的参数值和递增值量二者都可以由系统确定或由用户确定。
参考图5和图5A,这些现有技术图描绘了对油团(oil blob)和盐水的常规理解,其提供了用于示例性岩石样本部分的固定接触角惯例。(摘自A.R.Kovscek等人的DOE/BC—9200106发行DE92 001062,1992年9月,为美国能源部准备的“A Pore-Level Scenario forthe Development of Mixed-Wettability in Oil Reservoirs”。)
本说明书中描述的主题和功能操作的实施例可以在数字电子电路、有形体现的计算机软件或固件、计算机硬件(包括本说明书中公开的结构及其结构等同物)或它们中的一个或多个的组合中实现。本说明书中描述的主题的实施例可以实现为一个或多个计算机程序(即,在有形非暂时性程序载体上编码的一个或多个计算机程序指令模块,用于由数据处理装置执行或控制数据处理装置的操作)。计算机存储介质可以是机器可读存储设备、机器可读存储基板、随机或串行存取存储器设备、或它们中的一个或多个的组合。
计算机程序也可以被称为或描述为程序、软件、软件应用、模块、软件模块、脚本或代码,可以用任何形式的编程语言编写,包括编译或解释语言或者声明性或过程性语言,它可以以任何形式部署,包括作为独立程序或作为模块、组件、子例程或适用于计算环境的其他单元。计算机程序可以但不必对应于文件系统中的文件。程序可以存储在保存其他程序或数据的文件的一部分中(例如,存储在标记语言文档中的一个或多个脚本、专用于所讨论的程序的单个文件中、或者在多个协调文件(例如,存储一个或多个模块,子程序或代码部分的文件)中)。可以部署计算机程序,使得程序在一个计算机上或在位于一个站点上或分布在多个站点上并通过数据通信网络互连的多个计算机上执行。
适合于执行计算机程序的计算机可以基于通用或专用微处理器或两者,或任何其他类型的中央处理单元。一般而言,中央处理单元将从只读存储器或随机存取存储器或两者接收指令和数据。计算机的基本元件是用于执行指令的中央处理单元以及用于存储指令和数据的一个或多个存储器设备。一般而言,计算机还将包括用于存储数据的一个或多个大容量存储设备(例如,磁、磁光盘或光盘),或者操作耦合以从中接收数据或将数据传输到一个或多个大容量存储设备,或两者兼有,但是计算机不需要具有此类设备。
适用于存储计算机程序指令和数据的计算机可读介质包括介质和存储器设备上的所有形式的非易失性存储器,包括例如半导体存储器设备(例如,EPROM,EEPROM和闪存设备)、磁盘(例如,内部硬盘或可移动盘)、磁光盘、CD ROM和DVD-ROM盘。处理器和存储器可以由专用逻辑电路补充或并入专用逻辑电路中。
Claims (38)
1.一种计算机实现的方法,用于在包括水和油的至少两个流体存在的情况下对物理岩石样本的老化和润湿性更改进行模拟,所述方法包括:
由计算系统检索物理岩石样本的表示;
由所述计算系统模拟所述物理岩石样本的所述表示的孔隙空间中的表面位置的老化和润湿性更改;
由所述计算系统相对于先前的润湿性更改来计算所述物理岩石样本的所述表示的所述孔隙空间中的表面位置的局部曲率的改变,其中由所述计算系统计算局部曲率的所述改变包括:
由所述计算系统计算所述孔隙空间中的所述表面位置的局部毛细管压和局部分离压;
对于所述表面位置中的一些表面位置,通过将所述毛细管压与基于局部曲率的压力和局部分离压之和进行比较,确定是否将发生水-油膜破裂;以及
基于计算出的局部曲率的所述改变,由所述计算系统确定与所述物理岩石样本的所述表示对应的所述表面位置的经修改的润湿性。
2.如权利要求1所述的方法,还包括:
根据计算出的局部曲率,确定是否将发生水膜破裂;以及
对于所述表面位置中的具有水膜破裂的至少一些表面位置,基于局部表面特性、局部流体化学成分和老化时间来确定润湿性更改的程度。
3.如权利要求1所述的方法,其中,产生还包括:
将润湿性模拟所需的输入参数指派给计算机模拟环境。
4.如权利要求1所述的方法,其中,所述岩石样本的所述表示被建模为三角网格,并且流体由体元表示,体元被指派有属性。
5.如权利要求4所述的方法,其中,体元能够被切割以使所述流体与所述三角网格相符。
6.如权利要求1所述的方法,还包括:
根据用户指定的条件,执行排水模拟以建立流体相分布。
7.如权利要求6所述的方法,其中,所述用户指定的条件选自包括实验室条件和储层条件的代表的组。
8.如权利要求7所述的方法,其中,执行排水模拟还包括
使用流体计算技术来模拟排水,以针对用户指定的条件建立通过所述岩石样本表示的流体相分布。
9.如权利要求1所述的方法,其中,执行排水模拟还包括
在所述岩石样本表示中的每个表面处执行局部曲率计算;以及
应用所述局部曲率计算的结果来计算每个表面位置处的所述局部毛细管压。
10.如权利要求1所述的方法,其中,在执行确定润湿性更改的数值表示的迭代之后,该方法还包括:
允许流体相由于经修改的表面润湿性而重新分布;以及
重复对润湿性更改的数值表示的确定。
11.如权利要求1所述的方法,其中,该方法以每个表面位置处计算出的接触角的形式提供表面润湿性更改的数值预测以及所得到的润湿性空间分布。
12.如权利要求1所述的方法,其中,离散的表面元和一个或多个流体元连通,并进行取决于矿物类型的局部特性、流体特性、流体成分、毛细管压和分离压的局部润湿性确定;并且其中润湿性更改被迭代地执行,然后进行流体相重新分布,其中润湿性更改被重复进行直到所述计算系统检测到收敛,所述收敛是通过后续重复没有显著改变所确定的接触角的值而测量得出的。
13.如权利要求1所述的方法,其中,模拟包括:
模拟通过与气井或油井相邻的储层岩石的多相流动行为,并且其中确定多相流动行为包括用于确定所述物理岩石样本的润湿性更改的所述方法。
14.如权利要求1所述的方法,其中迭代地执行润湿性更改,然后进行流体相重新分布,重复进行润湿性更改,直到所述计算系统检测到通过后续重复测量得到的收敛,而不会明显改变确定的接触角的值。
15.一种计算机系统,包括:
一个或多个处理器设备;
存储器,耦合到所述一个或多个处理器设备;
存储装置,存储可执行的计算机指令,所述计算机指令用于在包括水和油的至少两个流体存在的情况下对物理岩石样本的老化和润湿性更改进行模拟,所述指令将所述一个或多个处理器配置为:
检索物理岩石样本的表示;
模拟所述物理岩石样本的所述表示的孔隙空间中的表面位置的老化和润湿性更改;
相对于先前的润湿性更改来计算所述物理岩石样本的所述表示的所述孔隙空间中的表面位置的局部曲率的改变,其中由所述计算机系统计算局部曲率的所述改变包括:
由所述计算机系统计算所述孔隙空间中的所述表面位置的局部毛细管压和局部分离压;
对于所述表面位置中的一些表面位置,通过将所述毛细管压与基于局部曲率的压力和局部分离压之和进行比较,
确定是否将发生水-油膜破裂;以及
基于计算出的局部曲率的所述改变,由所述计算机系统确定与所述物理岩石样本的所述表示对应的所述表面位置的经修改的润湿性。
16.如权利要求15所述的计算机系统,其中,所述计算机系统还被配置为:
根据计算出的局部曲率,确定是否将发生水膜破裂;以及
对于所述表面位置中的具有水膜破裂的至少一些表面位置,基于局部表面特性、局部流体化学成分和老化时间来确定润湿性更改的程度。
17.如权利要求15所述的计算机系统,其中,用于产生的指令还将所述一个或多个处理器配置为:
将润湿性模拟所需的输入参数指派给计算机模拟环境。
18.如权利要求15所述的计算机系统,其中,所述岩石样本的所述表示被建模为三角网格,并且流体由体元表示,体元被指派有属性。
19.如权利要求18所述的计算机系统,其中,体元能够被切割以使所述流体与所述三角网格相符。
20.如权利要求15所述的计算机系统,还将所述一个或多个处理器配置为:
根据用户指定的条件,执行排水模拟以建立流体相分布。
21.如权利要求20所述的计算机系统,其中,所述用户指定的条件选自包括实验室条件和储层条件的代表的组。
22.如权利要求21所述的计算机系统,其中,用于执行排水模拟的指令还将所述一个或多个处理器配置为:
使用流体计算技术来模拟排水,以针对用户指定的条件建立通过所述岩石样本表示的流体相分布。
23.如权利要求15所述的计算机系统,其中,用于执行排水模拟的指令还将所述一个或多个处理器配置为:
在所述岩石样本表示中的每个表面处执行局部曲率计算;以及
应用所述局部曲率计算的结果来计算每个表面位置处的所述局部毛细管压。
24.如权利要求15所述的计算机系统,其中,在执行确定润湿性更改的数值表示的迭代之后,所述指令还将所述一个或多个处理器配置为:
允许流体相由于经修改的表面润湿性而重新分布;以及
重复对润湿性更改的数值表示的确定。
25.如权利要求15所述的计算机系统,其中,所述计算机系统以每个表面位置处计算出的接触角的形式提供表面润湿性更改的数值预测以及所得到的润湿性空间分布。
26.如权利要求15所述的计算机系统,其中,离散的表面元和一个或多个流体元连通,并进行取决于矿物类型的局部特性、流体特性、流体成分、毛细管压和分离压的局部润湿性确定。
27.如权利要求15所述的计算机系统,其中,润湿性更改被迭代地执行,然后进行流体相重新分布,其中润湿性更改被重复进行直到所述计算机系统检测到收敛,所述收敛是通过后续重复没有显著改变所确定的接触角的值而测量得出的。
28.一种计算机可读非暂态存储设备,存储可执行的计算机指令,以在包括水和油的至少两个流体存在的情况下对物理岩石样本的老化和润湿性更改进行模拟,所述指令用于使计算系统:
检索物理岩石样本的表示;
模拟所述物理岩石样本的所述表示的孔隙空间中的表面位置的老化和润湿性更改;
相对于先前的润湿性更改来计算所述物理岩石样本的所述表示的所述孔隙空间中的表面位置的局部曲率的改变,其中由所述计算系统计算局部曲率的所述改变包括:
由所述计算系统计算所述孔隙空间中的所述表面位置的局部毛细管压和局部分离压;
对于所述表面位置中的一些表面位置,通过将所述毛细管压与基于局部曲率的压力和局部分离压之和进行比较,确定是否将发生水-油膜破裂;以及
基于计算出的局部曲率的所述改变,由所述计算系统确定与所述物理岩石样本的所述表示对应的所述表面位置的经修改的润湿性。
29.如权利要求28所述的计算机可读非暂态存储设备,其中,所述计算机可读非暂态存储设备还包括用于以下操作的指令:
根据计算出的局部曲率,确定是否将发生水膜破裂;以及
对于所述表面位置中的具有水膜破裂的至少一些表面位置,基于局部表面特性、局部流体化学成分和老化时间来确定润湿性更改的程度。
30.如权利要求28所述的计算机可读非暂态存储设备,其中,用于产生的指令还包括用于以下操作的指令:
将润湿性模拟所需的输入参数指派给计算机模拟环境。
31.如权利要求28所述的计算机可读非暂态存储设备,其中,所述岩石样本的所述表示被建模为三角网格,并且流体由体元表示,体元被指派有属性。
32.如权利要求31所述的计算机可读非暂态存储设备,其中,体元能够被切割以使所述流体与所述三角网格相符。
33.如权利要求28所述的计算机可读非暂态存储设备,还包括用于以下操作的指令:
接收用户指定的条件;以及
根据接收到的用户指定的条件,执行排水模拟以建立流体相分布。
34.如权利要求28所述的计算机可读非暂态存储设备,其中,用于执行排水模拟的指令还包括用于以下操作的指令:
在所述岩石样本表示中的每个表面处执行局部曲率计算;以及
应用所述局部曲率计算的结果来计算每个表面位置处的所述局部毛细管压。
35.如权利要求28所述的计算机可读非暂态存储设备,其中,在执行确定润湿性更改的数值表示的迭代之后,所述计算系统还被配置为:
允许流体相由于经修改的表面润湿性而重新分布;以及
重复对润湿性更改的数值表示的确定。
36.如权利要求28所述的计算机可读非暂态存储设备,其中,所述指令使所述计算系统以每个表面位置处计算出的接触角的形式提供表面润湿性更改的数值预测以及所得到的润湿性空间分布。
37.如权利要求28所述的计算机可读非暂态存储设备,其中,离散的表面元和一个或多个流体元连通,并进行取决于矿物类型的局部特性、流体特性、流体成分、毛细管压和分离压的局部润湿性确定。
38.如权利要求28所述的计算机可读非暂态存储设备,其中,润湿性更改被迭代地执行,然后进行流体相重新分布,其中润湿性更改被重复进行直到所述计算系统检测到收敛,所述收敛是通过后续重复没有显著改变所确定的接触角的值而测量得出的。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201862615701P | 2018-01-10 | 2018-01-10 | |
US62/615,701 | 2018-01-10 | ||
PCT/US2019/013079 WO2019140108A1 (en) | 2018-01-10 | 2019-01-10 | Determining fluid flow characteristics of porous mediums |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111936841A CN111936841A (zh) | 2020-11-13 |
CN111936841B true CN111936841B (zh) | 2023-11-17 |
Family
ID=67139489
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201980012848.7A Active CN111936841B (zh) | 2018-01-10 | 2019-01-10 | 确定多孔介质的流体流特点 |
Country Status (5)
Country | Link |
---|---|
US (1) | US11714040B2 (zh) |
EP (1) | EP3737928A4 (zh) |
JP (1) | JP7360390B2 (zh) |
CN (1) | CN111936841B (zh) |
WO (1) | WO2019140108A1 (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110573842B (zh) | 2017-01-26 | 2021-05-11 | 达索系统西姆利亚公司 | 基于流体占据时间的多相流可视化图 |
US11530598B2 (en) | 2018-08-21 | 2022-12-20 | Dassault Systemes Simulia Corp. | Determination of oil removed by gas via miscible displacement in reservoir rock |
US11613984B2 (en) | 2019-09-04 | 2023-03-28 | Dassault Systemes Simulia Corp. | Determination of hydrocarbon mobilization potential for enhanced oil recovery |
US11847391B2 (en) | 2020-06-29 | 2023-12-19 | Dassault Systemes Simulia Corp. | Computer system for simulating physical processes using surface algorithm |
US11907625B2 (en) | 2020-12-29 | 2024-02-20 | Dassault Systemes Americas Corp. | Computer simulation of multi-phase and multi-component fluid flows including physics of under-resolved porous structures |
US12002159B2 (en) | 2021-06-09 | 2024-06-04 | Dassault Systemes Americas Corp. | Prediction of mechanical properties of sedimentary rocks based on a grain to grain parametric cohesive contact model |
CN114088585B (zh) * | 2021-11-15 | 2022-08-12 | 中国石油大学(北京) | 一种页岩润湿角预测模型建立的方法、装置、设备及介质 |
US12044121B2 (en) | 2021-11-23 | 2024-07-23 | Saudi Arabian Oil Company | In-situ wettability measurement using multiple partitioning tracers |
CN115406806B (zh) * | 2022-11-02 | 2023-04-11 | 山东省中地易采石油技术有限责任公司 | 一种致密油气藏两相渗流实验的计量装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102822446A (zh) * | 2010-02-12 | 2012-12-12 | 英国石油勘探运作有限公司 | 用于预测注入到含油储层的微生物的影响的方法和系统 |
CN104849765A (zh) * | 2009-12-16 | 2015-08-19 | 英国石油勘探运作有限公司 | 用于测量岩石润湿性的方法 |
CN105393110A (zh) * | 2013-08-06 | 2016-03-09 | Bp北美公司 | 在模拟应力和应变条件下的岩石物理属性的基于图像的直接数值模拟 |
WO2016159807A1 (ru) * | 2015-03-27 | 2016-10-06 | Шлюмберже Текнолоджи Корпорейшн | Способ определения равновесной смачиваемости образца горной породы |
Family Cites Families (102)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3818757A (en) | 1972-05-05 | 1974-06-25 | Saratoga Sys Inc | Dual path ultrasonic fluid flow metering system and method |
GB2066466A (en) | 1979-12-03 | 1981-07-08 | Water Res Centre | Ultrasonic Flow Measurement |
US4498134A (en) | 1982-01-26 | 1985-02-05 | Hughes Aircraft Company | Segregator functional plane for use in a modular array processor |
US4507726A (en) | 1982-01-26 | 1985-03-26 | Hughes Aircraft Company | Array processor architecture utilizing modular elemental processors |
US4729098A (en) | 1985-06-05 | 1988-03-01 | General Electric Company | System and method employing nonlinear interpolation for the display of surface structures contained within the interior region of a solid body |
US4831519A (en) | 1985-12-12 | 1989-05-16 | Itt Corporation | Cellular array processor with variable nesting depth vector control by selective enabling of left and right neighboring processor cells |
US4809202A (en) | 1985-12-27 | 1989-02-28 | Thinking Machines Corporation | Method and apparatus for simulating systems described by partial differential equations |
JPH0648486B2 (ja) | 1986-10-08 | 1994-06-22 | 日本電気株式会社 | ベクトルデ−タ処理装置 |
US4933895A (en) | 1987-07-10 | 1990-06-12 | Hughes Aircraft Company | Cellular array having data dependent processing capabilities |
JP2635617B2 (ja) | 1987-09-29 | 1997-07-30 | 株式会社東芝 | 半導体素子特性評価用の直交格子点の発生方法 |
US4989166A (en) | 1988-07-12 | 1991-01-29 | Hitachi, Ltd. | Method for synthesizing analysis model and flow analysis system |
US5038302A (en) | 1988-07-26 | 1991-08-06 | The Research Foundation Of State University Of New York | Method of converting continuous three-dimensional geometrical representations into discrete three-dimensional voxel-based representations within a three-dimensional voxel-based system |
US5159690A (en) | 1988-09-30 | 1992-10-27 | Massachusetts Institute Of Technology | Multidimensional cellular data array processing system which separately permutes stored data elements and applies transformation rules to permuted elements |
JPH03276071A (ja) | 1990-03-27 | 1991-12-06 | Yoshiomi Kondo | 流体及び電磁流体の物理量予測方法 |
DE69122754T3 (de) | 1990-05-09 | 2000-10-19 | Memtec Ltd., South Windsor | Verfahren zur herstellung von pvdf hohlfasermembranen |
US5377129A (en) | 1990-07-12 | 1994-12-27 | Massachusetts Institute Of Technology | Particle interaction processing system |
JPH05508498A (ja) | 1990-07-12 | 1993-11-25 | マサチユセツツ・インスチチユート・オブ・テクノロジー | 流体力学並列プロセッサー |
IE69192B1 (en) | 1990-12-21 | 1996-08-21 | Hitachi Europ Ltd | A method of generating partial differential equations for simulation a simulation method and a method of generating simulation programs |
US5069065A (en) * | 1991-01-16 | 1991-12-03 | Mobil Oil Corporation | Method for measuring wettability of porous rock |
US5442733A (en) | 1992-03-20 | 1995-08-15 | The Research Foundation Of State University Of New York | Method and apparatus for generating realistic images using a discrete representation |
JP3426647B2 (ja) | 1992-06-24 | 2003-07-14 | 日本電信電話株式会社 | 3次元トポグラフィシミュレーションのための一般化されたソリッドモデリング |
US5361385A (en) | 1992-08-26 | 1994-11-01 | Reuven Bakalash | Parallel computing system for volumetric modeling, data processing and visualization |
US5424963A (en) | 1992-11-25 | 1995-06-13 | Photon Research Associates, Inc. | Molecular dynamics simulation method and apparatus |
US5848260A (en) | 1993-12-10 | 1998-12-08 | Exa Corporation | Computer system for simulating physical processes |
US5606517A (en) | 1994-06-08 | 1997-02-25 | Exa Corporation | Viscosity reduction in physical process simulation |
US5548694A (en) | 1995-01-31 | 1996-08-20 | Mitsubishi Electric Information Technology Center America, Inc. | Collision avoidance system for voxel-based object representation |
US5640335A (en) | 1995-03-23 | 1997-06-17 | Exa Corporation | Collision operators in physical process simulation |
US5605517A (en) | 1995-11-21 | 1997-02-25 | General Motors Corporation | Multi-speed power transmission |
US5953239A (en) | 1997-12-29 | 1999-09-14 | Exa Corporation | Computer simulation of physical processes |
US7209873B1 (en) | 1999-09-03 | 2007-04-24 | Kliegel James R | Method of analyzing and computing anisotropic turbulent flows |
US6318156B1 (en) | 1999-10-28 | 2001-11-20 | Micro Motion, Inc. | Multiphase flow measurement system |
US6928399B1 (en) | 1999-12-03 | 2005-08-09 | Exxonmobil Upstream Research Company | Method and program for simulating a physical system using object-oriented programming |
US6516080B1 (en) | 2000-04-05 | 2003-02-04 | The Board Of Trustees Of The Leland Stanford Junior University | Numerical method of estimating physical properties of three-dimensional porous media |
FR2811760B1 (fr) | 2000-07-17 | 2002-09-13 | Inst Francais Du Petrole | Methode pour modeliser des deplacements de fluides dans un milieu poreux tenant compte d'effets d'hysteresis |
US6915245B1 (en) | 2000-09-14 | 2005-07-05 | General Atomics | Method of simulating a fluid by advection of a time-weighted equilibrium distribution function |
FR2844355B1 (fr) * | 2002-09-11 | 2005-03-25 | Inst Francais Du Petrole | Methode de mesure de la mouillabilite de roches par resonnance magnetique nucleaire |
US7379852B2 (en) | 2004-02-18 | 2008-05-27 | Chevron U.S.A. Inc. | N-phase interface tracking method utilizing unique enumeration of microgrid cells |
EP1729151B1 (en) | 2005-05-27 | 2011-06-29 | Services Petroliers Schlumberger | Method and apparatus for measuring the wettability of geological 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 |
US20070016389A1 (en) | 2005-06-24 | 2007-01-18 | Cetin Ozgen | Method and system for accelerating and improving the history matching of a reservoir simulation model |
US7590514B1 (en) | 2006-05-24 | 2009-09-15 | Livermore Software Technology Corporation | Method and system for simulating inflation of an enclosed volume with flexible boundary using gas particles |
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 |
US20080012853A1 (en) | 2006-07-14 | 2008-01-17 | William Geiger | Generating mesh from implicit surface |
US7558714B2 (en) | 2006-08-10 | 2009-07-07 | Exa Corporation | Computer simulation of physical processes |
CA2660877A1 (en) | 2006-08-17 | 2008-02-21 | University Of Massachusetts | Stabilized emulsions, methods of preparation, and related reduced fat foods |
US7983886B2 (en) | 2007-09-10 | 2011-07-19 | Chevron U.S.A. Inc. | Methods for performing simulation of surfactant flooding of a hydrocarbon reservoir |
EP2289046A4 (en) | 2008-05-23 | 2016-03-30 | Fei Co | IMAGE DATA PROCESSING |
BRPI0911949A2 (pt) | 2008-06-25 | 2019-09-24 | Prad Research And Development Limited | método para determinar a permeabilidade de uma formação subterrânea e programa de computador armazenado numa mídia legível por computador |
GB2462261A (en) | 2008-07-28 | 2010-02-03 | Fujitsu Ltd | Method, apparatus and computer program for simulating behaviou r of thermodynamic systems |
JP5715568B2 (ja) | 2008-10-10 | 2015-05-07 | ヴェロシス,インク. | マイクロチャネルプロセス技術を使用するプロセスおよび装置 |
US8155377B2 (en) | 2008-11-24 | 2012-04-10 | Ingrain, Inc. | Method for determining rock physics relationships using computer tomographic images thereof |
KR101199475B1 (ko) * | 2008-12-22 | 2012-11-09 | 한국전자통신연구원 | 3차원 모델 생성 방법 및 장치 |
CN101556687A (zh) | 2009-04-23 | 2009-10-14 | 上海交通大学 | 微细通道内两相流的三维可视化系统及方法 |
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 |
CN102472649A (zh) | 2009-07-07 | 2012-05-23 | 西门子公司 | 用于测量多相流体流的设备和方法 |
JP2011191848A (ja) | 2010-03-12 | 2011-09-29 | Toyota Central R&D Labs Inc | 拡散現象解析方法、拡散現象解析装置およびプログラム |
US10083254B2 (en) | 2010-06-15 | 2018-09-25 | Exxonmobil Upstream Research Company | Method and system for stabilizing formulation methods |
US8412501B2 (en) | 2010-06-16 | 2013-04-02 | Foroil | Production simulator for simulating a mature hydrocarbon field |
US8757871B2 (en) | 2010-08-16 | 2014-06-24 | The Board Of Trustees Of The University Of Illinois | Particle dynamics microscopy using temperature jump and probe anticorrelation/correlation techniques |
CN102383783B (zh) | 2010-09-03 | 2015-06-17 | 中国石油化工股份有限公司 | 一种分析缝洞型油藏孔洞间油水流动特征的方法 |
CA2807300C (en) | 2010-09-20 | 2017-01-03 | Exxonmobil Upstream Research Company | Flexible and adaptive formulations for complex reservoir simulations |
US8768628B2 (en) * | 2010-10-20 | 2014-07-01 | Shawket Ghedan | Rise in core wettability characterization method |
CA2814669A1 (en) | 2010-11-23 | 2012-05-31 | Exxonmobil Upstream Research Company | Variable discretization method for flow simulation on complex geological models |
US8645070B2 (en) | 2010-11-24 | 2014-02-04 | Chevron U.S.A. Inc. | System and method for estimating fluid distribution in a subterranean reservoir |
US8583411B2 (en) | 2011-01-10 | 2013-11-12 | Saudi Arabian Oil Company | Scalable simulation of multiphase flow in a fractured subterranean reservoir as multiple interacting continua |
US9051825B2 (en) | 2011-01-26 | 2015-06-09 | Schlumberger Technology Corporation | Visualizing fluid flow in subsurface reservoirs |
CN102174888B (zh) * | 2011-03-03 | 2012-05-23 | 康志勇 | 一种基于储集岩岩性参数的地层数据处理方法 |
US9507047B1 (en) | 2011-05-10 | 2016-11-29 | Ingrain, Inc. | Method and system for integrating logging tool data and digital rock physics to estimate rock formation properties |
CA2840942C (en) | 2011-07-12 | 2017-05-02 | Ingrain, Inc. | Method for simulating fractional multi-phase/multi-component flow through porous media |
US9037440B2 (en) | 2011-11-09 | 2015-05-19 | Exa Corporation | Computer simulation of fluid flow and acoustic behavior |
US10360324B2 (en) | 2011-12-09 | 2019-07-23 | Dassault Systemes Simulia Corp. | Computer simulation of physical processes |
EP2795528A4 (en) | 2012-02-10 | 2016-06-29 | Landmark Graphics Corp | SYSTEMS AND METHOD FOR EVALUATING THE LIQUID BREAKDOWN TIMES AT PRODUCING BORING HOURS |
US20130317791A1 (en) | 2012-04-26 | 2013-11-28 | Conocophillips Company | Hydrodynamic slug flow model |
BR112014028302A2 (pt) | 2012-05-14 | 2019-09-24 | Landmark Graphics Corp | método e sistema de exibir uma representação gráfica dos parâmetros de produção de hidrocarbonetos |
RU2596692C2 (ru) | 2012-05-14 | 2016-09-10 | Лэндмарк Графикс Корпорейшн | Способ и система для оценки прогнозной добычи углеводородов |
US9140117B2 (en) | 2012-07-13 | 2015-09-22 | Ingrain, Inc. | Method for evaluating relative permeability for fractional multi-phase, multi-component fluid flow through porous media |
US9619592B2 (en) | 2012-08-31 | 2017-04-11 | Schlumberger Technology Corporation | Analysis of enhanced oil recovery (EOR) processes for naturally-fractured reservoirs |
US10550690B2 (en) * | 2013-05-16 | 2020-02-04 | Dassault Systemes Simulia Corp. | Mass exchange model for relative permeability simulation |
WO2015030837A1 (en) | 2013-08-27 | 2015-03-05 | Halliburton Energy Services, Inc. | Simulating fluid leak-off and flow-back in a fractured subterranean |
US9913616B2 (en) | 2013-07-30 | 2018-03-13 | Heartflow, Inc. | Method and system for modeling blood flow with boundary conditions for optimized diagnostic performance |
WO2015017648A2 (en) | 2013-07-31 | 2015-02-05 | Exa Corporation | Temperature coupling algorithm for hybrid thermal lattice boltzmann method |
AR097342A1 (es) | 2013-08-15 | 2016-03-09 | Univ Texas | Sistema alcalino tolerante a cationes minerales multivalentes para eor química |
US9377392B2 (en) | 2013-09-05 | 2016-06-28 | Proptester, Inc. | Methods and systems for testing fluids on crushed formation materials under conditions of stress |
US9910173B2 (en) | 2013-11-15 | 2018-03-06 | Schlumberger Technology Corporation | Saturation end-point adjustment |
US20150141303A1 (en) | 2013-11-18 | 2015-05-21 | The Board Of Regents Of The University Of Oklahoma | Binary and ternary surfactant blends for enhanced oil recovery in reservoir brines with extremely high total dissolved solids |
US9797225B2 (en) | 2013-11-27 | 2017-10-24 | Saudi Arabian Oil Company | Data compression of hydrocarbon reservoir simulation grids |
CN103698102B (zh) | 2013-12-17 | 2016-08-17 | 哈尔滨工程大学 | 一种管束间高温高压汽液两相流型实验装置及流型判别方法 |
WO2015142931A1 (en) | 2014-03-18 | 2015-09-24 | Crouse Bernd | Fluid blob tracking for evaluation of multiphase flow simulations |
FR3023316B1 (fr) | 2014-07-04 | 2016-08-19 | Ifp Energies Now | Procede d'exploitation d'un gisement petrolier a partir d'une technique de positionnement des puits a forer |
US10169498B2 (en) | 2015-04-14 | 2019-01-01 | California Institute Of Technology | Particle-based simulation through power diagrams |
GB2554013B (en) | 2015-05-01 | 2021-03-24 | Geoquest Systems Bv | Multiphase flow in porous media |
US10436003B2 (en) * | 2015-12-17 | 2019-10-08 | Baker Hughes, A Ge Company, Llc | Fluid blocking analysis and chemical evalution |
CN105910996A (zh) | 2016-04-07 | 2016-08-31 | 昆明理工大学 | 一种利用rgb颜色模型确定多相流混合时间的方法 |
WO2017222509A1 (en) | 2016-06-22 | 2017-12-28 | Schlumberger Technology Corporation | Visualizations of reservoir simulations with fracture networks |
CN110573842B (zh) | 2017-01-26 | 2021-05-11 | 达索系统西姆利亚公司 | 基于流体占据时间的多相流可视化图 |
US10691846B2 (en) * | 2017-03-01 | 2020-06-23 | International Business Machines Corporation | Capillary network simulations based on a low-dimensional representation of porous media |
EP3710555A1 (en) | 2017-11-14 | 2020-09-23 | Total S.A. | Surfactant composition for the reinjection of produced water |
US10830713B2 (en) | 2017-11-20 | 2020-11-10 | DigiM Solution LLC | System and methods for computing physical properties of materials using imaging data |
US20190203593A1 (en) * | 2017-12-29 | 2019-07-04 | Shawn Fullmer | Method and System for Modeling in a Subsurface Region |
US11530598B2 (en) | 2018-08-21 | 2022-12-20 | Dassault Systemes Simulia Corp. | Determination of oil removed by gas via miscible displacement in reservoir rock |
US11613984B2 (en) | 2019-09-04 | 2023-03-28 | Dassault Systemes Simulia Corp. | Determination of hydrocarbon mobilization potential for enhanced oil recovery |
US11847391B2 (en) | 2020-06-29 | 2023-12-19 | Dassault Systemes Simulia Corp. | Computer system for simulating physical processes using surface algorithm |
-
2019
- 2019-01-09 US US16/243,285 patent/US11714040B2/en active Active
- 2019-01-10 WO PCT/US2019/013079 patent/WO2019140108A1/en unknown
- 2019-01-10 CN CN201980012848.7A patent/CN111936841B/zh active Active
- 2019-01-10 JP JP2020538656A patent/JP7360390B2/ja active Active
- 2019-01-10 EP EP19738994.3A patent/EP3737928A4/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104849765A (zh) * | 2009-12-16 | 2015-08-19 | 英国石油勘探运作有限公司 | 用于测量岩石润湿性的方法 |
CN102822446A (zh) * | 2010-02-12 | 2012-12-12 | 英国石油勘探运作有限公司 | 用于预测注入到含油储层的微生物的影响的方法和系统 |
CN105393110A (zh) * | 2013-08-06 | 2016-03-09 | Bp北美公司 | 在模拟应力和应变条件下的岩石物理属性的基于图像的直接数值模拟 |
WO2016159807A1 (ru) * | 2015-03-27 | 2016-10-06 | Шлюмберже Текнолоджи Корпорейшн | Способ определения равновесной смачиваемости образца горной породы |
Also Published As
Publication number | Publication date |
---|---|
US11714040B2 (en) | 2023-08-01 |
JP7360390B2 (ja) | 2023-10-12 |
EP3737928A1 (en) | 2020-11-18 |
WO2019140108A1 (en) | 2019-07-18 |
EP3737928A4 (en) | 2021-03-03 |
CN111936841A (zh) | 2020-11-13 |
JP2021510419A (ja) | 2021-04-22 |
US20190212241A1 (en) | 2019-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111936841B (zh) | 确定多孔介质的流体流特点 | |
Jing et al. | Coal cleat reconstruction using micro-computed tomography imaging | |
US11163923B2 (en) | Automated upscaling of relative permeability and capillary pressure in multi-porosity systems | |
Alpak et al. | Prediction of fluid topology and relative permeability in imbibition in sandstone rock by direct numerical simulation | |
CN109891460B (zh) | 具有两个或更多个不混溶相的岩石的岩石物理性质的直接数值模拟 | |
EP3030888A1 (en) | Image-based direct numerical simulation of petrophysical properties under simulated stress and strain conditions | |
EP3077618B1 (en) | Tuning digital core analysis to laboratory results | |
Ayoub et al. | A new correlation for accurate prediction of oil formation volume factor at the bubble point pressure using Group Method of Data Handling approach | |
US20190203593A1 (en) | Method and System for Modeling in a Subsurface Region | |
JP2022104831A (ja) | 低解像度の多孔質構造の物理を含む多相多成分流体流のコンピュータシミュレーション | |
Amabeoku et al. | Evaluation and application of digital rock physics (DRP) for special core analysis in carbonate formations | |
Fayazi et al. | Estimation of pseudo relative permeability curves for a heterogeneous reservoir with a new automatic history matching algorithm | |
Ebeltoft et al. | Parameter based scal-analysing relative permeability for full field application | |
Cui et al. | Apparent permeability and representative size of shale: a numerical study on the effects of organic matter | |
Taylor et al. | Geometric and hydraulic void constrictions in granular media | |
WO2004083827A1 (en) | Method and system for predicting the apparent flow conductivity of a heterogeneous medium | |
French | Multiscale modeling of particle transport in petroleum reservoirs | |
EP4113117A1 (en) | System and method for correlating oil distribution during drainage and imbibition using machine learning | |
Foroughi et al. | Investigating the permeability–porosity relation of percolation-based porous media using the lattice Boltzmann method | |
Ringrose et al. | Upscaling flow properties | |
Blears | Coupled Pore-To-Continuum Multiscale Modeling of Dynamic Particle Filtration Processes in Porous Media | |
Assis et al. | Sensitivity analysis of hydraulic properties of naturally fractured reservoirs | |
Victor | Pore scale modeling of rock transport properties | |
Villegas et al. | Shape reconstruction from two-phase incompressible flow data using level sets | |
de Chalendar | Impact of Ostwald ripening on residually trapped carbon dioxide |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240409 Address after: Massachusetts USA Patentee after: Dassault Systemes Americas Corp. Country or region after: U.S.A. Address before: Rhode Island USA Patentee before: Dassault Systemes Simulia Corp. Country or region before: U.S.A. |