1 / 73

Multiple Removal with Local Plane Waves

Multiple Removal with Local Plane Waves. Dmitri Lokshtanov. Content. Motivation WE multiple suppression operator Fast 2D/3D WE approach for simple sea-floor 2D/3D WE approach for irregular sea-floor Conclusions. Motivation. Seismic processing and imaging - main challenges:

kmcmullen
Download Presentation

Multiple Removal with Local Plane Waves

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Multiple Removal with Local Plane Waves Dmitri Lokshtanov

  2. Content • Motivation • WE multiple suppression operator • Fast 2D/3D WE approach for simple sea-floor • 2D/3D WE approach for irregular sea-floor • Conclusions

  3. Motivation • Seismic processing and imaging - main challenges: • Velocity model building for sub-salt and sub-basalt imaging • Removal of multiples from strong irregular boundaries

  4. Near offset section (no AGC)

  5. Depth migration with water velocity

  6. Input shot gathers (no AGC)

  7. Multiple suppression • For multiples from complex boundaries the methods based on periodicity or kinematic discrimination usually don’t work or are not sufficient. • In such cases the main demultiple tools are based on the Surface Related Multiple Elimination (SRME) or Wave-Equation (WE) techniques.

  8. SRME(Berkhout, 1982; Verschuur, 1991) – advantages and limitations • Does not require any structural information. Predicts all free-surface multiples • As a rule becomes less efficient with increased level of interference of multiples of different orders • Requires the same dense sampling between sources as between receivers • Noise in data and poor sampling significantly degrade the prediction quality • Missing traces required by 3D SRME are reconstructed with least-square Fourier or Radon interpolation; residual NMO correction; DMO/inverse DMO; migration/demigration

  9. WE approach versus SRME • SRME is the method of preference for data from areas with deep sea-floor, especially when a thick package of strong reflectors is present below the sea-floor • WE approach is especially efficient when the main free-surface multiples are just ‘pure’ water-layer multiples and peg-legs. Gives usually better results than SRME when several orders of multiples are involved • 3D WE approach has less sampling problems than 3D SRME and it gives a flexiblility in methods for wavefield extrapolation depending on complexity of structure

  10. The operator Pg transforms the primary reflection event recorded at receiver 1 into the multiple event recorded at receiver 2 (Wiggins, 1988; Berryhill & Kim, 1986). The operator Pg transforms the primary reflection event recorded at receiver 1 into the multiple event recorded at receiver 2 (Wiggins, 1988; Berryhill & Kim, 1986).

  11. Principles of WE approach

  12. Adaptive subtraction of predicted multiples

  13. Wave-equation approach – main features • All predicted multiples are split into 3 terms, where each term requires the same amplitude correction • All source-side and receiver-side multiples of all orders are suppressed simultaneously in one consistent step • The prediction and the adaptive subtraction of multiples are performed in the same domain • Fast version (WEREM) for a simple sea-floor. Slower version for irregular sea-floor

  14. Why in the tau-p domain • Easier to apply antialiasing protection • No problems with muting of direct arrival • Easier to define ‘multiple’ zone of tau-p domain and mute it away • Estimated reflection coefficients are explicitly angle dependent

  15. 2D WEREM – prediction of multiples - 1

  16. 2D WEREM – prediction of multiples - 2

  17. Velocity model used to generate synthetic FD data

  18. Constant P sections (angle at the surface is about 3º) Input After Werem After Remul

  19. Constant P sections (angle at the surface is about 15º) Input After Werem After Remul

  20. Velocity model 2 used to generate synthetic FD data

  21. Constant P sections (angle at the surface is about 3º)  m. residual Input After Werem

  22. T. Shetland  T. Draupne  T. Brent Stack before multiple suppression Stack after Werem

  23. Constant P sections (angle at the surface is about 10º Input After Werem

  24. Stack before multiple suppression (left) and after Werem multiple suppression (right). The pink line shows the expected position of the first-order water-layer peg-leg from the Top Cretaceous (black line). The multiple period is about 140 msec.

  25. Constant P sections (angle at the surface is about 8º Input After Werem multiple suppression Difference

  26. raw stack

  27. stack after WEREM

  28. 500 m input WEREM

  29. 4000 m input WEREM

  30. Improving the results - local prediction / subtraction of multiples • Within the same prediction term, for the same CMP and the same p we have events reflected at different positions along the water bottom • Inconsistency between prediction and subtraction in case of rapid variation of sea-floor reflectivity • The problem is partly solved by applying adaptive subtraction in different time windows • Or by making prediction dependent on both p and offset (window)

  31. 3D WEREM – basic features • 3D data can be represented as a sum of plane waves with different vertical angles and azimuths from the source-side and receiver-side. • Current quasi 3D marine acquisition does not allow full 4D decomposition • Decomposition uniquely defines the direction of propagation from the receiver-side and is an integral over crossline slownesses from the source-side • The result of decomposition are used for exact prediction of multiples from the receiver-side and approximate prediction from the source-side • The approximation is that the crossline slowness from the source-side is the same as from the receiver-side (the same azimuth for 1D structures). The approximation allows us to mix data for flip flop shooting

  32. Input constant P section (small angles)

  33. Predicted multiples – R-side (small angles)

  34. Input constant P section (small angles)

  35. Constant P section – after prediction / subtraction (small angles)

  36. Difference (Input – 3D WEREM), small angles

  37. Input constant P section (larger angles)

  38. Constant P section – after 3D WEREM (larger angles)

  39. Werem - conclusions • Very efficient when the main assumptions are met: strongest multiples are water-layer multiples and peg-legs and the sea-floor is simple • Very fast - each predicted p trace is simply obtained as a sum of time-delayed input traces with the same p from the neighbour CMPs

  40. WE for irregular sea-floor • Kinematic prediction of multiples (extrapolation through the water layer) takes into account coupling between incident and reflected / scattered plane waves with different slownesses • Both multiple reflections and diffractions are predicted • The procedure starts from the Radon transformed CS gathers (no interleaving is required) • In 3D exact prediction from the receiver side; approximate prediction from the source side

  41. 2D prediction of multiples from the receiver side for irregular sea-floor

  42. 2D prediction of multiples from the source side for irregular sea-floor

  43. Velocity model for FD modelling

More Related