Nghiem, Long X. (Computer Modelling Group Ltd.) | Mirzabozorg, Arash (University of Calgary) | Chen, Zhangxin John (University of Calgary) | Hajizadeh, Yasin (Computer Modelling Group Ltd.) | Yang, Chaodong (Computer Modelling Group Inc.)
History matching of reservoir flow models based only on production data may not reveal deficiencies that affect future predictions. Incorporating saturation and temperature profile data that come from 4D seismic surveys in the history matching process can reduce the uncertainty of reservoir models for the prediction stage. We constructed a field reservoir model from which production history, saturation and temperature profile history were obtained. We started the history matching process with a base reservoir model, the petro-physical properties of which were substantially different than those of the field reservoir model. We propose a new methodology for matching the fluid and temperature profiles by adjusting reservoir petro-physical properties. In this methodology, some grid blocks in a reservoir model were selected judiciously to capture the overall saturation and temperature distribution profiles. In addition to well production data, we included the saturation and temperature profiles at these grid blocks as extra objective functions during the history matching process. The DECE optimization is used to reduce the objective function. We applied this method in a Steam Assisted Gravity Drainage (SAGD) process and matched the saturation and temperature profiles with an average error of less than 2%.
Excessive water production from unwanted zones in oil producing wells is one of the major challenges faced by the oil industry. The applicability of organically crosslinked polymer (OCP) systems as sealants for water shutoff treatments in temperatures up to 350°F is well documented. However, their effectiveness at temperatures above 350°F has not been evaluated. This paper presents experimental data from using an OCP system for water shutoff treatments at 400°F.
At temperatures around 400°F, crosslinking is expected to happen faster and can lead to premature gelation of the recipe before the entire treatment is in place. Thus, controlling the gelation time at such temperatures is extremely crucial. Optimizing the amount of retarder is essential to provide adequate time for placement of the treatment fluid. This paper provides gelation time data at temperatures between 350 and 400°F with different amounts of retarder. With an optimum amount of retarder, the OCP showed a gelation time of 1 hr 20 min.
This paper also describes the experimental setup used to study and determine the long-term stability of the OCP system at 400°F. Sand packs measuring 1-ft long were used for the test to simulate formation conditions. Once the optimized OCP recipe was gelled inside the sand pack, measurements were taken by gradually applying incremental differential pressure (?P) to evaluate the sealant at temperature, as well as the threshold ?P the system could withstand. Even after one month at 400°F, the OCP recipe was able to sustain a ?P of 950 psi over the sand pack.
The data indicates the applicability of this system as an effective conformance product to shut off water-producing zones over an extended period of time at 400°F.
Sanyal, Tirtharenu (Kuwait Oil Company) | Al-Hamad, Khairyah (KOC) | Jain, Anil Kumar (KOC) | Al-Haddad, Ali Abbas (KISR) | Kholosy, Sohib (KISR) | Ali, Mohammad A.J. (Kuwait Inst. Scientific Rsch.) | Abu Sennah, Heba Farag (Kuwait Oil Company)
Improved oil recovery for heavy oil reservoirs is becoming a new research study for Kuwaiti reservoirs. There are two mechanisms for improved oil recovery by thermal methods. The first method is to heat the oil to higher temperatures, and thereby, decrease its viscosity for improved mobility. The second mechanism is similar to water flooding, in which oil is displaced to the production wells. While more steam is needed for this method than for the cyclic method, it is typically more effective at recovering a larger portion of the oil.
Steam injection heats up the oil and reduce its viscosity for better mobility and higher sweep efficiency. During this process, the velocity of the moving oil increases with lower viscosity oil; and thus, the heated zone around the injection well will have high velocity. The increase of velocity in an unconsolidated formation is usually accompanied with sand movement in the reservoir creating a potential problem.
The objective of this study was to understand the effect of flowrate and viscosity on sand production in heavy oil reservoir that is subjected for thermal recovery process. The results would be useful for designing completion under steam injection where the viscosity of the oil is expected to change due to thermal operations.
A total of 21 representative core samples were selected from different wells in Kuwait. A reservoir condition core flooding system was used to flow oil into the core plugs and to examine sand production. Initially, the baseline liquid permeability was measured with low viscosity oil and low flowrate. Then, the flowrate was increased gradually and monitored to establish the value for sand movement for each plug sample. At the end of the test, the produced oil containing sand was filtered for sand content.
The result showed that sand production increased with higher viscosity oil and high flowrate. However, sand compaction at the injection face of the cores was more significant than sand production. In addition, high confining pressure contributes to additional sand production. The average critical velocity was estimated ranged from 18 to 257 ft/day for the 0.74 cp oil, 2 to 121 ft/day for the 16 cp oil, and 1 to 26 ft/day for the 684 cp oil.
Transverse fractures created from horizontal wells are a common choice in tight and shale gas reservoirs. Previous work has shown that proppant pack permeability reduction due to non-Darcy flow in a transverse fracture from a horizontal well causes significant reduction in the fracture performance when the gas formation permeability exceeds 0.5 md. There are other configurations and architectures such as aligning the well trajectory with the fracture, either by drilling horizontal wells in the direction that results in longitudinal fractures or by just sticking with drilling vertical wells. However, when drilling and fracturing costs are considered, productivity is not the only optimization consideration.
The field example illustrates a case when the apparent choice to use transverse fractures from horizontal wells proved to be suboptimal from the productivity perspective, but fundamental considering economics. Parametric studies for permeability ranging from 0.01 to 5 md illustrate the importance of economics in addition to physical performance. For similar reservoir characteristics, the optimum fractured well architecture varies considerably, and therefore an extensive reservoir engineering approach may be necessary beyond the well completions and/or current prejudices and inadequate understanding.
A live oil sample was subjected to a solid detection system (SDS) to measure asphaltene onset point (AOP) at 3850 psi, and asphaltene content of 1.3%. A high-resolution digital camera was used to measure asphaltene particle size distribution. The result showed that asphaltene particles were not uniform in size, but has a normal distribution of 100-120 µm. Asphaltene reversibility to dissolved back into the oil with increasing pressure was only 35% of the original deposition. Two core samples were examined for formation damage due to asphaltene deposition. A Low permeability core showed significant permeability reduction exceeding 50% of its baseline permeability, and the higher permeability core showed less permeability decline, even with the same asphaltene precipitation.
The field X is a brown heavy oil field producing under strong bottom water drive since the mid-1980. Production is from a combination of Amin aeolian and Al Khlata glacial reservoir sediments. At present, the development is focused on drilling horizontal infill wells. One of the biggest challenges is the unfavorable mobility contrast between the heavy oil and water causing early water breakthrough.
The Amin Formation, the primary reservoir, is characterized by a high net to gross ratio and an average porosity of 30 %. However the initial hydrocarbon saturation at the same porosity often varies by 20 % in different parts of the field. Furthermore, core measurements show an order of magnitude scatter in permeability at the same porosity, indicating the presence of different facies. In early studies these variations were attributed mainly to the grain size variations. A later petrographical study found that the abundance of clays and feldspars could also severely reduce permeability, but may retain high porosity.
In the current Study it was found that the rocks have variable radioactivity due to the presence of radioactive Potassium isotope associated with feldspars. A fare correlation was observed between the grain size and the content of feldspars from core. A novel approach to reservoir characterization integrating core and logs was developed leading to a major breakthrough in the reservoir characterization including:
• Enhanced permeability prediction using normalized Gamma Ray (GR) log as 3rd parameter;
• Facies identification using normalized Gamma Ray cut-off;
• Facies based Saturation-Height models.
This work is a good example of advances in reservoir characterization achieved by integrating core and log data. It results in better understanding of reservoir properties distribution, optimization of completions of new wells and improvement of further development scenarios. In particular, abnormally high gross production and high water cut in the north of the field is currently in line with new facies scheme.
The purpose of history matching is to achieve geological realizations calibrated to the historical performance of the reservoir. For complex geological structures it is usually intractable to run tens of thousands of full reservoir simulation to trace the most probable geological model. Hence the inadequacy of the history-matching results frequently leads to poor estimation of the true model and high uncertainty in production forecasting. Reduced-order modeling procedures, which have been applied in many application areas including reservoir simulation, represent a promising means for constructing efficient surrogate models. Nonlinear dimensionality reduction techniques allow for encapsulating the high-resolution complex geological description of reservoir into a low-dimensional subspace, which significantly reduces number of unknowns and provides an efficient way to construct a proxy model based on the the reduced-dimension parameters.
Polynomial Chaos Expansions (PCE) is a powerful tool to quantify uncertainty in dynamical system when there is probabilistic uncertainty in the system parameters. In reservoir simulation it has been shown to be more accurate and efficient compared to traditional experimental design (ED). PCEs have a significant advantage over other response surfaces as the convergence to the true probability distribution is proved when the order of the PCE is increased. Accordingly PCE proxy can be used as the pseudo-simulator to represent the surface responses of the uncertain variables. When the objective and constraints of a reservoir model is described by multivariate polynomial functions, there are very efficient algorithms to compute the global solutions. We have developed a workflow at which incorporates PCE to find the global minimum of the misfit surface and assess the uncertainty associated with. The accuracy of the PCE proxy increases with the additional trial runs of the reservoir simulator.
We conduct a two dimensional synthetic case study of a fluvial channel as well as a real field example to demonstrate the effectiveness of this approach. Kernel Principal Component Analysis (KPCA) is used to parameterize the complex geological structure. The study has revealed useful reservoir information and delivered more reliable production forecasts.
PCE-based history match enhances the quality and efficiency of the estimation of the most probable geological model and improve the confidence interval of production forecasts.
Cinar, Yildiray (The University of New South Wales) | Arns, Christoph (The University of New South Wales) | Dehghan Khalili, Ahmad (The University of New South Wales) | Yanici, Sefer (The University of New South Wales)
Resistivity measurements play a key role in hydrocarbon in place calculations for oil and gas reservoirs. They are a direct indi-cator of fluid saturation and connected pore space available in the formation. Carbonate rocks, which host around half of the world's hydrocarbons, exhibit a wide range of porosities with scales spanning from nanometres to centimetres. The often sig-nificant amount of microporosity displayed by Carbonate rocks emphasizes the necessity of an adequate characterization of their micro-features and their contribution to hydrocarbon in place. In this paper we examine upscaling methods to probe for-mation factor of a fully saturated carbonate sample using an X-ray CT based numerical approach and compare to experimental measurements.
Three-dimensional high-resolution X-ray CT enables the numerical calculation of petrophysical properties of interest at the pore scale with resolutions down to a few microns per voxel. For more complex and heterogeneous samples however, a direct calculation of petrophysical properties is not feasible, since the required resolution and a sufficient field of view cannot be obtained simultaneously. Thus an integration of measurements at different scale is required. In this study a carbonate sample of 38mm in diameter is first scanned using the X-ray CT method with a resolution of 26µm. After accompanying experimental measurements on the full plug, four 5mm plugs were drilled vertically from this sample and X-ray CT images of these plugs acquired at resolutions down to 2.74 µm. We calculate the porosity of the sample (macro- and micro-porosities) using the phase separation methods and then predict the formation factor of the sample at several scales using a Laplace solver. The formation factor is calculated by using a general value of m=2 as cementation factor for intermediate porosity voxels. We compare to experimental measurements of formation factor and porosity both at the small plug and full plug scale and find good agreement.
To assess the degree of uncertainty of the numerical estimate, we probe the extent of heterogeneity by investigating the size of a representative elementary volume (REV) for formation factor. We find that for the considered heterogeneous carbonate sam-ple, formation factor varies considerably over intervals less than a centimetre. Our results show that this variation could be explained by different cementation exponents applied at the micro-voxel scale, with the exemption of one plug, for which the cementation exponent would have to be unreasonably low. These cementation factors are derived by direct comparison be-tween numerical simulation and experiment. We conclude that for one plug an error in experimental measurement might have occurred. The numerical approach presented here therefore aids in quality control. Excluding this plug in the upscaling proce-dure improves the agreement with the experimental result for the whole core while still underestimating formation factor. Al-lowing for a constant m=2 in the simulation at the small scale and using directly the resulting relationship between porosity and formation factor in the upscaling process leads to an overestimation of formation factor.
Thread compound "dope?? in the vernacular, has been used routinely in assembling joints of casing and tubing. The practice in almost universal application in the oil and gas industry involves the manual application of the lubricant in a fashion that is rudimentary, non-systematic and unquantifiable. There is evidence presented in this paper that damage to the near-well zone and other unpleasant events may be associated with the thread compound.
This paper presents the results of both laboratory and field investigations quantifying the effects of the dope on near-well damage. During the assembly of tubing and casing a portion of the thread compound is exuded inside and outside the connection and gets access to the well fluids through the tubing and annular space. Studies presented here show that the dope forms a suspension which penetrates and damages the formation. The studies used standard fluid circulation velocities during typical completion operations.
To characterize and quantify the problem, core samples from the El Tordillo field, with different permeabilities were used. The samples were subjected to the circulation of the suspension created by the thread compound and the completion fluid, measuring the change in the core permeability. The work simulated the well conditions during water injection for water injection wells and during acid treatments for producer wells. A significant reduction in permeability, manifested by a fast and a very large increase in pressure, was measured, at the front face of the core sample. The same measurements showed a far smaller impact in the core body suggesting very minor penetration of dope particles.
This paper describes the laboratory and field work, with description of the test protocols, well conditions and laboratory emulation of field conditions that were used.
A multilateral (MLT) well with an advanced intelligent completion string was recently completed in the Middle East. The well was designed as a "stacked?? dual producer in the upper and lower reservoir, and was drilled using the latest geo-steering techniques to accurately place the wellbore in a highly faulted and geologically complex structure. Rotary-steerable drilling systems (RSS) were used in several of the hole sections, along with advanced logging-while-drilling (LWD) tools including multi-pole acoustic, azimuthal deep resistivity, and resistivity at bit. Encounters with unstable shale and faults made the drilling difficult, but the decisions made in real-time to navigate the well resulted in a very high percentage of net pay in both laterals.
This well combined TAML Level 4 multilateral (MLT) technology with passive inflow control devices in the laterals and an advanced intelligent completion system in the mainbore. The TAML Level 4 multilateral junction was cemented to isolate unstable shale above the reservoir and to provide zonal isolation from the lateral completions, which were compartmentalized into stages with proprietary swellable packers and inflow control devices (ICDs). The intelligent completion was run in the mainbore with two interval control valves (ICVs) and isolation ball valve (LV ICV) to manage the production from each of the two laterals independently. The ICVs and LV ICV are controlled hydraulically through four control lines to surface, which were run in a flat-pack with one electric line to control a downhole gauge package for each lateral. Finally, the well was configured to allow the installation of a large electric submersible pump (ESP) to be run inside the upper 9-5/8-in. production tubing.
This project required intensive planning and coordination for more than a year in advance, which made the project successful despite the difficult drilling conditions and resulted in very little NPT for wellbore construction operations. This paper will focus on the planning, execution and lessons learned from the project.
In the existing horizontal wells in the target sand reservoir of the target field, premature water breakthrough caused the water cut trend to increase within months of production. . This occurred because the reservoir has a very high permeability sands along with active faults containing high viscous reservoir fluids.
New technologies were required to overcome the issue, maximize reservoir contact and enhance a more uniform oil production from a single location. Introducing the smart TAML Level-4 MLT well design to this reservoir along with inflow control device (ICD), inflow control valve (ICV), isolation ball valve (LV ICV) and other downhole gauges proved to be the optimum solution. It also aided in managing the production and the reservoir proactively to achieve maximum oil recovery. Moreover, drilling several laterals from a single wellbore with the ability to control production from both laterals had a great economic advantage because of the optimized cost effective field management.