Optimal Waterflood Management Using Rate Control

2007 ◽  
Vol 10 (05) ◽  
pp. 539-551 ◽  
Author(s):  
Ahmed Alhuthali ◽  
Adedayo Oyerinde ◽  
Akhil Datta-Gupta

Summary Field-scale rate optimization problems often involve highly complex reservoir models, production-and-facilities related constraints, and a large number of unknowns. These factors make optimal reservoir management through rate- and flood-front control difficult without efficient optimization tools. Some aspects of the optimization problem have been studied before mainly using an optimal control theory. However, the applications to date have been rather limited to small problems because of the computation time and the complexities associated with the formulation and solution of adjoint equations. Field-scale rate optimization for maximizing waterflood sweep efficiency under realistic field conditions has remained largely unexplored. This paper proposes a practical and efficient approach for computing optimal injection and production rates, thereby managing the waterflood front to maximize sweep efficiency and delaying the arrival time to minimize water cycling. Our work relies on equalizing the arrival times of the waterflood front at all producers within selected subregions of a waterflood project. The arrival-time optimization has favorable quasilinear properties, and the optimization proceeds smoothly even if our initial conditions are far from the solution. Furthermore, the sensitivity of the arrival time with respect to injection and production rates can be calculated analytically using a single-flow simulation. This makes our approach computationally efficient and suitable for large-scale field applications. The arrival time optimization ensures appropriate rate allocation and flood-front management by delaying the water breakthrough at the producing wells. Several examples are presented to support the robustness and efficiency of the proposed optimization scheme. These include several 2D-synthetic examples for validation purposes and a 3D field application. In addition, we demonstrate the potential of the approach to optimize the flow profile along injection/production segments of horizontal-smart wells. Introduction Waterflooding is by far the most commonly used method to improve oil recovery after primary depletion. In spite of its many favorable characteristics, reservoir heterogeneity—particularly permeability contrast—can have an adverse impact on the performance of waterflooding. The presence of high-permeability streaks can severely reduce the sweep efficiency, leading to an early water arrival at the producers and bypassed oil. Also, an increased cost is associated with water recycling and handling. One approach to counteract the impact of heterogeneity and improve waterflood sweep efficiency is optimal rate allocation to the injectors and producers (Asheim 1988; Sudaryanto and Yortsos 2001; Brouwer et al. 2001; Brouwer and Jansen 2004; Grinestaff 1999; Grinestaff and Caffrey 2000). Through optimal rate control, we can manage the propagation of the flood front, delay water breakthrough at the producers, and also increase the recovery efficiency. Previous efforts to optimize waterflooding relied on optimal control theorem to allocate injection/production rates for fixed well configurations. Asheim (1988) investigated the optimization of waterflood based on maximizing net present value (NPV) for multiple vertical injectors and one producer where the rate profiles change throughout the optimization time. Sudaryanto and Yortsos (2001) used maximizing the displacement efficiency at water breakthrough as the objective for the optimization with two injectors and one producer. The optimal injection policy was found to be bang bang type. That is, the injectors were operated only at their extreme values—either at the maximum allowable injection rate or fully shut. The optimization then involved finding the switch time between the two injectors to ensure simultaneous water arrival at the producing well. Brouwer et al. (2001) studied the static optimization of waterflooding with two horizontal smart wells containing permanent downhole well-control valves and measurement equipment. The static optimization implies that the flow rates of the inflow-control valves (ICVs) along the well segments were kept constant during the waterflooding process until the water arrived at the producer. Various heuristic algorithms were utilized to minimize the impact of high-permeability streaks on the waterflood performance through rate control. The results indicated that the optimal rate allocation can be obtained by reducing the distribution of water-arrival times at various segments along the producer. Subsequently, Brouwer and Jansen (2004) extended their work to dynamic optimization of waterflooding with smart wells using the optimal control theory. The optimization was performed on one horizontal producer and one horizontal injector. Each well is equipped with 45 ICVs. The objective was to maximize the NPV, and it was achieved through changing the rate profile along the well segments throughout the optimization period. Both rate-constrained and bottomhole-pressure-constrained well conditions were studied.

2010 ◽  
Vol 13 (03) ◽  
pp. 406-422 ◽  
Author(s):  
Akhil Datta-Gupta ◽  
Ahmed H.H. Alhuthali ◽  
Bevan Yuen ◽  
Jerry Fontanilla

2009 ◽  
Author(s):  
Ahmed Humaid H. Alhuthali ◽  
Akhil Datta-Gupta ◽  
Bevan Bun Wo Yuen ◽  
Jerry Pasco Fontanilla

