Skip to content

Commit 4fa99d3

Browse files
authored
Merge pull request #6 from OPM/development
gas-water co2store
2 parents 92d3d15 + 2db7f5a commit 4fa99d3

File tree

16 files changed

+238
-153
lines changed

16 files changed

+238
-153
lines changed

docs/_sources/configuration_file.rst.txt

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -30,18 +30,18 @@ The following input lines in the configuration file are:
3030
3131
"""Set the model parameters"""
3232
spe11c master #Name of the spe case (spe11a, spe11b, or spe11c) and OPM Flow version (master or release)
33-
complete #Name of the co2 model (immiscible or complete)
33+
complete gaswater #Name of the co2 model (immiscible or complete) and co2store implementation (gaswater or gasoil [oil properties are set to water internally in OPM flow])
3434
cartesian #Type of grid (cartesian, tensor, or corner-point)
35-
8400. 5000. 1200. #Length, width, and depth [m]
35+
8400 5000 1200 #Length, width, and depth [m]
3636
84 #If cartesian, number of x cells [-]; otherwise, variable array of x-refinment
3737
50 #If cartesian, number of y cells [-]; otherwise, variable array of y-refinment [-] (for spe11c)
38-
12 #If cartesian, number of z cells [-]; if tensor, variable array of z-refinment; if corner-point, fix array of z-refinment (17 entries)
39-
70. 36.25 #Temperature bottom and top rig [C]
40-
19620000. #Pressure on the top [Pa]
38+
12 #If cartesian, number of z cells [-]; if tensor, variable array of z-refinment; if corner-point, fix array of z-refinment (18 entries)
39+
70 36.25 #Temperature bottom and top rig [C]
40+
19620000 0.1 #Pressure on the top [Pa] and multiplier for the permeability in the z direction [-]
4141
1e-9 2e-8 0 #Diffusion (in liquid and gas) [m^2/s] and dispersion [m] (disperison only available in Flow master)
42-
8.5e-1 2500. #Rock specific heat and density (for spe11b/c)
43-
5e4 1. #Pore volume on lateral boundaries and width of buffer cell [m] (for spe11b/c)
44-
150. 10. #Elevation of the parabola and back [m] (for spe11c)
42+
8.5e-1 2500 #Rock specific heat and density (for spe11b/c)
43+
0 5e4 1 #Added pore volume on top boundary (for spe11a [if 0, free flow bc]), pore volume on lateral boundaries, and width of buffer cell [m] (for spe11b/c)
44+
150 10 #Elevation of the parabola and back [m] (for spe11c)
4545
4646
In line 5 you specify if you are using OPM Flow from the master branch or from the latest stable release (OPM-flow 2023.10 release).
4747
This since there are continuous changues in the OPM master branch (e.g., implementation of mechanical dispersion). Then we
@@ -52,7 +52,7 @@ with the defined number of elements in lines 9 to 11. The tensor grid allows to
5252
is first divided with the number of entries in the array, and after it divides each of these elements by the assigned number in
5353
the array entry. The corner-point mode generates a grid where the x and y direction are defined as in the array mode, but the
5454
cell faces in the z-direction follows the lines as defined in the `lines_coordinates.geo <https://github.com/OPM/pyopmspe11/blob/main/src/pyopmspe11/reference_mesh/lines_coordinates.geo>`_ script,
55-
resulting in 17 levels. Then, the z-refinment in each of these levels is set. See the configuration files in the `tests <https://github.com/OPM/pyopmspe11/blob/main/tests>`_ and
55+
resulting in 18 levels. Then, the z-refinment in each of these levels is set. See the configuration files in the `tests <https://github.com/OPM/pyopmspe11/blob/main/tests>`_ and
5656
`examples <https://github.com/OPM/pyopmspe11/blob/main/examples>`_ folder for the setting of these grids.
5757

5858
.. figure:: figs/cartesian.png
@@ -61,7 +61,7 @@ resulting in 17 levels. Then, the z-refinment in each of these levels is set. Se
6161

6262
Examples of a cartesian (top), tensor (middle), and corner-point (bottom) grids. The cartesian is generated using 170, 100, and 120
6363
elements, the tensor grid by setting the x-array to 50,100,20; the y-array to 5,10,70,10,5; and the z-array to 5,10,20,20,20,20,20,5; and
64-
the corner-point grid using the same xy-arrays as in the tensor grid and for the z-refinment 4,4,3,3,5,3,7,5,13,5,7,5,3,9,21,21,2.
64+
the corner-point grid using the same xy-arrays as in the tensor grid and for the z-refinment 4,4,3,3,5,3,7,5,13,5,4,3,5,3,9,21,21,2.
6565

