## Publications

150 results found

Masihi M, Shams R, King PR, 2022, Pore level characterization of Micro-CT images using percolation theory, *Journal of Petroleum Science and Engineering*, Vol: 211, Pages: 110113-110113, ISSN: 0920-4105

Flow through porous media depends strongly on the spatial distribution of the geological heterogeneities which appear on all length scales. We lack precise information about heterogeneity distribution on various scales, from pore level to reservoir scale. However, some sources provide suitable information. At pore scale, for example, the micro-CT images show considerable insights into pore space structures and play valuable role in porous media characterization. The consequence of all geological heterogeneities is a great deal of uncertainty in dynamic performance of porous media which can be investigated using percolation theory. The main percolation quantities include the connected pore fraction; , the backbone fraction; , the dangling ends fraction; , and the effective permeability; . By finite size scaling within percolation theory, these quantities (i.e. ) become some functions of total pore fraction; , and the system dimensionless length; . In this work we examine the functional forms of percolation quantities on two-dimensional micro-CT images and develop new correlations for such quantities in three-dimensions. The results show good agreement on micro-CT samples with different pore size distribution and wide level of connectivity.

Adeyemi B, Ghanbarian B, Winter CL,
et al., 2022, Determining effective permeability at reservoir scale: Application of critical path analysis, *Advances in Water Resources*, Vol: 159, Pages: 1-12, ISSN: 0309-1708

Determining the effective permeability (keff) of geological formations has broad applications to site remediation, aquifer discharge or recharge, hydrocarbon production, and enhanced oil recovery. The objectives of this study are: (1) to explore an approach to estimating keff at the reservoir scale using the critical path analysis (CPA), (2) to evaluate the accuracy of this new approach by comparing the estimated keff to the numerically simulated effective permeability, and (3) to compare the performance of CPA estimates of keff with estimates by three other models i.e., perturbation theory (PT), effective-medium approximation (EMA), and renormalization group theory (RGT). We construct two- and three-dimensional random (uncorrelated) geologic formations based on permeability measurements from the Borden site and assume that the permeability distribution conforms to the log-normal probability density function over a wide range of means and standard deviations. Comparing keff estimated via CPA to keff values derived from numerical flow simulations indicates that CPA provides accurate estimations in both two and three dimensions over a wide range of heterogeneity levels, similar to RGT. Inter-model comparisons show that although PT and EMA provide reasonable keff estimations in rather homogeneous formations, they substantially overestimate the effective permeability in highly heterogeneous formations.

Ladipo L, Blunt MJ, King PR, 2022, Crossflow effects on low salinity displacement in stratified heterogeneity, *Journal of Petroleum Science and Engineering*, Vol: 208, Pages: 1-26, ISSN: 0920-4105

Crossflow is a major factor affecting recovery efficiency in heterogeneous permeable media. In typical water-oil displacements, viscous-dominated crossflow improves oil recovery efficiency relative to no-crossflow depending on the shock-front and/or the mobility ratio across the displacement front. Its impact is not yet fully understood for augmented or engineered waterfloods such as controlled/low salinity waterflooding (LSWF). This is critical in such a flood with two distinct displacement shock-fronts – unlike a standard waterflood – that are potentially influenced by mixing of the brines which further complicates the crossflow behaviour. This paper presents a comprehensive treatment of crossflow effects on recovery or displacement efficiency along stratified media of contrasting properties during LSWF considering physical dispersion.We define dimensionless numbers to characterize no-crossflow, viscous- and gravity-dominated crossflow regimes for different mobility-ratios. In two-dimensional numerical simulations, we explore the influence of property contrasts and mobility-ratios across the two distinct shock-fronts on the viscous crossflow behaviour in a LSWF. The sensitivity of viscous crossflow recovery and (low-salinity) engineered-water sweep efficiency to mobility-ratios is evaluated at different performance times relative to no-crossflow displacement.Viscous crossflow (VC) is found to be relevant in water-oil displacements for permeability contrasts less than or equal to 1000, but less important for EOR low-salinity displacement once the permeability contrast exceeds 50. For the mobility-ratio cases considered, VC is unidirectional – from the fast to slow layer – only when the mobility-ratios across the two distinct shock-fronts are both favourable. Unlike a typical water-oil displacement, the weak dependency of late-time recovery efficiency on VC is observed to be a function of the mobility-ratio and dispersion in LSWF. An unfavour

