Investigation of spatial patterns of ground-water exchange with lakes, using a three-dimensional numerical model

Hydrogeologic variables controlling groundwater exchange with inflow and flow-through lakes were simulated using a three-dimensional numerical model (MODFLOW) to investigate and quantify spatial patterns of lake bed seepage and hydraulic head distributions in the porous medium surrounding the lakes. Also, the total annual inflow and outflow were calculated as a percentage of lake volume for flow-through lake simulations. The general exponential decline of seepage rates with distance offshore was best demonstrated at lower anisotropy ratio (i.e., Kh/Kv = 1, 10), with increasing deviation from the exponential pattern as anisotropy was increased to 100 and 1000. 2-D vertical section models constructed for comparison with 3-D models showed that groundwater heads and seepages were higher in 3-D simulations. Addition of low conductivity lake sediments decreased seepage rates nearshore and increased seepage rates offshore in inflow lakes, and increased the area of groundwater inseepage on the beds of flow-through lakes. Introduction of heterogeneity into the medium decreased the water table and seepage ratesnearshore, and increased seepage rates offshore in inflow lakes. A laterally restricted aquifer located at the downgradient side of the flow-through lake increased the area of outseepage. Recharge rate, lake depth and lake bed slope had relatively little effect on the spatial patterns of seepage rates and groundwater exchange with lakes.

Estudo paramétrico da injeção de polímeros em reservatórios de petróleo

Developing an efficient methodology for oil recovery is extremely important . Within the range of enh anced oil recovery, known as EOR, the injection of polymer solutions becomes effective in controlling the mobility of displacing fluid . This method consists of adding polymers to the injection water to increase its viscosity, so that more water diffuses in to the porous medium and increasing the sweep efficiency in the reservoir. This work is studied by numerical simulation , application of the injection polymer solution in a homogeneous reservoir , semisynthetic with similar characteristics to the reservoirs of the Brazilian Northeast , numerical simulations were performed using thermal simulator STARS from CMG (Computer Modelling Group ). The study aimed to analyze the influence of some parameters on the behavior of reservoir oil production, with the response to cumulative production. Simulations were performed to analyze the influence of water injection, polymer solution and alternating injection of water banks and polymer solution, comparing the results for each simulated condition. The primary outcomes were: oil viscosity, percentage of injected polymer, polymer viscosity and flow rate of water injection. The evaluation of the influence of variables consisted of a complete experimental design followed a Pareto analysis for the purpose of pointing out which va riables would be most influential on the response represented b y the cumulative oil production . It was found that all variables significantly influenced the recovery of oil and the injection of polymer solution on an ongoing basis is more efficient for the cumulative production compared to oil recovery by continuous water injection. The primary recovery show ed low levels of oil production , water injection significantly improves the pro duction of oil in the reservoir , but the injection of polymer solution em erges as a new methodology to increase the production of oil, increasing the life of the well and possible reduction of water produced.

Estudo das forças capilar, viscosa e gravitacional através de grupos adimensionais na simulação da drenagem gravitacional assistida por gás (GAGD)

The distribution and mobilization of fluid in a porous medium depend on the capillary, gravity, and viscous forces. In oil field, the processes of enhanced oil recovery involve change and importance of these forces to increase the oil recovery factor. In the case of gas assisted gravity drainage (GAGD) process is important to understand the physical mechanisms to mobilize oil through the interaction of these forces. For this reason, several authors have developed physical models in laboratory and core floods of GAGD to study the performance of these forces through dimensionless groups. These models showed conclusive results. However, numerical simulation models have not been used for this type of study. Therefore, the objective of this work is to study the performance of capillary, viscous and gravity forces on GAGD process and its influence on the oil recovery factor through a 2D numerical simulation model. To analyze the interplay of these forces, dimensionless groups reported in the literature have been used such as Capillary Number (Nc), Bond number (Nb) and Gravity Number (Ng). This was done to determine the effectiveness of each force related to the other one. A comparison of the results obtained from the numerical simulation was also carried out with the results reported in the literature. The results showed that before breakthrough time, the lower is the injection flow rate, oil recovery is increased by capillary force, and after breakthrough time, the higher is the injection flow rate, oil recovery is increased by gravity force. A good relationship was found between the results obtained in this research with those published in the literature. The simulation results indicated that before the gas breakthrough, higher oil recoveries were obtained at lower Nc and Nb and, after the gas breakthrough, higher oil recoveries were obtained at lower Ng. The numerical models are consistent with the reported results in the literature