2021 ◽  
Vol 13 (15) ◽  
pp. 3014
Author(s):  
Feng Wang ◽  
Dongkai Yang ◽  
Guodong Zhang ◽  
Jin Xing ◽  
Bo Zhang ◽  
...  

Sea surface height can be measured with the delay between reflected and direct global navigation satellite system (GNSS) signals. The arrival time of a feature point, such as the waveform peak, the peak of the derivative waveform, and the fraction of the peak waveform is not the true arrival time of the specular signal; there is a bias between them. This paper aims to analyze and calibrate the bias to improve the accuracy of sea surface height measured by using the reflected signals of GPS CA, Galileo E1b and BeiDou B1I. First, the influencing factors of the delay bias, including the elevation angle, receiver height, wind speed, pseudorandom noise (PRN) code of GPS CA, Galileo E1b and BeiDou B1I, and the down-looking antenna pattern are explored based on the Z-V model. The results show that (1) with increasing elevation angle, receiver height, and wind speed, the delay bias tends to decrease; (2) the impact of the PRN code is uncoupled from the elevation angle, receiver height, and wind speed, so the delay biases of Galileo E1b and BeiDou B1I can be derived from that of GPS CA by multiplication by the constants 0.32 and 0.54, respectively; and (3) the influence of the down-looking antenna pattern on the delay bias is lower than 1 m, which is less than that of other factors; hence, the effect of the down-looking antenna pattern is ignored in this paper. Second, an analytical model and a neural network are proposed based on the assumption that the influence of all factors on the delay bias are uncoupled and coupled, respectively, to calibrate the delay bias. The results of the simulation and experiment show that compared to the meter-level bias before the calibration, the calibrated bias decreases the decimeter level. Based on the fact that the specular points of several satellites are visible to the down-looking antenna, the multi-observation method is proposed to calibrate the bias for the case of unknown wind speed, and the same calibration results can be obtained when the proper combination of satellites is selected.


Stroke ◽  
2021 ◽  
Vol 52 (Suppl_1) ◽  
Author(s):  
Brian Drumm ◽  
Paul Bentley ◽  
Zoe Brown ◽  
Lucio D’Anna ◽  
Tsering Dolkar ◽  
...  

Introduction: There are reports of changes in the numbers of stroke admissions and time intervals to receiving emergency treatments during the COVID-19 pandemic. We examined the impact of the COVID-19 pandemic on the stroke thrombolysis rate and delay to thrombolysis treatment in a regional stroke centre in London, UK. Methods: COVID-19 testing began at our hospital on 3 March 2020. Clinical data for all acute stroke admissions were routinely collected as part of a national Sentinel Stroke National Audit Programme (SSNAP) and all thrombolysis data were entered into our local thrombolysis database. We retrospectively extracted the relevant patient data for the period of March to May 2020 (COVID group) and compared to the same period in 2019 (pre-COVID group). Results: Compared with pre-COVID, there was a 17.5% fall in total stroke admissions (from 315 to 260) during COVID; but there were no significant differences in the demographics, stroke severity, proportions with known time of onset, or median onset-to-arrival time. The thrombolysis rates amongst ischemic strokes were not significantly different between the two groups (59/260=23% pre-COVID vs. 41/228=18% COVID, p=.19). For thrombolysis patients, their stroke severity and demographics were similar between the two both groups. Median onset-to-needle time was significantly longer by 22 minutes during COVID [127 (IQR 94-160) vs. 149 (IQR 110-124) minutes, p=.045]; this delay to treatment was almost entirely due to a longer median onset-to-arrival time by 16 minutes during COVID (p=.029). Favorable early neurological outcomes post-thrombolysis (defined as an improvement in NIHSS by ≥4 points at 24 hours) were similar (45% vs. 46%, p=.86). Conclusion: COVID-19 pandemic had a negative impact on prehospital delays which in turn significantly increased onset-to-needle time, but without affecting the chance of a favorable early neurological outcome. Our data highlight the need to maintain public awareness of taking immediate action when stroke symptoms occur during the COVID-19 pandemic.


2021 ◽  
Author(s):  
Tsubasa Onishi ◽  
Hongquan Chen ◽  
Jiang Xie ◽  
Shusei Tanaka ◽  
Dongjae Kam ◽  
...  