da Motta Pires PR, King PR, 2021, Interpolating qualitative time-lapse seismic interpretations, *SPE Journal*, Vol: 26, Pages: 3702-3718, ISSN: 1086-055X

Time-lapse seismic surveys, even if repeated at different calendar intervals, can only capture snapshots of a reservoir’s state. To obtaincontinuous information from these surveys, we propose a method to interpolate qualitative time-lapse seismic interpretations. Asinputs, we take the sequence of interpreted boundaries of an expanding time-lapse seismic anomaly. Our solution combines the fastmarching method (FMM) and the Pollock (1988)pathline-tracing algorithm to estimate the intermediate boundary positions at timesbetween acquisitions. For validation, we use a synthetic data set representing the scenario of a waterflooded sandstone reservoir. Theresults confirm that, despite small localized discrepancies in shape, at each time, the estimated boundary positions are in close agreement with the truth case. When time-lapse anomalies result primarily from changes in water saturations, the estimated boundary-arrivaltimes resemble the actual water-front-arrival times. Hence, combined with classical immiscible displacement theory, the estimatedarrival times can support a practical, model-independent, quantitative estimation of water saturations in the reservoir

Gago P, Konstantinou C, Biscontin G,
et al., 2020, Stress inhomogeneity effect on fluid-induced fracture behaviour into weakly consolidated granular systems, *Physical Review E: Statistical, Nonlinear, and Soft Matter Physics*, Vol: 102, ISSN: 1539-3755

We study the effect of stress inhomogeneity on the behavior of fluid-driven fracture development in weakly consolidated granular systems. Using numerical models we investigate the change in fracture growth rate and fracture pattern structure in unconsolidated granular packs (also referred to as soft-sands) as a function of the change in the confining stresses applied to the system. Soft-sands do not usually behave like brittle, linear elastic materials, and as a consequence, poroelastic models are often not applicable to describe their behavior. By making a distinction between “cohesive” and “compressive” grain-grain contact forces depending on their magnitude, we propose an expression that describes the fluid opening pressure as a function of the mean value and the standard deviation of the “compressive stress” distribution. We also show that the standard deviation of this distribution can be related with the extent to which fracture “branches” reach into the material.

Mittal S, Westbroek MJE, King PR,
et al., 2020, Path integral Monte Carlo method for the quantum anharmonic oscillator, *EUROPEAN JOURNAL OF PHYSICS*, Vol: 41, ISSN: 0143-0807

- Author Web Link
- Cite
- Citations: 2

Gago P, Raeini A, King P, 2020, A spatially resolved fluid-solid interaction model for dense granular packs/Soft-Sand., *Advances in Water Resources*, Vol: 136, Pages: 1-7, ISSN: 0309-1708

Fluid flow through dense granular packs or soft sands can be described as a Darcy’ s flow for low injection rates, as the friction between grain-grain and grain-walls dominate the solid system behaviour. For high injection rates, fluid forces can generate grain displacement forming flow channels or “fractures”, which in turn modify local properties within the system, such as permeability and stress distribution. Due to this kind of “self organized” behaviour, a spatially resolved model for these interactions is required to capture the dynamics of these systems. In this work, we present a resolved model based on the approach taken by the CFDEM open source project which uses LIGGGHTS – a discrete elements method (DEM)– to model the granular behaviour and OpenFoam finite volume library for computational fluid dynamics (CFD), to simulate the fluid behaviour. The capabilities provided by the DEM engine allows the properties of the solid phase, such as inter-grain cohesion and solid confinement stress to be controlled. In this work the original solver provided by the CFDEM project was modified so as to deal with dense granular packs more effectively. Advantages of the approach presented are that it does not require external “scaling parameters” to reproduce well known properties of porous materials and that it inherits the performance provided by the CFDEM project. The model is validated by reproducing the well-known properties of static porous materials, such as its permeability as a function of the system porosity, and by calculating the drag coefficient for a sphere resting inside a uniform flow. Finally, we present fracture patterns obtained when modelling water injection into a Hele-Shaw cell, filled with a dense granular pack.

Gago P, King P, Wieladek K, 2020, Fluid-induced fracture into weakly consolidated sand: Impact of confining stress on initialization pressure, *Physical Review E: Statistical, Nonlinear, and Soft Matter Physics*, Vol: 101, Pages: 012907-1-012907-6, ISSN: 1539-3755

