×

通过多层非线性消除增强的全隐式两相VT-flash组成流模拟。 (英语) Zbl 07524787号

概述:VT闪蒸计算具有可变的体积、温度和摩尔数规格,是传统PT闪蒸计算相行为建模的一种新的风险替代方法。迄今为止,VT闪光主要用作解决相平衡问题的独立计算,很少有工作将其应用于组分流动模拟。在本研究中,采用多层非线性消除方法,开发了一种改进的两相VT-闪速合成流算法。与之前的工作相比,新算法的鲁棒性和效率显著提高,因为非线性消除方法隐式地消除了局部较大的非线性,从而恢复了牛顿迭代的大步长。为了提高计算效率,使用自适应时间步长控制来调整时间步长。此外,使用改进的阴影区方法绕过了不必要的稳定性测试,该方法能够适应大时间步长下的实质性成分变化。通过大量的数值算例验证了所提出的VT-flash组合流算法的鲁棒性和有效性,该算法具有多层非线性消除功能。尽管收敛问题并没有完全解决(其根源在于相边界处的不可微平衡压力),但通过多层非线性消除,时间精细化的发生大大减少。研究还发现,多层非线性消除通常会略微增加牛顿迭代次数,但显著增大时间步长。因此,在多层非线性消除方法下,VT-闪蒸成分流模拟的整体计算效率得到了提高。

MSC公司:

