Flow and transport modelling of fractured aquifers based on a geostatistical approach [Elektronische Ressource] / von Anongnart Assteerawatt

Informations

• Publié par: universitat_stuttgart
• Publié le: 01 janvier 2008
• Nombre de lectures: 44
• Langue: English
• Poids de l'ouvrage: 8 Mo

Extrait

Flow and Transport Modelling of Fractured Aquifers
based on a Geostatistical Approach






Von der Fakultät Bau- und Umweltingenieurwissenschaften der
Universität Stuttgart zur Erlangung der Würde eines
Doktor-Ingenieurs (Dr.-Ing.) genehmigte Abhandlung



Vorgelegt von
Anongnart Assteerawatt
aus Bangkok, Thailand




Hauptberichter: Prof. Dr.-Ing Rainer Helmig
Mitberichter: Prof. Dr. rer. nat. Dr.-Ing. habil. András Bárdossy
Prof. Dr.-Ing Reinhard Hinkelmann

Tag der mündlichen Prüfung: 21. Juli 2008









Institut für Wasserbau der Universität Stuttgart
2008





Heft 176 Flow and Transport Modelling
of Fractured Aquifers based on
a Geostatistical Approach


von
Dr.-Ing.
Anongnart Assteerawatt













Eigenverlag des Instituts für Wasserbau der Universität Stuttgart D93 Flow and Transport Modelling of Fractured Aquifers based on
a Geostatistical Approach























Titelaufnahme der Deutschen Bibliothek


Assteerawatt, Anongnart:
Flow and Transport Modelling of Fractured Aquifers based on a Geostatistical
Approach / von Anongnart Assteerawatt. Institut für
Wasserbau, Universität Stuttgart. - Stuttgart: Inst. für Wasserbau, 2008

(Mitteilungen / Institut für Wasserbau, Universität Stuttgart: H. 176)
Zugl.: Stuttgart, Univ., Diss., 2008)
ISBN 978-3-933761-80-4

NE: Institut für Wasserbau <Stuttgart>: Mitteilungen


Gegen Vervielfältigung und Übersetzung bestehen keine Einwände, es wird lediglich
um Quellenangabe gebeten.




