Classifications

• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
  • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
  • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
  • G01V1/30—Analysis
  • G01V1/301—Analysis for determining seismic cross-sections or geostructures
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
  • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
  • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
  • G01V1/30—Analysis
  • G01V1/301—Analysis for determining seismic cross-sections or geostructures
  • G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
  • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
  • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
  • G01V1/282—Application of seismic models, synthetic seismograms
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
  • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
  • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
  • G01V1/30—Analysis
  • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
  • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
  • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
  • G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
  • G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
  • G01V20/00—Geomodelling in general
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
  • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
  • G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
  • G01V3/34—Transmitting data to recording or processing apparatus; Recording data
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06F—ELECTRIC DIGITAL DATA PROCESSING
  • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
  • G06F17/10—Complex mathematical operations
  • G06F17/15—Correlation function computation including computation of convolution operations
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06F—ELECTRIC DIGITAL DATA PROCESSING
  • G06F30/00—Computer-aided design [CAD]
  • G06F30/20—Design optimisation, verification or simulation
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
  • G01V2200/00—Details of seismic or acoustic prospecting or detecting in general
  • G01V2200/10—Miscellaneous details
  • G01V2200/14—Quality control
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
  • G01V2210/00—Details of seismic processing or analysis
  • G01V2210/60—Analysis
  • G01V2210/61—Analysis by combining or comparing a seismic data set with other data
  • G01V2210/616—Data from specific type of measurement
  • G01V2210/6169—Data from specific type of measurement using well-logging
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
  • G01V2210/00—Details of seismic processing or analysis
  • G01V2210/60—Analysis
  • G01V2210/66—Subsurface modeling
  • G01V2210/665—Subsurface modeling using geostatistical modeling
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
  • G01V2210/00—Details of seismic processing or analysis
  • G01V2210/60—Analysis
  • G01V2210/66—Subsurface modeling
  • G01V2210/667—Determining confidence or uncertainty in parameters

Definitions