This paper studies the process of fluid injection driven fractures in granular packs where particles are held together by external confining stresses and weak intergrain cohesion. We investigate the process of fracture formations in soft sand confined into a radial Hele-Shaw cell. Two main regimes are well known for fluid injection in soft sand. For low fluid injection pressures it behaves as a solid porous material while for high enough injection pressures grain rearrangement takes place. Grain rearrangements lead to the formation of fluid channels or “fractures,” the structure and geometry of which depend on the material and fluid properties. Due to macroscopic grain displacements and the predominant role of dissipative frictional forces in granular system dynamics, these materials do not behave as conventional brittle, linear elastic materials and the transition between these two regimes cannot usually be described using poroelastic models. In this work we investigate the change in the minimum fluid pressure required to start grain mobilization as a function of the confining stresses applied to the system using a spatially resolved computational fluid dynamics–discrete element method numerical model. We show that this change is proportional to the applied stress when the confining stresses can be regarded as uniformly distributed among the particles in the system. A preliminary analytical expression for this change is presented.

Ladipo L, Blunt MJ, King PR, 2020, A salinity cut-off method to control numerical dispersion in low-salinity waterflooding simulation, *Journal of Petroleum Science and Engineering*, Vol: 184, Pages: 1-19, ISSN: 0920-4105

Low-salinity or controlled salinity waterflooding (LSWF) is a promising enhanced oil recovery (EOR) technique. In simulations of this process, numerical dispersion smears saturation fronts, causing errors in the results. The objective of this work is to control these effects in LSWF simulation. We examine the impact of numerical dispersion on simulated LSWF performance. The low-salinity (LS) front is smeared even at unfeasibly fine grids. The velocities of the water fronts are altered. A numerical mixed zone forms around the interface between the injected and resident brines. This mimics a typical physical mixing effect.In reservoir simulation, threshold salinities are defined where the low salinity effect (LSE) is first encountered. It has been suggested that numerical dispersion effects can be corrected by imposing effective thresholds. We demonstrate that existing methods to evaluate these effective salinities do not accurately predict the salt front movement especially when dispersion is significant.We propose a simple simulation-based approach to evaluate the effective salinities based on the conservation of volumes of the resident and injected brines in the reference and upscaled solutions. After comparing analytical and corrected coarse-grid solutions in one-dimension, the effectiveness of the approach is demonstrated in multi-dimensional systems.A method is proposed to control the numerical mixed zone to replicate a physical longitudinal mixing effect. This method is demonstrated in one-dimension and does not require a fine-grid numerical solution as a benchmark. We investigate the effects of effective thresholds on the modeled transverse mixing or dispersion. A method to model transverse dispersion in simulations with an effective longitudinal component is suggested. This method is extended to the explicit modeling of physical dispersion in systems with transverse flows.We can now simply evaluate the effective salinities for a simulation grid; and control i

Ladipo L, Blunt M, King P, 2020, Optimizing low salinity waterflooding with controlled numerical influence of physical mixing considering uncertainty

Controlled/Low Salinity Waterflooding (LSWF) is an augmented waterflood with well-reported improved displacement efficiency compared with conventional waterfloods. Physical mixing or dispersion of the injected low-salinity (LS) brine with the formation high-salinity (HS) brine substantially reduces the low-salinity effect. Numerical dispersion often misrepresents this mixing in conventional LSWF-simulations, causing errors in the results. Uncertainty in the reservoir description further makes the evaluated performance questionable. Existing studies have suggested optimal amounts for the injected LS-brine to sustain its displacement stability during interwell flows with physical mixing, but with poor or no consideration of uncertainty. This work focuses on optimizing the injected LS-brine amount considering reported flow uncertainties while ensuring adequate correction of the erroneous influence of numerical dispersion on physical mixing. We investigate the impacts of flow uncertainties on the optimal LS slug-size. The sensitivity of the optimal slug-size to heterogeneity is examined under uncertainty. We evaluate how the interaction between physical mixing and geological heterogeneity influences slug integrity and performance. We propose an improved 'effective salinities' concept to evaluate appropriate effective salinities to characterize the desired representative physical mixing supressing the large numerical dispersion effects usually encountered in coarse-grid LSWF-simulations. This ensures reliable representation of physical dispersion in such grids. We consider different models with characterized levels of heterogeneity and essential variables that control the impact of mixing on LSWF performance based mainly on reported data. New indicators are defined to evaluate the displacement stability and performance of injected LS-brine thereby relating its technical and economic performance. Slug performance is evaluated at different injection times to examine the se

