EP2689275B1 - Method and apparatus for analyzing data in time-lapse geophysical surveys - Google Patents

Method and apparatus for analyzing data in time-lapse geophysical surveys Download PDF

Info

Publication number
EP2689275B1
EP2689275B1 EP11741397.1A EP11741397A EP2689275B1 EP 2689275 B1 EP2689275 B1 EP 2689275B1 EP 11741397 A EP11741397 A EP 11741397A EP 2689275 B1 EP2689275 B1 EP 2689275B1
Authority
EP
European Patent Office
Prior art keywords
data
survey
attribute
attribute value
cmp
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.)
Active
Application number
EP11741397.1A
Other languages
German (de)
French (fr)
Other versions
EP2689275A1 (en
Inventor
John GREGOR
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ion Geophysical Corp
Original Assignee
Ion Geophysical Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ion Geophysical Corp filed Critical Ion Geophysical Corp
Publication of EP2689275A1 publication Critical patent/EP2689275A1/en
Application granted granted Critical
Publication of EP2689275B1 publication Critical patent/EP2689275B1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/308Time lapse or 4D effects, e.g. production related effects to the formation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/34Displaying seismic recordings or visualisation of seismic data or attributes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/34Displaying seismic recordings or visualisation of seismic data or attributes
    • G01V1/345Visualisation of seismic data or attributes, e.g. in 3D cubes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • G01V1/3808Seismic data acquisition, e.g. survey design
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • G01V1/3817Positioning of seismic devices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • G01V1/3843Deployment of seismic devices, e.g. of streamers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/612Previously recorded data, e.g. time-lapse or 4D
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/70Other details related to processing
    • G01V2210/74Visualisation of seismic data

Definitions

  • This disclosure relates generally to geophysical exploration systems, and more particularly to data acquisition and analysis in subsequent geophysical surveys.
  • seismic exploration techniques utilizing, among other things, seismic and other wave exploration techniques to find oil and gas reservoirs within the Earth.
  • seismic exploration techniques often include controlling the emission of seismic energy into the Earth with a seismic source of energy (e.g. , dynamite, air guns, vibrators, etc.), and monitoring the Earth's response to the seismic source with a receiver (e.g. , a geophone, a hydrophone, etc.).
  • a seismic source of energy e.g. , dynamite, air guns, vibrators, etc.
  • a receiver e.g. , a geophone, a hydrophone, etc.
  • a follow-up, or monitor survey may be taken near the same location as the original, baseline survey.
  • the data that the monitor survey provides is called 4D data, because it provides a fourth dimension to the data being collected-time.
  • this 4D survey may allow scientists and engineers to see how much oil or gas has been removed from the reservoir, how much oil or gas remains in the reservoir, and/or how the oil or gas has moved within the reservoir since the baseline survey.
  • the monitor survey may involve firing the sources and recording the reflected seismic signals at the receivers.
  • some system of quality control may be established to facilitate coordinating data sets between the two surveys.
  • Repeatability is a term used to refer to the measure of how closely the monitor survey correlates to the baseline survey.
  • positional or geometric repeatability is a measure of how positionally or geometrically close a monitor survey's data is to a baseline survey's data. Even if the geometric repeatability of a monitor survey is good (for example, where sources and receivers are installed permanently on land or on the ocean floor), other factors may decrease the repeatability, and therefore the reliability, of the data obtained in the monitor survey.
  • a few aspects of a survey that may be of concern in terms of repeatability include: source position, receiver position, the common mid-point (CMP) position between the source and the receiver, offset or radial distance between the source and the receiver, azimuth, and components and combinations of these.
  • CMP common mid-point
  • a prior system for measuring the repeatability between surveys is disclosed in the International Patent Application published under number WO2010/054282 (CONOCOPHILLIPS COMPANY).
  • a first aspect of the present invention relates to a method for processing geophysical data as defined in independent claim 1, and optionally any of dependent claims 2 to 13.
  • the method includes matching a first set of data from a monitor survey with a second set of data from a baseline survey. It also includes generating an attribute value for each datum in the first set, storing each attribute value with each of the datum, and processing the first set of data based on the stored attribute values.
  • Another method of processing geophysical data may include correlating a first set of data from a first survey with a second set of data from a second survey. It also includes generating an attribute value for each datum in the first set, associating each generated attribute value with the datum from the first set of data, and processing the first set of data based on the attribute values.
  • a second aspect of the present invention relates to a tangible storage medium as defined in independent claim 14, and optionally dependent claim 15.
  • the medium may store a plurality of instructions, including instructions that correlate a first set of data from a first survey with a second set of data from a second survey, instructions that generate an attribute value for each datum in the first set, instructions that associate each generated attribute value with the datum from the first set of data from which the attribute value was generated, and instructions that process the first set of data based on the attribute values.
  • Another tangible storage medium may store a filtered set of CMPs, wherein the filtered set of CMPs is filtered from a first set of CMPs, and the filtering is based on a set of attribute values.
  • seismic exploration techniques utilizing, among other things, seismic and other wave exploration techniques to find oil and gas reservoirs within the Earth.
  • seismic exploration techniques often include controlling the emission of seismic energy into the Earth with a seismic source of energy (e.g. , dynamite, air guns, vibrators, etc.), and monitoring the Earth's response to the seismic source with a receiver (e.g. , a geophone, a hydrophone, etc.).
  • a seismic source of energy e.g. , dynamite, air guns, vibrators, etc.
  • a receiver e.g. , a geophone, a hydrophone, etc.
  • a follow-up, or monitor survey may be taken near the same location as the original, baseline survey.
  • the data that the monitor survey provides is called 4D data, because it provides a fourth dimension to the data being collected-time.
  • this 4D survey may allow scientists and engineers to see how much oil or gas has been removed from the reservoir, how much oil or gas remains in the reservoir, and/or how the oil or gas has moved within the reservoir since the baseline survey.
  • the monitor survey may involve firing the sources and recording the reflected seismic signals at the receivers.
  • some system of quality control may be established to facilitate coordinating data sets between the two surveys.
  • Repeatability is a term used to refer to the measure of how closely the monitor survey correlates to the baseline survey.
  • positional or geometric repeatability is a measure of how positionally or geometrically close a monitor survey's data is to a baseline survey's data. Even if the geometric repeatability of a monitor survey is good (for example, where sources and receivers are installed permanently on land or on the ocean floor), other factors may decrease the repeatability, and therefore the reliability, of the data obtained in the monitor survey.
  • a few aspects of a survey that may be of concern in terms of repeatability include: source position, receiver position, the common mid-point (CMP) position between the source and the receiver, offset or radial distance between the source and the receiver, azimuth, and components and combinations of these.
  • CMP common mid-point
  • systems and methods are disclosed that provide scientists and engineers the ability to analyze monitor survey geophysical data and monitor survey coverage that have been filtered based on one or more repeatability attributes.
  • Figure 1A shows a top view of a vessel 101 towing a source 102 and several receivers 103 on streamers behind the vessel 101.
  • Figure 1B shows a side-view of the vessel 101 shown in Figure 1A with the source 102 and receivers 103 being towed behind the vessel 101 just beneath the surface of the water.
  • the embodiment depicted in Figures 1A and 1B illustrates the source and receiver being towed by the same vessel, however, other possible combinations are possible.
  • either the source and/or receivers may be towed by separate vessels or may be implemented in land-based acquisition systems.
  • the source and/or receivers may be stationary while the other is towed behind the vessel.
  • the receivers 103 may be positioned deeper in the water, for example, by using streamer steering devices, such as the DigiFIN® brand steering device available from ION Geophysical, Inc.
  • streamer steering devices such as the DigiFIN® brand steering device available from ION Geophysical, Inc.
  • multiple sources may be used.
  • any type of source(s) or receiver(s) may be used, including for example, 1-, 2-, or 3-dimensional sources or receivers.
  • the source 102 may emit seismic energy (e.g., through an air gun), which may reflect off various portions of the Earth 104 and may be received back at the receivers 103 (as shown by the propagating seismic waves in Figure 1B ).
  • the signal received and processed at the receivers 103 may provide data that is useful in determining the composition of various portions of the Earth 104 proximate the location where the signal was reflected, which may include an oil and/or gas reservoir 105. If the amount of oil and/or gas in the reservoir 105 is depleted over time, then subsequent surveys conducted in substantially the same location as the first survey may indicate various properties of this depletion such as: decreasing pore pressures, migration of oil/water and/or gas/water contacts, drop in impedance, and so forth.
  • Figures 2A through 2D illustrate operations 200 that may utilize a predetermined metric or attribute as a criteria for binning coverage.
  • a predetermined metric or attribute may be used as a filter to bin coverage, for the sake of discussion, the following disclosure illustrates the operations 200 in the context of using a repeatability metric as a filter.
  • some embodiments of operations 200 may include binning a CMP if a repeatability attribute for that CMP falls within a specified range.
  • the operations 200 may analyze data from the full offset range (i.e. data from every offset step, with an offset step being one or more receiver intervals) or data from a subset of the full offset range (such as only one offset step from within each offset zone, or only one offset zone within the full offset range).
  • Operation 205 illustrates the sub-operations 206-210 associated with conducting a baseline survey.
  • a seismic source may fire shots or bursts of seismic energy into the Earth 104 in sub-operation 206. This energy reflects back and is received and measured as a signal at a receiver at a specific offset at sub-operation 207, with the offset being, in some embodiments, the radial range from the source to a specific receiver or group of receivers.
  • operation 205 may include sub-operations 208 and 209 that include calculating and storing CMPs at each offset step. This is illustrated in Figure 1A with box 106 indicating a first offset where each of the receivers in the box 106 has a different CMP vis-à-vis the source 102.
  • Sub-operation 209 may also store related context information with each CMP, including, for example acquisition and/or seismic attributes such as source ID, receiver ID, source position, receiver position, time, water depth, vessel speed, offset, azimuth, and so forth.
  • each of the receivers 103 in box 108 (which are at a different offset than the receivers in box 106), may have a different CMP.
  • Operation 205 may in some embodiments also include sub-operation 210 where a grid of rectangles referred to as "bins" may be defined. Note that although the bins are illustrated in operation 205 as the final sub-operation, they may be defined at any point in the baseline survey.
  • the width of a bin may be substantially the same as the spacing between CMP lines and the height of a bin may be substantially the same as the spacing between CMPs on a single CMP line. This is illustrated in Figure 1A with a bin 110 indicated by a dashed box. Of course, bins of other dimensions also may be defined with different heights or widths, depending on the resolution required.
  • a visual representation of the baseline survey may be generated after "binning" each CMP by associating data from each source receiver pair with the bin where its CMP is located (this operation is not shown in Figs. 2A through 2D ).
  • a visual representation may be a spatial coverage plot (similar to the plot shown Fig. 3 ) or any other representation that conveys coverage information, such as charts, lists, graphs, and so forth.
  • Within each bin there may be CMPs for different receiver offset steps. If a bin has at least one CMP from each offset step, it may be referred to as "full fold.” In other words, full fold coverage is achieved by obtaining a complete distribution of source to receiver offset steps over the entire offset range.
  • Each bin may include, for example a count of how many CMPs are present at each offset step, or how "full" the fold is. Each bin may additionally or alternatively include a percentage of how full the fold is. In some cases, there may be more than one CMP per offset step, referred to as "duplicates," and these duplicates may be included or excluded during binning.
  • a monitor survey may be conducted, as shown in operation 215.
  • the amount of time that elapses between the baseline and monitor surveys is one or more years.
  • the monitor survey in operation 215, along with its sub-operations 216-219, may be similar to the baseline survey: seismic shots may be fired in sub-operation 216, seismic energy may be reflected off the Earth and be received at one or more receiver(s) in sub-operation 217, the signal may be recorded, a CMP calculated for each source receiver pair in sub-operation 218, and the CMP for each offset may be stored in sub-operation 219.
  • the baseline survey CMPs may be binned and a visual representation of the coverage of the monitor survey may be generated (this operation is not shown in Figs. 2A through 2D ).
  • Each bin may include, for example, a count of how many CMPs are present at each offset, or any other useful information.
  • Figure 3 illustrates an unfiltered coverage plot for the monitor survey that may result from the sub-operations of operation 215.
  • the x and y axes of Figure 3 represent a bin grid, as defined in operation 205.
  • the shading in of each rectangular region represents how "full" that particular bin is, as shown in the key at the bottom of Figure 3 .
  • bins that are white are 61% or more full.
  • the shading scale can be configured to be any type of shading, including different colors and so forth.
  • the thresholds in the key and plot can be set to any percentage range; they may also be set to display the actual fold value rather than a percentage, or any other useful information.
  • flow may proceed to operation 225 where data from the monitor survey may be matched to data from the baseline survey and one or more metrics may be generated.
  • the sup-operations of operation 225 may include first, for each CMP of each shot fired in the monitor survey, determining a CMP from the baseline survey that is the "closest match" in sub-operation 226. In some embodiments, this match may be based upon the closest geometric match.
  • One example of geometric closeness may be measured by the geometric repeatability value Dsrc + Drx. Dsrc is the difference in source positions between the baseline survey CMP and the monitor survey CMP at the time of firing the seismic source respectively.
  • Drx is the difference in the receiver positions between the baseline survey CMP and the monitor survey CMP at the time of firing the seismic source.
  • the positions of the source and/or receiver may be with reference to a navigation system aboard the vessel 101, such as by using a Global Positioning System (GPS).
  • GPS Global Positioning System
  • Dsrc + Drx represents the sum of these two differences. Evaluating Dsrc + Drx values in order to determine the closest CMP match may ensure that the same azimuth of data is obtained so as to be referencing substantially the same area of the subsurface, which may allow, for example, the characteristics of the reservoir 105 to be consistently characterized over time.
  • each CMP in the monitor survey After each CMP in the monitor survey has been matched to a CMP in the baseline survey, the matched pairs of CMPs may be compared again in sub-operation 227 in order to generate an attribute value in sub-operation 228, and those attribute values may be stored with the CMP hits of the monitor survey in sub-operation 229.
  • the attribute value may be called a repeatability metric, or a repeatability attribute, if it represents a comparison between a baseline survey CMP and a monitor survey CMP; repeatability referring to how well the monitor survey 'repeated' the baseline survey. Accordingly, each CMP in the monitor survey may be individually characterized according to how it matches any one or more CMP(s) in the baseline survey.
  • flow may proceed from operation 225 (via 230) to operation 235 where the CMP hits of the monitor survey from operation 215, the associated attribute values determined in operation 225, and a user specified attribute parameter range or ranges (received in operation 234) all may be used to selectively bin the monitor survey data.
  • the monitor survey may be filtered such that only those CMPs whose attribute value is within the specified parameter range remain in the binned monitor survey data.
  • Figure 4 illustrates a sample filtered coverage visual representation with the filter being Dsrc + Drx less than or equal to 50 meters.
  • Figure 4 the x and y axes of Figure 4 are the same as Figure 3 , and represent the same area of the bin grid. Again, the shading of each rectangular region represents how "full" that particular bin is, with any CMPs whose Dsrc + Drx is greater than 50 meters being excluded. Hence, several bins that are 61% or more full in Figure 3 are significantly less full in Figure 4 after the filter has been applied.
  • Figures 3 and 4 only show coverage for the far-mid source receiver offsets, however, the operations 200 may be used to create a coverage plot (or attribute plot) over any desired offset range.
  • several different attributes may be combined, allowing the creation of several different filters or one compound filter. For example, a user may only want CMP hits with Dsrc + Drx being less than 50 meters and only hits that are from a certain receiver or streamer, or only hits from a certain offset or offset range.
  • a user may desire to create a filtered plot that only shows bins that have a certain percentage (e.g., 90%) of CMP hits with Dsrc + Drx less than 50 meters.
  • attribute value ranges may be used alone or in combination with one or more other attribute values filtering the data during binning.
  • any range of attribute values may be used, including by selecting only data above a certain range, below a certain range, or within an upper and lower bound.
  • flow may alternatively proceed from operation 225 (via 230) to operation 240.
  • operation 240 rather than generating a static, filtered 3D coverage of the monitor survey that can be viewed by a user (as is done in operation 235), the hybrid binning of operation 240 prepares data such that a user may later interactively filter coverage, or interactively create other visual representations based on data obtained during the monitor survey.
  • each matched CMP from the monitor survey may be binned together with its attribute value (as opposed to operation 235, where only CMPs falling within a user input range are binned without their attribute value).
  • the matched CMP hit may be binned with its attribute value by appending the attribute value to the file or files containing the information for each binned CMP.
  • the output of operation 240 may allow a user to either interactively filter coverage, as shown in operation 245, or to interactively filter attributes, as shown in operation 250, because, as mentioned, the attribute values upon which filtering will be done are stored with the binned CMPs for every source receiver pair in the monitor survey. Note that operation 240 may proceed generally from any survey and attribute generation, and is not limited to repeatability attributes.
  • a user may interactively input a desired attribute parameter range (for example, a maximum Dsrc + Drx value) and filter the 3D binned coverage by only including those CMPs whose attribute values are within the input parameter range.
  • a desired attribute parameter range for example, a maximum Dsrc + Drx value
  • a user may be able to interactively change the parameters, received in operation 244, and thereby specify how the data is combined.
  • a user may be able to change the filtering criteria. So, for example, a user may initially specify that Dsrc + Drx be less than 50 meters, and view the binned 3D coverage that has been filtered based on that specification.
  • a coverage plot may be created over any desired offset range (such as, for example, Dsrc + Drx being less than 50 meters but greater than 25 meters) and also, two or more user specified attributes may be combined, which may allow the creation of several different filters and/or one compound filter.
  • Figures 5A through 5C illustrate one possible set of images that may be output from operation 245.
  • Figure 5A shows an unfiltered coverage plot
  • Figure 5B shows a coverage plot filtered with Dsrc + Drx less than 50 meters
  • Figure 5C shows a coverage plot filtered with Dsrc + Drx less than 25 meters.
  • Figure 5 shows coverage plots for the far-mid offset range.
  • many other filters may be applied, such as, for example a filter of Dsrc + Drcx being less than 50 meters but greater than 25 meters.
  • a user can interactively input a desired offset range (received in operation 249), and view, instead of coverage, an attribute plot showing only attribute values from CMPs within the specified offset range.
  • a user may be able to specify that the plot show the minimum attribute value for each bin, the mean attribute value for each bin, the maximum attribute value for each bin, and so forth. If multiple attribute values were binned with the CMPs, a user may also be able to specify which attribute value to show for each bin.
  • Figure 6 illustrates an attribute plot showing the minimum value of Dsrc + Drx within each of the bins for the 1800-2000m offset range. This plot is generally similar to the plots illustrated in Figures 3, 4, and 5A through 5C , except instead of showing coverage (in percentage points), the plot in Figure 6 shows attribute values, in this case the minimum Dsrc + Drx value for each min specified in meters.
  • attribute data may be binned for monitor survey CMPs at a few offset steps (typically 3 or 4).
  • the output of this operation is a plot showing, for each bin, the value associated with that bin representing how "good" the repeatability was for one of the selected offset steps.
  • Repeatability attributes may include geometric values such as Dsrc (difference in source positions for the baseline and monitor surveys, which has radial, inline, and crossline components), Drx (difference in receiver positions for the baseline and monitor surveys, which also has radial, inline, and crossline components), Dcmp (difference in CMP positions for the baseline and monitor surveys), Doffset (difference in the radial offset from source to receiver for the baseline and monitor surveys), Dazimuth (difference in azimuth, the bearing from source to receiver), Dsrc + Drx (the combined difference of the source and receiver positions for the baseline and monitor surveys), Dcdp (the difference in the CDP positions for the baseline and monitor surveys, which has radial, inline, and crossline components), and so forth.
  • Dsrc difference in source positions for the baseline and monitor surveys, which has radial, inline, and crossline components
  • Drx difference in receiver positions for the baseline and monitor surveys, which also has radial, inline, and crossline components
  • Repeatability attributes also may include seismic values, such as static and dRMS (the difference in Root Mean Square amplitude between the monitor and baseline surveys) values, and so forth.
  • Other attributes that may be used for filtering coverage include seismic attributes (such as RMS, signal to noise ratio, and so forth) and acquisition attributes (such as water depth, streamer or receiver depth, feather, streamer separation, vessel speed, offset, azimuth, source ID, receiver ID, streamer ID, time, and so forth). If, for example, receiver ID is used as a filter, it may be possible to view the effect on coverage of removing one or more receivers interactively (such as while the survey is being conducted, or after the survey is complete while reviewing the data). Any of these attributes may be used alone or in combination with one another, in either filtered binning or hybrid binning, as described above.
  • the operations 200 described above may be implemented using various combinations of hardware and/or software.
  • Figure 7 illustrates an embodiment of the computer system 735 capable of storing and/or processing navigation and seismic data, such as to produce images according to the operations 200.
  • the computer system 735 may be a personal computer and/or a handheld electronic device aboard the vessel 101 (shown in Figures 1A and 1B ).
  • the computer system 735 may be an implementation of enterprise level computers, such as one or more blade-type servers within an enterprise in a land-based computer system.
  • a keyboard 740 and mouse 741 may be coupled to the computer system 735 via a system bus 748. The keyboard 740 and the mouse 741, in one example, may introduce user input to the computer system 735 and communicate that user input to a processor 743.
  • An input/output unit 749 (I/O) coupled to the system bus 748 represents such I/O elements as a printer, audio/video (A/V) I/O, etc.
  • Computer 735 also may include a video memory 744, a main memory 745 and a mass storage 742, all coupled to the system bus 748 along with the keyboard 740, the mouse 741 and the processor 743.
  • the mass storage 742 may include both fixed and removable media, such as magnetic, optical or magnetic optical storage systems and any other available mass storage technology.
  • the bus 748 may contain, for example, address lines for addressing the video memory 744 or the main memory 745.
  • the system bus 748 also may include a data bus for transferring data between and among the components, such as the processor 743, the main memory 745, the video memory 744 and the mass storage 742.
  • the video memory 744 may be a dual-ported video random access memory. One port of the video memory 744, in one example, is coupled to a video amplifier 746, which is used to drive a monitor 747.
  • the monitor 747 may be any type of monitor suitable for displaying graphic images, such as a cathode ray tube monitor (CRT), flat panel, or liquid crystal display (LCD) monitor or any other suitable data presentation device.
  • CTR cathode ray tube monitor
  • LCD liquid crystal display
  • the computer system includes a processor 743, which may be any suitable microprocessor or microcomputer.
  • the computer system 735 also may include a communication interface 750 coupled to the bus 748.
  • the communication interface 750 provides a two-way data communication coupling via a network link.
  • the communication interface 750 may be a satellite link, a local area network (LAN) card, a cable modem, and/or wireless interface.
  • the communication interface 750 sends and receives electrical, electromagnetic or optical signals that carry digital data streams representing various types of information.
  • Code received by the computer system 735 may be executed by the processor 743 as the code is received, and/or stored in the mass storage 742, or other non-volatile storage for later execution. In this manner, the computer system 735 may obtain program code in a variety of forms.
  • Program code may be embodied in any form of computer program product such as a medium configured to store or transport computer readable code or data, or in which computer readable code or data may be embedded. Examples of computer program products include CD-ROM discs, ROM cards, floppy disks, magnetic tapes, computer hard drives, servers on a network, and solid state memory devices.
  • the data processing system may execute operations that allow for the filtering using repeatability and other metrics.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (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)

Description

  • This Patent Cooperation Treaty patent application claims benefit under 35 U.S.C. § 119(e) to U.S. provisional application No. 61/466,595 , entitled "Method and apparatus for analyzing data in subsequent geophysical surveys" filed on March 23, 2011.
  • TECHNICAL FIELD
  • This disclosure relates generally to geophysical exploration systems, and more particularly to data acquisition and analysis in subsequent geophysical surveys.
  • BACKGROUND
  • Petrochemical products such as oil and gas, are ubiquitous in society and can be found in everything from gasoline to children's toys. Because of this, the demand for oil and gas remains high. In order to meet this high demand, it is important to locate oil and gas reserves in the Earth. Scientists and engineers conduct "surveys" utilizing, among other things, seismic and other wave exploration techniques to find oil and gas reservoirs within the Earth. These seismic exploration techniques often include controlling the emission of seismic energy into the Earth with a seismic source of energy (e.g., dynamite, air guns, vibrators, etc.), and monitoring the Earth's response to the seismic source with a receiver (e.g., a geophone, a hydrophone, etc.). By observing the reflected seismic signals detected by the receiver during the survey, the geophysical data pertaining to reflected signals may be acquired and these signals may indicate the composition of the Earth proximate the survey location.
  • As portions of the oil and gas are removed from the reservoir, a follow-up, or monitor survey may be taken near the same location as the original, baseline survey. Often-times, the data that the monitor survey provides is called 4D data, because it provides a fourth dimension to the data being collected-time. For example, this 4D survey may allow scientists and engineers to see how much oil or gas has been removed from the reservoir, how much oil or gas remains in the reservoir, and/or how the oil or gas has moved within the reservoir since the baseline survey. In cases where the sources and receivers for the baseline survey are installed on land or on the surface of the ocean floor and not removed between surveys, the monitor survey may involve firing the sources and recording the reflected seismic signals at the receivers. In other cases, such as where the surveys are conducted using streamers towed behind a vessel, or using transient sources and receivers on land, some system of quality control may be established to facilitate coordinating data sets between the two surveys.
  • "Repeatability" is a term used to refer to the measure of how closely the monitor survey correlates to the baseline survey. There are several different forms of repeatability. For example, positional or geometric repeatability is a measure of how positionally or geometrically close a monitor survey's data is to a baseline survey's data. Even if the geometric repeatability of a monitor survey is good (for example, where sources and receivers are installed permanently on land or on the ocean floor), other factors may decrease the repeatability, and therefore the reliability, of the data obtained in the monitor survey. A few aspects of a survey that may be of concern in terms of repeatability include: source position, receiver position, the common mid-point (CMP) position between the source and the receiver, offset or radial distance between the source and the receiver, azimuth, and components and combinations of these. In order to accurately determine the amount of oil or gas removed from a reservoir, scientists and engineers may need to know how closely the monitor survey correlates to one or more aspects of the baseline survey. A prior system for measuring the repeatability between surveys is disclosed in the International Patent Application published under number WO2010/054282 (CONOCOPHILLIPS COMPANY).
  • SUMMARY
  • A first aspect of the present invention relates to a method for processing geophysical data as defined in independent claim 1, and optionally any of dependent claims 2 to 13. The method includes matching a first set of data from a monitor survey with a second set of data from a baseline survey. It also includes generating an attribute value for each datum in the first set, storing each attribute value with each of the datum, and processing the first set of data based on the stored attribute values.
  • Another method of processing geophysical data may include correlating a first set of data from a first survey with a second set of data from a second survey. It also includes generating an attribute value for each datum in the first set, associating each generated attribute value with the datum from the first set of data, and processing the first set of data based on the attribute values.
  • A second aspect of the present invention relates to a tangible storage medium as defined in independent claim 14, and optionally dependent claim 15. The medium may store a plurality of instructions, including instructions that correlate a first set of data from a first survey with a second set of data from a second survey, instructions that generate an attribute value for each datum in the first set, instructions that associate each generated attribute value with the datum from the first set of data from which the attribute value was generated, and instructions that process the first set of data based on the attribute values.
  • Another tangible storage medium may store a filtered set of CMPs, wherein the filtered set of CMPs is filtered from a first set of CMPs, and the filtering is based on a set of attribute values.
  • BRIEF DESCRIPTION OF THE DRAWINGS
    • Figure 1A illustrates a top view of a vessel towing a seismic source and a plurality of seismic receivers positioned on streamers towed behind the vessel.
    • Figure 1B illustrates a side view of the vessel of Figure 1A.
    • Figures 2A through 2D illustrate a system and method of using one or more attribute value(s) as criteria for binning or filtering coverage.
    • Figure 3 is an unfiltered spatial coverage plot.
    • Figure 4 illustrates a spatial coverage plot filtered according to one embodiment.
    • Figures 5A through 5C illustrate spatial coverage plots that may be generated by interactively filtering coverage using one or more attribute value(s) as criteria.
    • Figure 6 illustrates a spatial attribute plot filtered according to another embodiment.
    • Figure 7 illustrates an embodiment of a computer system capable of storing and/or processing navigation and seismic data, such as to produce images according to the operations in Figs. 2A through 2D.
    DETAILED DESCRIPTION
  • Petrochemical products such as oil and gas, are ubiquitous in society and can be found in everything from gasoline to children's toys. Because of this, the demand for oil and gas remains high. In order to meet this high demand, it is important to locate oil and gas reserves in the Earth. Scientists and engineers conduct "surveys" utilizing, among other things, seismic and other wave exploration techniques to find oil and gas reservoirs within the Earth. These seismic exploration techniques often include controlling the emission of seismic energy into the Earth with a seismic source of energy (e.g., dynamite, air guns, vibrators, etc.), and monitoring the Earth's response to the seismic source with a receiver (e.g., a geophone, a hydrophone, etc.). By observing the reflected seismic signals detected by the receiver during the survey, the geophysical data pertaining to reflected signals may be acquired and these signals may indicate the composition of the Earth proximate the survey location.
  • As portions of the oil and gas are removed from the reservoir, a follow-up, or monitor survey may be taken near the same location as the original, baseline survey. Often-times, the data that the monitor survey provides is called 4D data, because it provides a fourth dimension to the data being collected-time. For example, this 4D survey may allow scientists and engineers to see how much oil or gas has been removed from the reservoir, how much oil or gas remains in the reservoir, and/or how the oil or gas has moved within the reservoir since the baseline survey. In cases where the sources and receivers for the baseline survey are installed on land or on the surface of the ocean floor and not removed between surveys, the monitor survey may involve firing the sources and recording the reflected seismic signals at the receivers. In other cases, such as where the surveys are conducted using streamers towed behind a vessel, or using transient sources and receivers on land, some system of quality control may be established to facilitate coordinating data sets between the two surveys.
  • "Repeatability" is a term used to refer to the measure of how closely the monitor survey correlates to the baseline survey. There are several different forms of repeatability. For example, positional or geometric repeatability is a measure of how positionally or geometrically close a monitor survey's data is to a baseline survey's data. Even if the geometric repeatability of a monitor survey is good (for example, where sources and receivers are installed permanently on land or on the ocean floor), other factors may decrease the repeatability, and therefore the reliability, of the data obtained in the monitor survey. A few aspects of a survey that may be of concern in terms of repeatability include: source position, receiver position, the common mid-point (CMP) position between the source and the receiver, offset or radial distance between the source and the receiver, azimuth, and components and combinations of these. In order to accurately determine the amount of oil or gas removed from a reservoir, scientists and engineers may need to know how closely the monitor survey correlates to one or more aspects of the baseline survey. Accordingly, systems and methods are disclosed that provide scientists and engineers the ability to analyze monitor survey geophysical data and monitor survey coverage that have been filtered based on one or more repeatability attributes.
  • Figure 1A shows a top view of a vessel 101 towing a source 102 and several receivers 103 on streamers behind the vessel 101. Figure 1B shows a side-view of the vessel 101 shown in Figure 1A with the source 102 and receivers 103 being towed behind the vessel 101 just beneath the surface of the water. For the sake of discussion, the embodiment depicted in Figures 1A and 1B illustrates the source and receiver being towed by the same vessel, however, other possible combinations are possible. For example, in other embodiments, either the source and/or receivers may be towed by separate vessels or may be implemented in land-based acquisition systems. In still other embodiments, the source and/or receivers may be stationary while the other is towed behind the vessel. In yet other embodiments, the receivers 103 may be positioned deeper in the water, for example, by using streamer steering devices, such as the DigiFIN® brand steering device available from ION Geophysical, Inc. In other embodiments, multiple sources may be used. Also, any type of source(s) or receiver(s) may be used, including for example, 1-, 2-, or 3-dimensional sources or receivers.
  • During operation, the source 102 may emit seismic energy (e.g., through an air gun), which may reflect off various portions of the Earth 104 and may be received back at the receivers 103 (as shown by the propagating seismic waves in Figure 1B). The signal received and processed at the receivers 103 may provide data that is useful in determining the composition of various portions of the Earth 104 proximate the location where the signal was reflected, which may include an oil and/or gas reservoir 105. If the amount of oil and/or gas in the reservoir 105 is depleted over time, then subsequent surveys conducted in substantially the same location as the first survey may indicate various properties of this depletion such as: decreasing pore pressures, migration of oil/water and/or gas/water contacts, drop in impedance, and so forth.
  • Figures 2A through 2D illustrate operations 200 that may utilize a predetermined metric or attribute as a criteria for binning coverage. Although any metric, or even multiple metrics, may be used as a filter to bin coverage, for the sake of discussion, the following disclosure illustrates the operations 200 in the context of using a repeatability metric as a filter. For example, some embodiments of operations 200 may include binning a CMP if a repeatability attribute for that CMP falls within a specified range. Also, the operations 200 may analyze data from the full offset range (i.e. data from every offset step, with an offset step being one or more receiver intervals) or data from a subset of the full offset range (such as only one offset step from within each offset zone, or only one offset zone within the full offset range).
  • Operations 200 will now be discussed in conjunction with the survey operations illustrated in Figures 1A and 1B. In Figure 2A, Operation 205 illustrates the sub-operations 206-210 associated with conducting a baseline survey. In general, during the baseline survey, a seismic source may fire shots or bursts of seismic energy into the Earth 104 in sub-operation 206. This energy reflects back and is received and measured as a signal at a receiver at a specific offset at sub-operation 207, with the offset being, in some embodiments, the radial range from the source to a specific receiver or group of receivers. In some embodiments, there may be several receivers at each offset and/or multiple offsets as indicated by boxes 106 and 108 in Figure 1A, in which each box contains several receivers and each box is located at a different offset.
  • For each source 102 - receiver 103 pair, operation 205 may include sub-operations 208 and 209 that include calculating and storing CMPs at each offset step. This is illustrated in Figure 1A with box 106 indicating a first offset where each of the receivers in the box 106 has a different CMP vis-à-vis the source 102. Sub-operation 209 may also store related context information with each CMP, including, for example acquisition and/or seismic attributes such as source ID, receiver ID, source position, receiver position, time, water depth, vessel speed, offset, azimuth, and so forth. Likewise, each of the receivers 103 in box 108 (which are at a different offset than the receivers in box 106), may have a different CMP. Operation 205 may in some embodiments also include sub-operation 210 where a grid of rectangles referred to as "bins" may be defined. Note that although the bins are illustrated in operation 205 as the final sub-operation, they may be defined at any point in the baseline survey. The width of a bin may be substantially the same as the spacing between CMP lines and the height of a bin may be substantially the same as the spacing between CMPs on a single CMP line. This is illustrated in Figure 1A with a bin 110 indicated by a dashed box. Of course, bins of other dimensions also may be defined with different heights or widths, depending on the resolution required.
  • After the baseline survey is complete, a visual representation of the baseline survey may be generated after "binning" each CMP by associating data from each source receiver pair with the bin where its CMP is located (this operation is not shown in Figs. 2A through 2D). A visual representation may be a spatial coverage plot (similar to the plot shown Fig. 3) or any other representation that conveys coverage information, such as charts, lists, graphs, and so forth. Within each bin, there may be CMPs for different receiver offset steps. If a bin has at least one CMP from each offset step, it may be referred to as "full fold." In other words, full fold coverage is achieved by obtaining a complete distribution of source to receiver offset steps over the entire offset range. Each bin may include, for example a count of how many CMPs are present at each offset step, or how "full" the fold is. Each bin may additionally or alternatively include a percentage of how full the fold is. In some cases, there may be more than one CMP per offset step, referred to as "duplicates," and these duplicates may be included or excluded during binning.
  • After a predetermined amount of time has elapsed, a monitor survey may be conducted, as shown in operation 215. In some embodiments, the amount of time that elapses between the baseline and monitor surveys is one or more years. The monitor survey in operation 215, along with its sub-operations 216-219, may be similar to the baseline survey: seismic shots may be fired in sub-operation 216, seismic energy may be reflected off the Earth and be received at one or more receiver(s) in sub-operation 217, the signal may be recorded, a CMP calculated for each source receiver pair in sub-operation 218, and the CMP for each offset may be stored in sub-operation 219. Then, similar to the binning and generating a visual representation of coverage described above, the baseline survey CMPs may be binned and a visual representation of the coverage of the monitor survey may be generated (this operation is not shown in Figs. 2A through 2D). Each bin may include, for example, a count of how many CMPs are present at each offset, or any other useful information. Figure 3 illustrates an unfiltered coverage plot for the monitor survey that may result from the sub-operations of operation 215. Referring to Figure 3, the x and y axes of Figure 3 represent a bin grid, as defined in operation 205. The shading in of each rectangular region represents how "full" that particular bin is, as shown in the key at the bottom of Figure 3. For example, bins that are white are 61% or more full. The shading scale can be configured to be any type of shading, including different colors and so forth. Likewise, the thresholds in the key and plot can be set to any percentage range; they may also be set to display the actual fold value rather than a percentage, or any other useful information.
  • Returning to Figure 2A, instead of generating an unfiltered coverage plot (like that shown in Figure 3), flow may proceed to operation 225 where data from the monitor survey may be matched to data from the baseline survey and one or more metrics may be generated. The sup-operations of operation 225 may include first, for each CMP of each shot fired in the monitor survey, determining a CMP from the baseline survey that is the "closest match" in sub-operation 226. In some embodiments, this match may be based upon the closest geometric match. One example of geometric closeness may be measured by the geometric repeatability value Dsrc + Drx. Dsrc is the difference in source positions between the baseline survey CMP and the monitor survey CMP at the time of firing the seismic source respectively. Likewise, Drx is the difference in the receiver positions between the baseline survey CMP and the monitor survey CMP at the time of firing the seismic source. In some embodiments, the positions of the source and/or receiver may be with reference to a navigation system aboard the vessel 101, such as by using a Global Positioning System (GPS). Dsrc + Drx represents the sum of these two differences. Evaluating Dsrc + Drx values in order to determine the closest CMP match may ensure that the same azimuth of data is obtained so as to be referencing substantially the same area of the subsurface, which may allow, for example, the characteristics of the reservoir 105 to be consistently characterized over time.
  • In the embodiments that utilize the repeatability metric Dsrc + Drx value to find a closest CMP match, all non-unique CMP hits from both monitor and baseline survey data may be analyzed to find a minimal Dsrc + Drx value for each bin. Some embodiments may perform this analysis at a single offset step while other embodiments may perform this analysis on multiple offset steps of interest. Furthermore, this analysis may involve exchanging the source and receiver positions in an effort to obtain different Dsrc + Drx values; this exchanging process is termed "reciprocity."
  • After each CMP in the monitor survey has been matched to a CMP in the baseline survey, the matched pairs of CMPs may be compared again in sub-operation 227 in order to generate an attribute value in sub-operation 228, and those attribute values may be stored with the CMP hits of the monitor survey in sub-operation 229. The attribute value may be called a repeatability metric, or a repeatability attribute, if it represents a comparison between a baseline survey CMP and a monitor survey CMP; repeatability referring to how well the monitor survey 'repeated' the baseline survey. Accordingly, each CMP in the monitor survey may be individually characterized according to how it matches any one or more CMP(s) in the baseline survey.
  • Next, as shown in Figure 2B, flow may proceed from operation 225 (via 230) to operation 235 where the CMP hits of the monitor survey from operation 215, the associated attribute values determined in operation 225, and a user specified attribute parameter range or ranges (received in operation 234) all may be used to selectively bin the monitor survey data. In this manner, the monitor survey may be filtered such that only those CMPs whose attribute value is within the specified parameter range remain in the binned monitor survey data. Figure 4 illustrates a sample filtered coverage visual representation with the filter being Dsrc + Drx less than or equal to 50 meters. Referring to Figure 4 in conjunction with Figure 3, the x and y axes of Figure 4 are the same as Figure 3, and represent the same area of the bin grid. Again, the shading of each rectangular region represents how "full" that particular bin is, with any CMPs whose Dsrc + Drx is greater than 50 meters being excluded. Hence, several bins that are 61% or more full in Figure 3 are significantly less full in Figure 4 after the filter has been applied.
  • It should be appreciated that that Figures 3 and 4 only show coverage for the far-mid source receiver offsets, however, the operations 200 may be used to create a coverage plot (or attribute plot) over any desired offset range. Also, as mentioned above, instead of just one user specified attribute parameter range, several different attributes may be combined, allowing the creation of several different filters or one compound filter. For example, a user may only want CMP hits with Dsrc + Drx being less than 50 meters and only hits that are from a certain receiver or streamer, or only hits from a certain offset or offset range. Alternatively, a user may desire to create a filtered plot that only shows bins that have a certain percentage (e.g., 90%) of CMP hits with Dsrc + Drx less than 50 meters. Any of several different attribute value ranges may be used alone or in combination with one or more other attribute values filtering the data during binning. Furthermore, any range of attribute values may be used, including by selecting only data above a certain range, below a certain range, or within an upper and lower bound.
  • Referring now to Figure 2C, flow may alternatively proceed from operation 225 (via 230) to operation 240. In operation 240, rather than generating a static, filtered 3D coverage of the monitor survey that can be viewed by a user (as is done in operation 235), the hybrid binning of operation 240 prepares data such that a user may later interactively filter coverage, or interactively create other visual representations based on data obtained during the monitor survey. In hybrid binning, each matched CMP from the monitor survey may be binned together with its attribute value (as opposed to operation 235, where only CMPs falling within a user input range are binned without their attribute value). The matched CMP hit may be binned with its attribute value by appending the attribute value to the file or files containing the information for each binned CMP. The output of operation 240 may allow a user to either interactively filter coverage, as shown in operation 245, or to interactively filter attributes, as shown in operation 250, because, as mentioned, the attribute values upon which filtering will be done are stored with the binned CMPs for every source receiver pair in the monitor survey. Note that operation 240 may proceed generally from any survey and attribute generation, and is not limited to repeatability attributes.
  • As described above, following hybrid binning in operation 240, flow may proceed to operation 245. Here, a user may interactively input a desired attribute parameter range (for example, a maximum Dsrc + Drx value) and filter the 3D binned coverage by only including those CMPs whose attribute values are within the input parameter range. A user may be able to interactively change the parameters, received in operation 244, and thereby specify how the data is combined. In other words, a user may be able to change the filtering criteria. So, for example, a user may initially specify that Dsrc + Drx be less than 50 meters, and view the binned 3D coverage that has been filtered based on that specification. Immediately following that, the user may be able to change the input to Dsrc + Drx being less than 25 meters. Without re-binning, the system 200 may be able to filter the coverage and visually display the binned 3D coverage filtered with the new input. Note that, as discussed above in connection with operation 235, a coverage plot may be created over any desired offset range (such as, for example, Dsrc + Drx being less than 50 meters but greater than 25 meters) and also, two or more user specified attributes may be combined, which may allow the creation of several different filters and/or one compound filter.
  • Figures 5A through 5C illustrate one possible set of images that may be output from operation 245. Specifically, Figure 5A shows an unfiltered coverage plot, Figure 5B shows a coverage plot filtered with Dsrc + Drx less than 50 meters, and Figure 5C shows a coverage plot filtered with Dsrc + Drx less than 25 meters. Note that Figure 5 shows coverage plots for the far-mid offset range. Also note that, as mentioned above, many other filters may be applied, such as, for example a filter of Dsrc + Drcx being less than 50 meters but greater than 25 meters.
  • Also, instead of proceeding to operation 245 following operation 240, flow may proceed to operation 250. There, a user can interactively input a desired offset range (received in operation 249), and view, instead of coverage, an attribute plot showing only attribute values from CMPs within the specified offset range. A user may be able to specify that the plot show the minimum attribute value for each bin, the mean attribute value for each bin, the maximum attribute value for each bin, and so forth. If multiple attribute values were binned with the CMPs, a user may also be able to specify which attribute value to show for each bin. Figure 6 illustrates an attribute plot showing the minimum value of Dsrc + Drx within each of the bins for the 1800-2000m offset range. This plot is generally similar to the plots illustrated in Figures 3, 4, and 5A through 5C, except instead of showing coverage (in percentage points), the plot in Figure 6 shows attribute values, in this case the minimum Dsrc + Drx value for each min specified in meters.
  • Referring now to Figure 2D, flow may alternately proceed from operation 225 (via 230) to operation 255, attribute binning. Here, attribute data may be binned for monitor survey CMPs at a few offset steps (typically 3 or 4). The output of this operation is a plot showing, for each bin, the value associated with that bin representing how "good" the repeatability was for one of the selected offset steps.
  • While the above discussion has specifically mentioned the use of repeatability attributes such as Dsrc + Drx for filtering, numerous CMP-based or shot-based attributes may be used for filtering including repeatability attributes, acquisition attributes, seismic attributes, and so forth. Repeatability attributes may include geometric values such as Dsrc (difference in source positions for the baseline and monitor surveys, which has radial, inline, and crossline components), Drx (difference in receiver positions for the baseline and monitor surveys, which also has radial, inline, and crossline components), Dcmp (difference in CMP positions for the baseline and monitor surveys), Doffset (difference in the radial offset from source to receiver for the baseline and monitor surveys), Dazimuth (difference in azimuth, the bearing from source to receiver), Dsrc + Drx (the combined difference of the source and receiver positions for the baseline and monitor surveys), Dcdp (the difference in the CDP positions for the baseline and monitor surveys, which has radial, inline, and crossline components), and so forth. Repeatability attributes also may include seismic values, such as static and dRMS (the difference in Root Mean Square amplitude between the monitor and baseline surveys) values, and so forth. Other attributes that may be used for filtering coverage include seismic attributes (such as RMS, signal to noise ratio, and so forth) and acquisition attributes (such as water depth, streamer or receiver depth, feather, streamer separation, vessel speed, offset, azimuth, source ID, receiver ID, streamer ID, time, and so forth). If, for example, receiver ID is used as a filter, it may be possible to view the effect on coverage of removing one or more receivers interactively (such as while the survey is being conducted, or after the survey is complete while reviewing the data). Any of these attributes may be used alone or in combination with one another, in either filtered binning or hybrid binning, as described above. The operations 200 described above may be implemented using various combinations of hardware and/or software.
  • Figure 7 illustrates an embodiment of the computer system 735 capable of storing and/or processing navigation and seismic data, such as to produce images according to the operations 200. In some embodiments, the computer system 735 may be a personal computer and/or a handheld electronic device aboard the vessel 101 (shown in Figures 1A and 1B). In other embodiments, the computer system 735 may be an implementation of enterprise level computers, such as one or more blade-type servers within an enterprise in a land-based computer system. A keyboard 740 and mouse 741 may be coupled to the computer system 735 via a system bus 748. The keyboard 740 and the mouse 741, in one example, may introduce user input to the computer system 735 and communicate that user input to a processor 743. Other suitable input devices may be used in addition to, or in place of, the mouse 741 and the keyboard 740. An input/output unit 749 (I/O) coupled to the system bus 748 represents such I/O elements as a printer, audio/video (A/V) I/O, etc.
  • Computer 735 also may include a video memory 744, a main memory 745 and a mass storage 742, all coupled to the system bus 748 along with the keyboard 740, the mouse 741 and the processor 743. The mass storage 742 may include both fixed and removable media, such as magnetic, optical or magnetic optical storage systems and any other available mass storage technology. The bus 748 may contain, for example, address lines for addressing the video memory 744 or the main memory 745.
  • The system bus 748 also may include a data bus for transferring data between and among the components, such as the processor 743, the main memory 745, the video memory 744 and the mass storage 742. The video memory 744 may be a dual-ported video random access memory. One port of the video memory 744, in one example, is coupled to a video amplifier 746, which is used to drive a monitor 747. The monitor 747 may be any type of monitor suitable for displaying graphic images, such as a cathode ray tube monitor (CRT), flat panel, or liquid crystal display (LCD) monitor or any other suitable data presentation device.
  • The computer system includes a processor 743, which may be any suitable microprocessor or microcomputer. The computer system 735 also may include a communication interface 750 coupled to the bus 748. The communication interface 750 provides a two-way data communication coupling via a network link. For example, the communication interface 750 may be a satellite link, a local area network (LAN) card, a cable modem, and/or wireless interface. In any such implementation, the communication interface 750 sends and receives electrical, electromagnetic or optical signals that carry digital data streams representing various types of information.
  • Code received by the computer system 735 may be executed by the processor 743 as the code is received, and/or stored in the mass storage 742, or other non-volatile storage for later execution. In this manner, the computer system 735 may obtain program code in a variety of forms. Program code may be embodied in any form of computer program product such as a medium configured to store or transport computer readable code or data, or in which computer readable code or data may be embedded. Examples of computer program products include CD-ROM discs, ROM cards, floppy disks, magnetic tapes, computer hard drives, servers on a network, and solid state memory devices. Regardless of the actual implementation of the computer system 735, the data processing system may execute operations that allow for the filtering using repeatability and other metrics.

Claims (15)

  1. A method of processing geophysical data, the method comprising the acts of:
    matching a first set of data from a monitor survey with a second set of data from a baseline survey (227), wherein the first and second sets of data comprise common mid-point positions, CMPs, between sources (102) and receivers (103);
    generating an attribute value for each datum in the first set (228), wherein the attribute value comprises a repeatability metric;
    storing each attribute value with each of said datum (229); and
    processing the first set of data based on the stored attribute values (230);
    characterised in that
    said repeatability metric is Dsrc + Drx, wherein Dsrc is a difference in source positions between a baseline survey CMP and a monitor survey CMP and Drx is a difference in receiver positions between the baseline survey CMP and the monitor survey CMP, each of which has radial, inline and crossline components.
  2. The method of claim 1, wherein the acts of matching and generating occur at substantially the same time.
  3. The method of claim 1, wherein the act of processing comprises:
    filtering the first set of data by removing data whose associated attribute lies outside a range of attribute values.
  4. The method of claim 1, wherein the act of processing comprises:
    filtering the first set of data by including data whose associated attribute lies within a range of attribute values;
    and, optionally, further comprises:
    binning, such that each datum within the filtered first set of data is associated with at least one of a plurality of bins.
  5. The method of claim 4, wherein each of the plurality of bins corresponds to a predefined boundary within a spatial grid.
  6. The method of claim 4, wherein the acts of filtering and binning happen at substantially the same time.
  7. The method of claim 4, further comprising the act of generating a visual representation of the binned data.
  8. The method of claim 1, wherein the monitor survey occurs after the baseline survey; and, optionally, wherein the geophysical data represents petrochemical deposits in the Earth.
  9. The method of claim 1, wherein the attribute value comprises a further repeatability metric selected from the group consisting of Dcmp, Doffset, Dazimuth, Dcdp, and dRMS.
  10. The method of claim 1, wherein the attribute value comprises a seismic attribute.
  11. The method of claim 1, wherein the attribute value comprises an acquisition attribute; and, optionally, wherein the acquisition attribute is a receiver number within an array of receivers.
  12. The method of claim 1, wherein the first and second sets of data also comprise context information for each CMP;
    and, optionally, wherein the context information comprises a source position, a receiver position, a receiver ID, a time, and an offset.
  13. The method of claim 1, wherein the difference between the first and second sets of data represents a depletion in petrochemical content over a period of time.
  14. A computer readable tangible storage medium storing a plurality of instructions that, when executed by a computer, cause the computer to:
    correlate a first set of data from a first survey with a second set of data from a second survey, wherein the first survey is a monitor survey and the second survey is a baseline survey and wherein the first and second sets of data comprise common mid-point positions, CMPs;
    generate an attribute value for each datum in the first set, wherein the attribute value comprises a repeatability metric;
    associate each generated attribute value with the datum from the first set of data from which the attribute value was generated;
    process the first set of data based on the attribute values; and
    store said attribute value with said datum from the first set of data;
    characterised in that
    said repeatability metric is Dsrc + Drx, wherein Dsrc is a difference in source positions between a baseline survey CMP and a monitor survey CMP and Drx is a difference in receiver positions between the baseline survey CMP and the monitor survey CMP, each of which has radial, inline and crossline components.
  15. The computer readable tangible storage medium of claim 14, wherein the processing instructions, when executed by the computer, further cause the computer to:
    filter out data from the first set whose associated attribute value lies outside a range of attribute values,
    optionally, wherein the filtering comprises: binning, such that each datum within the filtered first set of data is associated with at least one of a plurality of bins;
    further optionally, wherein each of the plurality of bins corresponds to a predefined boundary within a spatial grid;
    or further wherein the computer is further caused to generate a visual representation of the binned data, optionally wherein the visual representation is a coverage plot;
    and/or wherein the monitor survey occurs after the baseline survey;
    and/or wherein the first and second sets of data represent petrochemical deposits in the Earth;
    and/or wherein the attribute value comprises a further repeatability metric selected from the group consisting of Dcmp, Doffset, Dazimuth, Dcdp, and dRMS;
    and/or wherein the attribute value comprises a seismic attribute;
    and/or wherein the attribute value comprises an acquisition attribute;
    and/or wherein the first and second sets of data further comprise context information for all CMPs;
    and/or wherein the difference between the first and second sets of data represents a depletion in petrochemical content over a period of time.
EP11741397.1A 2011-03-23 2011-08-01 Method and apparatus for analyzing data in time-lapse geophysical surveys Active EP2689275B1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201161466595P 2011-03-23 2011-03-23
PCT/US2011/046086 WO2012138369A1 (en) 2011-03-23 2011-08-01 Method and apparatus for analyzing data in time -lapse geophysical surveys

Publications (2)

Publication Number Publication Date
EP2689275A1 EP2689275A1 (en) 2014-01-29
EP2689275B1 true EP2689275B1 (en) 2020-10-07

Family

ID=44629835

Family Applications (1)

Application Number Title Priority Date Filing Date
EP11741397.1A Active EP2689275B1 (en) 2011-03-23 2011-08-01 Method and apparatus for analyzing data in time-lapse geophysical surveys

Country Status (4)

Country Link
US (1) US10379240B2 (en)
EP (1) EP2689275B1 (en)
CN (1) CN103534614B (en)
WO (1) WO2012138369A1 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160161620A1 (en) * 2013-06-28 2016-06-09 Cgg Services Sa System and method for estimating repeatability using base data
CN105549083B (en) * 2015-12-25 2018-03-13 中国石油天然气集团公司 A kind of method and apparatus for determining seismic observation system repeatability and integrally measuring

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6438069B1 (en) 1996-09-13 2002-08-20 Pgs Data Processing, Inc. Method for time lapse reservoir monitoring
US5873051A (en) 1996-09-20 1999-02-16 Geoquest, A Division Of Schlumberger Technology Corporation Method and apparatus for processing at least two seismic data sets during a step to derive a third data set
US6831877B2 (en) 1999-11-08 2004-12-14 Vecta Technologies, Inc. System and method for determining an azimuth of a seismic energy source
US6925387B2 (en) 2003-08-14 2005-08-02 Westerngeco L.L.C. Method and apparatus for kinematically linking multiple seismic domains
US7027354B2 (en) * 2003-09-22 2006-04-11 4Th Wave Imaging Corp. Method of obtaining pore pressure and fluid saturation changes in subterranean reservoirs by forward modeling
US7092824B2 (en) 2003-10-20 2006-08-15 Ascend Geo Llp Methods and systems for interactive investigation of geophysical data
GB2420408B (en) 2004-11-19 2008-03-12 Geophysique Cie Gle Method for processing at least two sets of seismic data
EP1889200A4 (en) 2005-06-09 2017-08-16 Exxonmobil Upstream Research Company Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys
US7974153B2 (en) 2006-09-29 2011-07-05 Geokinetics Acquisition Company Inc. Three-dimensional seismic survey methods using a perturbation pattern to provide bin fractionation
US7636275B2 (en) 2007-02-06 2009-12-22 Conocophillips Company Direct time lapse inversion of seismic data
US8488409B2 (en) * 2007-05-17 2013-07-16 Westerngeco L.L.C. Acquiring azimuth rich seismic data in the marine environment using a regular sparse pattern of continuously curved sail lines
WO2009109663A2 (en) 2008-03-07 2009-09-11 Shell Internationale Research Maatschappij B.V. Method of marine time-lapse seismic surveying
US8339898B2 (en) * 2008-05-25 2012-12-25 Westerngeco L.L.C. Processing seismic data using combined regularization and 4D binning
CN101614826B (en) 2008-06-26 2012-01-04 王振华 Method and device for realizing binning homogenization in three-dimensional seismic data processing
US9207347B2 (en) 2008-07-05 2015-12-08 Westerngeco L.L.C. Determining the structure of a towed seismic spread element
US8717846B2 (en) 2008-11-10 2014-05-06 Conocophillips Company 4D seismic signal analysis

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
None *

Also Published As

Publication number Publication date
US10379240B2 (en) 2019-08-13
CN103534614B (en) 2017-08-04
WO2012138369A9 (en) 2020-05-07
EP2689275A1 (en) 2014-01-29
CN103534614A (en) 2014-01-22
WO2012138369A1 (en) 2012-10-11
US20140012508A1 (en) 2014-01-09

Similar Documents

Publication Publication Date Title
CA2712465C (en) Method for determining geological information related to a subsurface volume of interest
US20090132170A1 (en) Seismic data processing and visualization
US7917295B2 (en) Modeling and filtering coherent noise in seismic surveying
CN102147481B (en) The correction based on inclination angle of the data reconstruction in three-dimensional surface related multiple is predicted
EP2047300A2 (en) Time lapse marine seismic surveying employing interpolated multicomponent streamer pressure data
EP1879052A2 (en) Time lapse marine seismic surveying employing interpolated multicomponent streamer pressure data
AU2013214831B2 (en) Method and apparatus for processing seismic data
US20170219731A1 (en) Method of visualizing and interpreting wide azimuth profile (wap)
EP2249183B1 (en) Modeling and filtering coherent noise in seismic surveying
EP2689275B1 (en) Method and apparatus for analyzing data in time-lapse geophysical surveys
EP2689274B1 (en) Method and apparatus for analyzing data in subsequent geophysical surveys
Calvert et al. Preserving azimuthal velocity information: Experiences with cross-spread noise attenuation and offset vector tile preSTM
Khakimov et al. Using active and passive near-field hydrophones to image the near-surface in ultra-shallow waters offshore Abu Dhabi
US20080232195A1 (en) Apparatus and Method for Processing Geophysical Information
Al Ramadhan et al. Processing frequent 4D land seismic data with buried sensors for CO2 monitoring
Babkina et al. Utilizing Near-Field Hydrophones for Near-Surface Imaging Over Shallow Salt Structure Offshore Abu Dhabi
CA2888252A1 (en) Method of enhancing flat spots in three-dimensional seismic interpretation
Brink Result-oriented marine seismic acquisition
Buia et al. Coil Shooting on Tulip Discovery: Seismic Processing Challenges, Opportunities and Results

Legal Events

Date Code Title Description
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

17P Request for examination filed

Effective date: 20131022

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

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: EXAMINATION IS IN PROGRESS

17Q First examination report despatched

Effective date: 20170511

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

INTG Intention to grant announced

Effective date: 20200610

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

AK Designated contracting states

Kind code of ref document: B1

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

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

Ref country code: AT

Ref legal event code: REF

Ref document number: 1321757

Country of ref document: AT

Kind code of ref document: T

Effective date: 20201015

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602011068838

Country of ref document: DE

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: NL

Ref legal event code: MP

Effective date: 20201007

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 1321757

Country of ref document: AT

Kind code of ref document: T

Effective date: 20201007

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210108

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210208

Ref country code: NL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: NO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210107

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG4D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210107

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210207

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602011068838

Country of ref document: DE

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

26N No opposition filed

Effective date: 20210708

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

REG Reference to a national code

Ref country code: DE

Ref legal event code: R119

Ref document number: 602011068838

Country of ref document: DE

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

REG Reference to a national code

Ref country code: BE

Ref legal event code: MM

Effective date: 20210831

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20210831

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20210831

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210207

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20210801

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20210801

Ref country code: FR

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20210831

Ref country code: DE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20220301

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20210831

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20110801

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20201007

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20240826

Year of fee payment: 14