Micropores can constitute up to 100% of the total porosity within carbonate hosted hydrocarbon reservoirs, usually existing within micritic fabrics. There is, however, only a rudimentary understanding of the contribution that these pores make to reservoir performance and hydrocarbon recovery. To further our understanding, a flexible, object-based algorithm has been developed to produce 3D computational representations of end-point micritic fabrics. By methodically altering model parameters, the state-space of microporous carbonates is explored. Flow properties are quantified using lattice-Boltzmann and network modelling methods. In purely micritic fabrics, it has been observed that average pore radius has a positive correlation with single-phase permeability and results in decreasing residual oil saturations under both water-wet and 50% fractionally oil-wet states. Similarly, permeability increases by an order of magnitude (from 0.6md to 7.5md) within fabrics of varying total matrix porosity (from 18% to 35%) due to increasing pore size, but this has minimal effect on multi-phase flow. Increased pore size due to micrite rounding notably increases permeability in comparison to original rhombic fabrics with the same porosity, but again, multi-phase flow properties are unaffected. The wetting state of these fabrics, however, can strongly influence multi-phase flow; residual oil saturations vary from 30% for a water-wet state and up to 50% for an 80% oil wet fraction. flow when directly connected. Otherwise, micropores control single-phase permeability magnitude. Importantly in these fabrics, recovery is dependent on both wetting scenario and pore-network homogeneity; under water-wet imbibition, increasing proportions of microporosity yield lower residual oil saturations. Finally, in grain-based fabrics where mesopores form an independently connected pore network, micropores do not affect permeability, even when they constitute up to 50% of the total porosity. Through examination of these three styles of microporous carbonates, it is apparent that micropores can have a significant impact on flow and sweep characteristics in such fabrics.

Horizontal wells are drilled to achieve improved reservoir coverage, high production rates, and to overcome water coning problems, etc. Many of these wells often produce at rates much below the expected production rates. Low productivity of horizontal wells is attributed to various factors such as drilling induced formation damage, high completion skins, and variable formation properties along the length of the wellbore as in the case of heterogeneous carbonate reservoirs. Matrix acidizing is used to overcome the formation damage by injecting the acid into the carbonate rock to improve well performance. Designing the matrix acidizing treatments for horizontal wells is a challenging task because of the complex process. The estimation of acid distribution along wellbore is required to analyze that the zones needing stimulation are receiving enough acid. It is even more important in cases where the reservoir properties are varying along the length of the wellbore. A model is developed in this study to simulate the placement of injected acid in a long horizontal well and to predict the subsequent effect of the acid in creating wormholes, overcoming damage effects, and stimulating productivity. The model tracks the interface between the acid and the completion fluid in the wellbore, models transient flow in the reservoir during acid injection, considers frictional effects in the tubulars, and predicts the depth of penetration of acid as a function of the acid volume and injection rate at all locations along the completion. A computer program is developed implementing the developed model. The program is used to simulate hypothetical examples of acid placement in a long horizontal section. A real field example of using the model to history match actual treatment data from a North Sea chalk well is demonstrated. The model will help to optimize acid stimulation in horizontal wells.

Horizontal wells are drilled to achieve improved reservoir coverage, high production rates, and to overcome water coning problems, etc. Many of these wells often produce at rates much below the expected production rates. Low productivity of horizontal wells is attributed to various factors such as drilling induced formation damage, high completion skins, and variable formation properties along the length of the wellbore as in the case of heterogeneous carbonate reservoirs. Matrix acidizing is used to overcome the formation damage by injecting the acid into the carbonate rock to improve well performance. Designing the matrix acidizing treatments for horizontal wells is a challenging task because of the complex process. The estimation of acid distribution along wellbore is required to analyze that the zones needing stimulation are receiving enough acid. It is even more important in cases where the reservoir properties are varying along the length of the wellbore. A model is developed in this study to simulate the placement of injected acid in a long horizontal well and to predict the subsequent effect of the acid in creating wormholes, overcoming damage effects, and stimulating productivity. The model tracks the interface between the acid and the completion fluid in the wellbore, models transient flow in the reservoir during acid injection, considers frictional effects in the tubulars, and predicts the depth of penetration of acid as a function of the acid volume and injection rate at all locations along the completion. A computer program is developed implementing the developed model. The program is used to simulate hypothetical examples of acid placement in a long horizontal section. A real field example of using the model to history match actual treatment data from a North Sea chalk well is demonstrated. The model will help to optimize acid stimulation in horizontal wells.