Shokrollahzadeh Behbahani S, Masihi M, Ghazanfari MH,
et al., 2019, Effect of Characteristic Time on Scaling of Breakthrough Time Distribution for Two-Phase Displacement in Percolation Porous Media, *Transport in Porous Media*, Vol: 130, Pages: 889-902, ISSN: 0169-3913

Determining the time of breakthrough of injected water is important when assessing waterflood in an oil reservoir. Breakthrough time distribution for a passive tracer (for example water) in percolation porous media (near the percolation threshold) gives insights into the dynamic behavior of flow in geometrically complex systems. However, the application of such distribution to realistic two-phase displacements can be done based on scaling of all parameters. Here, we propose two new approaches for scaling of breakthrough time (characteristic times) in two-dimensional flow through percolation porous media. The first is based on the flow geometry, and the second uses the flow parameters of a representative homogenous model. We have tested the effectiveness of these two approaches using a large number of dynamic simulations. The results show significant improved distribution curves for the breakthrough (transit) time between an injector and a producer located in a heterogeneous porous medium in comparison with the previous scaling methods.

Zhou Y, Muggeridge AH, Berg CF,
et al., 2019, Effect of Layering on Incremental Oil Recovery From Tertiary Polymer Flooding, *SPE RESERVOIR EVALUATION & ENGINEERING*, Vol: 22, Pages: 941-951, ISSN: 1094-6470

- Author Web Link
- Cite
- Citations: 6

Westbroek MJE, Coche G-A, King PR,
et al., 2019, Pressure statistics from the path integral for Darcy flow through random porous media, *Journal of Physics A: Mathematical and Theoretical*, Vol: 52, ISSN: 1751-8113

The path integral for classical statistical dynamics is used to determine the properties of one-dimensional Darcy flow through a porous medium with a correlated stochastic permeability for several spatial correlation lengths. Pressure statistics are obtained from the numerical evaluation of the path integral by using the Markov chain Monte Carlo method. Comparisons between these pressure distributions and those calculated from the classic finite-volume method for the corresponding stochastic differential equation show excellent agreement for Dirichlet and Neumann boundary conditions. The evaluation of the variance of the pressure based on a continuum description of the medium provides an estimate of the effects of discretization. Log-normal and Gaussian fits to the pressure distributions as a function of position within the porous medium are discussed in relation to the spatial extent of the correlations of the permeability fluctuations.

Westbroek MJE, King PR, Vvedensky DD, et al., 2019, Pressure and flow statistics of Darcy flow from simulated annealing, Publisher: arXiv

The pressure and flow statistics of Darcy flow through a random permeablemedium are expressed in a form suitable for evaluation by the method ofsimulated annealing. There are several attractive aspects to using simulatedannealing: (i) any probability distribution can be used for the permeability,(ii) there is no need to invert the transmissibility matrix which, while not afactor for single-phase flow, offers distinct advantages for the case ofmultiphase flow, and (iii) the action used for simulated annealing is eminentlysuitable for coarse graining by integrating over the short-wavelength degreesof freedom. In this paper, we show that the pressure and flow statisticsobtained by simulated annealing are in excellent agreement with the moreconventional finite-volume calculations.

King PR, Masihi M, 2018, Percolation theory in reservoir engineering, ISBN: 9781786345233

This book aims to develop the ideas from fundamentals of percolation theory to practical reservoir engineering applications. Through a focus on field scale applications of percolation concepts to reservoir engineering problems, it offers an approximation method to determine many important reservoir parameters, such as effective permeability and reservoir connectivity and the physical analysis of some reservoir engineering properties. Starring with the concept of percolation theory, it then develops into methods to simple geological systems like sand-bodies and fractures. The accuracy and efficiency of the percolation concept for these is explained and further extended to more complex realistic models. Percolation Theory in Reservoir Engineering primarily focuses on larger reservoir scale flow and demonstrates methods that can be used to estimate large scale properties and their uncertainty, crucial for major development and investment decisions in hydrocarbon recovery.

Gago PA, King P, Muggeridge A, 2018, Fractal growth model for estimating breakthrough time and sweep efficiency when waterflooding geologically heterogeneous rocks, *Physical Review Applied*, Vol: 10, ISSN: 2331-7019