6666
.. warning::
6767
Dispersion seems to work fine for the spe11a case where there are no thermal effects nor water evaporation. However, for spe11b/c there are

docs/configuration_file.html

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -102,18 +102,18 @@ <h2>Reservoir-related parameters<a class="headerlink" href="#reservoir-related-p
102102
<p>The following input lines in the configuration file are:</p>
103103
<div class="highlight-python notranslate"><div class="highlight"><pre><span></span><span class="linenos"> 4</span><span class="sd">&quot;&quot;&quot;Set the model parameters&quot;&quot;&quot;</span>
104104
<span class="linenos"> 5</span><span class="n">spe11c</span> <span class="n">master</span> <span class="c1">#Name of the spe case (spe11a, spe11b, or spe11c) and OPM Flow version (master or release)</span>
105-
<span class="linenos"> 6</span><span class="n">complete</span> <span class="c1">#Name of the co2 model (immiscible or complete)</span>
105+
<span class="linenos"> 6</span><span class="n">complete</span> <span class="n">gaswater</span> <span class="c1">#Name of the co2 model (immiscible or complete) and co2store implementation (gaswater or gasoil [oil properties are set to water internally in OPM flow])</span>
106106
<span class="linenos"> 7</span><span class="n">cartesian</span> <span class="c1">#Type of grid (cartesian, tensor, or corner-point)</span>
107-
<span class="linenos"> 8</span><span class="mf">8400.</span> <span class="mf">5000.</span> <span class="mf">1200.</span> <span class="c1">#Length, width, and depth [m]</span>
107+
<span class="linenos"> 8</span><span class="mi">8400</span> <span class="mi">5000</span> <span class="mi">1200</span> <span class="c1">#Length, width, and depth [m]</span>
108108
<span class="linenos"> 9</span><span class="mi">84</span> <span class="c1">#If cartesian, number of x cells [-]; otherwise, variable array of x-refinment</span>
109109
<span class="linenos">10</span><span class="mi">50</span> <span class="c1">#If cartesian, number of y cells [-]; otherwise, variable array of y-refinment [-] (for spe11c)</span>
110-
<span class="linenos">11</span><span class="mi">12</span> <span class="c1">#If cartesian, number of z cells [-]; if tensor, variable array of z-refinment; if corner-point, fix array of z-refinment (17 entries)</span>
111-
<span class="linenos">12</span><span class="mf">70.</span> <span class="mf">36.25</span> <span class="c1">#Temperature bottom and top rig [C]</span>
112-
<span class="linenos">13</span><span class="mf">19620000.</span> <span class="c1">#Pressure on the top [Pa]</span>
110+
<span class="linenos">11</span><span class="mi">12</span> <span class="c1">#If cartesian, number of z cells [-]; if tensor, variable array of z-refinment; if corner-point, fix array of z-refinment (18 entries)</span>
111+
<span class="linenos">12</span><span class="mi">70</span> <span class="mf">36.25</span> <span class="c1">#Temperature bottom and top rig [C]</span>
112+
<span class="linenos">13</span><span class="mi">19620000</span> <span class="mf">0.1</span> <span class="c1">#Pressure on the top [Pa] and multiplier for the permeability in the z direction [-]</span>
113113
<span class="linenos">14</span><span class="mf">1e-9</span> <span class="mf">2e-8</span> <span class="mi">0</span> <span class="c1">#Diffusion (in liquid and gas) [m^2/s] and dispersion [m] (disperison only available in Flow master)</span>
114-
<span class="linenos">15</span><span class="mf">8.5e-1</span> <span class="mf">2500.</span> <span class="c1">#Rock specific heat and density (for spe11b/c)</span>
115-
<span class="linenos">16</span><span class="mf">5e4</span> <span class="mf">1.</span> <span class="c1">#Pore volume on lateral boundaries and width of buffer cell [m] (for spe11b/c)</span>
116-
<span class="linenos">17</span><span class="mf">150.</span> <span class="mf">10.</span> <span class="c1">#Elevation of the parabola and back [m] (for spe11c)</span>
114+
<span class="linenos">15</span><span class="mf">8.5e-1</span> <span class="mi">2500</span> <span class="c1">#Rock specific heat and density (for spe11b/c)</span>
115+
<span class="linenos">16</span><span class="mi">0</span> <span class="mf">5e4</span> <span class="mi">1</span> <span class="c1">#Added pore volume on top boundary (for spe11a [if 0, free flow bc]), pore volume on lateral boundaries, and width of buffer cell [m] (for spe11b/c)</span>
116+
<span class="linenos">17</span><span class="mi">150</span> <span class="mi">10</span> <span class="c1">#Elevation of the parabola and back [m] (for spe11c)</span>
117117
</pre></div>
118118
</div>
119119
<p>In line 5 you specify if you are using OPM Flow from the master branch or from the latest stable release (OPM-flow 2023.10 release).
@@ -125,7 +125,7 @@ <h2>Reservoir-related parameters<a class="headerlink" href="#reservoir-related-p
125125
is first divided with the number of entries in the array, and after it divides each of these elements by the assigned number in
126126
the array entry. The corner-point mode generates a grid where the x and y direction are defined as in the array mode, but the
127127
cell faces in the z-direction follows the lines as defined in the <a class="reference external" href="https://github.com/OPM/pyopmspe11/blob/main/src/pyopmspe11/reference_mesh/lines_coordinates.geo">lines_coordinates.geo</a> script,
128-
resulting in 17 levels. Then, the z-refinment in each of these levels is set. See the configuration files in the <a class="reference external" href="https://github.com/OPM/pyopmspe11/blob/main/tests">tests</a> and
128+
resulting in 18 levels. Then, the z-refinment in each of these levels is set. See the configuration files in the <a class="reference external" href="https://github.com/OPM/pyopmspe11/blob/main/tests">tests</a> and
129129
<a class="reference external" href="https://github.com/OPM/pyopmspe11/blob/main/examples">examples</a> folder for the setting of these grids.</p>
130130
<figure class="align-default">
131131
<img alt="_images/cartesian.png" src="_images/cartesian.png" />
@@ -138,7 +138,7 @@ <h2>Reservoir-related parameters<a class="headerlink" href="#reservoir-related-p
138138
<figcaption>
139139
<p><span class="caption-text">Examples of a cartesian (top), tensor (middle), and corner-point (bottom) grids. The cartesian is generated using 170, 100, and 120
140140
elements, the tensor grid by setting the x-array to 50,100,20; the y-array to 5,10,70,10,5; and the z-array to 5,10,20,20,20,20,20,5; and
141-
the corner-point grid using the same xy-arrays as in the tensor grid and for the z-refinment 4,4,3,3,5,3,7,5,13,5,7,5,3,9,21,21,2.</span><a class="headerlink" href="#id1" title="Link to this image"></a></p>
141+
the corner-point grid using the same xy-arrays as in the tensor grid and for the z-refinment 4,4,3,3,5,3,7,5,13,5,4,3,5,3,9,21,21,2.</span><a class="headerlink" href="#id1" title="Link to this image"></a></p>
142142
</figcaption>
143143
</figure>
144144
<div class="admonition warning">

