Seawater Breakthrough Monitoring and Reservoir-Model Improvement Using Natural Boron

2021 ◽  
Author(s):  
Yanqing Wang ◽  
Zhe Liu ◽  
Xiang Li ◽  
Shiqian Xu ◽  
Jun Lu

Abstract Natural geochemical data, which refer to the natural ion concentration in produced water, contain important reservoir information, but is seldomly exploited. Some ions were used as conservative tracers to obtain better knowledge of reservoir. However, using only conservative ions can limit the application of geochemical data as most ions are nonconservative and can either interact with formation rock or react with other ions. Besides, mistakenly using nonconservative ion as being conservative may cause unexpected results. In order to further explore the nonconservative natural geochemical information, the interaction between ion and rock matrix is integrated into the reservoir simulator to describe the nonconservative ion transport in porous media. Boron, which is a promising nonconservative ion, is used to demonstrate the application of nonconservative ion. Based on the new model, the boron concentration data together with water production rate and oil production rate are assimilated through ensemble smoother multiple data assimilation (ES-MDA) algorithm to improve the reservoir model. Results indicate that including nonconservative ion data in the history matching process not only yield additional improvement in permeability field, but also can predict the distribution of clay content, which can promote the accuracy of using boron data to determine injection water breakthrough percentage. However, mistakenly regarding nonconservative ion being conservative in the history matching workflow can deteriorate the accuracy of reservoir model.

2021 ◽  
Author(s):  
Yanqing Wang ◽  
Zhe Liu ◽  
Yuchen Zhang ◽  
Jun Lu

Abstract Geochemical data in produced water contain important reservoir information but are seldomly exploited, especially for the nonconservative chemicals. Some conservative chemical data have been integrated in history matching workflow to obtain better knowledge of reservoirs. However, assuming reservoir chemicals being conservative is impractical because most chemicals are involved in interactions with other chemicals or reservoir rock, and mistakenly regarding nonconservative chemicals as being conservative can cause large error. Nevertheless, once the interactions can be accurately described, nonconservative chemical data can be used to obtain more reservoir information. In this work, a new physicochemical model is proposed to describe the transport of natural nonconservative chemicals (barium and sulfate) in porous media. Both physical reactions, such as ion adsorption and desorption, and chemical reactions, such as barite deposition, are integrated. Based on the new model, the ensemble smoother with multiple data assimilations (ES-MDA) method is employed to update reservoir model parameters by assimilating oil production rate, water production rate, and chemical data (barium and sulfate concentration). Data assimilation results show that integrating geochemical data in ES-MDA algorithm yields additional improvements in estimation of permeability. Besides, clay content distribution, which is critical in injection water breakthrough percentage calculation, can be accurately estimated with relative root mean square error (rRMSE) being as small as 0.1. However, mistakenly regarding nonconservative chemicals as conservative can cause large errors in reservoir parameters estimation. Accurately modeling the chemical interactions is crucial for integrating chemical data in history matching algorithm.


2013 ◽  
Vol 28 (04) ◽  
pp. 369-375 ◽  
Author(s):  
Oscar Vazquez ◽  
Ross A. McCartney ◽  
Eric Mackay

2021 ◽  
Vol 5 (1) ◽  
pp. 119-131
Author(s):  
Frzan F. Ali ◽  
Maha R. Hamoudi ◽  
Akram H. Abdul Wahab

Water coning is the biggest production problem mechanism in Middle East oil fields, especially in the Kurdistan Region of Iraq. When water production starts to increase, the costs of operations increase. Water production from the coning phenomena results in a reduction in recovery factor from the reservoir. Understanding the key factors impacting this problem can lead to the implementation of efficient methods to prevent and mitigate water coning. The rate of success of any method relies mainly on the ability to identify the mechanism causing the water coning. This is because several reservoir parameters can affect water coning in both homogenous and heterogeneous reservoirs. The objective of this research is to identify the parameters contributing to water coning in both homogenous and heterogeneous reservoirs. A simulation model was created to demonstrate water coning in a single- vertical well in a radial cross-section model in a commercial reservoir simulator. The sensitivity analysis was conducted on a variety of properties separately for both homogenous and heterogeneous reservoirs. The results were categorized by time to water breakthrough, oil production rate and water oil ratio. The results of the simulation work led to a number of conclusions. Firstly, production rate, perforation interval thickness and perforation depth are the most effective parameters on water coning. Secondly, time of water breakthrough is not an adequate indicator on the economic performance of the well, as the water cut is also important. Thirdly, natural fractures have significant contribution on water coning, which leads to less oil production at the end of production time when compared to a conventional reservoir with similar properties.


SPE Journal ◽  
2021 ◽  
pp. 1-25
Author(s):  
Ahmed Adila ◽  
Emad W. Al-Shalabi ◽  
Waleed AlAmeri

