EP3769126A1 - System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions - Google Patents
System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversionsInfo
- Publication number
- EP3769126A1 EP3769126A1 EP19721762.3A EP19721762A EP3769126A1 EP 3769126 A1 EP3769126 A1 EP 3769126A1 EP 19721762 A EP19721762 A EP 19721762A EP 3769126 A1 EP3769126 A1 EP 3769126A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- volume
- seismic
- attribute
- interest
- reservoir
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 90
- 229930195733 hydrocarbon Natural products 0.000 title claims description 28
- 150000002430 hydrocarbons Chemical class 0.000 title claims description 27
- 238000004458 analytical method Methods 0.000 claims abstract description 33
- 230000000007 visual effect Effects 0.000 claims abstract description 5
- 239000012530 fluid Substances 0.000 claims description 29
- 239000004215 Carbon black (E152) Substances 0.000 claims description 19
- 239000000203 mixture Substances 0.000 claims description 18
- 238000003860 storage Methods 0.000 claims description 15
- 238000004519 manufacturing process Methods 0.000 claims description 12
- 239000012267 brine Substances 0.000 claims description 11
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 claims description 11
- 238000002347 injection Methods 0.000 claims description 10
- 239000007924 injection Substances 0.000 claims description 10
- 239000011148 porous material Substances 0.000 claims description 10
- 230000004044 response Effects 0.000 description 39
- 238000012360 testing method Methods 0.000 description 11
- 230000008859 change Effects 0.000 description 10
- 238000004891 communication Methods 0.000 description 10
- 238000009826 distribution Methods 0.000 description 9
- 238000005259 measurement Methods 0.000 description 9
- 239000004576 sand Substances 0.000 description 9
- 238000012545 processing Methods 0.000 description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 239000000463 material Substances 0.000 description 7
- 238000013459 approach Methods 0.000 description 6
- 230000008021 deposition Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000000704 physical effect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 239000004927 clay Substances 0.000 description 4
- 239000007787 solid Substances 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000010207 Bayesian analysis Methods 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 238000010793 Steam injection (oil industry) Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010494 dissociation reaction Methods 0.000 description 1
- 230000005593 dissociations Effects 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 125000001183 hydrocarbyl group Chemical group 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
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/308—Time lapse or 4D effects, e.g. production related effects to the formation
-
- 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/612—Previously recorded data, e.g. time-lapse or 4D
-
- 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/614—Synthetically generated data
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
- G01V2210/632—Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]
Definitions
- the present disclosure relates generally to methods and systems for probabilistic analysis of geologic features using seismic data and, in particular, methods and systems for assessing the probability of hydrocarbons in a subterranean reservoir based on seismic inversions generated from two or more seismic surveys performed at different times.
- 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 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.
- shots individual instances of the seismic source being activated
- One step in conventional processing of seismic reflection data involves adding multiple seismic traces that share a common mid-point, but have different source-receiver offsets. This is commonly called“stacking”. Stacking generally improves the signal to noise ratio, but can result in ambiguity surrounding the cause of the seismic amplitudes. For example, a high seismic amplitude could indicate either the presence of fluids or the presence of a particular lithology.
- AVO amplitude versus offset
- AVO angle
- AVO amplitude versus offset
- this type of AVA data may not be stacked thereby to preserve information that can be used to distinguish indicators of fluids from indicators of lithology.
- a hydrocarbon-bearing sand may generally have an increasingly negative seismic amplitude at further source- receiver offsets compared to a water-bearing sand which may be indicated by a decrease in positive seismic amplitude at further source-receiver offsets.
- the production of hydrocarbons causes changes in the elastic parameters of the earth. These changes may occur due to water displacing oil (or vice versa), water displacing gas (or vice versa), or gas displacing oil (or vice versa), within the reservoir interval. In other cases, the changes in the elastic parameters may occur due to enhanced hydrocarbon recovery operations, C02 injection, or clathrate dissociation from solid to gas.
- Time-lapse (4D) seismic data is acquired to compare seismic data at different times via two or more seismic surveys, a seismic survey at time one (Ti) and another seismic survey from time two (T 2 ), conducted months or years apart. The differences in the seismic responses for Ti and T 2 are at least partially due to fluid movement and/or pressure changes due to production or injection of water or gas.
- the seismic survey from Ti is referred to as the baseline survey
- the seismic survey from T 2 is referred to as the monitor survey.
- the seismic survey from Ti is an early monitor survey
- the seismic survey from T 2 is another monitor survey recorded at some time T 2 where T 2 is months or years after Ti .
- the above methods may however often be biased and may not truly represent the geologic features.
- conventional methods may fail where seismic data quality is low, such as where random and/or coherent noise is prevalent, or where seismic gathers are not flat.
- Project cost is dependent upon accurate prediction of the position of physical boundaries and fluid content within the Earth. 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. These decisions also include identifying locations for producing wells and injection wells, as well as how to adjust production rates or injection rates to optimize production over time.
- a computer-implemented method of reservoir property assessment in a subterranean volume of interest including receiving a first digital seismic attribute volume inverted from a seismic dataset recorded at a first time representative of a subsurface volume of interest, a second digital seismic attribute volume inverted from a seismic dataset recorded at a second time representative of the subsurface volume of interest, and a range of geological and geophysical parameters possible in the subsurface volume of interest; identifying at least one layer of interest; identifying at least one spatial area of interest for the at least one layer to define a reservoir volume of interest; computing an attribute difference volume from the first digital seismic attribute volume and the second digital seismic attribute volume within the reservoir volume of interest;
- 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 analyzing geologic features using seismic data, in accordance with some embodiments
- Figure 2 is a flowchart of one step from an embodiment
- Figures 3 - 6 are examples of other steps from various embodiments.
- Figure 7 illustrates steps and results from an embodiment
- Figure 8 is a block diagram illustrating a time-lapse fluid assessment system, in accordance with some embodiments.
- Figure 9 illustrates a flowchart of a method of analyzing geologic features using seismic attributes calculated by seismic inversions, in accordance with other embodiments.
- Figure 10 is a flowchart of one step from an embodiment.
- Like reference numerals refer to corresponding parts throughout the drawings.
- Described below are methods, systems, and computer readable storage media that provide a manner of geologic analysis using seismic data. These embodiments are designed to calculate probabilities of hydrocarbons (i.e. fluid property estimation) in subsurface geologic features and changes in those probabilities after production and/or injection.
- Probabilities of hydrocarbons i.e. fluid property estimation
- Industry standard techniques use deterministic estimation of the underlying geologic and geophysical parameters which contribute to the amplitude versus angle response utilizing forward modeling or inversion.
- the subsurface parameters of interest are the thickness, pore fluid (brine, oil, gas), hydrocarbon saturation, porosity, etc.
- the present method combines probabilistic AVA/AVO (amplitude versus angle/amplitude versus offset) and spatial summation of amplitude versus offset gathers with a Bayesian analysis to determine the range of geologic and geophysical parameters that will fit a user-selected range of measured field responses with selected areas.
- the probabilistic estimation builds a model space with a regular grid, then a singular bin is located for a given seismic trace and the property estimation is based on counting models in that singular bin.
- the present invention allows boxes based on the seismic data to be defined in the model space based on the probabilistic analysis from which the property estimation is done by counting models in the boxes.
- Seismic imaging of the subsurface is used to identify potential hydrocarbon reservoirs.
- Seismic data is acquired at a surface (e.g. the earth’s surface, ocean’s surface, or at the ocean bottom) as seismic traces which collectively make up the seismic dataset.
- the seismic dataset may be processed and imaged via a pre-stack method in order to analyze the seismic amplitude versus angle (AY A) or offset (AVO).
- Seismic surveys generally called a baseline survey and one or more monitor surveys, conducted at different times (months or years apart) are used to monitor changes in the subsurface and are processed and imaged to create images that will show differences in seismic amplitudes.
- the present invention includes embodiments of a method and system for assessing changes in reservoir properties over a period of time in a subterranean reservoir to determine the probability of hydrocarbons remaining after production and/or injection, in some embodiments estimating the probability of various saturation changes and/or pressure changes.
- Saturation changes are used to describe fluid changes in reservoirs that contain more than one type of fluid or gas or reservoirs where one fluid is partially replacing another as a result of hydrocarbon production or injection.
- Reservoir properties may include at least one of pore fluid content, porosity, brine composition, hydrocarbon composition, pressure, temperature, or any combination thereof.
- Determining the most probable changes in reservoir properties in a geologic feature and a range of possible changes in reservoir properties allows strategic planning around 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, selecting and purchasing appropriate completion and production equipment, and enhanced production strategies such as water or steam injection, as well as ultimately drilling into an optimum location to produce the hydrocarbons.
- FIG. 1 illustrates a flowchart of a method 100 for time-lapse geologic analysis of a subsurface volume of interest.
- At operation 10 at least two seismic datasets that were recorded at different times (i.e., baseline and monitor datasets), generally months or years apart, are received.
- the earlier seismic dataset is generally referred to as the baseline dataset and the subsequent datasets are monitor datasets.
- a seismic dataset includes a plurality of traces recorded at a plurality of seismic sensors. Due to changes in the reservoir properties caused by hydrocarbon production and/or injection, the seismic responses recorded at the two different times will be different in affected areas.
- Method 100 moves on to process the seismic datasets 11 using substantially similar processing flows to create digital seismic images.
- These datasets may be subjected to a number of seismic processing steps, such as deghosting, multiple removal, spectral shaping, and the like, before undergoing a pre-stack seismic imaging process.
- seismic processing steps such as deghosting, multiple removal, spectral shaping, and the like.
- the processing should preserve the seismic signal and reduce noise.
- the resultant digital seismic images may be, for example, a pre-stack seismic image, one or more seismic angle stacks, or one or more digital seismic horizon amplitude maps.
- the seismic horizon amplitude maps may have been computed at a series of angles (or summation of adjacent angles) in place of migrated seismic gathers.
- the seismic amplitude maps are computed by extracting the seismic amplitude from the migrated seismic gathers (either exact amplitude, or a computation of seismic amplitude at times around the horizon computed as average, absolute, rms, maximum, minimum, or other computational method) at the interpreted horizon time.
- the seismic horizons may be represented in time or depth, being optionally flattened.
- the seismic image and seismic horizons may be two-dimensional (2-D) (e.g., a horizontal dimension“x” and a time or depth dimension “z”) or three-dimensional (3-D) data sets (e.g., two perpendicular horizontal dimensions“x” and“y” and a time or depth dimension“z”).
- the seismic horizon may be representative of the top of a hydrocarbon reservoir (top sand) and/or the base of the hydrocarbon reservoir (base sand).
- the seismic images are interpreted to identify at least one spatial area on a seismic horizon that has differing amplitudes between the two seismic images.
- the seismic horizon should be representative of the reservoir that is being monitored. An example of this can be seen in Figure 3.
- Figure 3 shows map-view panels of a seismic horizon from a first time 30 and the same seismic horizon from a second time 32.
- a full range of seismic amplitude data has been stacked, which in this example embodiment is seismic amplitude data between angles 4° and 60°, as part of a data preprocessing step.
- the map indicates different regions of varying seismic amplitudes (indicated in differing shades) mostly correlating with the distribution of lithology, as well as liquids and gas, e.g., hydrocarbons.
- Interpretation of the seismic horizons shows that most of the amplitudes do not change between the two surveys.
- the amplitudes do change, so this area is selected for analysis by the rest of method 100.
- the amplitude difference in this area is calculated at operation 13.
- one or more areas of interest are identified on the seismic horizons.
- the areas of interest may be identified on a map view of the one or more seismic horizons, e.g., as polygons, wherein the map view may be colored (or shaded or contoured) to indicate the seismic amplitudes along the particular horizon.
- a probabilistic analysis is performed for the seismic amplitude versus angle (AVA) responses in at least one spatial area identified in the seismic image on at least one seismic horizon.
- An example of a method for doing this probabilistic analysis is shown in Figure 2 as method 200.
- This method may include, for example, using the method of US 2016/0209531, System and Method for Analyzing Geologic Features Using Seismic Data, which is incorporated herein in its entirety.
- a pre-stack seismic image contains multiple seismic horizons that represent seismic events identified or selected, in an embodiment, by a user as being of interest. These seismic horizons may represent a single thin lithology, such as a sand layer or a shale layer, or an interface within one or between two or more lithologies.
- each area of interest may encase a large number of seismic trace locations.
- a sufficient number of seismic trace locations may be on the order of thousands of trace locations.
- the statistical data ranges are influenced and determined by a range of geology enclosed in the selected area of interest (i.e. polygon) and noise.
- the range of geology may include, for example, changes in thickness, porosity, grain size, cementation, mineralogical composition, or the like.
- Statistical stability of the data is ensured by making the area of interest (polygon) sufficiently large to ensure that the noise is averaged out, as well as large enough to contain a representative sampling of the geology.
- the statistical data ranges may be represented by P50 and an upper and a lower probabilistic value for seismic amplitudes, each of the upper and lower values being similarly offset from the P50 value.
- the upper and lower probabilistic values may respectively be selected as a P10 and a P90 probabilistic value, a P20 and a P80 probabilistic value, a P30 and a P70 probabilistic value, or the like. These values are provided by way of example only and are not meant to be limiting.
- the P50 probabilistic value represents the underlying signal, while the upper and lower probabilistic values are indicative of a probabilistic range which represents the variable geology and/or noise.
- a variety of statistics may be computed from the aggregated seismic amplitudes, i.e. in addition, or alternatively, to the probabilistic values mentioned above.
- the statistical data ranges may include one or more of an average or mean (such as an average absolute amplitude), a mode, or a standard deviation such as RMS amplitude. It will be appreciated that other statistical measures may also be used.
- the use of many seismic trace locations from the areas of interest may assist in obtaining statistically significant data, in that the data may be more stable and distinct.
- angle stacks may be created by summing the seismic traces for each time or depth sample at two or more angles, e.g., adjacent angles.
- the angle stacks may be narrow, summing over a few adjacent angles, or broad, summing over many angles such as 10° - 20°.
- the ranges of angles summed over may overlap (e.g., 10° - 20° and 15° - 25°).
- a normalization based on the number of traces summed may be used in order to obtain an optimum presentation of the results.
- these narrow angle stacks may in some instances stabilize the trend of the AVA curves produced. It will however be appreciated that in many cases there may be no need for this type of stacking.
- the statistical data ranges may be based on other criteria such as the gradient or rate of change of the seismic amplitude response with angle or other industry-recognized measurements in the field (e.g., fanfar, grenv).
- FIG. 2 shows an embodiment of a method 200 for performing operation 14 of Figure 1.
- Operation 20 receives the baseline, monitor, amplitude difference datasets from the previous operations of method 100.
- Operation 21 of method 200 determines possible ranges of geological and geophysical parameters expected in the reservoir zone being analyzed that affect the seismic amplitude versus angle response.
- the expected ranges of geological and geophysical parameters are determined by the user based on nearby known information (e.g., previously drilled wells), estimated from theoretical equations, or other such information sources to provide results which may best characterize the expected geological and geophysical parameters expected in the reservoir zone.
- Geophysical parameters may include elastic properties such as P-wave velocity (Vp), S-wave velocity (Vs), and density.
- Geological parameters may include brine composition, hydrocarbon composition, pressure, temperature, porosity, reservoir thickness, mineralogical composition, and other factors. These determinations may be done by regional analysis, geologic inference or analogs, petrophysical analysis from analog well logs, or other means. In one
- Geological parameters may be determined, for example, for a situation in which there is advance knowledge of the deposition environment of the material. In this case, that knowledge may allow the user to determine information regarding what types of materials are likely to be present as well as what relationship various layers are likely to have.
- an eolian deposition environment would tend to include sandstones that are relatively free of clay and relatively well-sorted.
- deltaic sandstones would tend to be higher in clay content.
- the types of sandstone generated would depend, at least in part, on whether the region under investigation includes wind-deposited or river delta deposited material and could be further differentiated based on specifics of the deposition environment.
- Geophysical parameters may be determined, for example, where there is local information available, such as from well cores or well logs from nearby wells.
- the seismic amplitude versus angle (AVA) responses are calculated in at least one spatial area for each of the baseline, monitor, and amplitude difference datasets. This may be done, for example, using the method of US 2016/0209531, System and Method for Analyzing Geologic Features Using Seismic Data, which is incorporated herein in its entirety.
- a pre-stack seismic image contains multiple seismic horizons that represent seismic events identified or selected, in an embodiment, by a user as being of interest. These seismic horizons may represent a single thin lithology, such as a sand layer or a shale layer, or an interface within one or between two or more lithologies.
- AVA probabilities are calculated for the baseline, monitor, and amplitude difference AVA responses. These computations and calculations may be performed by reading seismic angle gathers, i.e. all of the seismic traces at a particular angle for an area of interest, identifying a time gate centered on the seismic horizon, and computing the aggregated amplitudes at each angle.
- the time gate has the effect of isolating a portion of each selected trace around a feature of interest in time. This process of computing the statistical data ranges for the seismic amplitudes in each of the areas of interest is
- the statistical data ranges may be represented by P50 and an upper and a lower probabilistic value for seismic amplitudes, each of the upper and lower values being similarly offset from the P50 value.
- the upper and lower probabilistic values may respectively be selected as a P10 and a P90 probabilistic value, a P20 and a P80 probabilistic value, a P30 and a P70 probabilistic value, or the like. These values are provided by way of example only and are not meant to be limiting.
- operation 23 proceeds to perform a full range of 2-layer or 3 -layer forward modeling with all combinations of the geological and geophysical parameters. This may be done, for example, using a method such as that described in US Patent 7,869,955, Subsurface Prediction Method and System, which is incorporated herein in its entirety.
- pseudo-wells including multiple types of synthetic well logs may be generated. Pseudo-wells may include physical properties such as Vp, Vs, density, porosity, shale volume (Vshale), saturation, pore fluid type or other properties.
- seismic models for the reservoir response first at conditions represented by the first seismic survey and then at a range of conditions representing expected changes in the reservoir properties that encompass the expected or measured properties represented by the time of the second seismic survey.
- These properties can be fluid saturation (brine, oil, gas), pressure, temperature, etc.
- These property changes should be represented by a number of discrete changes. In an embodiment, this may be a small number of discrete changes such as 2 - 5.
- the modeling of the reservoir at the initial state may include variations in reservoir thickness, porosity, and other properties.
- this operation may also construct a series of results of the amplitude difference between the reservoir properties corresponding to the first (i.e., baseline) seismic data set and the suspected discrete parameter changes represented by the second (i.e., monitor) seismic data set.
- the AVA probabilities are calculated at operation 25.
- the pseudo-wells may be generated using a partially random approach. Rather than using a simple stochastic approach, in which any particular physical model is equally likely, the generation of the pseudo-wells may be constrained by physical constraints. The constraining may take place prior to the generating, or alternately, purely stochastic pseudo- wells may be later constrained (e.g., by eliminating wells having characteristics outside the constraints). As will be appreciated, it is likely to be more efficient to first constrain, then generate, the wells, but either approach should be considered to be within the scope of the present invention.
- the forward modeling of operation 23 will produce modeled (i.e. synthetic) seismic gathers containing AVA effects for the various combinations of geological and geophysical parameters. Forward modeling may be done, for example, using some form of the Zoeppritz equation, full waveform modeling, or other such seismic modeling method that may be appropriate including that explained by US Patent 7,869,955. Then at operation 25, these synthetic seismic gathers are used to calculate the probability of various seismic amplitudes within the area of interest, thereby allowing the statistical data ranges of seismic amplitudes to be determined.
- Figure 4 shows AVA curves for three different fluid contents (brine/wet, fizz, and gas), including the P50 values (the triangle, star, and square symbols) with range bars indicating the P20 - P80 ranges.
- Fizz is generally considered to be a low saturation, non-commercial amount of hydrocarbon gas (1% to 15% gas saturation) contained in the rock pores along with formation water.
- the seismic amplitude responses should be determined separately for brine, low and high hydrocarbon saturation, and different hydrocarbon fluids.
- the measured response ranges may also be segregated by different geological assessment of the mineralogical composition of the reservoir and non reservoir rocks (i.e. facies) simulated in the forward modeling step. Other examples of the forward modeled responses can be seen in Figure 5.
- the different grayscale dots indicate amplitudes as very-far-stack vs. near-stack for different fluid contents. Boxes defining the baseline amplitudes (amplitudes at the earlier time), monitor amplitudes (amplitudes at the later time), and difference amplitudes are based on the AVA probabilities calculated in operation 24, calculated from the input digital seismic images, are shown.
- Figure 6 shows a similar plot of the baseline and monitor boxes but the forward modeled results have been simplified to the modeled fluid vector rather than the grayscale dots.
- Figure 6 shows a similar plot of the baseline and monitor boxes but the forward modeled results have been simplified to the modeled fluid vector rather than the grayscale dots.
- an ellipse or other such shape to represent the spatial distribution of the data about the central value.
- a mathematical distribution characterizing the distribution of the data around the P50 amplitude could be estimated and used from operation 24 and forward in the analysis.
- Figure 6 shows the box in two dimensions, the box (or ellipse or mathematical distribution) may be multi-dimensional. For example, if statistical data ranges are found for four different angles, the box would have four dimensions.
- Method 200 can now proceed to operation 26, estimating the probability for changes in pore fluid saturation based on comparison of the calculated AVA probabilities from operation 24 and the calculated modeled AVA probabilities from operation 25. This estimation is done by comparing the amplitude difference AVA responses and the baseline and monitor AVA responses. By way of example, this may be done by using a two-box or three-box test, to estimate the change in reservoir properties in each polygon separately by considering the difference in the measured seismic amplitude versus angle responses between the first (baseline) and second (monitor) seismic survey within a single spatial polygon.
- the successful seismic models must have a calculated response which fits all of the measured response ranges.
- the result of method 200 can be used to estimate changes in reservoir properties.
- Seismic models can be used to estimate other reservoir properties.
- the estimated reservoir properties may include pore fluid content, porosity, brine composition, hydrocarbon composition, pressure, temperature, or any combination thereof. These estimated reservoir properties are estimates of the average geology in the spatial area of interest. Reservoir properties such as porosity, thickness, and Vshale can be estimated from these seismic models. Statistical measurements can be computed and summarized.
- interval i.e. layer-based
- interval i.e. layer-based
- time thickness is a candidate for the seismic attribute, or average amplitude for an interval.
- Time thickness could be used as a seismic attribute to assure that the model and seismic thickness are broadly matched.
- time thickness for the seismic traces could be used to adjust the priors for the reservoir thickness of the models (greater than tuning thickness, less than tuning thickness, a mixture). Average amplitude for an interval could be useful for estimation of NTG (net to gross).
- FIG. 9 is a flowchart illustrating method 900 for estimating changes in reservoir properties from seismic attributes calculated by seismic inversion.
- At operation 90 of Figure 9 at least two seismic attribute volumes
- these seismic attribute volumes may include, by way of example and not limitation, one or more of AI, SI, p, El, lr and m, as calculated by seismic inversions of the seismic datasets recorded at the first and second times.
- AI seismic attribute volume
- AI multiple seismic attribute volumes
- the same type of seismic attribute volumes must be received for both the baseline and monitor times.
- a spatial area and layer is identified where the seismic attribute changes between the baseline and monitor volumes. This is similar to operation 12 of method 100 from Figure 1.
- the spatial area is selected on a given seismic layer (i.e. interval).
- this layer may be a reservoir layer such as a sand layer.
- the layer could be defined based on seismic amplitudes or seismic impedance.
- the layer top would be defined as being at the high-to-low impedance transition and the layer base would be at the low-to-high impedance transition.
- the layer thickness may be greater than one sample, so an average attribute value or other representative value (e.g., root-mean- square (RMS)) may be used to represent the seismic attribute for the layer at a particular location.
- RMS root-mean- square
- operation 92 calculates the difference between the baseline attribute volume and the monitor attribute volume. This is done as a point-by-point subtraction to produce an attribute difference volume.
- the attribute difference volume, baseline attribute volume, and monitor attribute volume for the spatial area are input to operation 93 which performs probabilistic attribute analysis.
- Figure 10 illustrates a method 1000 for performing probabilistic attribute analysis.
- Operation 1002 receives the attribute difference volume, baseline attribute volume, and monitor attribute volume for the spatial area.
- Operation 1006 performs statistical analysis of the volumes to calculate statistical data ranges for the attribute volumes. These computations and calculations may be performed by reading all of the attributes for an area of interest and within the layer and then determining statistical data ranges of the attributes.
- the statistical data ranges may be represented by P50 and an upper and a lower probabilistic value for seismic attributes, each of the upper and lower values being similarly offset from the P50 value.
- the upper and lower probabilistic values may respectively be selected as a P10 and a P90 probabilistic value, a P20 and a P80 probabilistic value, a P30 and a P70 probabilistic value, or the like. These values are provided by way of example only and are not meant to be limiting.
- Operation 1001 of method 1000 determines possible ranges of geological and geophysical parameters expected in the reservoir layer being analyzed that affect the seismic attributes.
- the expected ranges of geological and geophysical parameters are determined by the user based on nearby known information (e.g., previously drilled wells), estimated from theoretical equations, or other such information sources to provide results which may best characterize the expected geological and geophysical parameters expected in the reservoir zone.
- Geophysical parameters may include elastic properties such as P-wave velocity (Vp), S-wave velocity (Vs), and density.
- Geological parameters may include brine composition, hydrocarbon composition, pressure, temperature, porosity, reservoir thickness, mineralogical composition, and other factors.
- Geological parameters may be determined, for example, for a situation in which there is advance knowledge of the deposition environment of the material. In this case, that knowledge may allow the user to determine information regarding what types of materials are likely to be present as well as what relationship various layers are likely to have.
- an eolian deposition environment would tend to include sandstones that are relatively free of clay and relatively well-sorted.
- deltaic sandstones would tend to be higher in clay content.
- the types of sandstone generated would depend, at least in part, on whether the region under investigation includes wind-deposited or river delta deposited material and could be further differentiated based on specifics of the deposition environment.
- Geophysical parameters may be determined, for example, where there is local information available, such as from well cores or well logs from nearby wells.
- operation 1003 proceeds to generate synthetic seismic attributes with all combinations of the geological and geophysical parameters. This may be done, for example, using a method such as that described in US Patent 7,869,955, Subsurface Prediction Method and System, which is incorporated herein in its entirety.
- pseudo-wells including multiple types of synthetic well logs may be generated. Pseudo-wells may include physical properties such as Vp, Vs, AI, SI, p, El, lr and m, porosity, shale volume (Vshale), saturation, pore fluid type or other properties.
- seismic attributes are calculated for the reservoir first at conditions represented by the first seismic survey and then at a range of conditions representing expected changes in the reservoir properties that encompass the expected or measured properties represented by the time of the second seismic survey.
- These properties can be fluid saturation (brine, oil, gas), pressure, temperature, etc.
- These property changes should be represented by a number of discrete changes. In an embodiment, this may be a small number of discrete changes such as 2 - 5.
- the parameters of the reservoir at the initial state may include variations in reservoir thickness, porosity, and other properties.
- Operation 1003 will calculate the synthetic seismic attributes that were the output from seismic inversions, such as one or more of AI, SI, p, El, lr and m, and the difference between the reservoir properties corresponding to the first (i.e., baseline) seismic data set and the suspected discrete parameter changes represented by the second (i.e., monitor) seismic data set. Once the synthetic attributes have been calculated at operation 1003, the synthetic attributes are used in operation 1007 to determine the synthetic seismic attributes that were the output from seismic inversions, such as one or more of AI, SI, p, El, lr and m, and the difference between the reservoir properties corresponding to the first (i.e., baseline) seismic data set and the suspected discrete parameter changes represented by the second (i.e., monitor) seismic data set. Once the synthetic attributes have been calculated at operation 1003, the synthetic attributes are used in operation 1007 to determine
- the pseudo-wells may be generated using a partially random approach. Rather than using a simple stochastic approach, in which any particular physical model is equally likely, the generation of the pseudo-wells may be constrained by physical constraints. The constraining may take place prior to the generating, or alternately, purely stochastic pseudo- wells may be later constrained (e.g., by eliminating wells having characteristics outside the constraints). As will be appreciated, it is likely to be more efficient to first constrain, then generate, the wells, but either approach should be considered to be within the scope of the present invention. Additionally, if the seismic attribute volumes received in operation 90 of method 900 are bandlimited, the pseudowells can be filtered to match the bandwidth of the seismic attribute volumes.
- Operation 1003 will produce synthetic (i.e. modeled) seismic attributes for the various combinations of geological and geophysical parameters. Examples of the synthetic attributes can be plotted similarly to those seen in Figure 5.
- the different grayscale dots indicate attributes as average shear impedance (SI) vs. average acoustic impedance (AI) for different fluid contents.
- SI shear impedance
- AI average acoustic impedance
- “average” means the average value for the layer being analyzed, which is typically a reservoir layer.
- these plots could instead be AI vs. El or some other plot of at least two seismic attributes for each of the monitor and baseline datasets.
- Boxes defining the baseline attributes (attributes at the earlier time), monitor attributes (attributes at the later time), and difference attributes are based on the statistical data ranges calculated in operation 1006, calculated from the baseline, monitor, and difference attribute volumes, could be drawn.
- the boxes may be centered at the P50 values and the width of the boxes may be determined by the statistical ranges, such as but not limited to P80-P20.
- P50 values and statistical ranges for 2 attributes define a box.
- the boxes can have higher dimensions, where a N-dimensional box is defined by N inversion attributes.
- a box around the P50 attribute values to represent the range of possible models
- a mathematical distribution characterizing the distribution of the data around the P50 attribute value could be estimated and used from operation 1006 and forward in the analysis.
- the box (or ellipse or mathematical distribution) may be multi-dimensional. For example, if statistical data ranges are found for four different attributes, the box would have four dimensions.
- Method 1000 can now proceed to operation 1007, estimating the probability for changes in reservoir properties based on comparison of the calculated attribute statistical data ranges from operation 1006 and the calculated synthetic attributes from operation 1003. This estimation is done by comparing the attribute differences and the baseline and monitor attributes. By way of example, this may be done by using a two-box or three-box test, to estimate the change in reservoir properties in each polygon separately by considering the difference in the measured attributes between the first (baseline) and second (monitor) seismic survey within a single spatial polygon for the layer of interest.
- the successful seismic models must have a calculated attribute which fits all of the measured attribute ranges.
- modeled/synthetic seismic attribute volumes at each measured parameter and an extent determined by statistical measurements (e.g., P20 - P80, standard deviation, etc.). Any two or all 3 of the tests above can be used to determine the final subset of successful models. When doing only a two-box test, one of the tests should be using the attribute difference data set. Analyze the total number of successful synthetic attributes which fit either the two or three tests used. The probability of each reservoir property change (e.g. saturation case) at the time of the second (i.e., monitor) seismic survey according to this hypothesis is the number of successful attributes for that saturation case divided by the total number of successful attributes.
- the result of method 1000 can be used to estimate changes in reservoir properties in addition to saturation changes.
- seismic attributes can be used to estimate other reservoir properties.
- the estimated reservoir properties may include pore fluid content, porosity, brine composition, hydrocarbon composition, pressure, temperature, or any combination thereof. These estimated reservoir properties are estimates of the average geology in the spatial area of interest. Reservoir properties such as porosity, thickness, and Vshale can be estimated from these seismic attributes. Statistical measurements can be computed and summarized.
- Figure 7 is an example of the steps of method 100 using method 200.
- the structure maps show the baseline map 71 A and the monitor map 71B.
- Diagram 72 shows AVA probability curves created at operation 24 of method 200.
- Diagram 73 shows the AVA curves for the baseline seismic data, the monitor seismic data, and the amplitude difference between the baseline and monitor seismic data with angles selected for use in subsequent steps of method 200.
- Diagram 74 shows a box for known baseline oil, meaning that at the time of the baseline survey, the reservoir contained oil. Seismic models generated by method 200 at operations 23 and 25.
- Diagram 75 shows the three-box test described above.
- Diagram 76 shows the probabilities calculated as a result of method 100. Although these results are displayed graphically, this is not meant to be limiting. Other methods of presenting the results, such as in a spreadsheet format, are possible.
- FIG 8 is a block diagram illustrating a time-lapse reservoir property assessment 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 property assessment 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 property assessment 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 reservoir property 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.
- Time-lapse module 520 executes the operations of method 100.
- Time-lapse module 520 may include data sub-module 525, which handles the seismic image including seismic gathers 525-1 through 525-N. This seismic data is supplied by data sub-module 525 to other sub-modules.
- AVA (amplitude versus angle) sub-module 522 contains a set of instructions
- the forward modeling function sub-module 523 contains a set of instructions 523-1 and accepts metadata and parameters 523-2 that will enable it to execute parts of operation 14 of method 100.
- the fluid content sub-module 524 contains a set of instructions 524-1 and accepts metadata and parameters 524-2 that will enable it to execute at least operation 15 of method 100.
- 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.
- any of the seismic data or processed seismic data products may be transmitted via the communication interface(s)
- the network interface 508 may be stored in memory 506.
- time-lapse reservoir property assessment system 500 is designed for method 100, those of skill in the art will appreciate that a very similar system including an attribute sub-module rather than an AVA sub-module can be used for method 900.
- 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 8) 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. For ease of
- 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.
- 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)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Fluid Mechanics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Claims
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/928,130 US10718876B2 (en) | 2017-10-26 | 2018-03-22 | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions |
PCT/IB2019/052331 WO2019180669A1 (en) | 2018-03-22 | 2019-03-22 | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions |
Publications (1)
Publication Number | Publication Date |
---|---|
EP3769126A1 true EP3769126A1 (en) | 2021-01-27 |
Family
ID=66397321
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP19721762.3A Withdrawn EP3769126A1 (en) | 2018-03-22 | 2019-03-22 | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions |
Country Status (4)
Country | Link |
---|---|
EP (1) | EP3769126A1 (en) |
AU (1) | AU2019237361B2 (en) |
CA (1) | CA3092287A1 (en) |
WO (1) | WO2019180669A1 (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2588685B (en) * | 2019-11-04 | 2022-05-25 | Equinor Energy As | Hydrocarbon exploration method |
CN113376690B (en) * | 2020-03-09 | 2023-09-26 | 中国石油天然气股份有限公司 | Reservoir parameter prediction method and system |
CN112305614B (en) * | 2020-10-20 | 2024-03-29 | 中海石油(中国)有限公司 | Method and system for describing space spreading range of gas cloud area |
CN113640880B (en) * | 2021-08-18 | 2024-03-08 | 核工业二0八大队 | Sandstone type uranium ore seismic exploration method and device |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7636275B2 (en) * | 2007-02-06 | 2009-12-22 | Conocophillips Company | Direct time lapse inversion of seismic data |
US7869955B2 (en) * | 2008-01-30 | 2011-01-11 | Chevron U.S.A. Inc. | Subsurface prediction method and system |
US8379482B1 (en) * | 2009-01-13 | 2013-02-19 | Exxonmobil Upstream Research Company | Using seismic attributes for data alignment and seismic inversion in joint PP/PS seismic analysis |
US10215870B2 (en) * | 2014-05-08 | 2019-02-26 | Chevron U.S.A. Inc. | System and method for analyzing geologic features using seismic data |
US9817142B2 (en) | 2015-01-15 | 2017-11-14 | Chevron U.S.A. Inc. | System and method for analyzing geologic features using seismic data |
-
2019
- 2019-03-22 CA CA3092287A patent/CA3092287A1/en not_active Abandoned
- 2019-03-22 EP EP19721762.3A patent/EP3769126A1/en not_active Withdrawn
- 2019-03-22 AU AU2019237361A patent/AU2019237361B2/en active Active
- 2019-03-22 WO PCT/IB2019/052331 patent/WO2019180669A1/en active Application Filing
Also Published As
Publication number | Publication date |
---|---|
WO2019180669A1 (en) | 2019-09-26 |
CA3092287A1 (en) | 2019-09-26 |
AU2019237361B2 (en) | 2021-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2008349455B2 (en) | Subsurface prediction method and system | |
US10534100B2 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on time-lapse seismic data | |
AU2019237361B2 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions | |
CA3043310A1 (en) | Method for estimating petrophysical properties for single or multiple scenarios from several spectrally variable seismic and full wavefield inversion products | |
CA2940406C (en) | Characterizing a physical structure using a multidimensional noise model to attenuate noise data | |
US20190227185A1 (en) | System and method for predicting fault seal from seismic data | |
US10884149B2 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic data | |
US10634805B2 (en) | System and method for predicting fault seal from seismic data | |
EP3685194B1 (en) | System and method for deriving reservoir stresses from 4d seismic data | |
Naseer | Delineating the shallow‐marine stratigraphic traps of Lower‐Cretaceous incised valley sedimentation, Pakistan using post‐stack seismic colour inversion | |
EP4133310A1 (en) | System and method for seismic inversion | |
US10718876B2 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions | |
US20190146108A1 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic data | |
US10571585B2 (en) | System and method for time-lapsing seismic imaging | |
Wahbah et al. | Cross Functional Improvements in Modelling Complex Carbonate Reservoirs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: UNKNOWN |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE |
|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
|
17P | Request for examination filed |
Effective date: 20200825 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
AX | Request for extension of the european patent |
Extension state: BA ME |
|
RIN1 | Information on inventor provided before grant (corrected) |
Inventor name: MAGILL, JAMES R. Inventor name: BARTEL, DAVID C. |
|
DAV | Request for validation of the european patent (deleted) | ||
DAX | Request for extension of the european patent (deleted) | ||
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN |
|
18D | Application deemed to be withdrawn |
Effective date: 20221001 |
|
P01 | Opt-out of the competence of the unified patent court (upc) registered |
Effective date: 20230522 |