1,721,026 research outputs found
BlrInflowNetwork_Manohar and Mohan Kumar 2013.inp
No full textBangalore Inflow Network Data; used in : Manohar, U., & Mohan Kumar, M. S. (2013). Modeling equitable distribution of water: Dynamic inversion-based controller approach. Journal of Water Resources Planning and Management, 140(5), 607-619.</div
Going Beyond Counting First Authors in Author Co-citation Analysis
Get PDFThe present study examines one of the fundamental aspects of author co-citation analysis (ACA) - the way co-citation
counts are defined. Co-citation counting provides the data on which all subsequent statistical analyses and mappings
are based, and we compare ACA results based on two different types of co-citation counting - the traditional type that
only counts the first one among a cited work's authors on the one hand and a non-traditional type that takes into
account the first 5 authors of a cited work on the other hand. Results indicate that the picture produced through this non-traditional author co-citation counting contains more coherent author groups and is therefore considerably clearer. However, this picture represents fewer specialties in the research field being studied than that produced through the traditional first-author co-citation counting when the same number of top-ranked authors is selected and analyzed. Reasons for these effects are discussed
Variations on the Author
Get PDF“Variations on the Author” discusses two of Eduardo Coutinho’s recent films (Um Dia na Vida, from 2010, and Últimas Conversas, posthumously released in 2015) and their contribution to the general question of documentary authorship. The director’s filmography is characterized by a consistent yet self-effacing form of authorial self-inscription: Coutinho often features as an interviewer that rather than express opinions propels discourses; an interviewer that is good at listening. This mode of self-inscription characterizes him as an author who is not expressive but who is nonetheless markedly present on the screen. In Um Dia na Vida, however, Coutinho is completely absent form the image, while Últimas Conversas, on the contrary, includes a confessional prologue that moves the director from the margins to the center of his films. This article examines the ways in which these works stand out in the filmography of a director who offers new insights into the notion of cinematic authorship
Appropriate Similarity Measures for Author Cocitation Analysis
Get PDFWe provide a number of new insights into the methodological discussion about author cocitation analysis. We first argue that the use of the Pearson correlation for measuring the similarity between authors’ cocitation profiles is not very satisfactory. We then discuss what kind of similarity measures may be used as an alternative to the Pearson correlation. We consider three similarity measures in particular. One is the well-known cosine. The other two similarity measures have not been used before in the bibliometric literature. Finally, we show by means of an example that our findings have a high practical relevance.information science;Pearson correlation;cosine;similarity measure;author cocitation analysis
Water Balance Studies In A Small Experimental Forested Watershed, South India
No full textForested watersheds play a dominant role in the global hydrological cycle. Very few experimental observatories especially in tropical forested regions of India have been undertaken. This study has been initiated for this reason and to gain insights into functioning of the hydrological system in such climatic conditions. This study involves experimental setup of a watershed, it’s monitoring till date, modelling of the hydrological processes observed and the challenges in modelling components of the water balance in this watershed.
A Small Experimental Watershed of 4.3 Km2 was set up at Mule Hole, in South India along the Kerala-Karnataka State borders, and is situated inside the Bandipur National park. After an overview of watershed studies, review of literature related to forest watershed studies and processes in the first two chapters, Chapter 3 introduces the study area, Mule Hole Experimental Watershed and explains the methodology used to study this watershed. Model SWAT was used initially to simulate the water balance components. A brief description of the model, methodology adopted and discussion on the results obtained is presented in Chapter 4. The watershed initially modelled as an ungauged watershed using the default parameters in the model, simulated very high groundwater contribution to the runoff. The calibrated model although performed favourably for annual average values and monthly calibration, the daily calibration was unsatisfactory. An auxiliary study on quantification of actual and potential evapotranspiration (ET0) has been carried out in Chapter 5 . Ten methods including Penman-Montieth were compared and evaluated for efficacy of the methods. All methods except for Hargreaves method showed agreement with the Penman-Montieth for annual average values. Priestly-Taylor method was found be the best estimator in comparison with Penman-Montieth method, when used to estimate AET. Adjusted Hargreaves and FAO Blaney -Criddle method were found to be very useful when few or limited climatic data were available for estimation of Potential evapotranspiration. A multidisciplinary approach of estimating recharge consisting of chloride mass balance technique coupled with study of water table fluctuations and groundwater flow analytical modelling has been attempted in Chapter 6. Direct and localized recharge was estimated at 45 mm/yr and indirect recharge 30 mm/yr for the monitored years in the watershed. The low values of recharge rates implied an unexpected very high evapotranspiration rate. It may be inferred that in the absence of groundwater flow to the stream, the recharge joins groundwater flow as outflow of the hydrologic system. An integrated lumped model incorporating the regolith zone and the capability of the tree roots to access this store is presented in Chapter 7. The model was able to simulate the pattern of lag-time between water table rise was observed in shallow piezometers in comparison with hillslope piezometers. The patterns of water table variation among the different hillslope piezometers suggest that they are linked with local processes and not by a regional aquifer dynamics. This study shows that water uptake, combined with the spatial variability of regolith depth, can account for the variable lag time between drainage events and groundwater rise observed for the different piezometers. Chapter 8 discusses the results, conclusions derived from this study and possibility of further scope of studies
Multi-scale Modeling of Nanoparticle Transport in Porous Media : Pore Scale to Darcy Scale
Get PDFAccurate prediction of colloid deposition rates in porous media is essential in several applications. These include natural filtration of pathogenic microorganisms such as bacteria, viruses, and protozoa, transport and fate of colloid-associated transport of contaminants, deep bed and river bank filtration for water treatment, fate and transport of engineered nanoparticles released into the environment, and bioremediation of contaminated sites. Colloid transport in porous media is a multi-scale problem, with length scales spanning from the sub-pore scale, where the particle-soil interaction forces control the deposition, up to the Darcy scale, where the macroscopic equations governing particle transport are formulated. Colloid retention at the Darcy scale is due to the lumped effect of processes occurring at the pore scale. This requires the incorporation of the micro-scale physics into macroscopic models for a better understanding of colloid deposition in porous media. That can be achieved through pore-scale modeling and the subsequent upscaling to the Darcy scale. Colloid Filtration Theory (CFT), the most commonly used approach to describe colloid attachment onto the soil grains in the subsurface, is found to accurately predict the deposition rates of micron-sized particles under favorable conditions for deposition. But, CFT has been found to over predict particle deposition rates at low flow velocity conditions, typical of groundwater flow, and for nanoscale particles. Also, CFT is found to be inapplicable at typical environmental conditions, where conditions become unfavorable for deposition, due to factors not considered in CFT such as deposition in the secondary minimum of the interaction energy profile, grain surface roughness, surface charge heterogeneity of grains and colloids, and deposition at grain-to-grain contacts. To the best of our knowledge, mechanistic-based models for predicting colloid deposition rates under unfavorable conditions do not exist. Currently, fitting the colloid breakthrough curve (BTC), obtained from the laboratory column-or field-scale experiments, to the advection-dispersion-deposition model is used to estimate the values of deposition rate coefficients. Because of their small size (less than 100 nm), nanoparticles, a sub-class of colloids, may interact with the porous medium in a different way as compared to the larger colloids, resulting in different retention mechanisms for nanoparticles and micron-sized particles. This emphasizes the need to study nanoparticles separately from larger, micrometer-sized colloids to better understand nanoparticle retention mechanisms.
The work reported in this thesis contributes towards developing mathematical models to predict nanoparticle movement in porous media. A comprehensive mechanistic approach is employed by integrating pore-scale processes into Darcy-scale models through pore-network modeling to upscale nanoparticle transport in saturated porous media to the Darcy scale, and to develop correlation equations for the Darcy-scale deposition parameters in terms of various measurable parameters at Darcy scale. Further, a one-dimensional mathematical model to simulate the co-transport of viruses and colloids in partially saturated porous media is developed to understand the relative importance of various interactions on virus transport in porous media.
Pore-network modeling offers a valuable upscaling tool to express the macroscopic behavior by accounting for the relevant physics at the underlying pore scale. This is done by idealizing the pore space as an interconnected network of pore elements of different sizes and variably connected to each other, and simulating flow and transport through the network of pores, with the relevant physics implemented on a pore to pore basis (Raoof, 2011). By comparing the results of pore-network modeling with an appropriate mathematical model describing the macro-scale behavior, a relationship between the properties at the macro scale and those at the pore scale can be obtained. A three dimensional multi-directional pore-network model, PoreFlow, developed by Raoof et al. (2010, 2013) is employed in this thesis, which represents the porous medium as an interconnected network of cylindrical pore throats and spherical pore bodies, to upscale nanoparticle transport from pore scale to the Darcy scale. The first step in this procedure is to obtain relationships between adsorbed mass and aqueous mass for a single pore. A mathematical model is developed to simulate nanoparticle transport in a saturated cylindrical pore by solving the full transport equation, considering various processes such as advection, diffusion, hydrodynamic wall effects, and nanoparticle-collector surface interactions. The pore space is divided into three different regions: bulk, diffusion and potential regions, based on the dominant processes acting in each of these regions. In both bulk and diffusion regions, nanoparticle transport is governed by advection and diffusion. However, in the diffusion region, the diffusion is significantly reduced due to hydrodynamic wall effects. Nanoparticle-collector interaction forces dominate the transport in the potential region where deposition occurs. A sensitivity analysis of the model indicates that nanoparticle transport and deposition in a pore is significantly affected by various pore-scale parameters such as the nanoparticle and collector surface potentials, ionic strength of the solution, flow velocity, pore radius, and nanoparticle radius. The model is found to be more sensitive to all parameters under favorable conditions. It is found that the secondary minimum plays an important role in the deposition of small as well as large nanoparticles, and its contribution is found to increase as the favorability of the surface for adsorption decreases.
Correlation equations for average deposition rate coefficients of nanoparticles in a saturated cylindrical pore under unfavorable conditions are developed as a function of nine pore-scale parameters: the pore radius, nanoparticle radius, mean flow velocity, solution ionic strength, viscosity, temperature, solution dielectric constant, and nanoparticle and collector surface potentials. Advection-diffusion equations for nanoparticle transport are prescribed for the bulk and diffusion regions, while the interaction between the diffusion and potential regions is included as a boundary condition. This interaction is modeled as a first-order reversible kinetic adsorption. The expressions for the mass transfer rate coefficients between the diffusion and the potential regions are derived in terms of the interaction energy profile between the nanoparticle and the collector. The resulting equations are solved numerically for a range of values of pore-scale parameters. The nanoparticle concentration profile obtained for the cylindrical pore is averaged over a moving averaging volume within the pore in order to get the 1-D concentration field. The latter is fitted to the 1-D advection-dispersion equation with an equilibrium or kinetic adsorption model to determine the values of the average deposition rate coefficients. Pore-scale simulations are performed for three values of Péclet number, Pe = 0.05, 5 and 50. It is found that under unfavorable conditions, the nanoparticle deposition at pore scale is best described by an equilibrium model at low Péclet numbers (Pe = 0.05), and by a kinetic model at high Péclet numbers (Pe = 50). But, at an intermediate Pe (e.g., near Pe = 5), both equilibrium and kinetic models fit the 1-D concentration field. Correlation equations for the pore-averaged nanoparticle deposition rate coefficients under unfavorable conditions are derived by performing a multiple-linear regression analysis between the estimated deposition rate coefficients for a single pore and various pore-scale parameters. The correlation equations, which follow a power law relationship with nine pore-scale parameters, are found to be consistent with the column-scale and pore-scale experimental results, and qualitatively agree with CFT.
Nanoparticle transport is upscaled from pore to the Darcy scale in saturated porous media by incorporating the correlations equations for the pore-averaged deposition rate coefficients of nanoparticles in a cylindrical pore into a multi-directional pore-network model, PoreFlow (Raoof et al., 2013). Pore-network model simulations are performed for a range of parameter values, and nanoparticle BTCs are obtained from the pore-network model. Those curves are then modeled using 1-D advection-dispersion equation with a two-site first-order reversible deposition, with terms accounting for both equilibrium and kinetic sorption. Kinetic sorption is found to become important as the favorability of the surface for deposition decreases. Correlation equations for the Darcy¬scale deposition rate coefficients under unfavorable conditions are developed as a function of various measurable Darcy-scale parameters, including: porosity, mean pore throat radius, mean pore water velocity, nanoparticle radius, ionic strength, dielectric constant, viscosity, temperature, and surface potentials on the nanoparticle and grain surface. The correlation equations are found to be consistent with the observed trends from the column experiments available in the literature, and are in agreement with CFT for all parameters, except for the mean pore water velocity and nanoparticle radius. The Darcy-scale correlation equations contain multipliers whose values for a given set of experimental conditions need to be determined by comparing the values of the deposition rate coefficients predicted by the correlation equations against the estimated values of Darcy-scale deposition parameters obtained by fitting the BTCs from column or field experiments with 1-D advection-dispersion-deposition model. They account for the effect of factors not considered in this study, such as the physical and chemical heterogeneity of the grain surface and nanoparticles, flow stagnation points, grain-to-grain contacts, etc.
Colloids are abundant in the subsurface and have been observed to interact with a variety of contaminants, including viruses, thereby significantly influencing their transport. A mathematical model is developed to simulate the co-transport of viruses and colloids in partially saturated porous media under steady state flow conditions. The virus attachment to the mobile and immobile colloids is described using a linear reversible kinetic model. It is assumed that colloid transport is not affected by the presence of attached viruses on its surface, and hence, colloid transport is decoupled from virus transport. The governing equations are solved numerically using an alternating three-step operator splitting approach. The model is verified by fitting three sets of experimental data published in the literature: (1) Syngouna and Chrysikopoulos (2013) and (2) Walshe et al. (2010), both on the co-transport of viruses and clay colloids under saturated conditions, and (3) Syngouna and Chrysikopoulos (2015) for the co-transport of viruses and clay colloids under unsaturated conditions. The model results are found to be in good agreement with the observed BTCs under both saturated and unsaturated conditions.
Then, the developed model was used to simulate the co-transport of viruses and colloids in porous media under unsaturated conditions, with the aim of understanding the relative importance of various processes on the co-transport of viruses and colloids. The virus retention in porous media in the presence of colloids is greater under unsaturated conditions as compared to the saturated conditions due to: (1) virus attachment to the air-water interface (AWI), and (2) co-deposition of colloids with attached viruses on its surface to the AWI. A sensitivity analysis of the model to various parameters showed that virus attachment to AWI is the most sensitive parameter affecting the BTCs of both free viruses and total mobile viruses, and has a significant effect on all parts of the BTC. The free and the total mobile virus BTCs are mainly influenced by parameters describing virus attachment to the AWI, virus interactions with mobile and immobile colloids, virus attachment to solid-water interface (SWI), and colloid interactions with SWI and AWI. The virus BTC is relatively insensitive to parameters describing the maximum adsorption capacity of the AWI for colloids, inlet colloid concentration, virus detachment rate coefficient from the SWI, maximum adsorption capacity of the AWI for viruses, and inlet virus concentration
State and Parameter Estimation in Water Distribution Systems
No full textWater Distribution Systems (WDS) are one of the most important infrastructures in a society, which has a direct impact on human lives. Water authorities around the world are expected to supply the required amount of water at adequate pressure and potable quality to the consumers at all times. For achieving this task, proper system and resource management need to be carried out. Optimal or near optimal management of the WDS is rendered challenging due to various issues the system is plagued with. WDS around the world have to deal with many issues such as leakage, water quality degradation, demand uncertainty, inequitable supply of water among the consumers etc. For better system management, monitoring, modelling and control of the WDS is required. Currently, many sensors are available commercially which are capable of measuring the water quality and quantity parameters of the system at very fine frequency. These sensors make real time monitoring of the water distribution system easier. In case of a water quality or quantity event in the WDS, this real time feed of data might be able to give a preliminary warning to the water authorities. For development of a decision support system, this monitored data need to be coupled with modelling and analysis of the WDS. Efficiency of the decision support system will be enhanced if modelling and analysis is carried out in real time or near real time, as the time required for event mitigation will be reduced. Automatic control of appurtenances in the WDS helps in appropriate operation and execution in the system as well.
Data assimilation methods were used in many fields of science and engineering for parameter and state estimation of the system in real time. Different mathematical filters like Kalman Filter, Extended Kalman Filter, EnKF etc. and their variants are used in real time state estimation in the field of ocean sciences, remote sensing, ground water etc. Whereas, application of data assimilation techniques in WDS analysis is quite recent and limited. The non-linear nature of WDS, combined with dynamic operational conditions makes direct application of data assimilation techniques difficult in WDS
The primary objective of the present study is to apply Ensemble Kalman Filter (EnKF) based data assimilation techniques for water quality and quantity state estimation in WDS under different parameter uncertainties. EnKF and its variants are used in this study, as EnKF is one of the most widely used mathematical filter, and it is easy to apply it to WDS, as it does not make any assumptions of system linearity. In this study, EnKF and its variants were formulated
for hydraulic and water quality state estimation in WDS. EnKF and three different variants (Iterative or Dual- EnKF (I-EnKF), Non-Iterative Restart EnKF (NIR-EnKF) and Iterative Restart EnKF (IR-EnKF)) of it are used throughout the study.
The main objectives of this thesis are as follows:
i. Application of EnKF based data assimilation technique for hydraulic state estimation in WDS (Chapter 3).
ii. Application of EnKF based data assimilation for water quality state estimation in presence of source uncertainty (Chapter 4)
iii. Application of EnKF based data assimilation for water quality state estimation in presence of parameter uncertainty (Chapter 5)
iv. Application of EnKF based data assimilation and Modified- Fault Sensitivity Matrix (FSM) for leak detection and localization in WDS(Chapter 6).
v. Achieving equitable water supply for New Bangalore Inflow Network , using Dynamic Inversion-Proportional Integral Derivative (DI-PID) and Model Predictive Controller (MPC) technology (Chapter 7
Two Dimensional Numerical Modelling Of Variably Saturated Flows
No full textThe prediction of moisture and contaminant transport through unsaturated soil to ground water is becoming increasingly important in the fields of hydrology, agriculture and environmental engineering. Computer aided simulation techniques enables one to conduct a series of systematic numerical experiments to analyze flow phenomenon in subsurface hydrology under various physical and chemical processes. The flow movement depends upon medium characteristics, initial and boundary conditions, which reflect, physical processes occurring below the ground. To understand the effects of physical process an efficient and accurate model is needed. Thus the model developed must be able to handle varied initial and boundary conditions. In this regard, infiltration into a very dry soil becomes a very important problem of study.
Most of the earlier numerical models developed are concentrated on the development of an efficient algorithm or the modelling of a particular process which govern the flow in unsaturated or saturated-unsaturated homogeneous medium. Not much work has been done on the analysis of variably saturated flow in layered soil medium. Models to simulate unsaturated flow through dry soils, especially through layered soils with varied boundary conditions are very limited. Further, not much studies have been reported in the literature on the prediction of seepage face development and the phreatic surface movement in variably saturated media with layering. These aspects are very important in determining the flow field and the discharge from the domain. A detailed literature review covering above aspects has been made and is reported in this thesis.
In the present work, two dimensional numerical models to predict the movement of wetting front in unsaturated domain and the movement of the phreatic surface in homogeneous and layered porous media under various initial and boundary conditions are developed based on finite difference and finite volume techniques. These models can handle flow in both rectangular flow domains and radial flow domains. The initial condition settings include the handling of very dry soil medium without any transformation of the governing equation, handling of infiltration and constant head conditions at the boundaries under steady state as well as transient scenarios. The models are also able to handle various soil moisture characteristics which depict the nonlinear behaviour between hydraulic conductivity, moisture content and pressure head in a soil media.
A mixed form of the governing partial differential equation is used in the present study as it leads to better mass conservation. The finite difference model uses a central difference approximation for the space derivatives and an Eulerian backward difference approximation for the time derivative. The fully implicit formulation is solved iteratively using Strongly Implicit Procedure after making Picard approximation for the highly nonlinear coefficients. The process of infiltration into an initially dry soil leads to the development of a steep wetting front. As the finite volume technique is naturally an upwind method, it is expected to play a positive role in modelling such processes accurately. Hence, a finite volume model is also developed by approximating the convective part using a MUSCL approach and a fully implicit central difference method for the diffusive part of the governing equation.
The models developed are validated using both experimental data and numerical solutions for problems reported in the literature. The validation problems cover a wide range of physical scenarios such as: infiltration into a very dry soil, infiltration into a dry soil column with gravity drainage, development of water table mound, steady state drainage in a sand filled wedge shaped tank with seepage face development and transient seepage face development in a rectangular domain. Five test problems are used for the validation of the models. The developed models perform very well for the test problems considered, indicating the models' capability in handling such situations. The results obtained by using the present models for simulating flow through highly unsaturated (very dry) soils show that the models perform very well when compared with models which use transformation techniques to handle such problems. The performance of the present models in comparison with the experimental data and numerical models available in the literature show the suitability of the present models in handling such situations.
The present models are also used to analyse various types of unsaturated flow problems with varying initial and boundary conditions. The boundary conditions considered are no flow and /or free flow conditions along the left, right and bottom boundaries with infiltration condition along a part of the top boundary. For the various cases considered in the present study, infiltration rate varies from 5 cm/day to 50cm/day through an initially very dry soil of -15000 cm pressure head in homogeneous and layered soils. Different types of soil media considered vary from sandy loam, loam and clay with horizontal and vertical layering of these soils. A total number of 14 cases are analysed. The results are discussed in terms of pressure head variation in the flow domain along with moisture redistribution for all the cases under consideration. It is observed from these studies that the infiltration rate play an important role on the wetting front movement through layered soils depending on the type of layering and the boundary conditions considered. The soil properties of various layers affect the movement of wetting front by changing the direction of movement. Even though the wetting front movement is predominantly vertical, there is a tendency for the wetting front to move in the horizontal direction as it moves from a coarse soil to fine soil. It is also observed that the vertical layering of soils with different hydraulic conductivity helps in redirecting the flow towards the bottom boundary through the neighboring coarser layers.
As finite volume method is more suitable for simulating sharp fronts, it is expected to perform better than finite difference method for simulating infiltration into very dry soils. So, a comparison is made between the performance of these two models by using the above test problems. It is observed from these studies that the performance of both the models are same except that the finite volume method takes much more CPU time than the finite difference model. Considering the type of problems tested, it is observed that for modelling unsaturated and saturated-unsaturated flows, finite difference method is better in comparison to finite volume method. It may be due to the predominant diffusive characteristics of the governing equation even while modelling flow through very dry soils.
Proper estimation of the seepage height is an important aspect in finding the discharge through the porous medium. It is observed from the literature that the use of a saturated flow model in such situations can lead to an underestimation of the discharge through the porous medium. This effect is more important when dealing with small dimension problems. It is also observed that various parameters which govern the moisture movement through saturated-unsaturated regions affect proper estimation of the seepage face height and there by discharge. Various factors like effect of boundary conditions, type of soil layering, problem dimension and aspect ratio on seepage face development and the associated phreatic surface formation is studied in the present work. It is seen from the present study that the seepage face development is more in rectangular flow domain than in radial flow domain for both homogeneous and layered soils. It is also seen that the seepage face development in rectangular problems are more sensitive than radial flow problems for various factors considered. The seepage height is also influenced by the tail water level. It is seen from the present study that as the tail water level increases the seepage face reduces with no seepage face development for some of the cases studied. This influence is relatively less for radial flow problems. As the length of the domain increases the seepage height decreases. It is seen that for different cases with same horizontal dimension, as the height of the domain increases the seepage face height also increases. This phenomenon is observed for both homogeneous and layered soil medium. The influence of the aspect ratio, which is the ratio of the length to height of the domain indicate that as the aspect ratio increases the seepage height decreases. The type of the soil layering is observed to have a very strong influence on the seepage face development. The study for understanding the effect of soil layering on the development of seepage face and phreatic surface suggest that as the coarseness of the material increases, the phreatic surface become flatter and its steepness increases with the fineness of the soil.
The present model is also used for studying the transient phreatic surface movement and the seepage face development. This is studied for homogeneous and layered rectangular soil medium. The present study is used to understand the effect of specific storage on the phreatic surface movement and the seepage face development. The studies indicate that the influence of specific storage on the seepage face development is insignificant in homogeneous soils with only very little effect in the early time for longer domains. It is also observed that the influence of the specific storage is significant in the case of layered soils. This effect depends on the type of layering and the problem dimension and is observed to have influence for relatively longer period. This observation suggests the importance of specific storage on transient seepage face development. When the specific storage effect is considered the drainage of the soil become faster resulting in a faster decline of the phreatic surface with time. The influence of specific storage is also studied considering the problem dimension effect. It is seen that as the aspect ratio increases, the effect of specific storage on the phreatic surface development decreases. The studies with change in the upstream boundary condition from a constant head to a no flow condition indicate that the effect of specific storage has no significant influence on the phreatic surface development for both homogeneous and layered soils
Analytical and Numerical Modeling for Heat Transport in a Geothermal Reservoir due to Cold Water Injection
Get PDFGeothermal energy is the energy naturally present inside the earth crust. When a large
volume of hot water and steam is trapped in subsurface porous and permeable rock
structure and a convective circulating current is set up, it forms a geothermal reservoir. A geothermal system can be defined as - convective water in the upper crust of earth, which
transfers heat from a heat source (in the reservoir) to a heat sink, usually the free surface. A geothermal system is made up of three main elements: a heat source, a reservoir and a fluid, which is the carrier that transfers the heat.
As an alternative source of energy geothermal energy has been under attention of the researchers for quite some time. The reason behind this is the existence of several benefits like clean and renewable source of energy which has considerable environmental advantage, with no chemical pollutants or wastes are generated due to geothermal emissions, and the reliability of the power resource. Hence research has been directed in several directions like exploration of geothermal resources, modeling the characteristics of different types of geothermal reservoirs and technologies to extract energy from them. The target of these models has been the prediction of the production of the hot water and steam and thus the estimation of the electricity generating potential of a geothermal reservoir in future years.
In a geothermal power plant reinjection of the heat depleted water extracted from the
geothermal reservoir has been a common practice for quite some time. This started for safe wastewater disposal and later on the technology was employed to obtain higher efficiency of heat and energy extraction. In most of the cases a very small fraction of the thermal energy present in the reservoir can be recovered without the reinjection of geothermal fluid. Also maintaining the reservoir pressure is essential which gradually reduces due to continuous extraction of reservoir fluid without reinjection, especially for reservoirs with
low permeabilities. Although reinjection of cold-water has several benefits, the possibility of premature breakthrough of the cold-water front, from injection well zone to production well zone, reduces the efficiency of the reservoir operation drastically. Hence for maintaining the reservoir efficiency and longer life of the reservoir, the injectionproduction well scheme is to be properly designed and injection and extraction rates are to
be properly fixed.
Modeling of flow and heat transport in a geothermal reservoir due to reinjection of coldwater has been attempted by several researchers analytically, numerically and
experimentally. The analytical models which exist in this field deal mostly with a single injection well model injecting cold-water into a confined homogeneous porous-fractured geothermal reservoir. Often the thermal conductivity is neglected in the analytical study considering it to be negligible which is not always so, as proved in this study. Moreover heterogeneity in the reservoir is also a major factor which has not been considered in any such analytical study. In the field of numerical modeling there also exists a need of a general coupled three-dimensional thermo-hydrogeological model including all the modes of heat transport (advection and conduction), the heat loss to the confining rocks, the regional groundwater flow and the geothermal gradient. No study existing so far reported such a numerical model including those mentioned above.
The present study is concerned about modeling the non-isothermal flow and heat transport in a geothermal reservoir due to reinjection of heat depleted water into a geothermal reservoir. Analytical and numerical models are developed here for the transient temperature distributions and advancement of the thermal front in a geothermal reservoir which is generated due to the cold-water injection. First homogeneous geothermal aquifers are considered and later heterogeneities of different kinds are brought into picture. Threedimensional numerical models are developed using a software code DuMux which solves flow and heat transport problems in porous media and can handle both single and
multiphase flows. The results derived by the numerical models have been validated using
the results from the analytical models derived in this study.
Chapter 1 of the thesis gives a brief introduction about different types of geothermal reservoirs, followed by discussion on the governing differential equations, the conceptual model of a geothermal reservoir system, the efficiency of geothermal reservoirs, the modeling and simulation concepts (models construction, boundary conditions, model calibration etc.). Some problems related with geothermal reservoirs and geothermal power is also discussed. The scenario of India in the context having a huge geothermal power
potential is described and different potential geothermal sites have been pointed out.
In Chapter 2, the concept of reinjection of the heat depleted (cold) water into the
geothermal reservoir is introduced. Starting with a brief history of the geothermal
reinjection, the chapter describes the purpose and the need of reinjection of geothermal fluid giving examples of different geothermal fields over the world where reinjection has been in practice and benefitted by that. The chapter further discusses on the problems and obstacles faced by the geothermal projects resulting from the geothermal reinjection, most
important of which is the thermal-breakthrough and cooling of production wells. Lastly the problem of this thesis is discussed which is to model the transient temperature distribution and the movement of the cold-water thermal front generated due to the reinjection. The need of this modeling is elaborated which represents the motivation of taking up the problem of the thesis.
Chapter 3 describes an analytical model developed for the transient temperature in a
porous geothermal reservoir due to injection of cold-water. The reservoir is composed of a confined aquifer, sandwiched between rocks of different thermo-geological properties. The heat transport processes considered are advection, longitudinal conduction in the
geothermal aquifer, and the conductive heat transfer to the underlying and overlying rocks of different geological properties. The one-dimensional heat transfer equation has been solved using the Laplace transform with the assumption of constant density and thermal properties of both rock and fluid. Two simple solutions are derived afterwards, first neglecting the longitudinal conductive heat transport and then heat transport to confining rocks. The analytical solutions represent the transient temperature distribution in the geothermal aquifer and the confining rocks and model the movement of the cold-water thermal front in them. The results show that the heat transport to the confining rocks plays an influential role in the transient heat transport here. The influence of some parameters, e.g. the volumetric injection rate, the longitudinal thermal conductivity and the porosity of
the porous media, on the transient heat transport phenomenon is judged by observing the variation of the transient temperature distribution with different values of the parameters.
The effects of injection rate and thermal conductivity have been found to be high on the results.
Chapter 4 represents another analytical model for transient temperature distribution in a heterogeneous geothermal reservoir underlain and overlain by impermeable rocks due to injection of cold-water. The heterogeneity of the porous medium is expressed by the spatial variation of the flow velocity and the longitudinal effective thermal conductivity of the medium. Simpler solutions are also derived afterwards first neglecting the longitudinal conduction, then the heat loss to the confining rocks depending on the situation where the
contribution of them to the transient heat transport phenomenon in the porous media is
negligible. Solution for a homogeneous aquifer with constant values of the rock and fluid parameters is also derived with an aim to compare the results with that of the
heterogeneous one. The effect of heat loss to the confining rocks in this case is also
determined and the influence of some of the parameters involved, on the transient heat
transport phenomenon is assessed by observing the variation of the results with different magnitudes of those parameters. Results show that the heterogeneity plays a major role in controlling the cold-water thermal front movement. The transient temperature distribution in the geothermal reservoir depends on the type of heterogeneity. The heat loss to the confining rocks of the geothermal aquifer also has influence on the heat transport
phenomenon.
In Chapter 5 another analytical model is derived for a heterogeneous reservoir where the heterogeneous geothermal aquifer considered is a confined aquifer consisted of homogeneous layers of finite length and overlain and underlain by impermeable rock
media. All the different layers in the aquifer and the overlying and underlying rocks are of different thermo-hydrogeological properties. Results show that the advancement of the cold-water thermal front is highly influenced by the layered heterogeneity of the aquifer.
As the cold-water thermal front encounters layers of different thermo-hydrogeological
properties the movement of it changes accordingly. The analytical solution derived here has been compared with a numerical model developed by the multiphysics software code COMSOL which shows excellent agreement with each other. Lastly it is shown that
approximation of the properties of a geothermal aquifer by taking mean of the properties of all the layers present will lead to erroneous estimation of the temperature distribution.
Chapter 6 represents a coupled three-dimensional thermo-hydrogeological numerical
model for transient temperature distribution in a confined porous geothermal aquifer due
to cold-water injection. This 3D numerical model is developed for solving more practical problems which eliminate the assumptions taken into account in analytical models. The numerical modeling is performed using a software code DuMux as mentioned before.
Besides modeling the three-dimensional transient temperature distribution in the model domain, the chapter investigates the regional groundwater flow has been found to be a very important parameter to consider. The movement of the thermal front accelerates or decelerates depending on the direction of the flow. Influence of a few parameters involved in the study on the transient heat transport phenomenon in the geothermal reservoir domain, namely the injection rate, the permeability of the confining rocks and the thermal conductivity of the geothermal aquifer is also evaluated in this chapter. The models have been validated using analytical solutions derived in this thesis. The results are in very good agreement with each other.
In Chapter 7 the main conclusions drawn from the study have been enlisted and the scope
of further research is also pointed out
- …
