Methods and systems for modeling subsurfaces containing partial faults

Patent metadata

Field	Value
Publication number	US-11041969-B2
Application number	US-201916389041-A
Country	US
Kind code	B2
Filing date	Apr 19, 2019
Priority date	Jun 22, 2018
Publication date	Jun 22, 2021
Grant date	Jun 22, 2021

Abstract

Disclosed herein are geologic modeling methods and systems employing function-based representations of horizons intersected by partial faults. An illustrative method embodiment includes: (a) obtaining a seismic image volume; (b) identifying a horizon within the seismic image volume, said horizon being intersected by a partial fault; (c) deriving a function-based representation of the horizon, the representation being continuous except across the partial fault; (e) constructing a watertight subsurface model using the function-based representation; (f) assigning petrophysical parameter values to compartments of the watertight subsurface model; and, optionally, (g) storing or displaying the watertight subsurface model.

First claim

What is claimed is:

1. A subsurface modeling method that comprises:
obtaining a seismic image volume;
identifying a horizon within the seismic image volume, said horizon being intersected by a partial fault;
deriving a function-based representation of the horizon, the representation being continuous except across the partial fault, wherein said deriving employs a radial basis function method in which a weighted sum of radial basis functions is fitted to nodes on the horizon, subject to elimination of direct interaction between nodes on opposite sides of the partial fault and wherein said elimination is provided by (i) adjusting matrix entries in an equation to solve for weight coefficients, said adjusting producing an effect of a substantially increased distance value between nodes on opposite sides of the partial fault or (ii) adjusting matrix entries in an equation to solve for weight coefficients, said adjusting producing an effect of a substantially increased distance value between nodes on opposite sides of the partial fault;
constructing a watertight subsurface model using the function-based representation; and
assigning petrophysical parameter values to compartments of the watertight subsurface model.

2. The method of claim 1 , further comprising at least one of:
storing the watertight subsurface model;
using the stored watertight subsurface model to provide an interactive visualization of subsurface structures in a region corresponding to the seismic image volume;
simulating movement of reservoir fluids in said, region; and
providing an evaluation of one or more reservoir production strategies including well placements.

3. The method of claim 1 , wherein said elimination is provided by adjusting matrix entries in an equation to solve for weight coefficients, said adjusting producing an effect of a substantially increased distance value between nodes on opposite sides of the partial fault.

4. The method of claim 1 , wherein said elimination is provided by zeroing radial basis function contributions between nodes on opposite sides of the partial fault in a matrix equation.

5. A subsurface modeling system that comprises;
a memory having subsurface modeling software; and
one or more processors coupled to the memory to execute the subsurface modeling software, the software causing the one or more processors to perform operations including the steps of the method of claim 1 .

6. A non-transitory, machine-readable medium comprising code executable by a computer to implement the method of claim 1 .

7. A subsurface modeling method that comprises;
obtaining a seismic image volume;
identifying a horizon within the seismic image volume, said horizon being intersected by a partial fault;
deriving a function-based representation of the horizon, the representation being continuous except across the partial fault, wherein said deriving employs mesh vertices on or around the horizon having associated signed distance values to the horizon, said mesh vertices including a set of mesh vertices that, define mesh cells traversed by the partial fault, the set of mesh vertices further having extrapolated distance values that are extrapolated across the partial fault from associated distance values for vertices on an opposite side of the partial fault, and wherein said deriving applies a level set function method in which function values are interpolated from said associated and extrapolated distance values;
constructing a watertight subsurface model using the function-based representation; and
assigning petrophysical parameter values to compartments of the watertight subsurface model.

8. The method of claim 7 , wherein the associated and extrapolated distance values are signed based on which side of the horizon their vertices are on.

9. The method of claim 7 , wherein said constructing includes closing any horizon mismatch extending beyond the partial fault.

10. The method of claim 7 , further comprising at least one of:
storing the watertight subsurface model;
using the stored watertight subsurface model to provide an interactive visualization of subsurface structures in a region corresponding to the seismic image volume;
simulating movement of reservoir fluids in said region; and
providing an evaluation of one or more reservoir production strategies including well placements.

11. A subsurface modeling system that comprises:
a memory having subsurface modeling software; and
one or more processors coupled to the memory to execute the subsurface modeling software, the software causing the one or more processors to perform operations including the steps of the method of claim 7 .

12. A non-transitory, machine-readable medium comprising code executable by a computer to implement the method of claim 7 .

13. A subsurface modeling method that comprises:
obtaining a seismic image volume;
identifying a horizon within the seismic image volume, said horizon being intersected by a partial fault;
deriving a function-based representation of the horizon, the representation being continuous except across the partial fault, wherein said deriving includes:
creating one or more extensions of the partial fault to divide the horizon into two separate pieces;
determining for each of multiple mesh vertices on or around the horizon an associated distance value to the horizon; and
applying a level set function method to fit a first function to a first of the two separate pieces, a second function to a second of the two separate pieces, and a merge function to a set of mesh vertices that define mesh cells traversed by said one or more extensions;

constructing a watertight subsurface model using the function-based representation; and
assigning petrophysical parameter values to compartments of the watertight subsurface model.

14. The method of claim 13 , further comprising at least one of:
storing the watertight subsurface model;
using the stored watertight subsurface model to provide an interactive visualization of subsurface structures in a region corresponding to the seismic image volume;
simulating movement of reservoir fluids in said region; and
providing an evaluation of one or more reservoir production strategies including well placements.

15. A subsurface modeling system that comprises:
a memory having subsurface modeling software; and
one or more processors coupled to the memory to execute the subsurface modeling software, the software causing the one or more processors to preform operations including the steps of the method of claim 13 .

16. A non-transitory, machine-readable medium comprising code executable by a computer to implement the method of claim 13 .

Assignees

• Lin Baojiu
• Huang Hao
• Exxonmobil Upstream Res Co

Inventors

• Lin Baojiu
• Huang Hao

Classifications

• G01V2210/642

  Faults · CPC title

• G01V1/302 Primary

  in 3D data cubes · CPC title

• G06T17/05

  Geographic models · CPC title

• G01V1/282

  Application of seismic models, synthetic seismograms · CPC title

• G01V2210/66

  Subsurface modeling · CPC title