Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
Stochastic analysis of fluid flow and tracer pathways in crystalline fracture networks
KTH, School of Architecture and the Built Environment (ABE), Land and Water Resources Engineering. (Water Resources Engineering)
2010 (English)Doctoral thesis, comprehensive summary (Other academic)
Abstract [en]

Understanding groundwater flow systems and how these control transport is an essential part in assessing the suitability of subsurface environments as hosts for storage of toxic waste. Therefore it is important to be able to integrate knowledge obtained from field characterisation of the subsurface with methods which can be used to evaluate and predict possible impact on surrounding environments.In this thesis I investigate the characteristics of flow and transport in discrete fracture networks by analysing Eulerian and Lagrangian descriptions within a stochastic framework. The analysis is conducted through numerical flow and transport simulations configured according to available field data, combined with independent theoretical analytic and semi-analytic methods which are able to reveal insight to relevant constitutive properties. It is shown that numerical simulations conducted with the discrete fracture network approach can be both conditioned and confirmed against field measurable quantities, and the developed theoretical methods are evaluated against results obtained from simulation. Thereby, a methodology which can provide links between field measurable quantities and tracer discharge is presented, developed and evaluated. It is shown to be robust with respect to underlying assumptions used for flow configurations.In particular, a specific sampling algorithm for obtaining a Lagrangian description of transport based on a Eulerian description of flow is proposed, evaluated and shown to be robust for the cases considered, providing accurate replications. Also a generalisation of both the advection-dispersion solution and the one-sided stable distribution is shown to be able to evaluate advective transport quantities, and combined with a Lagrangian retention model it is shown to be a fairly accurate and robust method for upscaling distributions, enabling predictions of transport in terms of tracer discharge. Evaluation of transport is also conducted against the advective-dispersion assumption, where results indicate advective transport is generally non-Fickian for the fracture networks and domain scales considered, but not necessarily anomalous. Additionally, the impact certain model assumptions have on tracer discharge are analysed. For example, transport is evaluated for assumptions regarding injection mode, fracture network heterogeneity, relationship between aperture and transmissivity, relationship between transmissivity and size, as well as scale and modelling dimension. In relation to hydraulic testing and flow analysis, a method for conditioning fracture transmissivity from field measurements of flow by simulation is developed and evaluated against homogenisation assumptions commonly used in field applications. Results indicate the homogenisation assumption generally fails for current interpretations of field data.

Abstract [sv]

Miljökonsekvensbedömningar av toxiskt avfall i djupt bergförvar kräver engrundläggande förståelse av grundvattenströmning samt hur detta påverkartransportfenomet. Därför är det viktigt att kunna integrerafältundersökningsdata från berggrundsmätningar med metoder som kan användas föratt utvärdera och förutsäga potentiella konsekvenser på omgivningen.I denna avhandling undersöker jag flödes- och transportegenskaper i diskretaspricknätverk genom stokastisk analys av eulerska och lagrangeskafältbeskrivningar. Analysen sker genom en kombination av dels numeriska flödes-och transportsimuleringar som är konfigurerade enligt tillgänglig fältdata, samtdels med oberoende teoretiska analytiska och semi-analytiska metoder som gerdjupare insyn i relevanta konstitutiva egenskaper. Resultat visar att dennumeriska simuleringsmetoden för diskreta spricknätverk kan både konditionerastill fältdata och bestyrkas gentemot mätbara kvantiter. Detta är av betydelse dåde teoretiska metoderna i sin tur är främst evaluerade gentemotsimuleringsresultat. Därmed utvecklas en bestyrkt metodik som kansammanlänka och i viss mån omvandla fältdata till uppskattningar av mängdenspårämnen i ett utflöde. Resultat indikerar att denna metodik är robust avseendeflera antaganden som har används i simuleringskonfigurationen.En särskild urvalsalgoritm introduceras som kan erhålla en lagrangesktransportbeskrivning utifrån ett eulerskt strömningsfält. även denna utvärderasavseende vissa simuleringsantaganden och resultat tyder på att den är robust förde undersökta fallen. Vidare föreslås en viss generalisering av lösningen tillden advektiva-dispersionsekvationen samt av ensidigt stabila (one-sided stable)sannolikhetsfördelningar som metod för att prediktera advektiva kvantitetergenom upskalning av transportfördelningar i rummet. Denna modell kombineras meden tidigare utvecklad metod för transportretention för att uppskatta reaktivagenombrottsfördelningar. Således blir det möjligt att prediktera reaktivtransport d v s rumslig upskalning av genombrottstider för spårämnestransport.Metoden används också för att evaluera ett linärt dispersionsantagande, därresultat indikerar att även advektiv transport kan påvisa icke-linärt beteende.Transport i spricknätverk utvärderas bland annat för modellantaganden avseendeinjektionsmetod, heterogenitet i spricknätverk, konstitutiva relationer mellanapertur och transmissivitet samt mellan transmissivitet och spricklängd, ochmodelleringsskala samt dimension. Beträffande hydrauliska testmetoder ochflödesanalys introduceras en simuleringsmetod för att konditioneraspricktransmissivitet från flödesmätningar. Detta jämförs med etthomogeniseringsantaganden som inte sällan används i fältundersökningar för atttolka flödesmätningar till spricktransmissivitet, och resultat tyder på attdetta antagande kan betydligt undervärdera transmissivitet.

