软件X 12（2020）100533原始软件出版物SHEMAT-Suite：一个用于模拟多孔介质Johannes Kellera，2004，Volker Ratha，1，Johanna Bruckmanna，Darius Mottaghya，2，作者：Christoph Clausera，Andreas Wolfb，Ralf Seidlerc，H. Martin Bückerc，d，Norbert Klitzschaa德国亚琛工业大学应用地球物理和地球物理能源研究所，52074 AachenbFachgebiet Scientific Computing，Technische Universität Darmstadt，64283 Darmstadt，Germanyc德国耶拿弗里德里希·席勒大学计算机科学研究所，07737耶拿dMichael Stifel Center Jena for Data Driven and Simulation Science，07737 Jena，Germanyar t i cl e i nf o文章历史记录：收到2020年收到修订版2020年5月20日接受2020年MSC：35L2065M0665M3268N1976R0576R1076R5086-0486-0886A05保留字：流多孔介质中的热量和物质输运GeologicalEnergy参数估计不确定性量化a b st ra ctSHEMAT-Suite是一个有限差分开源代码，用于模拟多孔介质中的该程序是用Fortran-95语言编写的，来源于地热和水文地质领域的地球科学研究它包括：（1）输入和输出的多功能处理，（2）地下参数建模的模块化框架，（3）多级OpenMP并行化，（4）通过随机方法（蒙特卡罗，EnergyKalman滤波器）和基于自动微分的确定性贝叶斯方法进行参数估计和数据同化，用于计算前向码的精确（无截断误差）导数©2020作者由爱思唯尔公司出版这是CC BY许可下的开放获取文章（http://creativecommons.org/licenses/by/4.0/）中找到。*通讯作者。电子邮件地址：jkeller@eonerc.rwth-aachen.dewww.example.com Keller）。1 现在在：都柏林高等研究所，爱尔兰都柏林2 现在：能源技术学院，FH亚琛应用科学大学，52428 Jülich，德国。https://doi.org/10.1016/j.softx.2020.1005332352-7110/©2020作者。 由Elsevier B.V.出版。这是一篇开放获取的文章，使用CC BY许可证（http://creativecommons.org/licenses/by/4.0/）。可在ScienceDirect上获得目录列表SoftwareX期刊主页：www.elsevier.com/locate/softx2J. Keller，V.Rath，J.Bruckmann等人粤公网安备44010502000011号当前代码版本当前代码版本v9.00用于此代码版本的代码/存储库的永久链接https://github.com/ElsevierSoftwareX/SOFTX_2020_135Code Ocean compute capsule临时doi：https://doi.org/10.24433/CO.9420115.v1法律代码许可证MIT许可证使用git的代码版本控制系统使用Fortran 90/95、OpenMP、CMake的软件代码语言、工具和服务编译 要求， 操作 环境依赖性Linux OS，GNU Make，CMake，Util-linux，GNU编译器选择，BLAS，LAPACK，HDF5，Doxygen，Python，Tapenade如果可用，链接到开发人员文档/手册https://git.rwth-aachen.de/SHEMAT-Suite/SHEMAT-Suite-open/-/wikis/home问题支持电子邮件shemat-suite-support@eonerc.rwth-aachen.de1. 动机和意义耦合的地下流体流动和热传输问题的数值模拟通常在一些地热能、地下水、核废料储存和二氧化碳封存等。模拟越来越多地基于大型和复杂的地质模型，需要复杂和并行的软件包，理想情况下，具有广泛的基准测试和实际应用测试的历史。SHEMAT-Suite就是这样一个软件：一个多功能的，通用的开源代码，用于模 拟 地 质 油 藏 中 的 流 动 ， 热 量 和 物 质 输 运 。 从 概 念 上 讲 ，SHEMAT-Suite将地质储层视为多孔介质。该代码解决了瞬态或稳态，正向和反向耦合的问题，在1D，2D和3D。它可以处理广泛的时间尺度，因此，可以解决技术和地质过程。SHEMAT-Suite可以处理岩石和流体性质随温度、压力和组分浓度变化而产生的非线性。此外，SHEMAT-Suite还能够进行确定性和随机性逆模拟和数据模拟，从而实现参数估计和不确定性量化。这些都是重要的特征，因为地下水流和输运建模受制于数据稀缺和不同来源的巨大不确定性。其他软件包也可以处理SHEMAT-Suite可以解决的问题。有些是作为购买研究代码或商业软件包提供的，如TOUGH [1]，FEFLOW [2]和STOMP [3]，其他则作为开源计划分发，包括OpenGeoSys [4]，DUMUX [5]和MOD- FLOW [6]。SHEMAT-Suite通过强大的功能工具包补充了可用的开源软件包，包括：(a)多孔介质流动模拟的多级OpenMP并行化[7-(c)包括密度驱动流在内的物种迁移;（d）通过随机方法（MonteCarlo、集合卡尔曼滤波器[10]）以及基于自动微分（AD）的确定性贝叶斯反演进行不确定性量化、参数估计和数据同化，以计算前向模型的精确雅可比矩阵[11]。SHEMAT-Suite目前的形式是从SHEMAT（热量和质量输运模拟器[12，13]）演变而来的，SHEMAT是一种完全耦合的流动和热量输运模型，由Trescott和Larson的等温USGS三维地下水模型[14，15]开发。SHEMAT-Suite源自Fortran-77代码SHEMAT，已在Fortran-95中重写[11]，随后在众多研究和博士项目中使用各种附加功能进行重构和增强。这些能力小型模拟，用于大型模拟的双共轭梯度稳定（Bi-CGStab）方法的实现[18];(4)由于其求解器实现的模块化结构，易于在SHEMAT-Suite中集成的附加求解器和预处理器;（5）从库BLAS [19-HDF 53;（7）确定性贝叶斯反演[11];（8）自动通过软件Tapenade [24]对前向码进行微分，以计算确定性反演框架中的雅可比矩阵;（9）从geo-统计库GSLib [25]，用于蒙特卡罗框架内随机模拟模式(10)贝叶斯更新作为集合卡尔曼滤波器（EnKF）分析方程[26]。5SHEMAT-Suite的当前开源版本仅限于多孔介质然而，多相流体版本，其可以，例如，模拟CO2、蒸汽[27，28]和超临界水[29]正在开发中。重要的应用领域，如地下水管理和污染物输运模拟，反映了SHEMAT-Suite迄今为止的科学贡献，可分为四类：地热，古气候和水文地质学的地球科学，以及，此外，逆方法的发展。地热是最重要的一类地质，也是软件开发的最初动力。SHEMAT-Suite在地热方面的应用范围从大规模地热模拟[30-第二个应用领域涉及古气候相关研究，包括用于模拟多孔介质中冻结和解冻过程的冰模块（例如[34，35]）。此外，SHEMAT-Suite还被用于广泛的地质研究，包括自上一个冰河时代以来美国新泽西州大陆架海底地下水排放的模拟[36]和澳大利亚珀斯盆地对流单元的模拟[37]。最后，使用SHEMAT-Suite进行了算法的一般开发。实施并比较了各种EnKF方法[10]。EnKF是一种连续的、基于集合的数据同化算法，可用于大规模模拟中的SHEMAT-Suite使用与高性能计算不同的技术[38，39]定向的服从中心实验设计框架EFCOSS [40]，以增强储层性质估计[41]。包括（1）SHEMAT-Suite [7]的并行化方案;（2）在有限差分离散化[16，17]中通过隐式Picard迭代求解耦合偏微分方程组，SHEMAT-Suite中的离散化在具有可变单元尺寸的矩形网格上实现;（3）直接线性求解器，3 https://www.hdfgroup.org/downloads/hdf5/ [2020/02/19]。4 http://www.gslib.com/gslib_help/sgsim.html[2020/02/26]5 https://github.com/geirev/EnKF_analysis [2020/02/19]。6 https://www.eocoe.eu/[访问日期：2020/02/19]。J. Keller，V.Rath，J.Bruckmann等人粤公网安备44010502000011号3Fig. 1. SHEMAT-Suite的分支及其在git-repository中的依赖项。2. 软件描述SHEMAT-Suite是一个用于模拟流动、热量和物质输运的正向和反向数值代码。SHEMAT-Suite有三种编译模式：用于纯正演计算的正演模式，用于确定性反演计算的自动微分模式，以及用于地质统计模拟、随机数据同化和参数估计的随机模式。大多数用户与SHEMAT-Suite的交互仅限于输入文件。对于前进模式，只需要一个主输入文件。对于反向模式，需要在主输入文件以及单独的输入文件中2.1. 软件功能SHEMAT-Suite的源代码分为两个级别，分支级别和目录级别。在分支级别，使用版本控制软件git将源代码划分为分支。在软件开发中，git-branch通常用于开发软件，特别是所谓的对于SHEMAT-Suite，这个概念稍微扩展到四个主分支（图1）。①的人。四个主分支源自软件的三个核心功能：（1）用于流量、热量和溶解物质传输的解决方案的正向模式，（2）用于贝叶斯反演的AD模式（ad），其采用基于导数的优化方案，使用切线线性和通过自动微分（AD）生成的伴随代码，用于计算精确的雅可比矩阵，以及（3）随机模式（sm），用于蒙特卡罗模拟、数据同化和使用EnKF的参数估计。四个分支中的第一个，master，只能在forward模式下编译，并且针对纯forward计算进行了优化。第二个分支master-ad可以在forward模式和AD模式下编译。第三个分支master-sm可以在正向模式和随机模式下编译。第四个也是最后一个分支，master-all，可以在所有三种模式下编译。一般来说，用户总是可以使用master-all。然而，较小的代码可以在存储器和计算时间方面为特定任务提供更好的性能。由于特殊的分支结构，SHEMAT-Suite开发人员必须以正确的方式引入更改，分支，然后将这些更改合并到下游分支中。影响正向计算的变更应添加到master，变更仅影响AD/逆计算到master-ad，变更仅影响Monte-Carlo框架和随机计算的master-sm。原则上，分支master-all应该只通过合并接收对其他分支很少有具体的变化以使上游分支中的改变兼容。SHEMAT-Suite的开源文档与SHEMAT-Suite的源代码一起提供。SHEMAT-Suite-Wiki7包括（1）关于SHEMAT-Suite的一般信息，（2）首次用户的教程，以及（3）关于使用三个主要功能的说明，前向计算，随机反演和确定性反演。对于开发人员来说，可以从源代码头文件生成Doxy- gen文档. 此外，开源git-repositorySHEMAT-Suite_Models-open8包含各种功能的测试模型SHEMAT-Suite和一个Python脚本，用于根据以前代码版本的结果测试SHEMAT-Suite_Models-open的README提供了最新的测试模型列表，如果可用，还提供了它们的分析基准。预处理和后处理的脚本在git-repositorySHEMAT-Suite_Scripts-open中提供。9例如，此存储库包含一个Python脚本，用于将ASCII输入文件转换为可由SHEMAT-Suite读取的hdf 5文件最后，我们计划逐步转向软件开发的持续集成方法。这将简化代码开发工作流程，最大限度地减少错误的引入，从而确保软件质量。2.2. 软件构架本节解释了前一节中描述的分支中的源代码目录（见图11）。2，用于正向计算的源代码目录的草图）。master分支包含以下目录：/forward/包含了大部分的前向计算代码站，包括处理输入和输出，时间离散化循环和非线性迭代循环。此代码是并行使用共享内存并行编程范式的OpenMP针对多核工作站。/solve/包含线性方程组的解。与/forward/的分离确保了附加求解器和预处理器的简单实现，并且对于自动微分软件的使用至关重要。/props/包含定义流体和岩石属性的动态行为和耦合的属性模块。 属性模块可设计为使用恒定或压力和温度相关的流体和岩石属性。此外，本发明还可以考虑溶解物质浓度（通常是盐浓度）对流体性质（例如密度）的影响。模块化结构允许添加特定于站点或应用程序的属性关系。/user/包含用户定义的输入子程序、输出子程序、用户定义的Fortran模块和通用在Picard迭代开始时调用的用户定义子例程。这些子程序的完整集合称为SHEMAT-Suite的用户模块。用户模块提供了一个简单的，7 https://git.rwth-aachen.de/SHEMAT-Suite/SHEMAT-Suite-open/-/wikis/教程[访问日期2020/03/04]。8 https://git.rwth-aachen.de/SHEMAT-Suite/SHEMAT-Suite_Models-打开/[访问2020/03/16]。9 https://git.rwth-aachen.de/SHEMAT-Suite/shemat-suite_scripts-打开[访问2020/03/17]。4J. Keller，V.Rath，J.Bruckmann等人粤公网安备44010502000011号图二、SH EM AT - Su ite 中 的前向计算图 ， 突出显示了 各种源代码目录的作用。灵活的方式来添加和调整特定问题的功能。开源SHEMAT-Suite在none目录中提供了一个示例性的模板用户模块。目录none包含不带功能的子例程模板，这些子例程模板用作开发新用户模块的起点。/hdf 5/包含用于输入和输出的HDF 5库接口。/cmake/包含使用CMake工具的编译实用程序。10/doc/包含用于生成Doxygen文档的输入。分支master-sm和master-ad添加以下目录。/simul/（sm）包含用于地质统计模拟、MC计算和EnKF更新的子例程。/inverse/（ad）包含确定性反转代码。/mkAD/（ad）包含用于生成AD代码的脚本。3. 说明性示例本节介绍一个使用SHEMAT-Suite进行模拟的简单示例。对于众多的科学应用实例，感兴趣的读者可以参考第1节和第4节以及表1中引用的文献。此示例的完整输入文件10 https://cmake.org/[访问日期：2020/02/26]。并且对于说明SHEMAT-Suite的各种功能的附加示例，可以在测试模型库SHEMAT-Suite_Models-open中找到。113.1. 正演计算Theis模型的数值解[42]用于说明SHEMAT-Suite的使用。Theis模型是一种解决方案本文从渗流方程出发，计算了水平承压含水层中抽水井引起的不稳定水位下降。我们将SHEMAT-Suite的数值解与解析解进行比较[13，43]。Code Ocean上提供了一个计算胶囊，用于重现说明性示例。12输入文件1显示了Theis问题的输入文件。输入文件的第一部分提供运行的一般信息，例如标题和属性模块的规范，在本例中是const。这意味着在整个模拟过程中，诸如流体密度和粘度的流体参数以及诸如渗透率的岩石参数保持恒定。第二部分定义网格，从x-y-z方向的单元数量开始，然后是单元大小。在这种情况下，网格在远离井处较粗，而在中心处较细11 https://git.rwth-aachen.de/SHEMAT-Suite/SHEMAT-Suite_Models-打开/[访问2020/03/16].12 https://doi.org/10.24433/CO.9420115.v1[暂定doi].J. Keller，V.Rath，J.Bruckmann等人粤公网安备44010502000011号5−·图三. THEIS模型的结果。（a）3.5天后的水头分布，（b）从（a）开始穿过抽水井的水头分布的水平断面，以及(c)沿横断面的数值和解析计算水头之间的差异在抽水井附近的模型。在时间步长控制下，指定瞬态模拟，并且将影响时间步进的一些求解器参数设置为它们的缺省值。有关这些特殊参数的更多信息，请参见SHEMAT-Suite文档。时间单位设置为一天（86 400秒），模拟运行3.5天，平均分为100个时间步长。然后，为线性求解器方法BiCGStab设置精度和最大迭代次数。侧水头Dirichlet边界条件设置为20 m，模型的顶部和底部为无流边界。在模型的中心，抽水井由单元格（35， 35， 1）处的汇表示。该汇被指定为泵送速率为0的情况。001 m3 s-1，导致地下水位下降可以与分析的泰斯公式进行比较的单元格。所有4761个网格单元的初始水头设置为20 m。设定流体（水）和岩石基质的性质为水，这些属性包括密度（1000 kg m−3），压缩-粘度（5 10−8m s2kg−1）和粘度（10−3kg m−1s−1）。为岩石基质的性质包括孔隙度（10%）、渗透率（10−12 m2）和压缩性（10−8 m s2 kg−1）。其他参数-TER输入包括热参数和用于物质传输的参数由于该示例性示例仅针对地下水流，因此这些附加参数不会影响模拟。在SHEMAT-Suite文档中，用户可以找到一些示例输入文件的完整演示，包括Theis模型的这个文件最后，将输出文件类型设置为VTK和HDF5。结果示于图3 .第三章。水头场显示了3.5天后抽油井周围的水位降深SHEMAT-Suite 的 HDF 5 输出 使用 Python 和matplotlib 可视 化[44]。与Theis解析解的比较4. 影响SHEMAT-Suite为多个地球科学领域做出了贡献SHEMAT-Suite的重要代码开发总结见表1。SHEMAT-Suite最重要的贡献是在地热能领域，这是软件开发的原始动机举几个最近的例子，安装在钻孔热交换器处的温度传感器已被用于使用SHEMAT-Suite[33，45]估计地下水速度。在渗流区模型中，高压电缆的热输出使用SHEMAT-Suite [46，47]建模。对于深层地热储层，使用代码[32，37]研究了突出地震反射体对温度分布和对流单体形成的影响。此外，目前的一个研究项目使用SHEMAT-Suite预测墨西哥火山口系统的地热潜力[48]。此外，使用SHEMAT-Suite的模拟有助于其他地球科学领域。开发了冰的属性模块[34]，目前用于模拟多孔介质中的冻结和解冻过程。在当前的研究项目中，SHEMAT- Suite被用于模拟海底地下水排泄6J. Keller，V.Rath，J.Bruckmann等人粤公网安备44010502000011号输入文件1：THEIS模型的输入文件！ #=>RUNINFO<<<<< =#！#标题Theis#linfo23 2 3#runmode0#主动式头部道具数量= const#USER = none！ #=>GRID<<<<<=#！#网格69 69 1#delx3*160.0 d0 1*80.0 d0 1*40.0 d0 3*20.0 d0 3*10.0 d05*5.0 d0 5*2.0 d0 11*1.0 d0 2 * 1.0 d0 1.0 d011*1.0 d0 5*2.0 d0 5*5.0 d0 3*10.0 d0 3*20.0 d01*40.0 d0 1*80.0 d0 3*160.0 d0#Dely3*160.0 d0 1*80.0 d0 1*40.0 d0 3*20.0 d0 3*10.0 d05*5.0 d0 5*2.0 d0 11*1.0 d0 2 * 1.0 d0 1.0 d011*1.0 d0 5*2.0 d0 5*5.0 d0 3*10.0 d0 3*20.0 d01*40.0 d0 1*80.0 d0 3*160.0 d0#delz20. D0！ #=====>>> >TIMESTEPCONTROL<<<<< ====#！#timetepcontrol11.0 1.0 1.0 0.0#tunit86400#时间间隔，记录或数据 =10.0 3.5100线！ #=>求解器控件<<<<<=====#！#nlsolve50 0#lsolvef（线性求解器控制）1.d-1064 500#非线性迭代控制）1.0d-101.0 d0！#====>>> >BOUNDARYANDDINITIALCONDITIONS<<<<<=====#！#headbcn，records=1，direction=03535 1-0.001 d0 0#headbcd，simple=right，value=init#headbcd，simple=back，value=init#头初始化4761*20.0 d0#tempinit4761*20.0 d0！#=>PROPERTIES<<<<<=====#！#液压支柱一千0d05. 0D-84218。0d00. 65d01. 0D-3 02d0单位数量0.1 D01. D01. d01d-12 1.d-81. D01. d0 3.35 d00. D02300000. D010天0. D00. D0！#=>几何<<<<<=====#！#uindex4761*1！#=>输出<<<<<=====#！#文件输出：vtk hdf自上一个冰河时代以来，新泽西州大陆架上[36]。SHEMAT-Suite的现代含水层流量模拟通常是基于集合的，并依赖于高性能计算[49]。SHEMAT-Suite也对计算地球科学产生了影响。使用基于自动微分的SHEMAT-Suite进行三维逆参数估计[11]。蒙特卡洛技术被应用于荷兰海牙地热储层中预期地热能使用的不确定性量化[50]。专门为SHEMAT-Suite开发了一种并行化方案，并已实施[7集合卡尔曼滤波器已被实施并广泛测试用于随机渗透率估计[10，50]。最优实验设计已被用于寻找地热储层中新钻孔的最佳位置[41]。代码的开放源代码将影响这些方法在地球科学中的更广泛使用我们现在将列出一些重要的事态发展，目前-目前不是SHEMAT-Suite但是计划在将来与开源出版物合并已经使用SHEMAT-Suite作为起点开发了多相流模块[28]。此外，还实施了等效裂缝模型方法，该方法通过基于模拟有限差分法的放大来表示裂缝渗透率[54，56]。由于SHEMAT-Suite是一个多孔介质模型，它不能模拟离散的裂缝模型，因此不能显式地通过裂缝流动或传输。然而，在某些条件下，通过对各向异性渗透率进行建模或应用随机连续统概念，可以近似计算通过裂缝岩石的流量[49]。模拟无压含水层的基本方法已经实施。包括电动势模拟模块，可用于研究多孔介质中流体流动产生的自势[57另一个模块在SHEMAT-Suite中实现，用于化学流体-岩石相互作用[55]。在当前的研究项目EooE中，J. Keller，V.Rath，J.Bruckmann等人粤公网安备44010502000011号7--表1SHEMAT-Suite的重要代码开发新实施的功能基于自动微分的逆参数估计Rath et al. 2006 [11]Mottaghy和Rath，2006年[34]用于不确定性量化和降低的Monte Carlo技术Vogt et al. 2010 [50]钻孔热交换器模块bMottaghy和Dijkshoorn，2012 [51]共享内存并行化Wolf，2011 [7]基于集合卡尔曼滤波器的数据同化Vogt et al. 2012 [31]使用自动微分的多相流模块aBüsing et al. 2014 [27]分布式内存并行化Rostami和Bücker，2014 [38]平面热激活岩土系统的热传递模型aKürten等人。2014年[52、53]使用全渗透率张量的各向异性流动模块aChen等人。2016 [54]使用自动微分的超临界水/蒸汽模块aBüsing et al. 2016 [29]通过最佳实验设计进行储层表征的最佳钻孔定位cSeidler et al. 2016 [41]盐底辟附近多孔沉积岩中的岩盐沉淀模型aLi et al.2017年[55]使用精确雅可比矩阵求解非均质多孔介质中的有效两相流aBüsing，2020 [28]a开放源码软件包中没有的功能。b开放源码软件包提供的简化功能cSHEMAT-Suite功能开放源码，需要额外软件实现SHEMAT-Suite与现代高性能软件的接口，例如用于数据处理的便携式数据接口（PDI13），用于并行求解器方法的PETSc[60]和用于并行实施集合卡尔曼滤波器的并行数据同化框架（PDAF [61最后，多相流实现使用软件PETSc [60]，以实现MPI并行并适用于分布式存储器CPU集群。5. 结论为了解决流动、热量和物质输运中出现的正演和反演问题，用户可以选择SHEMAT-Suite，原因有很多，包括（1）软件的模块化特征，特别是在软件开发中使用git分支的新颖性;（2）开箱即用的自动微分（AD），用于计算反演中的精确雅可比矩阵;（3）其广泛的应用领域，例如浅层和深层地热能、水文地质学、工程学;（4）SHEMAT-Suite的功能，包括：（1）软件的模块化特征，特别是在软件开发中使用git分支的新颖性;（2）自动微分（AD），用于计算反演中的精确雅可比矩阵;（3）其广泛的应用领域，例如浅层和深层地热能、水文地质学、工程学;（4）SHEMAT-Suite的功能，包括：（4）高性能计算能力，即：OpenMP实现，以及（5）确定性和随机反演的能力。SHEMAT-Suite源代码的最近开源出版物具有许多优点：（1）它使SHEMAT-Suite可用于更大的代码用户和开发者社区，（2）它使使用SHEMAT-Suite的现有研究可重现，（3）它激发了新的有趣的开发和代码添加，以及（4）它启动了SHEMAT-Suite用户社区的形成。竞合利益作者声明，他们没有已知的竞争性财务利益或个人关系，可能会影响本文报告的工作致谢SHEMAT-Suite的重要开发及其作者列于表1。此外，除了本文作者和表1中列出的作者外，以下研究人员还使用和测试了SHEMAT-Suite ： Sascha Bodenburg ， Paromita Deb ， Anozie Ebigbo ，Dominique Knapp ， Jan Niederau ， Ariel Thomas ， Mathis vanWickeren。Gabriele Margaret在SHEMAT-Suite中实现了原始版本的集合卡尔曼滤波器。Michael Kühn和Heike Stöfen写道：13 https://pdi.julien-bigot.fr/0.6.0/[访问日期：2020/03/17].Shiyuan Li集成到SHEMAT-Suite中的用于模拟化学流体-岩石相互作用的原始SHEMAT-code。SHEMAT-Suite的分布式内存并行化主要由Ali Rostami开发。这项工作的财政支持是由欧盟委员会地平线2020研究和创新计划根据赠款协议第676629号（EoCoE）和第824158号（EoCoE-II）提供的。在过去的15年里，许多不同的资助组织支持了SHEMAT-Suite的工作，我们非常感谢欧盟（EU）、德国科学基金会（DFG）和德国联邦政府的一些部委：教育和研究（BMBF），环境，自然保护和核安全（BMU）经济事务和能源局（BMWi）和德国联邦辐射防护局（BfS）。引用[1]Finsterle S，Sonnenthal EL，Spycher N.使用TOUGH模拟器套件ComputGeosci 2014;65：2网址：//dx.doi.org/10.1016/j.cageo.2013.06.009网站。[2] Diersch H-J. FEFLOW：多孔和裂缝介质中流动、质量和热输运的有限元模拟。Springer-VerlagBerlinHeidelberg;2014，http://dx.doi.org/10.1007/978-3-642-38739-5。[3] White MD，Oostrom M，Rockhold ML，Rosing M.使用STOMP模拟器对汉福德工厂四氯化碳迁移进行可扩展建模包气带J 2008;7（2）：654-66。http://dx.doi.org/10.2136/vzj2007.0070网站。[4]Kolditz O ， Bauer S ， Bilke L ， Böttcher N ， Delfs JO ， Fischer T ， et al.OpenGeoSys：an open-source initiative for numerical simulation of thermo-hydro-mechanical/chemical（THM/C）processes in porous media. Environ.地球Sci.2012;67（2）：589-99.http://dx.doi.org/10.1007/s12665-012-1546-x.[5][10] Flemisch B，Darcis M，Erbertseder K，Faigle B，Lauser A，MosthafK，et al. DuMu x：DUNE for multi-phase，component，scale，physics，.多孔介质中的流动和输运。Adv Water Resour 2011;34（9）：1102网址：//dx.doi.org/10.1016/j.advwatres.2011.03.007网站。[6]Hughes JD，Langevin CD，Banta ER. MODFLOW 6框架的文档。Tech.方法2017年。http://dx.doi.org/10.3133/tm6a57网站。[7]沃尔夫AEinSoftwarekonzeptzurhierarchischendelisierungvonstochastischenunddeterministischenInversionsproblemenaufmodernenccNUMA-PlattformenunterNutzungautomatischerProgramtransformation-tion，博士论文，RWTH Aachen University; 2011.[8]Wolf A，Rath V，Bücker HM.地热模拟包的并行化：四种多核架构的案例研究。In：Bischof C，Bücker M，GiovannP，Joubert G，Lippert T，MohrB ， Peters F ， editors. 并 行 计 算 ： 体 系 结 构 、 算 法 和 应 用 。 Advances inparallel computing ， vol. 15 ， Amsterdam ， The Netherlands ： IOS Press;2008，p. 451-8，与NIC-导演共同出版，作为NIC系列第38卷，2007年。网址http://hdl.handle.net/2128/6498。[9]Bücker HM，Rasch A，Rath V，Wolf A.地热模拟正问题和反问题的半自动并 行 化 。In ： Proceedings of the 24th ACM symposium on appliedcomputing，Honolulu，Hawaii，USA，March 8-12 ，2009，Vol. 2. NewYork：ACM Press; 2009，p. 971-5. http://dx.doi.org/10.1145/1529282.1529495网站。8J. Keller，V.Rath，J.Bruckmann等人粤公网安备44010502000011号[10]Keller J，Hendricks Franssen H-J，Marvedig G.比较七种变体 集合卡尔曼滤波 器 ： 需 要 多 少 合 成 实 验 ？ Water Resour Res 2018;54 （ 9 ） ： 6299-318.http://dx.doi.org/10的网站。1029/2018 wr023374。[11] Rath V，Wolf A，Bücker HM. 基于自动微分的地下水流与传热耦合三维联合反演：灵敏度计算、验证与合成实例。Geophys J Int 2006;167（1）：453http://dx.doi.org/10.1111/j.1365-246x.2006.03074.x网站。[12]克劳泽角UntersuchungenzurTrennungderkonduktivenundkonvektivenAnteile im Wärmetransport in einem Sedimentbecken amBeispiel desOberrheintalgrabens. Düsseldorf：Werner Verlag;1988.[13]Clauser C编辑热含水层中反应流的数值模拟：SHEMAT和处理SHEMAT。Springer-VerlagBerlinHeidelbergNewYork;2003，http://dx.doi.org/10.1007/978-3-642-55684-5.[14] Trescott PC，Larson S.三维地下水流模拟有限差分模型文献。公开文件报告，美国地质调查局; 1976年，http://dx.doi.org/10.3133/ofr76591。[15]作者：J.用强隐式方法求解三维地下水流方程。J Hydrol 1977;35（1-2）：49-60. http：//dx.doi.org/10.1016/0022-1694（77）90076-2.[16]林奇博士.数值偏微分方程为环境科学家和工程师：第一个实用课程. NewYork，NY：Springer Science&Business Media;2005.[17]Huyakorn PS.地下水流计算方法。New York，NY：Academic Press;2012.[18]van der Vorst HA. Bi-CGSTAB：Bi-CG的一种快速、平滑收敛的变体，用于求 解 非 对 称 线 性 系 统 。 SIAM J Sci Stat Comput 1992;13 （ 2 ） ： 631-44.http://dx.doi.org/10.1137/0913035网站。[19]Lawson CL，Hanson RJ，Kincaid DR，Krogh FT.基本线性代数子程序，供Fortran 使 用 .ACMTransMathSoftware1979;5 （ 3 ） ： 308-23.http://dx.doi.org/10.1145/355841.355847网站。[20]Dongarra JJ，Du Croz J，Hammarling S，
我的内容管理
展开
