Transmission without reverberation by iterative incomplete time-reversal

Patent metadata

Field	Value
Publication number	US-2016334270-A1
Application number	US-201515111739-A
Country	US
Kind code	A1
Filing date	Jan 13, 2015
Priority date	Jan 14, 2014
Publication date	Nov 17, 2016
Grant date	—

Abstract

For wavefield analysis/processes, constructing f1 solutions and utilizing the f1 solutions of a medium using iterative incomplete time-reversal. Using sources convolved with the f1 solution, propagating a wavefield without internal reverberation in a medium. Many applications of the resulting reverberation-free wavefield or virtual data are also disclosed.

First claim

1 . A method for generating waves propagating into a medium bounded by an upper level and a bottom level, wherein at least one source is at or above an upper level and at least one receiver is at the upper level, wherein reflected responses by the medium are received by the at least one receiver, the method, comprising:
(a) obtaining an f1 solution for the medium; and (b) activating the at least one source according to the f1 solution to generate waves wherein the waves propagate into the medium, wherein the waves focus at a focal point without reverberations.

2 . The method of claim 1 , wherein obtaining an f1 solution of the medium comprises:
(a) injecting a single downgoing pulse by the at least one source into subsurface and recording the reflection responses by the at least one receiver; (b) time-reversing the recorded reflection response and re-injecting it back into the medium by the at least one source, and recording its reflection by the at least one receiver; (c) muting the recorded reflected time-reversed reflection response just before it re-focuses on the source pulse; and (d) time-reversing the muted reflected time-reversed reflection response and adding it just after the single downgoing pulse forming a new downgoing pulse.

3 . The method of claim 2 , wherein obtaining an f1 solution of the medium further comprises:
(e) repeating steps (a)-(d) until convergence or a predefined stopping criterion is reached, wherein in repeating the steps in step (a) of the repeated step a previous downgoing pulse is replaced with a new downgoing pulse as obtained in step (d).

4 . The method of claim 1 , wherein obtaining an f1 solution of the medium comprises:
performing by a computer an iterative incomplete time-reverse f1 solution construction comprising:
(a) injecting a single downgoing pulse into subsurface and recording the reflection responses;
(b) time-reversing the recorded reflection response and re-injecting it back into the medium, and recording its reflection;
(c) muting the recorded reflected time-reversed reflection response just before it re-focuses on the source pulse; and
(d) time-reversing the muted reflected time-reversed reflection response and adding it just after the single downgoing pulse forming a new downgoing pulse.

5 . The method of claim 4 , wherein obtaining an f1 solution of the medium further comprises:
(c) repeating steps (a)-(d) until convergence or a predefined stopping criterion is attained, wherein in each iteration the previous downgoing pulse is replaced with the new downgoing pulse as obtained in step (d).

6 . The method as in claim 1 , wherein the waves focused at the single focal point cause material changes to the medium at the focal point.

7 . The method as in claim 6 , wherein the medium at the focal point is biological tissues or earth formations.

8 . The method as in claim 1 , wherein the waves generated by the at least one source are symbol carriers.

9 . The method as in claim 1 , wherein the waves generated by the at least one source are seismic waves; and wherein the waves propagating within the medium have substantially no internal multiple reflections.

10 . A computer implemented method for processing measurements of reflected responses of waves propagating into a medium bounded by an upper level and a bottom level, wherein at least one source is at or above a upper level and at least one receiver is at the upper level, wherein reflected responses by the medium are received by the at least one receiver, the method, performed by a computer, comprising:
(a) obtaining an f1 solution of the medium; and (b) activating the at least one source according to the f1 solution to generate waves wherein the waves propagate into the medium, wherein the waves focus at a focal point without reverberations.

11 . The method as in claim 10 , wherein obtaining an f1 solution of the medium comprises:
(a) injecting a single downgoing pulse into subsurface and recording the reflection responses; (b) time-reversing the recorded reflection response and re-injecting it back into the medium, and recording its reflection; (c) muting the recorded reflected time-reversed reflection response just before it re-focuses on the source pulse; and (d) time-reversing the muted reflected time-reversed reflection response and adding it just after the single downgoing pulse forming a new downgoing pulse.

12 . The method of claim 11 , wherein obtaining an f1 solution of the medium further comprises:
(e) repeating steps (a)-(d) until convergence or a predefined stopping criterion, each time replacing the previous downgoing pulse with the new downgoing pulse obtained in step (d).

13 . The method as in claim 11 , wherein a subsurface object has a focal time t f , which is the time for a wave traveling from the upper level to a focal depth x f ,
wherein (a) injecting a single downgoing pulse into the subsurface is injected at time t=−t f ; wherein (c) muting the recorded reflected time-reversed reflection response just before it re-focuses on the source pulse is muted at time t=−t f ; wherein (d) time-reversing the muted reflected time-reversed reflection response and adding it just after the single downgoing pulse is at time t=−t f ; and further comprising muting the reflection response for t>t f , when primary reflections from all the interfaces above a desired focus have arrived.

14 . The method as in claim 13 , wherein the focal depth is within or below the object.

15 . The method as in claim 13 , wherein the object is in a 1D, 2D or 3D system.

16 . The method as in claim 13 ,
wherein the object is a geophysical target ( 65 ), a biological target ( 71 ), a remote-sensing target or a non-invasive investigation target.

17 . The method as in claim 10 , further comprising:
from the wavefields, retrieving true amplitude full-wave responses by receivers at desired locations within the medium or below the medium, wherein the retrieved true amplitude full wave responses are virtual data.

18 . The method as in claim 17 ,
wherein the desired locations are on an enclosing boundary of an object, the method further comprising: applying a seismic migration imaging process on the virtual data accounting for both transmission and reflection data; or applying a waveform inversion process on the virtual data accounting for both transmission and reflection data.

19 . The method as in claim 18 , the method further comprising:
reconstructing extended-image gathers with the virtual data accounting for both transmission and reflection data; or applying a migration velocity analysis or image domain waveform inversion with the virtual data accounting for both transmission and reflection data.

20 . The method as in claim 19 , the method further comprising:
selecting a subset of the object and the corresponding virtual data; and applying an imaging or inversion process for high-resolution targeted characterization on the selected subset.

21 . The method as in claim 20 , wherein second measurements of reflected responses recorded by receivers at a second time, the method further comprising:
processing the second measurements with the method of claim 1 to obtain a second virtual data at the second time; selecting a subset of the object and the corresponding second virtual data; applying an imaging or inversion process for high-resolution targeted characterization on the selected subset for the second time; and comparing the two h

Assignees

• Westerngeco Llc

Inventors

• Van Manen Dirk-Jan
• Pires De Vasconcelos Ivan

Classifications

• G01V2210/679

  Reverse-time modeling or coalescence modelling, i.e. starting from receivers · CPC title

• A61B8/0866

  involving foetal diagnosis; pre-natal or peri-natal diagnosis of the baby · CPC title

• G01V1/364

  Seismic filtering (G01V1/37 takes precedence) · CPC title

• A61B8/5207

  involving processing of raw data to produce diagnostic data, e.g. for generating an image · CPC title

• G01H3/00 Primary

  Measuring {characteristics of} vibrations by using a detector in a fluid (G01H7/00, G01H9/00, G01H11/00 take precedence) · CPC title