Transcript Mesh Generator for HYDROGEOCHEM & TOUGH2

The application of mesh generator
on
HYDROGEOCHEM & TOUGH2
Presenter:
Adviser:
Date:2010/12/09
Outline
What is mesh generator?
Mesh generator is applied to HYDROGEOCHEM 4.0
Mesh generator is applied to TOUGH2
Future work
2
3
1 00 150 20 0 25 0 30 0 35 0 4 00 4 50 50 0 550
27 0 1 0 00
27 0 0 0 00
26 9 9 0 00
Y (m )
DTM
50
26 9 8 0 00
26 9 7 0 00
26 9 6 0 00
26 9 5 0 00
21 8 0 00
22 0000
2 220 00
X (m )
22 4000
4
Mesh generator
8000 nodes
DTM
46
60
0 4
48
80
0 5
50
00
0 5
52
20
0 5
54
40
0 5
56
60
0 5
58
80
0
4
2697000
460 480 500 520 540 560 580
2697000
Y (m )
2696900
Y
Y (m
(m))
2696900
2696800
2696800
2696700
2696700
221500
221600
X (m )
221700
221800
221500
221600
X (m )
221700
221800
5
HYDROGEOCHEM 4.0
Pressure
head
Total head
DTM
P H (m )
T
20
40
60
Z (m )
80 100 120 140 160 180
2 .6 9 9 0 E + 0 6
50 100 150 200 250 300 350 400 450 500 550
2 .6 9 9 0 E + 0 6
B.C. TH=0 m
2 .6 9 8 0 E + 0 6
Y (m )
Y (m )
2 .6 9 8 0 E + 0 6
2 .6 9 7 0 E + 0 6
I.C. TH=0 m
2 .6 9 7 0 E + 0 6
2 .6 9 6 0 E + 0 6
2 .6 9 6 0 E + 0 6
B.C. TH=200 m
2 .6 9 5 0 E + 0 6
2 .2 0 0 0 E + 0 5
2 .2 2 0 0 E + 0 5
X (m )
2 .6 9 5 0 E + 0 6
2 .2 0 0 0 E + 0 5
2 .2 2 0 0 E + 0 5
X (m )
6
BetaD2
Volume
D1
TOUGH2 Input Formats
β
D1
1 E1 1
EE
L
E L E 1 E LEEL
D2
ELE 2
EL
LE
E 22
E
ELE 2
7
ELEME----1----*----2----*----3----*----4----*----5----*----6----*----7----*----8
CONNE----1----*----2----*----3----*----4----*----5----*----6----*----7----*----8
ELE01 ELE02
Name
MaterialPERMVolumeD1 area for D2
Area
X
COS(β)
Y
Z
0000000034
10.5000D+000.5000D+000.7000D+010.0000D+00
heat exchange
0000000001MATE10.7000D+010.1400D+02
00000
20.3500D+010.3000D+010.1000D+010.1000D+01
0.5000D+000.3500D+010.0000D+00
0000100035MATE10.6000D+010.1200D+02
00001
10.5000D+000.5000D+000.6000D+010.0000D+00
0.5000D+000.1000D+020.0000D+00
0000100002MATE10.6000D+010.1200D+02
00002
20.3000D+010.3000D+010.1000D+010.1000D+01
0.5000D+000.1600D+020.0000D+00
0000200036MATE10.2500D+010.5000D+01
00003
10.5000D+000.5000D+000.6000D+010.0000D+00
0.5000D+000.2025D+020.0000D+00
0000200003MATE10.1000D+010.2000D+01
00004
20.3000D+010.1250D+010.1000D+010.1000D+01
0.5000D+000.2200D+020.0000D+00
0000300037MATE10.2500D+010.5000D+01
00005
10.5000D+000.5000D+000.2500D+010.0000D+00
0.5000D+000.2375D+020.0000D+00
0000300004
20.1250D+010.5000D+000.1000D+010.1000D+01
*
0000400038
10.5000D+000.5000D+000.1000D+010.0000D+00
*
0000400005MATE10.1500D+040.3000D+04
00945
20.5000D+000.1250D+010.1000D+010.1000D+01
0.5750D+040.1495D+030.0000D+00
0000500039MATE20.1500D+040.3000D+04
00946
10.5000D+000.5000D+000.2500D+010.0000D+00
0.5750D+040.1525D+030.0000D+00
0000500006MATE10.1500D+040.3000D+04
00947
20.1250D+010.3000D+010.1000D+010.1000D+01
0.5750D+040.1555D+030.0000D+00
0000600040MATE10.3000D+040.6000D+04
00948
10.5000D+000.5000D+000.6000D+010.0000D+00
0.5750D+040.1600D+030.0000D+00
0000600007MATE10.6000D+040.1200D+05
00949
20.3000D+010.6000D+010.1000D+010.1000D+01
0.5750D+040.1690D+030.0000D+00
00950
MATE10.3000D+040.6000D+04
0.5750D+040.1780D+030.0000D+00
8
TOUGH2 Synthetic cases
B.C. P=10000 pa
Vertical
Horizontal
200
200
200
10 0 0 years
Y
Z (m )
Y
Z (m )
19 5
36
08
0 y ears
150
100
100
100
50
50
50
0
0
100
50
X (m )
0
0
100
50
X (m )
0
0
100
50
X (m )
100
I.C. P=10000 pa
200
200
22 2
68
09
0 years
1 0200
5 000 yyears
ears
100
y ears
1 *1 02^790 0y ears
150
100
3 .6* 1 0 ^ 26 years
9 * 10 ^ 1 6 y ears
150
Y
Z (m )
150
200
Y
Z (m )
50
50
X (m )
Y
Z (m )
Z (m )
100
100
4 .6 E + 0 5
4 .4 E + 0 5
4 .2 E + 0 5
4 .0 E + 0 5 5 0
3 .8 E + 0 5
3 .6 E + 0 5
3 .4 E + 0 5
0
3 .2 E + 0 5
0
3 .0 E + 0 5
2 .8 E + 0 5
2 .6 E + 0 5
2 .4 E + 0 5
2 .2 E + 0 5 2 0 0
2 .0 E + 0 5
1 .8 E + 0 5
1 .6 E + 0 5
150
1 .4 E + 0 5
1 .2 E + 0 5
1 .0 E + 0 5
19
2
602
0 y ears
150
Y
Z (m )
P
1 5 04 .8 E + 0 5
150
Y
Z (m )
150
12
2
36
08
0 y ears
200
Y
Z (m )
200
100
100
50
50
B.C. P=50000 pa
50
0
0
0
0
50
50
XX
(m )
5 0 100
X (m )
0
0
1 0 0 X50X(m )
100
0
0
50
X (m )
100
0
0
50
X (m )
100
9
TOUGH2 Synthetic cases
Vertical
200
1 0 0 0 y e ars
150
Z (m )
150
100
50
0
0
2 .9 E + 0 5
2 .8 E + 0 5
2 .7 E + 0 5
2 .6 E + 0 5
2 .5 E + 0 5
2 .4 E + 0 5
2 .3 E + 0 5
2 .2 E + 0 5
2 .1 E + 0 5
2 .0 E + 0 5
1 .9 E + 0 5
1 .8 E + 0 5
1 .7 E + 0 5
1 .6 E + 0 5
1 .5 E + 0 5
1 .4 E + 0 5
1 .3 E + 0 5
1 .2 E + 0 5
1 .1 E + 0 5
B.C. P=10000 pa
100
50
0
0
500
200
I.C. P=10000 pa
1000
X (m )
1500
2000
B.C. P=300005 .7pa
* 1 0 ^ 1 8 y ea rs
5 10500
Z (m )
P
Z (m )
200
1000
1500
X (m )
2000
100
50
0
0
500
1000
X (m )
1500
2000
10
TOUGH2 Synthetic cases
Vertical
200
1 0 ^ 38 y e a rs
Z (m )
150
200
100
B.C. P=10000 pa
50
150
P
100
500
1000
X (m )
1500
2000
I.C. P=10000 pa
50
0
0
0
44.5 E + 0 5
44.0 E + 0 5
33.5 E + 0 5
33.0 E + 0 5
22.5 E + 0 5
22.0 E + 0 5
11.5 E + 0 5
200
1 0 ^ 716 0y ey aers
a rs
150
Z (m )
Z (m )
0
5 10000
1000
X (m )
50
B.C. P=50000 pa
0
500
1000
X (m )
2000
1500
2000
200
12 0* ^1 170 0^y1ye0ae rs
ayrs
e a rs
150
Z (m )
7259
elements
0
1500
100
50
0
11
0
500
1000
X (m )
1500
2000
Illustrative Examples
1100m, 11MPa
Upper Shale Cap (impermeable)
150
Shale layers
Y (m)
30 meters
30 meters
100
30 meters
50
Hydrostatic pressure
30 meters
Domain: 6000m×180m
Injection well: (0,22)
Injection rate: 0.1585 kg/s
Porosity:
sand- 35%
shale- 10.25%
Permeability:
sand- 3×10-12m2
shale- 1×10-14m2
Temperature: 37 ℃
Salinity: 3.2 wt% NaCl
Pressure of upper boundary: 11MPa
30 meters
Injection Well
0
22 meters
Lower Shale Cap (impermeable)
0
2000
X (m)
4000
6000
12
SG: 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7
Inject 2 years
Z(m)
150
100
50
well point
0
0
1000
2000
3000
X(m)
4000
5000
13
Future work
14
XC O2 
C O 2 m ass
1000 g ( H 2 O m ass )  N aC l m ass  C O 2 m ass
8000
Inject 2 year
XCO2
Z(m)
7500
0.035
0.03
0.025
0.02
0.015
0.01
0.005
7000
6500
28000
30000
32000
X(m)
34000
36000
8000
Inject 20 year
Z(m)
7500
7000
6500
28000
30000
32000
X(m)
34000
36000
15
6500
6500
7500
7000
Z(m)
Z(m)
0.03
0.025
0.02
0.015
0.01
0.005
28000 6500
7000
30000
32000
34000
Future
28000 X(m)work
30000
32000
X(m)
36000
34000
6500
28000
28000
6500
30000
28000
28000
30000
32000
30000
X(m)
30000
34000
32000
X(m)
32000
X(m)36000 34000
32000
X(m)
34000
34000
XCO2a
0.035
0.03
0.025
0.02
0.015
0.01
0.005
XCO2a
0.035
0.03
0.025
0.02
0.015
16
Thanks for
your attention