RESERVOIRS

Field

[0001] The subject disclosure relates to the field of tracking flooding fronts of gas introduced into subterranean reservoirs. More specifically, the subject disclosure relates to techniques for tracking of non-uniform flooding fronts of gas injected in oil reservoirs.

Background

[0002] The use of C0 ₂ or gas as agents for a better recovery of oil in reservoirs has been used for many years. One difficulty is the occurrence of "fingering" resulting from viscous instability between the flooding fluid and the flooded fluid. Fingering is a complex non-linear mechanism difficult to estimate in real conditions and particularly when local reservoir heterogeneities that cannot be captured trigger the phenomenon.

[0003] The fundamental theory of fingering instability occurrence is well explained by the "Buckley-Leverett" model. See I. Brailovski et al. ("Fingering Instability in Water-Oil Displacement," Transport in Porous Media (2006) Vol. 63. pp 363-380, hereinafter, "Brailovski"). The exact shape of fingering when occurring is very difficult to predict.

[0004] Fig. 1 is a series of 2D images 100 of fingering modeled in a Heleshaw tank, from De Wit, A., Bertho, Y. and Martin, M., "Viscous fingering of miscible slices," Physics of Fluids, 17, 054114 (2005). In the series 100, a darker color fluid 110 is being pushed through a lighter color fluid 112, showing "Saffman-Taylor" instabilities and resulting in extended "fingers" such as finger 120, which can be of very long extension. The shapes of fingers are highly variable but as a rule they are shown to be "slender," more or less 2 dimensional (contained in a high permeability layer) as can be seen in the modeled image of Fig. 1. Fig. 2 depicts results of numerical modeling of such instabilities from Brailovski. From the image 200, it can be seen that the instabilities can be of complicated shapes. The extension of the fingers is generally radial and each finger is rather slender. [0005] Due to the difficulty in predicting the complex shapes of the flooding fronts of gas being injected into subterranean oil reservoirs, there is a need for techniques for monitoring such fronts.

Summary

[0006] This summary is provided to introduce a selection of concepts that are further described below in the detailed description. This summary is not intended to identify key or essential features of the claimed subject matter, nor is it intended to be used as an aid in limiting the scope of the claimed subject matter.

[0007] In accordance with some embodiments a system is provided for monitoring an injected gas, such as C0 ₂ gas injected into the subterranean reservoir through one or more injector wellbores. The system includes a downhole gravity measurement tool adapted to be deployed in a borehole and to make gravity measurements at a downhole location; and a processing system adapted and programmed to receive data representing the gravity measurements from the downhole gravity measurement tool and based at least in part thereon to monitor parameters associated with gas introduced into a subterranean reservoir. According to some embodiments the monitored parameters can be: distance from the downhole location to a front associated with the injected C0 ₂ gas, inclination angle associated with the injected C0 ₂ gas and a reference plane, or location of an interface associated with the injected C0 ₂ gas. According to some embodiments, the monitoring can also include identification of shape characteristics of an interface associated with the C0 ₂ gas, such as elongated finger- like shapes. According to some embodiments the monitoring of parameters associated with the introduced gas is based on horizontal components of the gravity measurements.

[0008] According to some embodiments, a method is provided for monitoring gas, such as C0 ₂ gas, injected into a subterranean reservoir. The method includes deploying a downhole gravity measurement tool in a borehole at a downhole location within the subterranean reservoir; making gravity measurements at the downhole location and generating therefrom data representing the measurements; and monitoring parameters associated with gas introduced into the subterranean reservoir based on the data representing the measurements. According to some embodiments a downhole gravity measurement tool is also deployed in one or more other boreholes and the data is compared and selecting based on measured gravity magnitude. According to some embodiments a plan for injecting the C0 ₂ gas is modified based on the identification of shape characteristics such as fingering.

[0009] Further features and advantages will become more readily apparent from the following detailed description when taken in conjunction with the accompanying Drawings.

Brief Description of the Drawings

[0010] The subject disclosure is further described in the detailed description which follows, in reference to the noted plurality of drawings by way of non-limiting examples of embodiments of the subject disclosure, in which like reference numerals represent similar parts throughout the several views of the drawings, and wherein:

[0011] Fig. 1 is a series of 2D images of fingering modeled in a Heleshaw tank as shown in prior art;

[0012] Fig. 2 depicts results of numerical modeling of instabilities as shown in prior art;

[0013] Fig. 3 illustrates a wellsite where gas is being injected and the flooding front is being tracked, according to some embodiments;

[0014] Figs. 4A and 4B are top-view, or horizontal plane diagrams illustrating injected gas flooding from a number of injections wells towards a number of observation wells at two different times, according to some embodiments;

