US20160349392A1 - Using seabed sensors and sea-surface reflections for structural imaging of a subsurface location in a geological formation - Google Patents
Using seabed sensors and sea-surface reflections for structural imaging of a subsurface location in a geological formation Download PDFInfo
- Publication number
- US20160349392A1 US20160349392A1 US15/166,902 US201615166902A US2016349392A1 US 20160349392 A1 US20160349392 A1 US 20160349392A1 US 201615166902 A US201615166902 A US 201615166902A US 2016349392 A1 US2016349392 A1 US 2016349392A1
- Authority
- US
- United States
- Prior art keywords
- wavefield
- location
- shot
- sensors
- seismic
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 38
- 230000015572 biosynthetic process Effects 0.000 title description 8
- 238000000034 method Methods 0.000 claims abstract description 62
- 230000008569 process Effects 0.000 claims abstract description 22
- 238000004519 manufacturing process Methods 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 10
- 238000013500 data storage Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 5
- 238000012544 monitoring process Methods 0.000 claims description 4
- 239000002245 particle Substances 0.000 claims description 3
- 230000009977 dual effect Effects 0.000 claims 1
- 238000004458 analytical method Methods 0.000 description 10
- 238000005755 formation reaction Methods 0.000 description 7
- 230000004044 response Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000003672 processing method Methods 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical group 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3843—Deployment of seismic devices, e.g. of streamers
- G01V1/3852—Deployment of seismic devices, e.g. of streamers to the seabed
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1427—Sea bed
-
- 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/612—Previously recorded data, e.g. time-lapse or 4D
- G01V2210/6122—Tracking reservoir changes over time, e.g. due to production
Definitions
- Geophysical surveys are often used for oil and gas exploration in geological formations, which may be located below marine environments.
- Various types of signal sources and geophysical sensors may be used in different types of geophysical surveys.
- Seismic geophysical surveys are based on the use of acoustic waves.
- a survey vessel may tow an acoustic source (e.g., an air gun or a marine vibrator) and a plurality of streamers along which a number of acoustic sensors (e.g., hydrophones and/or geophones) are located.
- acoustic sensors may be located on the seabed rather than being towed.
- Acoustic waves generated by the survey source may be transmitted to the Earth's crust and then reflected back and captured at the towed and/or seabed geophysical sensors.
- Data collected during a marine geophysical survey may be analyzed to locate hydrocarbon-bearing geological structures, and thus determine where deposits of oil and natural gas may be located.
- FIG. 1 is a diagram illustrating a geophysical survey system in accordance with an embodiment of the invention.
- FIG. 2 depicts an exemplary array of seabed sensors in accordance with an embodiment of the invention.
- FIG. 3 is a flow chart of a method using seabed sensors and sea surface reflections for structural imaging of a subsurface location in a geological formation in accordance with an embodiment of the invention.
- FIG. 4A depicts shot locations of a first seismic survey (A) in accordance with an embodiment of the invention.
- FIG. 4B depicts shot locations of a second seismic survey (B) in accordance with an embodiment of the invention.
- FIG. 4C depicts wavefields from the first and second seismic surveys (A and B) that are utilized in combination to image a subsurface location (P) in accordance with an embodiment of the invention.
- FIG. 4D depicts virtual wavefields after reversal of source and sensor definitions of FIG. 4C in accordance with an embodiment of the invention.
- FIG. 4E depicts the virtual wavefield utilized for imaging the subsurface location (P) after deconvolution or cross-correlation is applied to the virtual wavefields of FIG. 4D in accordance with an embodiment of the invention.
- FIG. 5A depicts shot locations of a first seismic survey (A) in accordance with an embodiment of the invention.
- FIG. 5B depicts shot locations of a second seismic survey (B) in the presence of a surface obstruction in accordance with an embodiment of the invention.
- FIG. 5C depicts wavefields from the first and second seismic surveys (A and B) that are utilized in combination to image a subsurface location (P) that may be underneath the obstruction in accordance with an embodiment of the invention.
- FIG. 5D depicts virtual wavefields after reversal of source and sensor definitions of FIG. 5C in accordance with an embodiment of the invention.
- FIG. 5E depicts the virtual wavefield utilized for imaging the subsurface location (P) after deconvolution or cross-correlation is applied to the virtual wavefields of FIG. 5D in accordance with an embodiment of the invention.
- FIG. 6A depicts shot locations of a first seismic survey (A) in accordance with an embodiment of the invention.
- FIG. 6B depicts shot locations of a second seismic survey (B) with different shot intervals than the first seismic survey (A) in accordance with an embodiment of the invention.
- FIG. 6C depicts wavefields from the first and second seismic surveys (A and B) that are utilized in combination to image a subsurface location (P) in accordance with an embodiment of the invention.
- FIG. 6D depicts virtual wavefields after reversal of source and sensor definitions of FIG. 6C in accordance with an embodiment of the invention.
- FIG. 6E depicts the virtual wavefield utilized for imaging the subsurface location (P) after deconvolution or cross-correlation is applied to the virtual wavefields of FIG. 6D in accordance with an embodiment of the invention.
- FIG. 7 is a flow chart of a method using seabed sensors and sea surface reflections for reservoir production monitoring in accordance with an embodiment of the invention.
- FIG. 8A depicts wavefields (A 0 and A 1 ) from a first seismic survey (A) at a first time that may be utilized in combination to image a subsurface location (P) at the first time in accordance with an embodiment of the invention.
- FIG. 8B depicts virtual wavefields after reversal of source and sensor definitions of FIG. 8A in accordance with an embodiment of the invention.
- FIG. 8C depicts the virtual wavefield utilized for imaging the subsurface location (P) at the first time in accordance with an embodiment of the invention.
- FIG. 8D depicts a wavefield (B 1 ) from a second seismic survey (B) at a second time that is utilized in combination with a wavefield (A 0 ) from the first seismic survey (A) to image the subsurface location (P) at the second time in accordance with an embodiment of the invention.
- FIG. 8E depicts virtual wavefields after reversal of source and sensor definitions of FIG. 8D in accordance with an embodiment of the invention.
- FIG. 8F depicts the virtual wavefield utilized for imaging the subsurface location (P) at the second time in accordance with an embodiment of the invention.
- FIG. 8G depicts a wavefield (C 1 ) from a third seismic survey (C) at a third time that is utilized in combination with a wavefield (A 0 ) from the first seismic survey (A) to image the subsurface location (P) at the third time in accordance with an embodiment of the invention.
- FIG. 8H depicts virtual wavefields after reversal of source and sensor definitions of FIG. 8G in accordance with an embodiment of the invention.
- FIG. 8I depicts the virtual wavefield utilized for imaging the subsurface location (P) at the third time in accordance with an embodiment of the invention.
- FIG. 9 shows a simplified example of a computer system which may be used in performing processing steps in accordance with an embodiment of the invention.
- Three-dimensional (3D) seismic surveying generally uses one or more acoustic (seismic) sources and an array of acoustic (seismic) sensors.
- a pattern of shots is performed by the source(s).
- each sensor in the array senses and records a seismic wavefield.
- the seismic data obtained includes multiple seismic wavefields for each of multiple shots.
- the seismic data may be processed to create images of a subsurface volume underneath the array of sensors.
- a 3D seismic survey design includes the positions of the shots and the locations of the sensors.
- the locations of the sensors are fixed.
- the positions of the shots may be designed to obtain seismic data covering the desired subsurface volume.
- Four-dimensional (4D) seismic surveying involves surveying the same geological formation at different points in time. For example, 4D surveying may be used to determine changes over time in fluid saturation, pressure, and so on. Because seabed sensors are fixed-position receivers that do not change position during the seismic acquisition, they may be redeployed at the same positions at a later time so as to be advantageously used in 4D marine seismic surveying to reduce variations between surveys.
- One problem or difficulty with marine seismic surveys may be caused by an obstruction that prevents positioning shots and/or sensors above a targeted sub-surface region.
- Such an obstruction makes it difficult to obtain seismic images of the sub-surface region under the obstruction in a conventional 3D or 4D seismic survey.
- a production platform may be located above the targeted sub-surface region.
- the production platform forms an obstruction that blocks shots from being taken at sea-surface positions in the area above the targeted sub-surface region.
- the production platform may also block sensors from being towed above the targeted region.
- the present disclosure presents an innovative technique to generate seismic images by combining seismic data from different surveys.
- the technique has various practical applications, such as imaging under an obstruction and efficient reservoir monitoring.
- FIG. 1 a diagram illustrating one embodiment of a Geophysical survey system 100 is shown.
- the system 100 includes a survey vessel 10 , signal source 32 , and geophysical sensors 22 .
- the survey vessel 10 may be configured to move along the surface 11 of body of water 15 such as a lake or the ocean, as shown.
- signal source 32 is configured to transmit an acoustic signal.
- the illustrated dashed line shows a direction of travel for a portion of an acoustic wave-front generated by signal source 32 .
- the acoustic signal proceeds through the sea floor 25 and is reflected by reservoir 120 , the sea surface (an air/water interface) 11 , and a reflector 27 located above reservoir 120 .
- the acoustic signal may be reflected from the sea surface multiple times (and reflected from the seabed in between the sea surface reflections).
- Survey vessel 10 may include equipment, shown generally at 12 and for convenience collectively referred to as “survey equipment.”
- Survey equipment 12 may include devices such as a data recording unit for making a record of signals generated by various geophysical sensors in the system 100 .
- Survey equipment 12 may also include navigation equipment, which may be configured to control, determine, and record, at selected times, the geodetic positions of: survey vessel 10 , signal source 32 , and/or other survey elements.
- Geodetic position may be determined using various devices, including global navigation satellite systems such as the global positioning system (GPS), for example.
- GPS global positioning system
- vessel 10 is configured to tow one or more seismic streamers (not shown).
- survey vessel 10 tows a single signal source 32 .
- survey vessel 10 may tow any appropriate number of signal sources, including as few as none or as many as six or more.
- the location of the signal sources may be centered behind survey vessel 10 or displaced from the center line, and may be at various distances relative to survey vessel 10 , including attached to the hull.
- Signal sources may be any type of signal source known in the art.
- Signal source 32 may include an array of multiple signal sources.
- signal source 32 may include a plurality of seismic sources (e.g., air guns).
- the term “signal source” may refer to a single signal source or to an array of signal sources.
- signal sources 32 are each coupled to survey vessel 10 at one end through winch 115 or a similar spooling device that enables changing the deployed length of a signal source cable.
- Survey equipment 12 may include signal source control equipment (not shown separately) for selectively operating and maneuvering signal source 32 .
- the specific configuration depicted in FIG. 1 is an example only. Other configurations may be used for the geophysical survey system.
- Geophysical sensors 22 may be referred to as seabed sensors and may be any type of geophysical sensor known in the art.
- the geophysical sensors may include pressure sensors (e.g., hydrophones) that measure pressure (P) as a scalar and particle-velocity sensors (e.g. geophones and accelerometers) that measure particle velocity (Z) as a vector.
- the pressure and particle-velocity sensors may be both included in a given seabed sensor.
- Other examples of geophysical sensors that may be used include pressure-time-gradient responsive seismic sensors, electrodes, magnetometers, temperature sensors, fiber optic sensors, or combinations of the foregoing.
- geophysical sensors 22 may measure, for example, seismic energy indicative of the response of various structures in the Earth's subsurface formation below the bottom 25 of body of water 15 to energy imparted into the subsurface formation by one or more of signal source 32 .
- Seismic energy may originate from signal source 32 , or an array of such signal sources, deployed in body of water 15 and towed by survey vessel 10 .
- geophysical sensors 22 are configured to store data until they are physically retrieved, at which point the stored data may be accessed.
- geophysical sensors 22 may be configured to communicate with survey equipment onshore or in survey vessel wirelessly or using transmission cables, in various embodiments.
- Geophysical sensors 22 may include relatively stable and/or synchronized clocks to timestamp recorded data.
- Geophysical sensors 22 may be comprised in ocean bottom cables (OBC) and linked together or may be ocean bottom nodes (OBN) and be individual recording systems (e.g., with a clock, battery, sensor, and data/time recorder integrated within each unit). In the illustrated embodiment, two of geophysical sensors 22 are located distance dl apart.
- Survey equipment 12 includes a computing system configured to, inter alia, process data from geophysical sensors 22 .
- a computing system at another location may perform geophysical analysis on data gathered by geophysical survey system 100 (e.g., on land after a survey has been conducted).
- the computing system may include or be configured to access a non-transitory storage medium having instructions stored thereon that are executable to cause the system to perform various operations described herein.
- the computing system may include one or more processors configured to execute the program instructions and one or more data storage devices for storing results in a tangible medium.
- the pattern of an array of seabed sensors for a survey may be determined in advance in order to generate data usable to perform the disclosed geophysical analysis techniques.
- the pattern may specify distances between sensors, sensor orientation, shot area, and other parameters.
- FIG. 2 depicts an exemplary array of seabed sensors in accordance with an embodiment of the invention.
- this exemplary array includes 380 sensors in a “raster-shaped” layout pattern.
- the numbers depicted indicate the seabed positions of sensors #1, #50, #100, #165, #235, #305 and #380.
- This layout pattern may be particularly useful in separated wavefield imaging by multiples (SWIM).
- FIG. 3 is a flow chart of a method 300 using seabed sensors and sea surface reflections for structural imaging of a subsurface location in a geological formation in accordance with an embodiment of the invention.
- the method 300 combines data from multiple surveys. The steps of FIG. 3 are discussed below in relation to the exemplary case illustrated in FIGS. 4A-4E .
- FIGS. 4A-4E illustrate the method steps of FIG. 3 used to acquire seismic data in multiple surveys (A and B) and process the acquired data in combination to image a subsurface location P in a subsurface region 29 .
- the subsurface region 29 may be difficult or impossible to image using the data from either survey alone due to the geometries of the shot locations relative to the sensor positions. Nevertheless, the subsurface region 29 may be imaged in three-dimensions (3D imaging) using the method described below for multiple subsurface points in the subsurface region 29 .
- a first marine seismic survey (survey A) may be performed.
- the performance of survey A involves generating seismic signals at a first set of positions and recording a first set of received signals using an array of seabed sensors.
- FIG. 4A An illustrative example of the first survey (survey A) is shown in FIG. 4A . Shown in FIG. 4A are: shot locations of survey A beneath the sea surface 11 , and seabed sensors 22 on the sea floor 25 .
- a second marine seismic survey (survey B) may be performed.
- the performance of survey B involves generating seismic signals at a second set of positions and recording a second set of received signals using the same array (location) of seabed sensors.
- FIG. 4B An illustrative example of the second survey (survey B) is shown in FIG. 4B . Shown in FIG. 4B are: shot locations of survey B beneath the sea surface 11 , and the array of seabed sensors 22 on the sea floor 25 . Note that, in accordance with the illustrated embodiment, the locations of the seabed sensors location 22 do not change between surveys A and B.
- up-going and down-going wavefields are separated for each survey. In an exemplary implementation, this may be accomplished by computing the summation and difference between the pressure (P) and particle velocity (Z) fields. Because the up-going wavefield is in-phase for both pressure and particle-velocity records, the up-going wavefield may be extracted by summing the two measurements (PZ-summation) after calibration. The down-going wavefield is canceled in the PZ-summation because it is out-of-phase on the P field compared to the Z field. The down-going wavefield may be extracted by determining the difference between the P field and the calibrated Z field.
- wavefield A from survey A and wavefield B from survey B are selected. Each selected wavefield originates at a shot location and ends at a sensor.
- the selected wavefield B originates at shot location B k from survey B and travels downward.
- Wavefield B reflects from subsurface point P in the subsurface region 29 so as to travel upward, then reflects from point r B of the sea surface 11 so as to travel downward.
- Wavefield B may then reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 again to travel downward until it is detected by a seabed sensor at seabed position S i .
- the selected wavefield A originates at shot location A j and travels downward.
- Wavefield A may reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position S i as that which detected wavefield B.
- the first sea-surface reflection point r B of wavefield B may be located at shot location A j from survey A.
- a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point r B of wavefield B and the shot location A j from survey A.
- the source and sensor definitions of wavefield A may be interchanged to obtain virtual wavefield A
- source and sensor definitions of wavefield B may be interchanged to obtain virtual wavefield B.
- the resultant virtual wavefields A and B are illustrated in FIG. 4D .
- virtual wavefield A originates at virtual source location 402 at location Si, travels upward to be reflected downward by the sea surface 11 , travels downward until reflected upward by the sea floor 25 , and travels upward to be detected by the first (intermediate) virtual sensor 404 at location A j .
- virtual wavefield B originates at virtual source location 402 at location Si, travels upward to be reflected downward by the sea surface 11 , travels downward until reflected upward by the sea floor 25 , travels upward to be reflected downward at reflection location r B , travels downward to be reflected upward at subsurface location P, and travels upward to be detected by the second (final) virtual sensor 406 at location B k .
- the subsurface location in the geological formation of interest may be imaged, either using a deconvolution or cross-correlation imaging process.
- the geophysical analysis equipment may be configured to deconvolve virtual wavefield B using virtual wavefield A to obtain deconvolved virtual wavefield B/A.
- the deconvolved virtual wavefield B/A may be effectively used as the imaging wavefield signal 410 .
- the imaging wavefield signal 410 may be considered as originating at a virtual source 412 at location A j (which is effectively the same as, or made to correspond with, r B ), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by the virtual sensor 414 at location B k .
- the geophysical analysis equipment may be configured to perform cross-correlation in order to image subsurface location P.
- a first data signal may be the data signal received at the first (intermediate) virtual sensor 404 at location A j
- a second data signal may be the data signal received at the second (final) virtual sensor 406 at location B k .
- the cross-correlation of the first data signal and the second data signal may provide the medium response (i.e., the Green's function) of the subsurface location P.
- the result is similar to imaging using a virtual source-receiver pair, where the virtual source 412 and virtual receiver 414 (shown in FIG. 4E ) are located at the positions of the virtual sensors 404 and 406 , respectively (shown in FIG. 4D ).
- the resultant image data and other pertinent data is stored in a tangible, non-volatile data storage medium. Tangible medium holding such image and other pertinent data may subsequently be provided as a geophysical data product.
- the series of FIGS. 4A-4E illustrates a first exemplary application of the method 300 where a subsurface region 29 may be difficult or impossible to image using the data from either survey alone.
- the series of FIGS. 5A-5E illustrates a second exemplary application of the method 300 where an obstruction at the sea surface blocks shot locations such that one survey has missing shots.
- FIG. 5A depicts an illustrative example of a first seismic survey (A) in accordance with an embodiment of the invention. Shown in FIG. 5A are: shot locations of survey A beneath the sea surface 11 , and seabed sensors 22 on the sea floor 25 .
- FIG. 5B depicts a second seismic survey (B) in the presence of a surface obstruction B in accordance with an embodiment of the invention.
- the surface obstruction B may be, for example, an offshore production platform, an offshore drilling platform, or a drillship. As depicted, the surface obstruction B blocks one or more shot locations that were planned to be taken during survey B.
- FIG. 5C depicts wavefields from the first and second seismic surveys (A and B) that are utilized in combination to image a subsurface location (P) that may be underneath the obstruction in accordance with an embodiment of the invention.
- FIG. 5C in the series of FIGS. 5A-5E corresponds to FIG. 4C in the series of FIGS. 4A-4E .
- the selected wavefield B originates at shot location B k from survey B and travels downward.
- Wavefield B reflects from subsurface point P in the subsurface region 29 so as to travel upward, then reflects from point r B of the sea surface 11 so as to travel downward.
- Wavefield B may then reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 again to travel downward until it is detected by a seabed sensor at the seabed position S i .
- the selected wavefield A originates at shot location A j and travels downward.
- Wavefield A may reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position S i as that which detected wavefield B.
- the first sea-surface reflection point r B of wavefield B may be located at shot location A j from survey A.
- a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point r B of wavefield B and the shot location A j from survey A.
- FIG. 5D depicts virtual wavefields after reversal of source and sensor definitions of FIG. 5C in accordance with an embodiment of the invention.
- FIG. 5D in the series of FIGS. 5A-5E corresponds to FIG. 4D in the series of FIGS. 4A-4E .
- virtual wavefield A originates at virtual source location 502 at location Si, travels upward to be reflected downward by the sea surface 11 , travels downward until reflected upward by the sea floor 25 , and travels upward to be detected by the first (intermediate) virtual sensor 504 at location A j .
- virtual wavefield B originates at virtual source location 502 at location Si, travels upward to be reflected downward by the sea surface 11 , travels downward until reflected upward by the sea floor 25 , travels upward to be reflected downward at reflection location r B , travels downward to be reflected upward at subsurface location P, and travels upward to be detected by the second (final) virtual sensor 506 at location B k .
- FIG. 5E depicts the virtual wavefield utilized for imaging the subsurface location (P) after deconvolution or cross-correlation is applied to the virtual wavefields of FIG. 5 D in accordance with an embodiment of the invention.
- FIG. 5E in the series of FIGS. 5A-5E corresponds to FIG. 4E in the series of FIGS. 4A-4E .
- the imaging wavefield signal 510 may be considered as originating at a virtual source 512 at location A j (which is effectively the same as, or made to correspond with, r B ), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by the virtual sensor 514 at location B k .
- the result is that subsurface location P in the subsurface region 29 may be imaged, despite the missing shot(s) due to the obstruction at the sea surface during survey B.
- FIGS. 6A-6E illustrates a third exemplary case where two surveys have different shot intervals from each other. Despite the different spacings between the shots, select data signals may be processed in combination as described below.
- FIG. 6A depicts a first seismic survey (A) in accordance with an embodiment of the invention. Shown in FIG. 6A are: shot locations of survey A beneath the sea surface 11 , and seabed sensors 22 on the sea floor 25 .
- FIG. 6B depicts a second seismic survey (B) with different shot intervals than the first seismic survey (A) in accordance with an embodiment of the invention.
- the shots in survey B are spaced farther apart than the shots in survey A.
- FIG. 6C depicts wavefields from the first and second seismic surveys (A and B) that are utilized in combination to image a subsurface location (P) in accordance with an embodiment of the invention.
- FIG. 6C in the series of FIGS. 6A-6E corresponds to FIG. 4C in the series of FIGS. 4A-4E .
- the selected wavefield B originates at shot location B k from survey B and travels downward.
- Wavefield B reflects from subsurface point P at a subsurface level 30 so as to travel upward, then reflects from point r B of the sea surface 11 so as to travel downward.
- Wavefield B may then reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 again to travel downward until it is detected by a seabed sensor at the seabed position S i .
- the selected wavefield A originates at shot location A j and travels downward.
- Wavefield A may reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position S i as that which detected wavefield B.
- the first sea-surface reflection point r B of wavefield B may be located at shot location A j from survey A.
- a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point r B of wavefield B and the shot location A j from survey A.
- FIG. 6D depicts virtual wavefields after reversal of source and sensor definitions of FIG. 6C in accordance with an embodiment of the invention.
- FIG. 6D in the series of FIGS. 6A-6E corresponds to FIG. 4D in the series of FIGS. 4A-4E .
- virtual wavefield A originates at virtual source location 602 at location Si, travels upward to be reflected downward by the sea surface 11 , travels downward until reflected upward by the sea floor 25 , and travels upward to be detected by the first (intermediate) virtual sensor 604 at location A j .
- virtual wavefield B originates at virtual source location 602 at location Si, travels upward to be reflected downward by the sea surface 11 , travels downward until reflected upward by the sea floor 25 , travels upward to be reflected downward at reflection location r B , travels downward to be reflected upward at subsurface location P, and travels upward to be detected by the second (final) virtual sensor 606 at location B k .
- FIG. 6E depicts the virtual wavefield utilized for imaging the subsurface location (P) after deconvolution or cross-correlation is applied to the virtual wavefields of FIG. 6D in accordance with an embodiment of the invention.
- FIG. 6E in the series of FIGS. 6A-6E corresponds to FIG. 4E in the series of FIGS. 4A-4E .
- the imaging wavefield signal 610 may be considered as originating at a virtual source 612 at location A j (which is effectively the same as, or made to correspond with, r B ), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by the virtual sensor 614 at location B k .
- the result is that subsurface location P at subsurface level 30 may be imaged using a combination of data from surveys A and B.
- FIG. 7 is a flow chart of a method 700 using seabed sensors and sea surface reflections for four-dimensional (4D) surveying for reservoir production monitoring in accordance with an embodiment of the invention.
- the steps 702 - 728 of FIG. 7 are discussed below in relation to the exemplary case illustrated in FIGS. 8A-8I .
- wavefields A 0 and A 1 from survey A are selected. As depicted in FIG. 8A , the two wavefields A 0 and A 1 originate at different shot locations and end at the same sensor. The selected wavefield A 1 originates at shot location A k from survey A and travels downward. Wavefield A 1 reflects from subsurface point P in the reservoir 31 so as to travel upward, then reflects from point r 1 of the sea surface 11 so as to travel downward. Wavefield A 1 may then reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 again to travel downward until it is detected by a seabed sensor at the seabed position S i .
- the selected wavefield A 0 originates at shot location A j and travels downward.
- Wavefield A 0 may reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position S i as that which detected wavefield A 1 .
- the first sea-surface reflection point r 1 of wavefield A 1 may be located at shot location A j from survey A 0 .
- a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point r 1 of wavefield A 1 and the shot location A j from survey A 0 .
- FIG. 8B depicts virtual wavefields after reversal of source and sensor definitions of FIG. 8A in accordance with an embodiment of the invention.
- virtual wavefield A 0 originates at virtual source location 802 at location Si, travels upward to be reflected downward by the sea surface 11 , travels downward until reflected upward by the sea floor 25 , and travels upward to be detected by the first (intermediate) virtual sensor 804 at location A j .
- virtual wavefield A 1 originates at virtual source location 802 at location Si, travels upward to be reflected downward by the sea surface 11 , travels downward until reflected upward by the sea floor 25 , travels upward to be reflected downward at reflection location r 1 , travels downward to be reflected upward at subsurface location P, and travels upward to be detected by the second (final) virtual sensor 806 -A at location A k .
- the subsurface location P in the reservoir 31 may be imaged, either using a deconvolution or cross-correlation imaging process.
- the resultant image data and other pertinent data may be stored in a tangible medium so as to be provided as a geophysical data product.
- the geophysical analysis equipment may be configured to deconvolve virtual wavefield A 1 using virtual wavefield A 0 to obtain deconvolved virtual wavefield A 1 /A 0 .
- the deconvolved virtual wavefield A 1 /A 0 may be effectively used as the imaging wavefield signal 810 .
- the imaging wavefield signal 810 may be considered as originating at a virtual source 812 at location A j (which is effectively the same as, or made to correspond with, r 1 ), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by the virtual sensor 814 at location A k .
- the geophysical analysis equipment may be configured to perform cross-correlation in order to image subsurface location P.
- a first data signal may be the virtual data signal received at the first virtual sensor at location A j
- a second data signal may be the virtual data signal received at the second virtual sensor at location A k .
- the cross-correlation of the first data signal and the second data signal may provide the medium response (i.e., the Green's function) of the subsurface location P.
- the result is similar to imaging using a virtual source-receiver pair, where the virtual source 812 and virtual receiver 814 (shown in FIG. 8C ) are located at the positions of the first and second virtual sensors, respectively.
- a second marine seismic survey (survey B) may be performed.
- the performance of survey B involves generating seismic signals at a second set of positions and recording a second set of received signals.
- Survey B may use an array of seabed sensors at the same positions as the array used in survey A.
- the seabed sensors in the array for survey B may be the same sensors as used for survey A, or the seabed sensors in the array for survey B may be redeployed into the same positions as the sensors used for survey A.
- survey B may use a sparser shot pattern than survey A.
- One exemplary shot location B k of the second survey (survey B) is shown in FIG. 8D .
- the shot location B k of survey B is at the same location, or effectively the same location, as the shot location A k of survey A.
- wavefields A 0 of survey A and B 1 from survey B are selected. As depicted in FIG. 8D , the two wavefields A 0 and B 1 originate at different shot locations and end at the same sensor position.
- the selected wavefield B 1 originates at shot location B k from survey B and travels downward.
- Wavefield B 1 reflects from subsurface point P in the reservoir 31 so as to travel upward, then reflects from point r 1 of the sea surface 11 so as to travel downward.
- Wavefield B 1 may then reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 again to travel downward until it is detected by a seabed sensor at the seabed position S i .
- the selected wavefield A 0 originates at shot location A j and travels downward.
- Wavefield A 0 may reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position S i as that which detected wavefield B 1 .
- the first sea-surface reflection point r 1 of wavefield B 1 may be located at shot location A j from survey A 0 .
- a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point r 1 of wavefield B 1 and the shot location A j from survey A 0 .
- FIG. 8E depicts virtual wavefields after reversal of source and sensor definitions of FIG. 8D in accordance with an embodiment of the invention.
- virtual wavefield A 0 originates at virtual source location 802 at location Si, travels upward to be reflected downward by the sea surface 11 , travels downward until reflected upward by the sea floor 25 , and travels upward to be detected by the first (intermediate) virtual sensor 804 at location A j .
- the subsurface location P in the reservoir 31 may be imaged, either using a deconvolution or cross-correlation imaging process. Comparison may then be made of the current image data of the subsurface location P generated in step 718 against the previous image data of the subsurface location P generated in step 708 . The comparison may involve generating difference data by taking a difference between the current and previous image data. The resultant image data, difference data, and other pertinent data may be stored in a tangible medium so as to be provided as a geophysical data product.
- the geophysical analysis equipment (on vessel or on shore), which uses a computer system to process the seismic data, may be configured to deconvolve virtual wavefield B 1 using virtual wavefield A 0 to obtain deconvolved virtual wavefield B 1 /A 0 .
- the deconvolved virtual wavefield B 1 /A 0 may be effectively used as the imaging wavefield signal 820 .
- the imaging wavefield signal 820 may be considered as originating at a virtual source 822 at location A j (which is effectively the same as, or made to correspond with, r 1 ), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by the virtual sensor 824 at location B k .
- the geophysical analysis equipment on vessel or on shore, which uses a computer system to process the seismic data, may be configured to perform cross-correlation in order to image subsurface location P.
- a first data signal may be the virtual data signal received at the first virtual sensor at location A j
- a second data signal may be the virtual data signal received at the second virtual sensor at location B k .
- the cross-correlation of the first data signal and the second data signal may provide the medium response (i.e., the Green's function) of the subsurface location P.
- the result is similar to imaging using a virtual source-receiver pair, where the virtual source 822 and virtual receiver 824 (shown in FIG. 8F ) are located at the positions of the first and second virtual sensors, respectively.
- a third marine seismic survey (survey C) may be performed.
- the performance of survey C involves generating seismic signals at a third set of positions and recording a third set of received signals.
- Survey C may use an array of seabed sensors at the same positions as the array used in survey A.
- the seabed sensors in the array for survey C may be the same sensors as used for survey A, or the seabed sensors in the array for survey C may be redeployed into the same positions as the sensors used for survey A.
- survey C may use a sparser shot pattern than survey A.
- One exemplary shot location C k of the third survey (survey C) is shown in FIG. 8G .
- the shot location C k of survey C is at the same location, or effectively the same location, as the shot location A k of survey A.
- wavefields A 0 of survey A and C 1 from survey C are selected. As depicted in FIG. 8G , the two wavefields A 0 and C 1 originate at different shot locations and end at the same sensor position.
- the selected wavefield C 1 originates at shot location C k from survey C and travels downward.
- Wavefield C 1 reflects from subsurface point P in the reservoir 31 so as to travel upward, then reflects from point r 1 of the sea surface 11 so as to travel downward.
- Wavefield C 1 may then reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 again to travel downward until it is detected by a seabed sensor at seabed position S i .
- the selected wavefield A 0 originates at shot location A j and travels downward.
- Wavefield A 0 may reflect from the sea floor 25 so as to travel upward, then reflect from the sea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position S i as that which detected wavefield C 1 .
- the first sea-surface reflection point r 1 of wavefield C 1 may be located at shot location A j from survey A 0 .
- a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point r 1 of wavefield C 1 and the shot location A j from survey A 0 .
- FIG. 8H depicts virtual wavefields after reversal of source and sensor definitions of FIG. 8G in accordance with an embodiment of the invention.
- virtual wavefield A 0 originates at virtual source location 802 at location Si, travels upward to be reflected downward by the sea surface 11 , travels downward until reflected upward by the sea floor 25 , and travels upward to be detected by the first (intermediate) virtual sensor 804 at location A j .
- the geophysical analysis equipment (on vessel or on shore), which uses a computer system to process the seismic data, may be configured to deconvolve virtual wavefield C 1 using virtual wavefield A 0 to obtain deconvolved virtual wavefield C 1 /A 0 .
- the deconvolved virtual wavefield C 1 /A 0 may be effectively used as the imaging wavefield signal 830 .
- the imaging wavefield signal 830 may be considered as originating at a virtual source 832 at location A j (which is effectively the same as, or made to correspond with, r 1 ), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by the virtual sensor 834 at location C k .
- the geophysical analysis equipment on vessel or on shore, which uses a computer system to process the seismic data, may be configured to perform cross-correlation in order to image subsurface location P.
- a first data signal may be the virtual data signal received at the first virtual sensor at location A j
- a second data signal may be the virtual data signal received at the second virtual sensor at location C k .
- the cross-correlation of the first data signal and the second data signal may provide the medium response (i.e., the Green's function) of the subsurface location P.
- the result is similar to imaging using a virtual source-receiver pair, where the virtual source 832 and virtual receiver 834 (shown in FIG. 8I ) are located at the positions of the first and second virtual sensors, respectively.
- the method 700 of FIG. 7 may continue with the performance and processing of further surveys (D, E, F, etc.). Each survey and the subsequent processing steps results in a “snapshot” image of the reservoir being monitored.
- the time gap between survey A and survey B may be typically of sufficient length for a change to be expected in the state of the reservoir 31 .
- the time gap between surveys A and B may be three months, six months, or a year, or longer.
- the time gap between survey B and survey C may be typically of sufficient length for a change to be expected in the state of the reservoir 31 .
- the time gap between surveys B and C may be three months, six months, or a year, or longer. And so on for further surveys to monitor a hydrocarbon reservoir over a span of years.
- the above-described method 700 of FIG. 7 may be used to efficiently monitor a production reservoir (i.e. perform a 4D survey).
- the above-described method 700 of FIG. 7 provides greater flexibility for subsequent shooting monitor surveys.
- the first survey may have a dense shot pattern for detailed 3D imaging
- the subsequent surveys may have a sparser shot pattern.
- the first survey may be recorded using a dense single shot survey
- the subsequent surveys may be recorded using cost-effective dual-source flip-flop acquisition.
- a time-lapse imaging sequence may be defined as: Base (A 1 /A 0 ); Monitor 1 (B 1 /A 0 ); Monitor 2 (C 1 /A 0 ), etc.
- FIG. 9 shows a simplified example of a computer apparatus 40 which may be used in performing processing steps in accordance with an embodiment of the invention.
- the computer apparatus 40 may be configured with executable instructions so as to perform the data processing methods described herein.
- This figure shows just one simplified example of a computer which may be used to perform the data processing methods described herein.
- Other types of computers may also be employed.
- a high-performance cluster of multi-processor server computers may be used in practice.
- the processing methods may also be performed using a cloud computing system.
- the simplified example of the computer apparatus 40 may include a processor 41 , such as those from the Intel Corporation of Santa Clara, Calif., for example.
- the computer apparatus. 40 may have a bus system 43 communicatively interconnecting its various components.
- the computer apparatus 40 may include one or more user input devices 42 (e.g., keyboard, mouse), a display monitor 44 (e.g., LCD, flat panel monitor, CRT), a computer network interface 45 (e.g., network adapter, modem), and a data storage system which may include one or more data storage devices 46 (e.g., hard drive, solid state memory, optical disk drive, USB memory) and a main memory 50 (e.g., RAM).
- data storage devices 46 e.g., hard drive, solid state memory, optical disk drive, USB memory
- main memory 50 e.g., RAM
- the main memory 50 includes executable code 52 and data 54 stored therein.
- the executable code 52 may comprise computer-readable program code (i.e., software) components which may be loaded from the data storage device 46 to the main memory 50 for execution by the processor 41 .
- the executable code 52 may be configured to perform computer-implemented steps in the methods described herein.
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Fluid Mechanics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Oceanography (AREA)
Abstract
Description
- This application claims the benefit of U.S. Provisional Application No. 62/168,988, filed on Jun. 1, 2015, which is incorporated by reference herein in its entirety.
- Geophysical surveys are often used for oil and gas exploration in geological formations, which may be located below marine environments. Various types of signal sources and geophysical sensors may be used in different types of geophysical surveys.
- Seismic geophysical surveys, for example, are based on the use of acoustic waves. In a marine seismic survey, a survey vessel may tow an acoustic source (e.g., an air gun or a marine vibrator) and a plurality of streamers along which a number of acoustic sensors (e.g., hydrophones and/or geophones) are located. Alternatively, acoustic sensors may be located on the seabed rather than being towed.
- Acoustic waves generated by the survey source may be transmitted to the Earth's crust and then reflected back and captured at the towed and/or seabed geophysical sensors. Data collected during a marine geophysical survey may be analyzed to locate hydrocarbon-bearing geological structures, and thus determine where deposits of oil and natural gas may be located.
-
FIG. 1 is a diagram illustrating a geophysical survey system in accordance with an embodiment of the invention. -
FIG. 2 depicts an exemplary array of seabed sensors in accordance with an embodiment of the invention. -
FIG. 3 is a flow chart of a method using seabed sensors and sea surface reflections for structural imaging of a subsurface location in a geological formation in accordance with an embodiment of the invention. -
FIG. 4A depicts shot locations of a first seismic survey (A) in accordance with an embodiment of the invention. -
FIG. 4B depicts shot locations of a second seismic survey (B) in accordance with an embodiment of the invention. -
FIG. 4C depicts wavefields from the first and second seismic surveys (A and B) that are utilized in combination to image a subsurface location (P) in accordance with an embodiment of the invention. -
FIG. 4D depicts virtual wavefields after reversal of source and sensor definitions ofFIG. 4C in accordance with an embodiment of the invention. -
FIG. 4E depicts the virtual wavefield utilized for imaging the subsurface location (P) after deconvolution or cross-correlation is applied to the virtual wavefields ofFIG. 4D in accordance with an embodiment of the invention. -
FIG. 5A depicts shot locations of a first seismic survey (A) in accordance with an embodiment of the invention. -
FIG. 5B depicts shot locations of a second seismic survey (B) in the presence of a surface obstruction in accordance with an embodiment of the invention. -
FIG. 5C depicts wavefields from the first and second seismic surveys (A and B) that are utilized in combination to image a subsurface location (P) that may be underneath the obstruction in accordance with an embodiment of the invention. -
FIG. 5D depicts virtual wavefields after reversal of source and sensor definitions ofFIG. 5C in accordance with an embodiment of the invention. -
FIG. 5E depicts the virtual wavefield utilized for imaging the subsurface location (P) after deconvolution or cross-correlation is applied to the virtual wavefields ofFIG. 5D in accordance with an embodiment of the invention. -
FIG. 6A depicts shot locations of a first seismic survey (A) in accordance with an embodiment of the invention. -
FIG. 6B depicts shot locations of a second seismic survey (B) with different shot intervals than the first seismic survey (A) in accordance with an embodiment of the invention. -
FIG. 6C depicts wavefields from the first and second seismic surveys (A and B) that are utilized in combination to image a subsurface location (P) in accordance with an embodiment of the invention. -
FIG. 6D depicts virtual wavefields after reversal of source and sensor definitions ofFIG. 6C in accordance with an embodiment of the invention. -
FIG. 6E depicts the virtual wavefield utilized for imaging the subsurface location (P) after deconvolution or cross-correlation is applied to the virtual wavefields ofFIG. 6D in accordance with an embodiment of the invention. -
FIG. 7 is a flow chart of a method using seabed sensors and sea surface reflections for reservoir production monitoring in accordance with an embodiment of the invention. -
FIG. 8A depicts wavefields (A0 and A1) from a first seismic survey (A) at a first time that may be utilized in combination to image a subsurface location (P) at the first time in accordance with an embodiment of the invention. -
FIG. 8B depicts virtual wavefields after reversal of source and sensor definitions ofFIG. 8A in accordance with an embodiment of the invention. -
FIG. 8C depicts the virtual wavefield utilized for imaging the subsurface location (P) at the first time in accordance with an embodiment of the invention. -
FIG. 8D depicts a wavefield (B1) from a second seismic survey (B) at a second time that is utilized in combination with a wavefield (A0) from the first seismic survey (A) to image the subsurface location (P) at the second time in accordance with an embodiment of the invention. -
FIG. 8E depicts virtual wavefields after reversal of source and sensor definitions ofFIG. 8D in accordance with an embodiment of the invention. -
FIG. 8F depicts the virtual wavefield utilized for imaging the subsurface location (P) at the second time in accordance with an embodiment of the invention. -
FIG. 8G depicts a wavefield (C1) from a third seismic survey (C) at a third time that is utilized in combination with a wavefield (A0) from the first seismic survey (A) to image the subsurface location (P) at the third time in accordance with an embodiment of the invention. -
FIG. 8H depicts virtual wavefields after reversal of source and sensor definitions ofFIG. 8G in accordance with an embodiment of the invention. -
FIG. 8I depicts the virtual wavefield utilized for imaging the subsurface location (P) at the third time in accordance with an embodiment of the invention. -
FIG. 9 shows a simplified example of a computer system which may be used in performing processing steps in accordance with an embodiment of the invention. - Three-dimensional (3D) seismic surveying generally uses one or more acoustic (seismic) sources and an array of acoustic (seismic) sensors. A pattern of shots is performed by the source(s). For each shot, each sensor in the array senses and records a seismic wavefield. The seismic data obtained includes multiple seismic wavefields for each of multiple shots. The seismic data may be processed to create images of a subsurface volume underneath the array of sensors.
- A 3D seismic survey design includes the positions of the shots and the locations of the sensors. In the case of marine seismic surveying using seabed sensors, the locations of the sensors are fixed. In that case, the positions of the shots (the shot survey) may be designed to obtain seismic data covering the desired subsurface volume.
- Four-dimensional (4D) seismic surveying involves surveying the same geological formation at different points in time. For example, 4D surveying may be used to determine changes over time in fluid saturation, pressure, and so on. Because seabed sensors are fixed-position receivers that do not change position during the seismic acquisition, they may be redeployed at the same positions at a later time so as to be advantageously used in 4D marine seismic surveying to reduce variations between surveys.
- One problem or difficulty with marine seismic surveys may be caused by an obstruction that prevents positioning shots and/or sensors above a targeted sub-surface region. Such an obstruction makes it difficult to obtain seismic images of the sub-surface region under the obstruction in a conventional 3D or 4D seismic survey.
- For example, a production platform may be located above the targeted sub-surface region. The production platform forms an obstruction that blocks shots from being taken at sea-surface positions in the area above the targeted sub-surface region. In the case of towed sensors, the production platform may also block sensors from being towed above the targeted region.
- The present disclosure presents an innovative technique to generate seismic images by combining seismic data from different surveys. The technique has various practical applications, such as imaging under an obstruction and efficient reservoir monitoring.
- Referring to
FIG. 1 , a diagram illustrating one embodiment of aGeophysical survey system 100 is shown. In the illustrated embodiment, thesystem 100 includes asurvey vessel 10, signalsource 32, andgeophysical sensors 22. - The
survey vessel 10 may be configured to move along thesurface 11 of body ofwater 15 such as a lake or the ocean, as shown. In the illustrated embodiment, signalsource 32 is configured to transmit an acoustic signal. The illustrated dashed line shows a direction of travel for a portion of an acoustic wave-front generated bysignal source 32. In the illustrated embodiment, the acoustic signal proceeds through thesea floor 25 and is reflected byreservoir 120, the sea surface (an air/water interface) 11, and areflector 27 located abovereservoir 120. In other embodiments, the acoustic signal may be reflected from the sea surface multiple times (and reflected from the seabed in between the sea surface reflections). -
Survey vessel 10 may include equipment, shown generally at 12 and for convenience collectively referred to as “survey equipment.”Survey equipment 12 may include devices such as a data recording unit for making a record of signals generated by various geophysical sensors in thesystem 100.Survey equipment 12 may also include navigation equipment, which may be configured to control, determine, and record, at selected times, the geodetic positions of:survey vessel 10, signalsource 32, and/or other survey elements. Geodetic position may be determined using various devices, including global navigation satellite systems such as the global positioning system (GPS), for example. In some embodiments,vessel 10 is configured to tow one or more seismic streamers (not shown). - In the
geophysical survey system 100 shown inFIG. 1 ,survey vessel 10 tows asingle signal source 32. In various embodiments,survey vessel 10 may tow any appropriate number of signal sources, including as few as none or as many as six or more. The location of the signal sources may be centered behindsurvey vessel 10 or displaced from the center line, and may be at various distances relative to surveyvessel 10, including attached to the hull. Signal sources may be any type of signal source known in the art. Signalsource 32 may include an array of multiple signal sources. For example, signalsource 32 may include a plurality of seismic sources (e.g., air guns). The term “signal source” may refer to a single signal source or to an array of signal sources. For example, in the illustrated embodiment,signal sources 32 are each coupled to surveyvessel 10 at one end throughwinch 115 or a similar spooling device that enables changing the deployed length of a signal source cable.Survey equipment 12 may include signal source control equipment (not shown separately) for selectively operating andmaneuvering signal source 32. The specific configuration depicted inFIG. 1 is an example only. Other configurations may be used for the geophysical survey system. -
Geophysical sensors 22 may be referred to as seabed sensors and may be any type of geophysical sensor known in the art. In an exemplary implementation, the geophysical sensors may include pressure sensors (e.g., hydrophones) that measure pressure (P) as a scalar and particle-velocity sensors (e.g. geophones and accelerometers) that measure particle velocity (Z) as a vector. The pressure and particle-velocity sensors may be both included in a given seabed sensor. Other examples of geophysical sensors that may be used include pressure-time-gradient responsive seismic sensors, electrodes, magnetometers, temperature sensors, fiber optic sensors, or combinations of the foregoing. In various implementations of the disclosure,geophysical sensors 22 may measure, for example, seismic energy indicative of the response of various structures in the Earth's subsurface formation below the bottom 25 of body ofwater 15 to energy imparted into the subsurface formation by one or more ofsignal source 32. Seismic energy, for example, may originate fromsignal source 32, or an array of such signal sources, deployed in body ofwater 15 and towed bysurvey vessel 10. - In some embodiments, geophysical sensors 22 (which may also be referred to as receivers) are configured to store data until they are physically retrieved, at which point the stored data may be accessed. In other embodiments,
geophysical sensors 22 may be configured to communicate with survey equipment onshore or in survey vessel wirelessly or using transmission cables, in various embodiments.Geophysical sensors 22 may include relatively stable and/or synchronized clocks to timestamp recorded data.Geophysical sensors 22 may be comprised in ocean bottom cables (OBC) and linked together or may be ocean bottom nodes (OBN) and be individual recording systems (e.g., with a clock, battery, sensor, and data/time recorder integrated within each unit). In the illustrated embodiment, two ofgeophysical sensors 22 are located distance dl apart. -
Survey equipment 12, in one embodiment, includes a computing system configured to, inter alia, process data fromgeophysical sensors 22. In other embodiments, a computing system at another location may perform geophysical analysis on data gathered by geophysical survey system 100 (e.g., on land after a survey has been conducted). The computing system may include or be configured to access a non-transitory storage medium having instructions stored thereon that are executable to cause the system to perform various operations described herein. The computing system may include one or more processors configured to execute the program instructions and one or more data storage devices for storing results in a tangible medium. - The pattern of an array of seabed sensors for a survey may be determined in advance in order to generate data usable to perform the disclosed geophysical analysis techniques. The pattern may specify distances between sensors, sensor orientation, shot area, and other parameters.
-
FIG. 2 depicts an exemplary array of seabed sensors in accordance with an embodiment of the invention. As depicted, this exemplary array includes 380 sensors in a “raster-shaped” layout pattern. For purposes of illustration, the numbers depicted indicate the seabed positions ofsensors # 1, #50, #100, #165, #235, #305 and #380. This layout pattern may be particularly useful in separated wavefield imaging by multiples (SWIM). - Exemplary Method Combining Data from Multiple Surveys for 3-Dimensional Structural Imaging
-
FIG. 3 is a flow chart of amethod 300 using seabed sensors and sea surface reflections for structural imaging of a subsurface location in a geological formation in accordance with an embodiment of the invention. Themethod 300 combines data from multiple surveys. The steps ofFIG. 3 are discussed below in relation to the exemplary case illustrated inFIGS. 4A-4E . - As described below,
FIGS. 4A-4E illustrate the method steps ofFIG. 3 used to acquire seismic data in multiple surveys (A and B) and process the acquired data in combination to image a subsurface location P in asubsurface region 29. Thesubsurface region 29 may be difficult or impossible to image using the data from either survey alone due to the geometries of the shot locations relative to the sensor positions. Nevertheless, thesubsurface region 29 may be imaged in three-dimensions (3D imaging) using the method described below for multiple subsurface points in thesubsurface region 29. - Per step 302-A, a first marine seismic survey (survey A) may be performed. The performance of survey A involves generating seismic signals at a first set of positions and recording a first set of received signals using an array of seabed sensors.
- An illustrative example of the first survey (survey A) is shown in
FIG. 4A . Shown inFIG. 4A are: shot locations of survey A beneath thesea surface 11, andseabed sensors 22 on thesea floor 25. - Per step 302-B, a second marine seismic survey (survey B) may be performed. The performance of survey B involves generating seismic signals at a second set of positions and recording a second set of received signals using the same array (location) of seabed sensors.
- An illustrative example of the second survey (survey B) is shown in
FIG. 4B . Shown inFIG. 4B are: shot locations of survey B beneath thesea surface 11, and the array ofseabed sensors 22 on thesea floor 25. Note that, in accordance with the illustrated embodiment, the locations of theseabed sensors location 22 do not change between surveys A and B. - Per
step 303, up-going and down-going wavefields are separated for each survey. In an exemplary implementation, this may be accomplished by computing the summation and difference between the pressure (P) and particle velocity (Z) fields. Because the up-going wavefield is in-phase for both pressure and particle-velocity records, the up-going wavefield may be extracted by summing the two measurements (PZ-summation) after calibration. The down-going wavefield is canceled in the PZ-summation because it is out-of-phase on the P field compared to the Z field. The down-going wavefield may be extracted by determining the difference between the P field and the calibrated Z field. - Per
step 304, wavefield A from survey A and wavefield B from survey B are selected. Each selected wavefield originates at a shot location and ends at a sensor. - As illustrated in
FIG. 4C , the selected wavefield B originates at shot location Bk from survey B and travels downward. Wavefield B reflects from subsurface point P in thesubsurface region 29 so as to travel upward, then reflects from point rB of thesea surface 11 so as to travel downward. Wavefield B may then reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 again to travel downward until it is detected by a seabed sensor at seabed position Si. - As further illustrated in
FIG. 4C , the selected wavefield A originates at shot location Aj and travels downward. Wavefield A may reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position Si as that which detected wavefield B. Note that the first sea-surface reflection point rB of wavefield B may be located at shot location Aj from survey A. Alternatively, a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point rB of wavefield B and the shot location Aj from survey A. - Per
step 306, the source and sensor definitions of wavefield A may be interchanged to obtain virtual wavefield A, and source and sensor definitions of wavefield B may be interchanged to obtain virtual wavefield B. The resultant virtual wavefields A and B are illustrated inFIG. 4D . - As depicted in
FIG. 4D , virtual wavefield A originates atvirtual source location 402 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, and travels upward to be detected by the first (intermediate) virtual sensor 404 at location Aj. As further depicted, virtual wavefield B originates atvirtual source location 402 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, travels upward to be reflected downward at reflection location rB, travels downward to be reflected upward at subsurface location P, and travels upward to be detected by the second (final) virtual sensor 406 at location Bk. - Per
step 308, the subsurface location in the geological formation of interest may be imaged, either using a deconvolution or cross-correlation imaging process. In one embodiment, the geophysical analysis equipment may be configured to deconvolve virtual wavefield B using virtual wavefield A to obtain deconvolved virtual wavefield B/A. The deconvolved virtual wavefield B/A may be effectively used as theimaging wavefield signal 410. - As depicted in
FIG. 4E , the imaging wavefield signal 410 may be considered as originating at avirtual source 412 at location Aj (which is effectively the same as, or made to correspond with, rB), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by thevirtual sensor 414 at location Bk. - In an alternate embodiment, the geophysical analysis equipment may be configured to perform cross-correlation in order to image subsurface location P. In this embodiment, a first data signal may be the data signal received at the first (intermediate) virtual sensor 404 at location Aj, and a second data signal may be the data signal received at the second (final) virtual sensor 406 at location Bk. The cross-correlation of the first data signal and the second data signal may provide the medium response (i.e., the Green's function) of the subsurface location P. The result is similar to imaging using a virtual source-receiver pair, where the
virtual source 412 and virtual receiver 414 (shown inFIG. 4E ) are located at the positions of the virtual sensors 404 and 406, respectively (shown inFIG. 4D ). - Per
step 310, the resultant image data and other pertinent data is stored in a tangible, non-volatile data storage medium. Tangible medium holding such image and other pertinent data may subsequently be provided as a geophysical data product. - As described above, the series of
FIGS. 4A-4E illustrates a first exemplary application of themethod 300 where asubsurface region 29 may be difficult or impossible to image using the data from either survey alone. Similarly, the series ofFIGS. 5A-5E illustrates a second exemplary application of themethod 300 where an obstruction at the sea surface blocks shot locations such that one survey has missing shots. - Like
FIG. 4A ,FIG. 5A depicts an illustrative example of a first seismic survey (A) in accordance with an embodiment of the invention. Shown inFIG. 5A are: shot locations of survey A beneath thesea surface 11, andseabed sensors 22 on thesea floor 25. -
FIG. 5B depicts a second seismic survey (B) in the presence of a surface obstruction B in accordance with an embodiment of the invention. The surface obstruction B may be, for example, an offshore production platform, an offshore drilling platform, or a drillship. As depicted, the surface obstruction B blocks one or more shot locations that were planned to be taken during survey B. -
FIG. 5C depicts wavefields from the first and second seismic surveys (A and B) that are utilized in combination to image a subsurface location (P) that may be underneath the obstruction in accordance with an embodiment of the invention.FIG. 5C in the series ofFIGS. 5A-5E corresponds toFIG. 4C in the series ofFIGS. 4A-4E . - As illustrated in
FIG. 5C , the selected wavefield B originates at shot location Bk from survey B and travels downward. Wavefield B reflects from subsurface point P in thesubsurface region 29 so as to travel upward, then reflects from point rB of thesea surface 11 so as to travel downward. Wavefield B may then reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 again to travel downward until it is detected by a seabed sensor at the seabed position Si. - As further illustrated in
FIG. 5C , the selected wavefield A originates at shot location Aj and travels downward. Wavefield A may reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position Si as that which detected wavefield B. Note that the first sea-surface reflection point rB of wavefield B may be located at shot location Aj from survey A. Alternatively, a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point rB of wavefield B and the shot location Aj from survey A. -
FIG. 5D depicts virtual wavefields after reversal of source and sensor definitions ofFIG. 5C in accordance with an embodiment of the invention.FIG. 5D in the series ofFIGS. 5A-5E corresponds toFIG. 4D in the series ofFIGS. 4A-4E . - As depicted in
FIG. 5D , virtual wavefield A originates atvirtual source location 502 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, and travels upward to be detected by the first (intermediate) virtual sensor 504 at location Aj. As further depicted, virtual wavefield B originates atvirtual source location 502 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, travels upward to be reflected downward at reflection location rB, travels downward to be reflected upward at subsurface location P, and travels upward to be detected by the second (final) virtual sensor 506 at location Bk. -
FIG. 5E depicts the virtual wavefield utilized for imaging the subsurface location (P) after deconvolution or cross-correlation is applied to the virtual wavefields of FIG. 5D in accordance with an embodiment of the invention.FIG. 5E in the series ofFIGS. 5A-5E corresponds toFIG. 4E in the series ofFIGS. 4A-4E . - As depicted in
FIG. 5E , the imaging wavefield signal 510 may be considered as originating at avirtual source 512 at location Aj (which is effectively the same as, or made to correspond with, rB), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by thevirtual sensor 514 at location Bk. The result is that subsurface location P in thesubsurface region 29 may be imaged, despite the missing shot(s) due to the obstruction at the sea surface during survey B. - Combining Data from Surveys with Different Shot Intervals
- Also similar to the series of
FIGS. 4A-4E , the series ofFIGS. 6A-6E illustrates a third exemplary case where two surveys have different shot intervals from each other. Despite the different spacings between the shots, select data signals may be processed in combination as described below. - Like
FIG. 4A ,FIG. 6A depicts a first seismic survey (A) in accordance with an embodiment of the invention. Shown inFIG. 6A are: shot locations of survey A beneath thesea surface 11, andseabed sensors 22 on thesea floor 25. -
FIG. 6B depicts a second seismic survey (B) with different shot intervals than the first seismic survey (A) in accordance with an embodiment of the invention. In the illustrated example, the shots in survey B are spaced farther apart than the shots in survey A. -
FIG. 6C depicts wavefields from the first and second seismic surveys (A and B) that are utilized in combination to image a subsurface location (P) in accordance with an embodiment of the invention.FIG. 6C in the series ofFIGS. 6A-6E corresponds toFIG. 4C in the series ofFIGS. 4A-4E . - As illustrated in
FIG. 6C , the selected wavefield B originates at shot location Bk from survey B and travels downward. Wavefield B reflects from subsurface point P at asubsurface level 30 so as to travel upward, then reflects from point rB of thesea surface 11 so as to travel downward. Wavefield B may then reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 again to travel downward until it is detected by a seabed sensor at the seabed position Si. - As further illustrated in
FIG. 6C , the selected wavefield A originates at shot location Aj and travels downward. Wavefield A may reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position Si as that which detected wavefield B. Note that the first sea-surface reflection point rB of wavefield B may be located at shot location Aj from survey A. Alternatively, a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point rB of wavefield B and the shot location Aj from survey A. -
FIG. 6D depicts virtual wavefields after reversal of source and sensor definitions ofFIG. 6C in accordance with an embodiment of the invention.FIG. 6D in the series ofFIGS. 6A-6E corresponds toFIG. 4D in the series ofFIGS. 4A-4E . - As depicted in
FIG. 6D , virtual wavefield A originates atvirtual source location 602 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, and travels upward to be detected by the first (intermediate) virtual sensor 604 at location Aj. As further depicted, virtual wavefield B originates atvirtual source location 602 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, travels upward to be reflected downward at reflection location rB, travels downward to be reflected upward at subsurface location P, and travels upward to be detected by the second (final) virtual sensor 606 at location Bk. -
FIG. 6E depicts the virtual wavefield utilized for imaging the subsurface location (P) after deconvolution or cross-correlation is applied to the virtual wavefields ofFIG. 6D in accordance with an embodiment of the invention.FIG. 6E in the series ofFIGS. 6A-6E corresponds toFIG. 4E in the series ofFIGS. 4A-4E . - As depicted in
FIG. 6E , the imaging wavefield signal 610 may be considered as originating at avirtual source 612 at location Aj (which is effectively the same as, or made to correspond with, rB), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by thevirtual sensor 614 at location Bk. The result is that subsurface location P atsubsurface level 30 may be imaged using a combination of data from surveys A and B. - Exemplary Method Combining Data from Multiple Surveys for 4-Dimensional Surveying
-
FIG. 7 is a flow chart of amethod 700 using seabed sensors and sea surface reflections for four-dimensional (4D) surveying for reservoir production monitoring in accordance with an embodiment of the invention. The steps 702-728 ofFIG. 7 are discussed below in relation to the exemplary case illustrated inFIGS. 8A-8I . - Per step 702, a first marine seismic survey (survey A) may be performed. The performance of survey A involves generating seismic signals at a first set of positions and recording a first set of received signals using an array of seabed sensors. An illustrative example of the first survey (survey A) is shown in
FIG. 8A . Shown inFIG. 8A are: shot locations of survey A beneath thesea surface 11, andseabed sensors 22 on thesea floor 25. - Per step 704, after separating up-going and down-going wavefields, wavefields A0 and A1 from survey A are selected. As depicted in
FIG. 8A , the two wavefields A0 and A1 originate at different shot locations and end at the same sensor. The selected wavefield A1 originates at shot location Ak from survey A and travels downward. Wavefield A1 reflects from subsurface point P in thereservoir 31 so as to travel upward, then reflects from point r1 of thesea surface 11 so as to travel downward. Wavefield A1 may then reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 again to travel downward until it is detected by a seabed sensor at the seabed position Si. The selected wavefield A0 originates at shot location Aj and travels downward. Wavefield A0 may reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position Si as that which detected wavefield A1. Note that the first sea-surface reflection point r1 of wavefield A1 may be located at shot location Aj from survey A0. Alternatively, a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point r1 of wavefield A1 and the shot location Aj from survey A0. - Per
step 706, the source and sensor definitions of wavefield A0 may be interchanged to obtain virtual wavefield A0, and source and sensor definitions of wavefield A1 may be interchanged to obtain virtual wavefield A1.FIG. 8B depicts virtual wavefields after reversal of source and sensor definitions ofFIG. 8A in accordance with an embodiment of the invention. As depicted inFIG. 8B , virtual wavefield A0 originates atvirtual source location 802 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, and travels upward to be detected by the first (intermediate) virtual sensor 804 at location Aj. As further depicted, virtual wavefield A1 originates atvirtual source location 802 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, travels upward to be reflected downward at reflection location r1, travels downward to be reflected upward at subsurface location P, and travels upward to be detected by the second (final) virtual sensor 806-A at location Ak. - Per step 708, the subsurface location P in the
reservoir 31 may be imaged, either using a deconvolution or cross-correlation imaging process. The resultant image data and other pertinent data may be stored in a tangible medium so as to be provided as a geophysical data product. - In one embodiment, the geophysical analysis equipment may be configured to deconvolve virtual wavefield A1 using virtual wavefield A0 to obtain deconvolved virtual wavefield A1/A0. The deconvolved virtual wavefield A1/A0 may be effectively used as the
imaging wavefield signal 810. As depicted inFIG. 8C , the imaging wavefield signal 810 may be considered as originating at avirtual source 812 at location Aj (which is effectively the same as, or made to correspond with, r1), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by thevirtual sensor 814 at location Ak. - In an alternate embodiment, the geophysical analysis equipment (on vessel or on shore) may be configured to perform cross-correlation in order to image subsurface location P. In this embodiment, a first data signal may be the virtual data signal received at the first virtual sensor at location Aj, and a second data signal may be the virtual data signal received at the second virtual sensor at location Ak. The cross-correlation of the first data signal and the second data signal may provide the medium response (i.e., the Green's function) of the subsurface location P. The result is similar to imaging using a virtual source-receiver pair, where the
virtual source 812 and virtual receiver 814 (shown inFIG. 8C ) are located at the positions of the first and second virtual sensors, respectively. - Per
step 712, after a first time gap, a second marine seismic survey (survey B) may be performed. The performance of survey B involves generating seismic signals at a second set of positions and recording a second set of received signals. Survey B may use an array of seabed sensors at the same positions as the array used in survey A. The seabed sensors in the array for survey B may be the same sensors as used for survey A, or the seabed sensors in the array for survey B may be redeployed into the same positions as the sensors used for survey A. In an exemplary implementation, survey B may use a sparser shot pattern than survey A. One exemplary shot location Bk of the second survey (survey B) is shown inFIG. 8D . The shot location Bk of survey B is at the same location, or effectively the same location, as the shot location Ak of survey A. - Per step 714, after separating up-going and down-going wavefields, wavefields A0 of survey A and B1 from survey B are selected. As depicted in
FIG. 8D , the two wavefields A0 and B1 originate at different shot locations and end at the same sensor position. The selected wavefield B1 originates at shot location Bk from survey B and travels downward. Wavefield B1 reflects from subsurface point P in thereservoir 31 so as to travel upward, then reflects from point r1 of thesea surface 11 so as to travel downward. Wavefield B1 may then reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 again to travel downward until it is detected by a seabed sensor at the seabed position Si. The selected wavefield A0 originates at shot location Aj and travels downward. Wavefield A0 may reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position Si as that which detected wavefield B1. Note that the first sea-surface reflection point r1 of wavefield B1 may be located at shot location Aj from survey A0. Alternatively, a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point r1 of wavefield B1 and the shot location Aj from survey A0. - Per
step 716, the source and sensor definitions of wavefield A0 may be interchanged to obtain virtual wavefield A0, and source and sensor definitions of wavefield B1 may be interchanged to obtain virtual wavefield B1.FIG. 8E depicts virtual wavefields after reversal of source and sensor definitions ofFIG. 8D in accordance with an embodiment of the invention. As depicted inFIG. 8E , virtual wavefield A0 originates atvirtual source location 802 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, and travels upward to be detected by the first (intermediate) virtual sensor 804 at location Aj. As further depicted, virtual wavefield B1 originates atvirtual source location 802 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, travels upward to be reflected downward at reflection location r1, travels downward to be reflected upward at subsurface location P, and travels upward to be detected by the second (final) virtual sensor 806-B at location Bk. - Per step 718, the subsurface location P in the
reservoir 31 may be imaged, either using a deconvolution or cross-correlation imaging process. Comparison may then be made of the current image data of the subsurface location P generated in step 718 against the previous image data of the subsurface location P generated in step 708. The comparison may involve generating difference data by taking a difference between the current and previous image data. The resultant image data, difference data, and other pertinent data may be stored in a tangible medium so as to be provided as a geophysical data product. - In one embodiment, the geophysical analysis equipment (on vessel or on shore), which uses a computer system to process the seismic data, may be configured to deconvolve virtual wavefield B1 using virtual wavefield A0 to obtain deconvolved virtual wavefield B1/A0. The deconvolved virtual wavefield B1/A0 may be effectively used as the
imaging wavefield signal 820. As depicted inFIG. 8F , the imaging wavefield signal 820 may be considered as originating at avirtual source 822 at location Aj (which is effectively the same as, or made to correspond with, r1), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by thevirtual sensor 824 at location Bk. - In an alternate embodiment, the geophysical analysis equipment (on vessel or on shore), which uses a computer system to process the seismic data, may be configured to perform cross-correlation in order to image subsurface location P. In this embodiment, a first data signal may be the virtual data signal received at the first virtual sensor at location Aj, and a second data signal may be the virtual data signal received at the second virtual sensor at location Bk. The cross-correlation of the first data signal and the second data signal may provide the medium response (i.e., the Green's function) of the subsurface location P. The result is similar to imaging using a virtual source-receiver pair, where the
virtual source 822 and virtual receiver 824 (shown inFIG. 8F ) are located at the positions of the first and second virtual sensors, respectively. - Per step 722, after a second time gap, a third marine seismic survey (survey C) may be performed. The performance of survey C involves generating seismic signals at a third set of positions and recording a third set of received signals. Survey C may use an array of seabed sensors at the same positions as the array used in survey A. The seabed sensors in the array for survey C may be the same sensors as used for survey A, or the seabed sensors in the array for survey C may be redeployed into the same positions as the sensors used for survey A. In an exemplary implementation, survey C may use a sparser shot pattern than survey A. One exemplary shot location Ck of the third survey (survey C) is shown in
FIG. 8G . The shot location Ck of survey C is at the same location, or effectively the same location, as the shot location Ak of survey A. - Per
step 724, after separating up-going and down-going wavefields, wavefields A0 of survey A and C1 from survey C are selected. As depicted inFIG. 8G , the two wavefields A0 and C1 originate at different shot locations and end at the same sensor position. The selected wavefield C1 originates at shot location Ck from survey C and travels downward. Wavefield C1 reflects from subsurface point P in thereservoir 31 so as to travel upward, then reflects from point r1 of thesea surface 11 so as to travel downward. Wavefield C1 may then reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 again to travel downward until it is detected by a seabed sensor at seabed position Si. The selected wavefield A0 originates at shot location Aj and travels downward. Wavefield A0 may reflect from thesea floor 25 so as to travel upward, then reflect from thesea surface 11 to travel downward until it is detected by the seabed sensor at the same seabed position Si as that which detected wavefield C1. Note that the first sea-surface reflection point r1 of wavefield C1 may be located at shot location Aj from survey A0. Alternatively, a wavefield interpolation process may be performed to obtain a correspondence between the first sea-surface reflection point r1 of wavefield C1 and the shot location Aj from survey A0. - Per
step 726, the source and sensor definitions of wavefield A0 may be interchanged to obtain virtual wavefield A0, and source and sensor definitions of wavefield C1 may be interchanged to obtain virtual wavefield C1.FIG. 8H depicts virtual wavefields after reversal of source and sensor definitions ofFIG. 8G in accordance with an embodiment of the invention. As depicted inFIG. 8H , virtual wavefield A0 originates atvirtual source location 802 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, and travels upward to be detected by the first (intermediate) virtual sensor 804 at location Aj. As further depicted, virtual wavefield C1 originates atvirtual source location 802 at location Si, travels upward to be reflected downward by thesea surface 11, travels downward until reflected upward by thesea floor 25, travels upward to be reflected downward at reflection location r1, travels downward to be reflected upward at subsurface location P, and travels upward to be detected by the second (final) virtual sensor 806-C at location Ck. - Per
step 728, the subsurface location P in thereservoir 31 may be imaged, either using a deconvolution or cross-correlation imaging process. Comparison may then be made of the current image data of the subsurface location P generated instep 728 against the previous image data of the subsurface location P generated in step 708 or 718. The comparison may involve generating difference data by taking a difference between the current and previous image data. The resultant image data, difference data, and other pertinent data may be stored in a tangible medium so as to be provided as a geophysical data product. - In one embodiment, the geophysical analysis equipment (on vessel or on shore), which uses a computer system to process the seismic data, may be configured to deconvolve virtual wavefield C1 using virtual wavefield A0 to obtain deconvolved virtual wavefield C1/A0. The deconvolved virtual wavefield C1/A0 may be effectively used as the
imaging wavefield signal 830. As depicted inFIG. 8I , the imaging wavefield signal 830 may be considered as originating at avirtual source 832 at location Aj (which is effectively the same as, or made to correspond with, r1), traveling downward to be reflected upward at the subsurface location P, and traveling upward to be detected by thevirtual sensor 834 at location Ck. - In an alternate embodiment, the geophysical analysis equipment (on vessel or on shore), which uses a computer system to process the seismic data, may be configured to perform cross-correlation in order to image subsurface location P. In this embodiment, a first data signal may be the virtual data signal received at the first virtual sensor at location Aj, and a second data signal may be the virtual data signal received at the second virtual sensor at location Ck. The cross-correlation of the first data signal and the second data signal may provide the medium response (i.e., the Green's function) of the subsurface location P. The result is similar to imaging using a virtual source-receiver pair, where the
virtual source 832 and virtual receiver 834 (shown inFIG. 8I ) are located at the positions of the first and second virtual sensors, respectively. - The
method 700 ofFIG. 7 may continue with the performance and processing of further surveys (D, E, F, etc.). Each survey and the subsequent processing steps results in a “snapshot” image of the reservoir being monitored. - Note that the time gap between survey A and survey B may be typically of sufficient length for a change to be expected in the state of the
reservoir 31. For example, the time gap between surveys A and B may be three months, six months, or a year, or longer. Similarly, the time gap between survey B and survey C may be typically of sufficient length for a change to be expected in the state of thereservoir 31. For example, the time gap between surveys B and C may be three months, six months, or a year, or longer. And so on for further surveys to monitor a hydrocarbon reservoir over a span of years. - The above-described
method 700 ofFIG. 7 may be used to efficiently monitor a production reservoir (i.e. perform a 4D survey). The above-describedmethod 700 ofFIG. 7 provides greater flexibility for subsequent shooting monitor surveys. - In one use case, while the first survey (survey A) may have a dense shot pattern for detailed 3D imaging, the subsequent surveys (survey B1, survey C, etc.) may have a sparser shot pattern. In a specific implementation, the first survey may be recorded using a dense single shot survey, while the subsequent surveys may be recorded using cost-effective dual-source flip-flop acquisition.
- By using the dense shooting geometry of the base-line survey A as the “forward extrapolated” wavefield A0 in combination of different recorded source wavefields A1, B1, C1, etc. as the “backward extrapolated” wavefield, the difference in the images may be used to provide a 4D signal corresponding to reservoir changes within the time gaps between acquisitions. A time-lapse imaging sequence may be defined as: Base (A1/A0); Monitor1 (B1/A0); Monitor2 (C1/A0), etc.
-
FIG. 9 shows a simplified example of acomputer apparatus 40 which may be used in performing processing steps in accordance with an embodiment of the invention. Thecomputer apparatus 40 may be configured with executable instructions so as to perform the data processing methods described herein. This figure shows just one simplified example of a computer which may be used to perform the data processing methods described herein. Other types of computers may also be employed. As a large amount of seismic data is typically processed, a high-performance cluster of multi-processor server computers may be used in practice. The processing methods may also be performed using a cloud computing system. - As depicted in
FIG. 9 , the simplified example of thecomputer apparatus 40 may include aprocessor 41, such as those from the Intel Corporation of Santa Clara, Calif., for example. The computer apparatus. 40 may have a bus system 43 communicatively interconnecting its various components. Thecomputer apparatus 40 may include one or more user input devices 42 (e.g., keyboard, mouse), a display monitor 44 (e.g., LCD, flat panel monitor, CRT), a computer network interface 45 (e.g., network adapter, modem), and a data storage system which may include one or more data storage devices 46 (e.g., hard drive, solid state memory, optical disk drive, USB memory) and a main memory 50 (e.g., RAM). - In the example shown in this figure, the
main memory 50 includesexecutable code 52 anddata 54 stored therein. Theexecutable code 52 may comprise computer-readable program code (i.e., software) components which may be loaded from thedata storage device 46 to themain memory 50 for execution by theprocessor 41. In particular, theexecutable code 52 may be configured to perform computer-implemented steps in the methods described herein. - Although specific embodiments have been described above, these embodiments are not intended to limit the scope of the present disclosure, even where only a single embodiment is described with respect to a particular feature. Examples of features provided in the disclosure are intended to be illustrative rather than restrictive unless stated otherwise. The above description is intended to cover such alternatives, modifications, and equivalents as would be apparent to a person skilled in the art having the benefit of this disclosure.
- The scope of the present disclosure includes any feature or combination of features disclosed herein (either explicitly or implicitly), or any generalization thereof, whether or not it mitigates any or all of the problems addressed herein. Accordingly, new claims may be formulated during prosecution of this application (or an application claiming priority thereto) to any such combination of features. In particular, with reference to the appended claims, features from dependent claims may be combined with those of the independent claims and features from respective independent claims may be combined in any appropriate manner and not merely in the specific combinations enumerated in the appended claims.
Claims (23)
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/166,902 US10317551B2 (en) | 2015-06-01 | 2016-05-27 | Using seabed sensors and sea-surface reflections for structural imaging of a subsurface location in a geological formation |
US16/391,782 US10962666B2 (en) | 2015-06-01 | 2019-04-23 | Using seabed sensors and sea-surface reflections for structural imaging of a subsurface location in a geological formation |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201562168988P | 2015-06-01 | 2015-06-01 | |
US15/166,902 US10317551B2 (en) | 2015-06-01 | 2016-05-27 | Using seabed sensors and sea-surface reflections for structural imaging of a subsurface location in a geological formation |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US16/391,782 Division US10962666B2 (en) | 2015-06-01 | 2019-04-23 | Using seabed sensors and sea-surface reflections for structural imaging of a subsurface location in a geological formation |
Publications (2)
Publication Number | Publication Date |
---|---|
US20160349392A1 true US20160349392A1 (en) | 2016-12-01 |
US10317551B2 US10317551B2 (en) | 2019-06-11 |
Family
ID=57398442
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/166,902 Active 2036-11-25 US10317551B2 (en) | 2015-06-01 | 2016-05-27 | Using seabed sensors and sea-surface reflections for structural imaging of a subsurface location in a geological formation |
US16/391,782 Active 2036-10-17 US10962666B2 (en) | 2015-06-01 | 2019-04-23 | Using seabed sensors and sea-surface reflections for structural imaging of a subsurface location in a geological formation |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US16/391,782 Active 2036-10-17 US10962666B2 (en) | 2015-06-01 | 2019-04-23 | Using seabed sensors and sea-surface reflections for structural imaging of a subsurface location in a geological formation |
Country Status (1)
Country | Link |
---|---|
US (2) | US10317551B2 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO20170155A1 (en) * | 2017-02-01 | 2018-05-07 | Polarcus Dmcc | Source array configuration for repeated marine seismic surveying of the same area |
US11650343B2 (en) | 2019-04-17 | 2023-05-16 | Pgs Geophysical As | Directional designature of marine seismic survey data |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4970697A (en) * | 1989-10-06 | 1990-11-13 | Amoco Corporation | Vertical marine seismic array |
US8559267B2 (en) * | 2006-10-26 | 2013-10-15 | Schlumberger Technology Corporation | Methods and apparatus of borehole seismic surveys |
US7426438B1 (en) * | 2007-03-16 | 2008-09-16 | Westerngeco L.L.C. | Technique to provide seismic data-based products |
US9146329B2 (en) * | 2008-12-17 | 2015-09-29 | Exxonmobil Upstream Research Company | System and method for reconstruction of time-lapse data |
US8756042B2 (en) * | 2010-05-19 | 2014-06-17 | Exxonmobile Upstream Research Company | Method and system for checkpointing during simulations |
US9360576B2 (en) * | 2012-04-27 | 2016-06-07 | Pgs Geophysical As | Methods and apparatus for generating deghosted seismic data |
US9007870B2 (en) | 2012-05-31 | 2015-04-14 | Pgs Geophysical As | Seismic surveying techniques with illumination areas identifiable from primary and higher-order reflections |
US9684085B2 (en) * | 2013-01-15 | 2017-06-20 | Cgg Services Sas | Wavefield modelling and 4D-binning for seismic surveys from different acquisition datums |
US9651695B2 (en) * | 2013-09-19 | 2017-05-16 | Pgs Geophysical As | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields |
US9689999B2 (en) | 2014-06-13 | 2017-06-27 | Pgs Geophysical As | Seismic imaging using higher-order reflections |
-
2016
- 2016-05-27 US US15/166,902 patent/US10317551B2/en active Active
-
2019
- 2019-04-23 US US16/391,782 patent/US10962666B2/en active Active
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO20170155A1 (en) * | 2017-02-01 | 2018-05-07 | Polarcus Dmcc | Source array configuration for repeated marine seismic surveying of the same area |
NO342332B1 (en) * | 2017-02-01 | 2018-05-07 | Polarcus Dmcc | Source array configuration for repeated marine seismic surveying of the same area |
WO2018143820A1 (en) * | 2017-02-01 | 2018-08-09 | Polarcus Dmcc | Source array configuration for repeated marine seismic surveying of the same area |
GB2571469A (en) * | 2017-02-01 | 2019-08-28 | Polarcus Dmcc | Source array configuration for repeated marine seismic surveying of the same area |
GB2571469B (en) * | 2017-02-01 | 2021-08-11 | Polarcus Dmcc | Source array configuration for repeated marine seismic surveying of the same area |
US11650343B2 (en) | 2019-04-17 | 2023-05-16 | Pgs Geophysical As | Directional designature of marine seismic survey data |
Also Published As
Publication number | Publication date |
---|---|
US10317551B2 (en) | 2019-06-11 |
US10962666B2 (en) | 2021-03-30 |
US20190250296A1 (en) | 2019-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9030910B2 (en) | Jointly interpolating and deghosting seismic data | |
EP2530491B1 (en) | Methods and apparatus for seismic exploration using pressure changes caused by sea-surface variations | |
US10520623B2 (en) | Methods and systems for marine survey acquisition | |
US11604297B2 (en) | Highly-sparse seabed acquisition designs adapted for imaging geological structure and/or monitoring reservoir production | |
US9689999B2 (en) | Seismic imaging using higher-order reflections | |
AU2009257824B2 (en) | System and technique to determine high order derivatives from seismic sensor data | |
US20130028048A1 (en) | Methods and apparatus for seismic imaging which accounts for sea-surface variations | |
US20140297190A1 (en) | Monitoring of source signature directivity in seismic systems | |
US10962666B2 (en) | Using seabed sensors and sea-surface reflections for structural imaging of a subsurface location in a geological formation | |
US20140172308A1 (en) | Repeatability indicator based on shot illumination for seismic acquisition | |
NO347442B1 (en) | A seismic sensor unit, a system for seismic surveying and a method for processing seismic data | |
US9798025B2 (en) | Regularization of multi-component seismic data | |
US20140165694A1 (en) | Methods and systems for quality control of seismic illumination maps | |
US20160363682A1 (en) | Efficient survey design to in-fill a pre-existing plural vessels seismic survey | |
US11287542B2 (en) | Seismic data processing including surface multiple modeling for ocean bottom cable systems | |
GB2539097A (en) | Highly-sparse seabed acquisition designs adapted for imaging geological structure and/or monitoring reservoir production | |
GB2529038A (en) | Seismic imaging using higher-order reflections |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: PGS GEOPHYSICAL AS, NORWAY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:LECERF, DIDIER B.P.;REEL/FRAME:038947/0946 Effective date: 20160617 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 4 |