Abstract: An unstructured grid model with actual well trajectory of individual multilateral wells of a subsurface reservoir is formed. Well trajectory data obtained during drilling of the wells and corresponding to well trajectory data stored as a structured grid model is provided as an input data set for unstructured grid simulation. The unstructured grid model may be formed in a computerized mainframe processor system, or by parallel reservoir simulation by processor nodes of a multicore processor of parallel processor nodes synchronized and under control of a master node.

Abstract: Orthogonal unstructured grids are automatically constructed for a field or reservoir model with two types of internal boundaries: complex wells and faults, or other discontinuities. The methodology is used to constructed simulation grids for reservoirs or fields which contains both complex fault planes and multi-lateral wells. A hierarchical grid point generation, prioritization, conflict point removal system is provided enabling the use of unconstrained Delaunay triangulation. High-quality orthogonal unstructured grids are produced with good convergence properties for reservoir simulation.

Abstract: Systems, methods, and computer-readable media are provided for a near-well unstructured grid model builder for generating a full-field unstructured grid for reservoir simulation. As described further below, the near-well unstructured grid model builder may include a workflow interface and a parallel unstructured grid model builder. The inputs to the near-well unstructured grid model builder may include existing well trajectory and completion data, future well data, a geological model, a structured grid simulation model, or any combination thereof. The near-well unstructured grid model builder may output a near-well unstructured grid having a specified grid resolution in regions of interest that include a well.

Abstract: A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.

Abstract: A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.

Abstract: A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.

Abstract: A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.

Abstract: A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.

Abstract: An unstructured grid model with actual well trajectory of individual multilateral wells of a subsurface reservoir is formed. Well trajectory data obtained during drilling of the wells and corresponding to well trajectory data stored as a structured grid model is provided as an input data set for unstructured grid simulation. The unstructured grid model may be formed in a computerized mainframe processor system, or by parallel reservoir simulation by processor nodes of a multicore processor of parallel processor nodes synchronized and under control of a master node.

Abstract: Orthogonal unstructured grids are automatically constructed for a field or reservoir model with two types of internal boundaries: complex wells and faults, or other discontinuities. The methodology is used to constructed simulation grids for reservoirs or fields which contains both complex fault planes and multi-lateral wells. A hierarchical grid point generation, prioritization, conflict point removal system is provided enabling the use of unconstrained Delaunay triangulation. High-quality orthogonal unstructured grids are produced with good convergence properties for reservoir simulation.

Abstract: Computer processing time and results are improved in fully-coupled fully-implicit well-reservoir simulation system using Jacobian matrix methodology. Approximate inverse preconditioners are provided which treat a well influence matrix at comparable accuracy and robustness to those for the grid-to-grid flow terms of system matrix. The methodology is highly parallelizable and the data processing can be performed faster, as fewer solver iterations are required to converge to the same acceptable tolerances.

Abstract: Computer processing time and results are improved in fully-coupled fully-implicit well-reservoir simulation system using Jacobian matrix methodology. Approximate inverse preconditioners are provided which treat a well influence matrix at comparable accuracy and robustness to those for the grid-to-grid flow terms of system matrix. The methodology is highly parallelizable and the data processing can be performed faster, as fewer solver iterations are required to converge to the same acceptable tolerances.

Abstract: Orthogonal unstructured grids are automatically constructed for a field or reservoir model with two types of internal boundaries: complex wells and faults, or other discontinuities. The methodology is used to constructed simulation grids for reservoirs or fields which contains both complex fault planes and multi-lateral wells. A hierarchical grid point generation, prioritization, conflict point removal system is provided enabling the use of unconstrained Delaunay triangulation. High-quality orthogonal unstructured grids are produced with good convergence properties for reservoir simulation.

Abstract: Unstructured grids are automatically constructed in a domain containing complex internal boundaries. Simulation grids are constructed for reservoirs or fields which contain complex fault planes. Reconciling among generated fault grid-points and other reservoir/field grid-points is performed, enabling the use of unconstrained Delaunay triangulation. High-quality orthogonal unstructured grids are provided with good convergence properties for reservoir simulation.

Abstract: 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.

Abstract: A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.

Abstract: A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.

Abstract: A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.

Abstract: A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.

Abstract: A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.