EP3351972A1 - Inversion iterative des donnees sismiques codees a base de la construction des donnees pseudo supersource - Google Patents

Inversion iterative des donnees sismiques codees a base de la construction des donnees pseudo supersource Download PDF

Info

Publication number
EP3351972A1
EP3351972A1 EP18154578.1A EP18154578A EP3351972A1 EP 3351972 A1 EP3351972 A1 EP 3351972A1 EP 18154578 A EP18154578 A EP 18154578A EP 3351972 A1 EP3351972 A1 EP 3351972A1
Authority
EP
European Patent Office
Prior art keywords
source
records
record
sources
simultaneous
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
EP18154578.1A
Other languages
German (de)
English (en)
Inventor
Christine E. Krohn
Partha S. Routh
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
ExxonMobil Upstream Research Co
Original Assignee
ExxonMobil Upstream Research Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ExxonMobil Upstream Research Co filed Critical ExxonMobil Upstream Research Co
Publication of EP3351972A1 publication Critical patent/EP3351972A1/fr
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/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/003—Seismic data acquisition in general, e.g. survey design
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/003—Seismic data acquisition in general, e.g. survey design
    • G01V1/005—Seismic data acquisition in general, e.g. survey design with exploration systems emitting special signals, e.g. frequency swept signals, pulse sequences or slip sweep arrangements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/24—Recording seismic data
    • G01V1/247—Digital recording of seismic data, e.g. in acquisition units or nodes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • 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
    • 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/364—Seismic filtering
    • G01V1/368—Inverse filtering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • G01V1/3808—Seismic data acquisition, e.g. survey design
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/10—Aspects of acoustic signal generation or detection
    • G01V2210/12—Signal generation
    • G01V2210/121—Active source
    • G01V2210/1212—Shot
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/10—Aspects of acoustic signal generation or detection
    • G01V2210/12—Signal generation
    • G01V2210/127—Cooperating multiple sources
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/10—Aspects of acoustic signal generation or detection
    • G01V2210/14—Signal detection
    • G01V2210/142—Receiver location
    • G01V2210/1423—Sea
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/30—Noise handling
    • G01V2210/32—Noise reduction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616—Data from specific type of measurement
    • G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
    • 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/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/66—Subsurface modeling

