Voici les éléments 1 - 3 sur 3
  • Publication
    Accès libre
    Inverse modeling of groundwater flow in the Rhône alluvial aquifer - impact of the Third Rhône correction
    Le principal objectif de cette thèse est l’évaluation quantitative de l’impact de la troisième correction du Rhône (PR3) sur l’aquifère alluvial de la plaine du Rhône. Le PR3 est un projet d’aménagement, dont le but est d’une part l’amélioration de la protection contre les inondations, et d’autre part la renaturalisation du lit de la rivière Rhône (Valais, Suisse). L’approche considérée pour améliorer la capacité de drainage du Rhône est l’élargissement du lit de la rivière, ce qui causera un abaissement du niveau de la rivière (d’environ 0.9m en moyenne).
    L’impact du PR3 sur les eaux souterraines a été estimé avec un modèle numérique d’écoulement permanent. Ce modèle a été calibré avec des données de charge hydraulique. Le modèle numérique permet d’intégrer les données existantes, et d’étudier ensuite le comportement possible du système lorsqu’il sera soumis à des conditions qui diffèrent de celles investiguées actuellement (état restauré). Quand le modèle est utilisé pour simuler l’état futur du système, ses prévisions contiennent nécessairement un certain degré d’incertitude. Deux sources d’incertitude ont été considérées dans ce travail : d’une part la non-unicité de la solution du problème inverse, qui consiste à estimer les paramètres du modèle sur la base des mesures disponibles de charge hydraulique; d’autre part l’incertitude concernant l’évolution de la conductivité hydraulique du lit du Rhône suite à la restauration. L’incertitude sur les prévisions découlant de la non-unicité du problème inverse a été traitée grâce à la méthode du “null Space Monte Carlo” (NSMC) implémentée dans PEST. L’incertitude induite par l’état futur incertain du lit du Rhône a été traitée avec une analyse de sensibilité des prévisions par rapport à la conductivité du lit de la rivière.
    Le modèle est basé sur une compréhension conceptuelle du système, acquise grâce à l’interprétation des données disponibles. Ces données suggèrent clairement que le système est fortement contrôlé par des paramètres à l’échelle régionale, et que les interactions de l’eau souterraine avec le Rhône, ainsi qu’avec les canaux, joue un rôle clé dans le fonctionnement de l’aquifère.
    Le modèle est cohérent avec les données mesurées sur le terrain (charge hydraulique, zones d’échanges entre eau de surface et nappe phréatique, flux d’échange dans les canaux). Un résultat inattendu du modèle est que les zones d’échanges importants entre eaux superficielles et souterraines se produit sur de brefs tronçon du Rhône et à proximité des gravières et des carrières dans le lit de la rivière.
    En supposant que les propriétés du fond du lit du Rhône vont revenir à un état stable et identique à l’état actuel (ce qui implique le développement d’une couche de colmatation), le modèle prévoit que, à terme, le niveau d’eau s’abaisse à cause du PR3. Les résultats de modélisation suggèrent aussi que les canaux de drainage vont considérablement atténuer l’influence de la restauration du Rhône, et réduisent de plus l’incertitude sur les prévisions.
    L’incertitude des prévisions découlant de la non-unicité du problème inverse est faibles (une gamme de 0.07m en moyenne entre les quantiles q90 et q10), ce qui indiquent que les prévisions du modèle calibré semblent très fiables dans la majeure partie du domaine. Toutefois, cette faible incertitude doit être interprétée avec précaution : l’étude d’un exemple synthétique simple a montré que dans un modèle non-linéaire, la méthode NSMC peut sous-estimer l’incertitude. De plus, l’incertitude sur les paramètres découlant de la non-unicité du problème inverse, et calculée avec la méthode NSMC, est considérée comme faible en comparaison de l’évolution incertaine des propriétés du lit de la rivière dans l’état restauré. Dans la partie amont du modèle, il a été remarqué que les prévisions dépendent fortement de la conductivité hydraulique du lit de la rivière à un tel point que, en dépit de l’abaissement du niveau d’eau de la rivière, un léger accroissement de la conductivité du lit de la rivière produit une hausse du niveau piézométrique dans la nappe., The main objective of this thesis is the quantitative evaluation of the impact of the Third Rhône correction (PR3) on groundwater within the Rhône alluvial aquifer. PR3 is a river restoration project which aims at improving flood protection and at upgrading the ecologically deficient streambed of the Rhône River (Valais, Switzerland). The approach taken to extend the drainage capacity of the Rhône River is a widening of the streambed, which causes at large a lowering of river stages (of about 0.9 m on average).
    The impact of PR3 on groundwater was estimated using a steady state numerical flow model that was calibrated against average hydraulic head. The numerical model allows integrating data and can be used to simulate the system behavior under conditions which differ from the ones of the investigated state. When the model is used to simulate future system behavior, its predictions contain necessarily a certain degree of uncertainty. Two sources of predictive uncertainty were explored herein. On the one hand the non-uniqueness of the solution of the inverse problem, i.e. the estimation of model parameters based on available measurements of hydraulic head; on the other hand the uncertainty as to how the streambed conductance of the Rhône River will be altered as a consequence of the restoration. Predictive uncertainty arising from the non-uniqueness of the inverse problem was assessed by the “null space Monte Carlo” (NSMC) method implemented in PEST; predictive uncertainty associated with the uncertain evolution of the streambed conductance was estimated by a sensitivity analysis of model predictions to streambed conductance.
    The model is based upon a conceptual understanding of the system gained through interpretation of available data. Those data clearly suggests that the system is strongly driven by regional controls and that the interaction of groundwater with the Rhône River, as well as with ditches, plays a key role.
    The model is consistent with field measurements (head measurements, patterns of SW-GW interaction, exchange flow at ditches). One unexpected result of the model is that high SW-GW exchange flows occur on short stretches of the Rhône River close to gravel pits and at in-stream quarries.
    Assuming that streambed properties will return to a stable state identical to the current state (involving the development of a clogging layer), the model predicts that the groundwater table will lower at large in response to the PR3 in the long-term. Modeling results suggest also that drainage ditches considerably attenuate the influence of the restoration, and that furthermore they reduce predictive uncertainty.
    The predictive uncertainty arising from the non-uniqueness of the inverse problem is small (q90–q10 inter-percentile range of 0.07 m on average), indicating that the predictions made using the calibrated model seem to be very reliable in most parts of the model. However, this low uncertainty has to be interpreted with care: a simple synthetic example showed that in the context of a nonlinear model, the NSMC methodology potentially underestimates uncertainty. Moreover, parameter uncertainties arising from the non-uniqueness of the inverse problem are considered as low as compared to the uncertain evolution of streambed properties in the restored state. In the upstream part of the model, predictions were found to depend strongly on streambed conductance, to the point that despite the lowering of river stages, already a slight increase of the streambed conductance produces a rising of the groundwater table., Das Hauptziel dieser Doktorarbeit ist es, den Einfluss der 3. Rhonekorrektion (PR3) auf das Grundwasser des alluvialen Rhonetal-Aquifers quantitativ abzuschätzen. PR3 ist ein Projekt zur Flussrestaurierung der Rhone (Wallis, Schweiz) mit den Hauptzielen, den Hochwasserschutz zu verbessern und den Fluss ökologisch aufzuwerten. Das Projekt verfolgt den Ansatz, die Abflusskapazität des Flusses mittels einer Aufweitung des Flussbettes zu erhöhen; generell bewirkt dies eine Absenkung des Wasserspiegels der Rhone (um ca. 0.9 m im Mittel).
    Der Einfluss des PR3 auf das Grundwasser wurde mit Hilfe eines stationären numerischen Grundwasserfliessmodelles abgeschätzt, welches unter Verwendung von Mittelwerten des hydraulischen Potentials im Aquifer kalibriert wurde. Das numerische Modell eignet sich gut zur Integration vorhandener Daten und kann dazu verwendet werden, ein System unter veränderten Bedingungen zu untersuchen. Modellvoraussagen über das Verhalten des Systems in der Zukunft sind notwendigerweise mit einer gewissen Unsicherheit behaftet. In dieser Arbeit wurden zwei mögliche Quellen, aus denen sich Unsicherheiten bezüglich der Modellvoraussagen ergeben, berücksichtigt: einerseits die Unmöglichkeit, die Parameter des Fliessmodelles aus den vorhandenen Messungen des hydraulischen Potentials im Aquifer eindeutig zu bestimmen; andererseits die Ungewissheit darüber, wie die Durchlässigkeit der Rhonesohle sich im Zuge des PR3 verändert. Die Schätzunsicherheit, welche aus der fehlenden Eindeutigkeit des Inversen Problems entsteht, wurde mittels der „null-space-Monte-Carlo-Methode“ (NSMC) in PEST bemessen; die Unsicherheit, welche sich aus der Abhängigkeit der Modellvorhersagen von der unsicheren Entwicklung der Durchlässigkeit der Flusssohle ergibt, wurde mittels einer Sensitivitätsanalyse abgeschätzt.
    Das Modell basiert auf einem konzeptuellen Verständnis des Systems, welches durch die Interpretation der vorhandenen Daten gewonnen wurde. Diese Daten deuten eindeutig an, dass das System zu einem bedeutenden Teil durch Prozesse auf regionalem Massstab gesteuert wird und dass die Wechselwirkung von Grundwasser mit der Rhone sowie mit Kanälen eine Schlüsselrolle einnimmt für das Funktionieren des Aquifers.
    Das Modell ist konsistent mit Feldmessungen (Messungen des hydraulischen Potentials, Muster der Wechselwirkung zwischen Oberflächengewässer und Grundwasser, Austauschraten an Kanälen). Ein unerwartetes Modellresultat ist das Auftreten von hohen Austauschraten zwischen der Rhone und dem Grundwasser auf kurzen Abschnitten der Rhone in der Nähe von Baggerseen und an Stellen mit Kiesentnahmen aus dem Rhonebett.
    Unter der Voraussetzung, dass die Eigenschaften des Rhonebettes langfristig den heutigen entsprechen (was die Entwicklung einer Kolmationsschicht mit einschliesst), sagt das Modell als Folge des PR3 langfristig eine Absenkung des Grundwasserspiegels voraus. Die Modellresultate deuten auch an, dass die Drainagekanäle die Auswirkungen der Restaurierung merklich abmindern und dass deren Präsenz ausserdem die Unsicherheit der Modellvoraussagen reduziert.
    Die Unsicherheit der Voraussagen, die aus der Schätzunsicherheit der Parameter resultiert, scheinen die Voraussagen des kalibrierten Modelles an den meisten Stellen des Modells sehr verlässlich (durchschnittlicher x90–x10 Interperzentilbereich von 0.07 m). Diese geringe Unsicherheit muss jedoch mit Vorsicht gedeutet werden: ein einfaches synthetisches Beispiel verdeutlichte, dass die NSMC-Methode in nichtlinearen Modellen die Unsicherheit potentiell unterschätzt. Ausserdem wurde die mit der NSMC-Methode ermittelte Schätzunsicherheit als klein eingestuft im Vergleich zur Unsicherheit, die im Zusammenhang mit der ungewissen Entwicklung der hydraulischen Durchlässigkeit der Rhonesohle im restaurierten Zustand besteht. Es zeigte sich, dass die Modellvorhersagen im obstromseitigen Teil des Modells stark von der Durchlässigkeit der Rhonesohle abhängen; so sehr, dass in diesem Bereich schon eine geringfügige Zunahme der Durchlässigkeit der Rhonesohle trotz einer projektbedingten Absenkung des Wasserspiegels der Rhone einen Anstieg des Grundwassers nach sich zieht.
  • Publication
    Accès libre
    Hydraulic testing of low-permeability formations: A case study in the granite of Cadalso de los Vidrios, Spain
    (2009)
    Mejías, Miguel
    ;
    ;
    In recent years, the study of low-permeability geological formations has undergone considerable development. This is mainly due to the use of natural geological barriers to confine waste disposals, preventing leaking water from bringing contaminants into contact with the biosphere and the groundwater resources.
    In that context, the Spanish Geological Survey (IGME, Spanish acronym), supported by the National Radioactive Waste Management Agency (ENRESA, Spanish acronym) and with the technical advice of the Swedish Nuclear Fuel and Waste Management Company (SKB) designed and built a Mobile Hydrogeological Unit (UMH, Spanish acronym) for low-permeability formations characterization, which has been operational since 1997 and has been used for different purposes. Among other possibilities, the UMH allows conducting: constant transient-state flow-rate injection tests, constant-head transient-state injection tests, pseudo-stationary state injection tests, pressure fall-off tests, slug tests and pulse tests.
    The main objectives of this article are to describe the hydraulic characterization methodology used by IGME to carry out hydraulic tests in low permeability environment, to compare different testing methods and to summarize the results that have been obtained when characterizing the leucogranites of Cadalso de los Vidrios in Spain. The study area presents an increasing interest for granite production inside the Community of Madrid. The petrological and structural characteristics of the granite rocks and the core-samples extracted from a 200 m deep borehole investigation are described. The packer tests are conducted with the Mobile Hydrogeological Unit.
    The tests are interpreted with the help of analytical solutions. The main software used is Hytool, an open source matlab toolbox that provides a library of analytical solutions and a set of routines to facilitate hydraulic tests interpretation.
    The results allow the elaboration of a comparative analysis of the applied hydraulic tests and to define the hydraulic conductivity optimum application interval most suited for each of the used methods.
    Hydraulic conductivity values obtained varies from 3.20•10−7 m/s, for the upper weathered layer, to 2.80•10−12 m/s, for the test section from 97.36 to 116.15 m depth.
    Finally, the hydraulic conductivity values obtained in this area are compared with other case studies of granite formations around the world.
  • Publication
    Accès libre
    Understanding diagnostic plots for well-test interpretation
    (2009) ; ;
    Mejias, Miguel
    Dans le domaine de l’interprétation des essais de pompage, un graphe de diagnostic est un graphe simultané du rabattement et de sa dérivée logarithmique en fonction du temps. Ce graphe est habituellement représenté en échelle bi-logarithmique. Les avantages principaux et les limites de cet outil sont discutés à l’aide de trois exemples hydrogéologiques réels. Des recommandations pratiques sont proposées afin d’utiliser ces graphes pour sélectionner un modèle conceptuel d’écoulement. Lorsque les mesures de rabattement sont bruitées, le calcul de la dérivée logarithmique présente des artefacts que l’on peut minimiser en échantillonnant le signal. Lorsque le débit de pompage varie, ou lors d’un test de récupération, les graphes de diagnostic peuvent aussi être utilisé à condition d’effectuer une déconvolution des données au préalable. Enfin, les effets d’une erreur d’estimation du temps de début d’un test sont discutés. Tous ces exemples montrent que les graphes de diagnostic constituent un outil extrêmement utile pour l’interprétation des essais hydrauliques., In well-test analysis, a diagnostic plot is a scatter plot of both drawdown and its logarithmic derivative versus time. It is usually plotted in log–log scale. The main advantages and limitations of the method are reviewed with the help of three hydrogeological field examples. Guidelines are provided for the selection of an appropriate conceptual model from a qualitative analysis of the log-derivative. It is shown how the noise on the drawdown measurements is amplified by the calculation of the derivative and it is proposed to sample the signal in order to minimize this effect. When the discharge rates are varying, or when recovery data have to be interpreted, the diagnostic plot can be used, provided that the data are pre-processed by a deconvolution technique. The effect of time shift errors is also discussed. All these examples show that diagnostic plots have some limitations but they are extremely helpful because they provide a unified approach for well-test interpretation and are applicable in a wide range of situations., En los análisis de ensayos de pozos, un gráfico de diagnóstico es un gráfico de dispersión de las depresiones y de sus derivadas logarítmica en función del tiempo. Usualmente se grafica en una escala log – log. Las principales ventajas y limitaciones del método se revisaron con la ayuda de tres ejemplos hidrogeológicos de campo. Se proveen pautas para la selección de un modelo conceptual apropiado a partir de un análisis cualitativo de la derivada logarítmica. Se muestra cómo el ruido en las mediciones de las depresiones resulta amplificado por el cálculo de las derivadas y se propone muestrear la señal con el objeto de minimizar este efecto. Cuando los ritmos de descarga varían, o cuando los datos de recuperación deben ser interpretados, el gráfico diagnóstico puede ser utilizado siempre que los datos sean pre-procesados por una técnica de deconvolución. Asimismo se discute el efecto de los errores de desplazamiento del tiempo. Todos estos ejemplos muestran que el gráfico diagnóstico tiene algunas limitaciones pero son extremadamente útiles porque proveen un enfoque unificado para la interpretación de ensayos de pozos y son aplicables en una amplia gama de situaciones., Nel dominio dell’interpretazione delle prove di pompaggio, un grafico diagnostico è un grafico simultaneo dell’abbassamento e della sua derivata logaritmica in funzione del tempo. Questo grafico è solitamente rappresentato in scala bi-logaritmica. I principali vantaggi ed i limiti di questo strumento sono discussi attraverso tre esempi idrogeologici reali. Delle raccomandazioni pratiche sono proposte al fine di utilizzare questi grafici per scegliere un modello concettuale di flusso. Quando le misure di abbassamento contengono del rumore, il calcolo della derivata logaritmica presenta degli artefatti che sono riducibili campionando il segnale. Quando la portata del pompaggio varia, o durante un test di recupero, i grafici diagnostici possono anche essere utilizzati a condizione di effettuare una preliminare deconvoluzione dei dati. Infine, gli effetti di un errore di stima del tempo di inizio di un test sono discussi. Tutti gli esempi mostrano che i grafici di diagnosi costituiscono uno strumento estremamente utile per l’interpretazione delle prove idrauliche.