System and method for automated reservoir targeting

• US 7,096,172 B2
• Filed: 07/18/2003
• Issued: 08/22/2006
• Est. Priority Date: 01/31/2003
• Status: Active Grant

• Current Assignee: Landmark Graphics Corporation (Halliburton Company)
• Original Assignee: Landmark Graphics Corporation (Halliburton Company)
• Inventors: McColpin, Glenn Robert, Colvin, Richard Daniel
• Primary Examiner(s): Rodriguez, Paul L.
• Assistant Examiner(s): Thangavelu, Kandasamy
• Application Number:
  US10/622,976
• Publication Number:
  US 20040153298A1
• Time in Patent Office: 1,131 Days
• Field of Search: 703/2, 703/10, 702/12, 702/13, 702/2, 702/14, 345/420, 345/762, 367/41
• US Class Current: 703/10
• CPC Class Codes: E21B 43/30 Specific pattern of wells, e.g. optimising the spacing of wells
• CPC Class Codes: G01V 11/00 Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Abstract

Methods and systems for reservoir targeting, the systems implementing methods comprising filtering a plurality of values in a three-dimensional model to eliminate values which are below a threshold to create a filtered three-dimensional model; developing a first matrix from the filtered three-dimensional model representing a two-dimensional model of the reservoir, wherein the first matrix includes a plurality of cell center locations, cell areas, and the plurality of values; developing a second matrix from the first matrix using a distance-weighted sum of the plurality of values; and selecting target locations from the second matrix based on the distance-weighted sum of the plurality of values.

Claims

1. A computer-implemented method of reservoir targeting, comprising:

(a) building, from seismic data, a three dimensional model of a reservoir;

(b) triangulating the three dimensional model of the reservoir to create a schematic model comprising a three dimensional grid of cells;

(c) filtering the schematic model to eliminate cells with attribute values below a threshold and/or do not meet other predetermined selection criteria;

(d) determining a set of contiguous cells for at least one selected X and Y location in the schematic model, wherein the set of contiguous cells is determined by moving a window of a plurality of cells to select a window of a plurality of cells that has a maximum value of a sum of values of an attribute of interest, wherein the sum is taken over the cells in the selected window, the plurality of cells having the maximum value of the sum being the most desirable cells for the at least one selected X and Y location;

(e) recording a center of location of the most desirable cells along with an area of the cells and the maximum value of the sum of values of the attribute of interest for the at least one selected X and Y location in a first matrix, wherein the first matrix is a two dimensional matrix;

(f) for each cell in the first matrix, taking each cell as a selected cell, calculating a distance weighted sum of values in the first matrix of all the cells within a multiple of a spacing radius from a center point of the selected cell, wherein a weight is selected to give more weight to cells located closer to the selected cell and less weight to cells located further from the selected cell and entering the distance weighted sum in a second matrix as an accumulated value for the selected cell;

(g) selecting the cell location having a maximum accumulated value in the second matrix as a target location;

(h) setting the value of the selected cell at the target location in the second matrix and the value of all the cells in the second matrix, within a multiple of the spacing radius, to zero in the first matrix;

(i) repeating steps (f) to (h) until a specified number of target locations are identified or there are no more cells with an accumulated value greater than zero.

2. The method of claim 1, wherein the accumulated value (AccumValue) associated with each center location is derived using the relationship:

AccumValue=(CumWeightedValue/CumWeight)*CumWeightedArea,

wherein

CumWeightedValue=Σ

cellvalue*weight,

CumWeightedArea=Σ

cellarea*weight, and

CumWeight=Σ

(SpacingRadius−

DistanceFromCell)/SpacingRadius,

where DistanceFromCell is defined as the actual distance from the cell for which the accumulated value is being calculated to a cell that is being taken into consideration and SpacingRadius is a user-defined value representing a reservoir draining radius for each target location.

