CN112861331A - 一种快速构建油气藏数值模拟器系数矩阵的方法 - Google Patents
一种快速构建油气藏数值模拟器系数矩阵的方法 Download PDFInfo
- Publication number
- CN112861331A CN112861331A CN202110116471.9A CN202110116471A CN112861331A CN 112861331 A CN112861331 A CN 112861331A CN 202110116471 A CN202110116471 A CN 202110116471A CN 112861331 A CN112861331 A CN 112861331A
- Authority
- CN
- China
- Prior art keywords
- grid
- oil
- hash table
- gas reservoir
- matrix
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种快速构建油气藏数值模拟器系数矩阵的方法,包括以下步骤:步骤1,获取油气藏基本数据;步骤2,建立单相油流渗流微分方程;步骤3,采用有限差分法对油气藏数学模型进行推导和离散;步骤4,对油气藏进行结构化网格剖分和网格编号;步骤5,建立哈希表存储网格坐标与编号关系;步骤6,对建立的哈希表结合每个键的长度进行重新排序,形成新的网格排序结构;步骤7,将哈希表按照对角线规律重新形成邻接图结构,组织网格间的连接关系,生成邻接表;步骤8,根据邻接表,组装渗流微分方程的矩阵方程,进行求解计算。本发明可以在一定程度上提高油气藏数值模拟器的求解速度,对于油气藏数值模拟器的编制具有实用价值。
Description
技术领域
本发明涉及油气藏数值模拟领域,具体涉及一种快速构建油气藏数值模拟器系数矩阵的方法。
背景技术
油气藏数值模拟源于石油工程师在不同开采条件下准确预测油气藏开发动态的需要。随着油气藏数值模拟技术所依据的数学模型变得越来越复杂,再加上油气藏的非均质性和岩石的各向异性引起的油气藏复杂性,流体性质和相对渗透率在平面上分布的变化,裂缝结构的复杂等,会造成网格模型的复杂,其构成的线性方程组也越来越大。往往完成一次算例,需要很多时间,对计算设备的性能要求也更高。而线性方程组的求解在一次模拟过程中耗时往往占了整个模拟时间的很大比重。因此,建立一种更加快速的方法或者更利于迭代消元的矩阵结构是有必要的。
发明内容
本发明主要是针对现有技术中的不足之处,本发明的目的是提供一种快速构建油气藏数值模拟器系数矩阵的方法。
为达到以上技术目的,本发明采用以下技术方案:
步骤1:获取计算所需的参数包括储层参数、网格参数及流体参数;
步骤2:根据物质守恒方程:
式中:ρ为流体密度;υx,υy,υz为x,y,z方向上流体流速;Ax,Ay,Az为x,y,z方向的流体流经截面积;qm为外部源汇项;Vb为控制单元总体积;αc为体积转换因子;φ为孔隙度;
及达西定律:
可推导出三维单相流体渗流方程:
式中:βc为传导率转换因子;Kx,Ky,Kz为x,y,z方向上的渗透率;μL为流体黏度;BL为流体体积系数;γL为流体重度;
步骤3:对步骤2中推导得到的非线性偏微分方程的空间交换项进行有限差分近似,得到其差分格式:
其中:
式中:i,j,k为直角坐标系中网格坐标;
上式中μ、B均为压力的函数,属于弱非线性项,采用中心加权取值,如下式所示:
然后,对每个网格的质量累积项进行线性化处理,使得
其中:
经过线性化处理后,构成形如AX=B的线性方程组。其中A为空间交换项差分后形成的邻接矩阵,B为累积项线性化后的常数数列,X为所求的未知量;
步骤4:将一个油藏划分为n个结构化网格,按总网格数n的网格进行自然排序,得出一个顺序排列的一维数列,代表网格编号;
步骤5:将一个油藏划分为n个结构化网格,按总网格数n的网格进行自然排序,得出一个顺序排列的一维数列,代表网格的顺序编号;以二维模型为例,nx表示模型的行数,ny表示列数,i表示第i行,j表示第j列,(i,j)表示某个网格所处的坐标位置,其自然排序的编号与坐标的关系为:
num=i×nx+j
步骤6:步骤5中得到的哈希表,每一个键对应的是矩阵的一条对角线,根据哈希表中每个键所对应数组的长度,按键的值从大到小,先对键值为偶数的数列,根据其数组长度,依次填入网格正常的排序,再用相同的步骤依次填入键值为奇数的数列,则可得到重新排序后的哈希表;
步骤7:将步骤6中重新排序后的哈希表以键的顺序为索引,将网格序号依次排入直角坐标系对角线中,即可构造出该方法中的网格序号的排列顺序,用于构造系数矩阵的连接关系,形成一个的邻接表;
步骤8:根据步骤7中的邻接表构造系数矩阵的连接关系,形成一个的邻接矩阵,该邻接矩阵则是线性方程组中A的形式,对应填入差分结果,结合直接法或迭代法进行求解,则可得到未知数的值。
有益效果:
与现有技术相比,本发明具有以下有益效果:
1、没有增加前处理中这一模块的增加运算时间,且网格数越大时会比原算法用时更少;
2、在求解线性方程组时计算速度比自然排序法构成的系数矩阵计算更快。
附图说明
图1为该方法下网格排序示意图(5*5);
图2为5*5网格形成的邻接矩阵二值图;
图3为3*3网格下的无向图;
图4为3*3网格下由邻接表推导出邻接矩阵示意图;
图5为不同网格数下两种方法前处理速度对比图;
图6为实施例1中51*51网格下生产200天的压力分布图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
实施例:
一种快速构建油气藏数值模拟器系数矩阵的方法,具体的计算方法包括以下步骤:
步骤1:获取计算所需的参数包括储层参数、网格参数及流体参数;
步骤2:根据物质守恒方程:
式中:ρ为流体密度;υx,υy,υz为x,y,z方向上流体流速;Ax,Ay,Az为x,y,z方向的流体流经截面积;qm为外部源汇项;Vb为控制单元总体积;αc为体积转换因子;φ为孔隙度;
及达西定律:
可推导出三维单相流体渗流方程:
式中:βc为传导率转换因子;Kx,Ky,Kz为x,y,z方向上的渗透率;μL为流体黏度;BL为流体体积系数;γL为流体重度;
步骤3:对步骤2中推导得到的非线性偏微分方程的空间交换项进行有限差分近似,得到其差分格式:
其中:
式中:i,j,k为直角坐标系中网格坐标;
上式中μ、B均为压力的函数,属于弱非线性项,采用中心加权取值,如下式所示:
然后,对每个网格的质量累积项进行线性化处理,使得
其中:
经过线性化处理后,构成形如AX=B的线性方程组。其中A为空间交换项差分后形成的邻接矩阵,B为累积项线性化后的常数数列,X为所求的未知量;
步骤4:将一个油藏划分为n个结构化网格,按总网格数n的网格进行自然排序,得出一个顺序排列的一维数列,代表网格编号;
步骤5:将一个油藏划分为n个结构化网格,按总网格数n的网格进行自然排序,得出一个顺序排列的一维数列,代表网格的顺序编号;以二维模型为例,nx表示模型的行数,ny表示列数,i表示第i行,j表示第j列,(i,j)表示某个网格所处的坐标位置,其自然排序的编号与坐标的关系为:
num=i×nx+j
步骤6:步骤5中得到的哈希表,每一个键对应的是矩阵的一条对角线,根据哈希表中每个键所对应数组的长度,按键的值从大到小,先对键值为偶数的数列,根据其数组长度,依次填入网格正常的排序,再用相同的步骤依次填入键值为奇数的数列,则可得到重新排序后的哈希表;
步骤7:将步骤6中重新排序后的哈希表以键的顺序为索引,将网格序号依次排入直角坐标系对角线中,即可构造出该方法中的网格序号的排列顺序,用于构造系数矩阵的连接关系,形成一个的邻接表;
步骤8:根据步骤7中的邻接表构造系数矩阵的连接关系,形成一个的邻接矩阵,该邻接矩阵则是线性方程组中A的形式,对应填入差分结果,结合直接法或迭代法进行求解,则可得到未知数的值。
实施例1:
本实施例以封闭地层的二维单相油藏为例,具体地层及流体数据如下:
地层原始压力,kPa | 岩石压缩系数,kPa-1 |
30000 | 1e-7 |
初始孔隙度,小数 | 原油压缩系数,无量纲 |
0.18 | 3.5e-6 |
原油粘度,Pa·s | 原油体积系数,无量纲 |
0.01 | 1.33 |
X方向绝对渗透率,D | Y方向绝对渗透率,D |
0.15 | 0.15 |
油藏X方向长度,m | 油藏Y方向长度,m |
500 | 500 |
网格数据如下表所示:
X方向网格数 | Y方向网格数 | 总网格数 |
51 | 51 | 2601 |
首先,输入上述油藏、流体及网格数据,并设置边界条件和模拟时间及时间步长。
接着,建立该油藏的二维单相渗流微分方程,并对其进行离散化,针对弱非线性项进行中心加权处理。对于强非线性项进行上游加权处理,得到其差分格式。
然后,对油藏进行网格剖分,按照网格划分方法对网格进行编号,并依据本发明提出的方法,按照对角线与横纵坐标的关系建立网格间的连接关系(为无向图),储存为邻接表,对应其差分格式的邻接矩阵。
最后,根据邻接矩阵组装矩阵,构成线性方程组,再利用时间步长及各参数值代入方程,进行迭代计算,直到模拟终止时间为止,即可得到各时间上油藏压力分布。
本发明所提供的一种快速构建油气藏数值模拟器系数矩阵的方法,在获取油气藏数值模拟的基本数据的基础上,基于物质平衡原理,结合达西定律和连续性方程,建立单相油流的基本渗流微分方程,并采用有限差分方法对油藏数学模型进行推导和离散,得到其差分格式。然后,对油气藏进行结构化网格剖分并进行网格编号,针对网格坐标与编号关系,建立哈希表存储其关系。按照所提出的方法对哈希表进行重新排序,形成新的网格排序结构。再将哈希表按照对角线规律重新形成邻接图结构,组织网格间的连接关系,生成邻接表。最后,根据生成的邻接表,组装渗流微分方程的矩阵方程,进行求解计算。本发明可以在一定程度上提高油气藏数值模拟器的求解速度,但编程实现该方法有较高的难度,这里提出了较为简易的实现方法,对于油气藏数值模拟器的编制具有实用价值。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (9)
1.一种快速构建油气藏数值模拟器系数矩阵的方法,其特征在于,所述的计算方法包括包括以下步骤:
步骤1:获取油气藏数值模拟的基本数据;
步骤2:基于物质平衡原理,结合达西定律和连续性方程,建立单相油流的基本渗流微分方程;
步骤3:采用有限差分方法对油藏数学模型进行推导和离散,得到其差分格式;
步骤4:对油气藏进行结构化网格剖分并进行网格编号;
步骤5:针对网格坐标与编号关系,建立哈希表存储其关系;
步骤6:根据步骤5中哈希表每个键所对应数组的长度,按键的值从小到大,先对键值为偶数的数列,根据其数组长度,依次填入网格正常的排序,再用同样的步骤依次填入键值为奇数的数列;
步骤7:将哈希表按照对角线规律重新排列,形成邻接图结构,并组织网格间的连接关系,生成邻接表;
步骤8:根据生成的邻接表,组装渗流微分方程的矩阵方程,进行求解计算。
2.根据权利要求1所述的方法,其特征在于,需要获取的参数包括储层参数、网格参数及流体参数。
5.根据权利要求1所述的方法,其特征在于,将一个油藏划分为n个结构化网格,按总网格数n的网格进行自然排序,得出一个顺序排列的一维数列,代表网格的顺序编号;以二维模型为例,其自然排序的编号与坐标的关系为:
num=i×nx+j
其中,nx表示模型的行数,ny表示列数,i表示第i行,j表示第j列,(i,j)表示某个网格所处的坐标位置。
6.根据权利要求1所述的方法,其特征在于,在步骤4的基础上,建立一个哈希表作为存储坐标与网格编号关系的容器,该哈希表的键为横纵坐标之和i+j,值为网格横纵坐标之和相等的序号组成的一个数列,即每一个i+j对应一条对角线上安放的网格序号所组成的数列,其中序号值为i×nx+j。并将每个键对应的数列都按从小到大的顺序排列。
7.根据权利要求1所述的方法,其特征在于,步骤5中得到的哈希表,每一个键对应的是矩阵的一条对角线,根据哈希表中每个键所对应数组的长度,按键的值从大到小,先对键值为偶数的数列,根据其数组长度,依次填入网格正常的排序,再用相同的步骤依次填入键值为奇数的数列,则可得到重新排序后的哈希表。
8.根据权利要求1所述的方法,其特征在于,将步骤6中重新排序后的哈希表以键的顺序为索引,将网格序号依次排入直角坐标系对角线中,即可构造出该方法中的网格序号排列顺序,用于构造系数矩阵的连接关系,形成一个邻接表。
9.根据权利要求1所述的方法,其特征在于,根据步骤7中的邻接表构造系数矩阵的连接关系,形成一个的邻接矩阵,该邻接矩阵则是线性方程组中A的形式,对应填入差分结果,结合直接法或迭代法进行求解,则可得到未知数的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110116471.9A CN112861331B (zh) | 2021-01-28 | 2021-01-28 | 一种快速构建油气藏数值模拟器系数矩阵的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110116471.9A CN112861331B (zh) | 2021-01-28 | 2021-01-28 | 一种快速构建油气藏数值模拟器系数矩阵的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112861331A true CN112861331A (zh) | 2021-05-28 |
CN112861331B CN112861331B (zh) | 2022-02-25 |
Family
ID=75987052
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110116471.9A Active CN112861331B (zh) | 2021-01-28 | 2021-01-28 | 一种快速构建油气藏数值模拟器系数矩阵的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112861331B (zh) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073772A (zh) * | 2011-01-18 | 2011-05-25 | 大庆油田有限责任公司 | 一种mdt压力数值测试模拟方法 |
WO2013143222A1 (zh) * | 2012-03-31 | 2013-10-03 | 中国科学技术大学 | 地震事件参数估计获取方法和***,地震事件搜索引擎 |
CN103559209A (zh) * | 2013-10-10 | 2014-02-05 | 河南大学 | 一种Voronoi Diagram与虚拟网格结合的高效空间最近邻查询方法 |
CN104112057A (zh) * | 2013-04-19 | 2014-10-22 | 中国石油化工股份有限公司 | 一种大尺度裂缝性油藏数值模拟方法 |
CN104335218A (zh) * | 2012-03-30 | 2015-02-04 | 爱迪德加拿大公司 | 使用基函数编码来保护可访问的*** |
CN105740521A (zh) * | 2016-01-25 | 2016-07-06 | 北京大学 | 油藏数值模拟***求解过程中的小网格剔除方法和装置 |
CN106294540A (zh) * | 2016-07-22 | 2017-01-04 | 长江大学 | 基于p‑stable局部敏感哈希检索数据样式的多点地质统计建模方法 |
CN107169227A (zh) * | 2017-06-09 | 2017-09-15 | 中国石油大学(华东) | 一种分段压裂水平井的粗网格模拟方法及*** |
CN109408838A (zh) * | 2017-08-16 | 2019-03-01 | 中国石油化工股份有限公司 | 一种对缝洞型油藏剩余油进行快速分析的方法及*** |
CN110321618A (zh) * | 2019-06-26 | 2019-10-11 | 中国石油大学胜利学院 | 一种低渗透油藏n2非混相驱数学模拟方法 |
CN111079335A (zh) * | 2019-12-20 | 2020-04-28 | 西南石油大学 | 一种嵌入式离散裂缝模型的数学推导和数值计算方法 |
CN111104766A (zh) * | 2019-12-16 | 2020-05-05 | 中国石油大学(华东) | 基于离散裂缝模型的油水两相非达西渗流数值模拟方法 |
WO2020224539A1 (zh) * | 2019-05-06 | 2020-11-12 | 西安华线石油科技有限公司 | 一种基于广义管流渗流耦合的流动模拟及瞬变井分析方法 |
CN111963163A (zh) * | 2020-09-16 | 2020-11-20 | 西南石油大学 | 一种基于气水两相流动页岩气藏双孔双渗模型构建及压力动态预测方法 |
-
2021
- 2021-01-28 CN CN202110116471.9A patent/CN112861331B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073772A (zh) * | 2011-01-18 | 2011-05-25 | 大庆油田有限责任公司 | 一种mdt压力数值测试模拟方法 |
CN104335218A (zh) * | 2012-03-30 | 2015-02-04 | 爱迪德加拿大公司 | 使用基函数编码来保护可访问的*** |
WO2013143222A1 (zh) * | 2012-03-31 | 2013-10-03 | 中国科学技术大学 | 地震事件参数估计获取方法和***,地震事件搜索引擎 |
CN104112057A (zh) * | 2013-04-19 | 2014-10-22 | 中国石油化工股份有限公司 | 一种大尺度裂缝性油藏数值模拟方法 |
CN103559209A (zh) * | 2013-10-10 | 2014-02-05 | 河南大学 | 一种Voronoi Diagram与虚拟网格结合的高效空间最近邻查询方法 |
CN105740521A (zh) * | 2016-01-25 | 2016-07-06 | 北京大学 | 油藏数值模拟***求解过程中的小网格剔除方法和装置 |
CN106294540A (zh) * | 2016-07-22 | 2017-01-04 | 长江大学 | 基于p‑stable局部敏感哈希检索数据样式的多点地质统计建模方法 |
CN107169227A (zh) * | 2017-06-09 | 2017-09-15 | 中国石油大学(华东) | 一种分段压裂水平井的粗网格模拟方法及*** |
CN109408838A (zh) * | 2017-08-16 | 2019-03-01 | 中国石油化工股份有限公司 | 一种对缝洞型油藏剩余油进行快速分析的方法及*** |
WO2020224539A1 (zh) * | 2019-05-06 | 2020-11-12 | 西安华线石油科技有限公司 | 一种基于广义管流渗流耦合的流动模拟及瞬变井分析方法 |
CN110321618A (zh) * | 2019-06-26 | 2019-10-11 | 中国石油大学胜利学院 | 一种低渗透油藏n2非混相驱数学模拟方法 |
CN111104766A (zh) * | 2019-12-16 | 2020-05-05 | 中国石油大学(华东) | 基于离散裂缝模型的油水两相非达西渗流数值模拟方法 |
CN111079335A (zh) * | 2019-12-20 | 2020-04-28 | 西南石油大学 | 一种嵌入式离散裂缝模型的数学推导和数值计算方法 |
CN111963163A (zh) * | 2020-09-16 | 2020-11-20 | 西南石油大学 | 一种基于气水两相流动页岩气藏双孔双渗模型构建及压力动态预测方法 |
Non-Patent Citations (5)
Title |
---|
姜瑞忠 等: "低渗透油藏非线性渗流新模型及试井分析", 《石油学报》 * |
张茂林 等: "三维四相多组分数学模型和数值模型的建立", 《断块油气田》 * |
王宝华 等: "大规模油藏数值模拟的块压缩存储及求解", 《石油勘探与开发》 * |
王英 等: "致密砂岩油藏注气非混相驱物质平衡及渗流机理研究", 《油气藏评价与开发》 * |
高大鹏 等: "精细分层注水约束的油藏数值模拟", 《计算物理》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112861331B (zh) | 2022-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Duddu et al. | A combined extended finite element and level set method for biofilm growth | |
CN111507024B (zh) | 基于gpu矩阵的离散元流固耦合数值模拟方法及*** | |
CN106991244B (zh) | 一种基于图论的裂隙网络连通性及渗流计算的方法 | |
CN111507026A (zh) | 一种模拟节点达西渗透流速的双重网格多尺度有限单元法 | |
CN114491831A (zh) | 一种基于断裂相场法的非均匀材料弥散裂纹j积分方法 | |
CN112100890B (zh) | 一种压裂水平缝起裂扩展全三维数学模型的计算方法 | |
CN112861331B (zh) | 一种快速构建油气藏数值模拟器系数矩阵的方法 | |
CN113011072B (zh) | 基于midas-pfc3d的离散元复杂模型识别方法 | |
CN107886573B (zh) | 一种复杂地质条件下边坡三维有限元网格生成方法 | |
CN111339688B (zh) | 基于大数据并行算法求解火箭仿真模型时域方程的方法 | |
Younesy et al. | Visualization of time-varying volumetric data using differential time-histogram table | |
CN103477345B (zh) | 多孔介质中流量的变量高保真模拟 | |
CN116822400A (zh) | 一种适合于大尺度平原河网一维非恒定流模拟方法 | |
CN115758911A (zh) | 融合点云残差网络和考虑滑移的流场和压力场预测方法 | |
CN104182558B (zh) | 一种缝洞野外露头水驱油数值模拟方法 | |
Nardean et al. | A novel and efficient preconditioner for solving Lagrange multipliers-based discretization schemes for reservoir simulations | |
CN115146408A (zh) | 一种基于三叉元结构的机床结构件正向设计方法 | |
CN115392032A (zh) | 一种gis-mpm无缝集成的动态三维地质模型构建方法 | |
CN108256143A (zh) | 一种飞行器表面多块结构网格的快速插值方法 | |
CN112307418A (zh) | 一种新型weno格式的高精度分数阶导数逼近方法 | |
CN112967147A (zh) | 一种考虑多尺度裂缝的基岩产量贡献率计算方法 | |
CN113656976B (zh) | 一种二维磁梯度张量快速数值模拟方法、装置和设备 | |
CN110765506A (zh) | 实体模型的多分辨率等几何拓扑优化方法 | |
CN110196983B (zh) | 一种基于配点理论的高维随机热传导问题谱分析方法 | |
CN118114540B (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 |