We describe a fast method for estimating flow through a porous medium with a heterogeneous permeability distribution. The main application is to contaminant transport in aquifers and recovery of oil by waterflooding, where such geological heterogeneities can result in regions of bypassed contaminants or oil. The extent of this bypassing is normally assessed by a numerical flow simulation that can take many hours of computer time. Ideally the impact of uncertainty in the geological description is then evaluated by the performing of many such simulations using different realizations of the permeability distribution. Obviously, a proper Monte Carlo evaluation may be impossible when the flow simulations are so computationally intensive. Consequently, methods from statistical mechanics, such as percolation theory and random walkers (such as diffusion-limited aggregation), have been proposed; however, these methods are limited to geological heterogeneities where the correlation lengths are smaller than the system size or to continuous permeability distributions. Here we describe a growth model that can be used to estimate the breakthrough time of the water (and hence the sweep efficiency) in most types of geologically heterogeneous rocks. We show how the model gives good estimates of the breakthrough time of water at the production well in a fraction of the time needed to perform a full flow simulation.

Westbroek MJE, Coche G-A, King PR,
et al., 2018, Evaluation of the path integral for flow through random porous media, *Physical Review E: Statistical, Nonlinear, and Soft Matter Physics*, Vol: 97, Pages: 1-5, ISSN: 1539-3755

We present a path integral formulation of Darcy's equation in one dimension with random permeability described by a correlated multivariate lognormal distribution. This path integral is evaluated with the Markov chain Monte Carlo method to obtain pressure distributions, which are shown to agree with the solutions of the corresponding stochastic differential equation for Dirichlet and Neumann boundary conditions. The extension of our approach to flow through random media in two and three dimensions is discussed.

Westbroek MJE, King PR, Vvedensky DD,
et al., 2018, User's guide to Monte Carlo methods for evaluating path integrals, *American Journal of Physics*, Vol: 86, Pages: 293-304, ISSN: 0002-9505

We give an introduction to the calculation of path integrals on a lattice, with the quantum harmonic oscillator as an example. In addition to providing an explicit computational setup and corresponding pseudocode, we pay particular attention to the existence of autocorrelations and the calculation of reliable errors. The over-relaxation technique is presented as a way to counter strong autocorrelations. The simulation methods can be extended to compute observables for path integrals in other settings.

Al-Shamma BR, Gosselin O, King PR, 2018, History matching using hybrid parameterisation and optimisation methods, 80th EAGE Conference and Exhibition 2018

Copyright © 2018, Society of Petroleum Engineers. Reservoir models are commonly used in the oil and gas industry to predict reservoir behaviour and forecast production in order to make important financial decision such as reserves estimations, infill well drilling, enhanced oil recovery schemes, etc. Conditioning reservoir models to dynamic production data is known as history matching, which is usually carried out in an attempt to enhance the predicted reservoir performance. Uncertainty quantification is also an important aspect of this task, and encompasses identifying multiple history matched models, which are constrained to a geological concept. History matching and uncertainty quantification can be accomplished by identifying and using efficient and speedy optimisation techniques. The assisted history matching practice usually includes two practices; the first of which is parameterisation, which consists of reducing the number of matching parameters in order to avoid adjusting too many variables with respect to the amount of production data available. A challenging situation results from over-parameterisation, in addition to an ill-posed formulation of the inverse problem. The second process involves optimisation, which aims at solving the inverse problem by reducing a misfit or objective function that defines the difference between simulated and production data. The main challenges of optimisation are local minima solutions and premature convergence. The success of optimisation is greatly dependent on the parameterisation strategy used. These algorithms that analyse various parameterisation methods, combined and examined with diverse optimisation algorithms lead us to suggest novel hybrid approaches addressing the two processes of assisted history matching. We propose a multistage combined parameterisation and optimisation history matching technique. Hybridisation of parameterisation and optimisation algorithms when designed in an optimum manner can combin

Al-Dhuwaihi A, King P, Muggeridge A, 2018, Upscaling for polymer flooding

