image
image
image
Česká geologická služba
Mathematical Modelling of the RESIBIL Project
3D geological model – methodology
(sources, method, results)
Skácelová Z., Mlčoch B.
Workshop Prague, March 4, 2019

image
image
image
Česká geologická služba
Creating of the 3D geological model
:
Depth of surface level
of individual geological layers:
Geology – individual strata (e.g. crystaline basement, Carboniferous,
Permian, Cretaceous, volcanic rocks)
Ingeneering geology (Cenozoic – Quternary)
Hydrogeology – aquifers and aquitards (e.g. Cretaceous: A -
Cenomanian, B - Lower Turonian, C - Middle and Upper Turonian),
Source
of data:
Geological maps
Boreholes database
Geophysics (electrical and electromagnetic sounding, seismic)
Definition of
faults
(simplified tectonic view)
Determining the
range
of layer (strata)

image
image
image
image
image
Česká geologická služba
3D modelling in CGS:
Boreholes
Database:
Geophysics
Map of scale 1 : 50 000
or 1 : 25 000
Outcrops
of individual geological
units
(altitude of relief)
Conversion to altitude
3D model of surface of individual geological units –
basement, aquifers
Geological maps:
Depth
of selected interfaces
+
Input file to modelling

image
image
image
image
Česká geologická služba
Boreholes processing.
Microsoft Access Application –
input
– data from CGS borehole database (.mdb format)

image
image
image
Česká geologická služba
Boreholes processing.
Klic
GDO
Vrt
nazev
x
y
z
Rozhra
ni
Typ
rozhra
ni
Hloubka
rozhrani
ve vrtu
Z
rozhrani
Mapa ZM
25
Horni
na
Přívl
astek
Přívlastek
mineralog
ický
Zrn
ito
st
Obl
ast
Regio
n
Jedn
otka
S
o
u
Vr
st
va
S
t
r
Pozná
mka
Klíč
GEO
Hloubka
vrtu
Autor
Datum
100965 L-10
-883730,8
-1006280,4 485,74
CMK
S
32
453,74
11-123 SVR
F
F2
F2B
33,3 Mlčoch 15.11.2011
101144 PJ-1
-886380
-1006650
510
CMK
S
7
503
11-123 SVR
KVR GRN
F
F2
F2C
50 Mlčoch 15.11.2011
110655 JD686
-851059,7
-1002297,3
495,6
CMK
S
43,8
451,8
11-212 120
BTT
3 F
F3
NEJ
Mlčoch 15.11.2011
112452 OD 543
-851609,3
-1002450,8
508,6
CMK
S
26
482,6
11-212 120
BTT
F
F3
NEJ
nepøesnì popsaný profil, rozhraní
Mlčoch mùž
15.11.2011
e být níž e
114927 M-135
-857136,43 -1008824,43 408,32
CMK
S
10
398,32
11-213 120
BTT
F
F3
NEJ
37,90 - 46,80 m aplit 66 Mlčoch 15.11.2011
122630 JE36
-853987,5
-1010626,8
390
CMK
S
40,6
349,4
11-214 120
B+A
3 F
F3
NEJ
porfyrický
Mlčoch 15.11.2011
125839 36H
-864602,2 -1012599,86 410,46
CMK
S
214
196,46
11-231 PRR
B+M
F
F2
F2B
252 Mlčoch 15.11.2011
110 TP 5
-774912,5
-974610,3 229,57
TR
S
182,8
46,77
02-322 porfyr
DVJ
F
F2
F2A
okatá
76,45 Mlčoch 11.10.2010
101 VR 24
-782817,38
-973696,73
324,5
TR
S
46,8
277,7
02-321 porfyr
DVJ
F
F2
F2A
164,59 Mlčoch 11.10.2010
55641 TH-33
-776803
-972036,4
311,3
TR
S
47,5
263,8
02-321 315
F
TR
47.50-51.00 m ztráta jádra
Mlčoch 11.10.2010
DIAMO 3367/65 -848003,07
-1003589,2 478,52
CMK
S
82,1
396,42
11-212 912
F
F2
Y
84,1 Mlčoch 15.11.2011
1 3061
-852786,5
-1007851
424,7
CMK
S
22,6
402,1
11-214 ELV
F
F3A
31,9 Mlčoch 15.11.2011
1 5051
-852753,2
-1007551,2
432
CMK
S
44,3
387,7
11-214 ELV
F
F3A
52 Mlčoch 15.11.2011
Microsoft Access Application –
output
– data in Excel table .xls
Storage
of processed data – database CGS Postgre SQL
Input
– .xls format
solid structure of table
not all atributes are required