Place, publisher, year, edition, pages
Stockholm: KTH , 2010. , xi, 62 p.
Series
Trita-LWR. PHD, ISSN 1650-8602 ; 1056
Keyword [en]
Groundwater; Fractured media; Transport; Retention; Macrodispersion; Lagrangian formulation; Stochastic modelling
National Category
Earth and Related Environmental Sciences Oceanography, Hydrology and Water Resources Other Environmental Engineering
Identifiers
URN: urn:nbn:se:kth:diva-11974ISBN: 978-91-7415-560-0 (print)OAI: oai:DiVA.org:kth-11974DiVA: diva2:291690
Public defence
2010-02-18, F3, Lindstedtsvägen 26, Kungl Tekniska högskolan, 10:00 (English)
Opponent
Supervisors
Available from: 2010-02-03 Created: 2010-02-02 Last updated: 2018-01-12
List of papers
1. Upscaling particle transport in discrete fracture networks: 1. Nonreactive tracers
Open this publication in new window or tab >>Upscaling particle transport in discrete fracture networks: 1. Nonreactive tracers
2007 (English)In: Water resources research, ISSN 0043-1397, E-ISSN 1944-7973, Vol. 43, no 10, W10428- p.Article in journal (Refereed) Published
Abstract [en]

We study tracer transport through discrete fracture networks and develop a methodology for upscaling particle breakthrough curves on the basis of fracture segment data. Our prime interest is to model the early arrival and peak of tracer breakthrough curves, i.e., to capture the bulk of the tracer mass arrival. This study is based on two-dimensional discrete fracture network simulations, combined with a truncated one-sided stable distribution as a model for upscaling particle transitions. Results indicate that this model can accurately capture the bulk mass and peak of the breakthrough distributions for an upscaled distance of at least 1 order of magnitude in terms of transport scale, which for our simulations is about 2 orders of magnitude greater than the mean fracture segment scale. We also introduce an accurate mapping algorithm for transforming Eulerian into Lagrangian flow statistics, without a priori knowledge of network connectivity.

Keyword
NONLOCAL REACTIVE TRANSPORT, POROUS-MEDIA, CHEMICAL HETEROGENEITY;, ANOMALOUS TRANSPORT, MULTIPLE FAMILIES, CRYSTALLINE ROCK, MASS-TRANSPORT, FLOW, DISPERSION, MODEL
National Category
Oceanography, Hydrology and Water Resources Other Environmental Engineering
Identifiers
urn:nbn:se:kth:diva-11967 (URN)10.1029/2006WR005334 (DOI)000250525500001 ()2-s2.0-36749058309 (Scopus ID)
Note
QC 20101005Available from: 2010-02-02 Created: 2010-02-02 Last updated: 2018-01-12Bibliographically approved
2. Upscaling particle transport in discrete fracture networks:  2. Reactive tracers
Open this publication in new window or tab >>Upscaling particle transport in discrete fracture networks:  2. Reactive tracers
2007 (English)In: Water resources research, ISSN 0043-1397, E-ISSN 1944-7973, Vol. 43, no 10, W10429- p.Article in journal (Refereed) Published
Abstract [en]

We study sorbing tracer transport through discrete fracture networks using astochastic Lagrangian framework, combined with the methodology for upscalingparticle breakthrough curves developed in the first part of this article series.Results indicate that this procedure can accurately predict expected normalizedtracer discharge for an upscaled distance of 1 order of magnitude in terms oftransport scale, which for our simulations is about 2 orders of magnitudegreater than the mean fracture segment scale. Specifically, we show theimportance of retaining the correlation between the water residence time τ andthe hydrodynamic control of retention β in order to make accurate tracerdischarge predictions. Also, we show that the extreme tails of τ and β distributions have essentially no impact on tracer discharge. These results areillustrated using the unlimited diffusion model, and for two hypotheticaltracers with properties designed to capture the behavior of many commonlyoccurring natural radionuclides.