Polymer flooding is a proven EOR/IOR process for viscous and light oil reservoirs alike. However, it results in the formation of two shocks front that require simulation models with fine grid blocks to represent field scale fluid movement. Therefore, upscaling is required to transfer such fluid behavior to coarser models. However, most upscaling methods are designed for waterflood only, while upscaling techniques for polymer flood are rarely discussed in the literature. In this paper, A new upscaling methodology specifically designed for polymer flooding is presented to address such impracticality. The methodology allows the average flow behavior to be captured, including the effects of small scale heterogeneity whilst compensating for the impact of increased numerical diffusion present in coarse grid models. The method is based on the pore volume weighted method for relative permeability pseudoization first derived by Emanuel and Cook (1974) for waterflooding but extends its implementation to model polymer specific parameters such as adsorption isotherm and viscosity-concentration function. The method is demonstrated on a series of simple reservoir models for a range of different aggregation ratios, showing overall improvement in the prediction of oil recovery, water cut, produced polymer concentration with time, and pressure response in the coarse grid models. This is demonstrated by comparing the predictions from the coarse grid, upscaled models with those from fine grid simulations and coarse grid simulations of the same model reservoir without the new upscaling methodology.

Al-Shamma BR, Gosselin O, King PR, 2018, History matching using hybrid parameterisation and optimisation methods

Reservoir models are commonly used in the oil and gas industry to predict reservoir behaviour and forecast production in order to make important financial decision such as reserves estimations, infill well drilling, enhanced oil recovery schemes, etc. Conditioning reservoir models to dynamic production data is known as history matching, which is usually carried out in an attempt to enhance the predicted reservoir performance. Uncertainty quantification is also an important aspect of this task, and encompasses identifying multiple history matched models, which are constrained to a geological concept. History matching and uncertainty quantification can be accomplished by identifying and using efficient and speedy optimisation techniques. The assisted history matching practice usually includes two practices; the first of which is parameterisation, which consists of reducing the number of matching parameters in order to avoid adjusting too many variables with respect to the amount of production data available. A challenging situation results from over-parameterisation, in addition to an ill-posed formulation of the inverse problem. The second process involves optimisation, which aims at solving the inverse problem by reducing a misfit or objective function that defines the difference between simulated and production data. The main challenges of optimisation are local minima solutions and premature convergence. The success of optimisation is greatly dependent on the parameterisation strategy used. These algorithms that analyse various parameterisation methods, combined and examined with diverse optimisation algorithms lead us to suggest novel hybrid approaches addressing the two processes of assisted history matching. We propose a multistage combined parameterisation and optimisation history matching technique. Hybridisation of parameterisation and optimisation algorithms when designed in an optimum manner can combine advantageous features of each method. This consisted

Masihi M, Gago P, King P, 2016, Estimation of the Effective Permeability of Heterogeneous Porous Media by Using Percolation Concepts, *Transport in Porous Media*, Vol: 114, Pages: 169-199, ISSN: 1573-1634

In this paper we present new methods to estimate the effective permeability (k_eff) of heterogeneous porous media with a wide distribution of permeabilities and various underlying structures, using percolation concepts. We first set a threshold permeability (k_th) on the permeability density function (pdf) and use standard algorithms from percolation theory to check whether the high permeable grid blocks (i.e. those with permeability higher than k_th) with occupied fraction of “p” first forms a cluster connecting two opposite sides of the system in the direction of the flow (high permeability flow pathway). Then we estimate the effective permeability of the heterogeneous porous media in different ways: a power law (k_eff=k_th p^m), a weighted power average (k_eff=[p.k_th^m+(1-p).k_g^m ]^(1/m) with k_g the geometric average of the permeability distribution) and a characteristic shape factor multiplied by the permeability threshold value. We found that the characteristic parameters (i.e. the exponent “m”) can be inferred either from the statistics and properties of percolation sub-networks at the threshold point (i.e. high and low permeable regions corresponding to those permeabilities above and below the threshold permeability value) or by comparing the system properties with an uncorrelated random field having the same permeability distribution. These physically based approaches do not need fitting to the experimental data of effective permeability measurements to estimate the model parameter (i.e. exponent m) as is usually necessary in empirical methods. We examine the order of accuracy of these methods on different layers of 10th SPE model and found very good estimates as compared to the values determined from the commercial flow simulators.

Westbroek MJE, King PR, Vvedensky DD, 2016, Path Integral Method for Flow through Random Porous Media