• the disclosed embodiments relate generally to techniques for managing hydrocarbon production from subsurface reservoirs and, in particular, to a method of analyzing seismic inversion property volumes to determine inversion quality and build a reservoir property model to be used to optimize hydrocarbon production.
• Seismic exploration involves surveying subterranean geological media for hydrocarbon deposits.
• a survey typically involves deploying seismic sources and seismic sensors at predetermined locations.
• the sources generate seismic waves, which propagate into the geological medium creating pressure changes and vibrations.
• Variations in physical properties of the geological medium give rise to changes in certain properties of the seismic waves, such as their direction of propagation and other properties.
• seismic waves Portions of the seismic waves reach the seismic sensors.
• Some seismic sensors are sensitive to pressure changes (e.g., hydrophones), others to particle motion (e.g., geophones), and industrial surveys may deploy one type of sensor or both.
• the sensors In response to the detected seismic waves, the sensors generate corresponding electrical signals, known as traces, and record them in storage media as seismic data.
• Seismic data will include a plurality of“shots” (individual instances of the seismic source being activated), each of which are associated with a plurality of traces recorded at the plurality of sensors. Seismic data can be inverted to determine reservoir properties such as impedance.
• Decisions include, but are not limited to, budgetary planning, obtaining mineral and lease rights, signing well commitments, permitting rig locations, designing well paths and drilling strategy, preventing subsurface integrity issues by planning proper casing and cementation strategies, and selecting and purchasing appropriate completion and production equipment.
• the seismic inversion quality property can be used to create a region in which geophysical constraints can be added to build a more deterministic reservoir property model and a region in which more statistically stochastic models must be built. Identification of the deterministic vs. stochastic regions is a key to optimizing reservoir management.
• a method of generating reservoir property models based on seismic inversion quality includes receiving a seismic inversion product volume and a seismic attribute volume representative of a subsurface volume of interest wherein the seismic inversion product volume and the seismic attribute volume are structured as a regular grid of cells; receiving, at the computer processor, well data from wells drilled in the subsurface volume of interest; identifying collocated cells in the seismic inversion product volume by comparing the seismic inversion product volume to the well data to match an equivalent seismic inversion property derived from the well data passing through that cell; creating attribute vectors from the seismic inversion product volume and the seismic attribute volume in each of the collocated cells and a range of neighboring cells; calculating a seismic inversion error magnitude property at the collocated cells from the difference between the seismic inversion product and the equivalent seismic inversion property derived from the well data; training a data analytics method with a preprocessor to predict the observed seismic inversion error magnitude property from the attribute vectors at all cells in the grid; verifying that the data analytics
• some embodiments provide a non-transitory computer readable storage medium storing one or more programs.
• the one or more programs comprise instructions, which when executed by a computer system with one or more processors and memory, cause the computer system to perform any of the methods provided herein.
• some embodiments provide a computer system.
• the computer system includes one or more processors, memory, and one or more programs.
• the one or more programs are stored in memory and configured to be executed by the one or more processors.
• the one or more programs include an operating system and instructions that when executed by the one or more processors cause the computer system to perform any of the methods provided herein.
• Figure 1 illustrates a flowchart of a method of building reservoir models, in accordance with some embodiments
• Figure 2 is an example of steps of a method of building reservoir models, in accordance with some embodiments.
• Figure 3 is another example of steps of a method of building reservoir models, in accordance with some embodiments.
• Figure 4 is a block diagram illustrating a reservoir modeling system, in accordance with some embodiments.
• Described below are methods, systems, and computer readable storage media that provide a manner of evaluating seismic inversion results in order to build reservoir property models. These embodiments are designed to be of particular use for calculating seismic inversion error with respect to sonic well logs in order to determine seismic inversion to build reservoir property models.
• the present invention includes embodiments of a method and system for generating reservoir property models based on seismic inversion products (i.e., impedance).
• seismic inversion products i.e., impedance
• FIG. 1 illustrates a flowchart of a method 100 for generating a reservoir property model.
• a seismic inversion product volume and a seismic attribute volume representative of a subsurface volume of interest are received.
• the seismic inversion product volume and the seismic attribute volume are structured as a regular grid of cells.
• the well data may include, for example, well logs such as sonic logs.
• the method determines attribute vectors from the seismic inversion product volume related to the well data. This may involve, by way of example and not limitation, identifying collocated cells in the seismic inversion product volume by comparing the seismic inversion product volume to the well data to match an equivalent seismic inversion property derived from the well data passing through that cell and creating attribute vectors from the seismic inversion product volume and the seismic attribute volume in each of the collocated cells and its neighboring cells.
• the equivalent seismic inversion property derived from the well data is the same property as the seismic inversion product, such as acoustic impedance (acoustic_velocity* density) which may be calculated from the sonic log and density log recorded in the well.
• the average collocated seismic values are found using the well-data vicinity radius.
• attribute vector is constructed from the collocated seismic values at the location and each of its neighbor locations (as defined by the U,V,W window size). If an additional optional property is chosen, the attribute vector has twice as many components, if three properties are chosen this gives three times as many components, etc. In this way, more than one seismic attribute and seismic inversion product can be included in the method.
• data analytics are used to determine the inversion quality based on the attribute vectors. This may include calculating a seismic inversion error magnitude property at the collocated cells from the difference between the seismic inversion product and the equivalent seismic inversion property derived from the well data; training a data analytics method such as k-nearest neighbors with a preprocessor such as principal component analysis (PCA) to predict the observed seismic inversion error magnitude property from the attribute vectors at all cells in the grid.
• PCA principal component analysis
• the PCA retains some principal components based on variance which is accomplished by partitioning the retained principal components into clusters based on characteristics in principal component space; identifying, for each cluster, data points that have known property values based on the location and the well data; and interpolating data points using an inverse distance method.
• the method may verify that the data analytics method accurately predicts the seismic inversion error magnitude using cross- validation such as jack-knife cross-validation. Then the method may generate an inversion quality volume using the prediction of seismic inversion error magnitude to assess inversion quality in each cell wherein a predicted high seismic inversion error magnitude indicates a low inversion quality for that cell.
• the inversion quality is used to build the reservoir property model.
• the seismic inversion product may be used to populate the reservoir property model.
• the reservoir property model is populated based on stochastic methods that do not rely on seismic data. Those of skill in the art are aware of many stochastic methods that may be used that do not require seismic data, seismic inversion products, or seismic attributes.
• the reservoir property model is used to optimize hydrocarbon production. This may be done, for example, by selecting locations to drill wells and/or planning and performing enhanced recovery techniques such as injection for water flooding, steam flooding, and the like.
• Figures 2 and 3 show intermediate results and a map of inversion quality as produced during the steps of method 100.
• method 100 is applied to generate the inversion quality volume of which a map-view slice is shown at the bottom.
• the inversion quality shows areal trends of around .8 quality, indicating that much of the seismic inversion product can be used for the reservoir model.
• the jack-knife cross-validation in the upper left shows a good linear correlation and the cumulative distribution function (CDF) plot of the inversion quality shown in the upper right shows the quality is high.
• Figure 3 shows the output of method 100 wherein much of the inversion quality is low.
• CDF cumulative distribution function
• FIG. 4 is a block diagram illustrating a reservoir model system 500, in accordance with some embodiments. While certain specific features are illustrated, those skilled in the art will appreciate from the present disclosure that various other features have not been illustrated for the sake of brevity and so as not to obscure more pertinent aspects of the embodiments disclosed herein.
• the reservoir modeling system 500 includes one or more processing units (CPUs) 502, one or more network interfaces 508 and/or other communications interfaces 503, memory 506, and one or more communication buses 504 for interconnecting these and various other components.
• the reservoir modeling system 500 also includes a user interface 505 (e.g., a display 505-1 and an input device 505-2).
• the communication buses 504 may include circuitry (sometimes called a chipset) that interconnects and controls communications between system components.
• Memory 506 includes high-speed random access memory, such as DRAM, SRAM, DDR RAM or other random access solid state memory devices; and may include non-volatile memory, such as one or more magnetic disk storage devices, optical disk storage devices, flash memory devices, or other non-volatile solid state storage devices. Memory 506 may optionally include one or more storage devices remotely located from the CPUs 502. Memory 506, including the non-volatile and volatile memory devices within memory 506, comprises a non-transitory computer readable storage medium and may store seismic data, velocity models, seismic images, and/or geologic structure information.
• memory 506 or the non-transitory computer readable storage medium of memory 506 stores the following programs, modules and data structures, or a subset thereof including an operating system 516, a network communication module 518, and a seismic imaging module 520.
• the operating system 516 includes procedures for handling various basic system services and for performing hardware dependent tasks.
• the network communication module 518 facilitates communication with other devices via the communication network interfaces 508 (wired or wireless) and one or more communication networks, such as the Internet, other wide area networks, local area networks, metropolitan area networks, and so on.
• communication network interfaces 508 wireless or wireless
• communication networks such as the Internet, other wide area networks, local area networks, metropolitan area networks, and so on.
• the quality module 520 executes the operations of method 100.
• Quality module 520 may include data sub-module 525, which handles the seismic dataset including seismic data 525-1 and well data 525-2. This data is supplied by data sub-module 525 to other sub-modules.
• Attribute vector sub-module 522 contains a set of instructions 522-1 and accepts metadata and parameters 522-2 that will enable it to execute operation 12 of method 100.
• the data analytics sub-module 523 contains a set of instructions 523-1 and accepts metadata and parameters 523-2 that will enable it to contribute to operations 14 and 16 of method 100. Although specific operations have been identified for the sub-modules discussed herein, this is not meant to be limiting. Each sub-module may be configured to execute operations identified as being a part of other sub-modules, and may contain other instructions, metadata, and parameters that allow it to execute other operations of use in processing seismic data and generate the seismic image. For example, any of the sub-modules may optionally be able to generate a display that would be sent to and shown on the user interface display 505-1. In addition, any of the data or processed data products may be transmitted via the
• communication interface(s) 503 or the network interface 508 may be stored in memory 506.
• Method 100 is, optionally, governed by instructions that are stored in computer memory or a non-transitory computer readable storage medium (e.g., memory 506 in Figure 4) and are executed by one or more processors (e.g., processors 502) of one or more computer systems.
• the computer readable storage medium may include a magnetic or optical disk storage device, solid state storage devices such as flash memory, or other non-volatile memory device or devices.
• the computer readable instructions stored on the computer readable storage medium may include one or more of: source code, assembly language code, object code, or another instruction format that is interpreted by one or more processors.
• some operations in each method may be combined and/or the order of some operations may be changed from the order shown in the figures.
• method 100 is described as being performed by a computer system, although in some embodiments, various operations of method 100 are distributed across separate computer systems.
• the phrase “if it is determined [that a stated condition precedent is true]” or “if [a stated condition precedent is true]” or “when [a stated condition precedent is true]” may be construed to mean “upon determining” or “in response to determining” or “in accordance with a determination” or “upon detecting” or “in response to detecting” that the stated condition precedent is true, depending on the context.
• stages that are not order dependent may be reordered and other stages may be combined or broken out. While some reordering or other groupings are specifically mentioned, others will be obvious to those of ordinary skill in the art and so do not present an exhaustive list of alternatives. Moreover, it should be recognized that the stages could be implemented in hardware, firmware, software or any combination thereof.

