CN113297777A - 碳酸盐岩油气藏酸化反应流多尺度数值模拟方法及系统 - Google Patents
碳酸盐岩油气藏酸化反应流多尺度数值模拟方法及系统 Download PDFInfo
- Publication number
- CN113297777A CN113297777A CN202110686452.XA CN202110686452A CN113297777A CN 113297777 A CN113297777 A CN 113297777A CN 202110686452 A CN202110686452 A CN 202110686452A CN 113297777 A CN113297777 A CN 113297777A
- Authority
- CN
- China
- Prior art keywords
- reaction
- scale
- numerical simulation
- flow
- carbonate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- 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
Abstract
本公开公开的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法及系统,包括:获取注入酸的物性信息和碳酸盐岩的物性信息;将获取的信息输入连续尺度反应‑流动模型中,获得酸化反应流数值模拟结果;其中,连续尺度反应‑流动模型包括尺度升级后微分方程组和描述渗透率、扩散系数随反应变化的关系式;尺度升级后微分方程组为对构建的酸岩反应流数学模型中的微尺度变量进行渐近展开获得;描述渗透率、扩散系数随反应变化的关系式通过对描述尺度升级后微分方程组中等效参数与微尺度变量关系的元胞问题进行求解获得。实现了对酸化反应流的准确和高效模拟。
Description
技术领域
本发明涉及油藏数值模拟技术领域,尤其涉及碳酸盐岩油气藏酸化反应流多尺度数值模拟方法及系统。
背景技术
本部分的陈述仅仅是提供了与本公开相关的背景技术信息,不必然构成在先技术。
碳酸盐岩油气藏分布广泛、储量丰富,酸化改造是碳酸盐岩油气井清除钻井污染,改善近井渗流能力,获得较高产能的常用措施,对于给定的油藏和酸液体系,找到最优的注入速度,使得在这一速度下酸溶蚀地层生成蚓孔,是决定酸化施工成功的关键。然而,由于尺度差异,岩心酸化实验得到的最优注入速度并不能直接应用于油藏尺度的酸化施工中,数值模拟可以方便地解决尺度限制的问题。
目前多采用连续尺度模型进行酸化数值模拟,但是该模型的应用还存在两个突出的问题:一方面,岩石的渗透率、孔隙半径以及孔隙比表面积与孔隙结构有关,由于酸岩反应持续地改变多孔介质的孔隙结构,因此岩石的物性参数也随着反应的进行而发生改变,如何将孔隙尺度上岩石物性参数随反应的变化等效成连续尺度模型的参数,是应用连续模型进行酸化优化设计的前提;另一方面,目前应用连续尺度模型进行酸化数值模拟时,只能计算岩心尺度的酸蚀现象,而将其应用于真实油藏尺度的酸化模拟时,给数值计算带来了很大挑战:为了准确计算碳酸盐岩酸化生成的溶蚀结构,数值模拟时用到的网格尺寸不能超过蚓孔直径对应的特征长度(约1~3cm)。而实际油藏酸化的目标区域是井筒周围米数量级的范围,所以对单个网格大小的限制条件导致了模拟油藏尺度酸化现象时需要用到海量的计算网格。应用常规数值算法求解这种规模网格数量的问题会导致数值计算工作量大、机时耗时长、对计算机内存要求高等问题,甚至在工程实际中无法被接受。
故发明人认为,现有的碳酸盐岩油气藏酸化数值模拟方法,不能保证数值模拟的准确性和高效性。
发明内容
本公开为了解决上述问题,提出了碳酸盐岩油气藏酸化反应流多尺度数值模拟方法及系统,既能保证酸化的数值模拟精度,又能有效减少计算量,提高数值模拟的速度,具备了数值模拟的准确性和高效性。
为实现上述目的,本公开采用如下技术方案:
第一方面,提出了碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,包括:
获取注入酸的物性信息和碳酸盐岩的物性信息;
将获取的信息输入连续尺度反应-流动模型中,获得酸化反应流数值模拟结果;
其中,连续尺度反应-流动模型包括尺度升级后微分方程组和描述渗透率、扩散系数随反应变化的关系式;尺度升级后微分方程组为对构建的酸岩反应流数学模型中的微尺度变量进行渐近展开获得;描述渗透率、扩散系数随反应变化的关系式通过对描述尺度升级后微分方程组中等效参数与微尺度变量关系的元胞问题进行求解获得。
第二方面,提出了碳酸盐岩油气藏酸化反应流多尺度数值模拟系统,包括:
信息获取模块,用于获取注入酸的物性信息和碳酸盐岩的物性信息;
酸化反应流数值模拟结果获取模块,用于将获取信息输入连续尺度反应-流动模型中,获得酸化反应流数值模拟结果;
其中,连续尺度反应-流动模型包括尺度升级后微分方程组和描述渗透率、扩散系数随反应变化的关系式;尺度升级后微分方程组为对构建的酸岩反应流数学模型中的微尺度变量进行渐近展开获得;描述渗透率、扩散系数随反应变化的关系式通过对描述尺度升级后微分方程组中等效参数与微尺度变量关系的元胞问题进行求解获得。
第三方面,提出了一种电子设备,包括存储器和处理器以及存储在存储器上并在处理器上运行的计算机指令,所述计算机指令被处理器运行时,完成碳酸盐岩油气藏酸化反应流多尺度数值模拟方法所述的步骤。
第四方面,提出了一种计算机可读存储介质,用于存储计算机指令,所述计算机指令被处理器执行时,完成碳酸盐岩油气藏酸化反应流多尺度数值模拟方法所述的步骤。
与现有技术相比,本公开的有益效果为:
1、本公开的连续尺度反应-流动模型为对酸岩反应流数学模型进行尺度升级得到,从而能够获得岩石物性参数随反应变化的关系式及等效系数的精确值,利用该连续尺度反应-流动模型进行酸化反应流数值模拟时,充分考虑了酸化过程中孔隙结构的变化,有效保证了数值模拟结果的准确性。
2、本公开应用模型降阶方法对连续尺度反应-流动模型进行数值求解,在保证计算精度的同时降低计算量,使得本公开进行酸化反应流数值模拟时,保证了数值模拟的准确性和高效性。
本发明附加方面的优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1为本公开实施例1公开方法的流程示意图;
图2为本公开实施例1公开的数字岩心图像;
图3为本公开实施例1公开的广义多尺度有限元方法的离线计算流程;
图4为本公开实施例1公开的广义多尺度有限元方法的在线计算流程;
图5为本公开实施例1公开的在不同注酸速度下生成的溶蚀模式(以孔隙度场表示)。
具体实施方式:
下面结合附图与实施例对本公开作进一步说明。
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在本公开中,术语如“上”、“下”、“左”、“右”、“前”、“后”、“竖直”、“水平”、“侧”、“底”等指示的方位或位置关系为基于附图所示的方位或位置关系,只是为了便于叙述本公开各部件或元件结构关系而确定的关系词,并非特指本公开中任一部件或元件,不能理解为对本公开的限制。
本公开中,术语如“固接”、“相连”、“连接”等应做广义理解,表示可以是固定连接,也可以是一体地连接或可拆卸连接;可以是直接相连,也可以通过中间媒介间接相连。对于本领域的相关科研或技术人员,可以根据具体情况确定上述术语在本公开中的具体含义,不能理解为对本公开的限制。
实施例1
在该实施例中,公开了碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,包括:
获取注入酸的物性信息和碳酸盐岩的物性信息;
将获取的信息输入连续尺度反应-流动模型中,获得酸化反应流数值模拟结果;
其中,连续尺度反应-流动模型包括尺度升级后微分方程组和描述渗透率、扩散系数随反应变化的关系式;尺度升级后微分方程组为对构建的酸岩反应流数学模型中的微尺度变量进行渐近展开获得;描述渗透率、扩散系数随反应变化的关系式通过对描述尺度升级后微分方程组中等效参数与微尺度变量关系的元胞问题进行求解获得。
进一步的,酸岩反应流数学模型基于碳酸盐岩的表征单元体,考虑流体在孔隙内的流动、溶质在液相中的传输以及液固界面随反应的变化构建。
进一步的,微尺度变量为与孔隙结构变化有关的变量,包括酸溶质浓度、流体压力、流体速度、水平集函数。
进一步的,采用扩展有限元方法对元胞问题进行求解,获得求解结果。
进一步的,对元胞问题进行求解的具体过程为:
对整体区域采用正交网格进行离散,对元胞问题中的水平集方程采用一阶迎风差分格式进行离散求解;
基于水平集方程求解的网格,通过添加斜对角面,将正交网格划分成四面体网格,在此网格上求解元胞节点自由度;
根据单元节点上水平集值,构建形函数;
基于单元的形函数项,构建单元刚度矩阵和荷载向量,在此基础上组装系统刚度矩阵和荷载向量,形成线性方程组;
对线性方程组进行求解,获得元胞问题的解。
进一步的,采用广义多尺度有限元方法对连续尺度反应-流动模型进行求解,获得酸化反应流数值模拟结果。
进一步的,注入酸的物性信息包括注入酸的粘度、浓度、酸溶质的扩散系数、酸盐反应速率;
碳酸盐岩的物性信息包括岩石初始渗透率、孔隙度、矿物成分。
对本实施例公开的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法进行详细说明。
应用数值模拟指导碳酸盐岩油气藏酸化施工时还存在两个突出的问题:一方面,岩石的渗透率、孔隙半径以及孔隙比表面积与孔隙结构有关,由于酸岩反应持续地改变多孔介质的孔隙结构,因此岩石的物性参数也随着反应的进行而发生改变,如何将孔隙尺度上岩石物性参数随反应的变化等效成连续尺度模型的参数,是应用连续模型进行酸化优化设计的前提;另一方面,目前应用连续尺度模型进行酸化数值模拟时,只能计算岩心尺度的酸蚀现象,而将其应用于真实油藏尺度的酸化模拟时,给数值计算带来了很大挑战:为了准确计算碳酸盐岩酸化生成的溶蚀结构,数值模拟时用到的网格尺寸不能超过蚓孔直径对应的特征长度(约1~3cm),而实际油藏酸化的目标区域是井筒周围米数量级的范围,所以对单个网格大小的限制条件导致了模拟油藏尺度酸化现象时需要用到海量的计算网格。应用常规数值算法求解这种规模网格数量的问题会导致数值计算工作量大、机时耗时长、对计算机内存要求高等问题,甚至在工程实际中无法被接受。
因此,本实施例公开了一种碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,将孔隙尺度上精确描述酸岩反应流的酸岩反应流数学模型进行尺度升级,获得连续尺度反应-流动模型中岩石物性参数随反应变化的关系式及等效系数的精确值,应用模型降阶方法对连续尺度反应-流动模型进行数值求解,在保证计算精度的同时降低计算量,建立一套对油藏尺度酸化反应流问题进行准确、高效求解的多尺度数值模拟方法。
本实施例公开的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,包括:
S1:获取注入酸的物性信息和碳酸盐岩的物性信息。
其中,注入酸的物性信息包括注入酸的粘度、浓度、酸溶质的扩散系数、酸盐反应速率等。
碳酸盐岩的物性信息包括岩石初始渗透率、孔隙度、矿物成分等。
S2:将S1获得的信息输入连续尺度反应-流动模型中,获得酸化反应流数值模拟结果。
其中,连续尺度反应-流动模型的构建过程,如图1所示,包括:
S21:构建碳酸盐岩多孔介质数字岩心,选定中心区域4×4×4个像素点组成的介质,计算其渗透率、比表面积;然后沿着x、y、z三个方向将选定的区域扩大,直到计算出的渗透率和比表面积值不再剧烈变化。将此区域作为表征单元体进行下述计算,如图2所示为表征单元体的数字岩心图像。
S22:基于S21得到的表征单元体,在三维数字岩心上,考虑流体在孔隙内的流动、溶质在液相中的传输、以及液固界面随反应的变化,建立微观尺度精确描述酸液在孔隙空间内流动并与岩石骨架发生化学反应的酸岩反应流数学模型。
具体的酸岩反应流数学模型为:
式中,u表示酸液的流动速度向量;P表示流体压力;μ表示注入酸粘度;C表示注入酸浓度;D表示酸溶质的扩散系数;t表示时间;n表示垂直于液固表面的单位法向量;Rm表示反应速率;表示岩石矿物体积分数;Vm表示矿物组分的摩尔体积;S表示液固表面的比表面积;kr表示反应速率常数;Cs表示酸液的饱和浓度。
S23:将酸溶质浓度、流体压力、流体速度、水平集函数等与孔隙结构变化有关的微尺度变量进行渐进展开,带入S22建立的酸岩反应流数学模型中,推导出尺度升级后微分方程组,即为连续尺度的控制方程,同时推导得到一组描述渐进展开项系数的微分方程组,即为元胞问题。
其中,尺度升级后微分方程组为:
式中,v=(vx,vy,vz)为达西速度向量;K为渗透率张量;μ为流体粘度;P为压力;φ为孔隙度;t为时间;Cf为酸在液相中的浓度;Cs为酸在液固表面的浓度;De为有效扩散张量;kc为传质系数;av为比表面积;ks为表面反应速度;ρs为岩石密度;α表示酸的溶解能力,定义为每摩尔酸反应能溶解的固体质量。
元胞问题为:
式中,j=1,2,3表示空间维度下标;Y表示元胞单元;下标l表示液相;Γ表示元胞单元中液相与固相的边界;ζ表示待求变量。
其中,液相与固相的界面通过水平集方程进行描述,方程为:
式中,L表示水平集,其意义为:
S24:对S23中的元胞问题进行求解,获得描述渗透率、扩散系数随反应变化的关系式。
具体为:应用扩展有限元方法对元胞问题进行数值求解,计算连续尺度模型中各等效参数在不同反应阶段对应的数值,基于数值计算的结果,拟合渗透率、扩散系数与孔隙度的关系式,获得描述渗透率、扩散系数随反应变化的关系式。
获得的渗透率及扩散系数随反应变化的关系式分别为:
D=6.23×10-10+0.667φ-1.102
式中,K为渗透率,单位为10-3μm;D为酸溶质的扩散系数,单位为m2/s。
应用扩展有限元方法对元胞问题进行数值求解的具体步骤包括:
(1)对整体区域采用正交网格进行离散,对元胞问题中的水平集方程采用一阶迎风差分格式进行离散求解。
(2)基于水平集方程求解的网格,通过添加斜对角面,将正交网格划分成四面体网格,在此网格上求解元胞节点自由度。
(3)根据单元节点上水平集值,构建形函数。如果单元与水平集零等值线不相交,采用标准有限元的形函数;如果单元与水平集零等值线相交,建立增强函数,构建形函数增强项。
(4)基于单元的形函数项,构建单元刚度矩阵和荷载向量,在此基础上组装系统刚度矩阵和荷载向量,形成线性方程组。
(5)对线性方程组进行求解,获得元胞问题的解。
S25:由S23获得的尺度升级后微分方程组和S24获得的渗透率及扩散系数随反应变化的关系式,加上相应的边界条件,获得连续尺度反应-流动模型。
具体的连续尺度反应-流动模型为:
将S1获取的信息输入S2构建的连续尺度反应-流动模型中,并输入不同的注酸速度,采用广义多尺度有限元方法对连续尺度反应-流动模型进行求解,获得不同注酸速度下的酸化反应流数值模拟结果,并根据不同注酸速度下的酸化反应流数值模拟结果,确定最优注酸速度。
采用广义多尺度有限元方法对连续尺度反应-流动模型进行求解,求解分为离线计算和在线计算两个阶段。
离线计算的算法流程如图3所示,具体步骤为:
①生成粗网格。对连续尺度反应-流动网格进行粗化,每个粗网格均包含多个细网格。
②构建样本空间用于计算离线空间。在每个粗网格上,采用调和延拓来构建样本空间。
③构建离线空间。对任意粗网格单元,采用谱分解方法对相应的样本空间进行正交化。选取从小到大排列在前n位(n为每个粗网格上的基函数个数,根据计算精度需求,n可为1个或多个,n越大精度越高,本算例中n=6)的特征值对应的特征向量来生成离线空间,使得样本空间中的任意元素都能通过离线空间中的n个向量近似地线性表出,即:
在线计算的算法流程如图4所示,具体步骤为:
①构建多尺度基函数空间。多尺度基函数空间的构建方法与离线空间的构建类似,将离线计算阶段构建的离线空间看作样本空间,同时考虑实际模型中边界条件、源汇项等参数的变化,采用谱分解方法进一步降低其维数可得到多尺度基函数空间。
②求解全局问题。给定物理问题的计算参数和初边值条件,包括注入酸的物性信息和碳酸盐岩的初始物性信息,将全局问题控制方程在基函数空间上投影,组装粗网格上的刚度矩阵,求解得到的线性方程组,即为待求方程在粗网格上的解,记为Pc。
③根据粗网格解与细网格解的对应关系,求得微分方程在细网格上的全局解Pf。
Pf=DTPc
式中,D为离线计算得到的线性表出系数dij组成的矩阵。
在输入连续尺度反应-流动模型中其余参数不变,仅改变酸的注入速度时,可以计算得到的不同注酸速度下生成的溶蚀模式,如图5所示。
本实施例公开的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,将孔隙尺度上精确描述酸岩反应流的酸岩反应流数学模型进行尺度升级,获得连续尺度反应-流动模型中岩石物性参数随反应变化的关系式及等效系数的精确值。应用模型降阶方法对连续尺度反应-流动模型进行数值求解,在保证计算精度的同时降低计算量,建立一套对油藏尺度酸化反应流问题进行准确、高效求解的多尺度数值模拟方法。
实施例2
在该实施例中,公开了碳酸盐岩油气藏酸化反应流多尺度数值模拟系统,包括:
信息获取模块,用于获取注入酸的物性信息和碳酸盐岩的物性信息;
酸化反应流数值模拟结果获取模块,用于将获取信息输入连续尺度反应-流动模型中,并输入不同的注酸速度,获得不同注酸速度下的酸化反应流数值模拟结果,根据不同注酸速度下的酸化反应流数值模拟结果,确定最优注酸速度;
其中,连续尺度反应-流动模型包括尺度升级后微分方程组和描述渗透率、扩散系数随反应变化的关系式;尺度升级后微分方程组为对构建的酸岩反应流数学模型中的微尺度变量进行渐近展开获得;描述渗透率、扩散系数随反应变化的关系式通过对描述尺度升级后微分方程组中等效参数与微尺度变量关系的元胞问题进行求解获得。
实施例3
在该实施例中,公开了一种电子设备,包括存储器和处理器以及存储在存储器上并在处理器上运行的计算机指令,所述计算机指令被处理器运行时,完成实施例1公开的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法所述的步骤。
实施例4
在该实施例中,公开了一种计算机可读存储介质,用于存储计算机指令,所述计算机指令被处理器执行时,完成实施例1公开的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法所述的步骤。
以上仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (10)
1.碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,其特征在于,包括:
获取注入酸的物性信息和碳酸盐岩的物性信息;
将获取的信息输入连续尺度反应-流动模型中,获得酸化反应流数值模拟结果;
其中,连续尺度反应-流动模型包括尺度升级后微分方程组和描述渗透率、扩散系数随反应变化的关系式;尺度升级后微分方程组为对构建的酸岩反应流数学模型中的微尺度变量进行渐近展开获得;描述渗透率、扩散系数随反应变化的关系式通过对描述尺度升级后微分方程组中等效参数与微尺度变量关系的元胞问题进行求解获得。
2.如权利要求1所述的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,其特征在于,酸岩反应流数学模型基于碳酸盐岩的表征单元体,考虑流体在孔隙内的流动、溶质在液相中的传输以及液固界面随反应的变化构建。
3.如权利要求1所述的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,其特征在于,微尺度变量为与孔隙结构变化有关的变量,包括酸溶质浓度、流体压力、流体速度、水平集函数。
4.如权利要求1所述的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,其特征在于,采用扩展有限元方法对元胞问题进行求解,获得求解结果。
5.如权利要求1所述的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,其特征在于,对元胞问题进行求解的具体过程为:
对整体区域采用正交网格进行离散,对元胞问题中的水平集方程采用一阶迎风差分格式进行离散求解;
基于水平集方程求解的网格,通过添加斜对角面,将正交网格划分成四面体网格,在此网格上求解元胞节点自由度;
根据单元节点上水平集值,构建形函数;
基于单元的形函数项,构建单元刚度矩阵和荷载向量,在此基础上组装系统刚度矩阵和荷载向量,形成线性方程组;
对线性方程组进行求解,获得元胞问题的解。
6.如权利要求1所述的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,其特征在于,采用广义多尺度有限元方法对连续尺度反应-流动模型进行求解,获得酸化反应流数值模拟结果。
7.如权利要求1所述的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法,其特征在于,注入酸的物性信息包括注入酸的粘度、浓度、酸溶质的扩散系数、酸盐反应速率;
碳酸盐岩的物性信息包括岩石初始渗透率、孔隙度、矿物成分。
8.碳酸盐岩油气藏酸化反应流多尺度数值模拟系统,其特征在于,包括:
信息获取模块,用于获取注入酸的物性信息和碳酸盐岩的物性信息;
酸化反应流数值模拟结果获取模块,用于将获取信息输入连续尺度反应-流动模型中,获得酸化反应流数值模拟结果;
其中,连续尺度反应-流动模型包括尺度升级后微分方程组和描述渗透率、扩散系数随反应变化的关系式;尺度升级后微分方程组为对构建的酸岩反应流数学模型中的微尺度变量进行渐近展开获得;描述渗透率、扩散系数随反应变化的关系式通过对描述尺度升级后微分方程组中等效参数与微尺度变量关系的元胞问题进行求解获得。
9.一种电子设备,其特征在于,包括存储器和处理器以及存储在存储器上并在处理器上运行的计算机指令,所述计算机指令被处理器运行时,完成权利要求1-7任一项所述的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法的步骤。
10.一种计算机可读存储介质,其特征在于,用于存储计算机指令,所述计算机指令被处理器执行时,完成权利要求1-7任一项所述的碳酸盐岩油气藏酸化反应流多尺度数值模拟方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110686452.XA CN113297777B (zh) | 2021-06-21 | 2021-06-21 | 碳酸盐岩油气藏酸化反应流多尺度数值模拟方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110686452.XA CN113297777B (zh) | 2021-06-21 | 2021-06-21 | 碳酸盐岩油气藏酸化反应流多尺度数值模拟方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113297777A true CN113297777A (zh) | 2021-08-24 |
CN113297777B CN113297777B (zh) | 2023-03-24 |
Family
ID=77328979
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110686452.XA Active CN113297777B (zh) | 2021-06-21 | 2021-06-21 | 碳酸盐岩油气藏酸化反应流多尺度数值模拟方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113297777B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114117791A (zh) * | 2021-11-26 | 2022-03-01 | 西安石油大学 | 一种碳酸盐岩酸化压裂数值模拟方法 |
CN114724641A (zh) * | 2022-06-07 | 2022-07-08 | 中国石油大学(华东) | 一种双矿物组分碳酸盐岩储层酸化数值模拟方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030225521A1 (en) * | 2002-05-31 | 2003-12-04 | Mohan Panga | Modeling, simulation and comparison of models for wormhole formation during matrix stimulation of carbonates |
CN110263434A (zh) * | 2019-06-20 | 2019-09-20 | 中国石油大学(华东) | 一种基于多尺度混合有限元的流动单元数值模拟方法 |
CN110826269A (zh) * | 2019-10-25 | 2020-02-21 | 衢州学院 | 基于不规则元胞的温度-结构场耦合拓扑优化设计方法 |
CN111553108A (zh) * | 2020-05-20 | 2020-08-18 | 中国石油大学(华东) | 一种页岩气藏流固耦合多尺度数值模拟方法 |
CN111814364A (zh) * | 2019-04-10 | 2020-10-23 | 中国石油化工股份有限公司 | 一种岩溶储层演化数值模拟方法 |
CN112069654A (zh) * | 2020-08-04 | 2020-12-11 | 中国石油大学(华东) | 一种碳酸盐岩酸化数值模拟方法 |
-
2021
- 2021-06-21 CN CN202110686452.XA patent/CN113297777B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030225521A1 (en) * | 2002-05-31 | 2003-12-04 | Mohan Panga | Modeling, simulation and comparison of models for wormhole formation during matrix stimulation of carbonates |
EP1509674A1 (en) * | 2002-05-31 | 2005-03-02 | Services Petroliers Schlumberger | Modeling, simulation and comparison of models for wormhole formation during matrix stimulation of carbonates |
CN1656299A (zh) * | 2002-05-31 | 2005-08-17 | 施蓝姆伯格技术公司 | 碳酸盐的基岩刺激中酸蚀孔洞形成的建模、模拟和模型比较 |
CN111814364A (zh) * | 2019-04-10 | 2020-10-23 | 中国石油化工股份有限公司 | 一种岩溶储层演化数值模拟方法 |
CN110263434A (zh) * | 2019-06-20 | 2019-09-20 | 中国石油大学(华东) | 一种基于多尺度混合有限元的流动单元数值模拟方法 |
CN110826269A (zh) * | 2019-10-25 | 2020-02-21 | 衢州学院 | 基于不规则元胞的温度-结构场耦合拓扑优化设计方法 |
CN111553108A (zh) * | 2020-05-20 | 2020-08-18 | 中国石油大学(华东) | 一种页岩气藏流固耦合多尺度数值模拟方法 |
CN112069654A (zh) * | 2020-08-04 | 2020-12-11 | 中国石油大学(华东) | 一种碳酸盐岩酸化数值模拟方法 |
Non-Patent Citations (1)
Title |
---|
PIYANG LIU ET AL.: "A general method for simulating reactive dissolution in carbonate rocks with arbitrary geometry", 《COMPUTATIONAL GEOSCIENCES》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114117791A (zh) * | 2021-11-26 | 2022-03-01 | 西安石油大学 | 一种碳酸盐岩酸化压裂数值模拟方法 |
CN114724641A (zh) * | 2022-06-07 | 2022-07-08 | 中国石油大学(华东) | 一种双矿物组分碳酸盐岩储层酸化数值模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113297777B (zh) | 2023-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113297777B (zh) | 碳酸盐岩油气藏酸化反应流多尺度数值模拟方法及系统 | |
CN102187342B (zh) | 多孔介质中多相流的基于间接误差的动态升级 | |
Hua et al. | Level set, phase-field, and immersed boundary methods for two-phase fluid flows | |
Loseille et al. | Boundary layer mesh generation and adaptivity | |
Liang | A simplified adaptive Cartesian grid system for solving the 2D shallow water equations | |
CN107038270B (zh) | 一种表面加工残余应力场引起的加工变形的计算方法 | |
Hørsholt et al. | Oil production optimization of black-oil models by integration of matlab and eclipse e300 | |
Durasov et al. | Debosh: Deep bayesian shape optimization | |
Capolei et al. | Oil reservoir production optimization using single shooting and ESDIRK methods | |
Burbulla et al. | A finite-volume moving-mesh method for two-phase flow in dynamically fracturing porous media | |
Yu et al. | TWO-GRID FINITE ELEMENT METHOD FOR THE STABILIZATION OF MIXED STOKES-DARCY MODEL. | |
Ewing | Upscaling of biological processes and multiphase flow in porous media | |
CN107977732A (zh) | 边水气藏水平井的见水时间预测方法及装置 | |
Van Heijst et al. | The spin-up of fluid in a rectangular container with sloping bottom | |
Goto et al. | Computing with vortices: Bridging fluid dynamics and its information-processing capability | |
CN111062156B (zh) | 基于网格收敛因子的网格划分方法和装置 | |
CN111339702A (zh) | 油藏数值模拟等效井筒半径计算方法 | |
CN112580861B (zh) | 一种针对非常规油气藏生产优化问题的控制方法及系统 | |
Miao et al. | Numerical calculation of channel dredging volume using 3D digital stratum model | |
Karlsen et al. | Semi-Godunov schemes for multiphase flows in porous media | |
Krafczyk et al. | Lattice‐gas simulations of two‐phase flow in porous media | |
CN117217131B (zh) | 一种油藏数值模拟方法、装置、设备和存储介质 | |
Tyvand et al. | Free-surface flow of a fluid body with an inner circular cylinder in impulsive motion | |
CN115758917A (zh) | 一种基于iwind-lr的水体除藻工程效益模拟方法和装置 | |
CN114970385A (zh) | 一种三相流体相变渗流的数值模拟方法和系统 |
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 |