Place, publisher, year, edition, pages
American Geophysical Union, 2007
Keyword
Stochastic transport, Upscaling, Retention, Fractured media
National Category
Oceanography, Hydrology and Water Resources Other Environmental Engineering
Identifiers
urn:nbn:se:kth:diva-11969 (URN)10.1029/2006WR005336 (DOI)000250525500002 ()
Note
QC 20101005Available from: 2010-02-02 Created: 2010-02-02 Last updated: 2018-01-12Bibliographically approved
3. Significance of injection modes and heterogeneity on spatial and temporal dispersion of advecting particles in two-dimensional discrete fracture networks
Open this publication in new window or tab >>Significance of injection modes and heterogeneity on spatial and temporal dispersion of advecting particles in two-dimensional discrete fracture networks
2009 (English)In: Advances in Water Resources, ISSN 0309-1708, E-ISSN 1872-9657, Vol. 32, no 5, 649-658 p.Article in journal (Refereed) Published
Abstract [en]

We investigate the impact of injection mode (flux and resident injection) and heterogeneity in hydraulic properties on dispersion of advecting particles in two-dimensional discrete fracture network models, using a Monte Carlo method. We find that the injection mode has a significant effect on dispersion: The resident injection mode exhibits anomalous features of transport whereas the flux injection mode tends to Gaussian transport; this observation is easily understood by considering phase diagrams where a limited number of particles entering low velocity fractures greatly increase macrodispersion. In spite of a sizeable portion of negative longitudinal velocities, it is shown that multiple crossings are negligible when quantifying longitudinal macrodispersion. A simple probabilistic expression of particle mass balance is shown to predict well the spatial distribution of advecting particles.

Place, publisher, year, edition, pages
Elsevier B.V., 2009
Keyword
Macrodispersion, Discrete fractured media, Advective transport, Injection mode, Multiple crossing, Anomalous transport, Gaussian transport, Stochastic simulation
National Category
Earth and Related Environmental Sciences Oceanography, Hydrology and Water Resources Other Environmental Engineering
Identifiers
urn:nbn:se:kth:diva-11970 (URN)10.1016/j.advwatres.2008.07.010 (DOI)000266253900003 ()2-s2.0-64849101017 (Scopus ID)
Note
QC 20100929Available from: 2010-02-02 Created: 2010-02-02 Last updated: 2018-01-12Bibliographically approved
4. Inference of field scale fracture transmissivities in crystalline rock using flow log measurements
Open this publication in new window or tab >>Inference of field scale fracture transmissivities in crystalline rock using flow log measurements
2010 (English)In: Water resources research, ISSN 0043-1397, E-ISSN 1944-7973, Vol. 46, no 11, W11502- p.Article in journal (Refereed) Published
Abstract [en]

Characterization of transmissivity for crystalline rock is conducted through simulation by conditioning against borehole flow rates obtained from high-resolution, in situ field measurements during extraction pumping. Full three-dimensional discrete fracture network simulations are carried out according to specifications obtained from site characterization data in a stochastic Monte Carlo setting. A novel method of conditioning is thereby introduced and applied using nonparametric comparison tests, which provide quantifiable measures of accuracy enabling evaluation of simulated results against field measurements. The assumption of a constitutive relationship ( perfect correlation) between fracture size and transmissivity is adopted. The method is evaluated against both single and multiple realizations, various domain size, and fracture length configurations and shown to be robust for the cases considered. When the introduced method of conditioning is applied, transmissivity parameterization can be inferred to a narrow range with a quantifiable accuracy in terms of a probability value. Results indicate that elementary interpretation of transmissivity based on homogenization of a porous medium will generally underestimate transmissivity. Further implications on advective transport for natural flow conditions are briefly evaluated, indicating advective breakthrough times can be overestimated up to a factor of about 10 in the median.

Keyword
Fractured media, Transmissivity, Pump tests, Conditioning, Stochastic simulation
National Category
Earth and Related Environmental Sciences Oceanography, Hydrology and Water Resources Other Environmental Engineering
Identifiers
urn:nbn:se:kth:diva-11971 (URN)10.1029/2009WR008367 (DOI)000283951400001 ()2-s2.0-78149306039 (Scopus ID)
Note
QC 20101208. Uppdaterad från submitted till published (20101208).Available from: 2010-02-02 Created: 2010-02-02 Last updated: 2018-01-12Bibliographically approved
5. Numerical and analytical modelling of advective travel times in realistic three-dimensional fracture networks
Open this publication in new window or tab >>Numerical and analytical modelling of advective travel times in realistic three-dimensional fracture networks
2011 (English)In: Water resources research, ISSN 0043-1397, E-ISSN 1944-7973, Vol. 47, no W02506, 16pp- p.Article in journal (Refereed) Published
Abstract [en]

