SOURCE SIGNATURE DETERMINATION AND MULTIPLE REFLECTION REDUCTION

Abstract: The signature of an energy source is determined from an inverse scattering Born series representing multiple reflected energy. The series comprises a polynomial in the inverse of the signature and has recorded data as the first term. The series is truncated, preferably to the first two terms to permit an analytical determination of the signature to be found. The value of the inverse signature which minimises the energy represented by the sub-series is found and this corresponds to the desired source signature. An iterative scheme may be adopted to improve the match to the actual signature so as to take into account the errors caused by truncating the scattering series.

Description: /20417 PCI7GB95/02934 SOURCE SIGNATURE DETERMINATION AND MULTIPLE REFLECTION REDUCTION The present invention relates to a method of determining a signature of a source of energy and to a method of reducing d e effects of multiple reflected energy. Such a method may be used in medical imaging and non-destructive evaluation, where a sample is irradiated with energy, for instance in the form of x-rays or ultrasound, in order to determine the internal structure of the sample non-invasively. Such a method may also be used with marine seismic reflection data obtained by means of seismic sources and receivers located in water and towed behind a seismic exploration vessel. These methods may be used during actual surveying and/or subsequently with recorded data from such surveys. During seismic surveying, a seismic source is repeatedly actuated and seismic receivers, such as hydrophones in marine seismic surveying, receive energy direct from the sources and reflected from various boundaries or interfaces. In the case of marine seismic surveying, energy propagates into the earth and is reflected back to the hydrophones from subterranean boundaries or interfaces, for instance between strata of different types. The recorded seismic data can be processed to reveal information about the structure of the earth in the area being surveyed. However, such reflections are contaminated by other reflection paths. For instance, energy from the sources is reflected at the sea surface directly to the hydrophones. Also, energy can be reflected more than once between the sources and the receivers. Such multiple reflections can take place within the earth. Also, energy initially travelling downwards from the sources can be reflected upwardly and then downwardly again from the surface of the sea before arriving at the hydrophones. Reflections of this type are referred to as "free-surface multiple reflections". Free- surface multiple reflections can be classified according to their order, which is equal to the number of reflections from the free-surface. Thus, first order free-surface reflections comprise energy initially travelling downwardly from the sources (as opposed to "ghosting" where energy travels upwardly and is reflected from the free- surface), is reflected upwardly from the sea bed or a boundary below the sea bed, and is then reflected downwardly from the free-surface to the hydrophones. Second order free-surface multiple reflections undergo two downward reflections from the sea-surface before being detected by the hydrophones, and so on. In order to remove or reduce the effects of multiple reflections including free- surface multiple reflections in seismic data, a good knowledge of the seismic source signature is required. The seismic source signature is the wave shape e.g. the pressure waveform with respect to time, of the seismic energy emitted by the seismic sources. Removal or reduction of the effects of multiple reflections in seismic data is performed before "stacking", the known process of summing seismic data related to each sub-region of the region being explored so as to improve the effective signal-to-noise ratio of the data. Good knowledge of the seismic source signature also enhances other data processing techniques such as deconvolution, migration and inversion for elastic parameters. Thus, good knowledge of the seismic source signature can significandy enhance the final data processing products. There are known techniques for determining the signature of a seismic source. One such known technique relies on a "convolution model" for seismic data, in which the model is defined as the convolution of the source signature with the impulse response of the earth (including the effects of reflections, refractions, multiples, and diffractions). This technique requires known information about some part of the sub- surface, such as the sea bottom, so as to allow estimation of the source signature to be treated as a linear inverse problem. However, such sub- surface information is often not available and the technique must then rely on statistical methods which require more severe assumptions, for instance about the phase of the signature or the "whiteness" of the reflectivity sequence. Accordingly, the practical use of this technique and its accuracy are limited. Other known techniques comprise actually measuring the signature. For instance, in very deep water, the direct wave from the sources to the hydrophone can be measured without interference from reflections from the sea bed and sub-strata. "Ghosting" reflections from the sea surface can be subtracted with a sufficient degree of accuracy to allow the source signature to be directly determined. However, it is generally not convenient to perform such measurements as this generally requires towing the sources and hydrophones to regions of deep water and then returning to a region to be surveyed. Further, the source signature may vary with time. For instance, where the or each source comprises an array of individual sources, such as air guns or water guns, one or more of the individual sources may cease operating during actual surveying, in which case the signature of the whole source will change and will therefore not be known with sufficient accuracy. Other techniques for measuring the source signature while performing seismic surveys are also known. However, many of these techniques require special data acquisition geometries i.e. special arrangements of sources and/or hydrophones, which are often not available. Examples of such techniques are disclosed in EP-A- 0 066 423, in Landr , M. and Sollie, R., 1992, Source signature determination by inversion: Geophysics, 57, 1633-1640, and in Weglein, A.B. and Secrest, B.G., 1990, Wavelet estimation for a multidimensional acoustic or elastic earth: Geophysics, 55, 902-913. According to a first aspect of the invention, there is provided a method of determining a signature of an energy source as defined in the appended Claim 1. According to a second aspect of the invention, there is provided a method of reducing the effects of multiple reflected energy as defined in the appended Claim 11. Preferred embodiments of the invention are defined in the other appended claims. It is thus possible to provide a technique which may be applied to marine seismic reflection data to allow source estimation for "prestack" data processing. This technique allows the source signature which permits reduction or removal of events due to first order free-surface reflections to be found. This technique makes use of the formulation of the relationship between free-surface reflections and the source signature as a scattering Born series as disclosed in Carvalho, P.M., Weglein, A.B. and Stolt, R.H., 1991, Examples of a nonlinear inversion method based on the T matrix of scattering theory: Application to multiple suppression: Mtg. Soc. Expl. Geophys., Expanded Abstracts, 1319-1322. The series is constructed exclusively with seismic data and the source signature without any knowledge of the subsurface other than the velocity of sea water. Thus, by using the scattering Born series for removing free-surface multiples, the source signature may be determined without any of the assumptions usually associated with known techniques based on the classical convolution model. Further, special acquisition geometries of sources and receivers are not required. The source signature may be determined directly from data recorded during seismic surveying. Further, the source signature thus determined may be used to reduce or remove the effects of free-surface multiple reflections so as to provide a substantial improvement in the quality of the data before stacking and when other processing techniques are applied. SUDSTITUTE SHEIZT (HULE 26) By using only the first two terms of the inverse scattering series, which is physically equivalent to dealing only with first order free-surface reflections, the relationship between free-surface reflections and the source signature has effectively been made linear as opposed to the nonlinear relationship which applies when higher order free-surface reflections and higher order terms of the inverse scattering series are considered. The mathematical problem of finding the source signature which minimises the "energy" of the seismic data (corresponding to substantially removing the effects of first order free-surface reflections) has an analytical solution and, in many cases, this solution permits the source signature to be determined to a desired degree of precision. However, this precision is limited by truncating the inverse scattering series and, if desired, a better approximation to the source signature may be obtained iteratively, for instance by finding a series of corrections to the source signature which progressively reduce the energy of the seismic data. Such iterations thus correct at least partially for the effects of truncation of the inverse scattering series. This iterative process may be stopped when a sufficient degree of precision has been achieved, for instance when further reductions in energy are not considered significant or when further corrections to the source signature are not considered significant. The invention will be further described, by way of example, with reference to the accompanying drawings, in which: Figure 1 is a schematic diagram illustrating marine seismic surveying; Figure 2 is a graph illustrating synthetic seismic reflection data as time against offset showing the effect of free-surface multiple reflections; Figure 3 is a graph similar to that of Figure 2 illustrating removal of the first order free-surface multiple reflections; Figure 4 is a graph similar to that of Figure 2 illustrating removal of first and second order free-surface multiple reflections; Figure 5 is a graph similar to that of Figure 2 illustrating removal of first, second, and third order free-surface multiple reflections; Figure 6 is a graph of amplitude against time showing the amplitude spectrum of the actual source signature and the initial estimated source signature corresponding to the data shown in Figure 2; Figure 7 is a graph of amplitude against frequency showing the frequency spectrum of the actual source signature and the initial estimated source signature corresponding to the data shown in Figure 2; Figure 8 is a graph of phase angle against frequency showing the phase spectrum of the actual source signature and d e initial estimated source signature corresponding to the data shown in Figure 2; Figure 9 is a graph similar to that shown in Figure 2 illustrating the seismic data after removal of free-surface multiple reflections using the initial estimated source signature shown in Figure 6; and Figures 10 to 13 correspond to Figures 6 to 9, respectively, and illustrate the estimated source signature and the effects of multiple removal after five iterations to reduce d e energy of the seismic data. Figure 14 is a flow-chart illustrating an example in accordance with the invention. For the purpose of illustration, the use of the present invention with marine seismic reflection data will be described in detail. However, these techniques may be applied to situations where multiple reflected energy can be represented by an inverse scattering series. Figure 1 illustrates a typical arrangement for performing marine seismic exploration. The drawing shows a section through the earth 1 below the sea 2 with the sea bed shown at 3 and die sea surface shown at 4. The sea surface 4 constitutes an interface between the air and the water and thus constitutes a free- surface in terms of seismic reflection data. An exploration vessel 5 tows a seismic source 6, for instance comprising an array of air guns. The vessel 5 also tows a seismic "streamer" 7, which is cable several hundreds or thousands of metres in length carrying hydrophones regularly spaced along d e cable and connected to data recording apparatus on board the vessel 5. The spacing between adjacent hydrophones is typically of die order of one or a few tens of metres. During exploration, the vessel 5 tows the source 6 and streamer 7 at relatively low speed along parallel lines above the sea bed 3 while actuating the source 6 at regular intervals. The seismic signals arriving at the hydrophones of the streamer 7 are recorded on board the vessel with or without data processing. The recorded seismic data are generally subjected to further processing elsewhere so as to reveal information about the structure of the earth 1 below the exploration region. The energy emitted by me source 6 arrives at the hydrophones of die streamer 7 via various different types of propagation paths. The propagation paths which are required for processing of the seismic data to reveal information about the structure of the earth are downward paths such as 8 followed by upward paths such as 9 to a typical hydrophone 10 following a single reflection, for instance at a boundary between different substrata illustrated at 11. However, these reflection data are contaminated with data arriving via other propagation paths. For instance, there is a direct path from the source 6 to each hydrophone, for instance shown as the path 12 to the hydrophone 13. Also, some energy initially travels upwardly from the source 6, for instance shown by the path 14, to be reflected downwardly by the sea surface 4 as shown by the path 15. The energy reflected downwardly by the sea surface 4 also travels to the sea bed 3 and into the earth 1 to be reflected back to the hydrophones. This gives rise to a virtual source shown in broken lines at 16. This effect is known as "ghosting". In general, direct waves following propagation paths such as 12 and virtual sources 16 resulting from reflection at the sea surface 4 do not cause a problem as their effects can relatively easily be subtracted from the seismic signals received by the hydrophones, such as 10 and 13. The seismic signals are also contaminated by multiple reflections and other effects such as refractions and diffractions. Multiple reflections can occur wiύiin the earth 1. However, the present invention is concerned with free-surface multiple reflections i.e. energy from the source 6 which initially travels downwardly but is received by the hydrophones after at least one downward reflection from the sea surface 4. The propagation of a typical free-surface multiple reflection is shown in Figure 1. Energy from the source 6 travels initially downwardly along a path 17. Part of the energy is reflected by the sea bed 3 to follow the path 18. This energy is substantially wholly reflected by the sea surface 4 along a path 19 and then again partially reflected upwardly by the sea bed 3 along a path 20 to arrive at the SUBSTITUTE SHEET (RULE 2G) hydrophone 10. Such free-surface multiple reflections are categorised by die number of downward reflections which take place at the free-surface or sea surface 4 following the initial downward propagation of die energy from the source 6. Thus, a first order free-surface multiple reflection undergoes one downward reflection at the surface 4 as shown in Figure 1. A second order free-surface multiple reflection undergoes two downward reflections at me surface 4, and so on. The present invention provides a technique for determining a signature of the source 6 and provides a technique for removing or reducing the effects of free-surface multiples, as will be described hereinafter. In order to illustrate this technique, the derivation of die method will be described with respect to a theoretical two-dimensional earth model comprising an inhomogeneous solid half-space overlain by a homogeneous fluid (water) layer. The sea surface 4 constitutes a free- surface because acoustic pressure vanishes in order to avoid infinite acceleration of the surface layer. The reflection coefficient at this interface is almost equal to -1. It is assumed mat the signal recorded from a typical hydrophone has had the effects of the direct wave and ghosting removed. Referring to horizontal (x) and vertical (z) spatial coordinates and a time coordinate t, die source position is denoted by (x s ,z s ) and die hydrophone (receiver) position is denoted by (x g ,z g ). The seismic data representing the pressure variation at (x g ,z g ) and at time t for a source located at (x.,z s ) is given by D 0 (x $ ,z,,t; x s ,z g ). By resetting the time variable to zero for each source and simplifying such that z g =z s =0, the seismic data may be equivalendy represented by D 0 (x s ,x g ,t). This is Fourier- transformed to die (w-k) domain with respect to x s ,x g ,t to give D 0 (k,,k g ,w), where k,,k., and w are the Fourier-transformed variables corresponding to x s ,x g , and t, respectively. To simplify notation, the Fourier-transformed data will be represented hereinafter by D 0 , and corresponding simplified notation will be used throughout. As disclosed by Carvalho, P.M., Weglein, A.B. and Stolt, R.H., 1991, Examples of a nonlinear inversion method based on die T matrix of scattering theory: Application to multiple suppression: Mtg. Soc. Expl. Geophys., Expanded Abstracts, 1319-1322, the scattering Born series for removing events due to free- surface reflections i.e. free-surface multiple reflections, from the seismic data D 0 can be written as: D p =D 0 +AD,+A 2 D 2 +... (1) where D p is the data without free-surface multiples and A=A(w) is the inverse of the Fourier-transformed source signature S(w), which is assumed to be only time dependent and not to vary with angle or source position i.e. A=l/S. The pressure fields D,,D 2> ... are given by oo D n (k g ,k s ,ω) = j dk. cos Θ.D σ (k g , k, ω).D n _ (k, k s , ω). (2) where n is a positive integer and The constant c is the seismic velocity of water and k is a generic horizontal wave number. The Born series in equation (1) for removing free-surface multiples is as follows. The series is constructed using the (Fourier-transformed) seismic data D 0 and the inverse source A only. The first term D 0 of the series is the actual data. The second term containing D, removes first order free-surface multiples. The next term including D 2 removes second order free-surface multiples, and so on. Thus, in order to remove die effects of free-surface multiple reflections, knowledge of die source signature, and hence its inverse A, is required. In order to determine the source signature (via its inverse A), it is necessary to find a solution to the equation (1), which is a polynomial in A. Equation (1) is a polynomial of high order and finding its solutions is therefore a difficult nonlinear problem. However, truncating of the polynomial to its first two terms corresponds to removal of first order free-surface multiples and yields a polynomial of first order which can be solved analytically. Thus, truncating equation (1) to its first two terms gives: D f = D σ + AD X + ε τ (4) where D f are the data from which the first order multiples have been removed and e ~ . describes the effects due to truncation, is small, and is nonlinearly related to the inverse source A. Removing the effects of first order free-surface multiples reduces e acoustic energy represented by die seismic data. When the source signature is correcdy determined, die reduction in energy will be maximised. Thus, it is required to find die source signature for which the energy of D p i.e. the data after the removal of first order multiples, is minimum. By using d e known least squares technique, the signature can be defined by that A which minimises: 5(A) = £> I 2 +μ| 2 (5) where: 12 D '/ = jdk g jdk jdω.D f .W D .D f * (6) and Λ|| 2 = σ 2 j dwj dw'.A(ω). W^ 1 (ω, ω').A * (ω f ) (7) The asterisk denotes complex conjugate. W D =W D (k, Jc g ,w) is a weighting function describing errors in the data and W A (w,w') describes a priori information about the source. The term || 2 is introduced to guarantee the stability of the solution. To simplify subsequent inversion formulae, the constant s 2 has been introduced into the definition of |j 2 . This criterion corresponds to finding die source which minimises the energy of data after the removal of the first order multiples. The source signature which permits the removal of the first order multiples will also permit the removal of higher order multiples so that, by finding die source signature in accordance with this technique, it may be used in equation (1) to remove free- surface multiples up to any desired order. The a priori information described by die weighting function W A (w,w') may be derived, for instance, from a deterministic estimate of the source signature, for instance using the technique disclosed in EP-A-0066423. Thus, a priori knowledge such as smoothness of the source spectrum may be incorporated through correlation between frequencies. In order to solve equation (5) so as to find the (inverse) source signature A, an initial analytical solution to equation (5) can be found (by ignoring the truncation errors (ε r )) by inversion to give: where N(ω) = \ dk g \ dk s .D 0 .W D .D; (9) and Q(ω) = \ dk g \ dk s .D ϊ .W D .D l * (10) and A <0> is die initial solution. In some circumstances, the initial solution will comprise a sufficiendy accurate and precise determination of die source signature so that it may men be substituted into equation (1) so as to provide data D p in which free-surface multiple reflections have been attenuated to a sufficient degree. Also, d e signature thus determined may be used for other purposes, for instance during subsequent data processing. However, where a more accurate determination of me source signature is required, this may be achieved by setting up an iterative scheme, of which A (0) is the starting solution, in order to accommodate die truncation errors which were ignored in die analytical solution represented by equation (8). Effectively, the nonlinear part of the problem is solved through the iterative linear solutions. Initialisation of die iterative scheme is as follows: and In general, A <0) permits a significant reduction of first order multiple energy through D f (0> . The iterations are set up to find die correction to A <0) which allows the removal of the remaining first order free-surface multiple energy. For each D 0 n) representing data containing the residual energy of die first order multiples, D 1 <n, can be found from equation (2). The correction δA (n) \o A ( " _l) can then be found by minimising This equation can be solved analytically in accordance with equation (8). D 0 is then replaced by D 0 (n) and D, is replaced by D, (n) . Finally, the source signature is updated in accordance witii: and die data are updated in accordance with D/ n) = D 0 (π) + ά4 (n) .zV π) (15) This iterative process may be stopped at any appropriate stage, for instance when a predetermined criterion is met. For instance, the iterations may be stopped when two successive solutions A (n) and A (n l> are sufficiently close (the absolute value of &\ {n) is less than a predetermined value). Alternatively, the iterations may be stopped when the reduction in energy represented by die data is less than a predetermined amount. When the source signature has been found to a required precision, it may be substituted into equation (1). This equation may be truncated at any desired term so as to remove or reduce the effects of free-surface multiples of orders up to and including die highest order of A in the truncated form of equation (1). In practice, the updated first order multiple term, D , in equation (15), can be kept identical to the initial first order multiple term, D^ because primary energy remains unchanged through iterations. The seismic reflection data required as die input for this technique are "shot gathers", for instance as described with reference to Figure 1, regularly sampled in space and time. It is assumed mat de-ghosting correction and subtraction of the direct waves have been performed on die input data. No assumptions of wavelet structure or phase are required in order for the technique to be performed. However, the source signature determined by this technique may be different from the actual source signature, depending on die effects of preprocessing of the data. The weighting function W A allows any a priori information about the source wavelet to be taken into account. For instance, in order to ensure that the source spectrum is smooth, this weighting function may be given as follows: (ω-ω')T W A (ω, ω') = sin c (17) The weighting function W D describes errors in the data. It may be used to ensure that certain parts of the data do not dominate the solution. For instance, this weighting function may be given by the following: • *.* » )={ffi_; where k- is about one tenth of the maximum wave number. Such a function rreedduucceess m diee iinnfflluueennccee ooff wwaavvee nniumbers in the neighbourhood k =0 (i.e. near the zero offset data) on the solution. The derivation of die present technique has been exemplified for a two dimensional eartii model. However, this can readily be modified to a three dimensional earth model with point sources. This simply requires that, in equations (9) to (16), the wave numbers k. and k β are replaced by vector wave numbers Further, the inverse source signature can be derived in die (x-w) domain. In tiiis case, equation (8) becomes: where N(ω) = j dx g j dx s .D 0 x s , x g , ω).D * (x s , x g ,ω) (21) and Q{ω) = J dx g J dx s .£> x s , x g , ω).W D (x, , x g , ω).D * (x_ , x g , ω) (22) The fields D 0 (x.,x g ,w) and D,(x s ,x g ,w) are the temporal Fourier transforms of D 0 (x s ,x g ,t) and D 0 (x s ,x g ,t), respectively. The weighting functions W D (x $ ,x g ,w) and W A (w,w ') must be suitably chosen for this domain. If, in addition to time variations, the source signature varies with angle or shot position, then d e technique described hereinbefore can be performed for A(w,k.) instead of A(w) in die (w -k) domain or for A(w,x $ ) instead of A(w) in die (x- w) domain. Figures 2 to 13 illustrate the use of this technique on synthetic data relating to an acoustic model of the earth comprising two homogeneous horizontal layers overlying a homogeneous half space. The deptii of the two horizontal interfaces are 120 and 285 metres. With increasing depth, d e velocities of the layers and the half space are 1500 meters per second, 2000 meters per second, and 2500 meters per second. A synthetic shot point gatiier was generated comprising 83 receivers witii a spacing of 12.5 meters. Figure 2 shows the syntiietic data as a graphical representation of recorded receiver data against receiver offset i.e. in the usual way for recorded seismic data. The data contain two primaries PI and P2 but several free-surface multiples. This corresponds to die first term DO in equation (1). Figure 3 shows the first two terms D 0 +AD, of equation (1). This shows the effect of removing first order free-surface multiples. Similarly, Figure 4 shows the first three terms D 0 +AD,+A D 2 , thus representing removal of the first and second order free-surface multiples. Figure 5 shows die first four terms D 0 +AD,+A 2 D 2 +A 3 D 3 i.e. with the first, second and tiiird order multiples removed. Figures 6 to 8 show the source signature or wavelet, the amplitude spectrum, and die phase spectrum, respectively, of die actual source signature (solid lines) and die initial estimated source signature (broken lines) before iteration to reduce die effects of truncation as described hereinbefore. Thus, most of the main characteristics of the source have been recovered by d e first estimation. In particular, the phase of the estimated source signature matches well with tiiat of the actual source signature. The result of free-surface multiple removal using this estimated source signature is shown in Figure 9. Some residual energy remains, for first order multiples in particular. Figures 10 to 13 correspond to Figures 8 to 9, respectively, but illustrate the result after five iterations. The match between the actual and estimated source signatures has significantly improved as shown in Figures 10 to 12. As shown in Figure 13, all of the free-surface multiple reflections have been successfully removed so mat only die primaries remain. It is thus possible to provide a technique which permits the signature of a seismic source to be determined without requiring any assumptions (about the source amplitude or phase or about the subsurface except the velocity of water in the marine survey) and without requiring any special configurations of hardware. Further, it is possible to provide a technique based on such source signature estimation to allow the effects of free- surface multiple reflections to be reduced or removed from seismic reflection data. The following situations generate multiple reflected energy in a manner similar to that of the marine seismic example described in detail hereinbefore. In such situations, the multiple energy can be used to estimate die source signature which can then be used for the removal or reduction of the multiple energy. (1) In vertical seismic profiling (VSP) applications, energy sources are located on the surface and receivers are located in die subsurface along a bore hole. Geophones closer to geological horizons can provide a more detailed stratigraphic picture. (2) Measurements can be derived from acoustic sources and receivers in a bore hole. Acoustic sources in a bore hole may be used to determine (a) die corrosion of the outer surface of the pipe, (b) die diird interface i.e. the bore hole wall, and (c) formation properties near the bore hole. These technologies are relevant to well development for drilling to production. Receivers are either coincident with die source or offset from it along the bore hole. The frequency range and radiation pattern of the source differ with different applications. In applications (b) and (c), reverberation of die wave in d e pipe can have a deleterious and disadvantageous effect. (3) In ground penetrating radar, electromagnetic reflection recorded on die surface of the earth is used to estimate near-surface properties. (4) In medical imaging and non-destructive material evaluation applications, a sample is irradiated with energy, such as x-rays or ultrasound, in order to determine the internal structure of the sample non-invasively.

Claims: CLAIMS 1. A method of determining a signature of an energy source from an inverse scattering series which represents multiple reflected energy and which comprises a polynomial in an inverse of the signature and has recorded data as die first term thereof, the method comprising: (a) selecting a sub-series comprising two terms of the inverse scattering series, the two terms comprising the first term and an Ntii order term where N is a positive integer, and (b) determining a value of the inverse signature which substantially minimises the energy represented by the sub-series. 2. A method as claimed in Claim 1, in which N is equal to 2 and die step (b) comprises analytically solving an expression for energy minimisation of the sub- series to determine die inverse signature. 3. A method as claimed in Claim 2, further comprising: (c) replacing the first term of the sub-series by the preceding value of die sub- series; (d) replacing the inverse signature in the sub-series by an inverse signature correction to form a modified sub-series; (e) determining die value of the inverse signature correction which substantially minimises me energy represented by die modified sub- series; and (f) adding die inverse signature correction to the inverse signature. 4. A method as claimed in Claim 3, in which the steps (c) to (d) are repeated until a predetermined criterion is achieved. 5. A method as claimed in Claim 4, in which me predetermined criterion is that die inverse signature correction is less than a predetermined value. 6. A method as claimed in Claim 4, in which the predetermined criterion is that die absolute value of die difference between consecutive values of die energy represented by the sub-series is less tiian a predefined value. 7. A metiiod as claimed in any one of the preceding claims, in which the recorded data comprise transformed recorded data. 8. A method as claimed in any one of the preceding claims, in which the recorded data comprise recorded data prcprocessed to remove or reduce direct wave and ghost effects. 9. A method as claimed in any one of the preceding claims, in which the recorded data comprise marine seismic data. 10. A metiiod as claimed in any one of die preceding claims, in which the recorded data comprise prestack seismic data. 1 1. A metiiod of reducing the effects of multiple reflected energy in recorded data, comprising: selecting the first M terms, where M is an integer greater than 1, of an inverse scattering series which represents free-surface multiple reflections and which comprises a polynomial in an inverse of an energy source signature; determining the energy source signature by a method as claimed in any one of die preceding claims; and substituting the determined signature into the selected first M terms to form data in which die effects of free-surface multiple reflections are reduced.

Inventors: IKELLE LUC THOMAS, ROBERTS GRAHAM, IKELLE, LUC, THOMAS, ROBERTS, GRAHAM
Applicants: GECO AS, IKELLE LUC THOMAS, PRAKLA SEISMOS GMBH, ROBERTS GRAHAM, GECO A.S., IKELLE, LUC, THOMAS, ROBERTS, GRAHAM
Application year: 1995
Main classification: G01V13
IPC codes: G01V1, G01V13, G10K11