WO2014110547A1 - Seismic data processing - Google Patents
Seismic data processing Download PDFInfo
- Publication number
- WO2014110547A1 WO2014110547A1 PCT/US2014/011407 US2014011407W WO2014110547A1 WO 2014110547 A1 WO2014110547 A1 WO 2014110547A1 US 2014011407 W US2014011407 W US 2014011407W WO 2014110547 A1 WO2014110547 A1 WO 2014110547A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- packets
- model
- seismic
- refraction
- received
- Prior art date
Links
- 238000012545 processing Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 107
- 238000003325 tomography Methods 0.000 claims abstract description 53
- 239000004215 Carbon black (E152) Substances 0.000 claims abstract description 16
- 229930195733 hydrocarbon Natural products 0.000 claims abstract description 16
- 150000002430 hydrocarbons Chemical class 0.000 claims abstract description 16
- 238000004519 manufacturing process Methods 0.000 claims abstract description 15
- 230000009189 diving Effects 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 abstract description 20
- 230000006870 function Effects 0.000 description 25
- 230000000875 corresponding effect Effects 0.000 description 24
- 230000010365 information processing Effects 0.000 description 17
- 230000008569 process Effects 0.000 description 12
- 230000003287 optical effect Effects 0.000 description 8
- 238000004891 communication Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 230000006855 networking Effects 0.000 description 3
- 230000002093 peripheral effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000005055 memory storage Effects 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 238000010793 Steam injection (oil industry) Methods 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005474 detonation Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000007723 transport mechanism Effects 0.000 description 1
- 230000005641 tunneling Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/32—Transforming one recording into another or one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
-
- 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Definitions
- Seismic exploration may utilize a seismic energy source to generate acoustic signals that propagate into the earth and partially reflect off subsurface seismic reflectors (e.g., interfaces between subsurface layers).
- the reflected signals are recorded by sensors (e.g., receivers or geophones located in seismic units) laid out in a seismic spread covering a region of the earth's surface.
- the recorded signals may then be processed to yield a seismic survey.
- seismic depth- migration projects may involve many iterations of model refinement before arriving at a final model. This refinement process is time-consuming and costly, because each iteration may include a complete remigration of the whole prestack seismic data volume, followed by interpretation of the changes in the image, and updating of the model to account for the changes.
- a method for processing seismic data corresponding to a region of interest may receive the seismic data.
- the method may separate the received seismic data into refraction packets and reflection packets.
- the method may receive a model for the region of interest.
- the method may update a portion of the received model using the refraction packets with refraction traveltime tomography.
- the method may use the updated model to facilitate hydrocarbon exploration or production.
- a method for processing seismic data corresponding to a region of interest may receive the seismic data.
- the method may decompose the received seismic data into seismic packets.
- the method may separate the seismic packets into refraction packets and reflection packets.
- the method may receive a model that describes the region of interest.
- the method may update a first portion of the received model using the refraction packets with refraction traveltime tomography.
- the method may update a second portion of the received model using the reflection packets with reflection tomography.
- the method may use an updated model based on the updated first portion and second portion of the received model to facilitate hydrocarbon exploration or production.
- a method for processing seismic data corresponding to a region of interest may separate the data into refraction packets and reflection packets.
- the method may receive a model that describes the region of interest.
- the method may update a first portion of the received model using the refraction packets with refraction traveltime tomography.
- the method may update a second portion of the received model using the reflection packets with reflection tomography.
- Figure 1 illustrates a diagrammatic view of marine seismic surveying in accordance with various implementations described herein.
- Figure 2 illustrates a flow diagram of a method for processing seismic data in accordance with various implementations described herein.
- Figure 3 illustrates ray tracing in accordance with various implementations described herein.
- Figure 4 illustrates a computer system in which the various technologies and techniques described herein may be incorporated and practiced.
- first, second, etc. may be used herein to describe various elements, these elements should not be limited by these terms. These terms are used to distinguish one element from another.
- a first object or block could be termed a second object or block, and, similarly, a second object or block could be termed a first object or block, without departing from the scope of the invention.
- the first object or block, and the second object or block are both objects or blocks, respectively, but they are not to be considered the same object or block.
- the term “if may be construed to mean “when” or “upon” or “in response to determining” or “in response to detecting,” depending on the context.
- the phrase “if it is determined” or “if [a stated condition or event] is detected” may be construed to mean “upon determining” or “in response to determining” or “upon detecting [the stated condition or event]” or “in response to detecting [the stated condition or event],” depending on the context.
- FIG. 1 illustrates a diagrammatic view of marine seismic surveying 10 in connection with implementations of various techniques described herein.
- a marine seismic acquisition system 10 may include a vessel 1 1 carrying control components and towing a plurality of seismic sources 16 and a plurality of streamers 18 equipped with seismic receivers 21 .
- the seismic sources 16 may include a single type of source, or different types.
- the sources may use any type of seismic generator, such as air guns, water guns, steam injection sources, controllable seismic sources, explosive sources such as dynamite or gas injection followed by detonation and the like.
- the streamers 18 may be towed by means of their respective lead-ins 20, which may be made from high strength steel or fiber-reinforced cables that convey electrical power, control, and data signals between the vessel 1 1 and the streamers 18.
- An individual streamer may include a plurality of seismic receivers 21 that may be distributed at spaced intervals along the streamer's length.
- the seismic receivers 21 may include hydrophone sensors as well as multi-component sensor devices, such as accelerometers.
- the streamers 18 may include a plurality of inline streamer steering devices (SSDs), also known as "birds.”
- SSDs may be distributed at appropriate intervals along the streamers 18 for controlling the streamers' depth and lateral movement.
- a single survey vessel may tow a single receiver array along individual sail lines, or a plurality of survey vessels may tow a plurality of receiver arrays along a corresponding plurality of the sail lines.
- the seismic sources 16 and the seismic streamers 18 may be deployed from the vessel 1 1 and towed slowly to traverse a region of interest.
- the seismic sources 16 may be periodically activated to emit seismic energy in the form of an acoustic or pressure wave through the water.
- the sources 16 may be activated individually or substantially simultaneously with other sources.
- the acoustic wave may result in one or more wavefields that travel coherently into the earth E underlying the water W.
- the wavefields strike interfaces 4 between earth formations, or strata, they may be reflected back through the earth E and water W along paths 5 to the various receivers 21 where the wavefields (e.g., pressure waves in the case of air gun sources) may be converted to electrical signals, digitized and transmitted to the integrated computer-based seismic navigation, source controller, and recording system in the vessel 1 1 via the streamers 18 and lead-ins 20.
- the wavefields e.g., pressure waves in the case of air gun sources
- the integrated computer-based seismic navigation, source controller, and recording system in the vessel 1 1 via the streamers 18 and lead-ins 20 Through analysis of these detected signals, it may be possible to determine the shape, position and lithology of the sub-sea formations, including those formations that may include hydrocarbon deposits.
- Figure 2 illustrates a flow diagram of a method for processing seismic data in accordance with various implementations described herein. It should be understood that while the operational flow diagram indicates a particular order of execution of the operations, in other implementations, the operations might be executed in a different order. Further, in some implementations, additional operations or blocks may be added to the method. Likewise, some operations or blocks may be omitted.
- seismic data are received for a region of interest (i.e., "the received seismic data").
- the seismic data may include a seismic volume obtained from a seismic survey.
- the seismic data may be from a common shot gather.
- the region of interest may include an area of the subsurface in the earth that may be of particular interest, such as for hydrocarbon production.
- an initial model may be received for the region of interest (i.e., "the received model").
- the initial model may be a velocity model or an anisotropic model that describes the region of interest.
- One such initial model may be obtained by using intercept times and local slopes based on a Herglotz-Wiechert inversion of the received seismic data from block 210.
- the received seismic data are decomposed into seismic packets.
- a seismic dataset may be viewed as a superposition of individual seismic attributes.
- a seismic packet may include data regarding these seismic attributes.
- the seismic packets may have the capability to recombine and approximately produce the original dataset.
- the decomposed seismic packets at block 230 may be recombined to produce, approximately, the received seismic data at block 210.
- a seismic packet may include beam components as well as direct and derived attributes of the received seismic data.
- Direct attributes may include seismic amplitudes, spatial locations of a beam center (e.g., x, y and z locations), acquisition source-receiver offsets, acquisition azimuths and coherency for seismic traces.
- Derived attributes may include the reflection angle, the reflection azimuth, a source-side time dip, a receiver-side time dip, a wavelet stretch and a beam spread of a seismic wave.
- the seismic packets may be locally coherent wave packets with compact support in both space and time, and expressed using the following equation:
- D (s> g>t) ⁇ 5 i (s,g,t ⁇ gi , , p ri ) Equation 1
- D(s,g,t) is the received seismic data
- S i ⁇ s, g,t ⁇ g i ,t i ,p ri ) is an i th seismic packet at center time t,, center source s,, and center receiver g, with a receiver-side time dip p r i (also called the time slope vector).
- the receiver-side time dip p ri may describe the direction that a wave approaches receiver g,.
- the receiver-side time dip p r i may include x-component p x or y-component p y of the approaching wave.
- Equation 1 the receiver-side time dip p r i may include x-component p x or y-component p y of the approaching wave.
- the seismic packets may be described with respect to the source wavelet and expressed using the following equation:
- the seismic packets from block 230 may be separated into refraction packets and reflection packets. For instance, refraction packets may be separated out from the seismic packets and sent to block 250. Reflection packets may be separated out of the seismic packets and sent to block 260.
- Both refraction packets and reflections packets may include the same type of information as the seismic packets from block 230.
- different types of packets may refer to different types of pressure waves that pass through the earth's subsurface.
- the refraction packets may describe information regarding a refractive wave from a source to a receiver, where the path of the wave is not changed by the wave's reflection at an impedance interface between the source and the receiver.
- refractive waves may include early-arrivals (e.g., diving waves, head waves or wide-angle reflections) that may undergo continuous refractions throughout the subsurface.
- Refracted waves may be strong and coherent while receiving little influence by multiples, shear or converted waves, which may make them very suitable for velocity analysis.
- refracted waves may be associated with long wave paths and large apertures, which may be used to obtain an initial model of long wavelength.
- the reflection packets may describe information regarding a reflected wave that travels from a source to a receiver, and where the reflecting wave's path is changed by one or more reflections at impedance interfaces between the source and the receiver.
- the reflected wave may exit an interface at the same angle as the angle of incidence, which may alter the reflected wave's path through the region of interest.
- refraction packets or reflection packets may be separated from the seismic packets using source to receiver traveltimes. For instance, a time window may be used to isolate the refraction packets from the rest of the seismic packets. As such, the time window may filter out seismic packets that arrive inside a predetermined time.
- the predetermined time may correspond to a designated time period when the majority of early-arrival waves have reached a corresponding receiver.
- the seismic packets that include acquisition times outside the designated time period may be determined to be a refraction packet, while those seismic packets with acquisition times inside the designated time period may be determined to be reflection packets or other types of data.
- source-receiver offsets may be used in the separation process at block 240 in a similar manner to the time window described above. For instance, a distance window may be used where seismic packets that include source-receiver offsets outside the distance window may be considered refraction packets and those inside the distance window may be considered reflection packets. [0034] At block 250, the refraction packets may be used to perform refraction traveltime tomography.
- the modeled traveltime t cal may be estimated by modeling or backprojecting a ray from the receiver at g using the receiver-side time dip p ri and searching for one point X along a corresponding ray path back to the source at s.
- Figure 3 illustrates an example of ray tracing that may be used in block 250.
- a ray path 350 between receivers 310 and point 330 may be modeled using the received model from block 220. Based on the time dip 335 recorded at receivers 310, the ray path 350 may be modeled backwards to the source 320 using the observed traveltime between the firing of the source 320 and the arrival of the seismic wave at the receivers 310. For instance, modeled ray path 360 using a calculated traveltime from the point 330 to the source 320 misses the source 320. Ray tracing may then be repeated until a ray path 340 is produced that reaches the source 320 from point 330.
- the early-arrival traveltime tomographic algorithm may take into account early-arrivals in addition to the first break arrivals, which may allow for no explicit selecting of first break arrivals. This may be advantageous as compared to known diving wave tomography.
- the refraction packets may be used to model a ray path from a source to a receiver.
- t cal is the total modeled traveltime
- t sx is the traveltime from the source to image point X
- t gx is the traveltime from the receiver to the image point X.
- t sx may be computed using a traveltime table and t gx may be computed by ray tracing one ray.
- a portion of the received model is updated based on the refraction traveltime tomography from block 250. Similar to other methods of tomography, updating the portion of the received model may include solving the refraction traveltime tomography problem by iteratively minimizing a cost function minimization criterion. This update may involve internal nonlinear iterations. As such, updating the received model may be performed iteratively by calculating new values for the received model and using those new values at block 250 if the updated model has not converged to a cost function minimization criterion at block 258. Thus, the gradient field for a cost function regarding the received model's parameters may be calculated and a conjugate gradient method may be used to update the received model.
- a cost function minimization criterion may designate when the error in an updated model approaches a specific degree of accuracy. When the updated model converges to a cost function minimization criterion, the updated model may be determined as accurate for the region of interest. If two refraction packets are for the same event, the packets' attributes may be highly correlated.
- a cost function minimization criterion may be the difference between a forward-modeled prediction of the region of interest and the received seismic data or the data in the refractions packets from block 240.
- One instance of a cost function may be shown by the following equation:
- Equation 4 the first order perturbation of the cost function shown in Equation 4 may be provided by the following equation:
- the process may return to block 250 to repeat refraction traveltime tomography using the updated model, an internal nonlinear iteration. If the updated model has converged to a corresponding cost function minimization criterion, the process may proceed to block 280 to generate a final model for the region of interest. In one implementation, the updated model may provide a shallow depth model for the region of interest. The process may also proceed to block 260 to further update the updated model using the reflection packets and reflection tomography.
- Blocks 260-278 may describe several different process flows in regard to the reflection packets from block 240.
- the process flow from blocks 250-258 may be performed simultaneously with the process flow involving reflection tomography in blocks 260-278 in order to generate a final model at block 280.
- the updated model produced by blocks 250-258 may be used for performing the reflection tomography at blocks 260-278.
- one or more blocks between 260-278 may be excluded from the process.
- the reflection packets from block 240 may be migrated into the depth domain (i.e., "migrated packets").
- the received seismic data at block 210 may include data in the time domain that corresponds to the acquisition time of the seismic data during a survey.
- the total traveltime image condition may be expressed by the following equation:
- t sjt is the traveltime from the source to an image point X
- t g is the traveltime from the image point to the receivers.
- a migrated packet may have compact support near the image point X.
- a seismic packet in the depth domain A i (x ⁇ X i , i ,0 i ⁇ i ,s i , g i ) may be determined using the following equation:
- r s i g denotes the total traveltime from source s, though the image point x i and then back to the receiver g h ⁇ ⁇ is the reflection angle, ⁇ ⁇ is the subsurface azimuth, v i represents the velocity at the packet center . along the normal direction n and ⁇ , ⁇ ⁇ ) denotes a 3D tapering function.
- subsurface information may be obtained regarding the migrated energy such as dips, azimuths and reflection angles, as well as associated surface information about the unmigrated data like source and receiver locations, source-side or receiver-side local time dips (e.g., dt/dsx, dt/dsy, dt/dgx, dt/dgy, dt/dhx and dt/dhy), coherencies and wavelets.
- source-side or receiver-side local time dips e.g., dt/dsx, dt/dsy, dt/dgx, dt/dgy, dt/dhx and dt/dhy
- the migrated packets are analyzed for whether a flat common image gather is produced. If a resulting common image gather from the migrated packets is flat, the migrated packets may proceed to block 280 for generating a final model. If the common image gather is not flat, the process may proceed to block 264, an internal nonlinear iteration.
- the migrated packets are sorted and binned into different groups based on common attributes. For instance, the migrated packets may be sorted according to offsets, reflection angles, or other packet attributes. Then, the sorted migrated packets may be binned into different groups such that binned packets are for the same event (e.g., for the same fault or the same reflection at a reflection interface). For instance, migrated packets may be binned into a group based on the following criteria: ⁇ - ⁇ ⁇ K - ) Equation 9
- i 3 ⁇ 4 and t h denote the traveltimes for two migrated packets for the same event with offsets ⁇ and h 2 , respectively.
- the source wavelets for migrated packets in the same group may be highly correlated, which may serve as a grouping criteria to avoid cycle skipping.
- Migrated packets within a group may also be selected such that the packets within a group are located within a predetermined distance from the corresponding event.
- User input may also be used to select or mask certain types of packets from the sorting or binning process, such as those packets associated with faults, steep events or a salt over-hang.
- a common image packet gather may be generated from the migrated packets.
- a common image packet gather may be formed using migrated packets that contribute to the same image location (i.e., the same location in the subsurface of the region of interest). If the migrated packets are resorted according to offsets or reflection angles, and the binned packets are reconstructed from the result, then a common image offset gather or a common image angle gathers may be obtained.
- two migrated packets that are in a common image packet gather may be represented as ⁇ ; and ⁇ ; .
- the two migrated packets may be compared in the form of a differential semblance as described by the following equation:
- Sij - ⁇ ⁇ ⁇ i ⁇ i ⁇ i , ⁇ P i , s i ,g i ) j (x ⁇ X j , n j ,e j ⁇ j , s j ,g j ) Equation 10
- Equation 10 may define the correlation between the two migrated packets.
- two migrated packets for the same event may be highly correlated.
- the differential semblance for the two migrated packets may be expressed by the following equation:
- Equation 1 1 Equation 1 1 where W i and W. are 3D tapering functions for the two migrated packets.
- the migrated packets or the reflection packets from block 250 may be used to perform reflection tomography. For instance, ray tracing may be done in the depth-domain similar to the ray path modeling done in the time-domain as described in regard to Figure 3 above.
- a portion of the received model from block 220, the updated model from block 258, or a model previously updated at block 274 may be updated based on the reflection tomography from block 270. Updating this portion may include solving the reflection tomography problem by minimizing a cost function minimization criterion in a similar manner to the approach used at block 254.
- Equation 12 Equation 12
- E l is the cost function for early-arrival traveltime tomography, such as the one expressed in Equation 4, ⁇ denotes a weight coefficient for early-arrival traveltime misfits, and s m) denotes the mismatch between any two migrated packets for I th event in k th common image gather.
- the indices of k and / are omitted during the following derivations since the migrated packets may be assumed to correspond for the same events in a common image packet gather.
- Equation 13 tapering function perturbations may be ignored, and w t and Wj are time-domain source wavelets for the two migrated packets described at block 268, respectively.
- Equation 15 may provide the Frechet derivative for a conjugated gradient solver.
- the process may return to block 270 to repeat reflection tomography with the updated model or to block 264 to repeat the sorting and binning of the migrated packets for generating the common image packet gathers. If the updated model has converged to a corresponding cost function minimization criterion, the process may proceed to block 280 to generate a final model for the region of interest.
- a final model for the region of interest may be generated from the updated model from block 258 or block 278.
- the final model may be an anisotropy model for the region of interest.
- the final model may also be a prestack depth migration (PSDM) volume or PSDM gather for the region of interest.
- PSDM prestack depth migration
- the final model may be used to facilitate hydrocarbon exploration or production in the region of interest.
- a method for processing seismic data corresponding to a region of interest may receive the seismic data.
- the method may separate the received seismic data into refraction packets and reflection packets.
- the method may receive a model for the region of interest.
- the method may update a first portion of the received model using the refraction packets with refraction traveltime tomography.
- the method may use the updated model to facilitate hydrocarbon exploration or production.
- the method may update a second portion of the received model using the reflection packets with reflection tomography.
- the method may also migrate the reflection packets into the depth domain.
- the received seismic data may also be from a common shot gather.
- the method may also decompose the received seismic data into a plurality of seismic packets, where the seismic packets may be separated into refraction packets and reflection packets. Updating the first portion of the received model may be performed iteratively.
- the received seismic data may include attributes such as a spatial location of a beam center, a spatial orientation of a beam, a source-receiver offset, a source-receiver azimuth, a reflection angle, a reflection azimuth, a wavelet identification, an amplitude of a seismic wave, a coherency for seismic traces, and a beam spread.
- One of the refraction packets may include the time dip of a seismic ray measured at a receiver.
- the method may also model one or more ray paths from the receiver to a source or vice versa using the time dip.
- the method may also model the traveltime of the ray paths from the receiver to the source or vice versa.
- the method may separate the received seismic data based on source to receiver traveltimes or source-receiver offsets.
- the refraction packets may describe information regarding a pressure wave that travels from a source to a receiver, where the path of the pressure wave is not changed by a reflection of the pressure wave at an impedance interface between the source and the receiver.
- the refraction packets may also describe information regarding a diving wave or a head wave.
- the reflection packets may describe information regarding a pressure wave that travels from a source to a receiver, and where the pressure wave is reflected from an impedance interface.
- the received model or the updated model may be a velocity model or an anisotropic model for the region of interest. Updating the first portion of the received model may include minimizing a cost function minimization criterion.
- the cost function minimization criterion may be a difference between a forward-modeled prediction and the received seismic data.
- an information processing apparatus for use in a computing system, and includes means for receiving seismic data corresponding to a region of interest.
- the information processing apparatus may also have means for separating the received seismic data into refraction packets and reflection packets.
- the information processing apparatus may also have means for receiving a model for the region of interest.
- the information processing apparatus may also have means for updating a portion of the received model using the refraction packets with refraction traveltime tomography.
- the information processing apparatus may also have means for using the updated model to facilitate hydrocarbon exploration or production.
- a computing system includes at least one processor, at least one memory, and one or more programs stored in the at least one memory, wherein the programs include instructions, which when executed by the at least one processor cause the computing system to receive seismic data corresponding to a region of interest.
- the programs may further include instructions to cause the computing system to separate the received seismic data into refraction packets and reflection packets.
- the programs may further include instructions to cause the computing system to receive a model for the region of interest.
- the programs may further include instructions to cause the computing system to update a portion of the received model using the refraction packets with refraction traveltime tomography.
- the programs may further include instructions to cause the computing system to use the updated model to facilitate hydrocarbon exploration or production.
- a computer readable storage medium which has stored therein one or more programs, the one or more programs including instructions, which when executed by a processor, cause the processor to receive seismic data corresponding to a region of interest.
- the programs may further include instructions, which cause the processor to separate the received seismic data into refraction packets and reflection packets.
- the programs may further include instructions, which cause the processor to receive a model for the region of interest.
- the programs may further include instructions, which cause the processor to update a portion of the received model using the refraction packets with refraction traveltime tomography.
- the programs may further include instructions, which cause the processor to use the updated model to facilitate hydrocarbon exploration or production.
- a method for processing seismic data corresponding to a region of interest may receive the seismic data.
- the method may decompose the received seismic data into a plurality of seismic packets.
- the method may separate the seismic packets into refraction packets and reflection packets.
- the method may receive a model that describes the region of interest.
- the method may update a first portion of the received model using the refraction packets with refraction traveltime tomography.
- the method may update a second portion of the received model using the reflection packets with reflection tomography.
- the method may use an updated model based on the updated first portion and second portion of the received model to facilitate hydrocarbon exploration or production.
- an information processing apparatus for use in a computing system, and includes means for receiving seismic data corresponding to a region of interest.
- the information processing apparatus may also have means for decomposing the received seismic data into a plurality of seismic packets.
- the information processing apparatus may also have means for separating the seismic packets into refraction packets and reflection packets.
- the information processing apparatus may also have means for receiving a model for the region of interest.
- the information processing apparatus may also have means for updating a first portion of the received model using the refraction packets with refraction traveltime tomography.
- the information processing apparatus may also have means for updating a second portion of the received model using the reflection packets with reflection tomography.
- the information processing apparatus may also have means for using an updated model based on the updated first portion and second portion of the received model to facilitate hydrocarbon exploration or production.
- a computing system includes at least one processor, at least one memory, and one or more programs stored in the at least one memory, wherein the programs include instructions, which when executed by the at least one processor cause the computing system to receive seismic data corresponding to a region of interest.
- the programs may further include instructions to cause the computing system to decompose the received seismic data into a plurality of seismic packets.
- the programs may further include instructions to cause the computing system to separate the seismic packets into refraction packets and reflection packets.
- the programs may further include instructions to cause the computing system to receive a model that describes the region of interest.
- the programs may further include instructions to cause the computing system to update a first portion of the received model using the refraction packets with refraction traveltime tomography.
- the programs may further include instructions to cause the computing system to update a second portion of the received model using the reflection packets with reflection tomography.
- the programs may further include instructions to cause the computing system to use an updated model based on the updated first portion and second portion of the received model to facilitate hydrocarbon exploration or production.
- a computer readable storage medium which has stored therein one or more programs, the one or more programs including instructions, which when executed by a processor, cause the processor to receive seismic data corresponding to a region of interest.
- the programs may further include instructions, which cause the processor to decompose the received seismic data into a plurality of seismic packets.
- the programs may further include instructions, which cause the processor to separate the seismic packets into refraction packets and reflection packets.
- the programs may further include instructions, which cause the processor to receive a model that describes the region of interest.
- the programs may further include instructions, which cause the processor to update a first portion of the received model using the refraction packets with refraction traveltime tomography.
- the programs may further include instructions, which cause the processor to update a second portion of the received model using the reflection packets with reflection tomography.
- the programs may further include instructions, which cause the processor to use an updated model based on the updated first portion and second portion of the received model to facilitate hydrocarbon exploration or production.
- a method for processing data corresponding to a region of interest may receive data corresponding to a region of interest.
- the method may separate the data into refraction packets and reflection packets.
- the method may receive a model that describes the region of interest.
- the method may update a first portion of the received model using the refraction packets with refraction traveltime tomography.
- the method may update a second portion of the received model using the reflection packets with reflection tomography.
- an information processing apparatus for use in a computing system, and includes means for receiving data corresponding to a region of interest.
- the information processing apparatus may also have means for separating the data into refraction packets and reflection packets.
- the information processing apparatus may also have means for receiving a model for the region of interest.
- the information processing apparatus may also have means for updating a first portion of the received model using the refraction packets with refraction traveltime tomography.
- the information processing apparatus may also have means for updating a second portion of the received model using the reflection packets with reflection tomography.
- a computing system includes at least one processor, at least one memory, and one or more programs stored in the at least one memory, wherein the programs include instructions, which when executed by the at least one processor cause the computing system to receive data corresponding to a region of interest.
- the programs may further include instructions to cause the computing system to separate the data into refraction packets and reflection packets.
- the programs may further include instructions to cause the computing system to receive a model that describes the region of interest.
- the programs may further include instructions to cause the computing system to update a first portion of the received model using the refraction packets with refraction traveltime tomography.
- the programs may further include instructions to cause the computing system to update a second portion of the received model using the reflection packets with reflection tomography.
- a computer readable storage medium which has stored therein one or more programs, the one or more programs including instructions, which when executed by a processor, cause the processor to receive data corresponding to a region of interest.
- the programs may further include instructions, which cause the processor to separate the data into refraction packets and reflection packets.
- the programs may further include instructions, which cause the processor to receive a model that describes the region of interest.
- the programs may further include instructions, which cause the processor to update a first portion of the received model using the refraction packets with refraction traveltime tomography.
- the programs may further include instructions, which cause the processor to update a second portion of the received model using the reflection packets with reflection tomography.
- Implementations of various technologies described herein may be operational with numerous general purpose or special purpose computing system environments or configurations.
- Examples of well known computing systems, environments, and/or configurations that may be suitable for use with the various technologies described herein include, but are not limited to, personal computers, server computers, hand-held or laptop devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, smartphones, smartwatches, personal wearable computing systems networked with other computing systems, tablet computers, and distributed computing environments that include any of the above systems or devices, and the like.
- program modules include routines, programs, objects, components, data structures, etc. that performs particular tasks or implement particular abstract data types. While program modules may execute on a single computing system, it should be appreciated that, in some implementations, program modules may be implemented on separate computing systems or devices adapted to communicate with one another. A program module may also be some combination of hardware and software where particular tasks performed by the program module may be done either through hardware, software, or both.
- the various technologies described herein may also be implemented in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network, e.g., by hardwired links, wireless links, or combinations thereof.
- the distributed computing environments may span multiple continents and multiple vessels, ships or boats.
- program modules may be located in both local and remote computer storage media including memory storage devices.
- FIG. 4 illustrates a schematic diagram of a computing system 400 in which the various technologies described herein may be incorporated and practiced.
- the computing system 400 may be a conventional desktop or a server computer, as described above, other computer system configurations may be used.
- the computing system 400 may include a central processing unit (CPU) 430, a system memory 426, a graphics processing unit (GPU) 431 and a system bus 428 that couples various system components including the system memory 426 to the CPU 430.
- CPU central processing unit
- GPU graphics processing unit
- the GPU 431 may be a microprocessor specifically designed to manipulate and implement computer graphics.
- the CPU 430 may offload work to the GPU 431 .
- the GPU 431 may have its own graphics memory, and/or may have access to a portion of the system memory 426.
- the GPU 431 may include one or more processing units, and the processing units may include one or more cores.
- the system bus 428 may be any of several types of bus structures, including a memory bus or memory controller, a peripheral bus, and a local bus using any of a variety of bus architectures.
- bus architectures include Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MCA) bus, Enhanced ISA (EISA) bus, Video Electronics Standards Association (VESA) local bus, and Peripheral Component Interconnect (PCI) bus also known as Mezzanine bus.
- the system memory 426 may include a read-only memory (ROM) 412 and a random access memory (RAM) 416.
- a basic input/output system (BIOS) 414 containing the basic routines that help transfer information between elements within the computing system 400, such as during start-up, may be stored in the ROM 412.
- the computing system 400 may further include a hard disk drive 450 for reading from and writing to a hard disk, a magnetic disk drive 452 for reading from and writing to a removable magnetic disk 456, and an optical disk drive 454 for reading from and writing to a removable optical disk 458, such as a CD ROM or other optical media.
- the hard disk drive 450, the magnetic disk drive 452, and the optical disk drive 454 may be connected to the system bus 428 by a hard disk drive interface 436, a magnetic disk drive interface 438, and an optical drive interface 440, respectively.
- the drives and their associated computer-readable media may provide nonvolatile storage of computer-readable instructions, data structures, program modules and other data for the computing system 400.
- computing system 400 is described herein as having a hard disk, a removable magnetic disk 456 and a removable optical disk 458, it should be appreciated by those skilled in the art that the computing system 400 may also include other types of computer-readable media that may be accessed by a computer.
- computer-readable media may include computer storage media and communication media.
- Computer storage media may include volatile and non-volatile, and removable and non-removable media implemented in any method or technology for storage of information, such as computer-readable instructions, data structures, program modules or other data.
- Computer storage media may further include RAM, ROM, erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other solid state memory technology, CD-ROM, digital versatile disks (DVD), or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by the computing system 400.
- Communication media may embody computer readable instructions, data structures, program modules or other data in a modulated data signal, such as a carrier wave or other transport mechanism and may include any information delivery media.
- modulated data signal may mean a signal that has one or more of its characteristics set or changed in such a manner as to encode information in the signal.
- communication media may include wired media such as a wired network or direct-wired connection, and wireless media such as acoustic, RF, infrared and other wireless media.
- the computing system 400 may also include a host adapter 433 that connects to a storage device 435 via a small computer system interface (SCSI) bus, a Fiber Channel bus, an eSATA bus, or using any other applicable computer bus interface. Combinations of any of the above may also be included within the scope of computer readable media.
- SCSI small computer system interface
- a number of program modules may be stored on the hard disk 450, magnetic disk 456, optical disk 458, ROM 412 or RAM 416, including an operating system 418, one or more application programs 420, program data 424, and a database system 448.
- the application programs 420 may include various mobile applications ("apps") and other applications configured to perform various methods and techniques described herein.
- the operating system 418 may be any suitable operating system that may control the operation of a networked personal or server computer, such as Windows® XP, Mac OS® X, Unix-variants (e.g., Linux® and BSD®), and the like.
- a user may enter commands and information into the computing system 400 through input devices such as a keyboard 462 and pointing device 460.
- Other input devices may include a microphone, joystick, game pad, satellite dish, scanner, or the like.
- These and other input devices may be connected to the CPU 430 through a serial port interface 442 coupled to system bus 428, but may be connected by other interfaces, such as a parallel port, game port or a universal serial bus (USB).
- a monitor 434 or other type of display device may also be connected to system bus 428 via an interface, such as a video adapter 432.
- the computing system 400 may further include other peripheral output devices such as speakers and printers.
- the computing system 400 may operate in a networked environment using logical connections to one or more remote computers 474.
- the logical connections may be any connection that is commonplace in offices, enterprise-wide computer networks, intranets, and the Internet, such as local area network (LAN) 476 and a wide area network (WAN) 466.
- the remote computers 474 may be another a computer, a server computer, a router, a network PC, a peer device or other common network node, and may include many of the elements describes above relative to the computing system 400.
- the remote computers 474 may also each include application programs 470 similar to that of the computer action function.
- the computing system 400 may be connected to the local network 476 through a network interface or adapter 444.
- the computing system 400 may include a router 464, wireless router or other means for establishing communication over a wide area network 466, such as the Internet.
- the router 464 which may be internal or external, may be connected to the system bus 428 via the serial port interface 442.
- program modules depicted relative to the computing system 400, or portions thereof, may be stored in a remote memory storage device 435. It will be appreciated that the network connections shown are merely examples and other means of establishing a communications link between the computers may be used.
- the network interface 444 may also utilize remote access technologies (e.g., Remote Access Service (RAS), Virtual Private Networking (VPN), Secure Socket Layer (SSL), Layer 2 Tunneling (L2T), or any other suitable protocol). These remote access technologies may be implemented in connection with the remote computers 474. [0084] It should be understood that the various technologies described herein may be implemented in connection with hardware, software or a combination of both.
- RAS Remote Access Service
- VPN Virtual Private Networking
- SSL Secure Socket Layer
- L2T Layer 2 Tunneling
- various technologies may take the form of program code (i.e., instructions) embodied in tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any other machine-readable storage medium wherein, when the program code is loaded into and executed by a machine, such as a computer, the machine becomes an apparatus for practicing the various technologies.
- the computing device may include a processor, a storage medium readable by the processor (including volatile and non-volatile memory and/or storage elements), at least one input device, and at least one output device.
- One or more programs that may implement or utilize the various technologies described herein may use an application programming interface (API), reusable controls, and the like.
- API application programming interface
- Such programs may be implemented in a high level procedural or object oriented programming language to communicate with a computer system.
- the program(s) may be implemented in assembly or machine language, if desired.
- the language may be a compiled or interpreted language, and combined with hardware implementations.
- the program code may execute entirely on a user's computing device, partly on the user's computing device, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or a server computer.
- processing techniques for collected data may also be used successfully with collected data types other than seismic data. While certain implementations have been disclosed in the context of seismic data collection and processing, those with skill in the art will recognize that one or more of the methods, techniques, and computing systems disclosed herein can be applied in many fields and contexts where data involving structures arrayed in a three-dimensional space and/or subsurface region of interest may be collected and processed, e.g., medical imaging techniques such as tomography, ultrasound, MRI and the like for human tissue; radar, sonar, and LIDAR imaging techniques; and other appropriate three-dimensional imaging problems.
- medical imaging techniques such as tomography, ultrasound, MRI and the like for human tissue
- radar, sonar, and LIDAR imaging techniques and other appropriate three-dimensional imaging problems.
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Described herein are implementations of various technologies for a method for processing seismic data corresponding to a region of interest. The method may receive the seismic data. The method may separate the received seismic data into refraction packets and reflection packets. The method may receive a model for the region of interest. The method may update a first portion of the received model using the refraction packets with refraction traveltime tomography. The method may use the updated model to facilitate hydrocarbon exploration or production.
Description
SEISMIC DATA PROCESSING CROSS-REFERENCE TO RELATED APPLICATIONS
[0001] This application claims benefit of United States provisional patent application serial number 61/75201 1 , entitled METHODS AND COMPUTING SYSTEMS FOR PROCESSING AND TRANSFORMING COLLECTED DATA FOR IMPROVED VISUALIZATION filed 14 January 2013 and United States non-provisional patent application serial number 14/153894, entitled SEISMIC DATA PROCESSING filed 13 January 2014, which are both herein incorporated by reference in their entirety.
BACKGROUND
[0002] This section is intended to provide background information to facilitate a better understanding of various technologies described herein. As the section's title implies, this is a discussion of related art. That such art is related in no way implies that it is prior art. The related art may or may not be prior art. It should therefore be understood that the statements in this section are to be read in this light, and applicant neither concedes nor acquiesces to the position that any given reference is prior art or analogous prior art.
[0003] Seismic exploration may utilize a seismic energy source to generate acoustic signals that propagate into the earth and partially reflect off subsurface seismic reflectors (e.g., interfaces between subsurface layers). The reflected signals are recorded by sensors (e.g., receivers or geophones located in seismic units) laid out in a seismic spread covering a region of the earth's surface. The recorded signals may then be processed to yield a seismic survey.
[0004] As seismic exploration becomes increasingly complex, the importance of analyzing seismic samples increases. In areas of complex geology, seismic depth- migration projects may involve many iterations of model refinement before arriving at a final model. This refinement process is time-consuming and costly, because each iteration may include a complete remigration of the whole prestack seismic data
volume, followed by interpretation of the changes in the image, and updating of the model to account for the changes.
[0005] Accordingly, there is a need for methods and computing systems that can employ more effective and accurate methods for identifying, isolating, and/or processing various aspects of seismic signals or other data that is collected from a subsurface region or other multi-dimensional space.
SUMMARY
[0006] In some implementations, a method for processing seismic data corresponding to a region of interest is provided. The method may receive the seismic data. The method may separate the received seismic data into refraction packets and reflection packets. The method may receive a model for the region of interest. The method may update a portion of the received model using the refraction packets with refraction traveltime tomography. The method may use the updated model to facilitate hydrocarbon exploration or production.
[0007] In some implementations, a method for processing seismic data corresponding to a region of interest is provided. The method may receive the seismic data. The method may decompose the received seismic data into seismic packets. The method may separate the seismic packets into refraction packets and reflection packets. The method may receive a model that describes the region of interest. The method may update a first portion of the received model using the refraction packets with refraction traveltime tomography. The method may update a second portion of the received model using the reflection packets with reflection tomography. The method may use an updated model based on the updated first portion and second portion of the received model to facilitate hydrocarbon exploration or production.
[0008] In some implementations, a method for processing seismic data corresponding to a region of interest is provided. The method may separate the data into refraction packets and reflection packets. The method may receive a model that
describes the region of interest. The method may update a first portion of the received model using the refraction packets with refraction traveltime tomography. The method may update a second portion of the received model using the reflection packets with reflection tomography.
[0009] The above referenced summary section is provided to introduce a selection of concepts that are further described below in the detailed description section. The summary is not intended to identify features of the claimed subject matter, nor is it intended to be used to limit the scope of the claimed subject matter. Furthermore, the claimed subject matter is not limited to implementations that solve any or most disadvantages noted in any part of this disclosure. Indeed, the systems, methods, processing procedures, techniques, and workflows disclosed herein may complement or replace conventional methods for identifying, isolating, and/or processing various aspects of seismic signals or other data that is collected from a subsurface region or other multi-dimensional space.
BRIEF DESCRIPTION OF THE DRAWINGS
[0010] Implementations of various technologies will hereafter be described with reference to the accompanying drawings. It should be understood, however, that the accompanying drawings illustrate various implementations described herein and are not meant to limit the scope of various technologies described herein.
[0011] Figure 1 illustrates a diagrammatic view of marine seismic surveying in accordance with various implementations described herein.
[0012] Figure 2 illustrates a flow diagram of a method for processing seismic data in accordance with various implementations described herein.
[0013] Figure 3 illustrates ray tracing in accordance with various implementations described herein.
[0014] Figure 4 illustrates a computer system in which the various technologies and techniques described herein may be incorporated and practiced.
DETAILED DESCRIPTION
[0015] The discussion below is directed to certain specific implementations. It is to be understood that the discussion below is for the purpose of enabling a person with ordinary skill in the art to make and use any subject matter defined now or later by the patent "claims" found in any issued patent herein.
[0016] Reference will now be made in detail to various implementations, examples of which are illustrated in the accompanying drawings and figures. In the following detailed description, numerous specific details are set forth in order to provide a thorough understanding of the claimed invention. However, it will be apparent to one of ordinary skill in the art that the claimed invention may be practiced without these specific details. In other instances, well known methods, procedures, components, circuits, and networks have not been described in detail so as not to unnecessarily obscure aspects of the claimed invention.
[0017] It will also be understood that, although the terms first, second, etc. may be used herein to describe various elements, these elements should not be limited by these terms. These terms are used to distinguish one element from another. For example, a first object or block could be termed a second object or block, and, similarly, a second object or block could be termed a first object or block, without departing from the scope of the invention. The first object or block, and the second object or block, are both objects or blocks, respectively, but they are not to be considered the same object or block.
[0018] The terminology used in the description herein is for the purpose of describing particular implementations and is not intended to limit the claimed invention. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will also be understood
that the term "and/or" as used herein refers to and encompasses any possible combinations of one or more of the associated listed items. It will be further understood that the terms "includes," "including," "comprises," and/or "comprising," when used in this specification, specify the presence of stated features, integers, blocks, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, blocks, operations, elements, components, and/or groups thereof.
[0019] As used herein, the term "if may be construed to mean "when" or "upon" or "in response to determining" or "in response to detecting," depending on the context. Similarly, the phrase "if it is determined" or "if [a stated condition or event] is detected" may be construed to mean "upon determining" or "in response to determining" or "upon detecting [the stated condition or event]" or "in response to detecting [the stated condition or event]," depending on the context.
[0020] Figure 1 illustrates a diagrammatic view of marine seismic surveying 10 in connection with implementations of various techniques described herein. A marine seismic acquisition system 10 may include a vessel 1 1 carrying control components and towing a plurality of seismic sources 16 and a plurality of streamers 18 equipped with seismic receivers 21 . The seismic sources 16 may include a single type of source, or different types. The sources may use any type of seismic generator, such as air guns, water guns, steam injection sources, controllable seismic sources, explosive sources such as dynamite or gas injection followed by detonation and the like. The streamers 18 may be towed by means of their respective lead-ins 20, which may be made from high strength steel or fiber-reinforced cables that convey electrical power, control, and data signals between the vessel 1 1 and the streamers 18. An individual streamer may include a plurality of seismic receivers 21 that may be distributed at spaced intervals along the streamer's length. The seismic receivers 21 may include hydrophone sensors as well as multi-component sensor devices, such as accelerometers. Further, the streamers 18 may include a plurality of inline streamer
steering devices (SSDs), also known as "birds." The SSDs may be distributed at appropriate intervals along the streamers 18 for controlling the streamers' depth and lateral movement. A single survey vessel may tow a single receiver array along individual sail lines, or a plurality of survey vessels may tow a plurality of receiver arrays along a corresponding plurality of the sail lines.
[0021] During acquisition, the seismic sources 16 and the seismic streamers 18 may be deployed from the vessel 1 1 and towed slowly to traverse a region of interest. The seismic sources 16 may be periodically activated to emit seismic energy in the form of an acoustic or pressure wave through the water. The sources 16 may be activated individually or substantially simultaneously with other sources. The acoustic wave may result in one or more wavefields that travel coherently into the earth E underlying the water W. As the wavefields strike interfaces 4 between earth formations, or strata, they may be reflected back through the earth E and water W along paths 5 to the various receivers 21 where the wavefields (e.g., pressure waves in the case of air gun sources) may be converted to electrical signals, digitized and transmitted to the integrated computer-based seismic navigation, source controller, and recording system in the vessel 1 1 via the streamers 18 and lead-ins 20. Through analysis of these detected signals, it may be possible to determine the shape, position and lithology of the sub-sea formations, including those formations that may include hydrocarbon deposits.
[0022] Figure 2 illustrates a flow diagram of a method for processing seismic data in accordance with various implementations described herein. It should be understood that while the operational flow diagram indicates a particular order of execution of the operations, in other implementations, the operations might be executed in a different order. Further, in some implementations, additional operations or blocks may be added to the method. Likewise, some operations or blocks may be omitted.
[0023] At block 210, seismic data are received for a region of interest (i.e., "the received seismic data"). For instance, the seismic data may include a seismic volume obtained from a seismic survey. In one implementation, the seismic data may be from
a common shot gather. The region of interest may include an area of the subsurface in the earth that may be of particular interest, such as for hydrocarbon production.
[0024] At block 220, an initial model may be received for the region of interest (i.e., "the received model"). For instance, the initial model may be a velocity model or an anisotropic model that describes the region of interest. One such initial model may be obtained by using intercept times and local slopes based on a Herglotz-Wiechert inversion of the received seismic data from block 210.
[0025] At block 230, the received seismic data are decomposed into seismic packets. Regarding decomposition, a seismic dataset may be viewed as a superposition of individual seismic attributes. As such, a seismic packet may include data regarding these seismic attributes. When a dataset is decomposed into a finite number of seismic packets, the seismic packets may have the capability to recombine and approximately produce the original dataset. In this regard, the decomposed seismic packets at block 230 may be recombined to produce, approximately, the received seismic data at block 210.
[0026] In regard to the seismic attributes, a seismic packet may include beam components as well as direct and derived attributes of the received seismic data. Direct attributes may include seismic amplitudes, spatial locations of a beam center (e.g., x, y and z locations), acquisition source-receiver offsets, acquisition azimuths and coherency for seismic traces. Derived attributes may include the reflection angle, the reflection azimuth, a source-side time dip, a receiver-side time dip, a wavelet stretch and a beam spread of a seismic wave.
[0027] Further, the seismic packets may be locally coherent wave packets with compact support in both space and time, and expressed using the following equation:
D(s> g>t) =∑ 5i (s,g,t \ gi , , pri ) Equation 1
where D(s,g,t) is the received seismic data, Si {s, g,t \ gi ,ti ,pri ) is an ith seismic packet at center time t,, center source s,, and center receiver g, with a receiver-side time dip pri (also called the time slope vector). The receiver-side time dip pri may describe the direction that a wave approaches receiver g,. As such, the receiver-side time dip pri may include x-component px or y-component py of the approaching wave. By virtue of reciprocity theorem one may also write Equation 1 in terms of the source-side time dip if the received seismic data is in common receiver gather order.
[0028] In another implementation, the seismic packets may be described with respect to the source wavelet and expressed using the following equation:
5i (s, g, t \ gi , t. , Pri ) = wi (t - 1. - Pri · (g - g. ))w{g - g. ) Equation 2 where w, is the source wavelet for the ith seismic packet, while W denotes a Gaussian weighting function for making the ith seismic packet compact within a local zone.
[0029] At block 240, the seismic packets from block 230 may be separated into refraction packets and reflection packets. For instance, refraction packets may be separated out from the seismic packets and sent to block 250. Reflection packets may be separated out of the seismic packets and sent to block 260.
[0030] Both refraction packets and reflections packets may include the same type of information as the seismic packets from block 230. However, different types of packets may refer to different types of pressure waves that pass through the earth's subsurface. The refraction packets may describe information regarding a refractive wave from a source to a receiver, where the path of the wave is not changed by the wave's reflection at an impedance interface between the source and the receiver. As such, refractive waves may include early-arrivals (e.g., diving waves, head waves or wide-angle reflections) that may undergo continuous refractions throughout the subsurface. Refracted waves may be strong and coherent while receiving little influence by multiples, shear or converted waves, which may make them very suitable for velocity
analysis. In another respect, refracted waves may be associated with long wave paths and large apertures, which may be used to obtain an initial model of long wavelength.
[0031] On the other hand, the reflection packets may describe information regarding a reflected wave that travels from a source to a receiver, and where the reflecting wave's path is changed by one or more reflections at impedance interfaces between the source and the receiver. During a reflection, the reflected wave may exit an interface at the same angle as the angle of incidence, which may alter the reflected wave's path through the region of interest.
[0032] Several methods may be used to separate the seismic packets produced at block 230 into refraction packets and reflection packets. In one implementation, refraction packets or reflection packets may be separated from the seismic packets using source to receiver traveltimes. For instance, a time window may be used to isolate the refraction packets from the rest of the seismic packets. As such, the time window may filter out seismic packets that arrive inside a predetermined time. The predetermined time may correspond to a designated time period when the majority of early-arrival waves have reached a corresponding receiver. The seismic packets that include acquisition times outside the designated time period may be determined to be a refraction packet, while those seismic packets with acquisition times inside the designated time period may be determined to be reflection packets or other types of data.
[0033] In another implementation, source-receiver offsets may be used in the separation process at block 240 in a similar manner to the time window described above. For instance, a distance window may be used where seismic packets that include source-receiver offsets outside the distance window may be considered refraction packets and those inside the distance window may be considered reflection packets.
[0034] At block 250, the refraction packets may be used to perform refraction traveltime tomography. For instance, using a refraction packet for a source located at s and receiver located at g, a total observed traveltime tobs, and a receiver-side time dip Pn (i.e., for packet i, a time dip in the x-direction and/or y-direction), the modeled traveltime tcal may be estimated by modeling or backprojecting a ray from the receiver at g using the receiver-side time dip pri and searching for one point X along a corresponding ray path back to the source at s.
[0035] Figure 3 illustrates an example of ray tracing that may be used in block 250. In an early-arrival traveltime tomographic algorithm, a ray path 350 between receivers 310 and point 330 may be modeled using the received model from block 220. Based on the time dip 335 recorded at receivers 310, the ray path 350 may be modeled backwards to the source 320 using the observed traveltime between the firing of the source 320 and the arrival of the seismic wave at the receivers 310. For instance, modeled ray path 360 using a calculated traveltime from the point 330 to the source 320 misses the source 320. Ray tracing may then be repeated until a ray path 340 is produced that reaches the source 320 from point 330. In one implementation, the early-arrival traveltime tomographic algorithm may take into account early-arrivals in addition to the first break arrivals, which may allow for no explicit selecting of first break arrivals. This may be advantageous as compared to known diving wave tomography.
[0036] Returning to block 250, the refraction packets may be used to model a ray path from a source to a receiver. Searching for the one point X may be performed so that a seismic ray may take a minimum amount of traveltime from the source s, through an image point X (e.g., the point 330 in Figure 3), and to the receiver g according to Fermat's principle expressed by the following equation being stationary: tad = t8t + tsx Equation 3
[0037] where tcal is the total modeled traveltime, tsx is the traveltime from the source to image point X, and tgx is the traveltime from the receiver to the image point X. tsx
may be computed using a traveltime table and tgx may be computed by ray tracing one ray.
[0038] At block 254, a portion of the received model is updated based on the refraction traveltime tomography from block 250. Similar to other methods of tomography, updating the portion of the received model may include solving the refraction traveltime tomography problem by iteratively minimizing a cost function minimization criterion. This update may involve internal nonlinear iterations. As such, updating the received model may be performed iteratively by calculating new values for the received model and using those new values at block 250 if the updated model has not converged to a cost function minimization criterion at block 258. Thus, the gradient field for a cost function regarding the received model's parameters may be calculated and a conjugate gradient method may be used to update the received model.
[0039] A cost function minimization criterion may designate when the error in an updated model approaches a specific degree of accuracy. When the updated model converges to a cost function minimization criterion, the updated model may be determined as accurate for the region of interest. If two refraction packets are for the same event, the packets' attributes may be highly correlated. In one implementation, a cost function minimization criterion may be the difference between a forward-modeled prediction of the region of interest and the received seismic data or the data in the refractions packets from block 240. One instance of a cost function may be shown by the following equation:
Equation 4
where ; and ; denote the observed and calculated traveltimes for the i wave arrivals for source s and receiver g, respectively. In updating a model using refraction traveltime tomography, the model's values with respect to one or more modeled ray paths may be adjusted in order to reduce Εγ (ιη) . Furthermore, the first order
perturbation of the cost function shown in Equation 4 may be provided by the following equation:
«1W= -∑∑∑ °bS - ^fc Equation 5
[0040] where 5tc" denotes the perturbation of the total traveltime for the calculated early-arrival ray path for the refraction packet. The model parameter perturbation dm may be determined using the following equation: dfal
< =∑ ¾/ Equation 6 dfal
where—^-denotes the traveltime changes corresponding to the change of the model
('l)l :
parameter m at subsurface grid node j along the calculated early-arrival ray path.
[0041] At block 258, it is determined whether the updated model has converged to a corresponding cost function minimization criterion. If the updated model from block 254 has not converged to a cost function minimization criterion, the process may return to block 250 to repeat refraction traveltime tomography using the updated model, an internal nonlinear iteration. If the updated model has converged to a corresponding cost function minimization criterion, the process may proceed to block 280 to generate a final model for the region of interest. In one implementation, the updated model may provide a shallow depth model for the region of interest. The process may also proceed to block 260 to further update the updated model using the reflection packets and reflection tomography.
[0042] Blocks 260-278 may describe several different process flows in regard to the reflection packets from block 240. For instance, in one implementation, the process flow from blocks 250-258 may be performed simultaneously with the process flow involving reflection tomography in blocks 260-278 in order to generate a final model at
block 280. In another implementation, the updated model produced by blocks 250-258 may be used for performing the reflection tomography at blocks 260-278. In another implementation, one or more blocks between 260-278 may be excluded from the process.
[0043] At block 260, the reflection packets from block 240 may be migrated into the depth domain (i.e., "migrated packets"). For instance, the received seismic data at block 210 may include data in the time domain that corresponds to the acquisition time of the seismic data during a survey. Through ray tracing, seismic packets may be migrated into the depth-domain by shooting a ray from receiver position gi with the receiver-side time dip pri , and searching for the source-side time dip such that the image point Xi = (xi , yi , zi ) satisfies a total traveltime image condition th which may be similar to Equation 3. The total traveltime image condition may be expressed by the following equation:
¼ + ^ , = ^ ' Equation 7 where tsjt is the traveltime from the source to an image point X tg is the traveltime from the image point to the receivers. Similar to seismic packets in the time-domain, a migrated packet may have compact support near the image point X. . Further, a seismic packet in the depth domain Ai (x \ Xi , i,0i^i,si , gi ) may be determined using the following equation:
Ai{x \ xi,ni,ei,(pi,si,g)-- W; Equation 8
where rs i g denotes the total traveltime from source s, though the image point xi and then back to the receiver gh θί is the reflection angle, φί is the subsurface azimuth, vi
represents the velocity at the packet center . along the normal direction n and Ψ^Χ, Χ^ η^) denotes a 3D tapering function.
[0044] Once a seismic packet is migrated into the depth domain, subsurface information may be obtained regarding the migrated energy such as dips, azimuths and reflection angles, as well as associated surface information about the unmigrated data like source and receiver locations, source-side or receiver-side local time dips (e.g., dt/dsx, dt/dsy, dt/dgx, dt/dgy, dt/dhx and dt/dhy), coherencies and wavelets.
[0045] In one implementation, at block 260, the migrated packets are analyzed for whether a flat common image gather is produced. If a resulting common image gather from the migrated packets is flat, the migrated packets may proceed to block 280 for generating a final model. If the common image gather is not flat, the process may proceed to block 264, an internal nonlinear iteration.
[0046] At block 264, the migrated packets are sorted and binned into different groups based on common attributes. For instance, the migrated packets may be sorted according to offsets, reflection angles, or other packet attributes. Then, the sorted migrated packets may be binned into different groups such that binned packets are for the same event (e.g., for the same fault or the same reflection at a reflection interface). For instance, migrated packets may be binned into a group based on the following criteria: ^ - ^ ^K - ) Equation 9
an where i¾ and th denote the traveltimes for two migrated packets for the same event with offsets ^ and h2 , respectively. The source wavelets for migrated packets in the same group may be highly correlated, which may serve as a grouping criteria to avoid cycle skipping. Migrated packets within a group may also be selected such that the packets within a group are located within a predetermined distance from the
corresponding event. User input may also be used to select or mask certain types of packets from the sorting or binning process, such as those packets associated with faults, steep events or a salt over-hang.
[0047] At block 268, a common image packet gather may be generated from the migrated packets. In some implementations, a common image packet gather may be formed using migrated packets that contribute to the same image location (i.e., the same location in the subsurface of the region of interest). If the migrated packets are resorted according to offsets or reflection angles, and the binned packets are reconstructed from the result, then a common image offset gather or a common image angle gathers may be obtained.
[0048] In one implementation, two migrated packets that are in a common image packet gather may be represented as Δ; and Δ; . The two migrated packets may be compared in the form of a differential semblance as described by the following equation:
Sij = -∑∑∑^ \ ^i ^i ^i , <Pi , si,gi ) j (x \ Xj , nj,ej ^j , sj ,gj ) Equation 10
y where the summation is over a local 3D space that may be in close proximity to the two migrated packets. Equation 10 may define the correlation between the two migrated packets. When a model accurately describes the region of interest, two migrated packets for the same event may be highly correlated. As such, the differential semblance for the two migrated packets may be expressed by the following equation:
Equation 1 1 where Wi and W. are 3D tapering functions for the two migrated packets.
[0049] At block 270, the migrated packets or the reflection packets from block 250 may be used to perform reflection tomography. For instance, ray tracing may be done in the depth-domain similar to the ray path modeling done in the time-domain as described in regard to Figure 3 above.
[0050] At block 274, a portion of the received model from block 220, the updated model from block 258, or a model previously updated at block 274 may be updated based on the reflection tomography from block 270. Updating this portion may include solving the reflection tomography problem by minimizing a cost function minimization criterion in a similar manner to the approach used at block 254.
[0051] In one implementation, the cost functions for the refraction traveltime tomography at block 250 and the reflection tomography at block 260 may be integrated into the following equation: Equation 12
where El is the cost function for early-arrival traveltime tomography, such as the one expressed in Equation 4, λ denotes a weight coefficient for early-arrival traveltime misfits, and sm) denotes the mismatch between any two migrated packets for Ith event in kth common image gather. The indices of k and / are omitted during the following derivations since the migrated packets may be assumed to correspond for the same events in a common image packet gather. By inserting Equation 1 1 into Equation 12, the first-order perturbation of the cost function for both the reflection tomography and refraction traveltime tomography may be obtained, and may be expressed by the following equation: δΕ =∑∑∑<¾¼ + δΕι = "ΣΣΣΣΣΣ^ ^ + W )+ λδΕι Equation 13 k I k I x y z
[0052] In Equation 13, tapering function perturbations may be ignored, and wt and Wj are time-domain source wavelets for the two migrated packets described at block 268, respectively.
[0053] If perturbations on a packet center location, a normal direction and the velocity along the normal direction are ignored, the perturbation of a source wavelet for a migrated packet may be expressed by the following equation:
[0054] where w' denotes the time derivative for the source wavelet, and δτ denotes the perturbation of the total traveltime for the ray path for the packet. The total perturbation for two migrated packets may then be expressed using the following equation:
[0055] where Equation 15 may provide the Frechet derivative for a conjugated gradient solver.
[0056] At block 278, it may be determined whether the updated model from block 274 has converged to a cost function minimization criterion. If the updated model has not converged to a cost function minimization criterion, the process may return to block 270 to repeat reflection tomography with the updated model or to block 264 to repeat the sorting and binning of the migrated packets for generating the common image packet gathers. If the updated model has converged to a corresponding cost function
minimization criterion, the process may proceed to block 280 to generate a final model for the region of interest.
[0057] At block 280, a final model for the region of interest may be generated from the updated model from block 258 or block 278. For instance, the final model may be an anisotropy model for the region of interest. The final model may also be a prestack depth migration (PSDM) volume or PSDM gather for the region of interest.
[0058] At block 290, the final model may be used to facilitate hydrocarbon exploration or production in the region of interest.
[0059] In some implementations, a method for processing seismic data corresponding to a region of interest is provided. The method may receive the seismic data. The method may separate the received seismic data into refraction packets and reflection packets. The method may receive a model for the region of interest. The method may update a first portion of the received model using the refraction packets with refraction traveltime tomography. The method may use the updated model to facilitate hydrocarbon exploration or production.
[0060] In some implementations, the method may update a second portion of the received model using the reflection packets with reflection tomography. The method may also migrate the reflection packets into the depth domain. The received seismic data may also be from a common shot gather. The method may also decompose the received seismic data into a plurality of seismic packets, where the seismic packets may be separated into refraction packets and reflection packets. Updating the first portion of the received model may be performed iteratively. The received seismic data may include attributes such as a spatial location of a beam center, a spatial orientation of a beam, a source-receiver offset, a source-receiver azimuth, a reflection angle, a reflection azimuth, a wavelet identification, an amplitude of a seismic wave, a coherency for seismic traces, and a beam spread. One of the refraction packets may include the time dip of a seismic ray measured at a receiver. The method may also
model one or more ray paths from the receiver to a source or vice versa using the time dip. The method may also model the traveltime of the ray paths from the receiver to the source or vice versa. The method may separate the received seismic data based on source to receiver traveltimes or source-receiver offsets. The refraction packets may describe information regarding a pressure wave that travels from a source to a receiver, where the path of the pressure wave is not changed by a reflection of the pressure wave at an impedance interface between the source and the receiver. The refraction packets may also describe information regarding a diving wave or a head wave. The reflection packets may describe information regarding a pressure wave that travels from a source to a receiver, and where the pressure wave is reflected from an impedance interface. The received model or the updated model may be a velocity model or an anisotropic model for the region of interest. Updating the first portion of the received model may include minimizing a cost function minimization criterion. The cost function minimization criterion may be a difference between a forward-modeled prediction and the received seismic data.
[0061] In some implementations, an information processing apparatus for use in a computing system is provided, and includes means for receiving seismic data corresponding to a region of interest. The information processing apparatus may also have means for separating the received seismic data into refraction packets and reflection packets. The information processing apparatus may also have means for receiving a model for the region of interest. The information processing apparatus may also have means for updating a portion of the received model using the refraction packets with refraction traveltime tomography. The information processing apparatus may also have means for using the updated model to facilitate hydrocarbon exploration or production.
[0062] In some implementations, a computing system is provided that includes at least one processor, at least one memory, and one or more programs stored in the at least one memory, wherein the programs include instructions, which when executed by
the at least one processor cause the computing system to receive seismic data corresponding to a region of interest. The programs may further include instructions to cause the computing system to separate the received seismic data into refraction packets and reflection packets. The programs may further include instructions to cause the computing system to receive a model for the region of interest. The programs may further include instructions to cause the computing system to update a portion of the received model using the refraction packets with refraction traveltime tomography. The programs may further include instructions to cause the computing system to use the updated model to facilitate hydrocarbon exploration or production.
[0063] In some implementations, a computer readable storage medium is provided, which has stored therein one or more programs, the one or more programs including instructions, which when executed by a processor, cause the processor to receive seismic data corresponding to a region of interest. The programs may further include instructions, which cause the processor to separate the received seismic data into refraction packets and reflection packets. The programs may further include instructions, which cause the processor to receive a model for the region of interest. The programs may further include instructions, which cause the processor to update a portion of the received model using the refraction packets with refraction traveltime tomography. The programs may further include instructions, which cause the processor to use the updated model to facilitate hydrocarbon exploration or production.
[0064] In some implementations, a method for processing seismic data corresponding to a region of interest is provided. The method may receive the seismic data. The method may decompose the received seismic data into a plurality of seismic packets. The method may separate the seismic packets into refraction packets and reflection packets. The method may receive a model that describes the region of interest. The method may update a first portion of the received model using the refraction packets with refraction traveltime tomography. The method may update a second portion of the received model using the reflection packets with reflection
tomography. The method may use an updated model based on the updated first portion and second portion of the received model to facilitate hydrocarbon exploration or production.
[0065] In some implementations, an information processing apparatus for use in a computing system is provided, and includes means for receiving seismic data corresponding to a region of interest. The information processing apparatus may also have means for decomposing the received seismic data into a plurality of seismic packets. The information processing apparatus may also have means for separating the seismic packets into refraction packets and reflection packets. The information processing apparatus may also have means for receiving a model for the region of interest. The information processing apparatus may also have means for updating a first portion of the received model using the refraction packets with refraction traveltime tomography. The information processing apparatus may also have means for updating a second portion of the received model using the reflection packets with reflection tomography. The information processing apparatus may also have means for using an updated model based on the updated first portion and second portion of the received model to facilitate hydrocarbon exploration or production.
[0066] In some implementations, a computing system is provided that includes at least one processor, at least one memory, and one or more programs stored in the at least one memory, wherein the programs include instructions, which when executed by the at least one processor cause the computing system to receive seismic data corresponding to a region of interest. The programs may further include instructions to cause the computing system to decompose the received seismic data into a plurality of seismic packets. The programs may further include instructions to cause the computing system to separate the seismic packets into refraction packets and reflection packets. The programs may further include instructions to cause the computing system to receive a model that describes the region of interest. The programs may further include instructions to cause the computing system to update a first portion of the
received model using the refraction packets with refraction traveltime tomography. The programs may further include instructions to cause the computing system to update a second portion of the received model using the reflection packets with reflection tomography. The programs may further include instructions to cause the computing system to use an updated model based on the updated first portion and second portion of the received model to facilitate hydrocarbon exploration or production.
[0067] In some implementations, a computer readable storage medium is provided, which has stored therein one or more programs, the one or more programs including instructions, which when executed by a processor, cause the processor to receive seismic data corresponding to a region of interest. The programs may further include instructions, which cause the processor to decompose the received seismic data into a plurality of seismic packets. The programs may further include instructions, which cause the processor to separate the seismic packets into refraction packets and reflection packets. The programs may further include instructions, which cause the processor to receive a model that describes the region of interest. The programs may further include instructions, which cause the processor to update a first portion of the received model using the refraction packets with refraction traveltime tomography. The programs may further include instructions, which cause the processor to update a second portion of the received model using the reflection packets with reflection tomography. The programs may further include instructions, which cause the processor to use an updated model based on the updated first portion and second portion of the received model to facilitate hydrocarbon exploration or production.
[0068] In some implementations, a method for processing data corresponding to a region of interest is provided. The method may receive data corresponding to a region of interest. The method may separate the data into refraction packets and reflection packets. The method may receive a model that describes the region of interest. The method may update a first portion of the received model using the refraction packets
with refraction traveltime tomography. The method may update a second portion of the received model using the reflection packets with reflection tomography.
[0069] In some implementations, an information processing apparatus for use in a computing system is provided, and includes means for receiving data corresponding to a region of interest. The information processing apparatus may also have means for separating the data into refraction packets and reflection packets. The information processing apparatus may also have means for receiving a model for the region of interest. The information processing apparatus may also have means for updating a first portion of the received model using the refraction packets with refraction traveltime tomography. The information processing apparatus may also have means for updating a second portion of the received model using the reflection packets with reflection tomography.
[0070] In some implementations, a computing system is provided that includes at least one processor, at least one memory, and one or more programs stored in the at least one memory, wherein the programs include instructions, which when executed by the at least one processor cause the computing system to receive data corresponding to a region of interest. The programs may further include instructions to cause the computing system to separate the data into refraction packets and reflection packets. The programs may further include instructions to cause the computing system to receive a model that describes the region of interest. The programs may further include instructions to cause the computing system to update a first portion of the received model using the refraction packets with refraction traveltime tomography. The programs may further include instructions to cause the computing system to update a second portion of the received model using the reflection packets with reflection tomography.
[0071] In some implementations, a computer readable storage medium is provided, which has stored therein one or more programs, the one or more programs including instructions, which when executed by a processor, cause the processor to receive data
corresponding to a region of interest. The programs may further include instructions, which cause the processor to separate the data into refraction packets and reflection packets. The programs may further include instructions, which cause the processor to receive a model that describes the region of interest. The programs may further include instructions, which cause the processor to update a first portion of the received model using the refraction packets with refraction traveltime tomography. The programs may further include instructions, which cause the processor to update a second portion of the received model using the reflection packets with reflection tomography.
Computing System
[0072] Implementations of various technologies described herein may be operational with numerous general purpose or special purpose computing system environments or configurations. Examples of well known computing systems, environments, and/or configurations that may be suitable for use with the various technologies described herein include, but are not limited to, personal computers, server computers, hand-held or laptop devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, smartphones, smartwatches, personal wearable computing systems networked with other computing systems, tablet computers, and distributed computing environments that include any of the above systems or devices, and the like.
[0073] The various technologies described herein may be implemented in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that performs particular tasks or implement particular abstract data types. While program modules may execute on a single computing system, it should be appreciated that, in some implementations, program modules may be implemented on separate computing systems or devices adapted to communicate with one another. A program module may also be some combination of
hardware and software where particular tasks performed by the program module may be done either through hardware, software, or both.
[0074] The various technologies described herein may also be implemented in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network, e.g., by hardwired links, wireless links, or combinations thereof. The distributed computing environments may span multiple continents and multiple vessels, ships or boats. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
[0075] Figure 4 illustrates a schematic diagram of a computing system 400 in which the various technologies described herein may be incorporated and practiced. Although the computing system 400 may be a conventional desktop or a server computer, as described above, other computer system configurations may be used.
[0076] The computing system 400 may include a central processing unit (CPU) 430, a system memory 426, a graphics processing unit (GPU) 431 and a system bus 428 that couples various system components including the system memory 426 to the CPU 430. Although one CPU is illustrated in Figure 4, it should be understood that in some implementations the computing system 400 may include more than one CPU. The GPU 431 may be a microprocessor specifically designed to manipulate and implement computer graphics. The CPU 430 may offload work to the GPU 431 . The GPU 431 may have its own graphics memory, and/or may have access to a portion of the system memory 426. As with the CPU 430, the GPU 431 may include one or more processing units, and the processing units may include one or more cores. The system bus 428 may be any of several types of bus structures, including a memory bus or memory controller, a peripheral bus, and a local bus using any of a variety of bus architectures. By way of example, and not limitation, such architectures include Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MCA) bus, Enhanced ISA (EISA) bus, Video Electronics Standards Association (VESA) local bus, and Peripheral
Component Interconnect (PCI) bus also known as Mezzanine bus. The system memory 426 may include a read-only memory (ROM) 412 and a random access memory (RAM) 416. A basic input/output system (BIOS) 414, containing the basic routines that help transfer information between elements within the computing system 400, such as during start-up, may be stored in the ROM 412.
[0077] The computing system 400 may further include a hard disk drive 450 for reading from and writing to a hard disk, a magnetic disk drive 452 for reading from and writing to a removable magnetic disk 456, and an optical disk drive 454 for reading from and writing to a removable optical disk 458, such as a CD ROM or other optical media. The hard disk drive 450, the magnetic disk drive 452, and the optical disk drive 454 may be connected to the system bus 428 by a hard disk drive interface 436, a magnetic disk drive interface 438, and an optical drive interface 440, respectively. The drives and their associated computer-readable media may provide nonvolatile storage of computer-readable instructions, data structures, program modules and other data for the computing system 400.
[0078] Although the computing system 400 is described herein as having a hard disk, a removable magnetic disk 456 and a removable optical disk 458, it should be appreciated by those skilled in the art that the computing system 400 may also include other types of computer-readable media that may be accessed by a computer. For example, such computer-readable media may include computer storage media and communication media. Computer storage media may include volatile and non-volatile, and removable and non-removable media implemented in any method or technology for storage of information, such as computer-readable instructions, data structures, program modules or other data. Computer storage media may further include RAM, ROM, erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other solid state memory technology, CD-ROM, digital versatile disks (DVD), or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage
devices, or any other medium which can be used to store the desired information and which can be accessed by the computing system 400. Communication media may embody computer readable instructions, data structures, program modules or other data in a modulated data signal, such as a carrier wave or other transport mechanism and may include any information delivery media. The term "modulated data signal" may mean a signal that has one or more of its characteristics set or changed in such a manner as to encode information in the signal. By way of example, and not limitation, communication media may include wired media such as a wired network or direct-wired connection, and wireless media such as acoustic, RF, infrared and other wireless media. The computing system 400 may also include a host adapter 433 that connects to a storage device 435 via a small computer system interface (SCSI) bus, a Fiber Channel bus, an eSATA bus, or using any other applicable computer bus interface. Combinations of any of the above may also be included within the scope of computer readable media.
[0079] A number of program modules may be stored on the hard disk 450, magnetic disk 456, optical disk 458, ROM 412 or RAM 416, including an operating system 418, one or more application programs 420, program data 424, and a database system 448. The application programs 420 may include various mobile applications ("apps") and other applications configured to perform various methods and techniques described herein. The operating system 418 may be any suitable operating system that may control the operation of a networked personal or server computer, such as Windows® XP, Mac OS® X, Unix-variants (e.g., Linux® and BSD®), and the like.
[0080] A user may enter commands and information into the computing system 400 through input devices such as a keyboard 462 and pointing device 460. Other input devices may include a microphone, joystick, game pad, satellite dish, scanner, or the like. These and other input devices may be connected to the CPU 430 through a serial port interface 442 coupled to system bus 428, but may be connected by other interfaces, such as a parallel port, game port or a universal serial bus (USB). A monitor
434 or other type of display device may also be connected to system bus 428 via an interface, such as a video adapter 432. In addition to the monitor 434, the computing system 400 may further include other peripheral output devices such as speakers and printers.
[0081] Further, the computing system 400 may operate in a networked environment using logical connections to one or more remote computers 474. The logical connections may be any connection that is commonplace in offices, enterprise-wide computer networks, intranets, and the Internet, such as local area network (LAN) 476 and a wide area network (WAN) 466. The remote computers 474 may be another a computer, a server computer, a router, a network PC, a peer device or other common network node, and may include many of the elements describes above relative to the computing system 400. The remote computers 474 may also each include application programs 470 similar to that of the computer action function.
[0082] When using a LAN networking environment, the computing system 400 may be connected to the local network 476 through a network interface or adapter 444. When used in a WAN networking environment, the computing system 400 may include a router 464, wireless router or other means for establishing communication over a wide area network 466, such as the Internet. The router 464, which may be internal or external, may be connected to the system bus 428 via the serial port interface 442. In a networked environment, program modules depicted relative to the computing system 400, or portions thereof, may be stored in a remote memory storage device 435. It will be appreciated that the network connections shown are merely examples and other means of establishing a communications link between the computers may be used.
[0083] The network interface 444 may also utilize remote access technologies (e.g., Remote Access Service (RAS), Virtual Private Networking (VPN), Secure Socket Layer (SSL), Layer 2 Tunneling (L2T), or any other suitable protocol). These remote access technologies may be implemented in connection with the remote computers 474.
[0084] It should be understood that the various technologies described herein may be implemented in connection with hardware, software or a combination of both. Thus, various technologies, or certain aspects or portions thereof, may take the form of program code (i.e., instructions) embodied in tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any other machine-readable storage medium wherein, when the program code is loaded into and executed by a machine, such as a computer, the machine becomes an apparatus for practicing the various technologies. In the case of program code execution on programmable computers, the computing device may include a processor, a storage medium readable by the processor (including volatile and non-volatile memory and/or storage elements), at least one input device, and at least one output device. One or more programs that may implement or utilize the various technologies described herein may use an application programming interface (API), reusable controls, and the like. Such programs may be implemented in a high level procedural or object oriented programming language to communicate with a computer system. However, the program(s) may be implemented in assembly or machine language, if desired. In any case, the language may be a compiled or interpreted language, and combined with hardware implementations. Also, the program code may execute entirely on a user's computing device, partly on the user's computing device, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or a server computer.
[0085] Those with skill in the art will appreciate that any of the listed architectures, features or standards discussed above with respect to the example computing system 400 may be omitted for use with a computing system used in accordance with the various embodiments disclosed herein because technology and standards continue to evolve over time.
[0086] Of course, many processing techniques for collected data, including one or more of the techniques and methods disclosed herein, may also be used successfully with collected data types other than seismic data. While certain implementations have
been disclosed in the context of seismic data collection and processing, those with skill in the art will recognize that one or more of the methods, techniques, and computing systems disclosed herein can be applied in many fields and contexts where data involving structures arrayed in a three-dimensional space and/or subsurface region of interest may be collected and processed, e.g., medical imaging techniques such as tomography, ultrasound, MRI and the like for human tissue; radar, sonar, and LIDAR imaging techniques; and other appropriate three-dimensional imaging problems.
[0087] Although the subject matter has been described in language specific to structural features and/or methodological acts, it is to be understood that the subject matter defined in the appended claims is not limited to the specific features or acts described above. Rather, the specific features and acts described above are disclosed as example forms of implementing the claims.
[0088] While the foregoing is directed to implementations of various technologies described herein, other and further implementations may be devised without departing from the basic scope thereof, which may be determined by the claims that follow. Although the subject matter has been described in language specific to structural features and/or methodological acts, it is to be understood that the subject matter defined in the appended claims is not limited to the specific features or acts described above. Rather, the specific features and acts described above are disclosed as example forms of implementing the claims.
Claims
1 . A method for processing seismic data corresponding to a region of interest, comprising:
receiving the seismic data;
separating the received seismic data into refraction packets and reflection packets;
receiving a model for the region of interest;
updating a first portion of the received model using the refraction packets with refraction traveltime tomography; and
using the updated model to facilitate hydrocarbon exploration or production.
2. The method of claim 1 , further comprising updating a second portion of the received model using the reflection packets with reflection tomography.
3. The method of claim 2, wherein updating the second portion of the received model comprises migrating the reflection packets into the depth domain.
4. The method of claim 1 , wherein the received seismic data is from a common shot gather.
5. The method of claim 1 , further comprising decomposing the received seismic data into a plurality of seismic packets, and wherein separating the received seismic data comprises separating the seismic packets into the refraction packets and the reflection packets.
6. The method of claim 1 , wherein updating the first portion of the received model is performed iteratively.
7. The method of claim 1 , wherein the received seismic data includes one or more attributes selected from the group consisting of:
a spatial location of a beam center;
a spatial orientation of a beam;
a source-receiver offset;
a source-receiver azimuth;
a reflection angle;
a reflection azimuth;
a wavelet identification;
an amplitude of a seismic wave;
a coherency for seismic traces; and
a beam spread.
8. The method of claim 1 , wherein at least one of the refraction packets comprises the time dip of a seismic ray measured at a receiver, and wherein updating the first portion of the received model comprises modeling one or more ray paths from the receiver to a source or vice versa using the time dip.
9. The method of claim 8, wherein modeling the ray paths comprises modeling the traveltime of the ray paths from the receiver to the source or vice versa.
10. The method of claim 1 , wherein separating the received seismic data comprises separating the received seismic data based on source to receiver traveltimes or source- receiver offsets.
1 1 . The method of claim 1 , wherein the refraction packets describe information regarding a pressure wave that travels from a source to a receiver, and wherein the path of the pressure wave is not changed by a reflection of the wave at an impedance interface between the source and the receiver.
12. The method of claim 1 , wherein the refraction packets describe information regarding a diving wave or a head wave.
13. The method of claim 1 , wherein the reflection packets describe information regarding a pressure wave that travels from a source to a receiver, and wherein the pressure wave is reflected from at least one impedance interface.
14. The method of claim 1 , wherein the received model or the updated model is a velocity model or an anisotropic model for the region of interest.
15. The method of claim 1 , wherein updating the first portion of the received model comprises minimizing a cost function minimization criterion.
16. The method of claim 15, wherein the cost function minimization criterion comprises a difference between a forward-modeled prediction and the received seismic data.
17. A method for processing seismic data corresponding to a region of interest, comprising:
receiving the seismic data;
decomposing the received seismic data into a plurality of seismic packets;
separating the seismic packets into refraction packets and reflection packets; receiving a model that describes the region of interest;
updating a first portion of the received model using the refraction packets with refraction traveltime tomography;
updating a second portion of the received model using the reflection packets with reflection tomography; and
using an updated model based on the updated first portion and second portion of the received model to facilitate hydrocarbon exploration or production.
18. The method of claim 17, wherein separating the seismic packets comprises separating the seismic packets based on source to receiver traveltimes or source- receiver offsets.
19. The method of claim 17, wherein at least one of the refraction packets comprises the time dip of a seismic ray measured at a receiver, and wherein updating the portion of the received model comprises modeling one or more ray paths from the receiver to a source or vice versa using the time dip.
20. A method, comprising:
receiving data corresponding to a region of interest;
separating the data into refraction packets and reflection packets;
receiving a model that describes the region of interest;
updating a first portion of the received model using the refraction packets with refraction traveltime tomography; and
updating a second portion of the received model using the reflection packets with reflection tomography.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP14737755.0A EP2943813A4 (en) | 2013-01-14 | 2014-01-14 | Seismic data processing |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201361752011P | 2013-01-14 | 2013-01-14 | |
US61/752,011 | 2013-01-14 | ||
US14/153,894 | 2014-01-13 | ||
US14/153,894 US20140200816A1 (en) | 2013-01-14 | 2014-01-13 | Seismic data processing |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2014110547A1 true WO2014110547A1 (en) | 2014-07-17 |
Family
ID=51165795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2014/011407 WO2014110547A1 (en) | 2013-01-14 | 2014-01-14 | Seismic data processing |
Country Status (3)
Country | Link |
---|---|
US (2) | US20140200816A1 (en) |
EP (1) | EP2943813A4 (en) |
WO (1) | WO2014110547A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105425286A (en) * | 2015-10-30 | 2016-03-23 | 中国石油天然气集团公司 | Earthquake time-travelling acquisition method and crosshole earthquake time-travelling tomography method based on the earthquake time-travelling acquisition method |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2534777B (en) * | 2013-11-19 | 2020-07-29 | Shell Int Research | User suite for interactive optimal stacking of wide azimuth seismic data |
CA2947847C (en) * | 2014-05-09 | 2018-08-14 | Exxonmobil Upstream Research Company | Efficient line search methods for multi-parameter full wavefield inversion |
US10422898B2 (en) | 2014-09-23 | 2019-09-24 | Westerngeco L.L.C. | Seismic data processing |
US10379236B2 (en) * | 2015-07-22 | 2019-08-13 | Conocophillips Company | WAVSEIS sourcing |
WO2017051242A1 (en) * | 2015-09-27 | 2017-03-30 | Schlumberger Technology Corporation | Work flow based acoustic processing system and method |
US10656294B2 (en) | 2017-05-17 | 2020-05-19 | Saudi Arabian Oil Company | Generating a velocity model using subsurface azimuth and reflection angle dependent full waveform inversion |
US10948617B2 (en) * | 2017-12-11 | 2021-03-16 | Saudi Arabian Oil Company | Generating a velocity model for a subsurface structure using refraction travel time tomography |
US11320557B2 (en) | 2020-03-30 | 2022-05-03 | Saudi Arabian Oil Company | Post-stack time domain image with broadened spectrum |
CN111856580B (en) * | 2020-07-27 | 2022-02-15 | 广州海洋地质调查局 | Ultra-long offset OBS data first-arrival wave energy enhancement method and processing terminal |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070189117A1 (en) * | 2003-08-23 | 2007-08-16 | Westerngeco L.Lc. | Multiple attenuation method |
US20070258324A1 (en) * | 2006-05-03 | 2007-11-08 | Baker Hughes Incorporated | Sub-Salt Reflection Tomography and Imaging by Walkaway VSP Survey |
US20070271041A1 (en) | 2006-05-17 | 2007-11-22 | Chengbin Peng | Diplet-Based Seismic Processing |
US20100312480A1 (en) * | 2009-04-24 | 2010-12-09 | Hansteen Fredrik | Method for monitoring fluid flow in a multi-layered system |
US20110134722A1 (en) * | 2009-12-07 | 2011-06-09 | Mediwound Ltd. | Simultaneous Joint Inversion of Surface Wave and Refraction Data |
US20120113749A1 (en) * | 2010-11-04 | 2012-05-10 | Pascal Edme | Computing a Calibration Term Based on Combining Divergence Data and Seismic Data |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3611278A (en) * | 1969-07-17 | 1971-10-05 | Mobil Oil Corp | Interval velocity determination |
US4110728A (en) * | 1976-09-24 | 1978-08-29 | Texaco Inc. | Seismic data processing apparatus for obtaining alpha values from seismic data |
ATE533073T1 (en) * | 2001-09-07 | 2011-11-15 | Inova Ltd | DEVICE AND METHOD FOR COLLECTION OF SEISMIC DATA |
AU2003265168A1 (en) * | 2003-09-16 | 2005-04-06 | Geosystem S.R.L. | Wide-offset-range pre-stack depth migration method for seismic exploration |
EP2667224B1 (en) * | 2006-09-04 | 2019-08-07 | Geosystem S.r.l. | Method for building velocity models for pre-stack depth migration via the simultaneous joint inversion of seismic, gravity and magnetotelluric data. |
US7970546B1 (en) * | 2007-06-22 | 2011-06-28 | Westerngeco L.L.C. | Diplet-based imaging of seismic data in shot or receiver records |
-
2014
- 2014-01-13 US US14/153,894 patent/US20140200816A1/en not_active Abandoned
- 2014-01-14 EP EP14737755.0A patent/EP2943813A4/en not_active Withdrawn
- 2014-01-14 WO PCT/US2014/011407 patent/WO2014110547A1/en active Application Filing
-
2015
- 2015-07-22 US US14/805,962 patent/US20150323688A1/en not_active Abandoned
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070189117A1 (en) * | 2003-08-23 | 2007-08-16 | Westerngeco L.Lc. | Multiple attenuation method |
US20070258324A1 (en) * | 2006-05-03 | 2007-11-08 | Baker Hughes Incorporated | Sub-Salt Reflection Tomography and Imaging by Walkaway VSP Survey |
US20070271041A1 (en) | 2006-05-17 | 2007-11-22 | Chengbin Peng | Diplet-Based Seismic Processing |
US20100312480A1 (en) * | 2009-04-24 | 2010-12-09 | Hansteen Fredrik | Method for monitoring fluid flow in a multi-layered system |
US20110134722A1 (en) * | 2009-12-07 | 2011-06-09 | Mediwound Ltd. | Simultaneous Joint Inversion of Surface Wave and Refraction Data |
US20120113749A1 (en) * | 2010-11-04 | 2012-05-10 | Pascal Edme | Computing a Calibration Term Based on Combining Divergence Data and Seismic Data |
Non-Patent Citations (2)
Title |
---|
M. C. TANIS ET AL.: "Diving-wave refraction tomography and reflection tomography for velocity model building", SEG TECHNICAL PROGRAM EXPANDED ABSTRACTS 2006, 1 January 2006 (2006-01-01), pages 3340 - 3344, XP055251816, DOI: doi:10.1190/1.2370225 |
See also references of EP2943813A4 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105425286A (en) * | 2015-10-30 | 2016-03-23 | 中国石油天然气集团公司 | Earthquake time-travelling acquisition method and crosshole earthquake time-travelling tomography method based on the earthquake time-travelling acquisition method |
Also Published As
Publication number | Publication date |
---|---|
EP2943813A1 (en) | 2015-11-18 |
US20140200816A1 (en) | 2014-07-17 |
EP2943813A4 (en) | 2016-04-06 |
US20150323688A1 (en) | 2015-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3607359B1 (en) | Generating common image gather using wave-field separation | |
US20150323688A1 (en) | Seismic data processing | |
US10317547B2 (en) | Noise model estimation in multimeasurement data | |
CA2801531C (en) | Method and system for determining source signatures after source ghost removal | |
US10520623B2 (en) | Methods and systems for marine survey acquisition | |
US9869783B2 (en) | Structure tensor constrained tomographic velocity analysis | |
EP3359982B1 (en) | Seismic sensor orientation | |
US11536863B2 (en) | Noise attenuation of multiple source seismic data | |
EP2981848B1 (en) | Joint inversion of geophysical attributes | |
US11092708B2 (en) | Processes and systems to enhance illumination and resolution of seismic images using multiple reflected wavefields | |
WO2015164388A1 (en) | Seismic data processing | |
Harding et al. | Velocity structure near IODP Hole U1309D, Atlantis Massif, from waveform inversion of streamer data and borehole measurements | |
Thiel et al. | Comparison of acoustic and elastic full‐waveform inversion of 2D towed‐streamer data in the presence of salt | |
US10132946B2 (en) | Methods and systems that combine wavefields associated with generalized source activation times and near-continuously recorded seismic data | |
EP3400464B1 (en) | Seismic data enhancement | |
EP3198308B1 (en) | Seismic data processing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 14737755 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2014737755 Country of ref document: EP |