Herausgegeben 2008 vom Eigenverlag des Instituts für Wasserbau
Druck: Document Center S. Kästl, Ostfildern Acknowledgement
First of all, I would like to thank my main supervisor Prof. Rainer Helmig for giv-
ing me an opportunity to do my research at the Institute of Hydraulic Engineering
(IWS), University Stuttgart. His valuable discussions, guidance, support and indi-
vidual characters have motivated and encouraged me throughout my time here. I
alsowouldliketoexpressmygrateful gratitudetomyco-supervisors Prof. Andra´s
Ba´rdossyandProf. ReinhardHinkelmann.
Furthermore, I would like to thank to Ha˚kon Hægland from the University of
Bergen for the collaboration in the development of the streamline tracing in a
fracture-matrix system. Iwould also liketo acknowledge Philipp Haasfor hishelp
onGermantranslation ofthesummaryofmywork.
ManythanksgoalsotoallcolleaguesatIWS,particularlytheHydrosyswhosefruit-
fulacademicdiscussionsaswellascheerfulnon-scientificactivitiescreateuniquely
special and memorable working environment. Special thanks to Prudence Law-
dayforproofreadingmydissertationandforhelponmyadministrativedocuments
here.
I extend my thanks to all my Thai friends for being accompanied throughout the
time in Germany. Finally, I amextremely thankful to myfamily, especiallymypar-
entsandmysisters,fortheirmentalsupportandendlessencouragement.Contents
List of Figures ii
List of Tables vii
Nomenclature viii
Abstract xi
Zusammenfassung xiii
1 Introduction 1
1.1 Motivation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Objectives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Structure oftheWork . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2 From Nature to Conceptual Models 5
2.1 GeometricalCharacteristicsofFractured Aquifers . . . . . . . . . . . 5
2.1.1 Fractureorientation . . . . . . . . . . . . . . . . . . . . . . . . 7
2.1.2 Fractureaperture . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.1.3 Fracturesizeandfracture trace . . . . . . . . . . . . . . . . . . 9
2.1.4 Fracturedensity . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.1.5 Fracturespace . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.2 ModelConceptsforFracture Aquifers . . . . . . . . . . . . . . . . . . 11
2.2.1 Scalesandrepresentativeelementaryvolume(REV) . . . . . . 11
2.2.2 Continuummodelapproach . . . . . . . . . . . . . . . . . . . 12
2.2.3 Discretemodelapproach . . . . . . . . . . . . . . . . . . . . . 13
2.2.4 Scalesofconsideration . . . . . . . . . . . . . . . . . . . . . . . 15
3 Stochastic Characteristics of Fracture Networks 18
3.1 Statistics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3.1.1 Univariatestatistics andbivariatestatistics . . . . . . . . . . . 18
3.1.2 Probabilitydistribution andrandomvariables . . . . . . . . . 19
3.2 Geostatistics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
3.2.1 Stationaryprocess . . . . . . . . . . . . . . . . . . . . . . . . . 23
3.2.2 Variogrammodels . . . . . . . . . . . . . . . . . . . . . . . . . 24ii Contents
3.3 StatisticalAnalysisofFractureCharacteristics . . . . . . . . . . . . . . 26
3.3.1 Fracture orientation . . . . . . . . . . . . . . . . . . . . . . . . 26
3.3.2 Fracture aperture . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3.3.3 Fracture sizeandfracturetrace . . . . . . . . . . . . . . . . . . 27
3.3.4 Fracture space . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.4 GeostatisticalAnalysisofFractureCharacteristics . . . . . . . . . . . 29
4 Generation of Fracture Network: Geostatistical Approach 30
4.1 GeneralConceptofFractureNetworkGeneration . . . . . . . . . . . 30
4.2 DescriptionoftheStudySite . . . . . . . . . . . . . . . . . . . . . . . . 31
4.3 AnalysisofSpatialCharacterisisticsofPliezhausenBlock . . . . . . . 32
4.4 GeostatisticalFracture Generation . . . . . . . . . . . . . . . . . . . . 42
4.4.1 Globaloptimizationtechnique: simulatedannealing . . . . . 42
4.4.2 Indicator-fieldgeneration . . . . . . . . . . . . . . . . . . . . . 44
4.4.3 Fracture-network generation . . . . . . . . . . . . . . . . . . . 46
5 Mathematical and Numerical Model 51
5.1 FormulationofFlowandTransportEquations . . . . . . . . . . . . . 51
5.1.1 Flowprocess . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
5.1.2 Transportprocess . . . . . . . . . . . . . . . . . . . . . . . . . . 55
5.2 DiscretizationofFlowandTransportEquations. . . . . . . . . . . . . 57
5.3 StreamlineTracing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
5.4 FluxRecoveryforStreamlineSimulation . . . . . . . . . . . . . . . . . 64
6 Comparison of Geostatistical and Statistical Fracture Generation 68
6.1 DescriptionofStudyCases. . . . . . . . . . . . . . . . . . . . . . . . . 68
6.2 Structure ofFracture Network . . . . . . . . . . . . . . . . . . . . . . . 76
6.3 FlowandTransportinFracture-Matrix System . . . . . . . . . . . . . 81
6.3.1 Modelset-up . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
6.3.2 Parameterscharacterizingsystembehavior . . . . . . . . . . . 83
6.3.3 Meshsizeanalysis . . . . . . . . . . . . . . . . . . . . . . . . . 84
6.3.4 Comparisonofstudycases . . . . . . . . . . . . . . . . . . . . 88
7 Comparison of Streamline Simulation and Advective-Diffusive Transport 94
7.1 MassTransportinStreamlines . . . . . . . . . . . . . . . . . . . . . . . 94
7.2 ComparisonStudyinPreliminaryTestCases . . . . . . . . . . . . . . 96
7.2.1 GroupA:Homogeneousandheterogeneousdomains . . . . . 96
7.2.2 GroupB:Singlefracture . . . . . . . . . . . . . . . . . . . . . . 100
7.2.3 GroupC:Systematicallydistributedfractures . . . . . . . . . 102
7.3 ComparisonStudyinComplexFractureNetworks . . . . . . . . . . . 104
8 Conclusion and Outlook 108List of Figures
1 Verfahren einer modifizierten Scanline-Technik zur Erzeugung eines
IndikatorfeldesauseinerOberfla¨chenkartierung. . . . . . . . . . . . . xvi
2 IllustrationderNachbarschafts-Richtungen mitdenZiffern0bis7. . xvii
3 Die gescannte Ansicht das Feldversuchsblocks von Nordwesten (a)
undSu¨dosten(b)(modifiziertnachDietrichetal.[2005]). . . . . . . . xviii
4 StandardisierteexperimentelleVariogrammederfu¨nfAußenseiten. . xix
5 Nachbarschaftsparameterderfu¨nfAußenseiten. . . . . . . . . . . . . xx
6 Kluftzellen-Dichtederfu¨nfAußenseiten. . . . . . . . . . . . . . . . . xx
2.1 Determination of a fracture orientation: azimuth, strike and dip
(modifiedafterDietrichetal.[2005]). . . . . . . . . . . . . . . . . . . . 7
2.2 Determination of polar coordinates of a unit normal vector of a frac-
turesurface. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.3 Fracture-orientation-distribution diagram from a Pliezhausen block
(afterDietrichetal.[2005]). . . . . . . . . . . . . . . . . . . . . . . . . 8
2.4 Twoopposite fracture surfacesofasingle fracture intersecting acore
sample of ca. 20 cm diameter (in cooperation with the Institute for
RoboticsandProcessControl, TechnicalUniversityofBraunschweig,
Germany)(afterSilberhorn-Hemminger[2002]). . . . . . . . . . . . . 9
2.5 Illustrationoffracturespacesfromasca