Análise da eficácia de um refrigerador doméstico com condensador modificado: revestimento com meio poroso dissipativo

The development of home refrigerators generally are compact and economic reasons for using simplified configuration. The thermodynamic coefficient of performance ( COP ) is limited mainly in the condenser design for reasons of size and arrangement ( layout ) of the project ( design ) and climatic characteristics of the region where it will operate. It is noteworthy that this latter limitation is very significant when it comes to a country of continental size like Brazil with diverse climatic conditions. The COP of the cycle depends crucially on the ability of heat dissipated in the condenser. So in hot climates like the northeast, north, and west-central dispel ability is highly attenuated compared to the south and southeast regions with tropical or subtropical climates when compared with other regions. The dissipation in compact capacitors for applications in domestic refrigeration has been the focus of several studies, that due to its impact on reducing costs and power consumption, and better use of the space occupied by the components of refrigeration systems. This space should be kept to a minimum to allow an increase in the useful storage volume of refrigerator without changing the external dimensions of the product. Due to its low cost manufacturing, wire on tube condensers continue to be the most advantageous option for domestic refrigeration. Traditionally, these heat exchangers are designed to operate under natural convection. Not always, the benefits of greater compactness of capacitors for forced outweigh the burden of pumping air through the external heat exchanger. In this work we propose an improvement in convective condenser changing it to a transfer mechanism combined in series with conductive pipes and wire to a moist convective porous medium and the porous medium to the environment. The porous media used in the coating was composed of a gypsum plaster impregnated fiber about a mesh of natural cellulosic molded tubular wire mesh about the original structure of the condenser , and then dried and calcined to greater adherence and increased porosity. The proposed configuration was installed in domestic refrigeration system ( trough ) and tested under the same conditions of the original configuration . Was also evaluated in the dry condition and humidified drip water under natural and forced with an electro - fan ( fan coil ) convection. Assays were performed for the same 134- refrigerant charge e under the same thermal cooling load. The performance was evaluated in various configurations, showing an improvement of about 72 % compared with the original configuration proposed in humidification and natural convection.

Fluxo de fluído através de um meio poroso fractal desordenado. Análise das tensões de cisalhamento e efeito de escala na estimativa das forças viscosas

In this work we have investigated some aspects of the two-dimensional flow of a viscous Newtonian fluid through a disordered porous medium modeled by a random fractal system similar to the Sierpinski carpet. This fractal is formed by obstacles of various sizes, whose distribution function follows a power law. They are randomly disposed in a rectangular channel. The velocity field and other details of fluid dynamics are obtained by solving numerically of the Navier-Stokes and continuity equations at the pore level, where occurs actually the flow of fluids in porous media. The results of numerical simulations allowed us to analyze the distribution of shear stresses developed in the solid-fluid interfaces, and find algebraic relations between the viscous forces or of friction with the geometric parameters of the model, including its fractal dimension. Based on the numerical results, we proposed scaling relations involving the relevant parameters of the phenomenon, allowing quantifying the fractions of these forces with respect to size classes of obstacles. Finally, it was also possible to make inferences about the fluctuations in the form of the distribution of viscous stresses developed on the surface of obstacles.

Modelagem multiescala de fenômenos eletrocinéticos em meios porosos carregados eletricamente: aplicação a meios porosos argilosos

In this master thesis, we propose a multiscale mathematical and computational model for electrokinetic phenomena in porous media electrically charged. We consider a porous medium rigid and incompressible saturated by an electrolyte solution containing four monovalent ionic solutes completely diluted in the aqueous solvent. Initially we developed the modeling electrical double layer how objective to compute the electrical potential, surface density of electrical charges and considering two chemical reactions, we propose a 2-pK model for calculating the chemical adsorption occurring in the domain of electrical double layer. Having the nanoscopic model, we deduce a model in the microscale, where the electrochemical adsorption of ions, protonation/ deprotonation reactions and zeta potential obtained in the nanoscale, are incorporated through the conditions of interface uid/solid of the Stokes problem and transportation of ions, modeled by equations of Nernst-Planck. Using the homogenization technique of periodic structures, we develop a model in macroscopic scale with respective cells problems for the e ective macroscopic parameters of equations. Finally, we propose several numerical simulations of the multiscale model for uid ow and transport of reactive ionic solute in a saturated aqueous solution of kaolinite. Using nanoscopic model we propose some numerical simulations of electrochemical adsorption phenomena in the electrical double layer. Making use of the nite element method discretize the macroscopic model and propose some numerical simulations in basic and acid system aiming to quantify the transport of ionic solutes in porous media electrically charged.

