Multi-shooting approach to seismic modeling and acquisition

• US 6,327,537 B1
• Filed: 07/19/1999
• Issued: 12/04/2001
• Est. Priority Date: 07/19/1999
• Status: Expired due to Fees

• Current Assignee: Luc T. Ikelle
• Original Assignee: Luc T. Ikelle
• Inventors: Ikelle, Luc T.
• Primary Examiner(s): McElheny, Jr., Donald E.
• Application Number:
  US09/356,571
• Time in Patent Office: 869 Days
• Field of Search: 702/14, 702/16, 702/17, 702/18, 367/73
• US Class Current: 702/14
• CPC Class Codes: G01V 1/37 specially adapted for seismic systems using continuous agitation of the ground , e.g. using pulse compression of frequency swept signals for enhancement of received signals
• CPC Class Codes: G01V 11/00 Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Abstract

A multi-shooting approach to seismic modeling and acquisition where several shot gathers can be generated simultaneously. The method is called a multi-shooting approach to seismic modeling and acquisition. A multi-shooting modeling method is disclosed which may be carried out on a computer system in either an explicit or an implicit manner. A multi-shooting acquisition method is also disclosed. The finite difference technique has the ability to generate elastic waves from several locations and at different time intervals simultaneously. The multi-shooting approach is based on the property of the finite difference technique. The decoding solution exploits the fact that specific time delays can be imposed between the various shots and that two closely space shooting points produce nearly identical responses.

Claims

1. A method for simulating a multi-shooting technique where a plurality of shots are generated simultaneously and recorded in a single shot gather, said method comprising the steps:

(a) let x denote a generic point in a three dimensional medium;

(b) let x_s⁽ⁱ⁾ denote a generic source position;

(c) let F_P(x,t) denote a field matrix of the wave motion at position x and time t;

(d) solve the equation $\left(D_{I,P}<br><br>\left(x,t\right)+M_{I,P}<br><br>\left(x\right)\right)<br><br>F_P<br><br>\left(x,t\right)=\underset{i=1}{\overset{i=n}{\sum <br><br>}}<br><br>\text{\hspace{1em}<br><br>}<br><br>Q_I^{\left(i\right)}<br><br>\left(x,t-{\mathrm{\tau <br><br>}}_i,x_s^{\left(i\right)}\right)$

wherein D_I,P(x,t) is a block off-diagonal spatial differentiation operator matrix which is specific for each type of wave motion under consideration, wherein M_I,P(x) is a medium matrix that is representative for the properties of the media in which the wave propagates, and wherein Q_I(x,t) is a volume source density matrix that is representative for the action of the volume sources that generate the wavefield, wherein I and P denote the pertaining matrix elements and the summation convention for repeated subscripts applies, where n shots are generated simultaneously and recorded in a single shot gather F_P(x,t), and wherein τ

_i describes time delays between various shots; and

(e) decoding the multi-shooting data into single shots.

2. The method for simulating a multi-shooting technique according to claim 1, wherein said decoding step is carried out utilizing a seismic time division multiple access technique.

3. The method for simulating a multi-shooting technique according to claim 1, wherein said decoding step is carried out utilizing a code division multiple access technique.

4. A method for simulating a multi-shooting technique where a plurality of shots are generated simultaneously and recorded in a single shot gather, said method comprising the steps:

(a) let x denote a generic point in a three dimensional medium;

(b) let x_s⁽ⁱ⁾ denote a generic source position;

(c) let F_P(x,t) denote a field matrix of the wave motion at position x and time t;

(d) solve the equation $\left(D_{I,P}<br><br>\left(x,t\right)+M_{I,P}<br><br>\left(x\right)\right)<br><br>F_P<br><br>\left(x,t\right)=\underset{i=1}{\overset{i=n}{\sum <br><br>}}<br><br>\text{\hspace{1em}<br><br>}<br><br>Q_I^{\left(i\right)}<br><br>\left(x,t-{\mathrm{\tau <br><br>}}_i,x_s^{\left(i\right)}\right)$

wherein D_I,P(x,t) is a block off-diagonal spatial differentiation operator matrix which is specific for each type of wave motion under consideration, wherein M_I,P(x) is a medium matrix that is representative for the properties of the media in which the wave propagates, and wherein Q_I(x,t) is a volume source density matrix that is representative for the action of the volume sources that generate the wavefield, wherein I and P denote the pertaining matrix elements and the summation convention for repeated subscripts applies, where n shots are generated simultaneously and recorded in a single shot gather F_P(x,t), and wherein τ

_i describes time delays between various shots;

(e) positioning a plurality of remote recording units for receiving and recording shooting signals as a function of time during a recording period;

(f) positioning a plurality of shooting sources at different physical locations;

(g) triggering each shooting source during the recording period; and

,

(h) decoding the multi-shooting data into single shots.

5. The method for simulating a multi-shooting technique according to claim 4, wherein said decoding step is carried out utilizing a seismic time division multiple access technique.

6. The method for simulating a multi-shooting technique according to claim 4, wherein said decoding step is carried out utilizing a code division multiple access technique.