×
索姆·谢斯;亚瑟·蒙科尔;拉米·尤尼斯

多孔介质中多相多组分流动全隐式模拟的局部线性系统。 （英文） Zbl 1434.86018号

计算。地质科学。 24，第2期，743-759（2020年）.
摘要：在求解完全隐式油藏模拟时间步长的过程中，考虑到计算机精度，经常会发现计算出的牛顿更新可能非常稀疏。这种稀疏性可以高达95%，并且可以在一次迭代到下一次迭代之间变化很大。在最近的工作中，在求解整个线性系统之前，开发了一个数学上合理的框架来预测稀疏模式。该理论仅限于多维非均匀环境中的一般标量非线性对流-扩散-反应问题。该理论被应用于减小在两相流的连续隐式时间步长期间计算的线性系统的尺寸。结果表明，线性化计算和线性求解过程可以局部化95%，同时保持精确的牛顿收敛性和最终解。受该方法取得巨大成功的启发，本文开发了算法扩展，以设计全隐式耦合多分量问题的定位算法。我们提出、应用并测试了一种新的算法，以解决从基于状态方程的合成模拟器中获得的双曲方程组。在SPE 10地质模型中，当应用于涉及六种热力学物种的各种全隐式流动和多组分输运问题时，根据系统中的局部性水平，观察到的计算工作量减少了6到49倍。我们将该算法应用于几种有无重力和毛细现象的注入和耗尽场景，以研究该方法对潜在异质性和复杂性的适应性和鲁棒性。我们证明，该算法能够实现高效和稳健的全分辨率全隐式模拟，而不会出现自适应离散化方法引入的误差或半隐式方法的稳定性问题。

引用于2文件

MSC公司：

86-10 地球物理相关问题的数学建模或模拟
86-08 地球物理问题的计算方法
76S05号 多孔介质中的流动；过滤；渗流
76T06型 液-液双组分流动

关键词：

局部线性系统;全隐式仿真;多孔介质

软件：

IPARS公司;纽顿图书馆
PDF格式BibTeX公司 XML格式引用
\textit{S.Sheth}等人，计算。地质科学。24，第2号，743--759（2020；Zbl 1434.86018）
全文： 内政部

参考文献：

[1] Allgower，E。；Böhmer，K。；Potra，F。；Rheinboldt，W.，算子方程及其离散化的网格独立性原理，印度工业大学。申请。数学。数字。分析。，23, 1, 160-169 (1986) ·Zbl 0591.65043号
[2] Axelsson，O.：迭代求解方法。剑桥大学出版社（1994）·Zbl 0795.65014号
[3] Aziz，K.，Settari，A.：油藏模拟。爱思唯尔应用科学（1979）
[4] 伯杰，M。；Colella，P.，《冲击流体动力学的局部自适应网格细化》，J.Compute。物理。，82, 1, 64-84 (1989) ·Zbl 0665.76070号 ·doi:10.1016/0021-9991（89）90035-1
[5] 蔡，X。；凯斯，D。；Marcinkowski，L.，《非线性加性schwarz预条件及其在计算流体动力学中的应用》，国际数值杂志。方法流体，40，12，1463-1470（2002）·Zbl 1025.76040号 ·doi:10.1002/fld.404
[6] Cao，H.、Tchelepi，H.，Wallis，J.、Yardumian，H.:油藏模拟的并行可扩展非结构化cpr型线性解算器。参加：SPE年度技术会议和展览（2005年）
[7] Chen，K.：矩阵预处理技术及其应用。剑桥大学出版社（2005）·Zbl 1079.65057号
[8] Chen，Z.，Huan，G.，Ma，Y.：多孔介质中多相流的计算方法。工业应用数学计算科学与工程学会（2006）·Zbl 1092.76001号
[9] 克里斯蒂，M。；Blunt，M.，《第十个spe比较解决方案项目：放大技术的比较》，Soc.Petroleum Eng.Reservoir Evalu。工程师，4，4308-317（2001）·doi:10.2118/72469-PA
[10] Coats，K.，《关于IMPES和一些基于IMPES的模拟模型的注释》，石油工程学会，5，245-251（2000）
[11] Dawson，C。；克利，H。；惠勒，M。；Woodward，C.，用预处理Newton-Krylov解算器求解两相流的一种并行、隐式、以单元为中心的方法，Compute。地质科学。，1, 3-4, 215-249 (1997) ·Zbl 0941.76062号 ·doi:10.1023/A:1011521413158
[12] Deufhard，P.：《非线性问题的牛顿方法：仿射不变性和自适应算法》，第35卷。施普林格（2011）·Zbl 1226.65043号
[13] Engstler，C。；Lubich，C.，不同时间尺度微分方程的多速率外推方法，计算，58，2，173-185（1997）·Zbl 0891.65083号 ·doi:10.1007/BF02684438
[14] 齿轮，C。；Wells，D.，多速率线性多步方法，BIT Numer。数学。，24, 4, 484-502 (1984) ·Zbl 0555.65046号 ·doi:10.1007/BF01934907
[15] Gerritsen，M。；Lambers，J.，非均质地层中地下水流模拟的局部-全局放大和网格适应性集成，计算。地质科学。，12, 2, 193-208 (2008) ·Zbl 1159.76362号 ·doi:10.1007/s10596-007-9078-2
[16] Hornung，R。；Trangenstein，J.，多孔介质流动的自适应网格细化和多级迭代，J.Compute。物理。，136,2522-545（1997年）·Zbl 0898.76074号 ·doi:10.1006/jcph.1997.5779
[17] 黄，F。；Cai，X.，不可压缩Navier-Stokes方程的并行非线性加法Schwarz预处理不精确牛顿算法，J.Compute。物理。，204, 2, 666-691 (2005) ·Zbl 1267.76083号 ·doi:10.1016/j.jcp.2004.10.025
[18] 凯利，C。；Sachs，E.，无限维问题牛顿类方法的网格独立性，J.积分Equ。申请。，3, 4, 549-573 (1991) ·Zbl 0756.65085号 ·doi:10.1216/jiea/1181075649
[19] Killough，J.E.，Kossack，C.A.：第五个比较解决方案项目：混相洪水模拟器的评估。摘自：SPE油藏模拟研讨会（1987年）
[20] Lacroix，S。；瓦西列夫斯基，Y。；Wheeler，M.，隐式并行精确油藏模拟器（ipars）中的解耦预处理器，编号。线性代数应用。，8, 8, 537-549 (2003) ·Zbl 1071.76583号 ·doi:10.1002/nla.264
[21] Moncorgé，A.，Tchelepi，H.，Jenny，P.：使用自然变量进行合成模拟的顺序全隐式公式。计算物理杂志（即将出版）。arXiv:1804.09751（2018）·兹比尔1415.76725
[22] Möyner，O.，Tchelepi，H.：等温状态方程组分问题的质量守恒序贯隐式多尺度方法。SPE期刊（即将出版）（2018年）
[23] 缪勒，S。；Stiriba，Y.，采用局部变化时间步长的守恒定律的完全自适应多尺度方案，J.Sci。计算。，30, 3, 493-531 (2007) ·Zbl 1110.76037号 ·数字对象标识代码：10.1007/s10915-006-9102-z
[24] 谢斯，S。；Younis，R.，全隐式两相流模拟中牛顿更新的局部计算，Procedia Compute。科学。，80, 1392-1403 (2016) ·doi:10.1016/j.procs.2016.05.445
[25] 谢斯，S。；Younis，R.，两相流和输运顺序隐式模拟中的局部线性系统，《石油工程学报》，22，5，1542-1569（2017）
[26] Sheth，S.，Younis，R.：通用全分辨率隐式油藏模拟的局部解算器。在：SPE油藏模拟会议（2017）
[27] 斯科格斯塔德，J。；Keilegavlen，E。；Nordbotten，J.，多孔介质非线性流动问题的区域分解策略，J.Compute。物理。，234, 439-451 (2012) ·doi:10.1016/j.jp.2012.10.001
[28] Sun，S.，Keyes，D.，Liu，L.：采用加性schwarz预处理非精确牛顿法的全隐式两相油藏模拟。参加：石油工程师学会油藏特征描述和模拟会议及展览（2013年）
[29] Trangenstein，J.，多孔介质流动的多尺度迭代技术和自适应网格细化，高级水资源。，25, 8, 1175-1213 (2002) ·doi:10.1016/S0309-1708（02）00053-2
[30] Wallis，J.、Kendall，R.、Little，T.等人：共轭残差法的约束残差加速度。摘自：SPE-13536-MS SPE油藏模拟研讨会论文集（1985）
[31] Younis，R.，关于一般两相流和输运现象时空局部性的一个尖锐分析界限，Procedia Compute。科学。，18, 473-480 (2013) ·doi:10.1016/j.procs.2013.05.211
[32] 尤尼斯，R。；Tchelepi，H。；Aziz，K.，《始终收敛的自适应局部连续-牛顿方法-非线性解算器》，SPE J.，15，2，526-544（2010）·doi:10.2118/119147-PA
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。