docs/searchindex.js

Lines changed: 1 addition & 1 deletion
Some generated files are not rendered by default. Learn more about customizing how changed files appear on GitHub.

docs/text/configuration_file.rst

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -30,18 +30,18 @@ The following input lines in the configuration file are:
3030
3131
"""Set the model parameters"""
3232
spe11c master #Name of the spe case (spe11a, spe11b, or spe11c) and OPM Flow version (master or release)
33-
complete #Name of the co2 model (immiscible or complete)
33+
complete gaswater #Name of the co2 model (immiscible or complete) and co2store implementation (gaswater or gasoil [oil properties are set to water internally in OPM flow])
3434
cartesian #Type of grid (cartesian, tensor, or corner-point)
35-
8400. 5000. 1200. #Length, width, and depth [m]
35+
8400 5000 1200 #Length, width, and depth [m]
3636
84 #If cartesian, number of x cells [-]; otherwise, variable array of x-refinment
3737
50 #If cartesian, number of y cells [-]; otherwise, variable array of y-refinment [-] (for spe11c)
38-
12 #If cartesian, number of z cells [-]; if tensor, variable array of z-refinment; if corner-point, fix array of z-refinment (17 entries)
39-
70. 36.25 #Temperature bottom and top rig [C]
40-
19620000. #Pressure on the top [Pa]
38+
12 #If cartesian, number of z cells [-]; if tensor, variable array of z-refinment; if corner-point, fix array of z-refinment (18 entries)
39+
70 36.25 #Temperature bottom and top rig [C]
40+
19620000 0.1 #Pressure on the top [Pa] and multiplier for the permeability in the z direction [-]
4141
1e-9 2e-8 0 #Diffusion (in liquid and gas) [m^2/s] and dispersion [m] (disperison only available in Flow master)
42-
8.5e-1 2500. #Rock specific heat and density (for spe11b/c)
43-
5e4 1. #Pore volume on lateral boundaries and width of buffer cell [m] (for spe11b/c)
44-
150. 10. #Elevation of the parabola and back [m] (for spe11c)
42+
8.5e-1 2500 #Rock specific heat and density (for spe11b/c)
43+
0 5e4 1 #Added pore volume on top boundary (for spe11a [if 0, free flow bc]), pore volume on lateral boundaries, and width of buffer cell [m] (for spe11b/c)
44+
150 10 #Elevation of the parabola and back [m] (for spe11c)
4545
4646
In line 5 you specify if you are using OPM Flow from the master branch or from the latest stable release (OPM-flow 2023.10 release).
4747
This since there are continuous changues in the OPM master branch (e.g., implementation of mechanical dispersion). Then we
@@ -52,7 +52,7 @@ with the defined number of elements in lines 9 to 11. The tensor grid allows to
5252
is first divided with the number of entries in the array, and after it divides each of these elements by the assigned number in
5353
the array entry. The corner-point mode generates a grid where the x and y direction are defined as in the array mode, but the
5454
cell faces in the z-direction follows the lines as defined in the `lines_coordinates.geo <https://github.com/OPM/pyopmspe11/blob/main/src/pyopmspe11/reference_mesh/lines_coordinates.geo>`_ script,
55-
resulting in 17 levels. Then, the z-refinment in each of these levels is set. See the configuration files in the `tests <https://github.com/OPM/pyopmspe11/blob/main/tests>`_ and
55+
resulting in 18 levels. Then, the z-refinment in each of these levels is set. See the configuration files in the `tests <https://github.com/OPM/pyopmspe11/blob/main/tests>`_ and
5656
`examples <https://github.com/OPM/pyopmspe11/blob/main/examples>`_ folder for the setting of these grids.
5757

