• 专利附录
  • 专利说明
  • 权利要求
  • 类似技术
  • 同族专利
  • 引用文献
  • 法律状态
  • 优先权
    • 专利摘要:
    A method for generating an image of a subsurface feature, comprises providing seismic data containing information about the feature, comprising i) a first dataset comprising a first up-going wave-field and a first down-going wave-field, ii) a second dataset comprising a second up-going wave-field and a second down-going wave-field, wherein the second dataset is collected at a time that is later than the first dataset by a time interval, creating an up-going 4D difference dataset U by subtracting one of the first and second up-going wave-fields from the other and creating a down-going 4D difference dataset D by subtracting one of the first and second down-going wave-fields from the other, creating a weighting function W that is a function of the similarity of the difference datasets, and creating an image of the feature by generating a 4D similarity stack, where the 4D similarity stack is defined as (U+D)*W/2.
    公开号:AU2012362876A1
    申请号:U2012362876
    申请日:2012-12-17
    公开日:2014-07-03
    发明作者:Paul James Hatchell；Kees HORNMAN；Jorge Louis Lopez
    申请人:Shell Internationale Research Maatschappij BV；
    IPC主号:G01V1-28
    专利说明:
    WO 2013/101518 PCT/US2012/070079 SUPPRESSING 4D-NOISE BY WEIGHTED STACKING OF SIMULTANEOUSLY ACQUIRED WAVE-FIELDS RELATED CASES 5 This application claims priority from U.S. application Serial No. 61/580,693, filed on 28 December 2011 and incorporated herein by reference. FIELD OF THE INVENTION [0001] The invention relates to seismic signal processing. More particularly, the invention relates to a method for producing two independent but simultaneous images of the subsurface, 10 e.g., using up- and down-going wave-fields, and using a weighted stack of the up-going and down-going 4D differences based on their similarity to produce a reduced-noise image. BACKGROUND OF THE INVENTION [0002] In the field of seismic data processing, "4D-noise" is an expression used to describe non-repeatable noises recorded during seismic surveys that do not disappear when time-lapse 15 seismic signals are compared. In deepwater seismic applications, 4D-noise can result from many sources, including imperfect repeats of acquisition geometry, differences in multiple energy due to changes in water velocities, seismic interference noise from other vessels, swell noise, etc. Likewise, 4D-noise onshore can result from seasonal changes in the shallow subsurface, time-varying receiver ghosts, etc. 20 [0003] In deepwater environments where seismic recordings are made using ocean bottom sensors such as nodes (OBN) or cables (OBC), it is possible to produce two independent images of the subsurface using the conventional up-going wave-field and the down-going wave-field, respectively, via the well-known technique of mirror migration. It is useful to think about how the 4D-signal and 4D-noise compare within these two types of images. The 25 4D-signals that result from changes in reservoir properties should be the same on both the up going and down-going data once the signals have been corrected for differences in wavelets and effective fold. In contrast, 4D-noise behaves quite differently. For example, 4D noise due to multiple contamination as a result of differences in water velocities will arrive at different times relative to a horizon of interest on the up- and down-going wave-fields. 30 Similar issues arise in any instance in which it is desirable to generate a single image from two independent but simultaneous images of the subsurface, namely the signal in both images will be the same but the noise in both images will be different. 1 WO 2013/101518 PCT/US2012/070079 [0004] One way of take advantage of the difference in noise signals is to stack the up- and down-going wave-fields, which effectively doubles the fold of the survey and statistically suppresses noises by a factor of 1if/ 7 . This method does not always give the desired level of enhancement, however and it remains desirable to provide an improved technique. 5 SUMMARY OF THE INVENTION [0005] The present invention allows reduced-noise seismic images to be generated from 4D seismic signals. In some embodiments, 4D-noise is suppressed by stacking the up-going and down-going 4D differences after applying a weighting function that is based on how similar 10 the differences are to one another. [0006] According to some embodiments, a method for generating an image of a subsurface feature, comprises the steps of a) providing seismic data containing information about the feature, comprising a first dataset comprising a first up-going wave-field and a first down going wave-field (U 1 , D 1 ) and a second dataset comprising a second up-going wave-field and 15 a second down-going wave-field (U 2 , D 2 ), wherein the second dataset is collected at a time that is later than the first dataset by a time interval, b) creating an up-going 4D difference dataset U by subtracting one of the first and second up-going wave-fields U 1 , U 2 from the other and creating a down-going 4D difference dataset D by subtracting one of the first and second down-going wave-fields DI, D 2 from the other, c) creating a weighting function W that 20 is a function of the similarity of the up-going 4D difference dataset and the down-going 4D difference dataset D; and e) creating an image of the feature by generating a 4D similarity stack, where the 4D similarity stack is defined as (U+D)*W/2. [0007] The weighting function W may be selected from the functions of the group consisting of NRMSD, (1-NRMSD(U,D)/NRMSDcutoff , NRMSD 2 , average difference of absolute 25 values, manual weighting, and powers of the foregoing functions. In preferred embodiments, the weighting function W is defined as W = (1-NRMSD(U,D)/2), where NRMSD(U,D) =2 * RMS(U-D)/[RMS(U)+RMS(D)]. [0008] The seismic data may be gathered using at least one source on the water surface and at least one receiver on the seabed. 30 [0009] The method may further include the step of aligning U 1 , D 1 , U 2 and D 2 with each other and/or aligning U 1 with U 2 and aligning D 1 with D 2 and may also include the step of aligning the up-going 4D difference dataset U with the down-going 4D difference dataset D. 2 WO 2013/101518 PCT/US2012/070079 [0010] In some embodiments, the method includes the step of scaling the data in at least two datasets so as to equalize the energy in each of a plurality of gates. The weighting function may also include a gate and the gate that is used to scale the data may be larger than the gate used in the weighting function. In particular, the weighting function may be an RMS 5 calculation that includes a gate, which may be less than 2 wavelengths, and the gate used to scale the data may be larger than the RMS gate. Thus, further, the RMS gate may selected to be less than 10 m. [0011] The time interval over which the present method is applied may be at least one month or at least one year. 10 [0012] According to other embodiments, a method for generating an image of a subsurface feature, comprises the steps of a) providing seismic data containing information about the feature, comprising a primary wave-field and a secondary wave-field (A 1 , B 1 ), b) creating a weighting function W that is a function of the similarity of the primary wave-field A 1 and the secondary wave-field B 1 ; and c) creating an image of the feature by generating a similarity 15 stack, where the similarity stack is defined as (A 1
    +B
    1 )*W/2. The method may further include the step of repeating steps a) - c) after a time interval and calculating a time-lapse dataset by comparing the resulting two similarity stacks. The method may further include the step of scaling the primary wave-field and the secondary wave-field so as to equalize the energy in each of a plurality of gates. In some embodiments, the weighting function may include a gate 20 and the gate used to scale the data may be larger than the gate used in the weighting function. The weighting function may be an RMS calculation that includes a RMS gate that is less than 2 wavelengths, wherein the gate used to scale the data may be larger than the RMS gate. [0013] The primary wave-field may be time-aligned with the secondary wave-field. [0014] The seismic data may be VSP data gathered using at least one source on the earth's 25 surface and at least one receiver in a borehole or the seismic data may be cross-well data gathered using at least one source in a first borehole and at least one receiver in a second borehole. Alternatively, the seismic data may be gathered using at least one source and at least two receivers that are located at different depths and the seismic data at each receiver may contain an up-going wave-field and a down-going wave-field comprising a receiver ghost. 30 [0015] In still other embodiments, the seismic data may be gathered using at least two receiver types, including but not limited to pressure sensors, strain sensors, velocity sensors, and acceleration sensors. 3 WO 2013/101518 PCT/US2012/070079 BRIEF DESCRIPTION OF THE DRAWINGS [0016] For a more detailed understanding of the invention, reference is made to the accompanying wherein: Figure 1 is a schematic illustration of a marine seismic system in which the present 5 invention can be used; Figure 2 is a schematic illustration of a VSP seismic system in which the present invention can be used; and Figure 3 is a schematic illustration of a cross-well seismic system in which the present invention can be used. 10 [0017] While some preferred embodiments of the invention are described below with respect to a marine seismic system, it will be understood that the invention is applicable in any seismic system. DETAILED DESCRIPTION OF A PREFERRED EMBODIMENTS 15 [0018] Referring initially to Figure 1, a marine seismic system 10 in which the present invention is advantageous includes at least one acoustic source 12 on the water's surface 14 and at least one sensing node 16 on or near the seafloor 18. A subsurface feature of interest, illustrated as horizon 22, lies in the subsurface 20 beneath seafloor 18. [0019] In system 10, there are multiple possible acoustic paths between each source 12 and 20 sensing node 16. In particular, some paths will result in an up-going wave arriving at node 16, as illustrated at 24, and some paths will result in a down-going wave arriving at node 16, as illustrated at 26. Each down-going wave-field has an additional reflection 28 at the air water interface where the reflection coefficient is close to -1.0. The down-going wave-field has two significant advantages over the up-going wave-field, namely, the down-going wave 25 field can image the near-seafloor environment and the down-going wave-field has a larger illumination area, resulting in increased fold at depth. [0020] Data for either the up-going or down-going wave-field will include information about both horizon 22 (and/or other subsurface features) and noise present at the time of sensing. This is evident in Reverse-Time Migrated (RTM) depth images of the baseline seismic data. 30 In preferred embodiments, the up-going data are time-aligned and/or depth-aligned to the down-going image so as to remove residual time shifts and are scaled to equalize the rms 4 WO 2013/101518 PCT/US2012/070079 energy, as is known in the art. With no further processing, the up-going and down-going images appear very similar. [0021] Using RTM data sets collected at different points in time, time-lapse (4D) differences can be generated by subtracting one data set from another and typically this is done after time 5 alignment of the baseline and monitor seismic data. 4D differences can be generated using up- or down-going migrations. When 4D differences for up- and down-going wave-fields are aligned and compared, a number of 4D signals stand out well above the noise levels. If the afore-mentioned scaling to match the energy levels with the higher-fold down-going wave field has been performed, the shallow portions of the up-going wave-field will contain more 10 noise than the deeper portions. [0022] Regardless, both the up- and down-going 4D difference sections will contain both 4D signal and 4D-noise. The 4D signals result from changes in reservoir properties over time and should be nearly identical on the up-going and down-going data, with only minor differences expected due to changes in offset distributions and raypath-transmission-imaging issues. In 15 contrast, the 4D-noise is expected to be completely different for the up- and down-going wave-fields. In other words, although noises from identical sources (e.g. multiples, swell noise, etc.) are present in both wave-fields, they will arrive at different times and positions. [0023] One way to suppress 4D-noise is to simply stack the up- and down-going 4D differences. Stacking is a good way to reduce noise and for purely random noise this generates 20 a statistical improvement of (1/1V). Visual inspection of the resulting image reveals that a stack of the 4D-differences looks somewhat less noisy than either of the inputs in the deeper portions, but is noisier in the shallow zones as compared with the down-going wave-field alone because of the noisier up-going data. [0024] According to preferred embodiments of the invention, an alternative method of 25 combining these datasets entails creating a weighted stack of the up- and down-going 4D differences using a weighting function based on the similarity of the two input datasets. By choosing the weighting function to be near 0.0 when the inputs are dissimilar and near 1.0 when they are similar, the 4D-noise can be muted. [0025] According to some preferred embodiments, a good weighting function can be 30 constructed based on the normalized root mean square (RMS) difference (NRMSD) of the 4D-difference data in a short gate that is preferable 0.5, 1, 1.5, or 2 wavelengths long. In one 5 WO 2013/101518 PCT/US2012/070079 preferred embodiment, described below, this gate is less than 200 ft, less than 100 ft, or in some instances less than 10 ft. [0026] In preferred embodiments where U and D represent the time-lapsed differences for the up- and down-going wave-fields, respectively, the NRMSD difference is: 5 NRMSD(U,D) = 2.* RMS(U-D)/[RMS(U)+RMS(D)]. Possible NRMSD values range from 0 (identical) to 2 (identical but with opposite sign) and can be used to generate a weighting function W = (1-NRMSD(U,D)/2). [0027] While NRMSD is one preferred type of weighting function, it will be understood that any suitable weighting function W that is a function of the similarity of the two datasets can 10 be used, including but not limited to the aforementioned W = (1-NRMSD(U,D)/2). Thus, W can alternatively be defined as (1-NRMSD(U,D)/NRMSD_cutoff , in which values of NRMS 2 above a predetermined cutoff value are set at zero, or can be based on NRMSD , or on an average difference of absolute values, manual weighting, and/or powers of any of the foregoing. 15 [0028] In embodiments where W 1 - NRMSD/NRMSDcuttoff, the cutoff value is set at or below 1 instead of allowing values between 1 and 2. This has the effect of suppressing noise more effectively, but at the expense of portions of the signal. In these embodiments, values of W less than zero are set to zero. [0029] It will be understood that scaling and/or time alignment is preferably performed on 20 each dataset as needed, preferably before calculating differences, stacking or other combinations of datasets. As is known in the art, scaling and alignment can be performed using matching filters. [0030] According to preferred embodiments, a "4D similarity stack" is defined as (U+D)*W/2. When the data are noisy, the weighting function W will be near 0 and the 4D 25 similarity stack will be small, effectively suppressing the 4D noises. In areas with true 4D signals, W will be near 1 and the signals will be preserved. An image of the 4D similarity stack derived using a 200 foot RMS gate showed a significant suppression of the 4D noise as compared to the non-weighted stack. [0031] To quantify the reduction in noise in the weighted stack, the NRMSD function can be 30 applied to the baseline and monitor data. By way of illustration only, in order to calculate this for the 4D similarity stack, an actual early down-going wave-field was used as the baseline data and a synthetic late (monitor) dataset was generated by subtracting the 4D similarity 6 WO 2013/101518 PCT/US2012/070079 stack from the baseline. Comparisons of NRMSD(baseline, monitor) as calculated for the up going, down-going and 4D similarity stack data in a 3000 foot gate centered at 15000 ft showed that NRMSD values are greatly reduced on the 4D similarity stack results and some of the underlying 4D signals stand out more from the noise. 5 [0032] As described, the similarity stack is effective at suppressing noise while preserving the desired 4D signal. In the following paragraphs, we compare theoretically what we expect in areas with and without 4D signal. Case 1: Regions with only random noise. [0033] In regions with only random noise, we expect NRMSD(U,D) -I2 and W - 0.3. The 10 4D similarity stack in this zone approximately equals 0.3*(U+D)/2, which reduces noise by a factor of 3 as compared with the conventional stack (U+D)/2. Case 2: Regions with a false 4D signal on either the up- or down-going data. [0034] In regions with large noise on one volume but not the other, NRMSD(U,D) - 2 and W 0. In these regions, the 4D similarity stack does an excellent job of removing false 4D 15 events. Case 3: Regions with both 4D signal and 4D noise. [0035] To simplify calculations for this case, we assume that 1) the level of 4D noise on the up and down going waves are comparable, and 2) 4D noise has zero correlation between up and down-going waves or with the 4D signal. With these assumptions, 1 12+2 *SNR 2 20 where SNR = RMS(4D-signal)/RMS(4D-noise) is the signal-to-noise ratio. For large SNR, W will be close to 1, but in noisy areas the weighting will also suppress the 4D signal. When SNR = 10, W = 0.93 and most of the 4D signal is preserved. In contrast, when SNR = 1, W = 0.5, so that in very noisy areas the 4D similarity stack will partially mute the signal. [0036] The method and examples described herein show that better time-lapse seismic data 25 can be obtained in deepwater environments where seismic recordings are made using ocean bottom sensors such as nodes (OBN) or cables (OBC). Specifically, it is possible to produce two independent images of the subsurface using up- and down-going wave-fields. The 4D signals should be nearly identical on these wave fields but the 4D noise differs greatly. A weighted stack of the up-going and down-going 4D differences based on their similarity will 30 suppress the 4D noises while preserving the signal. 7 WO 2013/101518 PCT/US2012/070079 [0037] More generally, whenever two independent images are acquired simultaneously, it can be assumed that the signal in both images will be the same but the noise in both images will be different. For deep sea applications these two images could refer to an up-going and a down-going wave-field, as discussed above, but the two independent-but-simultaneous images 5 could also refer to two images obtained in an acquisition with two (nearly) collocated detectors, for instance a pressure sensor and a velocity sensor, with two nearby detectors in an over/under configuration, or with two towed receiver streamers at different depths. Thus, the present concepts can be applied both offshore and onshore. While the foregoing description refers in some instances to up-going and down-going wave-fields, it will be understood that 10 any combination of such independent images could be used. [0038] In particular, the present invention may be applied to images generated with one or more seismic sources 12 on the surface and one or more receivers 16 in a borehole 30, as in the case of Vertical Seismic Profiling (VSP), illustrated schematically in Figure 2. The conventional reflection image is generated by migrating the up-going wave-field that arrives 15 at the receivers from reflections below. An image may also be generated using the down going wave-field by the technique of mirror or multiple migration. [0039] In still other embodiments, illustrated schematically in Figure 3, the present invention may be applied to images generated with one or more seismic sources 12 in a first borehole 30 and receivers 16 in another, nearby borehole 32, as in the case of cross-well seismic 20 surveying. For a given reflecting interface, one may generate an image using the up-going wave-field 24 which requires the seismic sources and receivers to be located above the interface in question. One may also generate an image using the down-going wave-field 26 which requires seismic sources and receivers to be located below the interface in question. Both situations occur in practice, as seismic sources and receivers are commonly placed above 25 and below the area of interest. [0040] In another embodiment, data may be collected at two or more (n) receivers at different depths or with different wave-fields. One or more weighting functions W can be derived for the aggregated signals belonging to the two or more receivers, respectively and these data can be summed using a multi-receiver (MR) similarity stack that is defined as MR = (R 1
    W
    1 + 30 R 2
    W
    2 + ... Rn Wn)/n, where Rn are the receiver data and the weighting functions Wn can be defined as functions of the similarity between each Rn and the other receivers. Two MR 8 WO 2013/101518 PCT/US2012/070079 similarity stacks based on data collected at different times can then be used to generated a noise-reduced 4D image. [0041] In another embodiment, data may be collected at two or more onshore receivers at different depths. In these instances, it may be desirable to enhance the up-going reflected 5 waves and to attenuate the down-going waves that are reflected from the free-surface, known as receiver ghosts. This can be done after aligning the up-going data, by calculating a weighting function based on the similarity of the two receivers' data and generating a weighted stack as described above. For shallow receivers, the time delay between up-going and down-going waves is may be less than the seismic wavelet length. In addition, seasonal 10 variations in near surface properties can cause seasonal variations in the receiver ghost, affecting the time lapse signal of the reservoir. Accordingly, enhanced suppression of the down-going waves, arriving shortly after the up-going waves, will improve the time lapse signal. [0042] In other variations, the seismic data may comprise VSP data that has been gathered 15 using at least one source on the earth's surface and at least one receiver in a borehole, or the seismic data may be cross-well data that has been gathered using at least one source in a first borehole and at least one receiver in a second borehole. Alternatively, the seismic data may be gathered using at least one source and at least two receivers that are located at different depths and the seismic data at each receiver may contain an up-going wave-field and a down-going 20 wave-field comprising a receiver ghost. In still other variations, the seismic data may be gathered using at least two receiver types, including but not limited to pressure sensors, strain sensors, velocity sensors, and acceleration sensors. [0043] It will be understood that other variations on the invention described above can be made without departing from the scope of the invention, which is set out in the claims that 25 follow. It will be understood that the sequential recitation of steps in the claims is not intended to require that the steps be performed sequentially unless expressly recited or logically required. 9
    权利要求:
    Claims (16)
    [1] 2. The method according to clairn I wherein the weighting function W is selected from 20 functions of the group consisting of NRMSD, (1-NRMSD(U,D)/NRMSD_cutoff, NRMSD 2 , average difference of absolute values, manual weighting, and powers of the foregoing functions.
    [2] 3. The method according to claim 2 wherein the weighting function W is defined as 25 W = (1-NRMSD(U,D)/2), where NRMSD(U,D) -2 * RMS(U-D)/[RMS(U)+RMS(D).
    [3] 4- The method according to claim I wherein the seismic data is gathered using at least one source on the water surface and at least one receiver on the seabed. 30 10 WO 2013/101518 PCT/US2012/070079
    [4] 5. The method according to claim 1, further including the step of aligning U 1 , DI, U 2 and D 2 (with each other.
    [5] 6. The method according to claim 1, further including the step of aligning U, with U 2 and 5 aligning D, with D2.
    [6] 7. The method according to claim 1, further including the step of aligning the up-going 4D difference dataset U with the down-going 4D difference dataset D. 10 8. The method according to claim 1, further including the step of scaling the data in at least two datasets so as to equalize the energy in each of a plurality of gates.
    [7] 9. The method according to claim 8 wherein the weighting function includes a gate and wherein the gate used to scale the data is larger than the gate used in the weighting function. 15
    [8] 10. The method according to claim 8 wherein the weighting function is an RMS calculation that includes a gate, wherein the RMS gate is less than 2 wavelengths, and wherein the gate used to scale the data is larger than the RMS gate. 20 I. The method according to claim 10 wherein the RMS gate is selected to be less than 10 M.
    [9] 12. The method according to claim I wherein the time interval is at least one month. 25 13. The method according to claim 1 wherein the time interval is at least one year.
    [10] 14. A method for generating an image of a subsurface nature, comprising the steps of: a) providing seismic data containing information about the feature, comprising a primary wave-field and a secondary wave-field (A,, BI); 30 b) creating a weighting function W that is a function of the similarity of the primary wave-field A, and the secondary wave-field BI; and c) creating an image of the feature by generating a similarity stack, where the similarity stack is defined as (A, +Bl)*Wt2. 11 WO 2013/101518 PCT/US2012/070079
    [11] 15. The method according to claim 14, further including the step of repeating steps a) - c) after a time interval and calculating a time-lapse dataset by comparing the resulting two similarity stacks. 5 16. The method according to claim 14 further including the step of scaling the primary wave-field and the secondary wave-field so as to equalize the energy in each of a plurality of gates.
    [12] 17. The method according to claim 16 wherein the weighting function includes a gate and 10 wherein the gate used to scale the data is larger than the gate used in the weighting function.
    [13] 18. The method according to claim 16 wherein the weighting function is an RMS calculation that includes a gate, wherein the RMS gate is less than 2 wavelengths, and wherein the gate used to scale the data is larger than the RMS gate. 15
    [14] 19. The method according to claim 14 further including the step of time-aligning the primary wave-field and the secondary wave-field.
    [15] 20. The method according to claim 14 wherein the seismic data is selected from the 20 groups consisting of VSP data gathered using at least one source on the earth's surface and at least one receiver in a borehole, cross-well data gathered using at )east one source in a first borehole and at least one receiver in a second borehole, data gathered using at least one source and at least two receivers that are located at different depths. 25 21. The method according to claim 14 wherein the seismic data at each receiver contains an up-going wave-field and a down-going wave-field comprising a receiver ghost.
    [16] 22. The method according to claim 14 wherein the seismic data is gathered using at least two receiver types, selected from the group consisting of pressure sensors, strain sensors, 30 velocity sensors, and acceleration sensors. 12
    类似技术:
    公开号 | 公开日 | 专利标题
    US20060190181A1|2006-08-24|Method for combining seismic data sets
    Yu et al.2017|Attenuation of noise and simultaneous source interference using wavelet denoising
    Moldoveanu et al.2012|Broadband seismic methods for towed-streamer acquisition
    Abbad et al.2009|Automatic nonhyperbolic velocity analysis
    US20120215453A1|2012-08-23|Device and method for multi-dimensional coherency driven denoising data
    AU2012362876B2|2015-01-15|Suppressing 4D-noise by weighted stacking of simultaneously acquired wave-fields
    Ravasi et al.2015|Vector-acoustic reverse time migration of Volve ocean-bottom cable data set without up/down decomposed wavefields
    Al-Hagan et al.2014|Iterative supervirtual refraction interferometry
    Majdański et al.2011|Attenuation of free-surface multiples by up/down deconvolution for marine towed-streamer data
    Jones2013|Tutorial: The seismic response to strong vertical velocity change
    Zhao et al.2018|Wavelet-crosscorrelation-based interferometric redatuming in 4D seismic
    Lester et al.2012|Multiple attenuation in crustal-scale imaging: examples from the TAIGER marine reflection data set
    Górszczyk et al.2019|Crustal-scale depth imaging via joint full-waveform inversion of ocean-bottom seismometer data and pre-stack depth migration of multichannel seismic data: a case study from the eastern Nankai Trough
    Staring et al.2021|Robust estimation of primaries by sparse inversion and Marchenko equation-based workflow for multiple suppression in the case of a shallow water layer and a complex overburden: A 2D case study in the Arabian Gulf
    Biondi et al.2014|Nonstretch normal moveout through iterative partial correction and deconvolution
    Dondurur et al.2012|Swell noise suppression by Wiener prediction filter
    GB2499883A|2013-09-04|Estimating time-shifts per event and per trace in seismic data
    Hatchell et al.2012|Suppressing 4D-noise by weighted stacking of up-going and down-going migrations
    AU2011248910B2|2014-10-16|System and method for compensating time and offset varying near-surface effects in seismic data background
    Burnstad et al.2012|Successful imaging of land hydrophone and dual sensor data in a dry desert environment
    Hatchell et al.2012|Suppressing 4D-noise by weighted stacking of up-going and down-going wave-fields
    Al-Anboori et al.2005|Approximate separation of pure-mode and converted waves in 3-C reflection seismics by τ-p transform
    Zhao et al.2017|Wavelet crosscorrelation and its application to virtual source denoising
    Whitmore2016|Extending the 3D primary image with multiples and mirrors
    Van Renterghem et al.2017|Land and ocean-bottom spatial gradient-based seismic wavefield separation
    同族专利:
    公开号 | 公开日
    CA2862250C|2021-03-02|
    GB2510785B|2017-07-12|
    BR112014015891A8|2021-03-16|
    CN104024890A|2014-09-03|
    WO2013101518A1|2013-07-04|
    GB2510785A|2014-08-13|
    CA3099540A1|2013-07-04|
    AU2012362876B2|2015-01-15|
    CA2862250A1|2013-07-04|
    US20150168575A1|2015-06-18|
    US9581709B2|2017-02-28|
    GB201410479D0|2014-07-30|
    CN104024890B|2016-08-31|
    引用文献:
    公开号 | 申请日 | 公开日 | 申请人 | 专利标题
    US5638338A|1995-10-16|1997-06-10|Peterson; Fred M.|Seismic processing apparatus and method|
    AT339699T|2000-12-07|2006-10-15|Westerngeco Seismic Holdings|METHOD AND DEVICE FOR PROCESSING SEISMATIC DATA|
    GB0222524D0|2002-09-27|2002-11-06|Westerngeco Seismic Holdings|Calibrating a seismic sensor|
    GB2420408B|2004-11-19|2008-03-12|Geophysique Cie Gle|Method for processing at least two sets of seismic data|
    US7952960B2|2006-10-03|2011-05-31|Bp Corporation North America Inc.|Seismic imaging with natural Green's functions derived from VSP data|
    US8259531B2|2008-03-31|2012-09-04|Statoil Asa|Method for reflection time shift matching a first and a second set of seismic reflection data|
    US8339898B2|2008-05-25|2012-12-25|Westerngeco L.L.C.|Processing seismic data using combined regularization and 4D binning|EP3039459A2|2013-08-27|2016-07-06|CGG Services SA|4d noise suppression|
    EP2846175B1|2013-09-06|2022-02-23|Services Pétroliers Schlumberger|Seismic survey analysis|
    EP3101450B1|2015-06-04|2021-04-14|SpotLight|Quick 4d detection seismic survey|
    CN105549083B|2015-12-25|2018-03-13|中国石油天然气集团公司|A kind of method and apparatus for determining seismic observation system repeatability and integrally measuring|
    法律状态:
    2015-05-14| FGA| Letters patent sealed or granted (standard patent)|
    优先权:
    申请号 | 申请日 | 专利标题
    US201161580693P| true| 2011-12-28|2011-12-28||
    US61/580,693||2011-12-28||
    PCT/US2012/070079|WO2013101518A1|2011-12-28|2012-12-17|Suppressing 4d-noise by weighted stacking of simultaneously acquired wave-fields|
    [返回顶部]