Modelagem multiescala de fenômenos eletrocinéticos em meios porosos carregados eletricamente: aplicação a meios porosos argilosos

In this master thesis, we propose a multiscale mathematical and computational model for electrokinetic phenomena in porous media electrically charged. We consider a porous medium rigid and incompressible saturated by an electrolyte solution containing four monovalent ionic solutes completely diluted in the aqueous solvent. Initially we developed the modeling electrical double layer how objective to compute the electrical potential, surface density of electrical charges and considering two chemical reactions, we propose a 2-pK model for calculating the chemical adsorption occurring in the domain of electrical double layer. Having the nanoscopic model, we deduce a model in the microscale, where the electrochemical adsorption of ions, protonation/ deprotonation reactions and zeta potential obtained in the nanoscale, are incorporated through the conditions of interface uid/solid of the Stokes problem and transportation of ions, modeled by equations of Nernst-Planck. Using the homogenization technique of periodic structures, we develop a model in macroscopic scale with respective cells problems for the e ective macroscopic parameters of equations. Finally, we propose several numerical simulations of the multiscale model for uid ow and transport of reactive ionic solute in a saturated aqueous solution of kaolinite. Using nanoscopic model we propose some numerical simulations of electrochemical adsorption phenomena in the electrical double layer. Making use of the nite element method discretize the macroscopic model and propose some numerical simulations in basic and acid system aiming to quantify the transport of ionic solutes in porous media electrically charged.

Evaluation of Gaseous Fluid Dispersion Through Packed Bed Using Radiotracers Technique

In this study, it was developed a methodology for the determination of the dispersion of a gaseous tracer in porous media using the radiotracer technique. In order to evaluate several porous media, a cylindrical filter was constructed in PVC and connected to a system with constant flow. Inside this unit silica crystals (16-20) mesh was used as porous media and CH3Br (Methyl Bromide) marked with 82Br was used as radiotracer. An instantaneous pulse of tracer was applied in the system entrance and registered by two NaI (3x3)” scintillation detectors located one before and the other after the filter. The curves produced by the radioactive cloud and recorded by the detector were analyzed statistically using the weighted moment method. The mathematical model one considered as great dispersion of tracer was used to evaluate the flow conditions inside the filter system. The results show us that the weight moment method associated with radiotracer techniques is useful to evaluated an industrial filter and allows to measure the residence time distribution, τ, and the axial dispersion, DAB, gas in a porous medium.

Simulação numérica da convecção mista em cavidade preenchida com meio poroso heterogêneo e homogêneo

In this work is presented mixed convection heat transfer inside a lid-driven cavity heated from below and filled with heterogeneous and homogeneous porous medium. In the heterogeneous approach, the solid domain is represented by heat conductive equally spaced blocks; the fluid phase surrounds the blocks being limited by the cavity walls. The homogeneous or pore-continuum approach is characterized by the cavity porosity and permeability. Generalized mass, momentum and energy conservation equations are obtained in dimensionless form to represent both the continuum and the pore-continuum models. The numerical solution is obtained via the finite volume method. QUICK interpolation scheme is set for numerical treatment of the advection terms and SIMPLE algorithm is applied for pressure-velocity coupling. Aiming the laminar regime, the flow parameters are kept in the range of 102≤Re≤103 and 103≤Ra≤106 for both the heterogeneous and homogeneous approaches. In the tested configurations for the continuous model, 9, 16, 36, and 64 blocks are considered for each combination of Re and Ra being the microscopic porosity set as constant φ=0,64 . For the pore-continuum model the Darcy number (Da) is set according to the number of blocks in the heterogeneous cavity and the φ. Numerical results of the comparative study between the microscopic and macroscopic approaches are presented. As a result, average Nusselt number equations for the continuum and the pore continuum models as a function of Ra and Re are obtained.

A new method for conditioning stochastic groundwater flow models in fractured media