In this study, the optimum conditions for wormhole propagation in horizontal well carbonate acidizing was investigated numerically using a horizontal well acidizing simulator. The factors that affect the optimum conditions are rock mineralogy, acid concentration, temperature and acid flux in the formation. The work concentrated on the investigation of the acid flux. Analytical equations for injection rate schedule for different wormhole models. In carbonate acidizing, the existence of the optimum injection rate for wormhole propagation has been confirmed by many researchers for highly reactive acid/rock systems in linear core-flood experiments. There is, however, no reliable technique to translate the laboratory results to the field applications. It has also been observed that for radial flow regime in field acidizing treatments, there is no single value of acid injection rate for the optimum wormhole propagation. In addition, the optimum conditions are more difficult to achieve in matrix acidizing long horizontal wells. Therefore, the most efficient acid stimulation is only achieved with continuously increasing acid injection rates to always maintain the wormhole generation at the tip of the wormhole at its optimum conditions. Examples of acid treatments with the increasing rate schedules were compared to those of the single optimum injection rate and the maximum allowable rate. The comparison study showed that the increasing rate treatments had the longest wormhole penetration and, therefore, the least negative skin factor for the same amount of acid injected into the formations. A parametric study was conducted for the parameters that have the most significant effects on the wormhole propagation conditions such as injected acid volume, horizontal well length, acid concentration, and reservoir heterogeneity. The results showed that the optimum injection rate per unit length increases with increasing injected acid volume. And it was constant for scenarios with different lateral lengths for a given system of rock/ acid and injected volume. The study also indicated that for higher acid concentration the optimum injection rate was lower. It does exist for heterogeneous permeability formations. Field treatment data for horizontal wells in Middle East carbonate reservoirs were also analyzed for the validation of the numerical acidizing simulator.

Estimation of reservoir fracture parameters, fracture orientation and density, from seismic data is often difficult because of one important question: Is observed anisotropy caused by the reservoir interval or by the effect of the lithologic unit or multiple units above the reservoir? Often hydrocarbon reservoirs represent a small portion of the seismic section, and reservoir anisotropic parameter inversion can be easily obscured by the presence of an anisotropic overburden. In this study, I show examples where we can clearly observe imprints of overburden anisotropic layers on the seismic response of the target zone. Then I present a simple method to remove the effect of anisotropic overburden to recover reservoir fracture parameters. It involves analyzing amplitude variation with offset and azimuth (AVOA) for the top of reservoir reflector and for a reflector below the reservoir. Seismic CMP gathers are transformed to delay-time vs. slowness (tau-p) domain. We then calculate the ratio of the amplitudes of reflections at the reservoir top and from the reflector beneath the reservoir. The ratios of these amplitudes are then used to isolate the effect of the reservoir interval and remove the transmission effect of the overburden. The methodology is tested on two sets of models - one containing a fractured reservoir with isotropic overburden and the other containing a fractured reservoir with anisotropic overburden. Conventional analysis in the x-t domain indicates that the anisotropic overburden has completely obscured the anisotropic signature of the reservoir zone. When the new methodology is applied, the overburden effect is significantly reduced. The methodology is also applied to an actual PP surface reflection (Rpp) 3D dataset over a reservoir in the Arabian Peninsula. Ellipse-fitting technique was applied to invert for two Fracture parameters: (1) Fracture density and (2) fracture direction. Fracture density inversion results indicate increased fracturing in the anticline structure hinge zone. Fracture orientation inversion results agree with Formation MicroImaging (FMI) borehole logs showing a WNW-ESE trend. This newly developed amplitude ratio method is suitable for quantitative estimation of fracture parameters including normal and tangential “weaknesses” (ΔN and ΔT respectively). Initially, inversion of conventional AVOA for ΔN and ΔT parameters indicates that the ΔN parameter is reliably estimated given an accurate background isotropic parameter estimation derived from borehole logging data. While ΔN parameter inversion is successful, inversion for ΔT parameter from Rpp information is not, presumably due to the dependence of ΔT estimation on many medium parameters for accurate prediction. The ΔN parameter is then successfully recovered when applied to the amplitude ratio values derived from synthetic data. It is important to recognize that ΔN parameter is directly proportional to fracture density and high ΔN values can be attributed to high crack density values. The ΔN parameter inversion is also applied to the amplitude ratios derived from real seismic data. This inversion requires fracture azimuth data input that is obtained from the fracture direction inversion using ellipse-fitting technique. The background Vp/Vs ratio. / text

