The Use of Capillary Tube Networks in Reservoir Performance Studies: 1. Equal-Viscosity Miscible Displacements

1971 ◽  
Vol 11 (02) ◽  
pp. 99-112 ◽  
Author(s):  
Ralph Simon ◽  
F.J. Kelsey

Abstract This paper concerns the use of network principles to study displacement phenomena in porous media. The information presented is for equal-viscosity, equal-density miscible displacements. The paper explains the reasons for using an interconnected network of capillary tubes to model the interconnected network of pores in a reservoir rock. A method is presented for defining the heterogeneity of a presented for defining the heterogeneity of a network of tubes based on tube-size and tube-location distribution functions. A technique is described for constructing a network whose heterogeneity models the heterogeneity of pores in a reservoir rock. The use of networks to provide information which can be used in the solution of reservoir engineering problems is illustrated with example calculations of the effect of heterogeneity on fingering, breakthrough, and selective plugging in linear systems, and the effect of heterogeneity on areal sweep efficiency in a five-spot pattern. Introduction Oil in a reservoir is contained in an interconnected three-dimensional network of pores. Direct evidence of the nature of this network of pores comes from examination of petrographic thin sections and three dimensional Scanning Electron Microscope (SEM) pictures of the pores. The SEM pictures show that the pores in a reservoir rock are channels through which flow can occur. These channels have highly irregular configurations so irregular that it is not practical at this time to calculate flow behavior through individual channels or through the interconnected network of the channels. It is practical, however, to use a computer to calculate flow behavior in an interconnected network of capillary tubes and several investigators have studied the problem of using a network of tubes to model a network of pores. pores. Fatt pioneered the idea of using a network of cubes model for reservoir engineering studies. He demonstrated that capillary pressures, relative permeabilities, and electrical resistivities permeabilities, and electrical resistivities calculated for a network model have the same characteristics as those measured for real pores in reservoir rocks. From this, Fatt concluded that the network of tubes is a valid model of real porous media. Rose reinforced Fatt's conclusion and showed that computers can be used to study the displacement characteristics of networks and to obtain results "…which can be supposed to have a direct bearing on the mechanics of petroleum recovery…" This paper takes two steps beyond the work of Fatt and Rose. First, it describes a technique for constructing a network whose heterogeneity models the heterogeneity of natural pores. This is done by matching calculated equal-viscosity miscible displacement behavior in the network with measured behavior in a laboratory core. Second, it illustrates the use of the network model for calculating the effects of heterogeneity on fingering, breakthrough, and plugging in linear systems and areal sweep efficiency in a five-spot pattern. The networks used in the studies in this paper consist of several hundred interconnected capillary tubes of different sizes. Four different types of connections or configurations were investigated and are shown below. These configurations are discussed in detail later in the paper. SPEJ P. 99

1968 ◽  
Vol 8 (01) ◽  
pp. 52-62 ◽  
Author(s):  
K.S. Lee ◽  
E.L. Claridge

Abstract Areal sweep efficiency of oil displacement by enhanced-viscosity water exhibiting pseudoplastic behavior was measured in a Hele-Shaw model representing one-quarter of a five-spot pattern. The pseudoplasticity of polymer solutions and the velocity distribution in the five-spot pattern produced a condition under which the mobility ratio between the displacing and the displaced fluid could not be assigned a single value. Instead, the movement of the displacement front is governed by local mobility ratios which are also time dependent. The areal sweep at breakthrough with polymer solutions was poorer than the sweep obtained with Newtonian fluids of comparable viscosity. However, the areal sweep and 1 PV throughput was greatly improved as compared to flood water without polymer. It was also demonstrated that, even after the oil-cut had declined to a low value during a regular waterflood, switching to polymer flood efficiently swept out the oil remaining in the model. Introduction The behavior of fluid displacements in isotropic porous media for various patterns of injection and production wells has been extensively investigated. These investigations all concerned Newtonian fluids, i.e., the viscosity of each fluid was constant regardless of flow rate. The generally unfavorable influence on areal sweep efficiency of higher mobility of the displacing fluid as compared to the mobility of the displaced fluid has been established for both miscible and immiscible fluids. The principle was also established that a close correspondence exists between miscible and immiscible flood front behavior, although oil recovery in a waterflood at unfavorable mobility ratio may be less than that observed in a miscible displacement at the same mobility ratio. This is true even when oil recovery is expressed on the basis of movable oil. The reason is that oil saturation only slowly achieves its final value behind the waterflood front in accordance with the Buckley-Leverett simultaneous flow relations. It is convenient to use miscible displacements for laboratory simulation of waterflood frontal advance since the interfacial tension forces which are negligible in proportion to viscous forces on a reservoir scales are thus made nonoperative in the laboratory model. For miscible displacements, the Hele-Shaw type of model adequately represents a porous medium so long as the appropriate scaling rules are observed in its design and operation. During simulation of waterflood front behavior in the laboratory by using miscible displacements, the behavior of connate water may ordinarily be disregarded since it is usually indistinguishable from flood water in this process. However, when the flood water is deliberately thickened to improve the mobility ratio between water and oil, the effect on the sweep efficiency due to generation of a connate water bank during the process must be considered. In a uniform porous medium, such a bank is generated and efficiently displaced by injection of thickened water. The oil originally in-place at the start of the waterflood is then displaced by connate water followed by thickened water. If the flood water must be thickened to obtain a favorable mobility ratio, the mobility of the oil phase is appreciably less than that of the connate water. Hence, the oil phase is inefficiently displaced by the connate water bank, and a considerable proportion of the oil comes in contact with and is displaced by the thickened waterflood front. SPEJ P. 52ˆ


