×
杰弗里·海曼。；马修·斯威尼。；卡尔·W·盖博。；丹尼尔·斯维亚茨基；康斯坦丁·利普尼科夫；J.David Moulton

在符合多维网格上使用模拟有限差分的三维离散裂缝矩阵模型中的流动和传输。 （英语） Zbl 07561073号

J.计算。物理学。 466，文章ID 111396，19 p.（2022）.
总结：我们提出了一个综合的工作流程，用离散裂缝矩阵方法模拟裂隙多孔介质中的单相流动和传输。该工作流程有三个主要部分：（1）三维断裂网络及其周围的协调网格生成方法，（2）使用二阶模拟有限差分方法离散控制方程，以及（3）高性能计算环境中数值方法的实现。本发明还提供了一种方法，用于创建断裂网络周围体积的一致Delaunay四面体化，其中断裂网格的三角形单元是体积网格中的面，该方法解决了常见的病理情况，并且降低了网格质量。我们的开源地下模拟器使用了一个过程内核层次结构（每个物理过程一个内核），允许断裂域和矩阵域的强耦合和弱耦合。我们基于流动和传输的解析解以及数值收敛性提供验证测试。我们还对该方法在复杂裂隙网络中的应用进行了多次阐述。在第一个示例中，我们通过考虑断裂网络作为流动屏障、主要通道或提供与周围基质相同阻力的两种情况，证明了该方法的鲁棒性。在第二个测试中，呈现了通过包含257个裂缝的三维随机生成网络的流动和传输。

引用于三文件

MSC公司：

65新元 偏微分方程边值问题的数值方法
76平方米 多孔介质中的流动；过滤；渗流
7.6亿 流体力学基本方法

关键词：

离散断裂矩阵；计算几何；模拟有限差分；高性能计算；裂隙介质中的流动与运移；Delaunay三角测量

软件：

拉格里特；dfn工程；PFLOTRAN公司
PDF格式BibTeX公司 XML格式引用
\textit{J.D.Hyman}等人，J.Compute。物理学。466，文章ID 111396，第19页（2022；Zbl 07561073）
全文： DOI程序 arXiv公司

参考文献：