Pore type variations account for complex velocity-porosity relationship and intensive permeability heterogeneity and consequently low oil and gas recovery in carbonate reservoir. However, it is a challenge for geologist and geophysicist to quantitatively estimate the influences of pore type complexity on velocity variation at a given porosity and porosity-permeability relationship. A new rock physics-based integrated approach in this study was proposed to quantitatively characterize the diversity of pore types and its influences on wave propagation in carbonate reservoir. Based on above knowledge, permeability prediction accuracy from petrophysical data can be improved compared to conventional approach. Two carbonate reservoirs with different reservoir features, one is a shallow carbonate reservoir with average high porosity (>10%) and another one is a supper-deep carbonate reservoir with average low porosity (<5%), are used to test the proposed approach. Paleokarst is a major event to complicate carbonate reservoir pore structure. Because of limited data and lack of appropriate study methods, it is a difficulty to characterize subsurface paleokarst 3D distribution and estimate its influences on reservoir heterogeneity. A method by integrated seismic characterization is applied to delineate a complex subsurface paleokarst system in the Upper San Andres Formation, Permian basin, West Texas. Meanwhile, the complex paleokarst system is explained by using a carbonate platform hydrological model, similar to modern marine hydrological environments within carbonate islands. How to evaluate carbonate reservoir permeability heterogeneity from 3D seismic data has been a dream for reservoir geoscientists, which is a key factor to optimize reservoir development strategy and enhance reservoir recovery. A two-step seismic inversions approach by integrating angle-stack seismic data and rock physics model is proposed to characterize pore-types complexity and further to identify the relative high permeability gas-bearing zones in low porosity reservoir (< 5%) using ChangXing super-deep carbonate reservoir as an example. Compared to the conventional permeability calculation method by best-fit function between porosity and permeability, the results in this study demonstrate that gas zones and non-gas zones in low porosity reservoir can be differentiated by using above integrated permeability characterization method.

Characterization of carbonate reservoirs by flow units is a practical way of reservoir zonation. This study represents a petrophysical-based method that uses well loggings and core plug data to delineate flow units within the most productive carbonate reservoir of Derdere Formation in Y field, Southeast Turkey. Derdere Formation is composed of limestones and dolomites. Logs from the 5 wells are the starting point for the reservoir characterization. The general geologic framework obtained from the logs point out for discriminations within the formation. 58 representative core plug data from 4 different wells are utilized to better understand the petrophysical framework of the formation. The plots correlating petrophysical parameters and the frequency histograms suggest the presence of distinctive reservoir trends. These discriminations are also represented in Winland porosity-permeability crossplots resulted in clusters for different port-sizes that are responsible for different flow characteristics. Although the correlation between core plug porosity and air permeability yields a good correlation coefficient, the formation has to be studied within units due to differences in port-sizes and reservoir process speed. Linear regression and multiple regression analyses are used for the study of each unit. The results are performed using STATGRAPH Version Plus 5.1 statistical software. The permeability models are constructed and their reliabilities are compared by the regression coefficients for predictions in un-cored sections. As a result of this study, 4 different units are determined in the Derdere Formation by using well logging data, and core plug analyses with the help of geostatistical methods. The predicted permeabilities for each unit show good correlations with the calculated ones from core plugs. Highly reliable future estimations can be based on the derived methods.

