Enhanced two point flux approximation scheme for reservoir simulation

Patent metadata

Field	Value
Publication number	US-11536873-B2
Application number	US-201716328270-A
Country	US
Kind code	B2
Filing date	Sep 28, 2017
Priority date	Sep 28, 2016
Publication date	Dec 27, 2022
Grant date	Dec 27, 2022

Abstract

A method for performing a modified two point flux approximation scheme is disclosed. The method includes: obtaining a first pressure value for a first neighbor cell and a second pressure value for a second neighbor cell, where the first neighbor cell has a first value of a reservoir property and the second neighbor cell as a second value of the reservoir property; determining a first weight using the first pressure value and a second weight using the second pressure value; calculating a third value of the reservoir property as a weighted average of the first value and the second value; and applying the third value to the first neighbor cell.

First claim

What is claimed is:

1. A method comprising:
obtaining a plurality of field measurements;
generating a digital reservoir model based on the plurality of field measurements and a gridding operation that assigns the plurality of field measurements to a plurality of cells defined for the digital reservoir model;
executing, using the digital reservoir model, a plurality of simulations to generate a plurality of simulation results, wherein each of the plurality of simulations comprises a corresponding set of inputs that are different amongst the plurality of simulations, and wherein, for each simulation in the plurality of simulations, executing comprises:
obtaining a first pressure value for a first neighbor cell of the plurality of cells and a second pressure value for a second neighbor cell of the plurality of cells, wherein the first neighbor cell has a first value of a reservoir property of the digital reservoir model and the second neighbor cell has a second value of the reservoir property;
determining a first weight using the first pressure value and a second weight using the second pressure value;
testing whether a difference between the first pressure value and the second pressure value satisfies a threshold,
generating a simulation result of the plurality of simulation results by:
responsive to the threshold being satisfied, calculating a third value of the reservoir property as a weighted average of the first value of the reservoir property and the second value of the reservoir property, and applying the third value of the reservoir property to the first neighbor cell; or
responsive to the threshold failing to be satisfied, calculating a fourth value of the reservoir property using a two point flux approximation scheme, and applying the fourth value to the first neighbor cell; and

selecting, from the corresponding set of inputs of the plurality of simulations, a selected set of inputs that correspond to a particular simulation result of the plurality of simulation results, wherein the particular simulation result matches one or more criterion, and wherein selecting is performed by an analyzer that takes, as input, at least the plurality of simulation results, and generates, as output, the selected set of inputs.

2. The method of claim 1 , further comprising:
applying the third value of the reservoir property to the second neighbor cell.

3. The method of claim 1 , further comprising:
testing whether a flow direction between the first neighbor cell and the second neighbor cell has changed between at least two iterations,
wherein the calculating the third value of the reservoir property as a weighted average is performed based on the flow direction changing.

4. The method of claim 1 , wherein determining the first weight using the first pressure value and the second weight using the second pressure value comprises calculating the first weight and the second weight as a function of the first pressure value and the second pressure value.

5. The method of claim 1 , wherein the set of inputs comprise one of a property of an injection fluid and a location of a wellbore.

6. The method of claim 1 , wherein the one or more criterion comprise at least one of a highest production, a lowest risk, a highest revenue expectation, a production within a second threshold, a cost within a third threshold, and a risk within a fourth threshold.

7. The method of claim 1 , further comprising: performing a field operation based on the selected set of inputs.

8. The method of claim 1 , wherein the first weight and the second weight comprise phase pressures between the first neighbor cell and the second neighbor cell.

9. A non-transitory computer readable medium comprising computer readable program code which, when executed on a computer, performs a computer-implemented method comprising:
obtaining a plurality of field measurements;
generating a digital reservoir model based on the plurality of field measurements and a gridding operation that assigns the plurality of field measurements to a plurality of cells defined for the digital reservoir model;
executing, using the digital reservoir model, a plurality of simulations to generate a plurality of simulation results, wherein each of the plurality of simulations comprises a corresponding set of inputs that are different amongst the plurality of simulations, and wherein, for each simulation in the plurality of simulations, executing comprises:
obtaining a first pressure value for a first neighbor cell of the plurality of cells and a second pressure value for a second neighbor cell of the plurality of cells, wherein the first neighbor cell has a first value of a reservoir property of the digital reservoir model and the second neighbor cell has a second value of the reservoir property;
determining a first weight using the first pressure value and a second weight using the second pressure value;
testing whether a difference between the first pressure value and the second pressure value satisfies a threshold,
generating a simulation result of the plurality of simulation results by:
responsive to the threshold being satisfied, calculating a third value of the reservoir property as a weighted average of the first value of the reservoir property and the second value of the reservoir property, and applying the third value of the reservoir property to the first neighbor cell; or
responsive to the threshold failing to be satisfied, calculating a fourth value of the reservoir property using a two point flux approximation scheme, and applying the fourth value to the first neighbor cell; and

selecting, from the corresponding set of inputs of the plurality of simulations, a selected set of inputs that correspond to a particular simulation result of the plurality of simulation results, wherein the particular simulation result matches one or more criterion, and wherein selecting is performed by an analyzer that takes, as input, at least the plurality of simulation results, and generates, as output, the selected set of inputs.

10. The non-transitory computer readable medium of claim 9 , further comprising computer readable program code for:
applying the third value of the reservoir property to the second neighbor cell.

11. The non-transitory computer readable medium of claim 9 , further comprising computer readable program code for:
testing whether a flow direction between the first neighbor cell and the second neighbor cell has changed between at least two iterations,
wherein the calculating the third value of the reservoir property as a weighted average is performed based on the flow direction changing.

12. The non-transitory computer readable medium of claim 9 , wherein determining the first weight using the first pressure value and the second weight using the second pressure value comprises calculating the first weight and the second weight as a function of the first pressure value and the second pressure value.

13. The non-transitory computer readable medium of claim 9 , wherein:
the set of inputs comprise one of a property of an injection fluid and a location of a wellbore; and
the one or more criterion comprise at least one of a highest production, a lowest risk, a highest revenue expectation, a production within a second threshold, a cost within a third threshold, and a risk within a fourth threshold.

14. The non-transitory computer readable medium of claim 9 , further comprising computer readable program code for:
performing a field operation based on the selected set of inputs.

15. The non-transitory computer readable medium of claim 9 , wherein the first weig

Assignees

• Schlumberger Technology Corp

Inventors

• Jonsthovel Tom
• Stone Terry Wayne

Classifications

• G01V3/18

  specially adapted for well-logging · CPC title

• G01V11/00

  Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00 · CPC title

• G06F30/20

  Design optimisation, verification or simulation (optimisation, verification or simulation of circuit designs G06F30/30) · CPC title

• G06F17/13

  Differential equations (using digital differential analysers G06F7/64) · CPC title

• G01V99/005 Primary

  Physics · mapped topic