×

非结构网格上的高阶FE-FV方法,用于非均质多孔介质中的输运和变粘度两相流。 (英语) Zbl 1349.76262号

小结:本文提出了在可渗透多孔介质中同时传输和吸附粘性物质的两相流的高阶数值模拟方法。提出的数值方案解决了模拟这一过程中的三大挑战。首先,组分传输与作用在载体相运动上的粘性力和毛细管力强烈耦合。毛细管部分的离散化特别困难,因为它对流动的影响产生非线性抛物线守恒方程。这些都适用于非线性有限元(FE),而毛细管对组分传输的贡献是一阶双曲线,经典FE不适用。我们通过对粘性部分使用带显式时间步长的有限体积(FVs)和对毛细管部分使用带隐式时间步法的有限元-有限体积(FEFV)组合格式,有效地解决了这一问题。其次,组件会经历流体动力弥散,区分数值弥散和物理弥散至关重要。我们为保持数值色散较低的相通量和分量通量开发了高阶公式,并将其与隐式有限元结合,以便充分考虑色散张量的非线性。第三,地下渗透介质显示出强烈的空间异质性,系数在许多数量级上都有变化,并且几何复杂性使得使用非结构化网格至关重要。在这项工作中,我们使用了节点中心FV,它将其解决流的能力与FE的灵活性结合在一起。文中给出了日益复杂的数值例子,证明了该方法的收敛性和鲁棒性,并证明了其对高度非均质、几何复杂的裂隙多孔介质的通用性。

MSC公司:

76米10 有限元方法在流体力学问题中的应用
76个M12 有限体积法在流体力学问题中的应用
76S05号 多孔介质中的流动;过滤;渗流
76Txx型 多相多组分流动
PDF格式BibTeX公司 XML格式引用
全文: DOI程序

参考文献:

[1] 胡贝尔,R。;Helmig,R.,《非均质多孔介质中的多相流:经典有限元方法与基于隐式压力-显式饱和度的混合有限元-有限体积法》,国际期刊Numer。液体方法,29,8,899-920(1999)·Zbl 0938.76053号
[2] 盖革,S。;罗伯茨,S。;Matthäi,S。;Zoppou,C。;Burri,A.,《将有限元和有限体积方法相结合,在高度非均质和结构复杂的地质介质中进行高效多相流模拟》,《地质流体》,4284-299(2004)
[3] Matthäi,S.K。;盖革,S。;罗伯茨,S.G。;Paluszny,A。;Belayneh,M。;伯里,A。;Mezentsev,A。;卢,H。;库姆·D·。;德里斯纳,T。;Heinrich,C.A.,《结构复杂油藏中多相流体流动的数值模拟》,(Jolley,S.J.;etal.,structurally complex reservoirs,vol.292(2007),Geol。社会规范出版物),405至429
[4] Durlowsky,L.,一种基于三角形的混合有限元-有限体积法,用于模拟穿过多孔介质的两相流,J.Compute。物理。,105, 252 (1993) ·Zbl 0768.76046号
[5] 蒙特亚古多,J。;Firoozabadi,A.,《裂缝介质注水模拟的控制体积模型:结合基质异质性和储层润湿性影响》,SPE J.,12,355-366(2007)
[6] Reichenberger,V。;雅各布斯,H。;Bastian,P。;Helmig,R.,《裂隙多孔介质中两相流的混合维有限体积法》,水资源部。,29, 1020-1036 (2006)
[7] 马丁·V。;贾夫雷,J。;Roberts,J.,《将裂缝和屏障建模为多孔介质中流动的界面》,SIAM J.Sci。计算。,26, 1667-1691 (2005) ·Zbl 1083.76058号
[8] Lee,S。;M.拉夫。;Jensen,C.,《天然裂缝地层中多长度尺度流动的分层建模》,《水资源》。决议,37,443-455(2001)
[9] Hoteit,H。;Firoozabadi,A.,不同毛细压力下非均质渗透介质中两相流的数值模拟,Adv.Water Resour。,31, 56-73 (2008)
[10] Niessner,J。;Helmig,R.,非均质多孔介质中三相三组分过程的多尺度建模,Adv.Water Resour。,30, 2309-2325 (2007)
[11] Strang,G.,《关于差分格式的构造和比较》,SIAM J.Numer。分析。,5, 506-517 (1968) ·Zbl 0184.38503号
[12] Nocedal,J。;Wright,S.,《数值优化》(1999),Springer Verlag:纽约Springer Verlag出版社·Zbl 0930.65067号
[13] Mallison,B。;Gerritsen先生。;Jessen,K。;Orr,F.,两相多组分流的高阶迎风格式,SPE J.,79691,297-311(2005)
[14] Mikys˘ka,J。;Firoozabadi,A.,稳健高效合成模拟的高阶方法实现,J.Compute。物理。,229, 2898-2913 (2010) ·Zbl 1307.76053号
[15] Karimi-Fard,M。;龚,B。;Durlowsky,L.,《根据详细的断裂特征生成粗尺度连续流模型》,《水资源》。第42号决议(2006年)
[16] 盖革,S。;Matthäi,S.K。;Niessner,J。;Helmig,R.,《裂缝性多孔介质中三组分三相流的Black-oil模拟》,SPE J.,14,338-354(2009)
[17] 摩尔加特,J。;孙,S。;Firoozabadi,A.,使用高阶有限元方法对重力三相流进行组合建模,《水资源》。第47号决议(2011年)
[18] 摩尔加特,J。;Firoozabadi,A.,非结构化和各向异性介质中Fickian扩散的高阶组分建模,高级水资源。,33, 951-968 (2010)
[19] Datta-Gupta,A。;Yoon,S。;Vasco,D.W。;Pope,G.A.,《井间示踪剂测试分区的反向建模:流线方法》,《水资源》。研究,38,10-1029(2002)
[20] 汗,S.A。;Pope,G.A。;Trangenstein,J.A.,《污染物清理问题和提高石油采收率的胶束/聚合物物理特性模型》,运输。《多孔医学》,24,35-79(1996)
[21] 华盛顿州西部。;Harwell,J.H.,《表面活性剂和地下修复》,环境。科学。技术。,26, 2324-2330 (1992)
[22] Sorbie,K.S。;Mackay,E.J.,注水中注入盐水、共生盐水和含水层盐水的混合及其与油田结垢的相关性,J.Petrol。科学。工程,27,85-106(2000)
[23] Mackay,E.,《预测原位硫酸盐垢沉积及其对生成离子浓度的影响》,《化学》。工程研究设计。,81, 326-332 (2003)
[24] 奥斯塔,T。;钢绞线,S。;马德兰,M。;Puntevold,T。;Korsnes,R.,《白垩中的海水:提高采收率和压实液》,SPE Res.Eval。工程,11648-654(2008)
[26] Javadpour,F.,《地质构造中的(CO_2)注入:从孔隙尺度过程确定宏观系数》,Transp。《多孔医学》,79,87-105(2009)
[27] 徐,T。;Sonnenthal,E。;Spycher,N。;Pruess,K.,TOUGHREACT-可变饱和地质介质中非等温多相反应地球化学运移的模拟程序:地热注入和(CO_2)地质封存的应用,计算。地质科学。,32, 145-165 (2006)
[28] Aavatsmark,I.,四边形网格多点通量近似介绍,计算。地质科学。,6405-432(2002年)·Zbl 1094.76550号
[29] Lie,K。;克罗格斯塔德,S。;利加登一世。;纳特维格,J。;尼尔森,H。;Skaflestad,B.,复杂网格上一致离散化的开源matlab实现,计算。地质科学。,1-26(2011年)
[30] 爱德华兹,M。;Zheng,H.,在结构化和非结构化网格上具有高度各向异性张量的拟正darcy-flux近似的双族,J.Compute。物理。,229, 594-625 (2010) ·Zbl 1253.76091号
[31] Edwards,M.,《基于流动网格的非结构化、控制体积分布、全张量有限体积格式》,计算。地质科学。,6, 433-452 (2002) ·Zbl 1036.76034号
[32] 拉明,S。;Edwards,M.,《高度畸变非结构化网格上多孔介质流动的高分辨率对流方案》,国际期刊数值。方法工程,761139-1158(2008)·Zbl 1195.76266号
[33] Edwards,M.,《二维结构化和非结构化网格上的高分辨率双曲线耦合椭圆含通量CVD格式》,国际期刊Numer。《液体方法》,51,1059-1077(2006)·Zbl 1158.76363号
[34] Edwards,M.,《三维结构化和非结构化网格上的高分辨率双曲线耦合椭圆含通量CVD方案》,国际期刊Numer。《液体方法》,51,1079-1095(2006)·Zbl 1158.76364号
[36] Sandwe,T。;贝尔雷,I。;Nordbotten,J.,《离散断裂矩阵模拟的有效多点通量近似方法》,J.Compute。物理。,231, 3784-3800 (2012) ·Zbl 1402.76131号
[37] Karimi-Fard,M。;杜洛夫斯基,L。;Aziz,K.,《适用于通用油藏模拟器的高效离散裂缝模型》,SPE J.,9,227-236(2004)
[38] Niessner,J。;赫尔米格,R。;雅各布斯,H。;Roberts,J.,非均质多孔介质两相流牛顿迭代中的界面条件和线性化方案,Adv.Water Resour。,28, 671-687 (2005)
[39] 布伦特,M.J。;刘凯。;Thiele,M.R.,预测水库流量的广义流线方法,Pet。地质科学。,2, 259-269 (1996)
[40] 金·M·J。;Datta-Gupta,A.,《流线模拟:当前观点》,《现场》,22,91-140(1998)
[41] Bear,J.,多孔介质中流体的动力学(1972),多佛:纽约州多佛·兹比尔1191.76001
[42] 马斯卡特,M.,《石油生产的物理原理》(1949年),麦格劳希尔出版社:纽约麦格劳·希尔出版社
[44] Van Genuchten,M.,预测非饱和土壤导水率的封闭式方程,土壤科学。《美国社会杂志》,44,892-898(1980)
[45] 施密德,K。;盖革,S。;Sorbie,K.,非混相两相流中吸附分散溶质的共电流和逆流吸胀分析解,计算。地质科学。,2012年6月16日·Zbl 1348.76171号
[46] Pope,G.A.,《分数流理论在提高原油采收率中的应用》,SPE J.,20,3,191-205(1980)
[47] 赫伯特,A。;Jackson,C。;Lever,D.,地下水流动和溶质运移与流体密度密切相关,水资源。Res.,241781-1795(1988)
[48] Nield,D。;Bejan,A.,多孔介质中的对流(2006),Springer Verlag·Zbl 1256.76004号
[49] Peaceman,D.W.,《数值油藏模拟基础》(1977年),爱思唯尔:爱思唯尔阿姆斯特丹
[50] 阿齐兹,K。;Settari,A.,《油藏模拟》(1979),应用科学出版社:英国伦敦应用科学出版社
[51] Gerritsen先生。;Durlowsky,L.,《油藏流体流动建模》,年。流体力学版次。,37, 211-238 (2005) ·兹比尔1082.76107
[52] Nayagum,D。;Schäfer,G。;Mosé,R.,使用混合和不连续有限元模拟多孔介质中的两相不可压缩流动,计算。地质科学。,8, 49-73 (2004) ·Zbl 1221.76119号
[53] Paluszny,A。;Matthäi,S。;Hohmeyer,M.,《复杂地质结构的混合有限元有限体积离散化和断裂岩石上演示的新模拟工作流》,地质流体,7186-208(2007)
[54] Durlowsky,L.,《混合和控制体积有限元逼近达西速度和相关量的精度》,《水资源》。研究,30965-973(1994)
[55] Ciarlet,P.G.,《椭圆问题的有限元方法》(1979年),北荷兰出版公司:北荷兰出版社,阿姆斯特丹,纽约,牛津
[56] Stüben,K.,《代数多重网格综述》,J.Compute。申请。数学。,128, 281-309 (2001) ·Zbl 0979.65111号
[57] Krechel,A。;Stüben,K.,基于子域分块的并行代数多重网格,并行计算。,27, 1009-1031 (2001) ·Zbl 0971.68215号
[58] Matthäi,S。;尼克·H。;疼痛,C。;Neuweiler,I.,《裂隙岩石中溶质运移的模拟:一种允许大时间步长的高阶精确有限元有限体积方法》,Transp。《多孔医学》,83,2,1-30(2009)
[60] 西格尔,P。;莫斯,R。;Ackerer,P。;Jaffre,J.,使用不连续和混合有限元组合求解对流-扩散方程,国际期刊数值。液体方法,24595-613(1997)·Zbl 0894.76041号
[61] 布伦特,M。;Rubin,B.,《石油油藏模拟的隐式流量限制方案》,J.Compute。物理。,102, 194-210 (1992) ·Zbl 0775.76109号
[62] 西莉亚,医学硕士。;布卢塔斯,E.T。;Zarba,L.R.,非饱和流动方程的一般质量守恒数值解,水资源。Res.,26,1483-1496(1990)
[63] 詹妮·P。;Tchelepi,H。;Lee,S.,带S型通量函数双曲守恒律的无条件收敛非线性解算器,J.Compute。物理。,228, 7497-7512 (2009) ·Zbl 1391.76553号
[64] Johansen,T。;Winther,R.,《模拟聚合物驱油的双曲守恒律方程组的黎曼问题的解》,SIAM J.Math。分析。,19, 541-566 (1988) ·Zbl 0658.35061号
[65] LeVeque,R.,《守恒定律的数值方法》(1992),Birkhäuser·Zbl 0847.65053号
[66] 巴菲特,T。;Clain,S.,非结构化网格的单斜和多斜MUSCL方法,J.Compute。物理。,229, 3745-3776 (2010) ·Zbl 1189.65204号
[67] Hubbard,M.E.,非结构网格上MUSCL型有限体积格式的多维坡度限制器,J.Compute。物理。,155, 54-74 (1999) ·Zbl 0934.65109号
[68] Leonard,B.,《不连续性对流建模的简单高精度分辨率程序》,国际数值杂志。方法工程,81291-1318(1988)·Zbl 0667.76125号
[69] Van Leer,B.,《走向最终保守差分格式》。II: 单调性和守恒结合在二阶格式中,J.Comput。物理。,14, 361-370 (1974) ·Zbl 0276.65055号
[71] Harten,A.,双曲守恒律的高分辨率格式,J.Comput。物理。,49, 357-393 (1983) ·Zbl 0565.65050号
[72] Sweby,P.,双曲守恒律使用通量限制器的高分辨率方案,SIAM J.Numer。分析。,21, 995-1011 (1984) ·Zbl 0565.65048号
[73] Arminjin,P。;Dervieux,A.,在二维任意有限元网格上构建TVD-类人工粘度,J.Compute。物理。,106, 176-198 (1993) ·Zbl 0771.65062号
[74] Jameson,A.,《气体动力学数值方案的分析和设计》。1:人工扩散、迎风偏置、限制器及其对精度和多重网格收敛的影响,国际计算杂志。流体动力学。,4, 171-218 (1995)
[75] 巴顿,P。;兰伯特,C。;Causon,D.,《非结构单元的正守恒高分辨率对流格式》,国际J·数值。方法工程,39,1821-1838(1996)·Zbl 0884.76048号
[76] 施密德,K.S。;盖革,S。;Sorbie,K.S.,《非混相两相流中溶质的同向和逆流渗透和分散的半解析解》,《水资源》。第47号决议(2011年)
[77] 斯皮瓦克,A。;价格,H。;Settari,A.,用变分法求解多维两相非混溶流动方程,SPE J.,17,27-41(1977)
[78] 库姆·D·。;Matthäi,S。;盖革,S。;Driesner,T.,求解复杂地质介质中流体流动的并行FE-FV方案,计算。地质科学。,34, 1697-1707 (2008)
[79] 尼克·H。;Matthái,S.,非均质介质中溶质运移嵌入不连续性的混合有限元有限体积法,渗流区J,10,299-312(2011)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。