METHODOLOGY FOR HPHISO QUALITY MAP FILTERED BY TRANSMISSIBILITY AND SCAN QUALITY FOR POST-PROCESSING OF OIL RESERVOIR FLOW SIMULATIONS

Abstract

The invention presents a methodology for evaluating the drainage efficiency of the proposed drainage mesh, by identifying poorly drained regions for mesh optimization and/or implementation of complementary projects, in addition to qualitatively evaluating the drainage effect of a field with different number of wells and units and, finally, comparatively evaluating the drainage quality from different reservoirs, thus observing opportunities for optimizing production management, such as activating different intelligent completion valves, for example. To achieve this objective, it applies a set of filters to the mobile HPHISO map in order to identify the regions that have the greatest potential to be explored and that will bring greater recovery gains for a given oil field.

Claims

1. A METHOD FOR THE HPHISO QUALITY MAP FILTERED BY TRANSMISSIBILITY, characterized in that it comprises following steps: a. performing scale cuts, and inform the top and bottom layers of the reservoir modeled in finite volumes, b. calculating the porous volume above the COA using the equations ${\mathrm{BVPacimaCOA}}_{\mathrm{ij}}=\{\begin{array}{c}0,{\mathrm{Sw}}_{\mathrm{ijk}}\left(t_i\right)>0,90\\ {\mathrm{BLOCKVOL}}_{\mathrm{ijk}},{\mathrm{Sw}}_{\mathrm{ijk}}\left(t_i\right)<0,90\end{array}{\mathrm{BLOCKVOL}}_{\mathrm{ijk}}=H_{\mathrm{ijk}}*{\mathrm{NTG}}_{\mathrm{ijk}}*{\mathrm{PHI}}_{\mathrm{ijk}}$  and c. calculating the transmissibility applied in each cell of the model
KH.sub.ijk=H.sub.ijk*√{square root over (PERMI.sub.ijk*PERMJ.sub.ijk)}, d. applying transmissibility filter on each model cell ${\mathrm{filtroKHcc}}_{\mathrm{ijk}}=\{\begin{array}{cc}1,& {\mathrm{KH}}_{\mathrm{ijk}}>\mathrm{khCorte}\\ 0,& {\mathrm{KH}}_{\mathrm{ijk}}<\mathrm{khCorte}\end{array}$ e. calculating the weighted average of the transmissibility filter by the porous volume by means of the equation ${\mathrm{filtroHK}}_{\mathrm{ij}}=\frac{{.Math.}_k{\mathrm{filtroKHcc}}_{\mathrm{ijk}}*{\mathrm{BVPacimaCOA}}_{\mathrm{ijk}}}{{.Math.}_k{\mathrm{BVPacimaCOA}}_{\mathrm{ijk}}}$ applied to kBase and kTopo limits

2. THE METHOD FOR THE HPHISO QUALITY MAP FILTERED BY TRANSMISSIBILITY according to claim 1, characterized in that, if the cell transmissibility is lower than a given limit, defined by the user, all the oil is zeroed for the filtered HPHISO map.

3. THE METHOD FOR THE HPHISO QUALITY MAP FILTERED BY TRANSMISSIBILITY according to claim 1, characterized in that the filter is applied to each cell of the grid.

4. THE METHOD FOR THE HPHISO QUALITY MAP FILTERED BY TRANSMISSIBILITY according to claim 3, characterized in that it is later integrated in the depth, generating an average value along the reservoir defined by the user.

5. THE METHOD FOR THE HPHISO QUALITY MAP FILTERED BY TRANSMISSIBILITY according to claim 1, characterized in that the user is also responsible for defining a recovery factor value where there is no fluid flow.

6. A METHOD FOR THE HPHISO QUALITY MAP FILTERED BY SCAN QUALITY, characterized in that it comprises the following steps: a. performing scale cuts, and inform the top and bottom layers of the reservoir modeled in finite volumes, b. calculating the porous volume above the COA using the equations ${\mathrm{BVPacimaCOA}}_{\mathrm{ij}}=\{\begin{array}{c}0,{\mathrm{Sw}}_{\mathrm{ijk}}\left(t_i\right)>0,90\\ {\mathrm{BLOCKVOL}}_{\mathrm{ijk}},{\mathrm{Sw}}_{\mathrm{ijk}}\left(t_i\right)<0,90\end{array}{\mathrm{BLOCKVOL}}_{\mathrm{ijk}}=H_{\mathrm{ijk}}*{\mathrm{NTG}}_{\mathrm{ijk}}*{\mathrm{PHI}}_{\mathrm{ijk}}$  and c. calculating the FR of each cell ${\mathrm{FRcc}}_{\mathrm{ijk}}=\frac{{\mathrm{So}}_{\mathrm{ijk}}\left(\mathrm{ti}\right)-{\mathrm{So}}_{\mathrm{ijk}}}{{\mathrm{So}}_{\mathrm{ijk}}\left(\mathrm{ti}\right)}$ d. applying a weighted average in each cell by the pore volume calculated in step ii using the equation ${\mathrm{FRmedio}}_{\mathrm{ij}}=\frac{{.Math.}_k{\mathrm{FRcc}}_{\mathrm{ijk}}*{\mathrm{BVPacimaCOA}}_{\mathrm{ijk}}}{{.Math.}_k{\mathrm{BVPacimaCOA}}_{\mathrm{ijk}}}$ e. applying the scan quality filter on the FRmedio weighted property ${\mathrm{filtroFR}}_{\mathrm{ij}}=-\left(\frac{1}{\mathrm{frCorte}}\right)*{\mathrm{FRmedio}}_{\mathrm{ij}}+1$

7. A METHOD FOR THE HPHISO QUALITY MAP FILTERED BY SCAN QUALITY, characterized in that it penalizes the regions that present a level of oil recovery above a certain economic cut-off value in the recovery, theoretical recovery factor theoretical FR, selected from a fractional flow fw.

8. THE METHOD FOR THE HPHISO QUALITY MAP FILTERED BY SCAN QUALITY according to claim 6, characterized in that it is estimated from a decreasing linear relationship assuming 1 as the value for regions with null FR and 0 for regions with FR greater than the limit defined by the user.

9. THE METHOD FOR THE HPHISO QUALITY MAP FILTERED BY TRANSMISSIBILITY AND SCAN QUALITY according to claims 1 and 6, characterized in that it comprises the following steps: a. by means of the porous volume, calculate the mobile HPHISO given by the equation
HPHISOm.sub.ij=Σ.sub.kH.sub.ij*PHI.sub.ij*SOm.sub.ij b. by means of the transmissibility filter method and the scan quality filter method and the mobile HPHISO, calculate the quality map of the filtered HPHISO using the equation
HPHISOm.sub.ij.sup.filtrado=HPHISOm.sub.ij*filtroFR.sub.ij*filtroKH.sub.ij

10. THE METHOD FOR THE HPHISO QUALITY MAP FILTERED BY TRANSMISSIBILITY AND SCAN QUALITY according to claim 10, characterized in that it generates the HPHISO Quality Map filtered as Indicators of Poorly Drained Regions and Potential for Mesh Optimization.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0009] The present invention will be described in more detail below, with reference to the attached figures which, in a schematic way and not limiting the inventive scope, represent examples of its embodiment. From FIGS. 1 to 12, there is the workflow of the proposed methodology, as shown below:

[0010] FIG. 1 illustrates the permo-porous quality filter (transmissibility).

[0011] FIG. 2 illustrates the theoretical FR vs. fractional flow, identifying economic mobile oil saturation.

[0012] FIG. 3 illustrates the scan quality filter.

[0013] FIG. 4 illustrates the porous volume map above the COA.

[0014] FIG. 5 illustrates the transmissibility map (equation 3).

[0015] FIG. 6 illustrates the map of the applied transmissibility filter (filtroKH) (equation 4).

[0016] FIG. 7 illustrates the averaged map of the transmissibility filter (equation 5).

[0017] FIG. 8 illustrates the FR map per cell (equation 6).

[0018] FIG. 9 illustrates the average FR map along column k (given by equation 7).

[0019] FIG. 10 illustrates the scan quality filter map (filtroFR) (equation 8).

[0020] FIG. 11 illustrates the mobile HPHISO map (equation 9).

[0021] FIG. 12 illustrates the Filtered mobile HPHISO map, which is the final product of this invention (equation 10).

[0022] From FIG. 13 onwards, the application scenarios of the proposed methodology are presented against the traditional approach to present the value that the same can add.

[0023] FIG. 13 illustrates the mobile HPHISO map before the start of production in the Pituba synthetic field.

[0024] FIG. 14 illustrates the drainage mesh over a mobile HPHISO map of FIG. 13.

[0025] FIG. 15 illustrates the drainage mesh on a mobile HPHISO map after 27 years of field production from the drainage mesh of FIG. 14.

[0026] FIG. 16 illustrates the drainage mesh with the addition of 3 mesh weighting wells positioned in the best locations according to the traditional HPHISO methodology—Palt1, Palt2 and Palt3.

[0027] FIG. 17 illustrates the gain in production curve and FR due to the addition of the 3 weighting wells according to the traditional methodology (FIG. 16).

[0028] FIG. 18 illustrates the original drainage mesh over the mobile HPHISO map filtered after 27 years of production.

[0029] FIG. 19 illustrates the drainage mesh with the addition of 3 mesh weighting wells positioned in the best locations according to the methodology proposed by this invention—Ptest1, Ptest2 and Ptest3.

[0030] FIG. 20 illustrates the production curve gain and FR due to the addition of the 3 weighting wells according to the proposed methodology (FIG. 19) against the original mesh and the mesh of FIG. 16.

[0031] FIG. 21 illustrates the map of mobile HPHISO filtered from the weighted mesh of FIG. 19 after 27 years of production, indicating that there are still poorly drained areas that could be future weighting opportunities.

[0032] FIG. 22 illustrates the entire flow of activities of the proposed methodology.

DETAILED DESCRIPTION OF THE INVENTION

[0033] The main inventive character of the methodology in question is to bring the concept of adding filters over the traditional HPHISO map in order to purge regions that appear to have a high volume of residual oil and/or areas that present high porosities, but low permeabilities, reducing the drainage capacity of the same. The invented filters were reservoir transmissibility filters and scan quality filters. Such filters will work by reducing the weight of regions of remaining oil identified as origin 1 and 2, highlighting more the regions of origin 3. This invention can be applied to any type of reservoir, whether carbonate or siliciclastic reservoirs.

[0034] The transmissibility filter was proposed to de-characterize the remaining oil in 2 scenarios that can arise for any field that is being worked: (i) naturally porous reservoirs, but with a low permeability level, making the drainage of this region very impaired, and (ii) regions of the reservoir that are considered porous whose permeability is penalized during the historical adjustment process.

[0035] So, the filter is applied in a simple way, defining that if the transmissibility of the cell, given by equation (3), is lower than a given limit (defined by the user), all the oil is zeroed for the filtered HPHISO map.

[0036] The filter is applied to each grid cell and, later, the same is integrated in the depth, generating an average value along the reservoir defined by the user. Likewise, the user is also responsible for defining the value of kh such that there should be no fluid flow, FIG. 1.

[0037] The scan quality filter has the purpose of penalizing those regions that present a level of oil recovery above a certain cut, to be defined by the user. This filter is necessary because the microscopic recovery factor profile against the fractional flow of water has an asymptotic profile, so that from a certain level of recovery flow FR, the marginal growth of the same in relation to the injected volumes is very small, that is, to increase the FR between the critical value and the theoretical maximum FR, it would be necessary to inject an unfeasible number of porous volumes. When reaching the critical FR, there should be a minimum economic oil saturation, but higher than the remaining oil saturation.

[0038] As a comparative basis with the Pituba field, the graph in FIG. 2 indicates the evolution of the theoretical recovery flow (FRt) with the fractional flow of water (fw). Note that, while fw is less than 80%, FRt is constant and equal to 30%. From that point on (considered fw at breakthrough), the FRt grows slowly (negative second derivative), until reaching FRt=57%, when the fw reaches 100%. This point is considered the maximum FRt. That is, from the critical point on, more and more water is produced for each unit of oil produced, which would significantly reduce the potential for increased recovery by increasing the drainage rate in that region.

[0039] It will be up to the user to define the acceptable limit value for the theoretical cut-off recovery factor (FRt). Since in post-processing it is only possible to evaluate the absolute recovery factor (absolute FR) within this cut-off value, the scan efficiency will be considered. That is, if the scan efficiency is naturally low, the observed FR will also be, so that the filter will apply high values, as can be seen in the graph of FIG. 3. This factor will depend, among others, on the type of rock, mineralogy, model scale, reservoir wettability, etc.

[0040] The scan quality filter (or filtroFR) is calculated from the cell-by-cell FR estimate (equation 6). Since there is a significant variability in the FR by region of the reservoir, due to heterogeneities, the average of the FR along the column (equation 7) is applied to, finally, apply the filter, given by equation 8. It uses the FR limit as input. Note that it may be necessary to truncate the maximum and minimum values of the FR filter to 1.0 and 0.0, respectively.

[0041] The entire logical flow of the methodology being proposed is presented below. The same is consolidated in FIG. 21.

[0042] Four basic information will be provided for the methodology: [0043] a. Transmissibility Cut (khCorte); [0044] b. Recovery Cut (frCorte); [0045] c. Top Layer of the Evaluated Reservoir (kTopo); [0046] d. Evaluated Reservoir Base Layer (kBase).

[0047] Calculate and distribute in the field the porous volume property above the COA (BLOCKPVOL above the COA):

[00001]$\begin{array}{cc}{\mathrm{BVPacimaCOA}}_{\mathrm{if}}=\{\begin{array}{cc}0,& {\mathrm{Sw}}_{\mathrm{ijk}}\left(t_i\right)>0,90\\ {\mathrm{BLOCKVOL}}_{\mathrm{ijk}},& {\mathrm{SW}}_{\mathrm{ijk}}\left(t_i\right)<0,90\end{array}& \left(1\right)\end{array}$$\begin{array}{cc}{\mathrm{BLOCKVOL}}_{\mathrm{ijk}}=H_{\mathrm{ijk}}*{\mathrm{NTG}}_{\mathrm{ijk}}*{\mathrm{PHI}}_{\mathrm{ijk}}& \left(2\right)\end{array}$

[0048] Calculate transmissibility in each model cell:

KH.sub.ijk=H.sub.ijk*√{square root over (PERMI.sub.ijk*PERMJ.sub.ijk)}  (3)

[0049] Apply transmissibility filter in each model cell:

[00002]$\begin{array}{cc}{\mathrm{filtroKHcc}}_{\mathrm{ijk}}=\{\begin{array}{cc}1,& {\mathrm{KH}}_{\mathrm{ijk}}>\mathrm{khCorte}\\ 0,& {\mathrm{KH}}_{\mathrm{ijk}}<\mathrm{khCorte}\end{array}& \left(4\right)\end{array}$

[0050] Calculate weighted average of the transmissibility filter calculated by (4) by BVPacimaCOA between kBase and kTopo (along the k):

[00003]$\begin{array}{cc}{\mathrm{filtroKH}}_{\mathrm{ij}}=\frac{{.Math.}_k{\mathrm{filtroKHcc}}_{\mathrm{ijk}}*{\mathrm{BVPacimaCOA}}_{\mathrm{ijk}}}{{.Math.}_k{\mathrm{BVPacimaCOA}}_{\mathrm{ijk}}}& \left(5\right)\end{array}$

[0051] Calculate FR cell-by-cell within the reservoir:

[00004]$\begin{array}{cc}{\mathrm{FRcc}}_{\mathrm{ijk}}=\frac{{\mathrm{So}}_{\mathrm{ijk}}\left(\mathrm{ti}\right)-{\mathrm{So}}_{\mathrm{ijk}}}{{\mathrm{So}}_{\mathrm{ijk}}\left(\mathrm{ti}\right)}& \left(6\right)\end{array}$

[0052] Calculate the weighted average of the FRcc property by BVPacimaCOA between kBase and kTopo (along the k):

[00005]$\begin{array}{cc}{\mathrm{FRmedio}}_{\mathrm{ij}}=\frac{{.Math.}_k{\mathrm{FRcc}}_{\mathrm{ijk}}*{\mathrm{BVPacimaCOA}}_{\mathrm{ijk}}}{{.Math.}_k{\mathrm{BVPacimaCOA}}_{\mathrm{ijk}}}& \left(7\right)\end{array}$

[0053] Calculate scan filter applied on the weighted property FRmedio:

[00006]$\begin{array}{cc}{\mathrm{filtroFR}}_{\mathrm{ij}}=-\left(\frac{1}{\mathrm{frCorte}}\right)*{\mathrm{FRmedio}}_{\mathrm{ij}}+1& \left(8\right)\end{array}$

[0054] Calculate the mobile HPHISO integrated along the column between kBase and kTop (along the k):

[00007]$\begin{array}{cc}{\mathrm{HPHISOm}}_{\mathrm{ij}}=\underset{k}{.Math.}{H}_{\mathrm{if}}*{\mathrm{PHI}}_{\mathrm{ij}}*{\mathrm{SOm}}_{\mathrm{ij}}& \left(9\right)\end{array}$

[0055] Product between filtroKH, filtroFR and HPHISOm, generating the Quality Map as Indicators of Poorly Drained Regions and Potential for Mesh Optimization, which is the invention itself:

HPHISOm.sub.ij.sup.filtrado=HPHISOm.sub.ij*filtroFR.sub.ij*filtroKH.sub.ij  (10)

[0056] To exemplify the presented method, a synthetic flow model, called Pituba, was used. The same was worked with a drainage mesh with 8 producers and 7 injectors, and an edge injection concept (which would tend to maximize the recovery).

[0057] The model is composed of 85,905 cells, between active and inactive, with only 23 layers, each of no more than 7 meters. The total VOIP of the model is 705.9 MM bbl (112.23 MM m.sup.3)) (776.5 MM BOE (4,751 MM kJ)), with medium viscosity oil (between 1.75 and 2.0 cp) and low solubility ratio (about 100 m.sup.3/m.sup.3). The field represented is also shallow, with original pressure at 305 kgf/cm.sup.2 (29,910 MPa).

[0058] The basal concept of the original drainage mesh was to use 8 producers distributed along the structural high and 7 injectors positioned on the flanks, providing a peripheral injection. The mesh in question is represented in FIG. 14. The productive life period of the field is assumed to be a total of 27 years.

[0059] The application of the proposed method, as presented above, would have the ability to identify only regions with good potential for production and that were not explored with the current mesh, as shown in FIG. 18.

[0060] By means of FIG. 18, it is clear that the entire region more to the NE of the field is much better drained than the region more to the SW, although this is not so clear from the mobile HPHISO map. Thus, this region would emerge as a potential region for locating other wells (complementary projects, for example) and thus maximizing the recovery of the field. An exercise was carried out to increase the mesh for 11 producers, keeping the 7 original injectors in 2 different ways: (1) by the traditional method —HPHISO map after 27 years of production, and (2) applying the proposed methodology of the filters.

[0061] The first approach is shown in FIG. 16, with additional producers positioned and named “Palt1”, “Palt2” and Palt3″ from the traditional mobile HPHISO map. This mesh will be called weighted mesh 1 and has the 3 wells in question basically positioned in the region with the greatest original volume of oil in the field.

[0062] As illustrated in FIG. 17, the variation of the field production curve between the original mesh and the weighted mesh 1 generated an increase in FR of the order of 0.5% (18.8% to 19.3%), that is, 0.17% per well, which is an insignificant result given the mobile HPHISO not yet drained by the original mesh (FIG. 16).

[0063] By observing the filtered mobile HPHISO map (proposed methodology), however, there can be observed that other regions not highlighted in FIG. 16 appear as areas of greater potential for positioning new wells. From this map (FIG. 18), 3 more producing wells were allocated, called “Ptest1”, “Ptest2” and “Ptest3”. This mesh was called Weighted Mesh 2.

[0064] Contrary to the results of the weighted mesh 1, the recovery gains in the fields were very expressive, with the 3 wells in question guaranteeing an increase of 1.8% in FR in relation to the original (18.8% to 20.6%) and 1.3% in relation to the scope of weighted mesh 1 (19.3% to 20.6%). When compared to the weighted mesh 1, the result is even more expressive because, with the same investment, a result of +0.5% FR or +1.8% FR can be obtained, depending on the use of the post-processing methodology that is used. In terms of updated production, the additional 3 wells according to the proposed methodology generated an increase of 14%, or 8 MM bbl (1,272 MM m.sup.3) updated.

[0065] The filtered HPHISO map of FIG. 19, indicating the quality of field drainage after the above-mentioned changes follows below. Note that, despite the mesh weighting, there are still regions of poorly drained oil, which would have the potential to still further increase the FR of the field. In any case, the evolution of the drainage quality of the field is noticeable, and, although the mobile HPHISO map indicates a large volume of mobile oil not yet drained, the proposed quality map indicates that the field drainage is efficient and there are few new optimization opportunities.

[0066] As indicated in the map in FIG. 20, the 3 wells are located in regions of greater mobile HPHISO, without, however, evaluating the drainage level of the regions of the deposit. See below the 3 new wells on the filtered mobile HPHISO map. The same clearly indicates that the scan in these positions was already very high, so their positioning in these locations should not bring gain in the recovery.

[0067] The perception is confirmed by means of the evaluation results of this new mesh, of 11 producers and 7 injectors, with the 3 additional wells positioned in these locations. Although the mesh weighting produces an anticipation of production, as in the Weighted Mesh 1 example, there is not the same gain in recovery as seen in the previous case. This was already expected because, as shown in the filtered map, the wells are at points where the HPHISO filtered after 27 years of production practically no longer indicates the presence of oil.

[0068] The following table summarizes the gains associated with increasing the scope with locations defined by the traditional method (mobile HPHISO) and the proposed method (filtered mobile HPHISO).

TABLE-US-00001 TABLE 1 Gains associated with different approaches. FR Updated Np Original Mesh 18.8%   58 MM bbl (MM m.sup.3) ) Weighted Mesh 1- 19.3% 62 MM bbl (+4 MM) Traditional Method (+0.5%)  (9,857 MM m.sup.3 (+0,636 MM m.sup.3) ) Weighted Mesh 2- 20.6% 66 MM bbl (+8 MM) Proposed Method (+1.8%) (10,493 MM m.sup.3 (+1,272 MM m.sup.3) )