One of the key problems in modelling flow in oil reservoirs is our lack of precise knowledge of the variations in flow properties across the field. At best we can infer the statistics of these variations from field observations. The challenge is to determine the statistics of the flow itself (flow rates, pressures etc.) from the statistics of the permeability variations. Conventional simulations are computationally very expensive unless smart sampling techniques or surrogate models are used. In this paper we demonstrate the use of a path integral formulation for this problem. To demonstrate how this methods works, we start with the one dimensional Darcy flow problem: q(x)=-K(x)dp(x)/dx where p(x) is the pressure, q(x) is the flow rate and K(x) is the rock permeability. The randomness of the porous medium is modelled by regarding K as a stochastic quantity which is assumed to follow Gaussian statistics. Because of the randomly varying rock structure, there is a variety of conceivable pressure realisations p(x). The path integral Z is an integral over all realisations with an appropriate probability measure. Once Z is evaluated, either analytically, or by standard Monte Carlo methods, any observable of interest, including pressure correlations can be easily obtained.

Westbroek MJE, King PR, Vvedensky DD, 2016, Path Integral Method for Flow through Random Porous Media

One of the key problems in modelling flow in oil reservoirs is our lack of precise knowledge of the variations in flow properties across the field. At best we can infer the statistics of these variations from field observations. The challenge is to determine the statistics of the flow itself (flow rates, pressures etc.) from the statistics of the permeability variations. Conventional simulations are computationally very expensive unless smart sampling techniques or surrogate models are used. In this paper we demonstrate the use of a path integral formulation for this problem. To demonstrate how this methods works, we start with the one dimensional Darcy flow problem: q(x)=-K(x)dp(x)/dx where p(x) is the pressure, q(x) is the flow rate and K(x) is the rock permeability. The randomness of the porous medium is modelled by regarding K as a stochastic quantity which is assumed to follow Gaussian statistics. Because of the randomly varying rock structure, there is a variety of conceivable pressure realisations p(x). The path integral Z is an integral over all realisations with an appropriate probability measure. Once Z is evaluated, either analytically, or by standard Monte Carlo methods, any observable of interest, including pressure correlations can be easily obtained.

Sadeghnejad S, Masihi M, King PR, et al., 2016, Study the effect of connectivity between two wells on secondary recovery efficiency using percolation approach

Estimating available hydrocarbon to be produced during secondary oil recovery is an ongoing activity in field development. The primary plan is normally scheduled during early stage of field's life through master development plan studies. During this period, due to the lake of certain data, estimation of the field efficiency is usually based on rules of thumb and not detailed field characterization. Hence, there is a great motivation to produce simpler physically-based methodologies. The minimum necessity inputs of percolation approach make it a useful tool for foration performance prediction. This approach enables us to attain a better assessment of the efficiency of secondary recovery methods at early production time. The main contribution of this study is to establish a continuum percolation model based on Monte Carlo simulation that can estimate the connectivity of good sands between two wells. In the classical percolation, the connectivity is considered between two lines and two faces of the system in 2-And 3-D; whereas, hydrocarbon production is achieved through wells with the shape of lines (e.g., vertical, horizontal, or deviated wells). In addition, the results showed that not implementation of the correct geometry of wells can alter the estimated results from the percolation approach.

Sadeghnejad S, Masihi M, King PR, 2016, Study the connectivity of good sands between two wells represented by two points using percolation theory

One of the major applications of percolation theory in petroleum engineering is investigation of connectivity in complex formations. Production normally is achieved through a heterogeneous porous media. Proper assessment of connectivity of formation considering its heterogeneity is important in formation evaluation. Percolation assumes that heterogeneity can be simplified to either permeable or impermeable rocktypes. Considering this, the system outcome (e.g., prediction of recovery) can be easily described by simple mathematical relationships which are entirely independent of small-scale details of formation. The main contribution of this work is to use a continuum percolation approach to estimate the connectivity of permeable sands via two points (P2P) representing two injection and production wells and comparing the results with the conventional line to line connectivity (L2L) in previous studies. In particular, the percolation exponents will be investigated both in P2P and L2L and their connectivity curves will be compared. For this purpose, an object-based technique based on Monte-Carlo simulation is used to model the spatial distribution of isotropic sandbodies in 2-D. The results showed that proper modelling of the shape of wells is a critical issue that can alter the obtained results associated with the amount of connected hydrocarbon when one uses the percolation approach.

Gago PA, King PR, Muggeridge AH, 2016, Fast estimation of effective permeability and sweep efficiency of waterflooding in geologically heterogeneous reservoirs