[0015] Fig. 5 is a top-down, horizontal view diagram showing the position of a C0 ₂ density anomaly at two different times, according to some embodiments; and

[0016] Fig. 6 is a flow chart showing processes for tracking a non-uniform injected gas flooding front, according to some embodiments. Detailed Description

[0017] Specific details are given in the following description to provide a thorough understanding of the embodiments. However, it will be understood by one of ordinary skill in the art that the embodiments may be practiced without these specific details. For example, systems, processes, and other elements in the invention may be shown as components in block diagram form in order not to obscure the embodiments in unnecessary detail. In other instances, well-known processes, structures, and techniques may be shown without unnecessary detail in order to avoid obscuring the embodiments. Further, like reference numbers and designations in the various drawings indicate like elements.

[0018] Also, it is noted that individual embodiments may be described as a process which is depicted as a flowchart, a flow diagram, a data flow diagram, a structure diagram, or a block diagram. Although a flowchart may describe the operations as a sequential process, many of the operations can be performed in parallel or concurrently. In addition, the order of the operations may be re-arranged. A process may be terminated when its operations are completed, but could have additional steps not discussed or included in a figure. Furthermore, not all operations in any particularly described process may occur in each embodiment. A process may correspond to a method, a function, a procedure, a subroutine, a subprogram, etc. When a process corresponds to a function, its termination corresponds to a return of the function to the calling function or the main function.

[0019] Furthermore, embodiments of the invention may be implemented, at least in part, either manually or automatically. Manual or automatic implementations may be executed, or at least assisted, through the use of machines, hardware, software, firmware, middleware, microcode, hardware description languages or any combination thereof. When implemented in software, firmware, middleware or microcode, the program code or code segments to perform the required tasks may be stored in a machine readable medium. A processor(s) may perform the required tasks.

[0020] Due to the difficulty in predicting non-uniform flooding patterns in a reservoir, the detection of such patterns relies instead on sufficient information to image the fingering system. Critical C0 ₂ has a density around 468kg/m , representing a good contrast for deep density measurements such as Gravimetry survey.

[0021] According to some embodiments, gravimetric tomography survey data in time lapses is used to characterize slender structures as fingering structures passing in the vicinity of observation wells. According to some embodiments, a characterization that the detected front flooding is not uniform but rather comprises slender structures is used as an input in the design and strategy used to carry out the injection and flooding.

[0022] According to some embodiments, gravity measurements are used to monitor the volumetric sweep of C0 ₂ and estimate the fingering phenomenon if any. The gravity measurements have been found to be effective due to the high density contrast between C0 ₂ and brine or oil which provides better and higher sensitivity for the gravity responses. See E. Gasperikova and G. M. Hoversten. Gravity monitoring of C0 ₂ movement during sequestration: Model studies. Geophysics Vol. 73. No. 6 (2008), which is incorporated by reference herein. According to some embodiments, an analytical forward model is developed and used. A strategy to monitor the movement of the C0 ₂ flooding front and to predict the fingering phenomenon is used to provide a warning tool to improve the sweep efficiency, according to some embodiments.

[0023] As a supercritical fluid, C0 ₂ exhibits different properties depending on pressure and temperature conditions. It has a high mobility ratio. Under miscible conditions, it has an affinity for oil and mixes to form a low viscosity. For the EOR process to be successful, the C0 ₂ should first be in contact with the oil, and then be directed to effectively sweep the reservoir. In heterogeneous oil reservoirs, the fingering phenomenon of C0 ₂ might be formed along high permeability zones as shown in Figs. 1, 4A and 4B. See, e.g., R. Berenblyum, G. Calderon, L. Kollbotn, L. M. Surguchev. Modeling C02 injection: IOR potential after waterflooding, SPE 113436 (2008). Such phenomenon reduces the sweep efficiency and causes poor flood performance. [0024] To monitor the volumetric sweep of C0 ₂ flooding, to predict whether there is any fingering phenomenon and to provide early warning of any potential loss of containment, the gravity measurement is a reasonable candidate to be used in C0 ₂ monitoring application. C0 ₂

3 is less dense and more compressible than brine or oil - its density is around 0.46 g/cm . The density contrast between C0 ₂ and brine or oil is high, and consequently the gravity measurements have a good sensitivity to density change and will have a stronger signal. According to some embodiments, such data are acquired by the borehole gravity meter (BHGM), which runs on wireline cable for reservoir monitoring. The BHGM has been found to provide a large depth of investigation and good vertical resolution at typical reservoir depths. The BHGM could be used in many applications as has been described in several studies. See: J. L. Brady and D.S. Wolcott, Gravity Methods: Useful Techniques for Reservoir Surveillance. (May 1993). SPE 26095; J. L. Brady et al. Improved Production Log Interpretation in Horizontal Wells Using a Combination of Pulsed Neutron Logs, Quantitative Temperature Log Analysis, Time Lapse LWD Resistivity Logs and Borehole Gravity. (May 1998). SPE 46222; K. Hadj-Sassi and J.-M. Donadille, Three-dimensional inversion of borehole gravity measurements for reservoir fluid monitoring. SPE 136928. April 2010; and K. Schultz. Monitoring fluid movement with the borehole gravity meter. Geophysics, Vol. 54, No. 10 (Oct 1989), pp. 1267-1273, 8 Figures).