Author(s):  
Jiehai Zhang ◽  
Arun Muley ◽  
Joseph B. Borghese ◽  
Raj M. Manglik

Enhanced heat transfer characteristics of low Reynolds number airflows in three-dimensional sinusoidal wavy plate-fin channels are investigated. For the computational simulation, steady state, constant property, periodically developed, laminar forced convection is considered with the channel surface at the uniform heat flux condition; the wavy-fin is modeled by its two asymptotic limits of 100% and zero fin efficiency. The governing equations are solved numerically using finite-volume techniques for a non-orthogonal, non-staggered grid. Computational results for velocity and temperature distribution, isothermal Fanning friction factor f and Colburn factor j are presented for airflow rates in the range of 10 ≤ Re ≤ 1500. The numerical results are further compared with experimental data, with excellent agreement, for two different wavy-fin geometries. The influence of fin density on the flow behavior and the enhanced convection heat transfer are highlighted. Depending on the flow rate, a complex flow structure is observed, which is characterized by the generation, spatial growth and dissipation of vortices in the trough region of the wavy channel. The thermal boundary layers on the fin surface are periodically disrupted, resulting in high local heat fluxes. The overall heat transfer performance is improved considerably, compared to the straight channel with the same cross-section, with a relatively smaller increase in the associated pressure drop penalty.


Author(s):  
Zhigang Li ◽  
Jun Li ◽  
Liming Song ◽  
Qing Gao ◽  
Xin Yan ◽  
...  

The modern gas turbine is widely applied in the aviation propulsion and power generation. The rim seal is usually designed at the periphery of the wheel-space and prevented the hot gas ingestion in modern gas turbines. The high sealing effectiveness of rim seal can improve the aerodynamic performance of gas turbines and avoid of the disc overheating. Effect of outer fin axial gap of radial rim seal on the sealing effectiveness and fluid dynamics was numerically investigated in this work. The sealing effectiveness and fluid dynamics of radial rim seal with three different outer fin axial gaps was conducted at different coolant flow rates using three-dimensional Reynolds-Averaged Navier-Stokes (RANS) and SST turbulent model solutions. The accuracy of the presented numerical approach for the prediction of the sealing performance of the turbine rim seal was demonstrated. The obtained results show that the sealing effectiveness of radial rim seal increases with increase of coolant flow rate at the fixed axial outer fin gap. The sealing effectiveness increases with decrease of the axial outer fin gap at the fixed coolant flow rate. Furthermore, at the fixed coolant flow rate, the hot gas ingestion increases with the increase of the axial outer fin gap. This flow behavior intensifies the interaction between the hot gas and coolant flow at the clearance of radial rim seal. The preswirl coefficient in the wheel-space cavity is also illustrated to analyze the flow dynamics of radial rim seal at different axial outer fin gaps.


2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Zhenglun Alan Wei ◽  
Zhongquan Charlie Zheng ◽  
Xiaofan Yang

