Modeling and simulating faults in subterranean formations

Patent metadata

Field	Value
Publication number	US-12099159-B2
Application number	US-202017594264-A
Country	US
Kind code	B2
Filing date	Apr 15, 2020
Priority date	Apr 15, 2019
Publication date	Sep 24, 2024
Grant date	Sep 24, 2024

Abstract

Methods for modeling and simulating fractured subterranean volumes include a method including obtaining geological data representing a subterranean volume, generating a structural model thereof in depositional space and in structural space. The method includes selecting a first cell and a second cell in the model, the first and second cells being juxtaposed in geological space and defining a fault face where the first and second cells are intersected by a fault, identifying a first point on the fault face, and calculating slip curves. Respective slip curves originate at the point and extend across the fault in geological space to a respective second point of a plurality of second points. The second points are co-located with the first point in the depositional space. The method includes calculating fault rock properties at the first point based on the slip curves and adjusting the model to include the fault rock properties.

First claim

What is claimed is:

1. A method, comprising:
obtaining, via one or more processors, geological data representing a subterranean volume, wherein the geological data relates to seismic data, a well log, a core sample, or a combination thereof, acquired from one or more components;
generating, via the one or more processors, a structural model of the subterranean volume in a depositional space and in a geological space, wherein the structural model comprises a grid of cells;
selecting, via the one or more processors, a first cell and a second cell of the grid of cells, the first and second cells being juxtaposed in the geological space and defining a fault face where the first and second cells are intersected by a fault;
identifying, via the one or more processors, a first point in the geological space that is on the fault face of the first and second cells in the geological space, wherein the first point in the geological space is a single point and is represented as a plurality of first points in the depositional space;
identifying, via the one or more processors, a plurality of second points in the depositional space, wherein respective second points of the plurality of second points in the depositional space are co-located with respective first points of the plurality of first points in the depositional space;
transforming, via the one or more processors, the structural model from the depositional space to the geological space to identify the plurality of second points in the geological space based on the plurality of second points in the depositional space;
calculating, via the one or more processors, a plurality of slip curves in the geological space, wherein respective slip curves originate at the first point in the geological space and extend across the fault in the geological space to respective second points of the plurality of second points in the geological space;
calculating, via the one or more processors, one or more fault rock properties at the first point based on clay properties encountered along the slip curves;
adjusting, via the one or more processors, the structural model to include the one or more fault rock properties to improve accuracy of the structural model to obtain a modified structural model having an improved geological representation of a reservoir; and
controlling, via the one or more processors, sensing, drilling, injecting, extracting, or a combination thereof, based on the modified structural model having the improved geological representation of the reservoir.

2. The method of claim 1 , wherein the one or more fault rock properties comprise fault thickness, fault permeability, or a combination thereof.

3. The method of claim 1 , wherein one of the plurality of first points is on an upthrown side of the fault, and another one of the plurality of first points is on a downthrown side of the fault.

4. The method of claim 3 , wherein the respective first points of the plurality of first points are offset across the fault from the respective second points of the plurality of second points in the geological space.

5. The method of claim 1 , wherein the plurality of slip curves comprises a downthrown slip curve that extends from the first point to a downthrown side of the fault, and an upthrown slip curve that extends from the first point to an upthrown side of the fault.

6. The method of claim 1 , wherein at least one of the one or more fault rock properties are calculated based on an average of the clay properties encountered along the slip curves.

7. The method of claim 1 , wherein generating the structural model comprises:
generating the structural model, in the geological space, of the subterranean volume using the geological data;
transforming the structural model from the geological space to the depositional space;
defining the grid of cells for the structural model in the depositional space; and
transforming the structural model, including the grid, to the geological space to define a depogrid.

8. The method of claim 1 , wherein the geological space comprises two horizontal axes and a depth axis, and wherein the depositional space comprises two horizontal axes and a depositional time axis.

9. The method of claim 1 , further comprising visualizing the modified structural model of the subterranean volume using a computer display.

10. The method of claim 1 , further comprising identifying a presence of the reservoir, or simulating a production fluid flow from the reservoir, using the modified structural model having the improved geological representation of the reservoir.

11. The method of claim 1 , further comprising:
calculating one or more fault-face properties for individual boundaries between adjacent cells that lie on the fault;
calculating transmissibility in the structural model using the one or more fault-face properties for the individual boundaries; and
adjusting the structural model by adjusting one or more of the fault-face properties of one of the individual boundaries.

12. A method, comprising:
obtaining geological data representing a subterranean volume, wherein the geological data relates to seismic data, a well log, a core sample, or a combination thereof, acquired from one or more components;
generating a static model of the subterranean volume, the static model comprising a grid of cells, wherein the static model is generated based on a structural model of the subterranean volume in a depositional space and a geological space, and faults in the subterranean volume are represented as boundaries between adjacent cells that lie on the faults;
calculating one or more fault-face properties for individual boundaries between the adjacent cells that lie on one of the faults based on a plurality of slip curves calculated in the geological space of the structural model;
calculating transmissibility in the static model using the one or more fault-face properties for the individual boundaries;
adjusting the static model by adjusting one or more of the fault-face properties of one of the individual boundaries to improve accuracy of the static model to generate a modified static model having an improved geological representation of a reservoir; and
controlling sensing, drilling, injecting, extracting, or a combination thereof, based on the modified static model having the improved geological representation of the reservoir.

13. The method of claim 12 , further comprising history matching using the modified static model, wherein history matching comprises determining uncertainty of one or more geological parameters to determine a sensitivity of the modified static model to modifications to the transmissibility.

14. The method of claim 13 , wherein the one or more geological parameters comprise static permeability, fault permeability and fault thickness multipliers, or a combination thereof.

15. The method of claim 12 , wherein the one or more fault-face properties comprise fault thickness and permeability.

16. The method of claim 15 , wherein the transmissibility is calculated based in part on the fault thickness and the permeability of the fault represented at the individual boundaries between the adjacent cells.

17. The method of claim 12 , wherein generating the static model comprises:
generating the structural model of the subterranean volume in the depositional space and in structural the geological space, wherein the structural model comprises a grid of cells;
selecting a first cell and a second cell of the grid of cells, the first and second cells being juxtaposed in the geological space and defining a fault fa

Assignees

• Schlumberger Technology Corp

Inventors

• Harris Simon
• Santoshini Samita
• Grenfell Stephen

Classifications

• G01V2210/642

  Faults · CPC title

• G01V1/301

  for determining seismic cross-sections or geostructures · CPC title

• E21B2200/20

  Computer models or simulations, e.g. for reservoirs under production, drill bits · CPC title

• G01V20/00 Primary

  Geomodelling in general · CPC title