5858
.. figure:: figs/cartesian.png
@@ -61,7 +61,7 @@ resulting in 17 levels. Then, the z-refinment in each of these levels is set. Se
6161

6262
Examples of a cartesian (top), tensor (middle), and corner-point (bottom) grids. The cartesian is generated using 170, 100, and 120
6363
elements, the tensor grid by setting the x-array to 50,100,20; the y-array to 5,10,70,10,5; and the z-array to 5,10,20,20,20,20,20,5; and
64-
the corner-point grid using the same xy-arrays as in the tensor grid and for the z-refinment 4,4,3,3,5,3,7,5,13,5,7,5,3,9,21,21,2.
64+
the corner-point grid using the same xy-arrays as in the tensor grid and for the z-refinment 4,4,3,3,5,3,7,5,13,5,4,3,5,3,9,21,21,2.
6565

6666
.. warning::
6767
Dispersion seems to work fine for the spe11a case where there are no thermal effects nor water evaporation. However, for spe11b/c there are

examples/spe11a.txt

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -3,17 +3,17 @@ flow --linear-solver=cprw --enable-tuning=true --enable-opm-rst-file=true --outp
33

44
"""Set the model parameters"""
55
spe11a master #Name of the spe case (spe11a, spe11b, or spe11c) and OPM Flow version (master or release)
6-
complete #Name of the co2 model (immiscible or complete)
6+
complete gasoil #Name of the co2 model (immiscible or complete) and co2store implementation (gaswater or gasoil [oil properties are set to water internally in OPM flow])
77
tensor #Type of grid (cartesian, tensor, or corner-point)
88
2.8 0.01 1.2 #Length, width, and depth [m]
99
8,15,20,20,15,15,15,6 #If cartesian, number of x cells [-]; otherwise, variable array of x-refinment
1010
1 #If cartesian, number of y cells [-]; otherwise, variable array of y-refinment [-] (for spe11c)
11-
4,6,6,6,6,6,6,6,6,6,6,6,6,6,6,5 #If cartesian, number of z cells [-]; if tensor, variable array of z-refinment; if corner-point, fix array of z-refinment (18 entries)
12-
20. 20. #Temperature bottom and top rig [C]
13-
110000. #Pressure on the top [Pa]
11+
4,6,6,6,6,6,6,6,6,6,6,6,6,6,6,5 #If cartesian, number of z cells [-]; if tensor, variable array of z-refinment; if corner-point, fix array of z-refinment (18 entries)
12+
20 20 #Temperature bottom and top rig [C]
13+
110000 1 #Pressure on the top [Pa] and multiplier for the permeability in the z direction [-]
1414
1e-9 1.6e-5 1e-2 #Diffusion (in liquid and gas) [m^2/s] and dispersion [m] (disperison only available in Flow master)
1515
0 0 #Rock specific heat and density (for spe11b/c)
16-
0 0 #Pore volume on lateral boundaries and width of buffer cell [m] (for spe11b/c)
16+
0 0 0 #Added pore volume on top boundary (for spe11a [if 0, free flow bc]), pore volume on lateral boundaries, and width of buffer cell [m] (for spe11b/c)
1717
0 0 #Elevation of the parabola [m] (for spe11c)
1818

