Results
- Showing results for:
- Reset all filters
Search results
-
Journal articleJacquemyn C, Jackson MD, Hampson GJ, et al., 2018,
Geometry, spatial arrangement and origin of carbonate grain-dominated, scour-fill and event-bed deposits: Late Jurassic Jubaila Formation and Arab-D Member, Saudi Arabia
, Sedimentology, Vol: 65, Pages: 1043-1066, ISSN: 0037-0746Outcrop analogues of the Late Jurassic lower Arab-D reservoir zone in Saudi Arabia expose a succession of fining-upward cycles deposited on a distal middle-ramp to outer-ramp setting. These cycles are interrupted by erosional scours that incise up to 1·8 m into underlying deposits and are infilled with intraclasts up to boulder size (1 m diameter). Scours of similar size and infill are not commonly observed on low-angle carbonate ramps. Outcrops have been used to characterize and quantify facies-body geometries and spatial relationships. The coarse grain size of scour-fills indicates scouring and boulder transport by debris or hyperconcentrated density flows strengthened by offshore-directed currents. Longitudinal and lateral flow transformation is invoked to produce the ‘pit and wing’ geometry of the scours. Scour pits and wings erode up to 1·8 m and 0·7 m deep, respectively, and are on average 50 m wide between wing tips. The flat bases of the scours and their lack of consistent aspect ratio indicate that erosion depth was limited by the presence of cemented firmgrounds in underlying cycles. Scours define slightly sinuous channels that are consistently oriented north–south, sub-parallel to the inferred regional depositional strike of the ramp, suggesting that local palaeobathymetry was more complex than commonly assumed. Weak lateral clustering of some scours indicates that they were underfilled and reoccupied by later scour incision and infill. Rudstone scour-fills required reworking of material from inner ramp by high-energy, offshore-directed flows, associated with storm action and the hydraulic gradient produced by coastal storm setup, to generate erosion and sustain transport of clasts that are generally associated with steeper slopes. Quantitative analysis indicates that these coarse-grained units have limited potential for correlation between wells as laterally continuous, highly permeable rese
- Abstract
- Open Access Link
- Cite
- Citations: 13
-
Journal articleDebbabi Y, Jackson M, Hampson G, et al., 2018,
Impact of the buoyancy–viscous force balance on two-phase flow in layered porous media.
, Transport in Porous Media, Vol: 2018, ISSN: 0169-3913Motivated by geological carbon storage and hydrocarbon recovery, the effect of buoyancy and viscous forces on the displacement of one fluid by a second immiscible fluid, along parallel and dipping layers of contrasting permeability, is characterized using five independent dimensionless numbers and a dimensionless storage or recovery efficiency. Application of simple dimensionless models shows that increased longitudinal buoyancy effects increase storage efficiency by reducing the distance between the leading edges of the injected phase in each layer and decreasing the residual displaced phase saturation behind the leading edge of the displacing phase. Increased transverse buoyancy crossflow increases storage efficiency if it competes with permeability layering effects, but reduces storage efficiency otherwise. When both longitudinal and transverse buoyancy effects are varied simultaneously, a purely geometrical dip angle group defines whether changes in storage efficiency are dominated by changes in the longitudinal or transverse buoyancy effects. In the limit of buoyancy-segregated flow, we report an equivalent, unidimensional flow model which allows rapid prediction of storage efficiency. The model presented accounts for both dip and layering, thereby generalizing earlier work which accounted for each of these but not both together. We suggest that the predicted storage efficiency can be used to compare and rank geostatistical realizations, and complements earlier heterogeneity measures which are applicable in the viscous limit.
-
Journal articleDebbabi Y, Stern D, Hampson GJ, et al., 2018,
Use of dimensionless scaling groups to interpret reservoir simulation results
, Journal of Petroleum Science and Engineering, Vol: 163, Pages: 270-282, ISSN: 1873-4715In conducting studies to make reservoir management decisions, it is important to efficiently interpret results of reservoir simulations. An understanding of how and why predicted reservoir performance changes with model parameters guides evaluation of production strategies as well as exploration of the impact of uncertainty in reservoir description. The aim of this work is to demonstrate the use of dimensionless scaling groups to interpret and qualitatively predict simulation results of multiphase flow in subsurface reservoirs with a large number of wells. Dimensionless scaling groups which quantify the balance between the forces causing fluid flow were computed between well pairs to rationalize simulation results. The data required to partition the model according to injector-producer pairs and estimate the scaling groups were obtained within minutes using simplified, single-phase numerical experiments. We show that scaling groups can be used to classify multiphase flow behaviours observed over the field into a small set of flow regimes characterized by the combination of their dominant forces. Changes in fluid distribution and reservoir performance with the model parameters can be analyzed in terms of changes in the force balance, and qualitatively predicted using the scaling groups. Predictions made using scaling groups may guide, and thereby reduce the use of, time-consuming multiphase flow simulations to optimize field development plans, to improve the calibration of reservoir models to production data and interpreted subsurface heterogeneity, and to assess the impact of reservoir uncertainties on production.
-
Journal articleMacAllister DJ, Jackson MD, Butler AP, et al., 2018,
Remote detection of saline intrusion in a coastal aquifer using borehole measurements of self potential
, Water Resources Research, Vol: 54, Pages: 1669-1687, ISSN: 0043-1397Two years of self‐potential (SP) measurements were made in a monitoring borehole in the coastal UK Chalk aquifer. The borehole SP data showed a persistent gradient with depth, and temporal variations with a tidal power spectrum consistent with ocean tides. No gradient with depth was observed at a second coastal monitoring borehole ca. 1 km further inland, and no gradient or tidal power spectrum were observed at an inland site ca. 80 km from the coast. Numerical modeling suggests that the SP gradient recorded in the coastal monitoring borehole is dominated by the exclusion‐diffusion potential, which arises from the concentration gradient across a saline front in close proximity to, but not intersecting, the base of the borehole. No such saline front is present at the two other monitoring sites. Modeling further suggests that the ocean tidal SP response in the borehole, measured prior to breakthrough of saline water, is dominated by the exclusion‐diffusion potential across the saline front, and that the SP fluctuations are due to the tidal movement of the remote front. The electrokinetic potential, caused by changes in hydraulic head across the tide, is one order of magnitude too small to explain the observed SP data. The results suggest that in coastal aquifers, the exclusion‐diffusion potential plays a dominant role in borehole SP when a saline front is nearby. The SP gradient with depth indicates the close proximity of the saline front to the borehole and changes in SP at the borehole reflect changes in the location of the saline front. Thus, SP monitoring can be used to facilitate more proactive management of abstraction and saline intrusion in coastal aquifers.
-
Journal articleZhang Z, Geiger S, Rood M, et al., 2017,
A Tracing Algorithm for Flow Diagnostics on Fully Unstructured Grids With Multipoint Flux Approximation
, SPE Journal, Vol: 22, Pages: 1946-1962, ISSN: 1930-0220Flow diagnostics is a common way to rank and cluster ensembles of reservoir models depending on their approximate dynamic behavior before beginning full-physics reservoir simulation. Traditionally, they have been performed on corner-point grids inherent to geocellular models. The rapid-reservoir-modeling (RRM) concept aims at fast and intuitive prototyping of geologically realistic reservoir models. In RRM, complex reservoir heterogeneities are modeled as discrete volumes bounded by surfaces that are sketched in real time. The resulting reservoir models are discretized by use of fully unstructured tetrahedral meshes where the grid conforms to the reservoir geometry, hence preserving the original geological structures that have been modeled.This paper presents a computationally efficient work flow for flow diagnostics on fully unstructured grids. The control-volume finite-element method (CVFEM) is used to solve the elliptic pressure equation. The flux field is a multipoint flux approximation (MPFA). A new tracing algorithm is developed on a reduced monotone acyclic graph for the hyperbolic transport equations of time of flight (TOF) and tracer distributions. An optimal reordering technique is used to deal with each control volume locally such that the hyperbolic equations can be computed in an efficient node-by-node manner. This reordering algorithm scales linearly with the number of unknowns.The results of these computations allow us to estimate swept-reservoir volumes, injector/producer pairs, well-allocation factors, flow capacity, storage capacity, and dynamic Lorenz coefficients, which all help approximate the dynamic reservoir behavior. The total central-processing-unit (CPU) time, including grid generation and flow diagnostics, is typically a few seconds for meshes with O (100,000) unknowns. Such fast calculations provide, for the first time, real-time feedback in the dynamic reservoir behavior while models are prototyped.
-
Journal articleSalinas P, Pavlidis D, Xie Z, et al., 2017,
A Discontinuous Control Volume Finite Element Method for Multi-Phase Flow in Heterogeneous Porous Media
, Journal of Computational Physics, Vol: 352, Pages: 602-614, ISSN: 0021-9991We present a new, high-order, control-volume-finite-element (CVFE) method for multiphase porous media flow with discontinuous 1st-order representation for pressure and discontinuous 2nd-order representation for velocity. The method has been implemented using unstructured tetrahedral meshes to discretize space. The method locally and globally conserves mass. However, unlike conventional CVFE formulations, the method presented here does not require the use of control volumes (CVs) that span the boundaries between domains with differing material properties. We demonstrate that the approach accurately preserves discontinuous saturation changes caused by permeability variations across such boundaries, allowing efficient simulation of flow in highly heterogeneous models. Moreover, accurate solutions are obtained at significantly lower computational cost than using conventional CVFE methods. We resolve a long-standing problem associated with the use of classical CVFE methods to model flow in highly heterogeneous porous media.
-
Journal articleDebbabi Y, Jackson MD, Hampson GJ, et al., 2017,
Capillary Heterogeneity Trapping and Crossflow in Layered Porous Media
, Transport in Porous Media, Vol: 120, Pages: 183-206, ISSN: 0169-3913We examine the effect of capillary and viscous forces on the displacement of one fluid by a second, immiscible fluid across and along parallel layers of contrasting porosity, and relative permeability, as well as previously explored contrasts in absolute permeability and capillary pressure. We consider displacements with wetting, intermediate-wetting and non-wetting injected phases. Flow is characterized using six independent dimensionless numbers and a dimensionless storage efficiency, which is numerically equivalent to the recovery efficiency. Results are directly applicable to geologic carbon storage and hydrocarbon production. We predict how the capillary–viscous force balance influences storage efficiency as a function of a small number of key dimensionless parameters, and provide a framework to support mechanistic interpretations of complex field or experimental data, and numerical model predictions, through the use of simple dimensionless models. When flow is directed across layers, we find that capillary heterogeneity traps the non-wetting phase, regardless of whether it is the injected or displaced phase. However, minimal trapping occurs when the injected phase is intermediate-wetting or when high-permeability layers contain a smaller moveable volume of fluid than low-permeability layers. A dimensionless capillary-to-viscous number defined using the layer thickness rather than the more commonly used system length is most relevant to predict capillary heterogeneity trapping. When flow is directed along layers, we show that, regardless of wettability, increasing capillary crossflow reduces the distance between the leading edges of the injected phase in each layer and increases storage efficiency. This may be counter-intuitive when the injected phase is non-wetting. Crossflow has a significant impact on storage efficiency only when high-permeability layers contain a smaller moveable volume of fluid than low-permeability layers. In that case, capillary he
-
Journal articleZhang J, Vinogradov J, Leinov E, et al., 2017,
Streaming potential during drainage and imbibition
, JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, Vol: 122, Pages: 4413-4435, ISSN: 2169-9313 -
Journal articleSalinas P, Pavlidis D, Xie Z, et al., 2017,
Improving the robustness of the control volume finite element method with application to multiphase porous media flow
, International Journal for Numerical Methods in Fluids, Vol: 85, Pages: 235-246, ISSN: 1097-0363Control volume finite element methods (CVFEMs) have been proposed to simulate flow in heterogeneous porous media because they are better able to capture complex geometries using unstructured meshes. However, producing good quality meshes in such models is nontrivial and may sometimes be impossible, especially when all or parts of the domains have very large aspect ratio. A novel CVFEM is proposed here that uses a control volume representation for pressure and yields significant improvements in the quality of the pressure matrix. The method is initially evaluated and then applied to a series of test cases using unstructured (triangular/tetrahedral) meshes, and numerical results are in good agreement with semianalytically obtained solutions. The convergence of the pressure matrix is then studied using complex, heterogeneous example problems. The results demonstrate that the new formulation yields a pressure matrix than can be solved efficiently even on highly distorted, tetrahedral meshes in models of heterogeneous porous media with large permeability contrasts. The new approach allows effective application of CVFEM in such models.
-
Journal articleDebbabi Y, Jackson MD, Hampson GJ, et al., 2017,
Viscous crossflow in layered porous media
, Transport in Porous Media, Vol: 117, Pages: 281-309, ISSN: 1573-1634We examine the effect of viscous forces on the displacement of one fluid by a second, immiscible fluid along parallel layers of contrasting porosity, absolute permeability and relative permeability. Flow is characterized using five dimensionless numbers and the dimensionless storage efficiency, so results are directly applicable, regardless of scale, to geologic carbon storage. The storage efficiency is numerically equivalent to the recovery efficiency, applicable to hydrocarbon production. We quantify the shock-front velocities at the leading edge of the displacing phase using asymptotic flow solutions obtained in the limits of no crossflow and equilibrium crossflow. The shock-front velocities can be used to identify a fast layer and a slow layer, although in some cases the shock-front velocities are identical even though the layers have contrasting properties. Three crossflow regimes are identified and defined with respect to the fast and slow shock-front mobility ratios, using both theoretical predictions and confirmation from numerical flow simulations. Previous studies have identified only two crossflow regimes. Contrasts in porosity and relative permeability exert a significant influence on contrasts in the shock-front velocities and on storage efficiency, in addition to previously examined contrasts in absolute permeability. Previous studies concluded that the maximum storage efficiency is obtained for unit permeability ratio; this is true only if there are no contrasts in porosity and relative permeability. The impact of crossflow on storage efficiency depends on the mobility ratio evaluated across the fast shock-front and on the time at which the efficiency is measured.
-
Journal articleMaes J, Muggeridge AH, Jackson MD, et al., 2017,
Scaling analysis of the In-Situ Upgrading of heavy oil and oil shale
, FUEL, Vol: 195, Pages: 299-313, ISSN: 0016-2361The In-Situ Upgrading (ISU) of heavy oil and oil shale is investigated. We develop a mathematical model for the process and identify the full set of dimensionless numbers describing the model. We demonstrate that for a model with nf fluid components (gas and oil), ns solid components and k chemical reactions, the model was represented by 9+k×(3+nf+ns-2)+8nf+2ns dimensionless numbers. We calculated a range of values for each dimensionless numbers from a literature study. Then, we perform a sensitivity analysis using Design of Experiments (DOE) and Response Surface Methodology (RSM) to identify the primary parameters controlling the production time and energy efficiency of the process. The Damköhler numbers, quantifying the ratio of chemical reaction rate to heat conduction rate for each reaction, are found to be the most important parameters of the study. They depend mostly on the activation energy of the reactions and of the heaters temperature. The reduced reaction enthalpies are also important parameters and should be evaluated accurately. We show that for the two test cases considered in this paper, the Damköhler numbers needed to be at least 10 for the process to be efficient. We demonstrate the existence of an optimal heater temperature for the process and obtain a correlation that can be used to estimate it using the minimum of the Damköhler numbers of all reactions.
-
Journal articleAl Mahrouqi D, Vinogradov J, Jackson MD, 2016,
Zeta potential of artificial and natural calcite in aqueous solution
, Advances in Colloid and Interface Science, Vol: 240, Pages: 60-76, ISSN: 0001-8686Despite the broad range of interest and applications, controls on calcite surface charge in aqueous solution, especiallyat conditions relevant to natural systems, remain poorly understood. The primary data source to understandcalcite surface charge comprises measurements of zeta potential. Here we collate and review previousmeasurements of zeta potential on natural and artificial calcite and carbonate as a resource for future studies,compare and contrast the results of these studies to determine key controls on zeta potential and where uncertaintiesremain, and report new measurements of zeta potential relevant to natural subsurface systems.The results show that the potential determining ions (PDIs) for the carbonate mineral surface are the lattice ionsCa2+, Mg2+ and CO32−. The zeta potential is controlled by the concentration-dependent adsorption of these ionswithin the Stern layer, primarily at the Outer Helmholtz Plane (OHP). Given this, the Iso-Electric Point (IEP) atwhich the zeta potential is zero should be expressed as pCa (or pMg). It should not be reported as pH, similarto most metal oxides.The pH does not directly control the zeta potential. Varying the pH whilst holding pCa constant yields constantzeta potential. The pH affects the zeta potential only by moderating the equilibrium pCa for a given CO2 partialpressure (pCO2). Experimental studies that appear to yield a systematic relationship between pH and zeta potentialare most likely observing the relationship between pCa and zeta potential, with pCa responding to the changein pH. New data presented here show a consistent linear relationship between equilibrium pH and equilibriumpCa or pMg irrespective of sample used or solution ionic strength. The surface charge of calcite is weakly dependenton pH, through protonation and deprotonation reactions that occur within a hydrolysis layer immediatelyadjacent to the mineral surface. The Point of Zero Charge (PZC) at which the surface charge is zero could b
-
Journal articleAl Mahrouqi D, Vinogradov J, Jackson MD, 2016,
Zeta potential of artificial and natural calcite in aqueous solution.
, Adv Colloid Interface Sci, Vol: 240, Pages: 60-76Despite the broad range of interest and applications, controls on calcite surface charge in aqueous solution, especially at conditions relevant to natural systems, remain poorly understood. The primary data source to understand calcite surface charge comprises measurements of zeta potential. Here we collate and review previous measurements of zeta potential on natural and artificial calcite and carbonate as a resource for future studies, compare and contrast the results of these studies to determine key controls on zeta potential and where uncertainties remain, and report new measurements of zeta potential relevant to natural subsurface systems. The results show that the potential determining ions (PDIs) for the carbonate mineral surface are the lattice ions Ca(2+), Mg(2+) and CO3(2-). The zeta potential is controlled by the concentration-dependent adsorption of these ions within the Stern layer, primarily at the Outer Helmholtz Plane (OHP). Given this, the Iso-Electric Point (IEP) at which the zeta potential is zero should be expressed as pCa (or pMg). It should not be reported as pH, similar to most metal oxides. The pH does not directly control the zeta potential. Varying the pH whilst holding pCa constant yields constant zeta potential. The pH affects the zeta potential only by moderating the equilibrium pCa for a given CO2 partial pressure (pCO2). Experimental studies that appear to yield a systematic relationship between pH and zeta potential are most likely observing the relationship between pCa and zeta potential, with pCa responding to the change in pH. New data presented here show a consistent linear relationship between equilibrium pH and equilibrium pCa or pMg irrespective of sample used or solution ionic strength. The surface charge of calcite is weakly dependent on pH, through protonation and deprotonation reactions that occur within a hydrolysis layer immediately adjacent to the mineral surface. The Point of Zero Charge (PZC) at which the surface char
-
Journal articleSalinas P, Pavlidis D, xie Z, et al., 2016,
Improving the convergence behaviour of a fixed-point-iteration solver for multiphase flow in porous media
, International Journal for Numerical Methods in Fluids, Vol: 84, Pages: 466-476, ISSN: 1097-0363A new method to admit large Courant numbers in the numerical simulation of multiphase flow is presented.The governing equations are discretised in time using an adaptive -method. However, the use of implicitdiscretisations does not guarantee convergence of the non-linear solver for large Courant numbers. In thiswork, a double-fixed point iteration method with backtracking is presented that improves both convergenceand convergence rate. Moreover, acceleration techniques are presented to yield a more robust non-linearsolver with increased effective convergence rate. The new method reduces the computational effort bystrengthening the coupling between saturation and velocity, obtaining an efficient backtracking parameter,using a modified version of Anderson’s acceleration and adding vanishing artificial diffusion.
-
Journal articleMacAllister DJ, Jackson MD, Butler AP, et al., 2016,
Tidal influence on self-potential measurements
, Journal of Geophysical Research. Solid Earth, Vol: 121, Pages: 8432-8452, ISSN: 2169-9313Long-term surface and borehole self-potential (SP) monitoring was conducted in the UK Chalk aquifer at two sites. The coastal site is ~1.7 km from the coast, and the inland site is ~80 km from the coast. At both sites, power spectral density analysis revealed that SP data contain the main ocean tidal periodic components. However, the principal lunar component (M2), the dominant ocean tidal component, was most significant at the coastal site. The M2 signal in surface-referenced SP data at the inland site was partly due to telluric currents caused by the geomagnetic ocean dynamo. Earth and/or atmospheric tides also contributed, as the SP power spectrum was not typical of a telluric electric field. The M2 component in borehole-referenced data at the inland site was below the significance level of the analysis method and was 2 orders of magnitude smaller than the M2 signal in borehole-referenced SP data at the coastal site. The tidal response of the SP data in the coastal borehole is, therefore, primarily driven by ocean tides. These cause changes in fluid pressure and chemical concentration gradients within the coastal aquifer, leading to time varying electrokinetic and exclusion-diffusion potentials. Borehole-referenced SP measurements could be used to characterize and monitor tidal processes in coastal aquifers such as the intrusion of seawater.
-
Journal articleMostaghimi P, Kamali F, Jackson MD, et al., 2016,
Adaptive mesh optimization for simulation of immiscible viscous fingering
, SPE Journal, Vol: 21, Pages: 2250-2259, ISSN: 1086-055XViscous fingering can be a major concern when waterflooding heavy-oil reservoirs. Most commercial reservoir simulators use low-order finite-volume/-difference methods on structured grids to resolve this phenomenon. However, this approach suffers from a significant numerical-dispersion error because of insufficient mesh resolution, which smears out some important features of the flow. We simulate immiscible incompressible two-phase displacements and propose the use of unstructured control-volume finite-element (CVFE) methods for capturing viscous fingering in porous media. Our approach uses anisotropic mesh adaptation where the mesh resolution is optimized on the basis of the evolving features of flow. The adaptive algorithm uses a metric tensor field dependent on solution-interpolation-error estimates to locally control the size and shape of elements in the metric. The mesh optimization generates an unstructured finer mesh in areas of the domain where flow properties change more quickly and a coarser mesh in other regions where properties do not vary so rapidly. We analyze the computational cost of mesh adaptivity on unstructured mesh and compare its results with those obtained by a commercial reservoir simulator on the basis of the finite-volume methods.
-
Journal articleJackson MD, Al-Mahrouqi D, Vinogradov J, 2016,
Zeta potential in oil-water-carbonate systems and its impact on oil recovery during controlled salinity water-flooding
, Scientific Reports, Vol: 6, ISSN: 2045-2322Laboratory experiments and field trials have shown that oil recovery from carbonate reservoirs can be increased by modifying the brine composition injected during recovery in a process termed controlled salinity water-flooding (CSW). However, CSW remains poorly understood and there is no method to predict the optimum CSW composition. This work demonstrates for the first time that improved oil recovery (IOR) during CSW is strongly correlated to changes in zeta potential at both the mineral-water and oil-water interfaces. We report experiments in which IOR during CSW occurs only when the change in brine composition induces a repulsive electrostatic force between the oil-brine and mineral-brine interfaces. The polarity of the zeta potential at both interfaces must be determined when designing the optimum CSW composition. A new experimental method is presented that allows this. Results also show for the first time that the zeta potential at the oil-water interface may be positive at conditions relevant to carbonate reservoirs. A key challenge for any model of CSW is to explain why IOR is not always observed. Here we suggest that failures using the conventional (dilution) approach to CSW may have been caused by a positively charged oil-water interface that had not been identified.
-
Journal articleAl-Mahrouqi, Vinogradov J, Jackson MD, 2016,
Temperature-dependence of the zeta potential in intact natural carbonates
, Geophysical Research Letters, Vol: 43, Pages: 11578-11587, ISSN: 1944-8007The zeta potential is a measure of the electrical charge on mineral surfaces and is an important control on subsurface geophysical monitoring, adsorption of polar species in aquifers, and rock wettability. We report the first measurements of zeta potential in intact, water-saturated, natural carbonate samples at temperatures up to 120°C. The zeta potential is negative and decreases in magnitude with increasing temperature at low ionic strength (0.01 M NaCl, comparable to potable water) but is independent of temperature at high ionic strength (0.5 M NaCl, comparable to seawater). The equilibrium calcium concentration resulting from carbonate dissolution also increases with increasing temperature at low ionic strength but is independent of temperature at high ionic strength. The temperature dependence of the zeta potential is correlated with the temperature dependence of the equilibrium calcium concentration and shows a Nernstian linear relationship. Our findings are applicable to many subsurface carbonate rocks at elevated temperature.
-
Journal articleAl Mahrouqi D, Vinogradov J, Jackson MD, 2016,
Temperature dependence of the zeta potential in intact natural carbonates
, GEOPHYSICAL RESEARCH LETTERS, Vol: 43, Pages: 11578-11587, ISSN: 0094-8276 -
Conference paperSalinas P, Pavlidis D, Xie Z, et al., 2016,
Dynamic unstructured mesh adaptivity for improved simulation of nearwellbore flow in reservoir scale models
, 15th European Conference on the Mathematics of Oil Recovery, Publisher: EAGEIt is well known that the pressure gradient into a production well increases with decreasing distanceto the well and may cause downwards coning of the gaswater interface, or upwards coning ofwateroil interface, into oil production wells; it can also cause downwards coning of the water table,or upwards coning of a saline interface, into water abstraction wells. To properly capture the localpressure drawdown into the well, and its effect on coning, requires high grid or mesh resolution innumerical models; moreover, the location of the well must be captured accurately. In conventionalsimulation models, the user must interact with the model to modify grid resolution around wells ofinterest, and the well location is approximated on a grid defined early in the modelling process.We report a new approach for improved simulation of nearwellbore flow in reservoirscale modelsthrough the use of dynamic unstructured adaptive meshing. The method is novel for two reasons.First, a fully unstructured tetrahedral mesh is used to discretize space, and the spatial location of thewell is specified via a line vector. Mesh nodes are placed along the line vector, so the geometry ofthe mesh conforms to the well trajectory. The well location is therefore accurately captured, and theapproach allows complex well trajectories and wells with many laterals to be modelled. Second,the mesh automatically adapts during a simulation to key solution fields of interest such as pressureand/or saturation, placing higher resolution where required to reduce an error metric based on theHessian of the field. This allows the local pressure drawdown and associated coning to be capturedwithout userdriven modification of the mesh. We demonstrate that the method has wideapplication in reservoirscale models of oil and gas fields, and regional models of groundwaterresources.
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.