[0025] Fig. 3 illustrates a wellsite where gas is being injected and the flooding front is being tracked, according to some embodiments. Wellsite 300 includes multiple injection wells, of which injection well 330 is shown, as well as multiple observation/monitoring wells, of which observation/monitoring well 310 is shown. Gas is being injected via the injection well 330 and other injection wells into the reservoir formation 336. The reservoir 336, for example, can be a heterogeneous oil reservoir. The injected gas, according to embodiments, is C0 ₂ gas and has flooded region 332 which includes a non-uniform front 334. The front 334 is being monitored using gravity measurements using a borehole gravity meter 332 that is part of a wireline toolstring 320. The toolstring 320 is deployed in well 310 via wire 312 from logging truck 314. Further details of a suitable borehole gravity meters is provided in U.S. Patent Application Publication Nos. 2011/0185806, and 2011/0191027, both of which are hereby incorporated by reference.

[0026] A data processing unit 350 is included, which according to some embodiments, is located within logging truck 314 and according to other embodiments is partially or fully located at other locations at the wellsite or one or more remote locations. The data processing unit 350 receives the gravity measurements from the gravity meter 322 and calculates therefrom, for example using the forward modeling techniques described herein, tracking information for the flooding front 334. The data processing unit 350 includes one or more central processing units 340, storage system 344, communications and input/output modules 340, a user display 346 and a user input system 348.

[0027] Figs. 4 A and 4B are top-view, or horizontal plane diagrams illustrating injected gas flooding from a number of injections wells towards a number of observation wells at two different times, according to some embodiments. Fig. 4A shows the C0 ₂ flooded zone 420 having a flooding front 422. The C0 ₂ is being injected at injection wells 412, 414, 416 and 418. The observation wells (which can also be production wells) 402, 404, 406 and 408 are used to track the movement of the flooded zone. Fig. 4B shows the C0 ₂ flooding at time T2=T1+At. During the time At, the additional areas 430 and 432 have become flooded with the C0 ₂ . The new flooding front 434 is shown on zone 430. Note that the flooding pattern shown in Fig. 4B has the appearance of viscous fingering.

[0028] The fingering phenomenon is understood to form following geological stratification mostly in the horizontal plane in the oil reservoir. In order to be sensitive to such phenomenon and then monitor its movement over time, the horizontal components of the gravity measurements - § _x and §>y - are estimated for this application. Such components of the gravity data have better potential for tracking this phenomenon, compared to the vertical component of the gravity field § _z . The horizontal components provide better coverage of the information on the horizontal layers. Their sensitivities to the horizontal movement of the C0 ₂ flooding are important compared to that given by the vertical component of the gravity

§ z . Hence, the horizontal components of the gravity measurements can be used to provide information on the fingering phenomenon of the C0 ₂ in an oil reservoir.

[0029] The horizontal components of the gravity field g are the gradients of the potential U in the x and j-directions, z being the vertical direction. Their expressions in the Cartesian coordinate system are described by the following equations:

[0030] U is the gravitational potential which respects the Laplace equation

( ² U(r _s ) = 0) outside of the source of the gravity field, G is the Newtonian gravitational

- 11 2 2

constant, G « 6.67 x 10 N*m /kg , the lengths r _s and r (in meters) respectively represent the observation location (measurement station) and the integration points and V is the body

2 2 2 3

volume. U is expressed in m /s , g is in m/s and p is in kg/m .

[0031] Fig. 5 is a top-down, horizontal view diagram showing the position of a C0 ₂ density anomaly at two different times, according to some embodiments. Injection wells 512, 514, 516 and 518 are shown as is an observation well 504 at station position x _s and y _s . The anomaly is shown at position 520 at time Tl and at position 530 at time T2 (where Τ2=Τ1+Δί). The viscous fingering of the C0 ₂ flooding of interest is described by the rectangular slab density anomaly as shown in Fig. 5. The C0 ₂ density anomaly is the distribution of the density contrast change over an interval of time [Tl, Τ2=Τ1+Δί]. The C0 ₂ movement over the time lapse At, in the x and y directions, is represented by Δχ and Ay, respectively. The distances between the finger front and the x and y axis, intersecting the observation well placed at (x _s , y _s ), are described by D _x and D _y , respectively. D _x o and D _y o (shown in the Fig. 5) represent the initial position of the C0 ₂ fingering front (C0 ₂ flooding) at Time Tl (D _x0 = ϋ _χ +Δχ and D _y0 = D _y +Ay).). The density contrast Δρ is assumed to be constant within the slab.