In the history of a reservoir, after the period of primary production, about 30 to 40%, of the original oil in place may be produced using a secondary recovery mechanism. Polymer injection, which is classified as a tertiary method, can be applied to the remaining oil in place. In this thesis, oil/water relative permeabilities, effect of polymer injection on end point relative permeabilities and residual oil saturations in heterogeneous carbonate reservoirs were investigated. Numereous core flood experiments were conducted on different heteroegneous carbonate cores taken from Midyat Formation. Before starting the displacement experiments, porosity, permeability and capillary pressure experiments were performed. The heterogeneity of the cores are depicted from thin sections. Besides the main aim stated above, effect of flow rate and fracture presence on end point relative permeability and on residual oil saturation and were investigated. According to the results of the displacement tests, end point hexane relative permeability increased when polymer solution was used as the displacing phase.Besides, end point hexane relative permeability increased with polymer injection and fracture presence.

Master of Science / Department of Geology / Matthew Totten / The Petersilie oil field in Ness County, Kansas produces out of the Mississippian System, a reservoir composed mainly of shallow water carbonates, at depths of around 4375 ft (1334 m). The lithology of the field ranges from limestone to dolomite, to interlaminated limestone-dolomite beds. Chert is commonly found throughout. Petersilie field lies to the west of the Central Kansas Uplift, and to the east of the Hugoton Embayment. The field saw much drilling activity in the 1960’s, when it reached a production peak of nearly 378,000 barrels of oil per year. Production declined swiftly after that until the late 1990’s, when waterflooding was successfully employed. In this study, a reservoir model was produced for the Mississippian as it occurs in Petersilie field using the Department of Energy’s EdBOAST reservoir modeling software, with the intent of providing a reference for future drilling activity in the Mississippian and determining reservoir characteristics that may have contributed to the effectiveness of waterflooding in this area. The reservoir model was checked by simulation with a companion reservoir simulator program, BOAST 98. Subsequent comparison of simulated and actual oil production curves demonstrates the reliability of well log and drill stem test data for the field and proves the reservoir model to be a good fit for the Mississippian in Petersilie. Production curve analysis of Petersilie indicates the field was an ideal candidate for waterflooding because it has a solution-gas drive mechanism. As the field approached depletion from primary recovery, oil saturations remained high. Petersilie also exhibits high porosity and good permeability. The BOAST software was found to be an effective and inexpensive means for understanding the Mississippian reservoir in central to south-central Kansas. It was determined that BOAST has potential for practical use by smaller independent oil companies targeting the Mississippian in Kansas.

The transport of components and phases plays a fundamental role in the success of an EOR process. Because many reservoirs have harsh conditions of salinity, temperature and rock heterogeneity, which limit process options, a robust system with flexibility is required. Systematic experimental study of formulations capable to transport surfactant as foam at 94°C, formulated in sea water, is presented. It includes methodology to conduct core floods in sand packs using foaming surfactants and to develop “surfactant blend ratio- salinity ratio maps” using equilibrium phase behavior to determine favorable conditions for oil recovery in such floods. Mathematical model able to reproduce the foam strength behavior observed in sand packs with the formulations studied is presented. Visualization of oil recovery mechanism from matrix is realized using a model system of micro-channels surrounded by glass beads to mimic matrix and fractures respectively. The observations illustrate how components may distribute within the matrix, thereby releasing oil into the fractures. The use of chemicals to minimize adsorption is required when surfactant adsorption is important. The presence of anhydrite may limit the use of sodium carbonate to reduce adsorption of carbonates. A methodology is presented to estimate the amount, if any, of anhydrite present in the reservoir. The method is based on brine software analysis of produced water compositions and inductively coupled plasma (ICP) analysis of core samples. X-ray powder diffraction (XRD) was used to verify the mineralogy of the rock. X-ray photoelectron spectroscopy (XPS) was used to obtain surface composition for comparison with bulk composition of the rock. Adsorption of surfactants was measured using dynamic and static adsorption experiments. Determining the flow properties of the rock samples via tracer analysis permitted the simulation of the dynamic adsorption process using a mathematical model that considers the distribution of adsorbed materials in the three different regions of pore space. Using this method allows one to predict adsorption in a reservoir via simulation.