image
image
image
image
image
image
Česká geologická služba
Geophysical data processing.
Electrical
and electromagnetic sounding
Gravity modelling
Reflection seismic
Refraction seismic
D
ata from CGS archives – technical reports with the coordinates of profiles, points
Processing – manual, determination of depth

image
image
image
image
Česká geologická služba
Geophysical data processing.
Storage of processed data – database CGS Postgre SQL
Input
– .xls format
solid structure of table
not all atributes are required
Output
– data fill in Excel table .xls – points of mesurement

image
image
image
image
image
image
Česká geologická služba
Outcrops – geological map.
Storage of processed data – database CGS Postgre SQL
Input
– .xls format
solid structure of table
Polygons of outcrops
(e.g. basement, granite
massif, Cretaceous –
Cenomanian, Turonian,
Coniacian, volcanic rocks) –
recent relief (ASCII format
.txt – contour line or grid)

image
image
image
image
image
Česká geologická služba
Storage of processed data
database CGS Postgre SQL
Search after:
ID
Coordinates
Name of boreholes
Geophysical method
Strata, rocks types etc.
Export:
xls. format
csv. format

image
image
image
image
image
image
Česká geologická služba
3D model in Surfer software
Method:
Use
Surfer
software
Processing all available data for each strata (boreholes, geophysisc) to input
file (.dat) including outcrops on the surface - relief (from geological maps).
Gridding methods:
a) triangulation with linear interpolation, b) kriging
Combination of surface
Upper Paleozoic
Recent relief
Crystalline basement
a)
b)

image
image
image
image
image
image
Česká geologická služba
3D model in Surfer software
Model accuracy – number of data and coverage of study area
Crystalline basement
Permian and Carboniferous
Quaternary

image
image
image
image
image
image
image
Česká geologická služba
3D model in Surfer software
Correction
of data error
Modelling without and with faults
Grid point density selection – precision of 3D
model
Surfer
– suitable SW for checking and removing
data errors, fast and user-friendly

image
image
image
image
Česká geologická služba
Outputs
: isolines of surface strata (grid)
Options presentation: e.g. Contour line map of surface, map of thickness,
3D model, cross-section
Thickness of Tertiary sediments in Cheb and
Sokolov basins
3D model in Surfer - results

image
image
image
image
Česká geologická služba
ResiBil – 3D geological model of Děčínský Sněžník pilot area – aquifers and aquitards
3D model in Surfer - results

image
image
image
image
Česká geologická služba
ResiBil study area – 3D model of basement
Lusatian thrust
Lusatian fault
Stráž fault
3D Model of aguifers A, B, BC – February 2019
3D model in Surfer - results

image
image
image
image
image
Česká geologická služba
Project REBILANCE – 3D model of aquifers in area 3
+ aquifer + aquitard D
aquifer + aquitard C,
BC, Cb, Ca
basement
+ aquifer + aquitard BC
+aquifer A
+aquifer+aquitard B
3D model in Surfer - results

image
image
image
image
image
image
Česká geologická služba
N
S
W
E
Lusatian Fault
Česká Kamenice Basin
Česká Kamenice Basin
Stráž Fault
The 3D geological model allows to
create depth cross-section:
a) cross-section of Cretaceous
Basin – aquifers A, B, BC, C
b) cross-section of ATC ryolite
3D model in Surfer - results
Ryolite body in Altenberg-Teplice
Caldera
a)
b)

image
image
image
image
image
image
Česká geologická služba
Lysá
Nymburk
Poděbrady
3D model in Surfer - results
Presentation of geological model for public:
a) Model of Cretaceous sediments on Upper Paleozoic to
determine of slope instability
b) Model of Quternary sediments (aquifer) in Labe area
with the volume calculation of groundwater horizons
678 839 132 m
3
a)

image
image
image
image
Česká geologická služba
Thanks for your attention.