Travel time distributions obtained from advective transport in multiple realisations of realistic discrete fracture network simulations are analysed using the truncated one-sided stable distribution, which has previously been shown to generalise both the advection-dispersion solution as well as one-sided stable distributions. Using this model it is shown that the Fickian assumption inherent in the advection-dispersion equation generally fails, despite that the first two moments of travel time essentially scale linearly with distance. It is also observed that the equally probable realisations drawn from the ensemble can produce a wide range of behaviour under the current configuration, such that Fickian conditions are almost obtained in some cases for increasing scales. Based on a small-scale calibration against particle breakthrough, the model is then shown to successfully predict limiting bounds of transport for a one order of magnitude increase in scale. Correlation in particle velocity is explicitly shown to be significant for scales close to the characteristic Lagrangian segment length. The network configuration is obtained from extensive site characterisation data at the Laxemar region in Sweden, and represents a block scale domain of reasonably sparse background fractures.

Keyword
Fractured media; Transport; Advection; Fickian dispersion; One-sided stable densities; Stochastic simulation
National Category
Earth and Related Environmental Sciences Oceanography, Hydrology and Water Resources Other Environmental Engineering
Identifiers
urn:nbn:se:kth:diva-11972 (URN)10.1029/2010WR009290 (DOI)000287039900001 ()2-s2.0-79851496854 (Scopus ID)
Note
QC 20110209. Uppdaterad från "accepted" till "published"(20110209).Available from: 2010-02-02 Created: 2010-02-02 Last updated: 2018-01-12Bibliographically approved
6. Transport and retention from single to multiple fractures in crystalline rock at Äspö (Sweden): 2. Fracture network simulations and generic retention model.
Open this publication in new window or tab >>Transport and retention from single to multiple fractures in crystalline rock at Äspö (Sweden): 2. Fracture network simulations and generic retention model.
2010 (English)In: Water resources research, ISSN 0043-1397, E-ISSN 1944-7973Article in journal (Refereed) Published
Abstract [en]

Hydrogeologic characterization of crystalline rock formations on the field scale is important for many applications but still presents a multitude of challenges [Neuman, 2005]. In this work we use comprehensive hydro-structural information and present a detailed simulation study of flow and advective transport in a discrete fracture network (DFN) that replicates the TRUE Block Scale rock volume at the Äspö Hard Rock Laboratory (Sweden). Simulated water residence time τ and hydrodynamic retention parameter β are used as independent constraints for estimating material retention properties as presented in part 1 of this series [Cvetkovic et al., 2009] whereas simulated mean water residence times are compared with observed values. We find that the DFN simulations reproduce water residence times reasonably well, indicating that the characterization data are sufficient and that the DFN model does capture dominant features of the flow paths analyzed. The empirical quadratic law that relates aperture and transmissivity seems to better reproduce calibrated mean water residence times than the theoretical cubic law, for the five flow paths. The active specific surface area (β/τ) [1/L] as inferred from simulations, is used for defining a generic retention model for the dominant rock type (”Äspö diorite”) that matches fairly well the entire range of calibrated retention parameters of the TRUE tests. The combination of Part 1 and this work provide a general, comprehensive methodology for evaluating tracer test results in crystalline rock where a comparable amount of information is available; critical to this methodology is that tracer tests are carried out using tracers with sufficiently different sorption affinities (of factor 10 to 100).

Place, publisher, year, edition, pages
American Geophysical Union, 2010
National Category
Engineering and Technology
Identifiers
urn:nbn:se:kth:diva-11973 (URN)10.1029/2009WR008030 (DOI)000277266300004 ()2-s2.0-77949637043 (Scopus ID)
Note
QC 20110103Available from: 2010-02-02 Created: 2010-02-02 Last updated: 2017-12-12Bibliographically approved

Open Access in DiVA

fulltext(5286 kB)574 downloads
File information
File name FULLTEXT01.pdfFile size 5286 kBChecksum SHA-512
f9eea3959acfd1c29ed3c31e0f65a7fbc5b0afa877777de73305909d66d310cdc1938f43288e0f4913ad3fffb3e7615a6bd90fe0fa05597ba82db53948a59a56
Type fulltextMimetype application/pdf

Search in DiVA

By author/editor
Frampton, Andrew
By organisation
Land and Water Resources Engineering
Earth and Related Environmental SciencesOceanography, Hydrology and Water ResourcesOther Environmental Engineering

Search outside of DiVA

GoogleGoogle Scholar
Total: 574 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

isbn
urn-nbn

Altmetric score

isbn
urn-nbn
Total: 557 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf