2021-05-09

Relative permeability calculation using an unsteady state approach

标题：用非稳态方法计算相对渗透率

In this example, an unsteady state (pressure difference between inlet and outlet) method will be used to simulate flow through a permeable material (in this case, a spherepack). Additionally, we show how to download and simulate on volumes hosted at the digitalrockportal.

简介：在本例中，将使用非稳态(入口和出口之间的压差)方法来模拟通过渗透性材料(在本例中是一个球包)的流动。此外，我们还展示了如何下载和模拟驻留在数字岩石门户上的例子。

Percolating path of a non-wetting fluid (rock and wetting fluid not shown).

非润湿流体的渗透路径(未显示岩石和润湿流体)。

1.Get geometry: Fetch the raw image from the digitalrockportal:

1.获取几何图形:从数字岩石门户获取原始图像:

fetch_image.sh

or manually by downloading it from here

或者从这里手动下载

2.Pre-processing: Run

2.预处理:

create_geom_4_2phase.m

to create a computationally efficient (and artifact free) domain for simulation

为模拟创建一个计算效率高(且没有工件)的域

3.Two-Phase simulation: Then, run

3.两相模拟

run2-phase.sh

To perform the simulation with the inputs specified in the 2-phase .xml-file at several capillary pressure increments. For more information about this procedure, please refer to the validation_YoungLaplace example.

使用两阶段.xml文件中指定的输入，在几个毛细管压力增量下执行模拟。有关此过程的更多信息，请参考validation_YoungLaplace示例。

4.Post-processing: Then, run

4.后处理

create_geoms_4_kr.m

to calculate saturations for the pressure-saturation curve and create geometries for running the relative permeability simulations. The Pc-Saturation curve can then be plot with the generated saturations and the pressure inputs for the simulation.

计算压力-饱和度曲线的饱和度，并为相对渗透率模拟创建几何图形。pc -饱和度曲线可以用生成的饱和度和模拟的压力输入来绘制。

5.One-Phase simulation: Finally, run

单向流

run1-phase.sh

to run the simulation for generating the relative permeability values with the inputs specified in the 1-phase .xml-file.

使用1阶段.xml文件中指定的输入运行模拟以生成相对渗透率值。

The capillary Pc-Saturation curve looks like this:

毛细管pc -饱和度曲线是这样的，示例见原网站。

Resulting Pc-Saturation curve

结果Pc-Saturation曲线

Additional Options

1.Handling different input file types

1.处理不同的输入文件类型

If you want to run simulations on your projects with different input file types like-

如果你想在你的项目上运行模拟与不同的输入文件类型，如-

Type 1	Raw image file (as shown above)
Type 2	MATLAB (.MAT) file
Type 3	image slices like .tiff/.png/.jpg
Type 4	spherical grain centre coordinates & radii

First, add your initial geometry to the input folder Then, update your specific file type in the script and run

首先，将初始的几何图形添加到输入文件夹中，然后更新脚本中的特定文件类型并运行

create_geom_4_2phase_more_file_types.m

2.Running the toolbox on supercomputing resources (eg. TACC)

2.在超级计算资源上运行工具箱(例如TACC)

 

1. Copy the toolbox and the created domain for simulation to your workspace on the supercomputer
2. Update details in the example script batchjob_4TACC.sh and make (compile) the (2-phase and 1-phase) codes.
3. Start the LBM simulation by running

1. 将工具箱和创建的用于模拟的域复制到超级计算机上的工作区中
2. 更新示例脚本batchjob_4TACC.sh中的详细信息，并生成(编译)(两阶段和一阶段)代码。

3. 通过运行启动LBM模拟

sbatch batchjob_4TACC.sh

Notes-

1. Make sure to test smaller cases first before running a full-fledged simulation on the supercomputer
2. Observe best practices of making/running jobs from appropriate nodes

1.在超级计算机上运行一个完整的模拟之前，请确保先测试较小的案例

2.观察从适当节点制造/运行作业的最佳实践

3.Percolation path calculation and visualization

3.渗流路径计算和可视化

You can measure the length and visualize the fluid path at first breakthrough after finishing the simulation. To do this, update details and run

在完成模拟后，您可以测量长度并在第一次突破时可视化流体路径。为此，请更新详细信息并运行

percolation_path.m

This can create 3D visualizations like the one shown below.

这可以创建如下所示的3D可视化，示例见原网站。

Percolating path of a non-wetting fluid (rock and wetting fluid not shown).

非润湿流体的渗透路径(未显示岩石和润湿流体)。

Note: This function may not work in OCTAVE

注意:此功能可能不能在OCTAVE中工作