Summary Engineered water injection (EWI) has gained popularity as an effective technique for enhancing oil recovery. Surfactant flooding is also a well-established chemical enhanced-oil-recovery (EOR) technique in the petroleum industry. The hybrid surfactant flooding/EWI (surfactant/EWI) technique has been studied experimentally and showed promising results. However, there are very limited numerical applications on the hybrid surfactant/EWI technique in carbonates in the literature. Also, the studies applied under harsh conditions of high temperature and high salinity are even fewer. In this study, a numerical-simulation model is developed and used to investigate the hybrid effect of surfactant/EWI in carbonates under harsh conditions. This developed model was validated by history matching a recently conducted surfactant coreflood in the secondary mode of injection. Oil recovery, pressure drop, and surfactant-concentration data were used. The surfactant-flooding model was then coupled with a geochemical model that captures different reactions involved during EWI. The geochemical reactions considered include aqueous, dissolution/precipitation, and ion-exchange reactions. The proposed model has been further validated by history matching another experimental data set. Furthermore, different simulation scenarios were considered, including waterflooding, surfactant flooding, EWI, and the hybrid surfactant/EWI technique. For the case of EWI, wettability alteration was considered as the main mechanism underlying incremental oil recovery. However, both wettability alteration and interfacial-tension (IFT) reduction mechanisms were considered for surfactant flooding depending on the type of surfactant used. The results showed that for the hybrid surfactant/EWI, wettability alteration is considered as the controlling mechanism where surfactant boosts oil-recovery rate through increasing oil relative permeability while EWI reduces residual oil. Moreover, the simulation runs showed that the hybrid surfactant/EWI is a promising technique for enhancing oil recovery from carbonates under harsh conditions. Also, hybrid surfactant/EWI results in a more water-wetting rock condition compared with that of EWI alone, which leads to lower injectivity, and hence, lower rate of propagation for ion-concentration waves. The hybrid surfactant/EWI outperformed other injection techniques followed by EWI, then surfactant flooding, and finally waterflooding. This work gives more insight into the application of hybrid surfactant/EWI on enhancing oil recovery from carbonates. The novelty is further highlighted through applying the hybrid surfactant/EWI technique under harsh conditions. In addition, the findings of this study can help in better understanding the mechanism behind enhancing oil recovery using the hybrid surfactant/EWI technique and the important parameters needed to model its effect on oil recovery.


2015 ◽  
Vol 18 (04) ◽  
pp. 564-576 ◽  
Author(s):  
O.. Vazquez ◽  
C.. Young ◽  
V.. Demyanov ◽  
D.. Arnold ◽  
A.. Fisher ◽  
...  

Summary Produced-water-chemistry (PWC) data are the main sources of information to monitor scale precipitation in oilfield operations. Chloride concentration is used to evaluate the seawater fraction of the total produced water per producing well and is included as an extra history-matching constraint to reevaluate a good conventionally history-matched (HM) reservoir model for the Janice field. Generally, PWC is not included in conventional history matching, and this approach shows the value of considering the nature of the seawater-injection front and the associated brine mixing between the distinctive formation water and injected seawater. Adding the extra constraint resulted in the reconceptualization of the reservoir geology between a key injector and two producers. The transmissibility of a shale layer is locally modified within a range of geologically consistent values. Also, a major lineament is identified which is interpreted as a northwest/southeast-trending fault, whereby the zero transmissibility of a secondary shale in the Middle Fulmar is locally adjusted to allow crossflow. Both uncertainties are consistent with the complex faulting known to exist in the region of the targeted wells. Other uncertainties that were carried forward to the assisted-history-matching phase included water allocation to the major seawater injectors; thermal fracture orientation of injectors; and the vertical and horizontal permeability ratio (Kv/Kh) of the Fulmar formation. Finally, a stochastic particle-swarm-optimization (PSO) algorithm is used to generate an ensemble of HM models with seawater fraction as an extra constraint in the misfit definition. Use of additional data in history matching has improved the original good HM solution. Field oil-production rate is interpreted as improved over a key period, and although no obvious improvement was observed in field water-production rate, seawater fraction in a number of wells was improved.


2021 ◽  
Author(s):  
Yanqing Wang ◽  
Xiang Li ◽  
Jun Lu

