Hydrocarbon migration and accumulation methods and systems
Methods for simulating hydrocarbon (HC) migration and accumulation in a subsurface formation are provided. The methods include determining a plurality of HC mass associated with a plurality of grid cells representing the subsurface formation. The methods also include determining a plurality of HC mass outflow magnitudes for one or more grid cells in the plurality of grid cells, the HC mass outflow magnitude for each of the one or more grid cells having an upper bound value based on the HC mass in that grid cell. The methods update the HC mass of the plurality of grid cells based on the plurality of HC mass outflow magnitudes. The methods also determine that a set of grid cells in the plurality of grid cells contain an excess mass of HC, and perform an accumulation process to model filling of a trap associated with the set of grid cells.
1. A computer-implemented method for simulating hydrocarbon (HC) migration and accumulation in a subsurface formation, comprising:
determining, with one or more hardware processors, a plurality of HC mass associated with a plurality of grid cells representing the subsurface formation;
determining, with the one or more hardware processors, a plurality of HC mass outflow magnitudes for one or more grid cells in the plurality of grid cells, the HC mass outflow magnitude for each of the one or more grid cells having an upper bound value based on the HC mass in that grid cell;
updating, with the one or more hardware processors, the HC mass of the plurality of grid cells based on the plurality of HC mass outflow magnitudes;
determining, with the one or more hardware processors, that a set of grid cells in the plurality of grid cells contain an excess mass of HC;
in response to determining that a set of grid cells contain an excess mass of HC, performing, with the one or more hardware processors, an accumulation process to model filling of a trap associated with the set of grid cells;
generating, with the one or more hardware processors, accumulation output data based on the accumulation process; and
generating, with the one or more hardware processors, at least one flow control equipment command based on the accumulation output data.
2. The computer-implemented method of claim 1 , wherein at least one grid cell of the trap is being filled, the at least one grid cell having a hydrocarbon potential that is less than hydrocarbon potential of neighboring grid cells adjacent to the at least one grid cell, the neighboring grid cells being outside the trap.
3. The computer-implemented method of claim 1 , wherein the HC mass of the plurality of grid cells is being updated based on the plurality of HC mass outflow magnitude determined for a time step, the time step being determined based on net efflux and volume of the one or more grid cells.
4. The computer-implemented method of claim 1 , wherein the magnitude of HC mass outflow of a grid cell is determined based on harmonic averaging of transmissibility of the grid cell and neighboring grid cells adjacent to the grid cell.
5. The computer-implemented method of claim 1 , wherein the HC mass of the plurality of grid cells is being updated based on the plurality of HC mass outflow magnitude determined for a time step, the time step being determined based on properties of one or more chemical components of the HC.
6. The computer-implemented method of claim 1 , wherein all grid cells in the trap have the same hydrocarbon potential.
7. The computer-implemented method of claim 1 , wherein the hydrocarbon potential is calculated based on a gravity-capillary equilibrium in the trap.
8. A non-transitory, computer-readable medium storing one or more instructions of a computer-implemented method for simulating hydrocarbon (HC) migration and accumulation in a subsurface formation, the one or more instructions executable by a computer system to perform operations comprising:
determining a plurality of HC mass associated with a plurality of grid cells representing the subsurface formation;
determining a plurality of HC mass outflow magnitudes for one or more grid cells in the plurality of grid cells, the HC mass outflow magnitude for each of the one or more grid cells having an upper bound value based on the HC mass in that grid cell;
updating the HC mass of the plurality of grid cells based on the plurality of HC mass outflow magnitudes;
determining that a set of grid cells in the plurality of grid cells contain an excess mass of HC;
in response to determining that a set of grid cells contain an excess mass of HC, performing an accumulation process to model filling of a trap associated with the set of grid cells;
generating accumulation output data based on the accumulation process; and
generating at least one flow control equipment command based on the accumulation output data.
9. The non-transitory, computer-readable medium of claim 8 , wherein at least one grid cell of the trap is being filled, the at least one grid cell having a hydrocarbon potential that is less than hydrocarbon potential of neighboring grid cells adjacent to the at least one grid cell, the neighboring grid cells being outside the trap.
10. The non-transitory, computer-readable medium of claim 8 , wherein the HC mass of the plurality of grid cells is being updated based on the plurality of HC mass outflow magnitude determined for a time step, the time step being determined based on net efflux and volume of the one or more grid cells.
11. The non-transitory, computer-readable medium of claim 8 , wherein the magnitude of HC mass outflow of a grid cell is determined based on harmonic averaging of transmissibility of the grid cell and neighboring grid cells adjacent to the grid cell.
12. The non-transitory, computer-readable medium of claim 8 , wherein the HC mass of the plurality of grid cells is being updated based on the plurality of HC mass outflow magnitude determined for a time step, the time step being determined based on properties of one or more chemical components of the HC.
13. The non-transitory, computer-readable medium of claim 8 , wherein all grid cells in the trap have the same hydrocarbon potential.
14. The non-transitory, computer-readable medium of claim 8 , wherein the hydrocarbon potential is calculated based on a gravity-capillary equilibrium in the trap.
15. A computer-implemented, parallel-processing system for simulating hydrocarbon (HC) migration and accumulation in a subsurface formation, the system comprising:
a computer memory; and
a hardware processor interoperably coupled with the computer memory and configured to perform operations comprising:
determining a plurality of HC mass associated with a plurality of grid cells representing the subsurface formation;
determining a plurality of HC mass outflow magnitudes for one or more grid cells in the plurality of grid cells, the HC mass outflow magnitude for each of the one or more grid cells having an upper bound value based on the HC mass in that grid cell;
updating the HC mass of the plurality of grid cells based on the plurality of HC mass outflow magnitudes;
determining that a set of grid cells in the plurality of grid cells contain an excess mass of HC;
in response to determining that a set of grid cells contain an excess mass of HC, performing an accumulation process to model filling of a trap associated with the set of grid cells:
generating accumulation output data based on the accumulation process; and
generating at least one flow control equipment command based on the accumulation output data.
16. The computer-implemented, parallel-processing system of claim 15 , wherein at least one grid cell of the trap is being filled, the at least one grid cell having a hydrocarbon potential that is less than hydrocarbon potential of neighboring grid cells adjacent to the at least one grid cell, the neighboring grid cells being outside the trap.
17. The computer-implemented, parallel-processing system of claim 15 , wherein the HC mass of the plurality of grid cells is being updated based on the plurality of HC mass outflow magnitude determined for a time step, the time step being determined based on net efflux and volume of the one or more grid cells.
18. The computer-implemented, parallel-processing system of claim 15 , wherein the magnitude of HC mass outflow of a grid cell is determined based on harmonic averaging of transmissibility of the grid cell and neighboring grid cells adjacent to the grid cell.
19. The computer-implemented, parallel-processing system of claim 15 , wherein the HC mass of the plurality of grid cells is being updated based on the plurality of HC mass outflow magnitude determined for a time step, the time step being determined based on properties of one or more chemical components of the HC.
20. The computer-implemented, parallel-processing system of claim 15 , wherein all grid cells in the trap have the same hydrocarbon potential calculated based on a gravity-capillary equilibrium in the trap.