Fault removal in geological models

• US 10,036,829 B2
• Filed: 08/23/2013
• Issued: 07/31/2018
• Est. Priority Date: 09/28/2012
• Status: Active Grant

• Current Assignee: Exxonmobil Upstream Research Company (Exxon Mobil Corporation)
• Original Assignee: Exxonmobil Upstream Research Company (Exxon Mobil Corporation)
• Inventors: Ghayour, Kaveh, Bi, Linfeng, Wu, Xiaohui
• Primary Examiner(s): Fernandez Rivas, Omar F
• Assistant Examiner(s): Mikowski, Justin C
• Application Number:
  US14/423,659
• Publication Number:
  US 20150293260A1
• Time in Patent Office: 1,803 Days
• CPC Class Codes: G01V 1/282 Application of seismic models, synthetic seismograms
• CPC Class Codes: G01V 20/00 Geomodelling in general
• CPC Class Codes: G01V 2210/66 Subsurface modeling
• CPC Class Codes: G06F 17/10 Complex mathematical operations function generation by table look-up G06F1/03; evaluation of elementary functions by calculation G06F7/544

Abstract

Method for transforming a discontinuous, faulted subsurface reservoir into a continuous, fault-free space where a complete geological model based on selected geological concepts can be built and updated efficiently. Faults are removed in reverse chronological order (62) to generate a pseudo-physical continuous layered model, which is populated with information according to the selected geological concept (68). The fault removal is posed as an optimal control problem where unknown rigid body transformations and relative displacements on fault surfaces are found such that deformation of the bounding horizons and within the volume near the fault surface are minimized (63). A boundary-element-method discretization in an infinite domain is used, with boundary data imposed only on fault surfaces. The data populated model may then be mapped back to the original faulted domain such that a one-to-one mapping between continuous and faulted spaces may be found to a desired tolerance (72).

Claims

1. A computer-implemented method for fault removal of a plurality faults in a subsurface geological model to populate the model with desired information, comprising:

(a) ordering the plurality of faults in reverse chronological order, and selecting a fault of the plurality of faults based on the reverse chronological order;

(b) removing the selected fault by iteratively solving, using a computer, an optimal control problem wherein Laplace'"'"'s equation is solved for an optimal set of rigid body transformations and boundary displacement vectors on the fault'"'"'s surface;

(c) repeating (b) as necessary to remove any remaining faults in the plurality of faults, one at a time, in reverse chronological order, resulting in a transformation of the geological model from faulted space to continuous space;

(d) populating said geological model with faults removed with selected grids or surfaces or physical property values;

(e) computing a mapping of the populated geological model back to the faulted space;

(f) storing the mapping and the populated geological model in the memory of the computer; and

(g) causing a well to be drilled based on the mapping for the production of hydrocarbons.

2. The method of claim 1, wherein the mapping back to the faulted state is performed in chronological order, one fault at a time.

3. The method of claim 1, wherein Laplace'"'"'s equation is solved by a boundary element method discretization in an infinite domain, with boundary data imposed on fault surfaces and not on bounding horizons.

4. The method of claim 1, wherein the optimal control problem comprises minimizing an objective function that measures deformation of bounding horizons and internal deformation near a fault surface.

5. The method of claim 1, wherein said optimal set of rigid body transformations and boundary displacement vectors on the fault'"'"'s surface describe subsurface horizons before the selected fault occurred.

6. The method of claim 1, wherein parameters of the rigid body transformations appear in boundary conditions used in the solving of Laplace'"'"'s equation.

7. The method of claim 1, wherein the mapping of the populated geological model back to the faulted space comprises determining an inverse of the transformation in (c).

8. The method of claim 7, wherein the inverse transformation is determined by iteratively solving an optimal control problem.

9. The method of claim 1, wherein removing the selected fault comprises matching two sides of the fault through an auxiliary conceptual surface that is determined in the solving of the optimal control problem, wherein both sides of the fault are free to move and meet at the auxiliary conceptual surface.

11. The method of claim 1, further comprising integrating rigid body transformations to minimize distortion.

12. The method of claim 11, wherein the rigid body transformations comprise one or more linear functions of the Cartesian coordinates.

13. The method of claim 1, further comprising displaying one of the mapping, the populated geological model, and any combination thereof.

10. A method for producing hydrocarbons from a subsurface region, comprising:

obtaining a seismic survey of the subsurface region;

using processed and interpreted data from the seismic survey to develop a geological model of the subsurface region, said geological model containing one or more faults;

using a method to populate the model with geophysical information, wherein the method comprises;

(a) ordering the plurality of faults in reverse chronological order, and selecting a fault of the plurality of faults based on the reverse chronological order;

(b) removing the selected fault by iteratively solving, using a computer, an optimal control problem wherein Laplace'"'"'s equation is solved for an optimal set of rigid body transformations and boundary displacement vectors on the fault'"'"'s surface;

(c) repeating (b) as necessary to remove any remaining faults in the plurality of faults, one at a time, in reverse chronological order, resulting in a transformation of the geological model from faulted space to continuous space;

(d) populating said geological model with faults removed with selected grids or surfaces or physical property values;

(e) computing a mapping of the populated geological model back to the faulted space; and

(f) storing the mapping and the populated geological model in the memory of the computer;

using the populated model in planning production of hydrocarbons from the subsurface region; and

producing hydrocarbons from the subsurface region in accordance with the plan.