Geological heterogeneity can adversely affect the macroscopic sweep efficiency when waterflooding oil reservoirs, however the exact distribution of permeability and porosity is generally not known. Engineers try to estimate the range of impacts heterogeneity might have on waterflood efficiency by creating multiple geological models and then simulating a waterflood through each of those realizations. Unfortunately each simulation can be computationally intensive meaning that it is generally not possible to obtain a statistically valid estimate of the expected sweep and the associated standard deviation. In this paper we show how the volume of unswept oil can be estimated rapidly (without flow simulations) from a geometrical characterization of the spatial permeability distribution. A "constriction" factor is defined which quantifies the effective cross-section area of the zones perpendicular to the principal flow direction. This is combined with a 'net-To-gross ratio' (which quantifies the fractional reservoir volume occupied by the zones that contribute to flow) to estimate effective permeability and the expected recovery factor for that realization. The method is tested using a range of realistic geological models, including SPE10 model 2 and its predictions are shown to agree well with values obtained using a well established commercial flow simulator.

Sadeghnejad S, Masihi M, King PR, 2016, Study the connectivity of good sands between two wells represented by two points using percolation theory

One of the major applications of percolation theory in petroleum engineering is investigation of connectivity in complex formations. Production normally is achieved through a heterogeneous porous media. Proper assessment of connectivity of formation considering its heterogeneity is important in formation evaluation. Percolation assumes that heterogeneity can be simplified to either permeable or impermeable rocktypes. Considering this, the system outcome (e.g., prediction of recovery) can be easily described by simple mathematical relationships which are entirely independent of small-scale details of formation. The main contribution of this work is to use a continuum percolation approach to estimate the connectivity of permeable sands via two points (P2P) representing two injection and production wells and comparing the results with the conventional line to line connectivity (L2L) in previous studies. In particular, the percolation exponents will be investigated both in P2P and L2L and their connectivity curves will be compared. For this purpose, an object-based technique based on Monte-Carlo simulation is used to model the spatial distribution of isotropic sandbodies in 2-D. The results showed that proper modelling of the shape of wells is a critical issue that can alter the obtained results associated with the amount of connected hydrocarbon when one uses the percolation approach.

Gago PA, King PR, Muggeridge AH, 2016, Fast estimation of effective permeability and sweep efficiency of waterflooding in geologically heterogeneous reservoirs

Geological heterogeneity can adversely affect the macroscopic sweep efficiency when waterflooding oil reservoirs, however the exact distribution of permeability and porosity is generally not known. Engineers try to estimate the range of impacts heterogeneity might have on waterflood efficiency by creating multiple geological models and then simulating a waterflood through each of those realizations. Unfortunately each simulation can be computationally intensive meaning that it is generally not possible to obtain a statistically valid estimate of the expected sweep and the associated standard deviation. In this paper we show how the volume of unswept oil can be estimated rapidly (without flow simulations) from a geometrical characterization of the spatial permeability distribution. A "constriction" factor is defined which quantifies the effective cross-section area of the zones perpendicular to the principal flow direction. This is combined with a 'net-To-gross ratio' (which quantifies the fractional reservoir volume occupied by the zones that contribute to flow) to estimate effective permeability and the expected recovery factor for that realization. The method is tested using a range of realistic geological models, including SPE10 model 2 and its predictions are shown to agree well with values obtained using a well established commercial flow simulator.

Sadeghnejad S, Masihi M, King PR, et al., 2016, Study the effect of connectivity between two wells on secondary recovery efficiency using percolation approach

Estimating available hydrocarbon to be produced during secondary oil recovery is an ongoing activity in field development. The primary plan is normally scheduled during early stage of field's life through master development plan studies. During this period, due to the lake of certain data, estimation of the field efficiency is usually based on rules of thumb and not detailed field characterization. Hence, there is a great motivation to produce simpler physically-based methodologies. The minimum necessity inputs of percolation approach make it a useful tool for foration performance prediction. This approach enables us to attain a better assessment of the efficiency of secondary recovery methods at early production time. The main contribution of this study is to establish a continuum percolation model based on Monte Carlo simulation that can estimate the connectivity of good sands between two wells. In the classical percolation, the connectivity is considered between two lines and two faces of the system in 2-And 3-D; whereas, hydrocarbon production is achieved through wells with the shape of lines (e.g., vertical, horizontal, or deviated wells). In addition, the results showed that not implementation of the correct geometry of wells can alter the estimated results from the percolation approach.

This data is extracted from the Web of Science and reproduced under a licence from Thomson Reuters. You may not copy or re-distribute this data in whole or in part without the written consent of the Science business of Thomson Reuters.