Abstract Seawater injection is widely used to maintain the offshore oil reservoir pressure and improve oil recovery. However, injecting seawater into reservoirs can cause many issues such as reservoir souring and scaling, which are tightly related to the seawater breakthrough percentage. Accurately calculating the seawater breakthrough percentage is important for estimating the severity of those problems and further developing corresponding strategies to solve those issues. The validation of using natural ion boron as tracer to calculate seawater breakthrough percentage was investigated. Boron can interact with clays, which can influence the accuracy in seawater breakthrough calculation using boron. Therefore, the interaction between boron and different clays at various conditions were studied, and Freundlich adsorption equation was used to describe the boron adsorption isotherms. Then boron adsorption isotherms were coupled into the reservoir simulator to investigate the boron transportation in porous media, and the results in turn were further analyzed to calculate the accurate seawater breakthrough percentage. Results indicated that boron adsorption by different clays varied. pH value of solution can significantly influence the amount of boron adsorbed. As a result, the boron concentration profile was delayed in coreflood test. The results of reservoir model fit perfectly with that of coreflood test, indicating the validation of boron reaction model. Based on the reservoir simulator results, boron concentration profile in produced water was successfully used to calculate seawater breakthrough percentage by considering the clay content distribution. However, the seawater breakthrough point cannot be calculated by boron as boron concentration is still at the formation level due to boron desorption.


2018 ◽  
Vol 6 (3) ◽  
Author(s):  
Anietie Okon ◽  
Dulu Appah ◽  
Julius Akpabio

In the Niger Delta, available correlations to predict water breakthrough time in thin oil rim reservoirs are based on generic reservoir models and/or experimental design approach. This approach had not considered the heterogeneity of the reservoir. Thus, the prediction of these available correlations for thin oil rim reservoirs in the Niger Delta is in doubt, considering the sensitive nature of developing thin oil rim reservoirs. Then, a correlation for water breakthrough time (tbt) was developed based on integrated reservoir model of thin oil rim reservoir in the Niger Delta. The obtained result indicated that the developed correlation predicted 1652.72 days compared to the actual Oilfield breakthrough time of 1653 days (about 4.53 years). Also, sensitivity study showed that the developed correlation and the integrated reservoir model predictions of oil production rate (qo), fractional well penetration (hp/h) and height above perforation-oil column (hap/h) on the water breakthrough time (tbt) were close and resulted in coefficient of determination (R2) of 0.9697, 0.8597 and 0.9553, respectively. Furthermore, the results depicted that water coning breakthrough time (tbt) depends directly on oil production rate (q) and well completion parameters: fractional well penetration (hp/h) and height above perforation (hap). Hence, to delay early water breakthrough in thin oil rim reservoirs, these completion parameters are consideration in vertical wells to achieve optimum oil recovery. Also, the developed correlation can be used as a quick and robust tool to predict water breakthrough time of thin oil rim reservoirs in the Niger Delta.


2021 ◽  
Vol 73 (04) ◽  
pp. 60-61
Author(s):  
Chris Carpenter

This article, written by JPT Technology Editor Chris Carpenter, contains highlights of paper SPE 199149, “Rate-Transient-Analysis-Assisted History Matching With a Combined Hydraulic Fracturing and Reservoir Simulator,” by Garrett Fowler, SPE, and Mark McClure, SPE, ResFrac, and Jeff Allen, Recoil Resources, prepared for the 2020 SPE Latin American and Caribbean Petroleum Engineering Conference, originally scheduled to be held in Bogota, Colombia, 17–19 March. The paper has not been peer reviewed. This paper presents a step-by-step work flow to facilitate history matching numerical simulation models of hydraulically fractured shale wells. Sensitivity analysis simulations are performed with a coupled hydraulic fracturing, geomechanics, and reservoir simulator. The results are used to develop what the authors term “motifs” that inform the history-matching process. Using intuition from these simulations, history matching can be expedited by changing matrix permeability, fracture conductivity, matrix-pressure-dependent permeability, boundary effects, and relative permeability. Introduction This article, written by JPT Technology Editor Chris Carpenter, contains highlights of paper SPE 199149, “Rate-Transient-Analysis-Assisted History Matching With a Combined Hydraulic Fracturing and Reservoir Simulator,” by Garrett Fowler, SPE, and Mark McClure, SPE, ResFrac, and Jeff Allen, Recoil Resources, prepared for the 2020 SPE Latin American and Caribbean Petroleum Engineering Conference, originally scheduled to be held in Bogota, Colombia, 17-19 March. The paper has not been peer reviewed. This paper presents a step-by-step work flow to facilitate history matching numerical simulation models of hydraulically fractured shale wells. Sensitivity analysis simulations are performed with a coupled hydraulic fracturing, geomechanics, and reservoir simulator. The results are used to develop what the authors term “motifs” that inform the history-matching process. Using intuition from these simulations, history matching can be expedited by changing matrix permeability, fracture conductivity, matrix-pressure-dependent permeability, boundary effects, and relative permeability. Introduction The concept of rate transient analysis (RTA) involves the use of rate and pressure trends of producing wells to estimate properties such as permeability and fracture surface area. While very useful, RTA is an analytical technique and has commensurate limitations. In the complete paper, different RTA motifs are generated using a simulator. Insights from these motif simulations are used to modify simulation parameters to expediate and inform the history- matching process. The simulation history-matching work flow presented includes the following steps: 1 - Set up a simulation model with geologic properties, wellbore and completion designs, and fracturing and production schedules 2 - Run an initial model 3 - Tune the fracture geometries (height and length) to heuristic data: microseismic, frac-hit data, distributed acoustic sensing, or other diagnostics 4 - Match instantaneous shut-in pressure (ISIP) and wellhead pressure (WHP) during injection 5 - Make RTA plots of the real and simulated production data 6 - Use the motifs presented in the paper to identify possible production mechanisms in the real data 7 - Adjust history-matching parameters in the simulation model based on the intuition gained from RTA of the real data 8 -Iterate Steps 5 through 7 to obtain a match in RTA trends 9 - Modify relative permeabilities as necessary to obtain correct oil, water, and gas proportions In this study, the authors used a commercial simulator that fully integrates hydraulic fracturing, wellbore, and reservoir simulation into a single modeling code. Matching Fracturing Data The complete paper focuses on matching production data, assisted by RTA, not specifically on the matching of fracturing data such as injection pressure and fracture geometry (Steps 3 and 4). Nevertheless, for completeness, these steps are very briefly summarized in this section. Effective fracture toughness is the most-important factor in determining fracture length. Field diagnostics suggest considerable variability in effective fracture toughness and fracture length. Typical half-lengths are between 500 and 2,000 ft. Laboratory-derived values of fracture toughness yield longer fractures (propagation of 2,000 ft or more from the wellbore). Significantly larger values of fracture toughness are needed to explain the shorter fracture length and higher net pressure values that are often observed. The authors use a scale- dependent fracture-toughness parameter to increase toughness as the fracture grows. This allows the simulator to match injection pressure data while simultaneously limiting fracture length. This scale-dependent toughness scaling parameter is the most-important parameter in determining fracture size.