[1] 国家研究委员会，《岩石断裂与流体流动：当代理解与应用》（1996），国家科学院出版社
[2] Neuman，S.，《通过断裂岩石量化流量和运输的趋势、前景和挑战》，《水文地质》。J.，13，1，124-147（2005）
[3] 范德夸克，J。；Sudicky，E.，《非水相液体的溶解和离散无裂缝多孔介质中水相污染物的迁移》，J.Contam。水文。，23，1-2，45-68（1996）
[4] Kueper，B.H。；McWhorter，D.B.，《断裂粘土和岩石中致密非水相液体的行为》，《地下水》，29，5，716-728（1991）
[5] 海曼，J。；吉姆·马丁内斯，J。；Viswanathan，H。；凯里·J。；波特，M。；Rougier，E。；Karra，S。；康（Kang，Q.）。；弗拉什，L。；Chen，L.，《理解水力压裂：多尺度问题》，Philos。事务处理。R.Soc.A，374，2078，第20150426条，pp.（2016）
[6] 米德尔顿，R。；凯里·J。；柯里埃，R。；海曼，J。；康（Kang，Q.）。；Karra，S。；吉姆·马丁内斯，J。；波特，M。；Viswanathan，H.，《页岩气和非水压裂液：超临界CO_2的机遇和挑战》，应用。能源，147500-509（2015）
[7] Hadgu，T。；Karra，S。；加里尼亚，E。；Makedonska，N。；海曼·J·D。；Klise，K。；Viswanathan，H.S。；Wang，Y.，离散裂缝网络和等效连续体模型的比较研究，用于模拟晶体宿主岩中假想核废料库远场的流动和传输，J.Hydrol。，553, 59-70 (2017)
[8] 乔伊斯，S。；哈特利，L。；Applegate，D。；霍克，J。；Jackson，P.，《温带气候条件下的多尺度地下水流动模拟，用于瑞典福斯马克拟建高级核废料处置场的安全评估》，Hydrogeol。J.，22，61233-1249（2014）
[9] 海曼·J·D。；吉梅内斯·马丁内斯，J。；盖博，C.W。；施陶夫，P.H。；Pawar，R.J.，《表征断裂盖层非均质性对超临界CO_2注入的影响》，Trans。多孔介质，1-21（2019）
[10] 詹金斯，C。；查德威克，A。；Hovorka，S.D.，《十年来监测和核查的最新进展》，国际J.Greenh出版社。天然气控制，40，312-349（2015）
[11] 贝尔雷，I。；加药器，F。；Keilegavlen，E.，《裂隙多孔介质中的流动：概念模型和离散方法综述》，Transp。多孔介质（2018年10月）
[12] Neuman，S.P。；Dener，J.S.，《使用可变尺度压力试验数据估算Oracle、亚利桑那州J.Hydrol附近断裂花岗岩的对数导水率协方差和分散性》。，102, 1-4, 475-501 (1988)
[13] Tsang，Y。；曾荫权（Tsang，C.）。；Hale，F.公司。；Dverstorp，B.，断裂介质随机连续介质模型中的示踪输运，《水资源》。决议，32，10，3077-3092（1996）
[14] Gerke，H。；Van Genuchten，M.T.，用于模拟结构性多孔介质中水和溶质优先运动的双重孔隙模型，water Resour。研究，29，2，305-319（1993）
[15] Lichtner，P。；Karra，S.，《多孔介质中多尺度多相多组分反应流建模：利用PFLOTRAN应用于CO_2封存和增强地热能》，（Al-Khoury，R.；Bundschuh，J.，《CO_2地质封存和压缩空气储能计算模型》（2014），CRC出版社），81-136
[16] 齐默尔曼，R.W。；陈，G。；哈古，T。；Bodvarsson，G.S.，《裂缝/基质流半解析处理的数值双重孔隙模型》，《水资源》。Res.，292127-2137（1993年）
[17] 卡卡斯，M.C。；Ledoux，E。；De Marsily，G。；Barbreau，A。；Calmels，P。；盖拉德，B。；Margritta，R.，用随机离散裂缝网络建模裂缝流：校准和验证：2。交通模型，水资源。研究，26，3，491-500（1990）
[18] Long，J。；雷默，J。；Wilson，C。；Witherspoon，P.，《不连续裂缝网络的多孔介质当量》，《水资源》。研究，18，3，645-658（1982）
[19] 梅洛特，J。；戴维·P。；勒戈克，R。；Darcel，C。；De Dreuzy，J.-R.，《随机分布和运动定义的离散裂缝网络模型中的连通性、渗透率和渠化》，《水资源》。研究，52，11，8526-8545（2016）
[20] Makedonska，N。；海曼，J.D.D。；Karra，S。；Painter，S.L。；盖博，C.W.W。；Viswanathan，H.S.，《评估内部孔径变化对千米级离散裂隙网络传输的影响》，水资源高级研究所。，94486-497（2016）
[21] Nordqvist，A.W。；Tsang，Y.W。；Tsang，C.F。；Dverstorp，B。；Andersson，J.，《裂隙岩石中流动和传输的可变孔径裂隙网络模型》，《水资源》。1703-1713年第28号、第6号决议（1992年）
[22] 艾哈迈德·R。；爱德华兹，M.G。；拉明，S。；Huisman，B.A。；Pal，M.，控制体积分布多点通量近似与低维断裂模型耦合，J.Compute。物理。，284462-489（2015年）·Zbl 1351.74073号
[23] 艾哈迈德·R。；爱德华兹，M.G。；拉明，S。；Huisman，B.A。；Pal，M.，三维控制体积分布多点通量近似与低维表面断裂模型耦合，J.Compute。物理。，303, 470-497 (2015) ·兹比尔1349.74343
[24] Antonietti，P.F。；Formaggia，L.公司。；斯科蒂，A。；维拉尼，M。；Verzott，N.，《裂隙多孔介质中流动的模拟有限差分近似》，ESAIM:Math。模型。数字。分析。，50, 3, 809-832 (2016) ·Zbl 1381.76231号
[25] Arraras，A。；Gaspar，F.J。；Portero，L。；Rodrigo，C.，裂隙多孔介质单相流动的混合多维几何多重网格方法，SIAM J.Sci。计算。，41,5，B1082-B1114（2019）·Zbl 1423.76327号
[26] 布恩，W.M。；Nordbotten，J.M。；Yotov，I.，裂缝性多孔介质中流动的稳健离散化，SIAM J.Numer。分析。，56, 4, 2203-2233 (2018) ·Zbl 1402.65149号
[27] Berrone，S。；A.鲍里奥。；Fidelibus，C。；Pieraccini，S。；Scialo，S。；Vicini，F.，离散裂缝矩阵模型中稳态流体流动的高级计算：fem-bem和vem-vem裂缝-块体耦合，GEM Int.J.Geomath。，9, 2, 377-399 (2018) ·Zbl 1419.65103号
[28] 弗莱米什，B。；Fumagalli，A。；Scotti，A.，《基于XFEM的裂缝性多孔介质中流动近似研究综述》，47-76（2016），Springer International Publishing：Springer国际出版公司Cham·Zbl 1371.76093号
[29] Fumagalli，A。；Scotti，A.，《具有非匹配网格的裂隙多孔介质中流动和传输的简化模型》（Cangiani，A.；Davidchack，R.L.；Georgoulis，E.；Gorban，A.N.；Levesley，J.；Tretyakov，M.V.，《数值数学与高级应用》，2011卷（2013），Springer:Springer Berlin，Heidelberg），499-507·Zbl 1273.76398号
[30] Odster，L.H。；Kvamsdal，T。；Larson，M.G.，多孔介质中流动和输运耦合问题的简单嵌入离散裂缝矩阵模型，计算机。方法应用。机械。工程，343，572-601（2019）·Zbl 1440.76147号
[31] 科佩尔，M。；马丁·V。；贾夫雷，J。；Roberts，J.E.，多孔介质中流动离散裂缝模型的拉格朗日乘子法，计算。地质科学。，23, 2, 239-253 (2019) ·Zbl 1414.76059号
[32] 曼祖尔，S。；爱德华兹，M.G。；Dogru，A.H。；Al-Shaalan，T.M.，《内部边界对齐非结构化网格生成和以单元为中心与以顶点为中心的CVD-MPFA性能》，计算。地质科学。，22, 1, 195-230 (2018) ·Zbl 1405.65160号
[33] Sandwe，T.H。；贝尔雷，I。；Nordbotten，J.M.，离散断裂矩阵模拟的有效多点通量近似方法，J.Compute。物理。，231, 9, 3784-3800 (2012) ·Zbl 1402.76131号
[34] 北施温克。；弗莱米什，B。；赫尔米格，R。；Wohlmuth，B.I.，《裂缝性多孔介质中的降维流动模型：交叉点和边界》，计算。地质科学。，19, 6, 1219-1230 (2015) ·Zbl 1391.76747号
[35] Svensk Kärnbränslehantering AB，安全评估数据报告SR-site（TR-10-52）。
[36] 哈特利，L。；Joyce，S.，《支持forsmark现场调查和安全评估的地下水流量建模方法和算法》，瑞典，J.Hydrol。，500, 200-216 (2013)
[37] C.P.杰克逊。；Hoch，A.R。；Todman，S.，非均匀连续多孔介质的自洽性，断裂介质的表征，水资源。研究，36，1，189-202（2000）
[38] 斯威尼，M.R。；盖博，C.W。；Karra，S。；Stauffer，P.H。；帕瓦尔·R·J。；Hyman，J.D.，《升级离散裂缝矩阵模型（UDFM）：裂缝多孔介质的八重定义连续表示法》，计算。地质科学。，1-18 (2019)
[39] Berrone，S。；Pieraccini，S。；Scialo，S.，离散裂缝网络流动的PDE约束优化公式，SIAM J.Sci。计算。，35、2、B487-B510（2013）·Zbl 1266.65188号
[40] 贝隆，S。；Pieraccini，S。；Scialó，S。；Vicini，F.，《大规模DFN流动模拟的并行求解器》，SIAM J.Sci。计算。，37、3、C285-C306（2015）·Zbl 1320.65167号
[41] 戴维·P。；勒高克（Le Goc，R.）。；Darcel，C.，《裂缝成核、增长和停止的模型以及裂缝密度和缩放的后果》，J.Geophys。固体地球研究，118，4，1393-1407（2013）
[42] 戴维·P。；勒戈克，R。；Darcel，C。；布尔·O。；de Dreuzy，J.R。；Munier，R.，《断裂尺度的可能通用模型及其对地壳流体力学的影响》，J.Geophys。固体地球研究，115，B10（2010）
[43] 德德鲁齐，J.-R。；Darcel，C。；戴维·P。；Bour，O.，根据幂律长度分布，裂缝中心的空间相关性对二维裂缝网络渗透率的影响，水资源。决议，40，1（2004）
[44] 德肖维茨，W。；Fidelibus，C.，通过边界元法推导三维离散裂隙网络的等效管网模拟，Water Resour。研究，35，9，2685-2691（1999）
[45] 德德鲁齐，J.-R。；梅赫斯特，Y。；Pichot，G.，裂缝尺度非均质性对三维离散裂缝网络流动特性的影响，地球物理杂志。固体地球研究，117，B11（2012）
[46] 埃赫尔，J。；德德鲁齐，J.-R。；Poirriez，B.，三维离散裂缝网络中的流动模拟，SIAM J.Sci。计算。，31, 4, 2688-2705 (2009) ·Zbl 1387.65124号
[47] 海曼·J·D。；盖博，C.W。；Painter，S.L。；Makedonska，N.，《随机生成的三维离散断裂网络的一致Delaunay三角剖分：网格化策略的特征抑制算法》，SIAM J.Sci。计算。，第36页，第4页，第1871-A1894页（2014年）·Zbl 1305.74082号
[48] 海曼·J·D。；Karra，S。；Makedonska，N。；盖博，C.W。；Painter，S.L。；Viswanathan，H.S.，DFN-works:模拟地下水流和输移的离散裂缝网络框架，计算。地质科学。，2015年10月19日，第84页
[49] 皮肖特，G。；埃赫尔，J。；de Dreuzy，J.-R.，求解随机离散裂隙网络中流动的广义混合砂浆方法，SIAM J.Sci。计算。，34、1、B86-B105（2012）·Zbl 1387.65121号
[50] 穆斯塔法，H。；Mustapha，K.，用优化网格模拟离散裂缝网络中流动的新方法，SIAM J.Sci。计算。，29, 1439 (2007) ·Zbl 1251.76056号
[51] 阀盖，E。；布尔·O。；Odling，东北部。；戴维·P。；主要，I。；Cowie，P。；Berkowitz，B.，地质介质中断裂系统的缩放，《地球物理学评论》。，39, 3, 347-383 (2001)
[52] 艾哈迈德·R。；爱德华兹，M.G。；拉明，S。；Huisman，文学学士。；Pal，M.，CVD-MPFA全压支持，耦合非结构离散断裂矩阵Darcy-flux近似，J.Compute。物理。，349, 265-299 (2017) ·Zbl 1380.76092号
[53] 博格达诺夫，I。；V.穆尔琴科。；Thovert，J.-F。；Adler，P.，《穿过裂缝性多孔介质的两相流》，Phys。E版，68，2，第026703条，pp.（2003）
[54] Helmig，R.，《地下多相流和输运过程：对水文系统建模的贡献》（1997），Springer-Verlag
[55] Karimi-Fard，M。；杜洛夫斯基，L.J。；Aziz，K.，适用于通用储层模拟器的高效离散裂缝模型，（SPE储层模拟研讨会（2003），石油工程师学会）
[56] Berrone，S。；Pieraccini，S。；Scialó，S.，《浸没交叉裂缝多孔介质中的流动模拟》，J.Compute。物理。，345, 768-791 (2017) ·Zbl 1378.76107号
[57] 北卡罗来纳州弗里。；马丁·V。；罗伯茨，J.E。；Saáda，A.，《将断裂建模为与非匹配网格的接口》，计算。地质科学。，16, 4, 1043-1060 (2012)
[58] 弗莱米什，B。；贝尔雷，I。；布恩，W。；Fumagalli，A。；北施温克。；斯科蒂，A。；斯蒂芬森，I。；Tatomir，A.，《裂缝性多孔介质中单相流动的基准》，Adv.Water Resour。，111, 239-258 (2018)
[59] 德尔普拉，M。；Fumagalli，A。；Scotti，A.，《完全耦合裂缝/整体达西流与XFEM的适定性》，SIAM J.Numer。分析。，55, 2, 785-811 (2017) ·Zbl 1362.76057号
[60] Nordbotten，J.M。；布恩，W.M。；Fumagalli，A。；Keilegavlen，E.，裂隙多孔介质中流动离散化的统一方法，计算。地质科学。，23, 2, 225-237 (2019) ·兹比尔1414.76064
[61] Abushaikha，A.S。；Terekhov，K.M.，《全张量渗透率通用地下油藏模拟的全隐式模拟有限差分格式》，J.Compute。物理。，406，第109194条pp.（2020）·Zbl 1453.76120号
[62] Lipnikov，K。；Manzini，G。；Shashkov，M.，《模拟有限差分法》，J.Compute。物理。，257, 1163-1227 (2014) ·兹比尔1352.65420
[63] 马丁·V。；贾夫雷，J。；Roberts，J.E.，《将裂缝和屏障建模为多孔介质中流动的界面》，SIAM J.Sci。计算。，26, 1667-1691 (2005) ·Zbl 1083.76058号
[64] Lipnikov，D.S.Konstantin；Shashkov，Mikhail，非结构化多面体网格上扩散问题的模拟有限差分离散化，J.Compute。物理。，2, 211, 473-491 (2006) ·Zbl 1120.65332号
[65] Barth，T.，非结构化网格的三维迎风Euler解算器，（AIAA第十届计算流体动力学会议，AIAA第10届计算流体力学会议，夏威夷州火奴鲁鲁，1991年6月24-27日（1991）），24-27
[66] 莫尔顿，J.D。；梅扎，J.C。；Buksas，M.W。；Day，M.，《多进程高性能计算模拟器Amanzi的高级设计》（2012年），美国能源部环境管理办公室：美国能源部华盛顿特区环境管理办公室，技术代表ASCEM-HPC-2011-03-1（LANL:LA-UR 12-22193）
[67] 库恩，E.T。；莫尔顿，J.D。；Painter，S.L.，《管理地表和近地表过程模拟的复杂性》，环境。模型。软质。，78, 134-149 (2016)
[68] Lipnikov，K。；莫尔顿，D。；Svyatskiy，D.，用Richards方程求解可变饱和流的无雅可比解算器的新预处理策略，高级水资源。，94, 11-22 (2016)
[69] 贝尔雷，I。；布恩，W.M。；弗莱米什，B。；Fumagalli，A。；Gläser博士。；Keilegavlen，E。；斯科蒂，A。；斯蒂芬森，I。；Tatomir，A。；Brenner，K.，《三维裂隙多孔介质中单相流动的验证基准》，Adv.Water Resour。，147，第103759条pp.（2021）
[70] 墨菲，M。；Mount，D.M。；Gable，C.W.，《协调Delaunay四面体化的点塑性策略》，国际计算杂志。地理。申请。，11, 06, 669-682 (2001) ·Zbl 1074.68646号
[71] Khamayseh，A。；Kuprat，A.，《非结构化网格的各向异性平滑和解决方案自适应》，国际期刊Numer。方法工程，39，18，3163-3174（1996）·Zbl 0873.73080号
[72] R.E.银行。；Smith，R.K.，使用后验误差估计进行网格平滑，SIAM J.Numer。分析。，34, 3, 979-997 (1997) ·兹伯利0873.65092
[73] LaGriT，洛斯阿拉莫斯网格工具箱（2013），（LaGriT）洛斯阿拉莫斯国家实验室
[74] Kuprat，A.，三维非结构化网格的自适应平滑技术（4 1996）
[75] 郭杰。；严博士。；Chen，L。；张，X。；O.Deussen。；Wonka，P.，通过最大泊松圆盘采样的四面体网格，计算。辅助Geom。设计。，43186-199（2016）·Zbl 1418.65189号
[76] Krotz，J。；斯威尼，M.R。；盖博，C.W。；海曼，J.D。；Restrepo，J.M.，网格离散断裂网络的可变分辨率泊松圆盘采样，J.Compute。申请。数学。，407，第114094条pp.（2022）·Zbl 1493.65256号
[77] Abelin，H。；Birgersson，L.公司。；莫雷诺，L。；维登，H。；Ågren，T.等人。；Neretnieks，I.，花岗岩中的大规模流动和示踪实验：2。结果和解释，水资源。研究，27，12，3119-3135（1991）
[78] Abelin，H。；内雷特尼克斯，I。；Tunbrant，S.公司。；Moreno，L.，《单个断裂中迁移的最终报告：实验结果和评估》，Nat.Genossenschaft fd Lagerung Radioaktiver-Abfälle（1985）
[79] Hyman，J.D.，《断裂网络中的流动通道：表征密度对优先流动路径形成的影响》，《水资源》。决议，56，9，条款e2020 WR027986 pp.（2020）
[80] 拉斯穆森，A。；Neretnieks，I.，《放射性核素在晶体岩石快速通道中的传输》，《水资源》。第22、8、1247-1256号决议（1986年）
[81] 海曼·J·D。；Dentz，M。；哈格伯格，A。；Kang，P.，《三维断裂网络中连接结构和传输特性》，J.Geophys。固体地球研究（2019年）
[82] 谢尔曼，T。；海曼·J·D。；Dentz，M。；Bolster，D.，《表征断裂密度对网络规模传输的影响》，J.Geophys。固体地球研究，125，1（2020）
[83] 德德鲁齐，J.-R。；戴维·P。；Bour，O.，长度和孔径幂律分布下二维随机裂隙网络的水力特性，水资源。决议，38，12（2002），12-1
[84] 海曼·J·D。；奥尔德里奇，G。；Viswanathan，H。；Makedonska，N。；Karra，S.，《裂缝尺寸和透射率相关性：根据截断的裂缝尺寸幂律分布对稀疏三维离散裂缝网络中输运模拟的影响》，水资源。决议，52，8，6472-6489（2016）
[85] 研磨棒，P。；Impey，M.，分形模拟多孔介质中的沟道和Fickian色散，《水资源》。研究，29，12，4077-4089（1993）
[86] 奥尔德里奇，G。；海曼·J·D。；Karra，S。；盖博，C.W。；Makedonska，N。；Viswanathan，H。；伍德林，J。；Hamann，B.，《使用流动拓扑图分析和可视化离散断裂网络》，IEEE Trans。视觉。计算。图表。，23, 8, 1896-1909 (2017)
[87] 海曼·J·D。；哈格伯格，A。；斯里尼瓦桑，G。；Mohd-Yusof，J。；Viswanathan，H.，使用基于图形的约简预测稀疏离散断裂网络的首次通过时间，Phys。E版，96，第013304条pp.（2017）
[88] 海曼·J·D。；哈格伯格，A。；Osthus，D。；Srinivasan，S。；Viswanathan，H。；Srinivasan，G.，《识别三维离散断裂网络中的主干：基于二分图的方法，多尺度模型》。同时。，16, 4, 1948-1968 (2018) ·Zbl 1457.76165号
[89] 威尔曼，T.P。；夏皮罗，A.M。；Hill，M.C.，《裂缝控制含水层中简化裂缝网络表示对惰性化学迁移的影响》，《水资源》。决议，45，1（2009年）
[90] Wood，A.T.，von Mises Fisher分布模拟，Commun。统计模拟。，23, 1, 157-164 (1994) ·Zbl 0825.62022号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。