Definitions

  • This disclosure relates generally to the field of geophysical prospecting and, more particularly, seismic data processing as used in hydrocarbon exploration. Specifically, the disclosure relates to a method for acquiring, at reduced seismic acquisition cost, data using simultaneous sources in the field, and then constructing pseudo source-records that better meet the requirements for using additional simultaneous sourcing for computer simulations or forward modeling as part of iterative inversion methods that update a subsurface model, such as FWI (Full Wavefield Inversion) or LSRTM (Least-Square Reverse Time Migration), with additional reduction in computational costs.
  • FWI Full Wavefield Inversion
  • LSRTM Least-Square Reverse Time Migration
  • Simultaneous sourcing is an emerging seismic acquisition method for reducing acquisition costs and improving spatial sampling.
  • surveys are acquired by locating a single point source or an array of point sources at a single source location, firing the sources at the same time and then recording the response for the time needed for the sources to finish firing followed by a listening time in which all returns from the subsurface target are recorded.
  • the firing of the sources can be repeated and multiple records can be recorded at the same location. Then, the source array is moved to another location, and the process is repeated.
  • the cost of acquiring seismic data by this sequential method is related to the time needed to record each individual source location and the number of such locations, and this cost often limits the ability to record data at fine sampling.
  • Simultaneous sourcing is most commonly used for vibroseis sources with long sweep functions, which can be easily encoded.
  • each individual vibrator can be driven by a sweep that differs in some manner from the sweeps for other vibrators within the array, for example using differences in the sweep phase, pseudorandom function, sweep rate, sweep frequency range, start time, etc.
  • Some methods require multiple sweeps and multiple records per location for separation. In the special case that the number of sweeps is greater than or equal to the number of vibrators, then the individual source records can be almost perfectly extracted from the multiple combined records by applying an inverse filter as described for the HFVS method in Sallas, et al. (U.S. Patent No. 5,721,710 ).
  • the cross talk noise is typically mitigated with an iterative data inversion and separation method ( Neelamani, et al., U.S. Patent No. 8,248,886 ) or by filtering ( Huo et al., U.S. Patent Publication No. 2012/0290214 ).
  • Simultaneous sourcing can also be used for impulsive sources but there are fewer and less powerful methods to encode impulsive sources. There is little cost saving benefit for use of simultaneous sourcing for land acquisition with dynamite, but use of simultaneous sourcing for airguns in marine acquisition can be beneficial, especially for wide-azimuth acquisition.
  • the use of random firing times for marine sources firing nearly simultaneously but located on different vessels was disclosed by Vaage (U.S. Patent No. 6,906,981 ). More recently, simultaneous sourcing has been proposed for multiple vessel shooting of wide-azimuth (WAZ) marine surveys ( Beasley et al., "A 3D simultaneous source field test processed using alternating projections: a new active separation method," Geophysical Prospecting 60, 591-601 (2012 )). Simultaneous sourcing is the only way that finely spaced (e.g. 25-m) source points, can be acquired in a single pass of the streamers. Without simultaneous sourcing, multiple passes are required and the survey takes much longer and costs are significantly
  • FIG. 1 We illustrate one configuration for a WAZ marine survey, in Fig. 1 to show the benefit of simultaneous sourcing.
  • the figure shows source line 123 , which is traversed by a source boat, and receiver line 121 , which is traversed by a boat pulling multiple streamers of hydrophones. Both boats move in parallel at the same speed, typically a minimum of 6 knots.
  • the first position of the boats are shown in black, and future positions are shown in grey.
  • the source boat is fired at position 103 while the receiver boat is at position 101 , and the response is recorded typically for about 10 s. During this 10 s, the boats are moving.
  • the source boat reaches the next shot point at 113 , typically 20-40 m from the previous shot, and the receiver boat reaches position 111.
  • the receiver boat can make 4 passes of receiver line 121 while a source boat traverses source lines 123-126 in sequence. This is an expensive option, but can yield a fine source sampling for each source line, for example, a 25-m source interval.
  • 4 source boats can be used, and the sources fired flip-flopping between lines. For example, a source can be fired at position 103 , then at position 114 , 135 , 136 , and then 143.
  • the source interval along each line is 200 m, much coarser than the fine-spaced survey. It is not possible to shoot and record at a finer shot spacing, because by the time the full record is acquired, the boats have moved tens of meters along the sail lines.
  • a finely spaced survey can be recorded with simultaneous shooting by firing all four sources within the same time record but with a small random delay or jitter in either the firing time or position. For example, the sources can be fired at positions 103, 104, 105, and 106 to form one record with overlapping source energy. Then for the next record, sources are fired at 113, 114, 115, and 116, etc.
  • the jitter is a form of encoding that allows the interference to be partially removed by filtering in processing. Since the boats are moving, a delay in firing time necessarily means a slight shift in the firing position around the nominal sourcing interval as determined by the speed of the vessel. Instead of requiring vessel-to-vessel time synchronization, it can be operationally simplier to implement random time delays by generating a "preplot" of sourcing positions along each line with random positional variations around the nominal source interval. During acquisition, each vessel shoots independently of the other vessels at the predermined sourcing positions. With this method, the exact firing position but not the firing time is predetermined, but the result is still randomization in time. In the current invention, the randomization of sourcing time or position is understood to be equivalent. In either case, it is important to determine the actual firing position and firing time and these values along with other sourcing characteristics comprise the encoding function.
  • Step 201 source records of length T record > T listen with multiple source excitations during the record are obtained. Some sort of field-encoding scheme such as jittered start times or position is used during the sourcing. T listen is the time needed for the energy to travel from the source to the target and then to the receiver.
  • Step 304 is a simple illustration of a source record for a single-source. In these and subsequent diagrams, the response of a single source is illustrated with a linear event 302 and a hyperbolic event 303.
  • the four sources fire with small time delays and a linear and hyperbolic event from each of the four sources interfere.
  • these 4 sources are at long crossline distances on the source lines 123-126 in Fig. 1 , but only one boat and source are shown in the cross-section view of 309.
  • the encoding functions including the source positions and start times, are determined for the sources that contribute to the records.
  • the source location and time projected onto the record window is indicated by the sunburst 301 for the single source.
  • the corresponding source positional variation is relatively small compared to the scale of the figure and is not illustrated in the diagram.
  • the simultaneous source record 309 is generated by sources at projected positions 305 , 306 , 307 , and 308 , each having a small time shift relative to each other. Then in Step 203 , the encoding function is used to extract individual source records, one for each source starting at the firing time of that source and continuing for the appropriate listening time T listen . For vibroseis data, this extraction can include the process of correlation by the particular sweep used for that source.
  • the single record 309 is copied 4 times and then shifted in time so that the record starts (zero time) at the firing time for each respective source.
  • the record 315 is made by copying record 309 and time shifting to the time of source 305.
  • Step 204 further processing methods are used to filter the interfering energy that is not desired on each source gather, or to use a sparse inversion scheme to improve the separation of the data, resulting in a separated seismogram for each source as if it has been recorded independently of the other sources.
  • Step 205 the separated source gathers can be conventionally imaged or inverted.
  • Fig. 2 shows a single long continuous record 401 with multiple source excitations illustrated by sunbursts. Unlike the marine streamer case, the sources are not fired at small intervals compared to the record length, and thus conditions are not met for use of the term areal array.
  • the initial pseudo-separation Step 203 involves extracting windows the size of the desired record length T listen starting at the firing time of one of the sources as shown in Fig. 4B .
  • a window corresponding to 412 starting at the source firing time 402 is copied and extracted to make source record 422. It has interfering energy from other sources at 433 and 424.
  • a window 413 is copied and extracted for source 403 to make source record 423
  • a window 414 is extracted for source 404 making source record 424 , and so on.
  • Simultaneous sourcing followed by source separation can also be used to assist with computationally-expensive seismic data simulation or forward modeling as described in Neelamani et al. (U.S. Patent No. 8,248,886 ).
  • Such forward modeling is a component of seismic imaging or seismic inversion with the output being an image of reflectivity or of formation properties such as the seismic velocity of the subsurface.
  • Forward modeling uses a detailed velocity model and computes the complex wavefields theoretically generated by each source. Considerable computer time can be saved by reducing the number of sources to be modeled at one time by using simultaneous sourcing with some sort of encoding scheme, and then separating the data into the individual source seismograms.
  • This method is identical to the field acquisition, but there are more choices of encoding schemes when done in the computer, and the specific encoded-sequence for a source is perfectly known.
  • One common encoding scheme is to use random scaling in which the output of each source is randomly multiplied by either plus or minus one. This scheme cannot be physically implemented in the field for impulse sources such as airguns or explosives.
  • simultaneous sourcing can be used to lower costs to acquire seismic data in the field or to simulate seismic data in the computer.
  • This involves recording one or more composite records containing interference from multiple sources.
  • This can be a short record with sources excited close together in time and forming a spatial source array. It also can be continuous long record with individual sources excited at random or fixed intervals.
  • the composite record must be separated into individual source gathers. Typically, this involves pseudo-separation by extracting a window around the firing-time of the sources and then using filtering or inversion operations to remove interference noise or crosstalk.
  • the number of records are the same or greater than the number of individual sources within a spatial array, the separation is quite good, but acquiring multiple records is expensive. With fewer records, there is a problem in that the separation is imperfect with some crosstalk noise remaining or important signal removed by the filtering or inversion.
  • Simultaneous sourcing is also used to save computational cost associated with imaging and inversion of seismic data.
  • individual seismic source gathers that were acquired sequentially, i.e. one source or source array shot at a time, are encoded in the computer and summed to form a simultaneous source record that is then used to form an image of seismic reflectivity or to determine subsurface properties.
  • Use of this method to increase the speed and reduce cost of conventional (non-iterative and does not improve a subsurface model) migration is disclosed by Ober et al. (U.S. Patent No. 6,021,094 ) and use of the method in inversion is disclosed by Krebs, et al. (U.S. Patent No. 8,121,823 ).
  • Crosstalk or interference between sources is also a problem for this use of simultaneous sourcing and such crosstalk manifests itself as noise in the imaging and inversion outputs.
  • the crosstalk can be minimized somewhat by optimizing the computer encoding functions, such as using random scaling instead of phase rotation, but the results may not be as good as the more computer-intensive sequential use of individual sources.
  • Simultaneous sourcing is particularly useful for inversion, such as full waveform inversion (FWI) and least-square reverse-time migration (LSRTM).
  • FWI full waveform inversion
  • LSRTM least-square reverse-time migration
  • FWI full waveform inversion
  • LSRTM least-square reverse-time migration
  • the sources can be re-encoded and re-summed every iteration and then used for a model update ( Krebs, U.S. Patent No. 8,121,823 ).
  • Each group of encoded and summed data may be called a realization of the data. The best results and reduced crosstalk are achieved when multiple realizations are used in the iterative process.
  • Step 501 a number of field records are obtained, each with the same spread extent L spread and record duration T listen .
  • the record duration T listen should include the time needed for seismic waves to travel from the source to the target and then to the receivers.
  • a single source or areal array can be used for each record.
  • the records are then computer encoded, preferably with a randomized encoding scheme in Step 502. For example, the records can be randomly multiplied by plus or minus 1 or phase rotated by a random factor. Then all the records in the sail line or swath or in the entire survey are summed or stacked, forming one simultaneous source record. This is called one realization of the data.
  • Step 503 the seismic response is simulated in the computer for traveltime T listen for all the sources at one time using the computer encoding scheme.
  • This step uses an initial or updated model.
  • the simulated and measured records are compared in Step 504 , and the comparison or misfit function is used to update the subsurface image or property model. If multiple iterations (Step 506) are needed, it is preferable to go back to Step 502 and re-encode the field records, making a second realization of the data. By changing the encoding each iteration, artifacts and residual noise are reduced.
  • simultaneous-sourcing for iterative inversion assumes that the receiver spread and record length are fixed, i.e. all receivers are recording for all sources for the same length of time so that the records can be summed together.
  • the computer is used to forward-model all the sources into all of the receivers as if they were initiated at the same time or nearly the same time. If the point source data are not recorded with a fixed spread, for example if different receiver locations are used to record different shots, then the forward-modeling case does not match the field data case. This can create problems in that the misfit function, the difference between the field and forward-modeled data, will be dominated by the missing energy between the forward modeling and measured data and will not be useful for updating the trial model.
  • Figure 6 illustrates the problem with acquiring data conventionally with a marine streamer and then using simultaneous sourcing to reduce the computation effort required in inversion.
  • a source is fired at position 602 and a record 601 is captured.
  • the boat then moves forward to position 604 and captures record 603 and then to location 606 for record 605.
  • the receivers are moving so the actual receivers are at different locations along the source line. If all the traces are arranged by their true positions along the sail lines, encoded and summed, a simultaneous source gather 610 is obtained. Then, if the three sources ( 622, 623, 624, corresponding to 602, 604, and 606 ) are simultaneously excited in the computer, the simulated record 612 is obtained.
  • the present invention uses simultaneous sourcing in the field in such a way as to overcome problems from non-fixed spreads and long recording times to yield a plurality of pseudo super-source records that can be computer encoded and stacked to make multiple realizations of the data that can be changed each iteration of the inversion.
  • This invention is a method for acquiring, at reduced seismic acquisition cost, data using simultaneous sources in the field, and then constructing pseudo source-records that better meet the requirements for using additional simultaneous sourcing for computer simulations or forward modeling as part of iterative inversion, such as FWI (Full Wavefield Inversion) or LSRTM (Least-Squares Reverse Time Migration), with additional reduction in computational costs.
  • FWI Full Wavefield Inversion
  • LSRTM Least-Squares Reverse Time Migration
  • the method can be used for marine streamer acquisition and other non-fixed spread geometries to acquire both positive and negative offsets and to mitigate the "missing data" problem for simultaneous-source FWI. It can also be used for land data to overcome issues with moving spreads and long continuous records, where a long continuous record means a data record too long to be effectively computer simulated.
  • a first embodiment of the invention is a method for performing simultaneous inversion (without separation) of multiple sources where the data being inverted are field data records generated by two or more interfering or overlapping sources. Steps of this method may include:
  • a second embodiment of the invention is an application of the first embodiment to data acquired under survey conditions in which the fixed-receiver assumption necessary for simultaneous-source inversion is not satisfied. Steps of this method may include:
  • the above-described first embodiment of the invention may be used without the additional features of the second embodiment, for example when processing data where all sources illuminate a full spread of receivers.
  • the updated or adjusted velocity model resulting from the present inventive method may be used to migrate the seismic data to generate an image of the subsurface, or for other seismic data processing and interpretation purposes relating to exploration for hydrocarbons.
  • Fig. 16 is a black-and-white reproduction of a color original.
  • the invention will be described in connection with example embodiments. However, to the extent that the following detailed description is specific to a particular embodiment or a particular use of the invention, this is intended to be illustrative only, and is not to be construed as limiting the scope of the invention. On the contrary, it is intended to cover all alternatives, modifications and equivalents that may be included within the scope of the invention, as defined by the appended claims.
  • the invention is first described in its basic form, then specific embodiments for marine and land data are described.
  • This invention uses simultaneous sourcing in the field in such a way as to enhance the ability to further use simultaneous sourcing in iterative inversion by reducing the effects of crosstalk noise and better approximating acquisition by a fixed spread of receivers.
  • the invention constructs, from acquisition records, what may be called pseudo super-source (or super-shot) records, each with the same duration and spatial extent.
  • Each pseudo super-source record contains recorded energy from multiple sources, each source energized with a field encoding scheme (e.g., random time shifts, random source positions, phase rotations, sweep function, or other method) and each record is constructed by the operations of windowing, time shifting, summing and appending the original field records.
  • a field encoding scheme e.g., random time shifts, random source positions, phase rotations, sweep function, or other method
  • the survey is acquired in a manner that allows these pseudo super-shot records to be constructed so that the sources can be properly simulated simultaneously in a computer.
  • seismic energy that would be generated by a synthetic source and recorded within a predetermined distance D source within the spread and time duration T source is represented within the measured pseudo super-shot record.
  • D source within the spread and time duration
  • T source is represented within the measured pseudo super-shot record.
  • the multiple super-shot records are then separately encoded in the computer, preferably with random scaling such as multiplying by randomly selected +1 or -1, and then summed and used for inversion or imaging.
  • the computer encoding scheme is changed in subsequent iterations of the inversion of the inversion or imaging.
  • Step 701 a seismic survey is planned including choosing the type of source and optimal number and spacing of desired sourcing locations and a field encoding scheme.
  • the field encoding scheme can include any parameter related to the source excitation including its location, firing time, frequency components, source phase, etc.
  • two parameters are selected: T source and D source .
  • T source is a time greater than or equal to T listen and, in addition, is a time long enough that the amplitude of that source is reduced to a level where it no longer provides substantial interference noise.
  • T source may be selected to be equal to T listen unless a long continuous record such as Fig. 12 has been acquired.
  • D source is related in that it is equal to or larger than the maximum offset distance of interest or equal to the distance by which the source amplitude is reduced so that it no longer provides substantial interference noise.
  • These parameters are used to guide the acquisition and subsequent generation of pseudo records such that the pseudo records can be accurately simulated within a distance D source and time T source of any source location and firing time. Ideally the parameters are as small as possible while still large enough to include returns from the target. Smaller parameters mean less stringent requirements for the acquisition and construction of pseudo-gathers.
  • Step 702 one or more field records are obtained that are generated with "simultaneous" sourcing so that energy from the different sources partially overlaps in time.
  • the sources do not have to be activated exactly simultaneously, and the small time shifts between them are one way of performing the field encoding referred to in Step 702.
  • a field record is typically all -- or a subset -- of the data recorded by the active receivers (moving or stationary) in one period of time, with a start time and a stop time and no gaps.
  • the field records can be discrete records of a fixed time duration or they can be a single, continuous time record.
  • the recording spread is moved during acquisition, then preferably some of the source points within the distance D source of the boundary of the first spread are repeated into the second spread with the same encoding scheme previously used so that all energy within the distance D source is recorded on both sets of spread positions so they can be appended together.
  • Step 703 a plurality of what may be called pseudo super-source records of fixed extent and duration are constructed.
  • the record extent would span the survey width, as if the survey had been recorded by a fixed spread of receivers the width of the survey, and the record duration would be at least as long as the time for seismic energy to propagate from the source to the target and to the receivers at the maximum useable distance or offset from the source.
  • the construction process can include operations such as extractions of various time windows and trace regions from the field records.
  • a pad of zero traces can be attached and a pad of time can be added before or after the windows.
  • An objective of the construction of a pseudo super-source record is that every receiver location within an offset distance D source from the location of the source has appropriate data, i.e. data that would have been recorded if there had been a fixed receiver spread when the source shot occurred. Typically, the data from every field record will appear in at least one pseudo super-source record.
  • the various windows can then be appended or summed together to form a pseudo super-shot record.
  • each shot that influences or contributes to the region of interest is identified along with its field encoding function, and start time relative to the zero time of the pseudo super-shot record.
  • the contributing or influential shots can be assumed to be those for which the source is excited within the distance D source and a time T source from the boundary of the region of interest. This information is combined with the computer-encoding function and used for the computer simulation step 706.
  • the different pseudo super-shot records are computer encoded, preferably, but not necessarily (any incoherent encoding scheme will work), by random scaling in which they are randomly multiplied by plus or minus one ( ⁇ 1). Then all the pseudo records are summed together to form one simultaneous source record. The computer is then used to compute the forward modeling simulation in one step for all the sources within the simultaneous record, which were identified in Step 704 , as if all the sources had been fired simultaneously or nearly simultaneously (Step 706 ) .
  • Step 706 a simultaneous-source record corresponding to a simultaneous-source measured record from Step 705
  • the simultaneous-source simulated record is generated using a combination of the computer encoding that was used in step 705 combined with the field encoding from step 704/702.
  • Step 707 the recorded records from 703 and the simulated records from 706 are compared over a region of interest, and the results are used to update the subsurface model. If more iterations of the imaging or inversion is needed as determined in Step 708, then preferably the computer encoding Step 705 is repeated with a different encoding function.
  • Step 701 acquisition is planned that modifies the acquisition geometry shown in Figure 1 , by locating a source at the rear of the streamer for every source in the front of the streamer as shown in Fig. 8 . As the source and streamer boats move forward, the rear sources follow the same source line as the front sources.
  • sources 803 and 807 track source line 801
  • sources 804 and 808 track source line 802.
  • This is different from conventional wide-azimuth geometries that may use a rear source ( Treadgold, et al., "Implementing A wide Azimuth Towed Streamer Field Trial, The What, Why, and Mostly How of WATS in Southern Green Canyon”, SEG Expanded Abstracts, 2901-2903 (2006 )), but locates the rear sources on different sail lines.
  • the near offset distances are the same for both front and rear sources.
  • the distance from source 803 to the nearest receiver 805 is the same as from source 807 to its nearest receiver 806.
  • All the sources fire within the same source interval but with different random time delays or random positions around the nominal source location, and a single record of fixed length is recorded as illustrated as 901 in Fig. 9 .
  • two source boats, 803 and 804 are at the front of the streamer and two source boats 807 and 808 fire at the rear of the streamer.
  • the position of the shot projected on the source record from two front sources are shown by the sunbursts in 903.
  • the boats move forward, all four sources firing every shot interval, for example every 25-m, with random time delays or random position jitter for each group of shots. This random time delay or position jitter is considered an encoding function.
  • source record 906 is recorded with the rear sources firing with the same pattern 904 as when this position was occupied by the front sources for record 903 at 901. And again when the rear sources reach position 907 they fire with the same encoding function and at the same location as when the front sources were at 905 ; see record 908. This pattern preferably continues for the complete sail line.
  • Step 703 a pseudo super-source record is constructed.
  • Each record that was recorded with identical rear and front sources at the same position are time aligned to match the source timing and appended. Traces may be padded (i.e., zeros added) at the end or beginning.
  • This pseudo super-source record now better approximates a fixed spread because both positive and negative offsets are recorded from each source position up to a distance of D source .
  • D source is naturally the streamer length. Now all these sources can be simultaneously simulated in the computer, for example by putting groups of sources at 922 , 932 , 942 and 952.
  • Step 703 additional pseudo super-source records are constructed, each having the same spatial extent and time duration as illustrated in Fig. 10 . If the first shot was at position 0, then it starts the first super-source record. At the next interval, e.g. 25-m from the first, the second super-source record is started - its shots are shifted one shot interval from the first. This continues until the full streamer length is used and the rear sources reach the first location of the front sources. For example, super source records 1001-1004 have source positions that are shifted by a source interval. Generally, the number of super-shot records that can be generated is related to the length of the streamer divided by the shot interval. Trace padding, for example at 1006 , may be needed to construct a fixed spread size for each super-record.
  • Step 704 the source location and encoding information, including time shifts, are determined for each pseudo super-source record relative to the boundaries of the pseudo record. For example, the start time of each source is adjusted by the time shifts used to form the pseudo record and is now relative to zero time of the pseudo record. This information will be used in step 706 , combined with the computer-encoding used in Step 705 , in simultaneously simulating the encoded pseudo records.
  • each pseudo super-source record containing many shots is encoded in the computer. Preferably this is done by randomly multiplying by +1 or -1. Alternatively, phase rotations or aother form of encoding can be used. Then, the encoded pseudo records are stacked or summed, as shown in the illustration of 1112 in Fig. 11 , to form one simultaneous source record for the full sail line.
  • Step 706 all the sources in the sail line are computer-simulated at one time using a combination of the field encoding determined in Step 704/702 and the computer encoding used in Step 705. Further savings in computational cost may be achieved by limiting the region of the model used in a single-sail line simulation. This simulation is illustrated with the sunbursts in 1116. This may involve extending or padding the modeling space by an additional region as indicated in 1113 , which allows the forward modeling to generate all the bits of energy recorded in the data window 1118.
  • Step 707 the measured simultaneous source record over the region of interest ( 1118 in the example) is compared to the simulated simultaneous source record and the result is used to update the subsurface model.
  • the computer encoding is preferably changed by repeating Step 705 with a different and encoding and thereby forming multiple realizations of the data that further reduce crosstalk and artifacts.
  • Step 701 Data acquisition on land or on the ocean bottom is considered next.
  • sensors are not moving as in a marine streamer, but the group of active receivers may change during the survey. This is often called rolling the spread.
  • Fig. 12 one record at one spread position is shown as 1212 and the second record corresponding to a different set of trace positions are shown in 1213. Often the receivers may record continuously making one long record for each spread position.
  • Step 701 the survey is planned, and a T source and D source identified. Sourcing may be activated randomly or in a pattern. Preferably, some of the energy generated from multiple source positions overlap in time and space to reduce acquisition time and cost.
  • source points within a distance of D source from the boundary are repeated for the new spread using the same relative time shifts. For example if sources 1201 , 1202 , and 1203 are within the distance D source from the boundary so that not all the energy is recorded within the distance D source from a source, these sources are repeated into the new spread position with the same relative time shift and the same encoding function previously used at that source position. For example 1221, 1222 and 1223 are a repeat of 1201, 1202 and 1203. This insures that the combined record has the captured the entire important signal, such as 1224, from each source.
  • Step 703 several pseudo super-source records of fixed size and duration are constructed.
  • the duration is longer than T listen but short enough to be efficiently simulated in the computer. Construction for this example may be iullustrated in two steps.
  • Fig. 13 the two records from Fig. 12 are combined to make a fixed spread. Because some shots near the boundary were repeated with the same relative timing and encoding function, we have insured that all the energy within a distance of D source and a time of T source are captured in the combined record.
  • isolated time windows are extracted preferably of length T source . There is not a need that the measured energy from any one source be isolated within the window, nor does any record need to start at the firing time of any source as with pseudo-deblending.
  • windows 1404 and 1406 are extracted from the record in Fig. 13 .
  • Step 704 the source locations, encoding function and firing time relative to zero time of the super source records 1401 and 1402 are identified.
  • the sources should be within a time of T source from above the top of the window or within a distance of D source from the boundaries of the record.
  • Step 705 the super-source records are computer encoded and summed, making a measured simultaneous source record as illustrated in Fig. 15 as 1501.
  • Step 706 the sourcing positions and both field and computer encoding functions are used to generate a simulated simultaneous source record as illustrated in 1502. This simulation is efficient, because all the sources are simulated simultaneously for a short period of time.
  • Step 707 the measured region of interest 1504 is compared to the simulated region of interest 1506 , and the result used to update a subsurface model. Because all the energy within the distance D source and T source from each source position is represented in both records, artifacts from the simultaneous record creation are avoided. With subsequent iterations of the imaging and inversion, preferably the super-source records would be combined with different encoding functions forming different realizations of all the data, reducing crosstalk noise.
  • Figure 16 shows a velocity model example in 2D. Due to patent law restrictions on the use of color, Fig. 16 is a black-and-white reproduction of a data display where velocity is quantitatively represented on a color scale.
  • the data were computer-simulated with a front and rear source and field encoding with random time delays up to 200 ms using the near-surface model shown in 1601.
  • Using the conventional method for simultaneous encoded-source inversion with the moving streamer data will not yield the correct results (not shown) because of the failure of the fixed-receiver assumption.
  • the super-source gathers were constructed as described in the marine embodiment, and a low-frequency inversion was performed using the present inventive method with the results shown at 1602. It can be seen that the model inferred by data inversion using the present inventive method compares very favorably with the "true" model 1601.

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Oceanography (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
EP18154578.1A 2013-08-23 2014-07-11 Inversion iterative des donnees sismiques codees a base de la construction des donnees pseudo supersource Withdrawn EP3351972A1 (fr)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201361869292P 2013-08-23 2013-08-23
PCT/US2014/046361 WO2015026451A2 (fr) 2013-08-23 2014-07-11 Recherche simultanée de sources pendant une acquisition sismique et une inversion sismique
EP14744730.4A EP3036566B1 (fr) 2013-08-23 2014-07-11 Recherche simultanée de sources pendant une acquisition sismique et une inversion sismique

Related Parent Applications (2)

Application Number Title Priority Date Filing Date
EP14744730.4A Division EP3036566B1 (fr) 2013-08-23 2014-07-11 Recherche simultanée de sources pendant une acquisition sismique et une inversion sismique
EP14744730.4A Division-Into EP3036566B1 (fr) 2013-08-23 2014-07-11 Recherche simultanée de sources pendant une acquisition sismique et une inversion sismique

Publications (1)

Publication Number Publication Date
EP3351972A1 true EP3351972A1 (fr) 2018-07-25

Family

ID=51230232

Family Applications (2)

Application Number Title Priority Date Filing Date
EP14744730.4A Not-in-force EP3036566B1 (fr) 2013-08-23 2014-07-11 Recherche simultanée de sources pendant une acquisition sismique et une inversion sismique
EP18154578.1A Withdrawn EP3351972A1 (fr) 2013-08-23 2014-07-11 Inversion iterative des donnees sismiques codees a base de la construction des donnees pseudo supersource

Family Applications Before (1)

Application Number Title Priority Date Filing Date
EP14744730.4A Not-in-force EP3036566B1 (fr) 2013-08-23 2014-07-11 Recherche simultanée de sources pendant une acquisition sismique et une inversion sismique

Country Status (9)

Country Link
US (1) US9772413B2 (fr)
EP (2) EP3036566B1 (fr)
AU (1) AU2014309376B2 (fr)
BR (1) BR112015030104A2 (fr)
CA (1) CA2913496C (fr)
DK (1) DK3036566T3 (fr)
ES (1) ES2671789T3 (fr)
MX (1) MX359179B (fr)
WO (1) WO2015026451A2 (fr)

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9188693B2 (en) * 2012-03-26 2015-11-17 Apache Corporation Method for acquiring marine seismic data
EP2856373A4 (fr) * 2012-05-24 2016-06-29 Exxonmobil Upstream Res Co Système et méthode de prédiction de la résistance d'une roche
US20150063064A1 (en) * 2013-09-03 2015-03-05 Pgs Geophysical As Methods and systems for attenuating noise in seismic data
US10036818B2 (en) 2013-09-06 2018-07-31 Exxonmobil Upstream Research Company Accelerating full wavefield inversion with nonstationary point-spread functions
EP3137926A4 (fr) 2014-04-28 2017-12-13 Westerngeco LLC Reconstruction de champs d'ondes
DE102014108789A1 (de) * 2014-06-24 2016-01-07 Byk-Gardner Gmbh Mehrstufiges Verfahren zur Untersuchung von Oberflächen sowie entsprechende Vorrichtung
US10838092B2 (en) * 2014-07-24 2020-11-17 Exxonmobil Upstream Research Company Estimating multiple subsurface parameters by cascaded inversion of wavefield components
US10073183B2 (en) 2014-10-20 2018-09-11 Pgs Geophysical As Methods and systems that attenuate noise in seismic data
US9784867B2 (en) 2015-04-01 2017-10-10 Schlumberger Technology Corporation Seismic data processing
US10928535B2 (en) 2015-05-01 2021-02-23 Reflection Marine Norge As Marine vibrator directive source survey
WO2016179206A1 (fr) 2015-05-05 2016-11-10 Schlumberger Technology Corporation Élimination d'effets d'acquisition dans des données sismiques marines
US10948615B2 (en) 2015-12-02 2021-03-16 Westerngeco L.L.C. Land seismic sensor spread with adjacent multicomponent seismic sensor pairs on average at least twenty meters apart
GB2547940A (en) * 2016-03-04 2017-09-06 Robertsson Johan Source separation method
WO2017218722A1 (fr) 2016-06-15 2017-12-21 Schlumberger Technology Corporation Systèmes et procédés pour atténuer le bruit dans des données sismiques et reconstruire des champs d'ondes sur la base des données sismiques
US10955576B2 (en) * 2016-08-19 2021-03-23 Halliburton Energy Services, Inc. Full waveform inversion of vertical seismic profile data for anisotropic velocities using pseudo-acoustic wave equations
GB2558630A (en) * 2017-01-12 2018-07-18 Seismic Apparition Gmbh Method for dealiasing data
WO2019008538A1 (fr) * 2017-07-06 2019-01-10 Chevron U.S.A. Inc. Système et procédé d'inversion de forme d'onde complète de données sismiques
WO2019017957A1 (fr) * 2017-07-20 2019-01-24 Halliburton Energy Services, Inc. Modélisation de dipôles destinée à des champs électriques et/ou magnétiques
US10788597B2 (en) 2017-12-11 2020-09-29 Saudi Arabian Oil Company Generating a reflectivity model of subsurface structures
WO2019152896A1 (fr) * 2018-02-02 2019-08-08 Fairfield Industries, Inc. Imagerie sismique avec une condition d'imagerie par décomposition temporelle
US11656377B2 (en) * 2018-03-30 2023-05-23 Cgg Services Sas Visco-acoustic full waveform inversion of velocity and Q
US11231516B2 (en) * 2018-05-15 2022-01-25 Exxonmobil Upstream Research Company Direct migration of simultaneous-source survey data
US11448790B2 (en) 2018-05-24 2022-09-20 King Abdullah University Of Science And Technology Method for partial differential equation inversion of data
US11372123B2 (en) 2019-10-07 2022-06-28 Exxonmobil Upstream Research Company Method for determining convergence in full wavefield inversion of 4D seismic data
US11914101B2 (en) 2020-01-31 2024-02-27 ExxonMobil Technology and Engineering Company Method for partitioning a search direction when using least squares reverse time migration
US11320557B2 (en) 2020-03-30 2022-05-03 Saudi Arabian Oil Company Post-stack time domain image with broadened spectrum
US11614555B2 (en) 2020-09-14 2023-03-28 China Petroleum & Chemical Corporation Method and system for connecting elements to sources and receivers during spectrum element method and finite element method seismic wave modeling
CN112415579A (zh) * 2020-11-03 2021-02-26 中国石油天然气集团有限公司 同时源随机激发方法、系统及装置
US20230194736A1 (en) * 2021-12-16 2023-06-22 Chevron U.S.A. Inc. System and method for robust seismic imaging
CN114296134B (zh) * 2021-12-24 2023-03-31 西安交通大学 一种深度卷积网络地震资料解混方法及系统
CN115060769B (zh) * 2022-06-07 2024-04-02 深圳大学 一种基于智能反演的隧道围岩裂隙及松动检测方法、系统

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5721710A (en) 1995-09-29 1998-02-24 Atlantic Richfield Company High fidelity vibratory source seismic method with source separation
US6021094A (en) 1998-12-03 2000-02-01 Sandia Corporation Method of migrating seismic records
US6906981B2 (en) 2002-07-17 2005-06-14 Pgs Americas, Inc. Method and system for acquiring marine seismic data using multiple seismic sources
WO2008042081A1 (fr) * 2006-09-28 2008-04-10 Exxonmobil Upstream Research Company Inversion itérative de données à partir de sources géophysiques simultanées
US20120073825A1 (en) 2010-09-27 2012-03-29 Routh Partha S Simultaneous Source Encoding and Source Separation As A Practical Solution For Full Wavefield Inversion
US20120143506A1 (en) 2010-12-01 2012-06-07 Routh Partha S Simultaneous Source Inversion for Marine Streamer Data With Cross-Correlation Objective Function
US8248886B2 (en) 2007-04-10 2012-08-21 Exxonmobil Upstream Research Company Separation and noise removal for multiple vibratory source seismic data
US20120215506A1 (en) 2011-02-18 2012-08-23 Rickett James E Waveform inversion by multiple shot-encoding for non-fixed spread geometries
US20120290214A1 (en) 2011-05-13 2012-11-15 Saudi Arabian Oil Company Coupled time-distance dependent swept frequency source acquisition design and data de-noising
US8437998B2 (en) 2010-09-27 2013-05-07 Exxonmobil Upstream Research Company Hybrid method for full waveform inversion using simultaneous and sequential source method
US20130238246A1 (en) 2012-03-08 2013-09-12 Jerome R. Krebs Orthogonal Source and Receiver Encoding

Family Cites Families (194)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3812457A (en) 1969-11-17 1974-05-21 Shell Oil Co Seismic exploration method
US3864667A (en) 1970-09-11 1975-02-04 Continental Oil Co Apparatus for surface wave parameter determination
US3984805A (en) 1973-10-18 1976-10-05 Daniel Silverman Parallel operation of seismic vibrators without phase control
US4168485A (en) 1974-08-12 1979-09-18 Continental Oil Company Simultaneous use of pseudo-random control signals in vibrational exploration methods
US4675851A (en) 1982-09-09 1987-06-23 Western Geophysical Co. Method for seismic exploration
US4545039A (en) 1982-09-09 1985-10-01 Western Geophysical Co. Of America Methods for seismic exploration
US4575830A (en) 1982-10-15 1986-03-11 Schlumberger Technology Corporation Indirect shearwave determination
US4594662A (en) 1982-11-12 1986-06-10 Schlumberger Technology Corporation Diffraction tomography systems and methods with fixed detector arrays
US4562540A (en) 1982-11-12 1985-12-31 Schlumberger Technology Corporation Diffraction tomography system and methods
FR2543306B1 (fr) 1983-03-23 1985-07-26 Elf Aquitaine Procede et dispositif pour l'optimisation des donnees sismiques
US4924390A (en) 1985-03-04 1990-05-08 Conoco, Inc. Method for determination of earth stratum elastic parameters using seismic energy
US4715020A (en) 1986-10-29 1987-12-22 Western Atlas International, Inc. Simultaneous performance of multiple seismic vibratory surveys
FR2589587B1 (fr) 1985-10-30 1988-02-05 Inst Francais Du Petrole Procede de prospection sismique marine utilisant un signal vibratoire code et dispositif pour sa mise en oeuvre
US4707812A (en) 1985-12-09 1987-11-17 Atlantic Richfield Company Method of suppressing vibration seismic signal correlation noise
US4823326A (en) 1986-07-21 1989-04-18 The Standard Oil Company Seismic data acquisition technique having superposed signals
US4686654A (en) 1986-07-31 1987-08-11 Western Geophysical Company Of America Method for generating orthogonal sweep signals
US4766574A (en) 1987-03-31 1988-08-23 Amoco Corporation Method for depth imaging multicomponent seismic data
US4953657A (en) 1987-11-30 1990-09-04 Halliburton Geophysical Services, Inc. Time delay source coding
US4969129A (en) 1989-09-20 1990-11-06 Texaco Inc. Coding seismic sources
US4982374A (en) 1989-10-23 1991-01-01 Halliburton Geophysical Services, Inc. Method of source coding and harmonic cancellation for vibrational geophysical survey sources
GB9011836D0 (en) 1990-05-25 1990-07-18 Mason Iain M Seismic surveying
US6005916A (en) 1992-10-14 1999-12-21 Techniscan, Inc. Apparatus and method for imaging with wavefields using inverse scattering techniques
US5469062A (en) 1994-03-11 1995-11-21 Baker Hughes, Inc. Multiple depths and frequencies for simultaneous inversion of electromagnetic borehole measurements
GB2322704B (en) 1994-07-07 1998-12-09 Geco As Method of Processing seismic data
US5583825A (en) 1994-09-02 1996-12-10 Exxon Production Research Company Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data
US5924049A (en) 1995-04-18 1999-07-13 Western Atlas International, Inc. Methods for acquiring and processing seismic data
DE69625978T2 (de) 1995-04-18 2003-11-13 Western Atlas Int Inc Gleichmässige unteroberflächendeckung bei anwesenheit von steilen neigungen
US5719821A (en) 1995-09-29 1998-02-17 Atlantic Richfield Company Method and apparatus for source separation of seismic vibratory signals
US5790473A (en) 1995-11-13 1998-08-04 Mobil Oil Corporation High fidelity vibratory source seismic method for use in vertical seismic profile data gathering with a plurality of vibratory seismic energy sources
US5715213A (en) 1995-11-13 1998-02-03 Mobil Oil Corporation High fidelity vibratory source seismic method using a plurality of vibrator sources
US5822269A (en) 1995-11-13 1998-10-13 Mobil Oil Corporation Method for separation of a plurality of vibratory seismic energy source signals
US5838634A (en) 1996-04-04 1998-11-17 Exxon Production Research Company Method of generating 3-D geologic models incorporating geologic and geophysical constraints
US5798982A (en) 1996-04-29 1998-08-25 The Trustees Of Columbia University In The City Of New York Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models
GB9612471D0 (en) 1996-06-14 1996-08-14 Geco As Method and apparatus for multiple seismic vibratory surveys
US5878372A (en) 1997-03-04 1999-03-02 Western Atlas International, Inc. Method for simultaneous inversion processing of well log data using a plurality of earth models
US6014342A (en) 1997-03-21 2000-01-11 Tomo Seis, Inc. Method of evaluating a subsurface region using gather sensitive data discrimination
US5999489A (en) 1997-03-21 1999-12-07 Tomoseis Inc. High vertical resolution crosswell seismic imaging
US5920828A (en) 1997-06-02 1999-07-06 Baker Hughes Incorporated Quality control seismic data processing system
FR2765692B1 (fr) 1997-07-04 1999-09-10 Inst Francais Du Petrole Methode pour modeliser en 3d l'impedance d'un milieu heterogene
GB2329043B (en) 1997-09-05 2000-04-26 Geco As Method of determining the response caused by model alterations in seismic simulations
US5999488A (en) 1998-04-27 1999-12-07 Phillips Petroleum Company Method and apparatus for migration by finite differences
US6219621B1 (en) 1998-06-30 2001-04-17 Exxonmobil Upstream Research Co. Sparse hyperbolic inversion of seismic data
US6388947B1 (en) 1998-09-14 2002-05-14 Tomoseis, Inc. Multi-crosswell profile 3D imaging and method
FR2784195B1 (fr) 1998-10-01 2000-11-17 Inst Francais Du Petrole Methode pour realiser en 3d avant sommation, une migration de donnees sismiques
US6574564B2 (en) 1998-10-01 2003-06-03 Institut Francais Du Petrole 3D prestack seismic data migration method
US6225803B1 (en) 1998-10-29 2001-05-01 Baker Hughes Incorporated NMR log processing using wavelet filter and iterative inversion
US6754588B2 (en) 1999-01-29 2004-06-22 Platte River Associates, Inc. Method of predicting three-dimensional stratigraphy using inverse optimization techniques
CA2362285C (fr) 1999-02-12 2005-06-14 Schlumberger Canada Limited Modelisation de zone souterraine a incertitude reduite
US6058073A (en) 1999-03-30 2000-05-02 Atlantic Richfield Company Elastic impedance estimation for inversion of far offset seismic sections
FR2792419B1 (fr) 1999-04-16 2001-09-07 Inst Francais Du Petrole Methode pour obtenir un modele optimal d'une caracteristique physique dans un milieu heterogene, tel que le sous-sol
GB9927395D0 (en) 1999-05-19 2000-01-19 Schlumberger Holdings Improved seismic data acquisition method
US6327537B1 (en) 1999-07-19 2001-12-04 Luc T. Ikelle Multi-shooting approach to seismic modeling and acquisition
FR2798197B1 (fr) 1999-09-02 2001-10-05 Inst Francais Du Petrole Methode pour former un modele d'une formation geologique, contraint par des donnees dynamiques et statiques
DK1746443T3 (en) 1999-10-22 2014-03-17 Fugro N V A method of calculating the elastic parameters and stone composition of subterranean formations using seismic data
US6480790B1 (en) 1999-10-29 2002-11-12 Exxonmobil Upstream Research Company Process for constructing three-dimensional geologic models having adjustable geologic interfaces
FR2800473B1 (fr) 1999-10-29 2001-11-30 Inst Francais Du Petrole Methode pour modeliser en 2d ou 3d un milieu heterogene tel que le sous-sol decrit par plusieurs parametres physiques
US6903999B2 (en) 2000-01-21 2005-06-07 Schlumberger Technology Corporation System and method for estimating seismic material properties
AU779802B2 (en) 2000-01-21 2005-02-10 Schlumberger Holdings Limited System and method for seismic wavefield separation
US6826486B1 (en) 2000-02-11 2004-11-30 Schlumberger Technology Corporation Methods and apparatus for predicting pore and fracture pressures of a subsurface formation
FR2805051B1 (fr) 2000-02-14 2002-12-06 Geophysique Cie Gle Methode de surveillance sismique d'une zone souterraine par utilisation simultanee de plusieurs sources vibrosismiques
GB2359363B (en) 2000-02-15 2002-04-03 Geco Prakla Processing simultaneous vibratory seismic data
US6687659B1 (en) 2000-03-24 2004-02-03 Conocophillips Company Method and apparatus for absorbing boundary conditions in numerical finite-difference acoustic applications
US6317695B1 (en) 2000-03-30 2001-11-13 Nutec Sciences, Inc. Seismic data processing method
EP1327162A2 (fr) 2000-10-17 2003-07-16 WesternGeco, L.L.C. Procede d'utilisation de balayages en cascade pour le codage de sources et la suppression d'harmoniques
US20020120429A1 (en) 2000-12-08 2002-08-29 Peter Ortoleva Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories
FR2818753B1 (fr) 2000-12-21 2003-03-21 Inst Francais Du Petrole Methode et dispositif de prospection sismique par emission simultanee de signaux sismisques obtenus en codant un signal par des sequences pseudo aleatoires
FR2821677B1 (fr) 2001-03-05 2004-04-30 Geophysique Cie Gle Perfectionnements aux procedes d'inversion tomographique d'evenements pointes sur les donnees sismiques migrees
US6751558B2 (en) 2001-03-13 2004-06-15 Conoco Inc. Method and process for prediction of subsurface fluid and rock pressures in the earth
US6927698B2 (en) 2001-08-27 2005-08-09 Larry G. Stolarczyk Shuttle-in receiver for radio-imaging underground geologic structures
US6545944B2 (en) 2001-05-30 2003-04-08 Westerngeco L.L.C. Method for acquiring and processing of data from two or more simultaneously fired sources
US6882958B2 (en) 2001-06-28 2005-04-19 National Instruments Corporation System and method for curve fitting using randomized techniques
GB2379013B (en) 2001-08-07 2005-04-20 Abb Offshore Systems Ltd Microseismic signal processing
US6593746B2 (en) 2001-08-27 2003-07-15 Larry G. Stolarczyk Method and system for radio-imaging underground geologic structures
US7672824B2 (en) 2001-12-10 2010-03-02 Westerngeco L.L.C. Method for shallow water flow detection
US7069149B2 (en) 2001-12-14 2006-06-27 Chevron U.S.A. Inc. Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume
US7330799B2 (en) 2001-12-21 2008-02-12 Société de commercialisation des produits de la recherche appliquée-Socpra Sciences et Génie s.e.c. Method and algorithm for using surface waves
US6842701B2 (en) 2002-02-25 2005-01-11 Westerngeco L.L.C. Method of noise removal for cascaded sweep data
GB2387226C (en) 2002-04-06 2008-05-12 Westerngeco Ltd A method of seismic surveying
FR2839368B1 (fr) 2002-05-06 2004-10-01 Total Fina Elf S A Methode de decimation de traces sismiques pilotee par le trajet sismique
US6832159B2 (en) 2002-07-11 2004-12-14 Schlumberger Technology Corporation Intelligent diagnosis of environmental influence on well logs with model-based inversion
FR2843202B1 (fr) 2002-08-05 2004-09-10 Inst Francais Du Petrole Methode pour former un modele representatif de la distribution d'une grandeur physique dans une zone souterraine, affranchi de l'effet de bruits correles entachant des donnees d'exploration
WO2004034088A2 (fr) 2002-10-04 2004-04-22 Paradigm Geophysical Corporation Procede et systeme permettant l'imagerie sismique en frequences limitees
GB2396448B (en) 2002-12-21 2005-03-02 Schlumberger Holdings System and method for representing and processing and modeling subterranean surfaces
US20040225483A1 (en) 2003-02-24 2004-11-11 Michal Okoniewski Fdtd hardware acceleration system
US6735527B1 (en) 2003-02-26 2004-05-11 Landmark Graphics Corporation 3-D prestack/poststack multiple prediction
US6999880B2 (en) 2003-03-18 2006-02-14 The Regents Of The University Of California Source-independent full waveform inversion of seismic data
US7184367B2 (en) 2003-03-27 2007-02-27 Exxonmobil Upstream Research Company Method to convert seismic traces into petrophysical property logs
US7436734B2 (en) 2003-04-01 2008-10-14 Exxonmobil Upstream Research Co. Shaped high frequency vibratory source
US7072767B2 (en) 2003-04-01 2006-07-04 Conocophillips Company Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data
NO322089B1 (no) 2003-04-09 2006-08-14 Norsar V Daglig Leder Fremgangsmate for simulering av lokale prestakk dypmigrerte seismiske bilder
GB2400438B (en) 2003-04-11 2005-06-01 Westerngeco Ltd Determination of waveguide parameters
US6970397B2 (en) 2003-07-09 2005-11-29 Gas Technology Institute Determination of fluid properties of earth formations using stochastic inversion
US6882938B2 (en) 2003-07-30 2005-04-19 Pgs Americas, Inc. Method for separating seismic signals from two or more distinct sources
GB2405473B (en) 2003-08-23 2005-10-05 Westerngeco Ltd Multiple attenuation method
US6944546B2 (en) 2003-10-01 2005-09-13 Halliburton Energy Services, Inc. Method and apparatus for inversion processing of well logging data in a selected pattern space
US6901333B2 (en) 2003-10-27 2005-05-31 Fugro N.V. Method and device for the generation and application of anisotropic elastic parameters
US7046581B2 (en) 2003-12-01 2006-05-16 Shell Oil Company Well-to-well tomography
US20050128874A1 (en) * 2003-12-15 2005-06-16 Chevron U.S.A. Inc. Methods for acquiring and processing seismic data from quasi-simultaneously activated translating energy sources
US7359283B2 (en) 2004-03-03 2008-04-15 Pgs Americas, Inc. System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers
US7791980B2 (en) 2004-05-21 2010-09-07 Westerngeco L.L.C. Interpolation and extrapolation method for seismic recordings
FR2872584B1 (fr) 2004-06-30 2006-08-11 Inst Francais Du Petrole Methode pour simuler le depot sedimentaire dans un bassin respectant les epaisseurs des sequences sedimentaires
EP1617309B1 (fr) 2004-07-15 2011-01-12 Fujitsu Limited Technique de simulation utilisant le raffinement de maillage espace-temps
US7646924B2 (en) 2004-08-09 2010-01-12 David Leigh Donoho Method and apparatus for compressed sensing
US7480206B2 (en) 2004-09-13 2009-01-20 Chevron U.S.A. Inc. Methods for earth modeling and seismic imaging using interactive and selective updating
FR2876458B1 (fr) 2004-10-08 2007-01-19 Geophysique Cie Gle Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples
GB2422433B (en) 2004-12-21 2008-03-19 Sondex Wireline Ltd Method and apparatus for determining the permeability of earth formations
US7373251B2 (en) 2004-12-22 2008-05-13 Marathon Oil Company Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data
US7230879B2 (en) 2005-02-12 2007-06-12 Chevron U.S.A. Inc. Method and apparatus for true relative amplitude correction of seismic data for normal moveout stretch effects
EP1859301B1 (fr) 2005-02-22 2013-07-17 Paradigm Geophysical Ltd. Suppressions multiples de migration en profondeur et dans le temps dans le domaine angulaire
US7840625B2 (en) 2005-04-07 2010-11-23 California Institute Of Technology Methods for performing fast discrete curvelet transforms of data
WO2006122146A2 (fr) 2005-05-10 2006-11-16 William Marsh Rice University Procede et appareil utilisant la technique du 'compressed sensing' distribue
US7405997B2 (en) 2005-08-11 2008-07-29 Conocophillips Company Method of accounting for wavelet stretch in seismic data
US20090164756A1 (en) 2005-10-18 2009-06-25 Tor Dokken Geological Response Data Imaging With Stream Processors
AU2006235820B2 (en) 2005-11-04 2008-10-23 Westerngeco Seismic Holdings Limited 3D pre-stack full waveform inversion
FR2895091B1 (fr) 2005-12-21 2008-02-22 Inst Francais Du Petrole Methode pour mettre a jour un modele geologique par des donnees sismiques
US7400552B2 (en) * 2006-01-19 2008-07-15 Westerngeco L.L.C. Methods and systems for efficiently acquiring towed streamer seismic surveys
GB2436626B (en) 2006-03-28 2008-08-06 Westerngeco Seismic Holdings Method of evaluating the interaction between a wavefield and a solid body
US7620534B2 (en) 2006-04-28 2009-11-17 Saudi Aramco Sound enabling computerized system for real time reservoir model calibration using field surveillance data
US20070274155A1 (en) 2006-05-25 2007-11-29 Ikelle Luc T Coding and Decoding: Seismic Data Modeling, Acquisition and Processing
US7725266B2 (en) 2006-05-31 2010-05-25 Bp Corporation North America Inc. System and method for 3D frequency domain waveform inversion based on 3D time-domain forward modeling
US7599798B2 (en) 2006-09-11 2009-10-06 Westerngeco L.L.C. Migrating composite seismic response data to produce a representation of a seismic volume
RU2008151147A (ru) 2006-12-07 2010-06-27 Каусел Оф Сайнтифик Энд Индастриал Рисерч (In) Способ вычисления точного импульсного отклика плоского акустического отражателя для точечного акустического источника при нулевом смещении
ATE543109T1 (de) 2007-01-20 2012-02-15 Spectraseis Ag Zeitumkehr-reservoir-lokalisierung
US7715986B2 (en) 2007-05-22 2010-05-11 Chevron U.S.A. Inc. Method for identifying and removing multiples for imaging with beams
US7974824B2 (en) 2007-06-29 2011-07-05 Westerngeco L. L. C. Seismic inversion of data containing surface-related multiples
JP2009063942A (ja) 2007-09-10 2009-03-26 Sumitomo Electric Ind Ltd 遠赤外線カメラ用レンズ、レンズユニット及び撮像装置
US20090070042A1 (en) 2007-09-11 2009-03-12 Richard Birchwood Joint inversion of borehole acoustic radial profiles for in situ stresses as well as third-order nonlinear dynamic moduli, linear dynamic elastic moduli, and static elastic moduli in an isotropically stressed reference state
US20090083006A1 (en) 2007-09-20 2009-03-26 Randall Mackie Methods and apparatus for three-dimensional inversion of electromagnetic data
WO2009067041A1 (fr) 2007-11-19 2009-05-28 Steklov Mathematical Institute Ras Procédé et système pour évaluer les propriétés caractéristiques de deux milieux en contact et de l'interface entre eux à partir d'ondes de surface mélangées se propageant le long de l'interface
US20090164186A1 (en) 2007-12-20 2009-06-25 Bhp Billiton Innovation Pty Ltd. Method for determining improved estimates of properties of a model
EP2238474A4 (fr) 2008-01-08 2018-06-20 Exxonmobil Upstream Research Company Inversion de forme spectrale et migration de données sismiques
US8577660B2 (en) 2008-01-23 2013-11-05 Schlumberger Technology Corporation Three-dimensional mechanical earth modeling
EP2260331B1 (fr) 2008-03-21 2017-10-11 Exxonmobil Upstream Research Company Procédé efficace pour inversion de données géophysiques
US8451684B2 (en) * 2008-03-28 2013-05-28 Exxonmobil Upstream Research Company Surface wave mitigation in spatially inhomogeneous media
EP2105765A1 (fr) 2008-03-28 2009-09-30 Schlumberger Holdings Limited Inversion simultanée de données d'induction pour la constante diélectrique et la conductivité électrique
US8275592B2 (en) 2008-04-07 2012-09-25 Westerngeco L.L.C. Joint inversion of time domain controlled source electromagnetic (TD-CSEM) data and further data
US8494777B2 (en) 2008-04-09 2013-07-23 Schlumberger Technology Corporation Continuous microseismic mapping for real-time 3D event detection and location
US8345510B2 (en) 2008-06-02 2013-01-01 Pgs Geophysical As Method for aquiring and processing marine seismic data to extract and constructively use the up-going and down-going wave-fields emitted by the source(s)
US20110182141A1 (en) 2008-08-14 2011-07-28 Schlumberger Technology Corporation Method and system for monitoring a logging tool position in a borehole
US8295124B2 (en) 2008-08-15 2012-10-23 Bp Corporation North America Inc. Method for separating independent simultaneous sources
US8559270B2 (en) 2008-08-15 2013-10-15 Bp Corporation North America Inc. Method for separating independent simultaneous sources
US20100054082A1 (en) 2008-08-29 2010-03-04 Acceleware Corp. Reverse-time depth migration with reduced memory requirements
US8296069B2 (en) 2008-10-06 2012-10-23 Bp Corporation North America Inc. Pseudo-analytical method for the solution of wave equations
US7616523B1 (en) 2008-10-22 2009-11-10 Pgs Geophysical As Method for combining pressure and motion seismic signals from streamers where sensors are not at a common depth
US9213119B2 (en) 2008-10-29 2015-12-15 Conocophillips Company Marine seismic acquisition
US20100118651A1 (en) 2008-11-10 2010-05-13 Chevron U.S.A. Inc. Method for generation of images related to a subsurface region of interest
US20100142316A1 (en) 2008-12-07 2010-06-10 Henk Keers Using waveform inversion to determine properties of a subsurface medium
US8095345B2 (en) 2009-01-20 2012-01-10 Chevron U.S.A. Inc Stochastic inversion of geophysical data for estimating earth model parameters
US9052410B2 (en) 2009-02-12 2015-06-09 Conocophillips Company Multiple seismic signal inversion
WO2010095859A2 (fr) 2009-02-17 2010-08-26 Shin Changsoo Dispositif et procédé pour l'imagerie de structure de subsurface
US8352190B2 (en) 2009-02-20 2013-01-08 Exxonmobil Upstream Research Company Method for analyzing multiple geophysical data sets
US9075163B2 (en) 2009-04-17 2015-07-07 Westerngeco L.L.C. Interferometric seismic data processing
US8176284B2 (en) 2009-08-11 2012-05-08 Texas Memory Systems, Inc. FLASH-based memory system with variable length page stripes including data protection information
US20110044127A1 (en) 2009-08-19 2011-02-24 Clement Kostov Removing free-surface effects from seismic data acquired in a towed survey
US8923093B2 (en) 2009-08-25 2014-12-30 Westerngeco L.L.C. Determining the quality of a seismic inversion
US20110131020A1 (en) 2009-09-09 2011-06-02 Conocophillips Company Dip guided full waveform inversion
GB2486121B (en) 2009-10-01 2014-08-13 Halliburton Energy Serv Inc Apparatus and methods of locating downhole anomalies
US9244181B2 (en) 2009-10-19 2016-01-26 Westerngeco L.L.C. Full-waveform inversion in the traveltime domain
US8861308B2 (en) 2009-12-07 2014-10-14 Westerngeco L.L.C. Simultaneous joint inversion of surface wave and refraction data
WO2011091216A2 (fr) 2010-01-22 2011-07-28 Schlumberger Canada Limited Évaluation en temps réel d'anisotropie et d'inclinaison de formation utilisant des mesures d'induction triaxiales
CA2787693A1 (fr) 2010-01-25 2011-07-28 CGGVeritas Services (U.S.) Inc. Procedes et systemes d'estimation de contraintes a l'aide de donnees sismiques
US20130098608A1 (en) 2010-01-29 2013-04-25 Robert Barnum Temporary field storage of gas to optimize field development
US8265875B2 (en) 2010-01-29 2012-09-11 Westerngeco L.L.C. Interpolation of periodic data
US8537638B2 (en) 2010-02-10 2013-09-17 Exxonmobil Upstream Research Company Methods for subsurface parameter estimation in full wavefield inversion and reverse-time migration
EP2545399A4 (fr) 2010-03-12 2017-10-25 CGG Veritas Services (U.S.) Inc. Procédés et systèmes pour réaliser une inversion élastique simultanée azimutale
US8680865B2 (en) 2010-03-19 2014-03-25 Schlumberger Technology Corporation Single well reservoir imaging apparatus and methods
US20110235464A1 (en) 2010-03-24 2011-09-29 John Brittan Method of imaging the earth's subsurface during marine seismic data acquisition
US8223587B2 (en) 2010-03-29 2012-07-17 Exxonmobil Upstream Research Company Full wavefield inversion using time varying filters
US9176244B2 (en) 2010-03-31 2015-11-03 Schlumberger Technology Corporation Data set inversion using source-receiver compression
US8576663B2 (en) 2010-04-30 2013-11-05 Schlumberger Technology Corporation Multicomponent seismic inversion of VSP data
KR101167715B1 (ko) 2010-04-30 2012-07-20 서울대학교산학협력단 복소 구배 최소자승법에의한 파형 역산을 이용한 지하 구조의 영상화 장치 및 방법
US8694299B2 (en) 2010-05-07 2014-04-08 Exxonmobil Upstream Research Company Artifact reduction in iterative inversion of geophysical data
US8756042B2 (en) 2010-05-19 2014-06-17 Exxonmobile Upstream Research Company Method and system for checkpointing during simulations
JP5812990B2 (ja) 2010-06-15 2015-11-17 電気化学工業株式会社 透光性硬質基板積層体の製造方法
US20110320180A1 (en) 2010-06-29 2011-12-29 Al-Saleh Saleh M Migration Velocity Analysis of Seismic Data Using Common Image Cube and Green's Functions
US8612188B2 (en) 2010-07-12 2013-12-17 The University Of Manchester Wave modelling
EA026517B1 (ru) * 2010-08-02 2017-04-28 Бп Корпорейшн Норт Америка Инк. Способ сейсмической разведки
EP2606452A4 (fr) 2010-08-16 2017-08-16 Exxonmobil Upstream Research Company Réduction de la dimensionnalité du problème de l'inversion conjointe
US20120051176A1 (en) 2010-08-31 2012-03-01 Chevron U.S.A. Inc. Reverse time migration back-scattering noise removal using decomposed wavefield directivity
WO2012047384A1 (fr) 2010-09-27 2012-04-12 Exxonmobil Upstream Research Company Procédé hybride d'inversion de forme d'onde complète à l'aide d'un procédé de source simultanée et séquentielle
WO2012041834A1 (fr) 2010-09-28 2012-04-05 Shell Internationale Research Maatschappij B.V. Estimation de modèle terrestre par une inversion de forme d'onde complète acoustique de données sismiques
US9134442B2 (en) 2010-12-16 2015-09-15 Bp Corporation North America Inc. Seismic acquisition using narrowband seismic sources
IT1404170B1 (it) * 2011-02-10 2013-11-15 Eni Spa Metodo di indagine sismica del sottosuolo
EP2691795A4 (fr) 2011-03-30 2015-12-09 Vitesse de convergence d'une inversion d'un champ d'onde complet utilisant une mise en forme spectrale
US20120275267A1 (en) 2011-04-26 2012-11-01 Ramesh Neelamani Seismic Data Processing
US20120316790A1 (en) 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for data inversion with phase extrapolation
US20120316844A1 (en) 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for data inversion with phase unwrapping
US20120316791A1 (en) 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for seismic data inversion by non-linear model update
US9075159B2 (en) 2011-06-08 2015-07-07 Chevron U.S.A., Inc. System and method for seismic data inversion
US9075162B2 (en) * 2011-11-10 2015-07-07 Pgs Geophysical As Method and system for separating seismic sources in marine simultaneous shooting acquisition
US9091788B2 (en) 2011-11-14 2015-07-28 Cggveritas Services Sa Device and method for de-blending simultaneous shooting data with apex shifted radon transform
US9103943B2 (en) 2011-11-28 2015-08-11 Fugro-Geoteam As Acquisition and processing of multi-source broadband marine seismic data
US9176930B2 (en) 2011-11-29 2015-11-03 Exxonmobil Upstream Research Company Methods for approximating hessian times vector operation in full wavefield inversion
US20130311149A1 (en) 2012-05-17 2013-11-21 Yaxun Tang Tomographically Enhanced Full Wavefield Inversion

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5721710A (en) 1995-09-29 1998-02-24 Atlantic Richfield Company High fidelity vibratory source seismic method with source separation
US6021094A (en) 1998-12-03 2000-02-01 Sandia Corporation Method of migrating seismic records
US6906981B2 (en) 2002-07-17 2005-06-14 Pgs Americas, Inc. Method and system for acquiring marine seismic data using multiple seismic sources
WO2008042081A1 (fr) * 2006-09-28 2008-04-10 Exxonmobil Upstream Research Company Inversion itérative de données à partir de sources géophysiques simultanées
US8121823B2 (en) 2006-09-28 2012-02-21 Exxonmobil Upstream Research Company Iterative inversion of data from simultaneous geophysical sources
US8248886B2 (en) 2007-04-10 2012-08-21 Exxonmobil Upstream Research Company Separation and noise removal for multiple vibratory source seismic data
US20120073825A1 (en) 2010-09-27 2012-03-29 Routh Partha S Simultaneous Source Encoding and Source Separation As A Practical Solution For Full Wavefield Inversion
US8437998B2 (en) 2010-09-27 2013-05-07 Exxonmobil Upstream Research Company Hybrid method for full waveform inversion using simultaneous and sequential source method
US20120143506A1 (en) 2010-12-01 2012-06-07 Routh Partha S Simultaneous Source Inversion for Marine Streamer Data With Cross-Correlation Objective Function
US20120215506A1 (en) 2011-02-18 2012-08-23 Rickett James E Waveform inversion by multiple shot-encoding for non-fixed spread geometries
US20120290214A1 (en) 2011-05-13 2012-11-15 Saudi Arabian Oil Company Coupled time-distance dependent swept frequency source acquisition design and data de-noising
US20130238246A1 (en) 2012-03-08 2013-09-12 Jerome R. Krebs Orthogonal Source and Receiver Encoding

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
BEASLEY ET AL.: "A 3D simultaneous source field test processed using alternating projections: a new active separation method", GEOPHYSICAL PROSPECTING, vol. 60, 2012, pages 591 - 601
HABER ET AL.: "An effective method for parameter estimation with PDE constraints with multiple right hand sides", PREPRINT - UBC, 2010, Retrieved from the Internet <URL:http://www.math.ubc.ca/-haber/pubs/PdeOptStochV5.pdf>
HAFEDH BEN-HADJ-ALI ET AL: "An efficient frequency-domain full waveform inversion method using simultaneous encoded sources", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS, US, vol. 76, no. 4, 1 July 2011 (2011-07-01), pages R109 - R124, XP001574430, ISSN: 0016-8033, [retrieved on 20110628], DOI: 10.1190/1.3581357 *
TREADGOLD ET AL.: "Implementing A wide Azimuth Towed Streamer Field Trial, The What, Why, and Mostly How of WATS in Southern Green Canyon", SEG EXPANDED ABSTRACTS, 2006, pages 2901 - 2903
VIRIEUX J ET AL: "An overview of full-waveform inversion in exploration geophysics", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS, US, vol. 74, no. Suppl. of 6, 1 November 2009 (2009-11-01), pages WCC1 - WCC26, XP001550475, ISSN: 0016-8033, DOI: 10.1190/1.3238367 *

Also Published As

Publication number Publication date
ES2671789T3 (es) 2018-06-08
CA2913496A1 (fr) 2015-02-26
AU2014309376B2 (en) 2016-11-17
CA2913496C (fr) 2018-08-14
EP3036566A2 (fr) 2016-06-29
BR112015030104A2 (pt) 2017-07-25
WO2015026451A2 (fr) 2015-02-26
DK3036566T3 (en) 2018-07-23
MX359179B (es) 2018-09-18
WO2015026451A3 (fr) 2015-04-30
US20150057938A1 (en) 2015-02-26
AU2014309376A1 (en) 2016-03-03
US9772413B2 (en) 2017-09-26
EP3036566B1 (fr) 2018-04-04

Similar Documents

Publication Publication Date Title
CA2913496C (fr) Recherche simultanee de sources pendant une acquisition sismique et une inversion sismique
EP3025170B1 (fr) Procédé de désignature de données sismiques acquises à l&#39;aide d&#39;une source mobile et système de traitement de données sismiques
EP2786176B1 (fr) Séparation de données de source simultanée
EP3120167B1 (fr) Reconstruction de données sismiques de sources impulsionnelles à partir de données de réseaux de canons à air à déclenchements répartis dans le temps
US8982663B2 (en) Subsurface imaging systems and methods with multi-source survey component segregation and redetermination
US9551800B2 (en) Device and method for deblending simultaneous shooting data using annihilation filter
US11269093B2 (en) Method and system for generating geophysical data
US10996359B2 (en) Removal of acquisition effects from marine seismic data
WO2011049921A2 (fr) Procédés destinés à traiter des données sismiques contaminées par une énergie cohérente rayonnée par plusieurs sources
Beasley et al. Wave equation receiver deghosting: A provocative example
US20150234066A1 (en) Method and device for de-blending seismic data using source signature
RU2570827C2 (ru) Гибридный способ для полноволновой инверсии с использованием способа одновременных и последовательных источников
Huff et al. Near offset reconstruction for marine seismic data using a convolutional neural network
NO20221017A1 (en) Seismic data recording and processing with different uncontaminated recording time lengths

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

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

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20180201

AC Divisional application: reference to earlier application

Ref document number: 3036566

Country of ref document: EP

Kind code of ref document: P

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

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

Free format text: STATUS: THE APPLICATION HAS BEEN WITHDRAWN

18W Application withdrawn

Effective date: 20200511