1,721,167 research outputs found
Development of a 3D groundwater model based on scarce data: the Wadi Kafrein catchment/Jordan
Full text linkIn this work, a 3D groundwater flow model integrating all important geological features of the hydrogeological system is developed to investigate hydrological processes in the Wadi Kafrein area of Jordan. A large amount of available geological and hydrological data is integrated to construct a 3D groundwater flow model for the Wadi Kafrein area. Using the newly developed mapping approach, the translation of the highly detailed geological formations to an unstructured finite element grids, can be accomplished with high precision. The existing data set for model calibration is scarce, which is a typical situation for many hydrogeological case studies. At first, the steady state calibration of the groundwater model is carried out based on the observation wells. Then, the time and space-dependent recharge from precipitation are applied at the top surface of the finite element model. The transient simulation is conducted during the period of 1996-2008 considering the abstraction rates of the production wells and discharge of the springs. The calculated water levels are close to the observed values. The difference is partly caused by return flows from irrigation and the groundwater inflow from the adjacent aquifers which are not taken into consideration so far. Since the Wadi Kafrein area is an important agricultural area in the semiarid region of the Lower Jordan Valley, the model developed in this study can be regarded as a useful tool for analyzing the hydrological processes and improving groundwater management practices elsewhere affected by similar geological and hydrogeological conditions
Beiträge zur rechnerischen Geotechnik : nicht-isothermischer Fluss in niedrig-porösen Medien
No full textTo better understand the coupling of thermal (T), hydraulic (H) and mechanical (M) processes (T-H-M processes) and their influence on the system behaviour, models allowing T-H-M coupling are developed. These models allow simulations in the near-field of the system. The modeling of non-isothermal thermo-hydraulic (TH) processes is important for applications such as geothermal energy generation, heat supported environmental remediation, and nuclear waste disposal. The work presented herein focuses on deep geological disposal of nuclear waste, and more specifically on the thermal and hydraulic processes in this application.
Thermal processes result directly from the heat radiation of the waste and include heat transport from the core to the bentonite buffer. Other processes of importance are vaporization and condensation associated with phase changes between the liquid and gaseous phases. Hydraulic processes of importance include water intrusion from the host rock to the buffer and eventually
to the core, as well as swelling and shrinking processes in the bentonite. Bentonite swells as a result of water intrusion from the host rock and dries as a result of the heat transport from the core.
The objective of the work is to formulate the processes mathematically and to integrate them into the object-oriented simulator GeoSys/RockFlow.
Buffer, host rock, and fluids in the gas and liquid phase form a multiphase-multicomponental system (porous medium). The TH model consists of a set of three balance equations. One balance equation for the water component, one balance equation for the air component and one energy balance equation. The three primary, or independent variables are gas pressure, water saturation, and temperature.
To solve these balance equations, equations describing the material modelled are necessary. Material properties include for example capillary pressure-saturation relationships, density equations, or viscosity calculations. For those processes, material parameters and state variables are highly non-linear and mostly functions of temperature, saturation, and pressure. Other than describing the material, the thermodynamic state of the system has to be described. This is achieved with equations of state, as for example functions for the calculation of liquid
pressure or mass fractions. When the material properties and the state functions are inserted into the balance equations, governing equations in the di.erential form are obtained. After numerical transformations, these equations are then solved by GeoSys/RockFlow.
The implementation allows phase changes between the fluid phases (gas and liquid) to occur explicitly. The model allows the simulation of processes in very low permeability clays with high capillary pressures. Examples for code validation are shown, where low permeability clay is desaturated.Um thermische (T), hydraulische (H) und mechanische (M) Prozesse, ihre Kopplungen und ihren Einfluss auf das Systemverhalten besser zu verstehen, werden T-H-M Modelle entwickelt. Diese Modelle lassen Modellierungen im Nahbereich des Systems zu. Die Simulierung von nicht-isothermen, thermo-hydraulischen (TH) Prozessen ist für Anwendungen wie Geothermie, Wäremeunterstützte Grundwassersanierung, und Atommüllentsorgung wichtig. Diese Doktorarbeit richtet sich spezifisch an den Anwendungen für Atommüllentsorgung aus, vor Allem an den thermischen und hydraulischen Prozessen in diesem Bereich.
Thermische Prozesse kommen durch die Wärmestrahlung des Abfalls zustande und zeichnen sich durch den Wärmetransport vom Kanister zum Bentonitpuffer aus. Andere wichtige thermische Prozesse sind Verdampfung und Kondensierung, die mit dem Phasenwechsel zwischen wässeriger- und gas Phase assoziiert sind. Wichtige hydraulische Prozesse sind das Eindringen von Wasser von dem Gastgestein in den Puffer und dann in den Kanister, sowie quellen und schrumpfen des Bentonitpuffers. Bentonit quillt durch das Eindringen von Wasser und trocknet durch den Wärmetransport des Abfalls.
Das Ziel der Arbeit ist die mathematische Formulierung der Prozesse und ihre Integrierung in den objekt-orientierten finite-elemente Code GeoSys/RockFlow. Puffer, Gastgestein und Fluide (in der wässerigen und gas Phase) formen gemeinsam ein mehrphasiges-mehrkomponenten System (poröses Medium). Das TH Modell benötigt drei Bilanzgleichungen, eine für die Wasserkomponente, eine für die Luftkomponente und eine für Energie. Die drei Primärvariablen sind Gasdruck, Wassersättigung und Temperatur. Um die Bilanzgleichungen zu lösen, werden Gleichungen, die das System beschreiben, benötigt. Dies sind Materialparameter, wie zum Beispiel Kapillardruck - Sättigungsbeziehungen oder Dichtegleichungen. Für die beschriebenen Prozesse sind diese Materialparameter und Zustandsgleichungen meist nicht-linear und meist Funktionen der Temperatur, der Sättigung und des Drucks. Nicht nur das Material, sondern auch der thermodynamische Zustand des Systems muß beschrieben werden. Dies wird mit Zustandsgleichungen erreicht, zum Beispiel Funktionen zur Berechnung des Wasserdrucks, oder der Massenfraktionen. Materialparameter und Zustandsgleichungen werden in die Bilanzgleichungen substituiert, es resultieren die drei Modellgleichungen in Differentialform. Diese Gleichungen werden nach Umformungen von GeoSys/RockFlow gelöst.
Die Implementierung läßt Phasenwechsel zwischen den Fluidphasen (wasser und gas) explizit zu. Das Modell ermöglicht Simulationen von sehr undurchlässigem Tonmaterial mit hohen Kapillardrücken. Beispiele der Modellvalidierung werden gezeigt, wo Ton durch hohe Temperaturen entsättigt wird
Experimentelle und numerische Untersuchungen zum Materialverhalten Geotechnischer Barrieren
No full textQuellfähige Tone spielen eine bedeutende Rolle bei aktuellen Konzepten zur Endlagerung hoch-radioaktiver Abfälle in kristallinem Wirtsgestein. Zum Schutz vor dem Austritt radioaktiver Substanzen wurden Multi-Barrieren-Konzepte entwickelt. Die Barrieren bestehen aus einem Kupferkanister, aus kompaktiertem Bentonit als Verfüllmaterial (die Geotechnische Barriere) sowie dem Wirtsgestein. Durch Korrosion der Kupferkanister oder Radiolyse kann bei hoch-radioaktiven Abfällen Wasserstoff entstehen. Durch die in diesem Konzept vorgesehene Aufstättigung des Verfüllmaterials mit Kluftwasser wird seine Permeabilität stark reduziert, so dass sich Wasserstoffgas in dem Raum zwischen Kanister und der Geotechnischen Barierre aufstauen kann. Durch diese Gasansammlung kann der Eindringdruck des Verfüllmaterials überschritten werden, und es kommt dann zu einer Migration der Gase durch die Geotechnische Barriere.
Zur Untersuchung der thermischen und hydraulischen Eigenschaften des Verfüllmaterials unter Bedingungen, die in untertätigen Endlagern für radioaktiven Abfall herrschen, wurden verschiedene Laborexperimente konzipiert. Alle Experimente wurde mit Natrium- (SPV Volclay) oder Kalzium-Bentonit (Calcigel), der mit Gesteinsmehl vermischt wurde, unter Verwendung von Wässern aus dem Felslabor Äspö durchgeführt. Kapillardruckkurven wurden für alle Sättigungsbereiche mit einer Kapillardruckmessdose und einem Thermohygrometer bestimmt. Hydraulische Säulenexperimente wurden mit einem speziell angefertigten Permeameter durchgeführt. Unter Anwendung der Darcy-Gleichung konnte die Permeabilität berechnet werden. Die thermischen und nicht-isothermalen Drainage Experimente wurden mit Hilfe inverser Modellierung mit dem Computerprogramm iTOUGH2 ausgewertet. Die thermischen Experimente wurden bis zur konstanten Temperaturverteilung durchgeführt. Für alle Zeitpunkte und an allen Messpunkten entsprach die simulierte Temperaturverteilung sehr gut den Messdaten. Die invers geschätzte gesättigte thermische Leitfähigkeit und gesättigte spezifische Wärmekapazität waren konsistent mit Ergebnissen aus empirischen Gleichungen. Anhand von nicht-isothermalen Drainage Experimenten wurden Drücke, Temperaturen und die drainierte Wassermenge gemessen. Mit Hilfe des automatisierten Kalibrierungsverfahrens im iTOUGH2 wurden die gesättigte Permeabilität, die thermische Leitfähigkeit, spezifische Wärmekapazität und der Gaseindringdruck der Bentonit/Gesteinsmehl-Mischung geschätzt. Diese Inversion ergab den gleichen Parametersatz zur Beschreibung der Kapillardruck–Sättigungsbeziehung wie die Messung mit der Kapillardruckmessdose und dem Thermohygrometer. Das nicht-isothermale Drainage Experiment diente des weitern zum Codevergleich zwischen TOUGH2 und RockFlow/RockMech sowie zur Validierung von RockFlow/RockMech.
Anhand unterschiedlichster experimenteller, numerischer und analytischer Verfahren konnten konsistente Materialparameter für das Verfüllmaterial bestimmt werden. Diese Materialparameter können nun als Basis zur Simulation weiterer thermisch-hydraulisch gekoppelter Prozesse in untertägigen Deponien für radioaktive Abfälle dienen. Die neu entwickelten Experimente in Kombination mit inverser Modellierung erlauben somit die Bestimmung von Materialparametern die für thermisch-hydraulische Prozesse unter Endlagerbedingungen von ausschlaggebender Bedeutung sind.Swelling clays play a major role in current concepts for the underground disposal of high-level nuclear waste in deep geological formations. In one of the multi-barrier concepts for preventing the escape of radioactive substances from a high-level nuclear waste repository, the barrier consists of a copper container, compacted bentonite as buffer and backfill (the geotechnical barrier), and the repository host rock. Corrosion of the copper canister and radiolysis both produce hydrogen. When the buffer and backfill are saturated with water and the permeability of the bentonite is reduced by swelling, any hydrogen that is produced can accumulate in the space between the container and the geotechnical barrier. This will result in pressures exceeding the entry pressure of the buffer and backfill, and passage of gas through the geotechnical barrier.
An experimental program was developed to investigate the thermal and hydraulic properties of the buffer and backfill under conditions expected to exist in a permanent repository for radioactive waste. All experiments were conducted with mixtures containing sodium- (SPV Volclay) or calcium-bentonite (Calcigel) and crushed rock, and we used water from the Äspö test site. Water retention curves were measured from low to high saturation using a pressure cell and a thermohygrometer. Hydraulic column experiments were carried out with a specially designed permeameter and Darcy’s law was applied to determine the permeability. The thermal and non-isothermal drainage laboratory experiments were analyzed with inverse modeling techniques using iTOUGH2. The thermal experiments were conducted until a time-invariant temperature distribution was reached. The simulated temperature distribution matched the measured data very well at all locations along the column and for all times. The inversely estimated thermal conductivity and specific heat were consistent with the predictions of the empirical relationships. Pressure, temperature and the drained water volume were measured with non-isothermal drainage experiments and jointly inverted to estimate absolute permeability, thermal conductivity, specific heat, and capillary strength parameters of the bentonite/crushed rock mixtures. Consistent capillary pressure curves were obtained with the inversion of transient data and the direct pressure cell method. One of the non-isothermal drainage experiments was treated as a benchmark between TOUGH2 and RockFlow/RockMech. The simulation of the experiment served also a validation for RockFlow/RockMech.
The parameters estimated using different experimental, numerical, and analytical procedures were consistent with one another, providing backfill material properties useful for the simulation of gas- and heat-generating nuclear waste repositories. The newly developed experimental setup in combination with inverse modeling allows the identification of key parameters governing hydraulic and thermal processes under repository conditions
Massentransport in Geklüfteten Medien - Der Übergang zum Anormalen Transport
No full textSubstantial understanding of the processes involved in fluid flow and tracer transport in geological systems can be gained from the use of numerical simulations. Regardless of the nature or scale of a geological system, the probability of encountering fractures in hard rock units is high. Therefore, modelling flow and transport in fractures and fracture networks is a key component of almost any complex hydrogeological simulation. Over the last few decades, with the increase in computational power and the availability of increasingly sophisticated and powerful specialized modelling software, numerical models have become increasingly commonplace, and are meanwhile regarded as vital tools in the wide field of geoscience.
This work portrays the development of a series of numerical models, with discrete fracture network geometry, based on the knowledge of an actual fracture network gained through prior experiments performed at the Grimsel Test Site (GTS) in the Swiss Alps. A three dimensional reconstruction of the actual fracture network geometry is developed from the extensive photographic data set provided by the Excavation Project (EP) and used as the basis for purely advective particle tracking models, as well as for fully advective-dispersive transport models. Realistic fracture apertures based on actual fracture measurements are generated through a geostatistical method and their influence upon flow and transport is studied. Various degrees of heterogeneity with regard to aperture distributions are generated and used in a Monte Carlo approach to examine their influence. Further, the retardation effects of matrix diffusion are also considered and incorporated into the models by means of a new semi-analytical method included in the FEM code Rockflow/Geosys V4. The resulting breakthrough curves are analysed and fitted with several analytical solutions, including advective dispersive transport with matrix diffusion and continuous time random walks (CTRW). It is shown that the standard advection dispersion equation (ADE) is not suitable for approximating or predicting breakthrough curves from heterogeneous models.
The development of the various models portrayed herein depended strongly on the ongoing advances made in the available modelling codes and in new and dedicated preprocessing tools. Several newly developed methods and tools were tested and applied during the project.
For the first time, a group of numerical simulations is presented, that incorporates complex fracture network geometry based on actual field data, geostatistically generated realistic apertures and the effect of matrix diffusion. This new level of realism provides the basis for new insights into the role of heterogeneity, matrix diffusion and the combination of both in the transition toward anomalous transport.Das Wissen über die an Fluss- und Transportprozessen beteiligten Mechanismen kann durch die Verwendung von numerischen Simulationen deutlich erweitert werden. Unabhängig der Natur oder der Grössenordnung des untersuchten Systems ist das Auftreten von Klüften und Spalten in Felsstrukturen sehr wahrscheinlich. Daher stellt das Modellieren von Fluss- und Transport in geklüfteten Medien eine wesentliche Komponente der meisten komplexen hydrogeologischen Simulationen dar. Mit der andauernden Entwicklung von immer Leistungsfähigeren Computern und dem Aufkommen von spezialisierter und leistungsfähiger Software zur numerischen Modellierung haben sich numerische Simulationen in den letzten Jahren und Jahrzehnten etabliert, und sind nun als wichtiges Werkzeug im weiten Feld der Geowissenschaften anerkannt.
Diese Arbeit stellt die Entwicklung und Anwendung einer Reihe solcher numerischer Modelle dar. Die Modelle basieren auf dem bestehenden Wissen über ein natürliches Kluftsystem in den Schweizer Alpen, welches durch frühere Experimente am Felslabor Grimsel untersucht wurde. Eine dreidimensionale Rekonstruktion der Kluftgeometrie wird entwickelt und als Basis für sowohl rein advective "Particle Tracking" Modelle, als auch für advektiv-dispersive Transportmodelle verwendet. Realistische Kluftöffnungsweiten, basierend auf gemessenen Weiten von untersuchten Klüften, werden mittels einer Geostatistischen Methode für die Modelle generiert. Der Einfluss dieser Heterogenität auf Fluss und Transport wird untersucht. Vergleichbare Öffnungsweitenverteilungen mit unterschiedlich starker Heterogenität werden generiert und deren Einfluss im Rahmen eines Monte Carlo Ansatzes Untersucht. Ebenso wird der retardierende Einfluss von Matrix Diffusion in den Modellen berücksichtigt und untersucht. Die resultierenden Durchbruchskurven werden mittels einer Auswahl an analytischen Methoden, unter anderem mit Matrix Diffusion und "Continuous Time Random Walks" (CTRW), angepasst und interpretiert. Es wird dabei gezeigt, dass die herkömliche Advektions-Dispersions Gleichung (ADE) nicht geeignet ist um Durchbrüche von heterogenen Modellen anzunähern oder vorauszusagen.
Die Entwicklung der vorgestellten Modelle wird ermöglicht durch die fortwährende Entwicklung im Bereich von numerischen Modellierungscodes, sowie durch neue und verbesserte Werkzeugen zur Erstellung und Verwaltung der Modellgeometrie. Etliche neue Methoden und Werkzeuge wurden im Rahmen der Arbeit getestet und verwendet.
Erstmals wird eine Gruppe von numerischen Simulationen vorgestellt, welche complexe, auf Felddaten beruhende Kluftnetzwerkgeometrien beinhalten, sowie geostatistisch generierte Öffnungsweiten und Matrix Diffusion. Dieser neue Grad an Realitätsnähe bildet die Basis für neue Erkenntnisse bezüglich der Rolle von Heterogenitäten, Matrix Diffusion und deren Kombination im Übergangsbereich zum anormalen Transport
Hydrogen energy storage in porous geological formations: Investigation of storage dimensioning, induced effects and monitoring methods
Full text linkPower production from renewable sources is fluctuating, governed by the availability of the sources such as wind or solar radiation. Increasing the share of renewable sources can therefore lead to shortages in power production. Large-scale energy storage in porous geological formations, e.g. using hydrogen (H2), produced during periods of ample power supply, could help mitigate offsets between supply and demand. In this thesis, hydrogen energy storage in porous geological formations is investigated using numerical scenario simulations with emphases on the quantification of storage dimensions, on induced hydraulic and thermal effects and on geophysical monitoring techniques. For this end, a coupling approach for the scientific open source simulator OpenGeoSys and the proprietary reservoir simulator ECLIPSE is developed. The coupled simulator combines the multiphase-multicomponent flow and transport of ECLIPSE with the heat and geochemical reactive transport of OpenGeoSys. Furthermore, the Joule-Thomson effect is implemented in OpenGeoSys for gases consisting of N2, O2, H2, CO2 and CH4. For the assessment of potential storage dimensions and induced hydraulic and thermal effects, a hypothetical storage operation is simulated at an existing structure in northern Germany. 25 realistic and heterogeneous parameter distributions are generated using a facies modelling approach. Additionally, the applicability of using seismic waveform inversion (FWI), electrical resistivity tomography (ERT) and gravity monitoring methods are investigated for H2 storage operations.Die Stromproduktion aus erneuerbaren Energien unterliegt Schwankungen bedingt durch die variable Verfügbarkeit der Energieträger. Durch einen hohen Anteil von erneuerbaren Energien an der Gesamtproduktion können so Engpässe in der Stromproduktion entstehen, welche durch Energiespeicher ausgeglichen werden müssen. Hierfür kann Wasserstoff, mit Hilfe von Überkapazitäten synthetisiert, als chemischer Energieträger innerhalb von porösen geologischen Formationen gespeichert werden. In der vorliegenden Arbeit werden potentielle Dimensionen einer Wasserstoffspeicherung, die induzierten hydraulischen und thermischen Effekte sowie die Anwendung von geophysikalischen Monitoringmethoden mit Hilfe von numerischen Szenariensimulation untersucht. Hierfür wurde zunächst eine Modellkopplung für den wissenschaftlichen Open-Source Simulator OpenGeoSys und den proprietären Reservoir-Simulator ECLIPSE (© Schlumberger) entwickelt. Der gekoppelte Simulator kombiniert die Mehrphasen-Mehrkomponenten Strömungssimulation von ECLIPSE mit dem Wärme- und Reaktiven-Transport-Modell von OpenGeoSys. Zusätzlich wurde das Wärmetransportmodell von OpenGeoSys durch den Joule-Thomson-Effekt für Gase, bestehend aus N2, O2, H2, CO2 und CH4, erweitert. Für die Abschätzung von möglichen Speicherdimensionen sowie der induzierten hydraulischen und thermischen Effekte wurde ein hypothetischer Speicherbetrieb, basierend auf einer existierenden geologischen Struktur in Norddeutschland, simuliert. Um die Variabilität der Verteilung geologischer Parameter abzubilden, wurden mit Hilfe eines Faziesmodells 25 realistische und gleich wahrscheinliche Parameterverteilungen erstellt. Des Weiteren wird in der vorliegenden Arbeit die Anwendbarkeit von seismischer Wellenform Inversion (FWI), Geoelektrik (ERT) und Gravimetrie für ein Monitoring eines H2-Speichers untersucht
Impact of aquifer heterogeneity on subsurface flow and salt transport at different scales: from a method determine parameters of heterogeneous permeability at local scale to a large-scale model for the sedimentary basin of Thuringia
Full text linkMore than 70% of the German population rely on groundwater as their daily drinking water supply. Hence, the understanding of processes like fluid flow and solute transport in porous media are of broad relevance. However, subsurface hydrology faces two problems. On the one hand, information about the subsurface structure and hydraulic properties is scarce. On the other hand, many aquifer parameters like the permeability which controls the flow velocity, are heterogeneously distributed in space. This thesis is dedicated to the question of how aquifer heterogeneity impacts on subsurface flow and transport processes at different scales. We first present a large-scale numerical model of the Thuringian basin in order to investigate the mechanisms of saltwater transport. The results underline the fact that permeability and its heterogeneous spatial distribution are decisive factors for the evolving flow and salt patterns. Our results show that the correlation structure and the degree of heterogeneity (variance) impact significantly on the amount of dissolved salt and the location where salt reaches near-surface regions. Thermally induced convection is not present. However, variations in fluid density due to dissolved salt can lead to significant changes in the distribution and amount of salt in the Thuringian Basin. We further focus on a method to determine the statistical parameters which describe the heterogeneous spatial distribution of permeability by analyzing well flow at local scale. This is of interest, because pumping tests are a widely used tool to infer on porous medium characteristics. We derive an analytical solution which describes the drawdown of a steady state pumping test in 3D heterogeneous anisotropic media by making use of upscaling theory. By combining the analytical solution with an inverse estimation strategy we determine the statistical parameters, like mean permeability, variance and correlation length from on-site pumping tests
Compressed air energy storage in porous geological formations: Investigation of storage characteristics and induced impacts
Full text linkWith the rapid growth of energy production from intermittent renewable sources like wind and solar power plants, large-scale energy storage options are required to compensate for fluctuating power generation on different time scales. Compressed air energy storage (CAES) in porous geological formations is seen as a promising underground storage option for balancing short-term diurnal fluctuations. CAES is a power-to-power energy storage, which converts electricity to mechanical energy, i.e. highly pressurized air, and stores it in the subsurface. This thesis aims at investigating the feasibility, dimensioning the storage capacity, and assessing the induced hydraulic, thermal and geochemical impacts of a large-scale CAES operation in porous formations. For this, scenario based numerical simulations are used
Experimental procedure for studying the degradation and alteration of limestone slabs applied on exterior cladding
Full text linkThe studied limestone is a well known Portuguese natural stone that occurs in Valverde—Alcanede region, located in the
Maciço Calcário Estremenho (center of Portugal). This stone is used in several exterior and interior applications such as
paving, cladding, masonry and decorations. Slabs made of the studied limestone were installed in ventilated facade with a
“kerf” anchoring system in a building located in Valencia (Spain). After approximately fve years, signs of degradation were
detected on the facades through changes in color, enhancment of cracks and fssures that caused instability and collapse on
several slabs. Limestone slabs comprise irregular patterns due to rock cutting across stylolites. These features represent ideal
paths for fuid circulation through the slab and can be identifed as vulnerability zones for exterior applications. Experimental
procedure included several laboratory analyses to study clays and the iron oxide contents. Microscopic petrography, XRD and
SEM observations were important to identify the interaction of the clay material in stylolites and fossil contours. Results show
the importance of establishing standard methods for selecting stone on cladding applications. From the results, it is possible
to understand that clay minerals have a strong infuence on the mechanical performance of this stone. Among other important
remarks, results show the importance of the development of standard procedures that acknowledge the structure and mineral
composition before setting these natural products as construction materials. Proper stone selection was found important to
avoid facade degradation, and hence contribute to prevent accidents and promote user safety as well as economic impacts.V. Pires gratefully acknowledges the support of Fundação
para a Ciência e Tecnologia (FCT) under the project UIDB/04449/2020
e UIDP/04449/202
- …
