埃里克·凯莱加夫伦;Runar Berge;Alessio Fumagalli先生;米歇尔·斯塔诺尼;伊瓦尔·斯特凡森;贾布里尔·瓦雷拉;英加·伯雷 PorePy：一个用于模拟裂缝性多孔介质中多物理过程的开源软件。 （英语） Zbl 1453.86001号 计算。地质科学。 25，编号1，243-265（2021）. 摘要：随着研究向日益先进的过程耦合和复杂的断裂网络方向发展，裂隙岩石动力学模型和专用数值方法的开发是一个活跃的研究领域。模拟模型中包含耦合过程受到裂缝高纵横比、裂缝网络复杂几何形状以及完全改变裂缝-岩石界面特征的过程的关键影响的挑战。本文对在模拟器中引入断裂的设计原则进行了一般性讨论，并定义了一个集成建模、离散化和计算机实现的框架。该框架是在开源模拟软件PorePy中实现的，该软件可以作为裂隙岩石多物理问题的柔性原型工具。基于裂缝及其交点作为低维对象的表示，我们讨论了混合维网格的数据结构、多物理问题的表示以及利用现有软件的离散化。我们进一步提出了一个蟒蛇在PorePy开源软件工具中实现这些概念，该工具旨在耦合模拟三维裂缝性储层中的流动和运移以及裂缝和储层的变形。我们通过多孔介质中流动、孔隙弹性和裂缝变形的基准进行验证。然后，通过非线性耦合流动和传输以及注入驱动裂缝变形的模拟，说明了框架的灵活性。所有结果都可以通过公开可用的模拟脚本进行复制。 引用于28文件 MSC公司： 86-04 地球物理相关问题的软件、源代码等 86-08 地球物理问题的计算方法 关键词：裂缝性油气藏;混合维几何;数值模拟;多重物理量;离散断裂矩阵模型;开源软件;可再生科学 软件：驼鹿;Firedrake公司;蟒蛇;杜穆克斯;孔隙比;三角形;TOUGH（艰难）;PFLOTRAN公司;dfn工程;Gms小时;TetGen公司;MRST公司;FEniCS公司;Matlab公司;ParaView视图 PDF格式BibTeX公司 XML格式引用 \textit{E.Keilegavlen}等人，计算。地质科学。25，编号1，243--265（2021；Zbl 1453.86001） 全文： 内政部 arXiv公司 参考文献： [1] Berkowitz，B.，《表征断裂地质介质中的流动和传输：综述》，《高级水资源》。，25, 8-12, 861-884 (2002) ·doi:10.1016/S0309-1708（02）00042-8 [2] 马丁·V。；贾夫雷，J。；Roberts，JE，《将裂缝和屏障建模为多孔介质中流动的界面》，SIAM J.Sci。计算。，26, 5, 1667-1691 (2005) ·Zbl 1083.76058号 ·doi:10.137/S1064827503429363 [3] 巴顿，N。；Bandis，S。；Bakhtar，K.，《岩石节理的强度、变形和导电性耦合》，国际岩石力学杂志。最小科学。地质力学。文章摘要。，22, 3, 121-140 (1985) ·doi:10.1016/0148-9062（85）93227-9 [4] 北卡罗来纳州弗里。；JE罗伯茨；Saada，A.，《将断裂建模为界面：Forchheimer断裂模型》，计算。地质科学。，12, 1, 91-104 (2008) ·兹比尔1138.76062 ·数字对象标识代码：10.1007/s10596-007-9062-x [5] Rutqvist，J。；Wu，Y-S；Tsang，C-F；Bodvarsson，G.，《裂隙多孔岩石中耦合多相流体流动、传热和变形分析的建模方法》，《国际岩石力学杂志》。最小科学。，39, 4, 429-442 (2002) ·doi:10.1016/S1365-1609（02）00022-9 [6] Burnell，J.等人：地热超级模型：下一代综合地球物理、化学和流动模拟建模工具。Proc World地热。恭喜。7 (2015) [7] Pruess，K.：TOUGH2：多相流体和热流通用数值模拟程序。报告LBL-29400（1991） [8] 通用电气哈蒙德；Lichtner，PC公司；Mills，RT，评估并行地下模拟器的性能：PFLOTRAN的一个示例：评估PFLOTRAN的并行性能，Water Resour。研究，50，1，208-228（2014）·doi:10.1002/2012WR013483 [9] 巴伦布拉特，GI；Zheltov，IP；Kochina，IN，裂隙岩石[地层]中均质液体渗流理论的基本概念，J.Appl。数学。机械。，24, 5, 1286-1303 (1960) ·Zbl 0104.21702号 ·doi:10.1016/0021-8928（60）90107-6 [10] Arbogast，T。；小道格拉斯；Hornung，U.，《通过均匀化理论推导单相流的双重孔隙度模型》，SIAM J.Math。分析。，21, 4, 823-836 (1990) ·兹比尔0698.76106 ·doi:10.1137/0521046 [11] Lemonnier，P。；Bourbiaux，B.，《天然裂缝性储层的模拟》。最新技术：第1部分-物理机制和模拟器公式，石油天然气科学。Technol公司。《L’Institut Fr.Pétrole》，第65、2、239-262页（2010年）·doi:10.2516/ogst/2009066 [12] Lemonnier，P。；Bourbiaux，B.，《天然裂缝性储层的模拟》。最新技术：第2部分-基质-断裂转移和数值研究的典型特征，石油天然气科学。Technol.-技术《L’Institut Fr.Pétrole》，第65、2、263-286页（2010年）·doi:10.2516/ogst/2009067 [13] 海曼，JD；Karra，S。；Makedonska，N。；山墙，CW；画家，SL；Viswanathan，HS，dfnWorks：模拟地下水流和输送的离散裂缝网络框架，计算。地质科学。，84, 10-19 (2015) ·doi:10.1016/j.cageo.2015.08.001 [14] 埃赫尔，J。；德鲁齐，J-R；Poirriez，B.，三维离散裂缝网络中的流动模拟，SIAM J.Sci。计算。，31, 4, 2688-2705 (2009) ·Zbl 1387.65124号 ·电话：10.1137/080729244 [15] Berrone，S。；Pieraccini，S。；Scialó，S.，《利用基于优化的扩展有限元方法模拟离散裂隙网络流动》，SIAM J.Sci。计算。，35、2、A908-A935（2013）·Zbl 1266.65187号 ·doi:10.1137/120882883 [16] 贝尔雷，I。；加药器，F。；Keilegavlen，E.，《裂隙多孔介质中的流动：概念模型和离散方法综述》，Transp。多孔介质，130，215-236（2018）·doi:10.1007/s11242-018-1171-6 [17] Noorishad，J。；Mehran，M.，求解裂隙多孔介质中瞬态传输方程的上游有限元方法，水资源。研究，18，3，588-596（1982）·doi:10.1029/WR018i003p00588 [18] 巴卡，RG；阿内特，RC；Langford，DW，《用有限元技术模拟裂隙多孔岩体中的流体流动》，国际J·数值。液体方法，4，4，337-348（1984）·Zbl 0579.76095号 ·doi:10.1002/fld.165040404 [19] Reichenberger，V。；雅各布斯，H。；Bastian，P。；Helmig，R.，《裂隙多孔介质中两相流的混合维有限体积法》，水资源部。，2006年7月29日，1020-1036·doi:10.1016/j.advwatres.2005.09.001 [20] 李，L。；Lee，SH，通过离散裂缝网络和均质介质对天然裂缝性油藏中的黑油进行高效现场模拟，SPE Reserve。评估。工程，11，4，750-758（2008）·doi:10.2118/103901-PA [21] Fumagalli，A。；斯科蒂，A。；Cangiani，A。；Davidchack，RL；Georgoulis，E。；AN Gorban；Levesley，J。；Tretyakov，MV，《具有非匹配网格的裂隙多孔介质中流动和传输的简化模型》，《数值数学与高级应用》2011年第499-507页（2013年），柏林：施普林格出版社，柏林·兹比尔1273.76398 ·doi:10.1007/978-3642-33134-353 [22] 弗莱米什，B。；Fumagalli，A。；斯科蒂，A。；文丘拉，G。；Benvenuti，E.，《基于XFEM的裂隙多孔介质流动近似方法综述》，离散化方法进展，47-76（2016），Cham:Springer International Publishing，Cham·Zbl 1371.76093号 ·doi:10.1007/978-3-319-41246-73 [23] 北施温克。；弗莱米什，B。；赫尔米格，R。；Wohlmuth，BI，《裂缝性多孔介质中的降维流动模型：交叉点和边界》，计算。地质科学。，19, 6, 1219-1230 (2015) ·Zbl 1391.76747号 ·doi:10.1007/s10596-015-9536-1 [24] 姜杰。；Younis，RM，一种用于裂缝性储层多相流的改进的基于投影的嵌入式离散裂缝模型（pEDFM），高级水利资源。，109, 267-289 (2017) ·doi:10.1016/j.advwatres.2017.09.017 [25] 弗莱米什，B。；Darcis，M。；Erbertseder，K。；费格尔，B。；Lauser，A。；Mosthaf，K。；缪兴，S。；努斯克，P。；Tatomir，A。；沃尔夫，M。；Helmig，R.，DuMux:DUNE用于多孔介质中的多相、组分、尺度、物理等流动和传输，高级水资源。，34, 9, 1102-1112 (2011) ·文件编号：10.1016/j.advwatres.2011.03.007 [26] 斯洛伐克马特哈伊；盖革，S。；罗伯茨，SG；Paluszny，A。；Belayneh，M。；Burri，A。；Mezentsev，A。；卢，H。；库姆·D·。；Driesner，T。；Heinrich，CA，《结构复杂油藏中多相流体流动的数值模拟》，Geol。Soc.伦敦。规范出版物。，292, 1, 405-429 (2007) ·doi:10.144/SP292.22文件 [27] D.加斯顿。；纽曼，C。；Hansen，G。；Lebrun-Grandié，D.，MOOSE:非线性方程耦合系统的并行计算框架，Nucl。工程设计。，239, 10, 1768-1778 (2009) ·doi:10.1016/j.nucingdes.2009.05.021 [28] Breede，K。；Dzebisashvili，K。；刘，X。；Falcone，G.，《增强型（或工程化）地热系统的系统回顾：过去、现在和未来》，《地热》。能源，1,1,4（2013）·doi:10.1186/2195-9706-1-4 [29] Wang，W。；Kolditz，O.，《多孔介质中热流体力学（THM）问题的面向对象有限元分析》，国际期刊编号。方法工程，69，1，162-201（2007）·Zbl 1176.74199号 ·doi:10.1002/nme.1770 [30] Březina，J。；Stebel，J。；Kozubek，T。；布莱赫塔，R。；谢斯特克，J。；罗兹洛日尼克，M。；C ermák，M.，多孔介质流动连续断裂模型的模型误差分析，科学与工程中的高性能计算，152-160（2016），Cham:Springer International Publishing，Cham·Zbl 1382.76240号 ·doi:10.1007/978-3-319-40361-811 [31] Lie，K.-A.：使用MATLAB/GNU倍频程的油藏模拟简介：MATLAB油藏模拟工具箱（MRST）用户指南，第1版。剑桥大学出版社（2019）·Zbl 1425.76001号 [32] 谎言，K-A；克罗格斯塔德，S。；利加登，IS；Natvig，JR；HM尼尔森；Skaflestad，B.，复杂网格上一致离散化的开源MATLAB实现，计算。地质科学。，16, 2, 297-322 (2012) ·Zbl 1348.86002号 ·doi:10.1007/s10596-011-9244-4 [33] Alns，M.等人：FEniCS项目1.5版。架构（architecture）。数字。柔和。3, (2015). doi:10.11588/ans.2015.100.20553 [34] Blatt，M.等人：分布式统一数值环境，2.4版。架构（architecture）。数字。柔和。2016年4月4日。doi:10.11588/ans.2016.100.26526 [35] Rathgeber，F。；哈姆，DA；米切尔，L。；兰格，M。；卢波里尼，F。；ATT麦克雷；GT Bercea；马卡尔，GR；Kelly，PHJ，Firedrake：通过组合抽象实现有限元方法自动化，ACM Trans。数学。软质。，43, 3, 1-27 (2016) ·Zbl 1396.65144号 ·doi:10.1145/2998441 [36] Boon，W.M.，Nordbotten，J.M.，Vatne，J.E.：混合维几何上的函数分析和外部微积分。Ann.Mat.（2020年）。doi:10.1007/s10231-020-01013-1·Zbl 1462.58001号 [37] Boon，W.M.，Nordbotten，J.M.：含薄夹杂物线性弹性的稳定混合有限元。arXiv公司。1903.01757, (2019) [38] Nordbotten，JM；布恩，WM；Fumagalli，A。；Keilegavlen，E.，裂隙多孔介质中流动离散化的统一方法，计算。地质科学。，23, 2, 225-237 (2019) ·兹比尔1414.76064 ·doi:10.1007/s10596-018-9778-9 [39] Geuzaine，C。；Remacle，J-F，Gmsh：一个具有内置预处理和后处理设施的三维有限元网格生成器，国际期刊Numer。方法工程，79，11，1309-1331（2009）·Zbl 1176.74181号 ·doi:10.1002/nme.2579 [40] Karimi-Fard，M.：适用于通用油藏模拟器的高效离散裂缝模型。SPE J.9（2），（2004）。网址：10.2118/88812-PA [41] Hui，M.-H.，Mallison，B.，Lim，K.-T.：模拟巨型碳酸盐岩储层裂缝的创新工作流程。程序。国际宠物。技术会议15（2008） [42] Berre，I.等人：三维裂隙多孔介质中单相流动的验证基准。arXiv公司。2002.07005, (2020) [43] Boon，W.M.，Nordbotten，J.M.，Yotov，I.：裂隙多孔介质中流动的稳健离散化。SIAM J.数字。分析。56(4), 2203-2233 (2018). doi:10.1137/17M1139102·Zbl 1402.65149号 [44] Quarteroni，A。；Valli，A.，偏微分方程的数值逼近（1997），柏林：施普林格出版社，柏林 [45] 加里波夫，TT；Karimi-Fard，M。；Tchelepi，HA，耦合流和地质力学的离散断裂模型，计算。地质科学。，20, 1, 149-160 (2016) ·Zbl 1392.76079号 ·doi:10.1007/s10596-015-9554-z [46] 尤卡尔，E。；Keilegavlen，E。；贝尔，I。；Nordbotten，JM，断裂介质变形的有限体积离散化，计算。地质科学。，22, 4, 993-1007 (2018) ·Zbl 1396.74099号 ·doi:10.1007/s10596-018-9734-8 [47] 麦克卢尔，MW；霍恩，RN，《增强地热系统中刺激机制的研究》，国际岩石力学杂志。最小科学。，72, 242-260 (2014) ·doi:10.1016/j.ijrmms.2014.07.011 [48] Coussy，O.，《多孔力学》（2003），威利：奇切斯特，威利·数字对象标识代码：10.1002/0470092718 [49] Berge，RL；贝尔雷，I。；凯莱加夫伦，E。；Nordbotten，JM；Wohlmuth，B.，多孔弹性介质的有限体积离散化，含由接触力学建模的裂缝，国际期刊数值。方法工程，121，4，644-663（2020）·doi:10.1002/nme.6238 [50] Dong，S。；曾磊。；Dowd，P。；徐，C。；Cao，H.，离散断裂网络中断裂交点检测的快速方法，计算。岩土工程。，98, 205-216 (2018) ·doi:10.1016/j.compgeo.2018.02.005 [51] Mallison，B.T.，Hui，M.H.，Narr，W.：离散裂缝建模工作流的实用网格算法。在英国牛津举行的第十二届欧洲石油开采数学会议上提出（2010年）。doi:10.3997/2214-4609.20144950 [52] 霍尔姆，R。；考夫曼，R。；海姆松，B-O；伊恩·E·。；Espedal，MS，复杂内部几何域的网格划分，数值。线性代数应用。，13, 9, 717-731 (2006) ·Zbl 1174.76363号 ·doi:10.1002/nla.505 [53] Berge，RL；克莱梅茨达尔；Lie，K-A，符合低维对象的非结构化Voronoi网格，计算。地质科学。，23, 1, 169-188 (2019) ·Zbl 1411.65155号 ·doi:10.1007/s10596-018-9790-0 [54] Shewchuk，JR，《三角形：设计2D质量的网格生成器和Delaunay三角测量器》，《应用计算几何：走向几何工程》，203-222（1996）·doi:10.1007/BFb0014497 [55] Si，H.，TetGen，基于Delaunay的优质四面体网格生成器，ACM Trans。数学。软质。，41, 2, 1-36 (2015) ·Zbl 1369.65157号 ·数字对象标识代码：10.1145/2629697 [56] Boffi，D。；布雷齐，F。；Fortin，M.，《混合有限元方法和应用》（2013），柏林：施普林格出版社，柏林·Zbl 1277.65092号 ·doi:10.1007/978-3-642-36519-5 [57] 达维加，LB；布雷齐，F。；马里尼，LD；Russo，A.，多边形网格上一般二阶椭圆问题的混合虚元方法，ESAIM Math。模型。数字。分析。，50, 3, 727-747 (2016) ·Zbl 1343.65134号 ·doi:10.1051/m2安/2015067 [58] 达维加，LB；布雷齐，F。；马里尼，LD；Russo，A.，H（div）和H（curl）协调虚拟元方法，数值。数学。，133, 2, 303-332 (2016) ·Zbl 1343.65133号 ·doi:10.1007/s00211-015-0746-1 [59] Nordbotten，JM，线性弹性以单元为中心的有限体积离散的收敛性，SIAM J.Numer。分析。，53, 6, 2605-2625 (2015) ·Zbl 1330.74171号 ·数字对象标识代码：10.1137/140972792 [60] Keilegavlen，E。；Nordbotten，JM，《弱对称弹性力学的有限体积法》，国际期刊数值。方法工程，112，8，939-962（2017）·doi:10.1002/nme.5538 [61] Nordbotten，JM，《Biot方程的稳定单元中心有限体积离散化》，SIAM J.Numer。分析。，54,2942-968（2016）·Zbl 1382.76187号 ·doi:10.1137/15M1014280 [62] Nordbotten，J.M.，Keilegavlen，E.：热弹性多点通量（MPFA）和应力（MPSA）有限体积法简介。arXiv公司。2001.01990, (2020) [63] Hüeber，S。；斯塔德勒，G。；Wohlmuth，BI，《库仑摩擦三维接触问题的原对偶主动集算法》，SIAM J.Sci。计算。，30, 2, 572-596 (2008) ·Zbl 1158.74045号 ·doi:10.1137/060671061 [64] Stefansson，I.、Berre，I.和Keilegavlen，E.：多孔介质中热流体力学过程和断裂接触力学的完全耦合数值模型。arXiv:2008.06289，（2020）·Zbl 1454.65152号 [65] Berge，R.L.，Berre，I.，Keilegavlen，E.，Nordbotten，J.M.：裂缝性多孔介质中的粘性指进。arXiv:1906.10472，（2019） [66] Budisa，A.，Boon，W.，Hu，X.：混合维辅助空间预条件子。arXiv:1910.04704，（2019）·Zbl 1452.65049号 [67] Budiša，a.，Hu，X.：裂缝性多孔介质中混合多维离散流动的块预条件。计算。地质科学。（2020年）。doi:10.1007/s10596-020-09984-z·Zbl 1460.76634号 [68] Ahrens，J.，Geveci，B.，Law，C.：ParaView：大型数据可视化的最终用户工具 [69] Fumagalli，A。；Keilegavlen，E。；Scialó，S.，《离散裂缝网络模拟中的一致、非一致和非匹配离散耦合》，J.Compute。物理。，376694-712（2019）·Zbl 1416.76306号 ·doi:10.1016/j.jcp.2018.09.048 [70] Fumagalli，A。；Keilegavlen，E.，离散裂缝矩阵模型的双重虚拟元方法，石油天然气科学。Technol.-技术新能源D’IFP版次。，74, 41 (2019) ·doi:10.2516/ogst/2019008 [71] 斯蒂芬森，I。；贝尔雷，I。；Keilegavlen，E.，裂隙多孔介质中流动的有限体积离散化，Transp。多孔介质，124,2439-462（2018）·doi:10.1007/s11242-018-1077-3 [72] 使用runscripts.doi:10.5281/zenodo.3374624实现PorePy。(2019) [73] 弗莱米什，B。；贝尔雷，I。；布恩，W。；Fumagalli，A。；北施温克。；斯科蒂，A。；斯蒂芬森，I。；Tatomir，A.，《裂缝性多孔介质中单相流动的基准》，Adv.Water Resour。，111, 239-258 (2018) ·doi:10.1016/j.advwatres.2017.10.036 [74] Mandel，J.，《土壤综合》（Consolidation des sols），《岩土工程》，第3期，第7期，第287-299页（1953年）·doi:10.1680/geot.1953.3.7.287 [75] 阿布斯莱曼，Y。；程，A-D；崔，L。；Detournay，E。；罗杰斯，J-C，《重新审视曼德尔问题》，《岩土工程》，46，2，187-195（1996）·doi:10.1680/geot.1996.46.2187 [76] 郑，AH-D；Detournay，E.，平面应变多孔弹性的直接边界元法，国际J数值。分析。方法地质力学。，12, 5, 551-572 (1988) ·Zbl 0662.73056号 ·doi:10.1002/nag.1610120508 [77] Mikelić，A。；王，B。；惠勒，MF，耦合流和地质力学迭代耦合的数值收敛研究，计算。地质科学。，18, 3-4, 325-341 (2014) ·Zbl 1386.76115号 ·doi:10.1007/s10596-013-9393-8 [78] 斯奈登，印第安纳州，傅里叶变换（1995），纽约：多佛出版社，纽约 [79] 克劳奇，SL；Starfield，AM，《固体力学中的边界元方法：在岩石力学和地质工程中的应用》（1983），伦敦：Allen&Unwin出版社，伦敦·Zbl 0528.73083号 ·数字对象标识代码：10.1115/1.3167130 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。