2021 ◽  
Author(s):  
Hui Zhao ◽  
Wei Liu ◽  
Xiang Rao ◽  
Guanglong Sheng ◽  
Huazhou Andy Li ◽  
...  

Abstract The data-driven interwell simulation model (INSIM) has been recognized as an effective tool for history matching and interwell-connectivity characterization of waterflooding reservoirs. INSIM-FT-3D (FT: front tracking) was recently developed to upgrade the applicationdimension of INSIM series data-driven models from two-dimensional (2D) to three-dimensional (3D). However, INSIM-FT-3D cannot accurately infer the dynamic change of well-connectivity and predict well's bottom-hole pressure (BHP). The main purpose of this study intends to expand the capability of INSIM-FT-3D to empower for the assimilation of BHPs, the reliable prediction of water breakthrough and the characterization of dynamic interwell-connectivities. The default setting of well index (WI) in INSIM-FT-3D based on Peaceman's equation does not yield accurate BHP estimates. We derive a WI that can honor the BHPs of a reference model composed of a set of 1D connections. When history matching BHPs of a 3D reservoir, we show that the derived WI is a better initial guess than that obtained from Peaceman's equation. We also develop a flow-path-tracking (FPT) algorithm to calculate the dynamic interwell properties (allocation factors and pore volumes (PVs)). Besides, we discuss the relationship between the INSIM-family methods and the traditional grid-based methods, which indicates that the INSIM-family methods can calculate the transmissibility of the connection between coarse-scale cells in a more accurate manner. As an improvement of INSIM-FT-3D, the newly proposed data-driven model is denoted as INSIM-FPT-3D. To verify the correctness of the derived WI, we present a 1D problem and a T-shaped synthetic reservoir simulation model as the reference models. BHPs and oil production rates are obtained as the observed data by running these two reference models with total injection/production-rate controls. An INSIM-FPT-3D model is created by specifying the transmissibilities and PVs that are the same as those in the reference model. By applying the derived WIs in INSIM-FPT-3D, the resulting BHPs and oil rates obtained agree well with the reference model without further model calibration. Applying INSIM-FPT-3D to a synthetic multi-layered reservoir shows that we obtain a reasonable match of both BHPs and oil rates with INSIM-FPT-3D. Compared with the FrontSim model, the INSIM-FPT-3D model after history matching is shown to match the dynamic PVs from FrontSim reasonably well and can correctly predict the timing of water breakthrough. By allowing for the assimilation of BHP data, we enable INSIM-FPT-3D to history match a green field with limited production history and forecast the timing of water breakthrough. The improved INSIM-FPT-3D leads to more accurate characterization of the interwell connectivities.


Sign in / Sign up

Export Citation Format

Share Document