A parallel implementation of an immersed-boundary (IB) method is presented for low Reynolds number flow simulations in a representative elementary volume (REV) of porous media that are composed of a periodic array of regularly arranged structures. The material of the structure in the REV can be solid (impermeable) or microporous (permeable). Flows both outside and inside the microporous media are computed simultaneously by using an IB method to solve a combination of the Navier–Stokes equation (outside the microporous medium) and the Zwikker–Kosten equation (inside the microporous medium). The numerical simulation is firstly validated using flow through the REVs of impermeable structures, including square rods, circular rods, cubes, and spheres. The resultant pressure gradient over the REVs is compared with analytical solutions of the Ergun equation or Darcy–Forchheimer law. The good agreements demonstrate the validity of the numerical method to simulate the macroscopic flow behavior in porous media. In addition, with the assistance of a scientific parallel computational library, PETSc, good parallel performances are achieved. Finally, the IB method is extended to simulate species transport by coupling with the REV flow simulation. The species sorption behaviors in an REV with impermeable/solid and permeable/microporous materials are then studied.


2005 ◽  
Vol 2005 (1) ◽  
pp. 77-89 ◽  
Author(s):  
W. Chon ◽  
R. S. Amano

When the airflow patterns inside a lawn mower deck are understood, the deck can be redesigned to be efficient and have an increased cutting ability. To learn more, a combination of computational and experimental studies was performed to investigate the effects of blade and housing designs on a flow pattern inside a1.1mwide corotating double-spindle lawn mower deck with side discharge. For the experimental portion of the study, air velocities inside the deck were measured using a laser Doppler velocimetry (LDV) system. A high-speed video camera was used to observe the flow pattern. Furthermore, noise levels were measured using a sound level meter. For the computational fluid dynamics (CFD) work, several arbitrary radial sections of a two-dimensional blade were selected to study flow computations. A three-dimensional, full deck model was also developed for realistic flow analysis. The computational results were then compared with the experimental results.


2019 ◽  
Vol 18 (02) ◽  
pp. 1850019
Author(s):  
Huiyuan Yu ◽  
Jiayi Zhu ◽  
Hongbo Ren ◽  
Shuxin Liu

Graphene-based aerogels with a three-dimensional interconnected network were fabricated via the hydrothermal self-assembly and thermal-annealing process. The aerogels were characterized by means of scanning electron microscopy and atomic absorption spectroscopy. The graphene-based aerogels showed highly porous structure and adsorption capacity for heavy metal ions. Thus, they would be the promising materials for removal of heavy metal ions from water.


2004 ◽  
Vol 126 (4) ◽  
pp. 692-699 ◽  
Author(s):  
Xiufang Gao ◽  
Bengt Sunde´n

The flow behavior in rib-roughened ducts is influenced by the inclination of ribs and the effect is investigated in the present study by Particle Image Velocimetry (PIV). The local flow structures between two adjacent ribs were measured. The Reynolds number was fixed at 5800. The flow field description was based on the PIV results in planes both parallel and perpendicular to the ribbed walls at various locations. The rib angle to the main flow direction was varied as 30 deg, 45 deg, 60 deg and 90 deg. The ribs induce three dimensional flow fields. The flow separation and reattachment between adjacent ribs are clearly observed. In addition, the inclined ribs are found to alter the spanwise distribution of the streamwise velocity component. The streamwise velocity component has its highest values at the upstream end of the ribs, and decreases continuously to its lowest values at the downstream end. Strong secondary flow motion occurs over the entire duct cross section for the inclined ribs. The flow structures between two consecutive ribs show that the fluid flows along the ribs from one end of the ribs to the other end, and then turns back at the transverse center. Downwash and upwash flows are observed at the upstream end and downstream end of the ribs, respectively.


Author(s):  
Ronghua Chen ◽  
Lie Chen ◽  
Wenxi Tian ◽  
Guanghui Su ◽  
Suizheng Qiu

In the typical boiling water reactor (BWR), each control rod guide tube supports four fuel assemblies via an orificed fuel support piece in which a channel is designed to be a potential corium relocation path from the core region to the lower head under severe accident conditions. In this study, the improved Moving Particle Semi-implicit (MPS) method was adopted to analyze the melt flow and ablation behavior in this region during a severe accident of BWR. A three-dimensional particle configuration was constructed for analyzing the melt flow behavior within the fuel support piece. Considering the symmetry of the fuel support piece, only one fourth of the fuel support was simulated. The eutectic reaction between Zr (the material of the corium) and stainless steel (the material of the fuel support piece) was taken into consideration. The typical melt flow and freezing behaviors within the fuel support piece were successfully reproduced by MPS method. In all the simulation cases, the melt discharged from the hole of the fuel support piece instead of plugging the fuel support piece. The results indicate that MPS method has the capacity to analyze the melt flow and solidification behavior in the fuel support piece.


Sign in / Sign up

Export Citation Format

Share Document