Sharif Digital Repository

In This work, the performance of a finite volume (FV) numerical method is examined for simulating naturally fractured hydrocarbon reservoir. The dual-porosity approach applied for fracture modeling, and the flow transfer formulated based on a generalized black oil model which is based on matrix formulation, leading to a system of equations prone to degeneracy. This FV high-resolution central scheme numerical method is not sensitive to the choice of CFL number like some of other central schemes and can conveniently handle degenerate equations appearing in the reservoir simulation governing equation. In this effort, one dimension multi phase flow simulated in naturally fractured reservoir and... 

One of the most important issues in the petroleum industry, is how to increase the oil production from hydrocarbon reservoirs. The optimum design of an enhanced oil recovery process needs the knowledge on the physical phenomena of the working flow in oil reservoirs. An important aspect of any oil recovery process is the effectiveness of process fluid in removing oil from the rock pores at microscopic scales, in a way that it determines the success or failure of the recovery process. Meanwhile, one of the most important factors on the fluid mobility and as a result the microscopic mobility efficiency is the pressure difference between different phases in narrow paths of porous media. This... 

In this thesis the oil reservoir simulation is carried out by using a combination of fast Multi-scale finite volume method and Streamline method considering gravity effects. In this simulation, the governing equations for incompressible fluid flow in porous media are solved using the IMPES (implicit pressure- explicit saturation) viewpoint. Multi- scale method is rapid (equivalent on coarse-scale methods) and accurate (equivalent on fine-scale methods) in solving the pressure field. Using the pressure field obtained through this way combined with Darcy equation make us able to extract the conservative velocity field. In this study, analytical and numerical Streamline methods are used to... 

Fractured reservoirs show a different behavior from common reservoirs because of the existence of a broad network of fractures. This phenomenon makes it necessary to apply special methods for fractured reservoirs in the procedure of reservoir simulation. Since twenty percent of petroleum content in the world is buried in fractured reservoirs, investigating these reservoirs is of great importance.
The first step in simulation of these kinds of reservoirs is to come up with a geometrical model which can be used to take the fracture network influence into account. In the course for reaching such an objective, various models have been developed which are based on specific assumptions and in... 

Magnetohydrodynamic (MHD) flow is concerned with the interaction between the magnetic field and electrically conducting fluid flow. In this investigation, the MHD flow in the duct of magnetohydrodynamic wave energy conversion (MWEC) generator was simulated by use of OpenFOAM software. In the MWEC generator the reciprocating motion of ocean waves lead to electrically conducting fluid flow in the duct of this generator. By applying magnetic field perpendicular to this flow, electric current generate in the fluid. This generator has a higher efficiency compared with other ocean wave energy conversion generators. Firstly, to evaluate accuracy and performance of mhdFOAM code of OpenFOAM, three... 

Innovations in laser, microwaves, and similar technologies have significantly advanced thermal treatments for diseases or even injuries concerning skin tissue. For a thorough understanding in the underlying mechanisms of bioheat transfer behavior of skin,a1D unsteady non-dimensional hyperbolic model of heat transfer through this tissue with metabolic heat generation which is subject to specific boundary conditions, is solved numerically using the finite difference method. A thermal shock is generated at the base of the tissue, which moves forward with a finite speed. A Numerical solution for a simple one-layer skin tissue is obtained. Then, the effects of various parameters, time step,... 

In this work, simulation of 3D three-phase naturally fractured reservoirs is studied. Both Dual Porosity and Dual Porosity-Dual Permeability models which are commonly used in commercial softwares such as ECLIPSE and CMG, are utilized. The governing equations include the mass (in Black-oil model) and linear momentum (in Darcy law form) equations, which are converted into pressure equations and a set of transport equations for the saturations and mole fractions of phases. These equations are solved once for the matrix domain and subsequently for the fracture domain. These two set of equations are coupled through transfer functions. The accuracy of computations significantly depends on how... 

In this study, bulk viscosity of suspensions is predicted using numerical simulation.The elaborated numerical method consists of the lattice Boltzmann (LB) method for fluid matrix simulation and the immersed boundary (IB) method for solid boundary implementation. In order to successfully implement the IB method in the frame-work of the LB method, two important challenges, which are body force modeling in the LB method and numerical velocity slip on solid boundaries, are discussed. As a first challenge, body force modeling in the LB method is studied by using a uni-fied framework which has been devised in this study. Using this unified framework,firstly, errors of different body force... 