3. The method of claim 1, wherein selecting target locations includes determining whether there are existing targets for the reservoir, and if existing targets are identified, eliminating possible targets within a predetermined distance from the existing targets before selecting new targets.

4. The method of claim 3, wherein selecting target locations includes an iterative process of selecting the targets based on a first preferred attribute value, eliminating other targets within a predetermined distance from an initial target, and selecting a next preferred attribute value for a next target location.

5. The method of claim 1, wherein selecting target locations includes an iterative process of selecting a target based on a first preferred attribute value, eliminating other targets within a predetermined distance from an initial target, and selecting a next preferred attribute value for a next target location.

6. The method of claim 1, wherein selecting target locations further comprises ranking the selected target locations and displaying a user-selected percentage of the ranked target locations.

7. A computer-readable medium having computer-executable instructions which when executed on a computer perform a process for reservoir targeting, the process comprising:

(a) building, from seismic data, a three dimensional model of a reservoir;

(b) triangulating the three dimensional model of the reservoir to create a schematic model comprising a three dimensional grid of cells;

(c) filtering the schematic model to eliminate cells with attribute values below a threshold and/or do not meet other predetermined selection criteria;

(d) determining a set of contiguous cells for at least one selected X and Y location in the schematic model, wherein the set of contiguous cells is determined by moving a window of a plurality of cells to select a window of a plurality of cells that has a maximum value of a sum of values of an attribute of interest, wherein the sum is taken over the cells in the selected window, the plurality of cells having the maximum value of the sum being the most desirable cells for the at least one selected X and Y location;

(e) recording a center of location of the most desirable cells along with an area of the cells and the maximum value of the sum of values of the attribute of interest for the at least one selected X and Y location in a first matrix, wherein the first matrix is a two dimensional matrix;

(f) for each cell in the first matrix, taking each cell as a selected cell, calculating a distance weighted sum of values in the first matrix of all the cells within a multiple of a spacing radius from a center point of the selected cell, wherein a weight is selected to give more weight to cells located closer to the selected cell and less weight to cells located further from the selected cell and entering the distance weighted sum in a second matrix as an accumulated value for the selected cell;

(g) selecting the cell location having a maximum accumulated value in the second matrix as a target location;

(h) setting the value of the selected cell at the target location in the second matrix and the value of all the cells in the second matrix, within a multiple of the spacing radius, to zero in the first matrix;

(i) repeating steps (f) to (h) until a specified number of target locations are identified or there are no more cells with an accumulated value greater than zero.

8. The computer-readable medium of claim 7, wherein developing the second accumulated value matrix further comprises deriving a accumulated value (AccumValue) associated with each center location using the relationship:

AccumValue=(CumWeightedValue/CumWeight)*CumWeightedArea,

wherein

CumWeightedValue=Σ

cellvalue*weight,

CumWeightedArea=Σ

cellarea*weight, and

CumWeight=Σ

(SpacingRadius−

DistanceFromCell)/SpacingRadius,

where DistanceFromCell is defined as the actual distance from the cell for which the accumulated value is being calculated to a cell that is being taken into consideration and SpacingRadius is a user-defined value representing a reservoir draining radius for each target location.

9. The computer-readable medium of claim 7, wherein selecting target locations includes determining whether there are existing targets for the reservoir, and if existing targets are identified, eliminating possible targets within a predetermined distance from the existing targets before selecting new targets.

10. The computer-readable medium of claim 9, wherein selecting target locations includes an iterative process of selecting the targets based on a first preferred attribute value, eliminating other targets within a predetermined distance from an initial target, and selecting a next preferred attribute value for a next target location.

11. The computer-readable medium of claim 7, wherein selecting target locations includes an iterative process of selecting a target based on a preferred attribute value, eliminating other targets within a predetermined distance from a initial target, and selecting a next preferred attribute value for a next target location.

12. The computer-readable medium of claim 7, wherein selecting target locations further comprises ranking the selected target locations and displaying a user-selected percentage of the ranked target locations.