1919
"""Set the saturation functions"""

examples/spe11b.txt

Lines changed: 11 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -3,17 +3,17 @@ flow --linear-solver=cprw --enable-tuning=true --enable-opm-rst-file=true --outp
33

44
"""Set the model parameters"""
55
spe11b master #Name of the spe case (spe11a, spe11b, or spe11c) and OPM Flow version (master or release)
6-
complete #Name of the co2 model (immiscible or complete)
6+
complete gaswater #Name of the co2 model (immiscible or complete) and co2store implementation (gaswater or gasoil [oil properties are set to water internally in OPM flow])
77
corner-point #Type of grid (cartesian, tensor, or corner-point)
8-
8400. 1. 1200. #Length, width, and depth [m]
8+
8400 1 1200 #Length, width, and depth [m]
99
4,40,20,20,20,10,10,4 #If cartesian, number of x cells [-]; otherwise, variable array of x-refinment
1010
1 #If cartesian, number of y cells [-]; otherwise, variable array of y-refinment [-] (for spe11c)
1111
2,2,1,1,2,1,3,2,6,2,3,2,2,8,4,10,10,1 #If cartesian, number of z cells [-]; if tensor, variable array of z-refinment; if corner-point, fix array of z-refinment (18 entries)
12-
70. 40. #Temperature bottom and top rig [C]
13-
19620000. #Pressure on the top [Pa]
12+
70 40 #Temperature bottom and top rig [C]
13+
19620000 0.1 #Pressure on the top [Pa] and multiplier for the permeability in the z direction [-]
1414
1e-9 2e-8 0 #Diffusion (in liquid and gas) [m^2/s] and dispersion [m] (disperison only available in Flow master)
15-
8.5e-1 2500. #Rock specific heat and density (for spe11b/c)
16-
5e4 1. #Pore volume on lateral boundaries and width of buffer cell [m] (for spe11b/c)
15+
8.5e-1 2500 #Rock specific heat and density (for spe11b/c)
16+
0 5e4 1 #Added pore volume on top boundary (for spe11a [if 0, free flow bc]), pore volume on lateral boundaries, and width of buffer cell [m] (for spe11b/c)
1717
0 0 #Elevation of the parabola and back [m] (for spe11c)
1818

1919
"""Set the saturation functions"""
@@ -49,8 +49,8 @@ PERM7 0 PORO7 0 THCONR7 2.00
4949

5050
"""Define the injection values ([hours] for spe11a; [years] for spe11b/c)"""
5151
"""injection time, time step size to write results, maximum solver time step, injected fluid (0 water, 1 co2) (well1), injection rate [kg/s] (well1), temperature [C] (well1), injected fluid (0 water, 1 co2) (well2), ..."""
52-
995 995 199 0 0 10 0 0 10
53-
5 5 5 0 0 10 0 0 10
54-
25 5 .1 1 0.035 10 0 0 10
55-
25 5 .1 1 0.035 10 1 0.035 10
56-
950 5 1 0 0 10 0 0 10
52+
995 995 1 0 0 10 0 0 10
53+
5 5 1 0 0 10 0 0 10
54+
25 5 .1 1 0.035 10 0 0 10
55+
25 5 .1 1 0.035 10 1 0.035 10
56+
950 5 1 0 0 10 0 0 10

0 commit comments

Comments
 (0)