Semi-elimination methodology for simulating high flow features in a reservoir
View Patent ↗A semi-elimination methodology for simulating high flow features in a reservoir and wells is disclosed. The reservoir and wells may be divided into a plurality of cells, including small cells in wells and the reservoir and bulk cells in the bulk of the reservoir, where the small cells are smaller (e.g., by pore volume) than the bulk cells. Processing of all of the cells, including all of the small cells, may be too computationally expensive, particularly when processing is iterative. In that regard, at least some of the small cells are partly processed in an iteration, such as for flow rates, compositions, or flow derivatives. After which, some or all of the small cells are eliminated from further processing in the iteration. In that way, high flow features in a reservoir and wells may be simulated effectively.
1. A method for producing hydrocarbons from a reservoir in a subsurface and wells in the reservoir, the method comprising:
producing a simulation model of the subsurface and wells in which the subsurface and wells are divided into a first set of cells and a second set of cells, at least one aspect of the first set of cells being different from the second set of cells;
obtaining data in the first set of cells from physical sensors;
obtaining data in the second set of cells from the physical sensors;
generating reservoir simulation data based on the data in the first set of cells, wherein generating the reservoir simulation data comprises processing, in an iteration of the reservoir simulation, at least some cells in the first set of cells for at least one of flow rates, compositions or flow derivatives;
eliminating some or all of the cells in the first set of cells from further processing in the iteration of the reservoir simulation, wherein the some or all of the cells in the first set of cells are eliminated from further processing based on a size of the some or all of the cells in the first set of cells;
processing, in the iteration of the reservoir simulation, cells in the second set of cells using the generated reservoir simulation data to create updated reservoir simulation data; and
extracting hydrocarbons using the updated reservoir simulation data to predict flow patterns in the reservoir.
2. The method of claim 1 , wherein the at least one aspect comprises size of cells.
3. The method of claim 2 , wherein the subsurface and wells are divided into small cells and bulk cells, the small cells being smaller than the bulk cells in pore volume;
wherein the first set of cells comprises the small cells; and
wherein the second set of cells comprises bulk cells in a bulk of the reservoir.
4. The method of claim 3 , wherein fewer than all of the small cells in the first set of cells are eliminated from further processing in the iteration of the reservoir simulation; and
wherein a remaining subset of the small cells in the first set of cells is processed, along with the bulk cells in the second set of cells, in the iteration of the reservoir simulation using the reservoir simulation data.
5. The method of claim 4 , wherein the reservoir simulation comprises a plurality of iterations; and
wherein for each respective iteration of the plurality of iterations, performing:
processing, in the respective iteration of the reservoir simulation, the at least some cells in the first set of cells for the at least one of flow rates, compositions, or flow derivatives in order to generate reservoir simulation data;
eliminating the at least some cells in the first set of cells from further processing in the respective iteration of the reservoir simulation;
processing, in the respective iteration of the reservoir simulation, cells in the second set of cells and the remaining subset of the small cells in the first set of cells using the reservoir simulation data; and
determining whether to continue iterating or terminate the reservoir simulation.
6. The method of claim 5 , further comprising selecting a subset from the set of small cells to process for the at least one of flow rates, compositions, or flow derivatives in order to generate reservoir simulation data; and
wherein a remainder of the small cells outside of the subset are not processed for the at least one of flow rates, compositions or flow derivatives in order to generate reservoir simulation data.
7. The method of claim 6 , wherein selecting the subset to process is based on location of the small cells relative to a boundary or partition for the reservoir simulation.
8. The method of claim 7 , wherein selecting the subset to process is based on whether a respective small cell is removed from a point for enforcing well control.
9. The method of claim 8 , further comprising processing, in the respective iteration of the reservoir simulation, the bulk cells and small cells outside the subset for elimination that are adjacent to at least one small cell in the subset for the at least one of flow rates, compositions or flow derivatives in order to generate reservoir simulation data.
10. The method of claim 9 , wherein eliminating some or all of the cells in the set of small cells from further processing in the respective iteration of the reservoir simulation comprises eliminating all of the small cells from further processing.
11. The method of claim 9 , wherein eliminating some or all of the cells in the set of small cells from further processing in the respective iteration of the reservoir simulation comprises eliminating some, but not all, of the small cells from further processing.
12. The method of claim 5 , wherein determining whether to continue iterating or terminate the reservoir simulation comprising determining whether to exit a Newton loop or to terminate the reservoir simulation.
13. The method of claim 5 , further comprising constructing a pressure basis;
wherein for each of the plurality of iterations, computation of one or more phase flow and flow derivative for at least some of the small cells is performed; and
wherein elimination of the at least some of the small cells in a linear system at each of the plurality of iterations.
14. The method of claim 13 , wherein for each of the plurality of iterations, multiphase flow and flow derivative computation for at least some of the small cells is performed.
15. The method of claim 14 , wherein determining multiphase flow and derivative computation for the at least some of the small cells comprises:
computing multiphase flow and flow derivatives between the bulk cells in the bulk and neighboring small cells which are eliminated;
computing multiphase flow and flow derivatives between remaining small cells and neighboring small cells which are eliminated;
assembling a linear system by eliminating the some or all of the small cells; and
solving the linear system.
16. The method of claim 15 , wherein determining multiphase flow and derivative computation for the at least some of the small cells comprises:
configuring a system of single phase steady state equations;
computing a matrix for expressing at least one of pressure or potential in the small cells as linear combinations of pressures or potentials in surrounding bulk cells in a bulk of the reservoir and remaining small cells; and
creating computational non-neighboring connections.
17. The method of claim 16 , further comprising:
computing fluid properties, capillary pressure, relative permeability and flow derivatives for cells remaining after elimination; and
computing the multiphase flow and the flow derivatives between neighboring cells in a remaining set of bulk and small cells.
18. The method of claim 17 , wherein, for a single phase, phase mass flow rates and flow derivatives for each connection between a respective cell in the bulk or in the remaining small cells and a respective eliminated small cell are computed.
19. The method of claim 17 , wherein, for multiple phases, total mass flow rates and flow derivatives for each connection between a respective cell in the bulk or in the remaining small cells and a respective eliminated small cell are computed.
20. The method of claim 19 , wherein, for the multiple phases:
pressure values in the at least some of the small cells are computed;
the at least some of the small cells are rearranged based on the pressure values; and
computing phase mass rates and derivatives for connections attached to the at least some of the small cells sequentially using upstream weighting for inflow from a respective remaining cell to a respective neighboring small cell which is eliminated, and partitioning for flow leaving the respective neighboring small cell.
21. The method of claim 1 , further comprising managing hydrocarbons in the subsurface based on the reservoir simulation.
22. The method of claim 1 , further comprising:
generating at least one of result data, tables, spreadsheets, or an image based on the reservoir simulation; and
outputting the at least one of the result data, the tables, the spreadsheets, or the image on a display.
23. The method of claim 1 , wherein the some or all of the cells in the first set of cells are small cells.