7.6亿 流体力学基本方法
76平方米 多孔介质中的流动;过滤;渗流
65新元 偏微分方程边值问题的数值方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Coats,K.H.,《油藏模拟》(1987年佩赫第48章),(《石油工程手册》(1987))
[2] 格里森,M.G。;Durlowsky,L.J.,《油藏流体流动建模》,年。Rev.流体机械。,37, 211-238 (2005) ·Zbl 1082.76107号
[3] Chen,Z.,《油藏模拟:石油开采中的数学技术》(2007),SIAM:SIAM Philadelphia,Pennsylvania·Zbl 1167.86001号
[4] Fussell,L.T。;Fussell,D.D.,《组分油藏模型的迭代技术》,SPE J.,19,211-220(1979)
[5] Coats,K.H.,状态方程组分模型,SPE J.,20,363-376(1980)
[6] Nghiem,L.X。;Fong,D.K。;Aziz,K.,《状态方程的成分建模》,SPE J.,21,687-698(1981)
[7] Young,L.C。;Stephenson,R.E.,《油藏模拟的广义组分方法》,SPE J.,23,727-742(1983)
[8] 王,P。;Stenby,E.H.,合成油藏模拟中的非迭代闪蒸计算算法,流体相平衡。,95, 93-108 (1994)
[9] 拉斯穆森,C.P。;Krejbjerg,K。;Michelsen,M.L。;Bjurström,K.E.,《利用合成、瞬态模拟应用提高闪光计算的计算速度》(SPE年度技术会议和展览,SPE年度科技会议和展览(2003年),科罗拉多州丹佛)
[10] Voskov,D.V。;Tchelepi,H.A.,用于流动模拟的成分空间参数化,(SPE油藏模拟研讨会,SPE油藏仿真研讨会,德克萨斯州休斯顿(2007))
[11] Voskov,D.V。;Tchelepi,H.A.,混相驱替模拟的成分空间参数化,Transp。多孔介质,75,111-128(2008)
[12] Voskov,D.V。;Tchelepi,H.A.,基于Tie-simplex的任意相数混合物热力学平衡计算数学框架,流体相平衡。,283, 1-11 (2009)
[13] Belkadi,A。;Yan,W。;Michelsen,M.L。;Stenby,E.H.,《成分模拟中加速闪光计算的两种方法的比较》(SPE油藏模拟研讨会(2011年),Woodlands:Woodlands Texas)
[14] Yan,B。;王,Y。;Nasrabadi,H。;基洛,J.E。;Wu,K.,使用合成空间加速合成模拟的闪光计算,J.Pet。科学。工程,159,1000-1008(2017)
[15] Wu,Y。;科维茨,C。;Sun,S。;Salama,A.,《加速两相组分流模拟中的闪蒸计算——稀疏网格的应用》,J.Compute。物理。,285, 88-99 (2015) ·兹比尔1351.76227
[16] Wu,Y。;Chen,Z.,高维稀疏网格在闪光计算中的应用:从理论到实现,流体相平衡。,464, 22-31 (2018)
[17] 加西米,M。;易卜拉欣,A。;Gildin,E.,《使用双线性近似技术进行油藏模拟的降阶建模》,(SPE拉丁美洲和加勒比石油工程会议,SPE拉美和加勒比石油工程师会议,委内瑞拉马拉开波(2014))
[18] 胡,X。;Wu,S。;吴晓红。;徐,J。;张春生。;张,S。;Zikatanov,L.,《预处理与油藏模拟应用相结合》,多尺度模型。模拟。,11, 507-521 (2013) ·Zbl 1426.76269号
[19] 刘,H。;王凯。;Chen,Z。;Jordan,K.E.,并行分布式系统上高度非均质油藏模拟的高效多级预处理程序,(SPE油藏模拟研讨会,SPE油藏仿真研讨会,德克萨斯州休斯顿(2015))
[20] Bui,Q。;Elman,H。;Moulton,D.,多孔介质中多相流的代数多重网格预处理,SIAM J.Sci。计算。,39,S662-S680(2017)·Zbl 1392.65024号
[21] Bui,Q.M。;Wang,L。;Osei-Kuffuor,D.,多孔介质中相变两相流的代数多重网格预处理器,高级水资源。,114, 19-28 (2018)
[22] 基洛,J.E。;Bhogeswara,R.,在分布式内存并行计算机上模拟组成油藏现象,J.Pet。技术。,43, 1-368 (1991)
[23] Han,C。;德尔沙德,M。;塞珀努里,K。;Pope,G.A.,《一个完全隐式、并行、成分化学驱模拟器》(SPE年度技术会议和展览,SPE年度科技会议和展览会,德克萨斯州达拉斯(2005))
[24] 王,Y。;Killough,J.E.,基于储层模型过度分解的并行/成分模拟负载平衡新方法,SPE J.,19,304-315(2013)
[25] Yang,H。;Sun,S。;李毅。;Yang,C.,并行计算机上油藏模拟的可扩展全隐式框架,计算。方法应用。机械。工程,330,334-350(2018)·Zbl 1439.76103号
[26] Yang,H。;Sun,S。;李毅。;Yang,C.,多组分可压缩流完全隐式互补公式的并行油藏模拟器,计算。物理学。社区。,244, 2-12 (2019) ·Zbl 07674826号
[27] Sammon,P.H.,用于成分模拟的动态网格细化和合并,(SPE储层模拟研讨会。SPE储层模拟研讨会,德克萨斯州休斯顿(2003))
[28] 医学博士Jackson。;戈麦斯,J.L。;莫斯塔吉米,P。;珀西瓦尔,J.R。;托利特,B.S。;帕夫利迪斯,D。;疼痛,C.C。;谢赫,A.H。;Muggeridge,A.H。;Blunt,M.J.,《使用曲面、自适应非结构化网格和控制体积单元方法进行流动模拟的油藏建模》(SPE油藏模拟研讨会(2013),Woodlands:Woodlands Texas)
[29] 俄勒冈州波利夫卡。;Mikyška,J.,《多孔介质中使用恒定体积闪蒸和通量计算的成分建模,无需相位识别》,J.Compute。物理。,272, 149-169 (2014) ·Zbl 1349.76825号
[30] 金德罗娃,T。;Mikyška,J.,给定体积、温度和摩尔数下多相平衡计算的通用算法,流体相平衡。,393, 7-25 (2015)
[31] Hoteit,H。;Firoozabadi,A.,《组合非连续Galerkin和混合方法的成分建模》,SPE J.,11,19-34(2006)
[32] Hoteit,H。;Firoozabadi,A.,用间断Galerkin和混合方法对无传递函数的离散断裂介质进行成分建模,SPE J.,11,341-352(2006)
[33] Druetta,P。;特西,P。;Persis,C.D。;Picchioni,F.,《石油开采过程和油藏模拟方法》,高级化学。工程科学。,6, 39-435 (2016)
[34] 俄勒冈州波利夫卡。;Mikyška,J.,《使用半隐式格式对多孔介质中两相流进行成分建模》,IAENG J.Appl。数学。,45, 218-226 (2015) ·Zbl 1512.76111号
[35] 蔡,X.C。;Li,X.,高局部非线性问题的基于限制可加性Schwarz非线性消除的非精确牛顿方法,SIAM J.Sci。计算。,33, 746-762 (2011) ·Zbl 1227.65045号
[36] 黄,F.-N。;林,H.-L。;Cai,X.C.,基于二级非线性消除的不精确牛顿方法预条件及其在激波管流计算中的应用,电子。事务处理。数字。分析。,37, 239-251 (2010) ·Zbl 1205.65180号
[37] Yang,H。;杨,C。;Sun,S.,多孔介质中两相流问题的带非线性消除的Active-set缩减空间方法,SIAM J.Sci。计算。,38,B593-B618(2016)·兹比尔1383.76384
[38] Yang,H。;Sun,S。;杨,C.,多孔介质两相流变分不等式解的非线性预处理半光滑牛顿方法,J.Compute。物理。,332, 1-20 (2017) ·兹比尔1378.76115
[39] Yang,H。;李毅。;Sun,S.,带毛细现象地下三相流的非线性预处理约束保持算法,计算。方法应用。机械。工程,367,第113140条pp.(2020)·Zbl 1442.76097号
[40] Yang,H。;Hwang,F.-N.,高度局部非线性多分量偏微分方程系统的自适应非线性消除预处理不精确牛顿算法,应用。数字。数学。,133, 100-115 (2018) ·Zbl 1397.65245号
[41] 罗,L。;Shiu,W.S。;陈,R。;Cai,X.C.,《人类动脉狭窄血流问题的非线性消除预处理不精确牛顿法》,J.Compute。物理。,399,第108926条pp.(2019)·Zbl 1453.65109号
[42] 罗,L。;蔡,X.C。;严,Z。;徐,L。;Keyes,D.E.,三维稳态不可压缩流动问题的多层非线性消除预处理不精确牛顿方法,SIAM J.Sci。计算。,42,B1404-B1428(2020)·Zbl 1456.76104号
[43] 罗伦茨,J。;布雷,B.G。;Clark,C.R.,根据储层流体的成分计算其粘度,J.Pet。技术。,16, 1-171 (1964)
[44] 彭大勇。;罗宾逊,D.B.,《一种新的双常数状态方程》,《工业工程化学》。芬达姆。,15, 59-64 (1976)
[45] 金德罗娃,T。;Mikyška,J.,在给定体积、温度和摩尔数下计算两相平衡的快速稳健算法,流体相平衡。,353, 101-114 (2013)
[46] Lie,K.-A.,《使用MATLAB/GNU Octave:MATLAB油藏模拟工具箱(MRST)用户指南》(2019年),剑桥大学出版社:英国剑桥大学出版社·Zbl 1425.76001号
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配,并且可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。