Many geological formations consist of crystalline rocks that have very low matrix permeability but allow flow through an interconnected network of fractures. Understanding the flow of groundwater through such rocks is important in considering disposal of radioactive waste in underground repositories. A specific area of interest is the conditioning of fracture transmissivities on measured values of pressure in these formations. This is the process where the values of fracture transmissivities in a model are adjusted to obtain a good fit of the calculated pressures to measured pressure values. While there are existing methods to condition transmissivity fields on transmissivity, pressure and flow measurements for a continuous porous medium there is little literature on conditioning fracture networks. Conditioning fracture transmissivities on pressure or flow values is a complex problem because the measurements are not linearly related to the fracture transmissivities and they are also dependent on all the fracture transmissivities in the network. We present a new method for conditioning fracture transmissivities on measured pressure values based on the calculation of certain basis vectors; each basis vector represents the change to the log transmissivity of the fractures in the network that results in a unit increase in the pressure at one measurement point whilst keeping the pressure at the remaining measurement points constant. The fracture transmissivities are updated by adding a linear combination of basis vectors and coefficients, where the coefficients are obtained by minimizing an error function. A mathematical summary of the method is given. This algorithm is implemented in the existing finite element code ConnectFlow developed and marketed by Serco Technical Services, which models groundwater flow in a fracture network. Results of the conditioning are shown for a number of simple test problems as well as for a realistic large scale test case.

Modelagem e previsão da perda de injetividade em poços canhoneados

Waterflooding is a technique largely applied in the oil industry. The injected water displaces oil to the producer wells and avoid reservoir pressure decline. However, suspended particles in the injected water may cause plugging of pore throats causing formation damage (permeability reduction) and injectivity decline during waterflooding. When injectivity decline occurs it is necessary to increase the injection pressure in order to maintain water flow injection. Therefore, a reliable prediction of injectivity decline is essential in waterflooding projects. In this dissertation, a simulator based on the traditional porous medium filtration model (including deep bed filtration and external filter cake formation) was developed and applied to predict injectivity decline in perforated wells (this prediction was made from history data). Experimental modeling and injectivity decline in open-hole wells is also discussed. The injectivity of modeling showed good agreement with field data, which can be used to support plan stimulation injection wells

Simulação de fluxo de fluidos em meios porosos desordenados uma análise de efeito de escala na estimativa da permeabilidade e do coeficiente de arrasto

The present study provides a methodology that gives a predictive character the computer simulations based on detailed models of the geometry of a porous medium. We using the software FLUENT to investigate the flow of a viscous Newtonian fluid through a random fractal medium which simplifies a two-dimensional disordered porous medium representing a petroleum reservoir. This fractal model is formed by obstacles of various sizes, whose size distribution function follows a power law where exponent is defined as the fractal dimension of fractionation Dff of the model characterizing the process of fragmentation these obstacles. They are randomly disposed in a rectangular channel. The modeling process incorporates modern concepts, scaling laws, to analyze the influence of heterogeneity found in the fields of the porosity and of the permeability in such a way as to characterize the medium in terms of their fractal properties. This procedure allows numerically analyze the measurements of permeability k and the drag coefficient Cd proposed relationships, like power law, for these properties on various modeling schemes. The purpose of this research is to study the variability provided by these heterogeneities where the velocity field and other details of viscous fluid dynamics are obtained by solving numerically the continuity and Navier-Stokes equations at pore level and observe how the fractal dimension of fractionation of the model can affect their hydrodynamic properties. This study were considered two classes of models, models with constant porosity, MPC, and models with varying porosity, MPV. The results have allowed us to find numerical relationship between the permeability, drag coefficient and the fractal dimension of fractionation of the medium. Based on these numerical results we have proposed scaling relations and algebraic expressions involving the relevant parameters of the phenomenon. In this study analytical equations were determined for Dff depending on the geometrical parameters of the models. We also found a relation between the permeability and the drag coefficient which is inversely proportional to one another. As for the difference in behavior it is most striking in the classes of models MPV. That is, the fact that the porosity vary in these models is an additional factor that plays a significant role in flow analysis. Finally, the results proved satisfactory and consistent, which demonstrates the effectiveness of the referred methodology for all applications analyzed in this study.

Porous, catalytically active palladium nanostructures by tuning nanoparticle interactions in an organic medium

We present a simple template-free method for the synthesis of interconnected hierarchical porous palladium nanostructures by controlling the aggregation of nanoparticles in organic media. The interaction between the nanoparticles is tuned by varying the dielectric constant of the medium consistent with DLVO calculations. The reaction products range from discrete nanoparticles to compact porous clusters with large specific surface areas. The nanoclusters exhibit hierarchical porosity and are found to exhibit excellent activity towards the reduction of 4-nitrophenol into 4-aminophenol and hydrogen oxidation. The method opens up possibilities for synthesizing porous clusters of other functional inorganics in organic media.