Landscapes

• Engineering & Computer Science (AREA)
• Physics & Mathematics (AREA)
• Remote Sensing (AREA)
• Life Sciences & Earth Sciences (AREA)
• General Physics & Mathematics (AREA)
• General Life Sciences & Earth Sciences (AREA)
• Geophysics (AREA)
• Environmental & Geological Engineering (AREA)
• Geology (AREA)
• Acoustics & Sound (AREA)
• Theoretical Computer Science (AREA)
• Mathematical Optimization (AREA)
• Data Mining & Analysis (AREA)
• Mathematical Physics (AREA)
• Pure & Applied Mathematics (AREA)
• Computational Mathematics (AREA)
• Mathematical Analysis (AREA)
• General Engineering & Computer Science (AREA)
• Computing Systems (AREA)
• Algebra (AREA)
• Databases & Information Systems (AREA)
• Software Systems (AREA)
• Computer Hardware Design (AREA)
• Evolutionary Computation (AREA)
• Geometry (AREA)
• Geophysics And Detection Of Objects (AREA)

Applications Claiming Priority (2)

Publications (2)

Family

ID=71608034

Family Applications (1)

Country Status (4)

Families Citing this family (10)

Family Cites Families (4)

• 2019
  • 2019-07-25 US US16/522,445 patent/US10908308B1/en active Active
• 2020
  • 2020-07-08 EP EP20740095.3A patent/EP4004607B1/de active Active
  • 2020-07-08 AU AU2020104491A patent/AU2020104491B4/en active Active
  • 2020-07-08 AU AU2020319015A patent/AU2020319015A1/en active Pending
  • 2020-07-08 WO PCT/IB2020/056406 patent/WO2021014252A1/en not_active Ceased

Also Published As

Similar Documents

Legal Events