[发明专利]油藏数值模拟井模型的处理方法及其相应系统在审
申请号: | 201810711096.0 | 申请日: | 2018-07-03 |
公开(公告)号: | CN109948171A | 公开(公告)日: | 2019-06-28 |
发明(设计)人: | 张地大;李健 | 申请(专利权)人: | 北京默凯斯能源技术有限公司 |
主分类号: | G06F17/50 | 分类号: | G06F17/50;G06F17/16 |
代理公司: | 北京中企讯知识产权代理有限公司 11677 | 代理人: | 熊亮 |
地址: | 100102 北京市朝阳区利*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 矩阵 网格 油藏数值模拟 矩阵方程 稀疏矩阵 原始油藏 迭代法 箭头形 井模型 修正 求解 快速收敛 相关参数 演化过程 油田数据 耦合 分裂 收敛 赋予 | ||
本发明提供一种油藏数值模拟井模型的处理方法及其相应系统。所述处理方法,包括获取原始油藏模拟箭头形稀疏矩阵;分裂原始油藏模拟箭头形稀疏矩阵形成网格方程矩阵、网格方程修正矩阵、井方程矩阵及井方程修正矩阵,并将前述分裂形成的矩阵进行耦合;赋予矩阵方程的y初始值,并对矩阵方程中的网格方程修正矩阵利用迭代法进行求解得x值为第一值;将所述第一值代入利用迭代法进行求解得y值为第二值;判断所述第一值、第二值是否同时收敛的步骤。采用所述处理方法,能够使处理结果快速收敛,能够更加快捷地与其他相关参数结合得到油田数据随时间的演化过程。
技术领域
本发明属于油田开采技术领域,具体涉及一种油藏数值模拟井模型的处理方法及其相应系统。
背景技术
工程实施前需要理论模拟计算的数据指引,通过模拟计算得到油田开采过程各种产量以及相应的干预措施,得到油田开采的经济效益最大化需要的一系列操作和数据。
目前通常的做法是采用有限差分法将油藏方程时间和空间项离散后得到七对角的线性矩阵,井方程则通过射孔与网格相耦合,将数据放在网格矩阵的左边和下边,构成一个更大的箭头型稀疏矩阵,通过求解此稀疏矩阵便可得到油田数据随时间的演化过程。
而井方程数据破坏网格方程的七对角结构,可能造成主对角不占优,增大网格方程的条件数,使得直接求解收敛较慢甚至在给定的迭代步数内达不到收敛精度要求,尤其是在井数量较大的情况下则收敛越不容易,从而使油田数据随时间的演化过程的获得过程尤其迟滞。
发明内容
因此,本发明要解决的技术问题在于提供一种油藏数值模拟井模型的处理方法及其相应系统,使得原始油藏模拟箭头形稀疏矩阵能更快的收敛,在结果收敛的情况下,能够更加快捷地与其他相关参数结合得到油田数据随时间的演化过程,进而能够更加快捷地得到油田开采的经济效益最大化需要的一系列操作和数据。
为了解决上述问题,本发明提供一种油藏数值模拟井模型的处理方法,包括如下步骤:
a.获取原始油藏模拟箭头形稀疏矩阵A’((m+n)×(m+n));
b.分裂所述原始油藏模拟箭头形稀疏矩阵A’((m+n)×(m+n)),形成网格方程矩阵A(n×n)、网格方程修正矩阵R(n×m)、井方程修正矩阵D(m×n)及井方程矩阵H(m×m),并将所述网格方程矩阵A(n×n)、网格方程修正矩阵R(n×m)通过公式Ax=b-Ry(1)进行耦合,将所述井方程修正矩阵D(m×n)及井方程矩阵H(m×m)通过公式Hy=c-Dx(3)进行耦合,其中,b、c为常数向量;
c.赋予y初始值,并代入所述公式(1),利用迭代法进行求解得x值为第一值;
d.将所述第一值代入所述公式(2)中,直接求解得y值为第二值;
e.判断所述第一值、第二值是否同时收敛,若收敛,则结束计算;若不收敛,则将第二值作为y的初始值,重复步骤c、步骤d,直到通过公式(1)、公式(2)得出的x值、y值同时收敛。
优选地,步骤c中的y初始值为0或上一时间步的解。
优选地,所述迭代法为BICGstab。
优选地,步骤a中的所述原始油藏模拟箭头形稀疏矩阵A’((m+n)×(m+n))通过有限差分法离散三维油藏模拟平衡方程获得。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京默凯斯能源技术有限公司,未经北京默凯斯能源技术有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810711096.0/2.html,转载请声明来源钻瓜专利网。