Abstract Streamline-based methods have proven to be effective for various subsurface flow and transport modeling problems. However, the applications are limited in dual-porosity and dual-permeability (DPDK) system due to the difficulty in describing interactions between matrix and fracture during streamline tracing. In this work, we present a robust streamline tracing algorithm for DPDK models and apply the new algorithm to rate allocation optimization in a waterflood reservoir. In the proposed method, streamlines are traced in both fracture and matrix domains. The inter-fluxes between fracture and matrix are described by switching streamlines from one domain to another using a probability computed based on the inter-fluxes. The approach is fundamentally similar to the existing streamline tracing technique and can be utilized in streamline-assisted applications, such as flow diagnostics, history matching, and production optimization. The proposed method is benchmarked with a finite-volume based approach where grid-based time-of-flight was obtained by solving the stationary transport equation. We first validated our method using simple examples. Visual time-of-flight comparisons as well as tracer concentration and allocation factors at wells show good agreement. Next, we applied the proposed method to field scale models to demonstrate the robustness. The results show that our method offers reduced numerical artifacts and better represents reservoir heterogeneity and well connectivity with sub-grid resolutions. The proposed method is then used for rate allocation optimization in DPDK models. A streamline-based gradient free algorithm is used to optimize net present value by adjusting both injection and production well rates under operational constraints. The results show that the optimized schedule offers significant improvement in recovery factor, net present value, and sweep efficiency compared to the base scenario using equal rate injection and production. The optimization algorithm is computationally efficient as it requires only a few forward reservoir simulations.


2018 ◽  
Vol 3 (2) ◽  
pp. 170-178
Author(s):  
Lidia Agustina Rumaal ◽  
Jehunias L. Tanesib ◽  
Jonshon Tarigan

Abstrak Telah dilakukan pemetaan daerah rawan tsunami berdasarkan estimasi waktu tiba gelombang dan tutupan lahan di Kabupaten Kupang Provinsi Nusa Tenggara Timur menggunakan aplikasi Penginderaan Jauh dan Sistem Informasi Geografi. Penelitian ini bertujuan untuk mengidentifikasi, memetakan daerah rawan tsunami dan tingkat kerawanannya menurut estimasi waktu tiba gelombang dan tutupan lahan sebagai upaya mitigasi dampak bencana tsunami terhadap kepadatan penduduk. Metode penelitian secara umum dibagi dalam empat tahap utama yaitu pembangunan basis data berupa pembuatan peta tutupan lahan, peta gempa dan peta batimetri. Analisis data kerawanan dari peta tutupan lahan dan etimasi waktu tiba gelombang, penyajian hasil data dalam bentuk tingkat kerawanan masing-masing peta dan analisis hasil penelitian berupa tingkat kerawanan secara kualitatif masing-masing daerah titik pantau menurut peta tutupan lahan maupun estimasi waktu tiba gelombang. Selain itu, dampak kerawanan tsunami diklasifikasikan menurut tingkat kepadatan penduduk untuk kebutuhan mitigasi sebagai berikut Kecamatan Kupang Timur, Kupang Barat, Sulamu, Amfoang Timur, Semau, Semau Selatan, Amfoang Utara, Amfoang Barat Daya, Amfoang Barat Laut dan Fatuleu Barat. Kata kunci : Peta rawan tsunami, Penginderaan Jauh, Sistem Informasi Geografi, Estimasi Waktu Tiba Gelombang  Abstract Mapping of hazard tsunami areas based on estimation of arrival time of wave and land cover in Kupang Regency of East Nusa Tenggara Province using remote sensing application and geographic information system has been done. The  aims of this research are to mapping the hazard tsunami area and tsunami vulnerability level in Kupang Regency East Nusa Tenggara according to the estimated arrival time of the wave and land cover as an effort to mitigate the impact of the tsunami disaster on population density. These generally devided into four main phase namely development of database in the form of land cover map , seismic maps and bathymetry maps, data analysis of research results in the form of qualitative vulnerability of each monitoring area according to land cover map and estimated wave arrival time. Presentation of data results in the form of vulnerability level of each map and analysis and results analysis of research the form of vulnerability level of each map and analysis of research results in the form of qualitative vulnerability of each monitoring area according to land cover map and estimated wave arrival time. And then, the impact of tsunami vulnerability is classified according to population density levels for mitigation needs as follows Kupang Timur, Kupang Barat, Sulamu, Amfoang Timur, Semau, Semau Selatan, Amfoang Utara, Amfoang Barat Daya, Amfoang Barat Laut and Fatuleu Barat. Keywords: Tsunami Hazard Map, Remote Sensing, Geographic Information System, Estimated Time of arrival Wave


Sign in / Sign up

Export Citation Format

Share Document