CN114021498B - 一种预测多相孔隙介质弹性模量的高效数值模拟方法 - Google Patents
一种预测多相孔隙介质弹性模量的高效数值模拟方法 Download PDFInfo
- Publication number
- CN114021498B CN114021498B CN202111303690.4A CN202111303690A CN114021498B CN 114021498 B CN114021498 B CN 114021498B CN 202111303690 A CN202111303690 A CN 202111303690A CN 114021498 B CN114021498 B CN 114021498B
- Authority
- CN
- China
- Prior art keywords
- fluid
- equation
- boundary
- numerical simulation
- modulus
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000004088 simulation Methods 0.000 title claims abstract description 23
- 239000011148 porous material Substances 0.000 title claims abstract description 15
- 239000012530 fluid Substances 0.000 claims abstract description 82
- 239000011435 rock Substances 0.000 claims abstract description 32
- 239000007787 solid Substances 0.000 claims abstract description 17
- 238000006073 displacement reaction Methods 0.000 claims abstract description 15
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 7
- 230000001133 acceleration Effects 0.000 claims description 8
- 230000000737 periodic effect Effects 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 230000010355 oscillation Effects 0.000 claims description 6
- 238000002591 computed tomography Methods 0.000 claims description 5
- 230000035699 permeability Effects 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 22
- 239000006185 dispersion Substances 0.000 abstract description 10
- 230000008878 coupling Effects 0.000 abstract description 2
- 238000010168 coupling process Methods 0.000 abstract description 2
- 238000005859 coupling reaction Methods 0.000 abstract description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 11
- 230000008569 process Effects 0.000 description 8
- 235000019738 Limestone Nutrition 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 3
- 239000003245 coal Substances 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 238000005315 distribution function Methods 0.000 description 3
- 238000011545 laboratory measurement Methods 0.000 description 3
- 239000006028 limestone Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003534 oscillatory effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种预测多相孔隙介质弹性模量的高效数值模拟方法,包括:通过流体分布图像获得流体分布特征;定义流体密度、粘度、体积模量分布特征;构建频率域Biot方程的位移‑压力形式;约束纵波边界条件;全局粗化网格剖分,局部网格精细剖分求解;使用有限元求解方程,获得应变固体位移场和流体压力场;计算纵波模量。本发明提供的预测多相孔隙介质弹性模量的高效数值模拟方法,相比于以前求解孔弹性Biot流‑固耦合位移方程,实现了对部分饱和岩石样品频散和衰减预测的有限元数值模拟方法,特别适用于模型中存在小非均质目标,可避免常规Biot固体位移‑流体位移方程计算效率低,四面体网格剖分内存占用量大,计算时间长的问题。
Description
技术领域
本发明涉及油气及煤层气勘探开发技术领域,具体涉及一种预测多相孔隙介质弹性模量的高效数值模拟方法。
背景技术
在油气或煤层气勘探开发过程中,评估储层岩石中多相流体对岩石声学属性的影响是地震资料成像、解释及油气或者煤层气资源储量评估的重要约束条件。获得岩石中流体分布与岩石声学属性关系的方法主要包括实验室直接测量法和数值模拟预测法:实验室测量法优点在于承认了岩石的各向异性及非均匀性特征,测量结果直接、准确,是对数值模拟方法的重要约束和验证。然而,较高的实验成本使实验室直接测量方法无法在研究者中广泛推广;数值模拟法相对于高成本的实验法而言,具有广义性,计算成本低等特点。因此,在过去数十年,众多研究者针对其展开了数量庞大的研究。
White等首先提出了流体分布会造成岩石介观尺度的非均质性,从而导致岩石的声波速度出现频散和衰减。除此之外,White等认为气-水交互的层状模型和均匀球状气泡分布模型等价于大尺度非均质体对地震波属性产生的影响。Dutta&Ode等将White的球状模型从二维情况拓展为三维情况。对这类模型的求解表明:在物理上,低频Biot纵波通过孔弹性介质的表现特征等价于通过更大尺度的粘弹性固体的表现特征。据此,一些研究者通过特定模型对孔弹性介质的地震响应进行了解析求解:Steven R.Pride和Berryman的双孔模型;Norris的层状模型;Johnson的任意几何形状的斑块饱和模型;Müller的随机多孔介质模型;Toms的随机离散和连续流体分布模型。
上述解析模型能很好的处理满足常规几何形状的非均质孔弹性模型。但是,对于较为复杂的非均质情况,人们通过数值模拟波传播过程进行研究。虽然数值模拟方法计算精度高,可以预测更加复杂的情况,但是,其在计算上昂贵,甚至有时无法完成计算。这是因为:需要非常精细的网格来表示非均质性;与地震波长相比,扩散长度非常小,因此在低频范围内,与流体压力平衡相关的扩散过程的分辨率是一个关键问题。为了克服限制,Masson&Pride提出了一种不同以往且有趣的数值计算方法。他们将随着时间变化的应力施加到非均质样品的边界上,并通过数值计算平均应力和应变场,确定其有效复模量。数值结果说明模拟过程的稳定性与岩石中流体粘度密切相关。然而,在时间域求解Biot方程时无法避免小的时间积分步长问题,同时需要将方程按照硬孔和非硬孔分开,这使得整个计算过程处于冗余状态。此外,对于频变模型,除了动态渗透率,方程中所有弹性模量均与频率相关,这将会在时间域求解时造成额外的卷积计算。相反,在频率域求解Biot方程则不需要考虑这类问题。在频率域可以避免卷积计算,使整个方程的求解过程变方便。Rubino等人还提出了一种升尺度数值模拟方法,以获得多相流体饱和多孔弹性介质的频散和衰减特征。该方法基于空间频率域中经典Biot方程的有限元解模拟振荡压缩性和剪切测试。虽然该数值方法大幅度提高了计算效率,但是在有限元计算过程中需要求取存在二阶导数的流体和固体位移方程,不利于有限元方法计算;除此之外,对于复杂的流体分布,尤其岩石模型中存在细小气泡的三维模型,常规网格剖分导致资源浪费巨大,最终常常无法计算,获得有效结果。
发明内容
为了克服现有技术方案的不足,本发明提供一种预测多相孔隙介质弹性模量的高效数值模拟方法,利用Biot位移-压力方程进行有限元算法高效求解非均质模型频散衰减特征。为了实现上述技术目的,本发明采用如下技术方案:
一种预测多相孔隙介质弹性模量的高效数值模拟方法,包括如下步骤:
S1:通过流体分布图像获得第i相流体分布范围Ωi,i=1,2…n;n为混合流体的最大相数;
S2:定义流体密度ρf为:ρf(x,y,z)=ρi,x,y,z∈Ωi,i=1,2…n,x,y,z为笛卡尔坐标系空间坐标;ρi是第i项流体密度;
S3:定义流体粘度为ηf:ηf=ηi,x,y,z∈Ωi,i=1,2…n;ηi是第i相流体粘度;
S4:定义流体体积模量Kf为:Kf=ki,x,y,z∈Ωi,i=1,2…n;
S5:构建频率域Biot方程的位移-压力形式:
其中,S=C:ε-αBpI;ρav为饱和岩石密度;ρav=ρd+φρf;ρd为岩石骨架密度,φ为孔隙度,ρc为流体复密度,τ为扭曲因子,κp为岩石渗透率;S为总应力,C为干燥岩石骨架柔度,ε为应变张量,αB为Biot-Willis系数,p为流体压力,I为单位张量矩阵,F为周期震荡的外力;ω为角频率,u为固体位移张量;
S6:约束边界条件:数值模拟的目标是实验室岩石样品,岩石样品轴向方向是z方向,垂直于z方向的平面上存在着正交的x、y方向,根据纵波模量的定义,横向应变εxx=εxy=εyy=εyx=0,及流体边界没有流体流出,因此边界条件为:
下边界为固定边界,需要同时满足两个条件:边界法向上流体加速度为0,固体位移是0,因此边界条件是:
上边界加载周期震荡应力,对应的边界条件是边界法向应力是σzz,流体流动的加速度是0,因此边界条件是:
S7:全局粗化网格剖分,局部网格精细剖分求解;
S8:使用有限元求解方程,获得应变固体位移场和流体压力场;
S9:计算纵波模量M:
其中,σzz是垂向应力,εzz是垂向应变;
则纵波速度Vp和纵波衰减Qp为:
优选地,步骤S1中流体分布图像通过CT扫描获得。
优选地,步骤S7具体包括:
S71:利用有限元算法使用粗化网格求解方程;
S72:计算所有网格元素上的偏微分方程的残差;
S73:估计所有网格元素的解中误差,其中第l方程的L2范数误差为h为局部网格尺寸,ρl为第1个偏微分方程的残差,qi为偏微分方程稳定性估计导数阶数,默认值是2,由于参数C是任意常数,计算出的误差实际上是误差指示因子;
S74:如果已进行了预设的优化次数或超过了最大有限元剖分数,则终止执行;
S76:重复步骤S71-S75。
与现有技术相比,本发明的有益效果:
本发明提供的预测多相孔隙介质弹性模量的高效数值模拟方法,相比于以前求解孔弹性Biot方程,实现对部分饱和岩石样品频散和衰减预测的有限元数值模拟方法,具有同等计算精度下,占据内存少,计算时间少的特点,特别适用于模型中存在小非均质目标,可避免常规四面体网格剖分占用大量内存,计算占用大量时间的问题。
附图说明
为了更清楚的说明本发明实施例或现有技术的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做简单的介绍,显而易见的,下面描述中的附图仅仅是本发明的一些实施例,对于本领域中的普通技术人员来说,在不付出创造性劳动的前提下,还可根据这些附图获得其他附图。
图1为气水两相流的双层模型;
图2为本发明实施例频散和衰减图;
图3为本发明实施例水/气的流体分布模式;
图4为本发明实施例流体体积模量、密度和粘度;
图5为本发明实施例横向边界条件加载位置;
图6为本发明实施例震源加载位置;
图7为本发明实施例固定边界条件加载位置;
图8为本发明实施例网格初步剖分;
图9为本发明实施例网格自适应剖分;
图10为本发明实施例在频率1000Hz时的流体压力分布图;
图11为本发明实施例体积模量随频率变化图;
图12为本发明实施例体积衰减随频率变化图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种预测多相孔隙介质弹性模量的高效数值模拟方法,包括如下步骤:
S1:通过流体分布图像获得第i相流体分布范围Ωi,流体分布图像通过CT扫描获得,i=1,2…n;n为混合流体的最大相数;
S2:定义流体密度ρf为:ρf(x,y,z)=ρi,x,y,z∈Ωi,i=1,2…n,x,y,z为笛卡尔坐标系空间坐标;ρi是第i项流体密度;
S3:定义流体粘度为ηf:ηf=ηi,x,y,z∈Ωi,i=1,2…n;ηi是第i相流体粘度;
S4:定义流体体积模量Kf为:Kf=ki,x,y,z∈Ωi,i=1,2…n;
S5:构建频率域Biot方程的位移-压力形式:
其中,ρav为饱和岩石密度;ρav=ρd+φρf;ρd为岩石骨架密度,φ为孔隙度,ρc为流体复密度,τ为扭曲因子,κp为岩石渗透率;S为总应力,C为干燥岩石骨架柔度,ε为应变张量,αB为Biot-Willis系数,p为流体压力,I为单位张量矩阵,F为周期震荡的外力;ω为角频率,u为固体位移张量;
S6:约束边界条件:数值模拟的目标是实验室岩石样品,岩石样品轴向方向是z方向,垂直于z方向的平面上存在着正交的x、y方向,根据纵波模量的定义,横向应变εxx=εxy=εyy=εyx=0,及流体边界没有流体流出,因此边界条件为:
下边界为固定边界,需要同时满足两个条件:边界法向上流体加速度为0,固体位移是0,因此边界条件是:
上边界加载周期震荡应力,对应的边界条件是边界法向应力是σzz,流体流动的加速度是0,因此边界条件是:
S7:全局粗化网格剖分,局部网格精细剖分求解,具体包括:
S71:利用有限元算法使用粗化网格求解方程;
S72:计算所有网格元素上的偏微分方程的残差;
S73:估计所有网格元素的解中误差,其中第1方程的L2范数误差为h为局部网格尺寸,ρl为第1个偏微分方程的残差,qi为偏微分方程稳定性估计导数阶数,默认值是2,由于参数C是任意常数,计算出的误差实际上是误差指示因子;
S74:如果已进行了预设的优化次数或超过了最大有限元剖分数,则终止执行;
S76:重复步骤S71-S75;
S8:使用有限元求解方程,获得应变固体位移场和流体压力场;
S9:计算纵波模量M:
其中,σzz为垂向应力,εzz为垂向应变;
则纵波速度Vp和纵波衰减Qp为:
本发明对实验中气/水部分饱和Indiana石灰岩中流体流动引起的频散和衰减进行预测,表1为气和水的流体物性参数,表2为干燥Indiana岩石物性参数,样品饱和度为88%,其中,步骤S1水/气的流体分布模式通过CT扫描获得,并显示在图3中。
样品中为气/水两相,按照步骤S2-S4,获得图4(a)所示的流体体积模量分布函数、图4(b)所示的流体密度分布函数、图4(c)所示的流体黏度分布函数。
按照步骤S6,约束边界条件为:
如图5所示,横向边界条件为:
如图6所示,上边界加载周期震荡应力,对应的边界条件是边界法向应力及流体流动的加速度是0,边界条件为:
如图7所示,下边界为固定边界,需要同时满足两个条件:边界法向上流体加速度为0,固体位移是0,因此边界条件是:
为求解方程,首先利用图8所示的粗化网格进行计算,根据残差构建图9所示的自适应网格。求解方程,获得Indiana灰岩在不同频率下的流体压力分布,图10展示了频率为1000Hz时的流体压力分布。根据步骤S9,获得频变体积模量和衰减,如图11和图12所示。
为说明本发明的效果,以附图1(a)层状模型为例,该模型上层为空气,下层为水,表1为对应的流体物性。表2为模型岩石物性信息。结合上述物性信息,利用常规有限元计算方法求解Biot流-固耦合位移方程,为保证分层界面的精度,采用细化四面体网格剖分模型,如图1(b)所示,计算结果获得的流体压场如图1(c)所示,体积模量频散和衰减变化特征如图2中实线所示,计算过程总耗时为2分32秒。
使用本发明方法计算时,首先将所需物性信息按照步骤S1-S4构建系数矩阵,结合步骤S5的Biot位移-压力方程,粗化网格剖分,如图1(d)所示。在获得初步残差后,利用自适应算法细化小目标地质体网格,如图1(e)所示,可以看到,在边界层处网格充分细化,但其他位置网格仍然为粗化网格。据此计算的流体压力场为图1(f)。所得流体压力场图1(c)和图1(f)差距不大,其对应的频散和衰减曲线为图2虚线所示,对比常规算法结果(图2实线),可以看到两者的计算结果基本一致,但本发明所使用的计算时间为32秒,比常规方法速度提高五倍之多。
表1流体物性参数表
流体属性 | 数值 |
水体积模量 | 2.25E9Pa |
气体积模量 | 1E5Pa |
水密度 | 1000kg/m<sup>3</sup> |
气密度 | 78kg/m<sup>3</sup> |
水黏度 | 0.001Pa·s |
气黏度 | 1.5E-4Pa·s |
表2岩石物性参数表
固体属性 | 数值 |
岩石排水体积模量 | 2.5E10Pa |
骨架颗粒模量 | 7.7E10Pa |
孔隙度 | 0.108 |
干骨架密度 | 2367.4kg/m<sup>3</sup> |
渗透率 | 9.8692E-18m<sup>2</sup> |
岩石排水剪切模量 | 1.52E10Pa |
扭曲因子 | 1 |
模型长度 | 0.4m |
模型高度 | 0.4m |
Claims (3)
1.一种预测多相孔隙介质弹性模量的高效数值模拟方法,其特征在于,包括如下步骤:
S5:构建频率域Biot方程的位移-压力形式:
其中,;;为饱和岩石密度;;为岩石骨架密度,为孔隙度,为流体复密度,为扭曲因子,为岩石渗透率;为总应力,C为干燥岩石骨架柔度,为应变张量,为Biot-Willis系数,为流体压力,I为单位张量矩阵,F为周期震荡的外力;为角频率,为固体位移张量;
下边界为固定边界,需要同时满足两个条件:边界法向上流体加速度为0,固体位移是0,因此边界条件是:
S7:全局粗化网格剖分,局部网格精细剖分求解;
S8:使用有限元求解方程,获得应变固体位移场和流体压力场;
S9:计算纵波模量M:
2.根据权利要求1所述的预测多相孔隙介质弹性模量的高效数值模拟方法,其特征在于,步骤S1中流体分布图像通过CT扫描获得。
3.根据权利要求1所述的预测多相孔隙介质弹性模量的高效数值模拟方法,其特征在于,步骤S7具体包括:
S71:利用有限元算法使用粗化网格求解方程;
S72:计算所有网格元素上的偏微分方程的残差;
S73:估计所有网格元素的解中误差,其中第l方程的L2范数误差为,h为局部网格尺寸,为第l个偏微分方程的残差,为偏微分方程稳定性估计导数阶数,默认值是2,由于参数C是任意常数,计算出的误差实际上是误差指示因子;
S74:如果已进行了预设的优化次数或超过了最大有限元剖分数,则终止执行;
S76:重复步骤S71-S75。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111303690.4A CN114021498B (zh) | 2021-11-05 | 2021-11-05 | 一种预测多相孔隙介质弹性模量的高效数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111303690.4A CN114021498B (zh) | 2021-11-05 | 2021-11-05 | 一种预测多相孔隙介质弹性模量的高效数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114021498A CN114021498A (zh) | 2022-02-08 |
CN114021498B true CN114021498B (zh) | 2022-10-11 |
Family
ID=80061076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111303690.4A Active CN114021498B (zh) | 2021-11-05 | 2021-11-05 | 一种预测多相孔隙介质弹性模量的高效数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114021498B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115308798B (zh) * | 2022-08-26 | 2023-04-07 | 中国矿业大学 | 一种高速低内存消耗的预测储层岩石弹性波速度的方法 |
CN115575505B (zh) * | 2022-10-10 | 2024-02-13 | 四川大学 | 一种应力作用条件下岩石纵波波速与衰减的计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105607123A (zh) * | 2016-01-15 | 2016-05-25 | 中国石油天然气股份有限公司 | 一种计算随机孔隙介质模型地震波特征信息的方法和装置 |
CN109116420A (zh) * | 2018-10-16 | 2019-01-01 | 河海大学 | 一种含裂隙的孔隙介质纵波速度与衰减预测方法 |
CN109471168A (zh) * | 2018-11-06 | 2019-03-15 | 河海大学 | 一种孔裂隙介质中纵波速度与衰减的预测方法 |
CN110263434A (zh) * | 2019-06-20 | 2019-09-20 | 中国石油大学(华东) | 一种基于多尺度混合有限元的流动单元数值模拟方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10571603B2 (en) * | 2014-02-19 | 2020-02-25 | Repsol, S.A. | Method implemented in a computer for the numerical simulation of a porous medium |
-
2021
- 2021-11-05 CN CN202111303690.4A patent/CN114021498B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105607123A (zh) * | 2016-01-15 | 2016-05-25 | 中国石油天然气股份有限公司 | 一种计算随机孔隙介质模型地震波特征信息的方法和装置 |
CN109116420A (zh) * | 2018-10-16 | 2019-01-01 | 河海大学 | 一种含裂隙的孔隙介质纵波速度与衰减预测方法 |
CN109471168A (zh) * | 2018-11-06 | 2019-03-15 | 河海大学 | 一种孔裂隙介质中纵波速度与衰减的预测方法 |
CN110263434A (zh) * | 2019-06-20 | 2019-09-20 | 中国石油大学(华东) | 一种基于多尺度混合有限元的流动单元数值模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114021498A (zh) | 2022-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114021498B (zh) | 一种预测多相孔隙介质弹性模量的高效数值模拟方法 | |
RU2428739C2 (ru) | Укрупнение сетки для моделей коллекторов путем повторного использования расчетов потока, полученных на основе геологических моделей | |
Manwart et al. | Lattice-Boltzmann and finite-difference simulations for the permeability for three-dimensional porous media | |
US6826520B1 (en) | Method of upscaling permeability for unstructured grids | |
US6662109B2 (en) | Method of constraining by dynamic production data a fine model representative of the distribution in the reservoir of a physical quantity characteristic of the subsoil structure | |
US9665537B2 (en) | Method for generating a fractured reservoir mesh with a limited number of nodes in the matrix medium | |
US8868390B2 (en) | Method of developing a petroleum reservoir from a flow model calibrated through pressure and saturation map scaling | |
US20130246031A1 (en) | Constructing Geologic Models From Geologic Concepts | |
US8364447B2 (en) | Method, program and computer system for conciliating hydrocarbon reservoir model data | |
Srisutthiyakorn et al. | What is the role of tortuosity in the Kozeny-Carman equation? | |
US11327201B2 (en) | Porosity prediction based on effective stress | |
Rezapour et al. | Upscaling of geological models of oil reservoirs with unstructured grids using lifting-based graph wavelet transforms | |
Allgeyer et al. | Numerical approximation of the 3d hydrostatic Navier–Stokes system with free surface | |
US8718992B2 (en) | Method for history matching of a geological model comprising a sub-seismic fault network | |
Nunna et al. | Novel diffuse source pressure transient upscaling | |
CN109613624A (zh) | 一种储层岩石声电性质联合模拟方法 | |
Jones et al. | Effect of image scaling and segmentation in digital rock characterisation | |
CN105301642A (zh) | 非均匀孔隙岩石及其固态有机质体积含量确定方法及装置 | |
CN115392090B (zh) | 基于三维数字岩心的地震波频散和衰减特征的预测方法 | |
US20220146706A1 (en) | Computer implemented method for manipulating a mesh for the discretization of a reservoir domain | |
Wendebourg et al. | Modeling compaction and isostatic compensation in SEDSIM for basin analysis and subsurface fluid flow | |
CN115793092A (zh) | 海域天然气水合物分布形态诊断方法、系统、设备及介质 | |
Shiriyev | Discrete fracture network modeling in a carbon dioxide flooded heavy oil reservoir | |
CN117238393A (zh) | 一种赋存水合物的沉积物的声学特征模拟方法和装置 | |
EP3918381A1 (en) | Hydrocarbon flow simulation |
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 |