[0032] To compute the integrals of the equations 1 and 2, we first assume that the C0 ₂ to be injected infinitely along the z-direction of the reservoir (we neglect the vertical effect of the vertical distribution of the C0 ₂ ). Thus, Equation 1 and 2 are integrated along the z direction, varying from - infinity to + infinity, resulting in the Equations 3 and 4 that no longer depend on z. We then take the limits of x integral from D _x to D _x o and the limits of y from D _y to D _y o. The horizontal attractions g and due to the C0 ₂ density anomaly are then expressed as follows:

DxO fDyO

g , ) = 2GApj _Di dxdy

Dy Eq. 3

(x - x _s ) + {y -y _s )

[0033] By integrating the term in brackets of the equation 3, with respect to x first then y, we get the x component of the gravity field §> ' x■

Similarly, the y component of the gravity field g _y is given by:

[0034] As we can see, the horizontal gravity measurements depend on D _x and D _y which represent the distances that separate the fingering front of C(¾ to the x-z and y-z planes that intersect the observation/monitoring well. The initial distances D _x o and D _y o are assumed to be known, as we may have information on the initial position of the C0 ₂ flooding at time Tl . Therefore, by using the forward modeling systems described by the Equations 5 and 6, the horizontal gravity measurements §> _x and §>y will directly provide: (1) the distances D _x and D _y from the observation well to the fingering front of C0 ₂ flooding; (2) the radial distance R

2 2 1/2

from the observation/monitoring well to the C0 ₂ finger front: R=(D _X + D _y ) ; and (3) the angle (inclination) between the x-z plane and the fingering front of the C0 ₂ flooding : a=tan ^Dy/Dx).

[0035] These parameters will provide an accurate estimation on the movement and the behavior of the fingering of C0 ₂ flooding that can be formed in the heterogeneous oil reservoir. Hence, we have an efficient method to track the C0 ₂ movement, which provides information to improve the problem of poor sweep efficiency of C0 ₂ in oil reservoirs.

[0036] Fig. 6 is a flow chart showing processes for tracking a non-uniform injected gas flooding front, according to some embodiments. In process 610, the gravity measurements

§ x and § _y are acquired at one or multi wells surrounding the C0 ₂ movement. This may be done using a borehole gravity meter. In process 612 the measurements § _x acquired at the multi wells are compared. In process 614 the measurements §> _y acquired at the multi wells are compared. In process 616 the maximum magnitude of the gravity response §, _x (max( §, _x )) is compared with that of g, _y (max( §, _y )) and the well position corresponding to the maximum magnitude of the horizontal gravity responses max( §, _x , §, _y ) is selected. In the case max( g )»max( §, _y ), the fingering phenomenon forms along the x-direction. Likewise, in the case max( §, _x )«max( §, _y ), the fingering phenomenon occurs along the y-direction.

[0037] In process 618, once the well is selected, the distances Dx and Dy are computed from that selected observation well to the fingering front of C0 ₂ flooding, using the developed forward modeling described by Equations 5 and 6. In process 620, the radial distance R is determined from the selected observation/monitoring well to the C0 ₂ finger

2 2 1/2

front: R=(Dx + Dy ) . In process 622 the angle (inclination) between the x-z plane and the fingering front of the C0 ₂ flooding: a=arctan(Dy/Dx) is determined.

[0038] Finally, according to some embodiments, in process 624, in cases where there is any appearance of the fingering phenomenon, the C0 ₂ injection pattern/strategy is modified in order to improve the sweep efficiency.

[0039] Although only a few example embodiments have been described in detail above, those skilled in the art will readily appreciate that many modifications are possible in the example embodiments without materially departing from this invention. Accordingly, all such modifications are intended to be included within the scope of this disclosure as defined in the following claims. In the claims, means-plus-function clauses are intended to cover the structures described herein as performing the recited function and not only structural equivalents, but also equivalent structures. Thus, although a nail and a screw may not be structural equivalents in that a nail employs a cylindrical surface to secure wood parts together, whereas a screw employs a helical surface, in the environment of fastening wood parts, a nail and screw may be equivalent structures. It is the express intention of the applicant not to invoke 35 U.S.C. § 1 12, paragraph 6 for any limitations of any of the claims herein, except for those in which the claim expressly uses the words "means for" together with an associated function.