In this study, simulation of two-phase water-oil compressible flow in naturally fractured reservoirs is studied. The studied reservoir has a 1D and horizontal ge-ometry that its properties such as porosity and permeability are homogeneous.Dual Porosity and Dual Porosity-Dual Permeability models are employed for simulation of reservoir. Reservoir descriptive equations are conservation of mass and momentum. Conservation of momentum equation is simplified by Darcy law and two-phase fluid flow is simulated with Black-oil model. Discretization of equations is done by finite volume method. In order to solve the equations of two-phase compressible fluid flow, the so-called IMPES (Implicit Pressure... 

In this project, a response of a single well in a gas condensate reservoir has been investigated by compositional modeling and implicit numerical method. The most significant goals of this paper are as follows: 1- Examination of the flow produced in a gas condensate reservoir in constant pressure or mass rate. 2- Investigation of variables affecting the productivity of reservoir. 3- The study of effect of injection and gas cycling to the well on the production. A three-zone model has been used in this study as if there is only gas phase in the third and second zones. Condensation occurs since the pressure drops below the dewpoint pressure in the second zone; however, gas phase is only... 

In this project، the immiscible two phase flow of oil-water in the absence of gravity in the two dimensional domain is investigated. The uncertainty quantification of pressure and saturation of fluid in the reservoir are calculated by the methods of statistical moment equation (presented by Tchelepi) and probabilistic collocation method combined with Karhunen Loeve expansion (presented by Heng Lee). Then the results are compared with Monte Carlo simulation. Matlab Reservoir Simulation Toolbox is used for flow simulation. The main problem based on Tchelepi’s work[56[ ، is a horizontal two dimensional problem that there is two injection and production wells at the two end point grids (a... 

Reservoir simulation is a useful tool for efficient reservoir management and protective exploitation of non-renewable hydrocarbon resources as the main supplier of development programs in the country. Due to the nature of underground hydrocarbon reservoirs, the simulation faces several challenges, both in rock and fluid parts. These challenges should be overcome in order to provide reliable results from simulations. One of the major challenges pertaining to rock, is large discontinuous variations accompanied by anisotropy of rock permeability fields which causes numerous problems when governing equations are solved numerically. A viable numerical procedure shall account for complex reservoir... 

The main goal of this study is simulation of the flow-like landslide phenomenon and investigation of the effect of bed material entrainment on its behavior. Debris flows are considered as one of the types of flow-like landslides, and flow of dry granular materials on horizontal and inclined beds as a simple model of these types of flows in nature is numerically simulated. Also, in order to investigate the effect of entrainment on dynamics of flow, a bed containing granular materials is installed on horizontal and inclined beds. The governing equations of the flow are solved in the continuum context using the Smoothed Particle Hydrodynamics (SPH) method and the depth-averaged assumptions are... 

It has been demonstrated that the key features of fracture networks can be identified with the help of high-resolution aerial images of outcrops (rock formations exposed at the Earth's surface) obtained by drones and photogrammetry. The aim of this study is to utilize these images to improve the geometrical modeling of Naturally Fractured Reservoirs (NFRs). Initially, an algorithm is developed for processing aerial images of outcrops, so that the pattern of fracture network can be obtained. After obtaining the fracture network, the Embedded Discrete Fracture Model (EDFM) implemented in the MATLAB Reservoir Simulation Toolbox (MRST) is employed to study two-phase (oil-water) flows in the... 

Measuring flow rate in various industries is an operational challenge that requires to be accurate and fast to avoid wasting energy and money. Various flowmeters are employed in different industries using different principles and methods. Sonar flowmeters are one of the acoustic methods of fluid flow measurement, which do this by listening to the sound of the flow via arrays of sensors. Sonar flow measurement can be achieved by finding either the convection velocity of eddies or the speed of sound in the fluid. The first method is suitable for measuring single-phase flows, which is the subject of this thesis. This method needs a detailed study of turbulence. In the present work, this is... 

In this research, a novel scheme for flux approximation in finite volume central schemes have been presented, evaluated and used in PARSI reservoir simulator to solve black oil equations with Trangenstein & Bell formulation with gravitational effects. Also, a numerical method, have been proposed to increase the ability of modeling of capillary phenomena to PARSI.Edward's dominant wave-capturing scheme (2005), with Harten-Hymen entropy fix (1983), has been used to solve flow equations with gravity. This scheme relies on the detection of the dominant wave in the system without recourse to characteristic decomposition and upwinding while avoiding the excessive numerical diffusion that is... 

The goal of this research is to introduce and improve the new multiscale finite volume method for simulation of multiphase flow in porous media. This method has recently presented as an efficient approach in computational fluids dynamics field. This thesis contains the followings: In the first step a brief introduction is presented about petroleum reservoir engineering and different methods of oil production. Then the importance of new simulation methods to anticipate the behavior of these reservoirs with detailed geographical information is discussed. In the second section black-oil formulation in petroleum reservoirs is descritized with the aid of finite volume method. Using general... 

In this work, the combination of a multi-scale finite volume and streamline methods is presented for subsurface flow modeling. The used multi-scale method in this work is the same with the traditional type of it in algorithm and formulation. However, because of using a new mesh structure for implementation of multi-scale finite volume method, the new method named by Staggered Mesh Multi-scale Finite Volume (SMMsFV) method. Using the staggered mesh has some advantages such as reducing the computational cost and increasing the accuracy of the multi-scale method. In the SMMsFV method, first the coarse grid and dual coarse grid are constructed on the underlying fine grid. Then, the basis and... 

Study of two-phase flowin porous media forms a basis for investigating various fluid systemsin the nature and industries. In recent years, more and more researchers use pore-scale mod-els to explore various physical phenomena occurring in complex media such as hydrocarbonreservoirs. After solving the flow equations in pore-scale, this approach can be utilized tocalculate some macroscopic values using particular forms of averaged flow field.In this research, Smoothed Particle Hydrodynamics (SPH) is used to simulate two-phaseflow in pore spaces of two- and three-dimensional porous media. The SPH method is a La-grangian meshless method that is suitable for flows with free surfaces, deformable... 

In this article, we investigate the transient pressure response and production rate of horizontal wells in a dual porosity-dual permeability naturally fractured reservoir.First, the point source solution is derived in Laplace transform space and then pressure distribution and horizontal well response is obtained using the principle of superposition. In our model, despite of dual porosity model, matrix system is directly connected to horizontal wellbore as well as fracture system and has an influence on type curves characteristics. We also showed that type curves characteristics are affected by horizontal well length, well eccentricity, external boundary conditions, permeability and storavity...