DE2419578C2 - Procedure for eliminating as far as possible the multiples in seismic records - Google Patents

Description

(a) Die Ursprungsspur wird entweder über ein Wiener-Filter, dessen Parameter (insbesondere Vorhersageoperator ft) an die örtlichen Verhältnisse (nämlich an Schallwellengeschwindigkeit im zu untersuchenden Medium und Abstand zwischen Sender und zur Datenspur gehörendem Empfänger) angepaßt wird, einer Dekonvolutoperation unterzogen, oder es wird ein anderes an sich bekanntes Verfahren angewendet, um einen ersten Schätzwert des mehrfach reflektierten Signals zu erhalten (F i g. 8, Blöcke SA, SB, SC); (a) The original track is either subjected to a deconvolute operation using a Wiener filter, the parameters of which (in particular the prediction operator ft) are adapted to the local conditions (namely to the speed of sound waves in the medium to be examined and the distance between the transmitter and the receiver belonging to the data track), or Another method, known per se, is used to obtain a first estimate of the multiply reflected signal (FIG. 8, blocks SA, SB, SC);

(b) der Schätzwert wird an geeigneten, insbesondere die Energieminima darstellenden Stellen zu Zeitsegmenten aufgeteilt (F i g. 8, Block SD)-, (b) the estimated value is divided into time segments at suitable points, in particular those representing the energy minima (FIG. 8, block SD) -,

(c) die Zeitsegmente des Schätzwertes werden mit den entsprechenden Z^rtsegmenten des Ursprungssignals untu- Anwendung unterschiedlicher Zeitverschiebungen "or den Schätzwert kreuzkorreliert (Block SE)-, (c) the time segments of the estimated value are cross-correlated with the corresponding Z ^ rtsegmente of the original signal and the application of different time shifts "or the estimated value (block SE) -,

(d) die Zeitverschiebung, die die beste Korrelation ergibt, wird ermittelt (Block SE)-, (d) the time shift that gives the best correlation is determined (block SE) -,

(e) der jedem Segment des Schätzwertes für maximale Korrelation zuzuordnende Korrekturfaktor wird bestimmt (Block SE)-, (e) the correction factor to be assigned to each segment of the estimated value for maximum correlation is determined (block SE) -,

(f) jedes Segment des Schätzwertes wird mit dem im Schritt (e) ermittelten zugehörigen KorreK-turfaktor multipliziert und um die im Schritt (d) ermittelte Zeitverschiebung verschoben (Block SF) und die so modifizierten Segmente von den entsprechenden Segmenten der Ursprungsspur abgezogen (Block 8G^ um eine Primärspur zu erhalten.(f) each segment of the estimated value is multiplied by the associated correction factor determined in step (e) and shifted by the time shift determined in step (d) (block SF) and the segments modified in this way are subtracted from the corresponding segments of the original track (block 8G ^ to get a primary track.

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daß aus der Primärspur ein verbesserter Schätzwert gewonnen und damit die Verfahrensschritte (b) bis (f) wiederholt werden. 2. The method according to claim 1, characterized in that an improved one from the primary track The estimated value is obtained and the process steps (b) to (f) are repeated.

3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, daß aus einer oder mehreren einzelnen Ursprungsspuren oder Primärspuren ein einziger verbesserter Schätzwert gewonnen und dieser Schätzwert auf Ursprungsdatenspuren angewendet wird (Verfahrensschritte (b) bis (f)), die selbst keinen ausreichend zuverlässigen Schätzwert liefern. 3. The method according to claim 1 or 2, characterized in that of one or more individual traces of origin or primary traces a single improved estimate obtained and this estimated value is applied to original data traces (method steps (b) to (f)), which themselves do not provide a sufficiently reliable estimate.

4. Verfahren nach einem der Ansprüche 1 bis 3. dadurch gekennzeichnet, daß zur Ausführung der Verfahrensschritte (c), (d) und (e) analog arbeitende Kreuzkorrelatoren dienen.4. The method according to any one of claims 1 to 3, characterized in that for executing the Process steps (c), (d) and (e) analogously operating cross-correlators are used.

5. Verfahren nach einem der Ansprüche 1 bis 3, dadurch gekennzeichnet, daß zur Ausführung der Verfahrensschritte (b) und (e) analog arbeitende Vergleicher dienen.5. The method according to any one of claims 1 to 3, characterized in that for carrying out the Process steps (b) and (e) analogously operating comparators are used.

6. Verfahren nach einem der Ansprüche 1 bis 3, dadurch gekennzeichnet, daß zur Ausführung des6. The method according to any one of claims 1 to 3, characterized in that to carry out the

Verfahrensschrittes (f) analog arbeitende Multiplikatoren und Addierer dienen.Process step (f) analog multipliers and adders are used.

7. Verfahren nach einem der Ansprüche 4 bis 6, dadurch gekennzeichnet, daß die Verfahrensschritte (b) bis (f) mittels eines Analogrechners arsgeführt7. The method according to any one of claims 4 to 6, characterized in that the method steps (b) to (f) are performed by means of an analog computer

ίο werden, der zum einen den Schätzwert und zum anderen die Ursprungsspur als Eingänge aufnimmt und folgendermaßen aufgebaut ist:ίο be, on the one hand the estimated value and on the other hand other records the original track as inputs and is structured as follows:

Der Eingangsanschluß für den Schätzwert ist mit einer Segmentierungseinrichtung (SD) verbunden, die den Schätzwert quadriert und das Umhüllungssignal (F i g. 7) an den Minima in Segmente aufteilt, außerdem mit einer Zeit- und Amplitudenvergleichseinrichtung (SE), die auch mit dem Eingangsanschluß für die Ursprungsspur sowie mit dem Ausgang der Segmentierungseinrichtung verbunden ist und unter Steuerung des Ausgangs der Segmentierungseinrichtung zwischen dem zugeführten Ursprungssignai und dem Schätzwert Korrelationen, Zeit- und Amplitudenvergleiche vornimmt und für jedes Segment einen Amplitudenkorrekturfaktor und eine Zeitverschiebung ermittelt, wobei diese Werte einer Korrektureinrichtung (SF) zugeführt werden, die auch mit dem Eingangsanschluß für den Schätzwert verbunden ist und am Schätzwert die ermitteltenThe input connection for the estimated value is connected to a segmentation device (SD) which squares the estimated value and divides the envelope signal (FIG. 7) into segments at the minima, as well as to a time and amplitude comparison device (SE), which is also connected to the The input connection for the original track and the output of the segmentation device is connected and, under control of the output of the segmentation device, it carries out correlations, time and amplitude comparisons between the supplied original signal and the estimated value and determines an amplitude correction factor and a time shift for each segment, these values being a correction device ( SF) , which is also connected to the input connection for the estimated value, and the determined value at the estimated value

Jf) Korrekturen vornimmt. Der Ausgang der Korrektureinrichtung und der Eingangsanschluß für die Ursprungsspur sind an einen Addierer (8C^ geführt, der die negierte korrigierte Schätzung und die Ursprungsspur addiert und die Summe als Primärspur an einen Ausgangsanschluß (S 2) liefert.Jf) making corrections. The output of the correction device and the input connection for the original track are fed to an adder (8C ^ which adds the negated corrected estimate and the original track and delivers the sum as a primary track to an output connection (S 2).

8. Verfahren nach einem der Ansprüche 1 bis 3, dadurch gekennzeichnet, daß die Verfahrensschritte (b) bis (f) mittels eines Digitalrechners ausgeführt werden.8. The method according to any one of claims 1 to 3, characterized in that the method steps (b) to (f) are carried out using a digital computer.

Die Erfindung betrifft ein Verfahren zur weitmöglichsten Beseitigung der Multiplen aus einer oder mehreren grafisch oder phonographisch reproduzierbaren, insbesondere digital abgetasteten und gespeicherten seismischen Datenspuren seismischer Wellen, die von einem Sender ausgehend mehrere vom Sender und voneinander einen Abstand aufweisenden Empfänger direkt und/oder nach mehrfacher Reflexion an horizontalen unterhalb der Oberfläche liegenden Schichten erreichen. The invention relates to a method for the greatest possible elimination of the multiples from one or more graphically or phonographically reproducible, in particular digitally scanned and stored seismic Data traces of seismic waves emanating from a transmitter several from the transmitter and from each other a spaced receiver directly and / or after multiple reflections at horizontal reach layers below the surface.

Ein derartiges Verfahren ist aus der US-PS 32 23 967 bereits bekannt.Such a method is from US-PS 32 23 967 already known.

Das in dieser Druckschrift dargestellte Verfahren arbeitet dadurch, daß mehrere Eingangssignalkurven zueinander in Beziehung gesetzt werden, um dadurch seismische Interferenzstörung^ zu eliminieren. Auch beim Erfindungsgegenstand wird mit zahlreichen Meßstationen gearbeitet, die ebenfalls bezüglich einer Signalquelle in bestimmter Weise angeordnet sind, wie es im Zusammenhang mit der Fig.6 noch näher erläutert werden wird. Die in dieser Figur dargestellten Spuren beruhen auf der Anordnung von 24 Empfangsstationen, die jeweils im Abstand von 50 m zueinander angeordnet sind, wobei der zum Sender nächstliegende Empfänger von diesem eine Entfernung von 300 mThe method shown in this document works by having several input signal curves can be related to each other to thereby eliminate seismic interference disturbance ^. Even When the subject invention is worked with numerous measuring stations, which also with respect to a Signal source are arranged in a certain way, as will be explained in more detail in connection with FIG will be explained. The tracks shown in this figure are based on the arrangement of 24 receiving stations, which are each arranged at a distance of 50 m from one another, with the one closest to the transmitter Receiver from this at a distance of 300 m

besitzt Die in F ί g. 6 unter a angegebenen Aufzeichnungsspuren stellen die Originalwiedergabe der Empfangssignale dieser Empfänger dar.possesses the in F ί g. 6 recording tracks indicated under a represent the original reproduction of the received signals of these receivers.

Für jede dieser einzelnen Ursprungsspuren wird nun in an sich bekannter Weise eine sog. »Schätzung« gewonnen, wobei die Art und Weise der Gewinnung dieser Schätzung an sich ohne Bedeutung ist Nur beispielsweise sei auf die Möglichkeit verwiesen, für diesen Zweck ein sog. Wiener-Filter zu verwenden, wie es in der Zeiu-hrift »Geophysics«, Bd. 34, Nr. 2, April 1969, Seiten 155-169 in großen Einzelheiten beschrieben wird. Auch die DE-OS 16 23 552 (sie ist analog zu der GB-PS 11 95 204 bzw. der US-PS 34 96 529) finden sich ebenfalls genaue Erläuterungen für Arbeitsweise und gerätetechnische Verwirklichung bestimmter Vorgänge, die als für den Erfindungsgegenstand bekannt unterstellt werden können und daher im folgenden nicht in größeren Einzelheiten erläutert werden.For each of these individual traces of origin, a so-called "estimate" is now made in a manner known per se. won, the manner in which this estimate was obtained is in itself of no importance. Only for example, reference is made to the possibility of using a so-called Wiener filter for this purpose, such as see it in the journal "Geophysics", Vol. 34, No. 2, April 1969, pages 155-169 is described in great detail. DE-OS 16 23 552 (it is analogous to GB-PS 11 95 204 and US-PS 34 96 529) there are also precise explanations for the method of operation and the implementation of certain processes in terms of equipment, which can be assumed to be known for the subject matter of the invention and therefore not in the following will be explained in greater detail.

Wie der bereits genannten Zitatstelle aus »Geophysics« entnommen werden kann, kann dieses an sich bekannte »Vorhersageverfahren« mittels verschiedenartiger Filteranordnungen verwirklicht werden, so z. B. mit dem inversen Filter der kleinsten Vielehen mit bestimmter Vorhersageentfernung, wobei durch Wahl dieser Vorhersageentfernung (bzw. Vorhersagezeit) das Maß der Auflösung festgelegt wird.As the already mentioned quotation from »Geophysics« can be inferred, this "prediction method", which is known per se, can be obtained by means of various types Filter assemblies are realized, so z. B. with the inverse filter of the smallest polygons with certain prediction distance, whereby by choosing this prediction distance (or prediction time) the The degree of resolution is determined.

Tastet man die zunächst analog vorliegenden seismischen Signale digital ab, kann ein Schätzverfahren angewendet werden, das, worauf bereits verwiesen wurde, als »Wiener-Filter-Verfahren« bezeichnet wird. Mit Hilfe von Wiener-Filtern können periodische Vorgänge mit festgelegten PeriodenJauern durch entsprechende Wahl der Parameter des Filters beseitigt werden. Mathematisch läßt sich der Vorgang durch eine Konvoluttransformation wiedergeben, jedoch braucht bei der technischen Verwirklichung nicht unbedingt ein Computer herangezogen zu werden, um diese mathematische Operation auch tatsächlich auszuführen, obwohl bei digital vorliegenden Werten sich diese Art der Berechnung anbietet. Statt dessen kann aber auch ein Analogverfahren mit Hilfe von entsprechendeii Analogschaltungen angewendet werden, oder aber auch, wenn dies auch recht unpraktisch und zeitraubend ist, eine Umsetzung von Hand mit zeichnerischen Mitteln.If the seismic signals, which are initially available in analog form, are digitally sampled, an estimation method can be used what has already been referred to is referred to as the "Wiener filter method". With the help of Wiener filters, periodic processes with defined periods can be carried out corresponding choice of the parameters of the filter can be eliminated. Mathematically, the process can be described by a Reproduce convolute transformation, but does not necessarily need a in the technical implementation Computers to be used to actually perform this mathematical operation, although this type of calculation is suitable for digital values. Instead, however, can also an analog process with the help of correspondingeii Analog circuits are used, or else, if this is also quite impractical and time-consuming is, a realization by hand with drawing means.

In der noch folgenden Beschreibung ist diese Verfahrensweise in F i g. 8 durch die Blöcke SA, SB und SC näher erläutert. Der Ausgang des Blockes SC stellt die aus der Ursprungssignalkurve gewonnene erste Schätzung dar, die dann zusammen mit der Ursprungskurve selbst im eigentlichen erfindungswesentlichen Teil, der die Blöcke SD bis SG umfaßt, verarbeitet werden, um eine besser auswertbare Spur zu erhalten, d. h. eine Spur, die weniger durch Mehrfachreflexionen (Multiple) gestört ist.In the description that follows, this procedure is shown in FIG. 8 explained in more detail by the blocks SA, SB and SC. The output of the block SC represents the first estimate obtained from the original signal curve, which is then processed together with the original curve itself in the actual essential part of the invention, which comprises the blocks SD to SG , in order to obtain a track that can be better evaluated, ie a track, which is less disturbed by multiple reflections.

Auch für die Zeit- und Amplitudenkorrekturen lassen sich bestimmte mathematische Ausdrücke angeben, durch die sich die entsprechende Korrektur zwar exakt definieren läßt, die aber bei der praktischen Ausführung des Verfahrens nicht unbedingt angewendet werden müssen, insbesondere dann nicht, wenn die Signale in analoger Form vorliegen. So beschreibt die DE-OS 16 23 552 sehr ausführlich analoge Korrelatoren, wobei diese Korrelatoren bestimmte Signalwerte am Ausgang liefern, die als »Korrelationskoeffizient« bezeichnet werden und - bei richtiger Eichung - zwischen O und 1 liegen, wobei O nicht vorhandene Korrelation und I vollständige Korrelation der miteinander korrelierten Größen bedeuten. Vollütändige Korrelation besagt, daß beide Signale bezüglien ihres Signalverlaufs identisch sind, wobei die Identität zeitgleich sein kann, so daß die Korrelation von 1 für eine Korrelationszeitverschie-Certain mathematical expressions can also be specified for the time and amplitude corrections, through which the corresponding correction can indeed be defined exactly, but in the practical implementation of the method do not necessarily have to be used, especially not if the signals in analogous form. For example, DE-OS 16 23 552 describes analog correlators in great detail, with these correlators deliver certain signal values at the output, which are referred to as the "correlation coefficient" and - with correct calibration - lie between O and 1, where O correlation does not exist and I mean complete correlation of the quantities correlated with one another. Complete correlation says that both signals are identical with respect to their signal course, the identity can be the same, so that the Correlation of 1 for a correlation time different

--, bung von O sich ergibt, oder sie kann bezüglich zeitverschobener Stgnalteile vorhanden sein, wobei sich die hohe Korrelation erst bei einer bestimmten Korrelationszeitverschiebung zeigt Bei seismischen Signalen liegen häufig derartige in engem Zusammen- -, exercise of O results, or it can be present with respect to time-shifted signal parts, the high correlation only showing up at a certain correlation time shift.

η hang stehenden Signale vor, die dann allerdings eben um ein bestimmtes Ausmaß zeitverschoben sind, bedingt durch die unterschiedlich lange Laufzeit Korreliert man diese beiden Signale miteinander, ergibt sich ein höchster Wert für die Korrelation — der allerdingsη hang present signals, which are then, however, time-shifted by a certain amount, conditionally Due to the different running times, if these two signals are correlated with one another, the result is a highest value for the correlation - but that one

ii nicht unbedingt 1 sein muß, da Störsignale und unterschiedliche durchlaufene Gebiete zu weiteren Abweichungen führen - bei einer Zeitverschiebung, die dem Laufzeitunterschied entsprichtii does not necessarily have to be 1, since interfering signals and different areas traversed lead to further deviations - with a time difference, the corresponds to the runtime difference

Sind in einem seismischen Signal mehrere Reflexio-If there are several reflection points in a seismic signal

JD nen vorhanden, die aber alle bestimmte Gemeinsamkeiten aufweisen und nur zueinander Teitlich verschoben sind, wie es bei Mehrfachreflexioriei' meist der Fall ist. kann man nun durch Korrelation der einzelnen reflektierten Anteile miteinander (siehe in F i g. 7 dieJD exist, but they all have certain similarities have and only partially shifted to each other are, as is usually the case with multiple reflections. one can now by correlating the individual reflected components with one another (see in FIG. 7 die

2ϊ einzelnen «Gebirgszüge« der Signalkurve) Einzelheiten innerhalb dieser Signalabschnitte genauer untersuchen.2ϊ individual "mountain ranges" of the signal curve) details examine more closely within these signal sections.

Der Korrelator, dem die zu korrelierenden SignaleThe correlator to which the signals to be correlated

zugeführt werden, liefert also sowohl bestimmte Korrelationswerte (Amplitudenwerte) als auch Zeitver-are supplied, so it provides both specific correlation values (amplitude values) and timing

jii schiebungswerte, die z. B. maximale Korrelationswerte zwischen den beiden Vorgängen repräsentieren.jii shift values that z. B. Maximum correlation values represent between the two processes.

Weitere Vorschläge zur richtigen Auswertung seismischer Aufzeichnungen, die durch Multiple überlagert werden, werden in den nachstehenden Veröffentlichun-Further suggestions for the correct evaluation of seismic records overlaid by multiple will be discussed in the following publications

ii gen beschrieben, auf die in der folgenden Beschreibung wieder Bezug genommen wird:ii genes described in the following description is referenced again:

A. Backus, M. M., 1959, »Water reverberations, their nature and elimination«: Geophysics, Bd. 24, S.A. Backus, M. M., 1959, "Water reverberations, their nature and elimination": Geophysics, Vol. 24, pp.

w 233-261: w 233-261:

B. Robinson, E. Α., 1967, »Predictive decomposition of time series with application to seismic exploration« : Geophysics, Bd. 32, S. 418 - 484;B. Robinson, E. Α., 1967, "Predictive decomposition of time series with application to seismic exploration ": Geophysics, Vol. 32, pp. 418-484;

C. Peacock, K. L. und Treitel, S.. 1969, »Predictive 4", deconvolution: theory and practice«: Geophysics, Bd. 34, S. 155-169;C. Peacock, K. L. and Treitel, S .. 1969, "Predictive 4", deconvolution: theory and practice ": Geophysics, Vol. 34, pp. 155-169;

D. Newman. P., 1973, »Divergence effects in a layered earth«: Geophysics, Bd. 38, Juni 73, S.481;D. Newman. P., 1973, "Divergence effects in a layered earth": Geophysics, Vol. 38, Jun. 73, p. 481;

E. Anstey, N. A. und Newman, P„ GB-PS 11 95 204 -mi (US-PS34 96 529).E. Anstey, N.A. and Newman, P "GB-PS 11 95 204 -mi (US-PS 34 96 529).

Aufgabe der vorliegenden Erfindung ist es, das Verfahren der eingangs genannten Art hinsichtlich der Dämpfung der Multiplen bei der Gewinnung seismiv. scher Aufzeichnungen noch weiter zu verbessern.The object of the present invention is to provide the method of the type mentioned with regard to Seismic attenuation of the multiples during extraction. shear records even further.

Gelöst wird dier.e Aufgabe durch die Kennzeichnen den Merkmale (a) bis (f) des Hauptanspruchs.The task is solved by the labeling features (a) to (f) of the main claim.

Erfindungsgemäß wird also zunächst im Verfahrensschritt (a) die Ursprungsspur entweder über ein «ι Wiener-Filter, dejsen Parameter (insbesondere Vorhersageoperator ft) an die örtlichen Verhältnisse (nämlich an Schallwellengeschwindigkeit im zu untersuchenden Medium im Abstand zwischen Sender und zur Datenspur gehörendem Empfänger) angepaßt wird, <r> einer Dekonvolutoperation unterzogen, oder auch ein anderes an sich bekannte,) Verfahren angewendet, um einen ersten Schätzwert des mehrfach reflektierten Signals zu erhalten, wobei auf die Blöcke SA, SB. 8CderAccording to the invention, the original track is first adapted to the local conditions (namely to the sound wave speed in the medium to be examined at the distance between the transmitter and the receiver belonging to the data track) using a Wiener filter, dejsen parameter (in particular the prediction operator ft) is, <r> subjected to a deconvolute operation, or another method known per se is used to obtain a first estimate of the multiply reflected signal, with the blocks SA, SB. 8Cder

Fig. 8 verwiesen werden kann. Dann wird gemäß den weiteren Verfahrensschritten die Ursprungssignalspur mit der «geschätzten« Signalspur, die bereits weniger Störungen enthält, aber möglicherweise die tatsächlichen Verhältnisse noch nicht ganz korrekt wiedergibt, abschnittweise korreliert, wobei diese Abschnitte nach bestimmten Zweckmäßigkeitserwägungen ausgewählt werden. Meistens wird man als Trennpunkt die Minima wählen, siehe z. B. die in F i g. 7 dargestellte quadrierte und mit einer Hüllkurve versehene, gemäß dem Wiener-Verfahren geschätzte Spur, die der untersten Spur (300 m Entfernung von der Signalquelle) der F i g. 6.6b. entspricht. Die Feststellung der Minima kann z. B. durch analoge elektronische Schaltungen erfolgen, alternativ natürlich auch durch digitale Anordnungen, falls die Spur in Form von Digitalwerten vorliegt. Der in Fig. 8 dargestellte Block SD enthält eine derartige Vorrichtung zur Feststellung der Minima, wobei dieser Block aufgrund der Minima dann auch die entsprechende Zeitauftrennung vornimmt und die zugehörigen Steuersignale an den Zeit- und Amplitudenvergleicher 8£ liefert, dem zum Zwecke des Vergleiches zum einen der Schätzwert, zum anderen der Originalwert der gerade untersuchten Spur zugeführt wird. In dieser Vergleichseinrichtung SE wird dann segmentweise eine Korrelation der beiden Signalwerte vorgenommen, wobei unterschiedliche Zeitverschiebungen angewendet werden. Zum Beispiel kann dieser Block SE analog arbeitende Korrelatoren aufweisen, die die beiden zu vergleichenden Signale kreuzkorreliercn und dabei z. B. schrittweise die Zeitverschiebung ändern, bis sich ein maximaler Wert für die Korrelation ergibt. Der zugehörige Zeitverschiebungswert wie auch das Verhältnis dieses maximalen Korrelationswertes zum Korrelationswert mit Zeitverschieoung 0 (entspricht dem in der späteren Beschreibung definierten Faktor K) werden dann benutzt, um den ersten Schätzwert entsprechend zu korrigieren, wie durch den Block SF der F i g. 8 angedeutet. Die so verbesserte Wiedergabe der Mehrfachreflexionen wird dann von der Ursprungskurve abgezogen, siehe Block SC, F i g. 8. um ein Ausgangssignal 52 zu erhalten, das im wesentlichen frei is. von den Mehrfachreflexionen.Fig. 8 can be referenced. Then, according to the further process steps, the original signal trace is correlated section by section with the “estimated” signal trace, which already contains less interference but may not yet reflect the actual conditions completely correctly, these sections being selected according to certain considerations of expediency . Most of the time you will choose the minimums as the dividing point , see e.g. B. the in F i g. 7, shown squared and provided with an envelope, estimated according to the Wiener method, which is the lowest track (300 m distance from the signal source) of FIG. 6. 6b. is equivalent to. The determination of the minima can, for. B. by analog electronic circuits, alternatively, of course, by digital arrangements, if the track is in the form of digital values. The block SD shown in Fig. 8 contains such a device for determining the minima, this block then also undertakes the appropriate time separation based on the minima and supplies the associated control signals to the time and amplitude comparator 8 £, the one for the purpose of comparison the estimated value and, on the other hand, the original value of the track being examined. In this comparison device SE , a correlation of the two signal values is then carried out in segments, with different time shifts being used. For example, this block SE can have analog working correlators which cross-correlate the two signals to be compared and thereby z. B. gradually change the time shift until there is a maximum value for the correlation. The associated time shift value as well as the ratio of this maximum correlation value to the correlation value with time shift 0 (corresponds to the factor K defined in the description below) are then used to correct the first estimated value accordingly, as indicated by the block SF in FIG. 8 indicated. The reproduction of the multiple reflections improved in this way is then subtracted from the original curve, see block SC, FIG. 8. To obtain an output signal 52 that is substantially free. from multiple reflections.

Auch diese Korrektur der Schätzung sowie üic Subtraktion der verbesserten Schätzung von der Originalkurve können analog oder auch digital vorgenommen werden, je nachdem in welcher Form die Signale vorliegen.This correction of the estimate as well as subtraction of the improved estimate from the Original curves can be made analog or digital, depending on the form in which the Signals are present.

Gemäß einer im Anspruch 2 beschriebenen vorteilhaften Weiterbildung des erfindungsgemäßen Verfahrens werden zum Zwecke der noch besseren Abschätzung bestimmte Verfahrensschritte wiederholt, nämlich die Schritte der Segmentierung. Kreuzkorrelation und Korrekturfaktorbestimmung sowie die Anwendung der Korrekturfaktoren auf die Ausgangsschätzwerte. According to an advantageous development of the method according to the invention described in claim 2, certain method steps, namely the segmentation steps, are repeated for the purpose of even better estimation. Cross-correlation and correction factor determination as well as the application of the correction factors to the initial estimated values.

Statt für jede Spur das Verfahren durchzuführen, kann es gemäß einer noch anderen Ausführungsform der Erfindung zweckmäßig sein, mehrere Spuren zusammenzufassen und aus diesen einen einzigen verbesserten Schätzwert zu gewinnen und diesen dann wieder auf die einzelnen Ursprungsspuren getrennt anzuwenden. Dieses Verfahren ist besonders dann anzuwenden, wenn der Störanteil bei den einzelnen Spuren verhältnismäßig hoch ist und erst durch Zusammenfassung mehrerer Spuren ein ausreichender Störabstand sich ergibtInstead of performing the method for each track , according to yet another embodiment of the invention it can be expedient to combine several tracks and to obtain a single improved estimated value from these and then to apply this again separately to the individual original tracks. This method is particularly useful if the noise component in the individual tracks is relatively high and a sufficient signal-to-noise ratio is only obtained by combining several tracks

Die weiteren Unteransprüche beschäftigen sich mit bestimmten gerätetechnischen Aspekten zur Ausführung des erfindungsgemäßen Verfahrens.The further subclaims deal with certain device-technical aspects for the execution of the method according to the invention.

Die Erfindung wird nachfolgend anhand von Ausführungsbeispielen näher erläutert, die in den Zeichnungen dargestellt sind.The invention is explained in more detail below on the basis of exemplary embodiments which are shown in the drawings are shown.

Es zeigtIt shows

Fig. 1 in der Form idealisierter Diagramme la, Xb und Ic einige der theoretischen Wege, längs deren sich ein seismisches Signal bei nur einer reflektierenden Schichtgrenze fortpflanzen kann, Fig. 1 in the form of idealized diagrams la, Xb and Ic some of the theoretical paths along which a seismic signal can propagate with only one reflective layer boundary,

F i g. 2 in der Form ähnlicher Diagramme 2a, 2b und 2c einige theoretische Wege längs deren sich ein seismisches Signal bei zwei reflektierenden Schichtgrenzen fortpflanzen kann, F i g. 2 in the form of similar diagrams 2a, 2b and 2c some theoretical paths along which a seismic signal can propagate at two reflective layer boundaries,

F i g. 3 eine idealisierte seismische Aufzeichnung, wie sie sich aus der mit Fig. I wiedergegebenen Situation ergeben könnte, wobei aufeinanderfolgende Signalkurven die Signale wiedergeben, die an Stellen unterschiedlicher Entfernung von der Quelle empfangen werden,F i g. 3 an idealized seismic recording, as it results from the situation reproduced with FIG could result, with successive signal curves reproduce the signals at different points Distance from the source are received,

Fig. 4 in gleicher Weise eine idealisierte seismische Aufzeichnung, wie sie sich äüä der mü Tig. 2 wiedergegebenen Situation ergeben kann, aus der jedoch Größen oder Ereignisse, wie sie entsprechend F i g. 1 auftreten können, ausgeschlossen worden sind,4 shows an idealized seismic recording in the same way as it appears from the mü Tig. 2 reproduced situation can result , from which, however, variables or events as they are in accordance with FIG. 1 can occur, have been excluded,

Fig.5 eine Zusammenfassung der Fig.3 und 4, die somit beide Arten von Mehrfachreflexionen enthält, wie sie sich aus den in F i g. 1 und 2 dargestellten Situationen ergeben können.5 shows a summary of FIGS. 3 and 4, which thus contains both types of multiple reflections, as can be seen from the FIGS. 1 and 2 can arise.

Fig. 6 und 7 bei einer Untersuchung gemäß dem erfin<2,s.igsgemäßen Verfahren tatsächlich gemachte Aufzeichnungen undFIGS. 6 and 7 actually made during an examination according to the method according to the invention Records and

F i g. 8 eine schematische Darstellung zur noch detaillierteren Erläuterung des erfindungsgemäßen Verfahrens.F i g. 8 is a schematic illustration to explain the invention in more detail Procedure.

Im einzelnen lassen Fig. 1 und 2 eine Oberfläche S (Wasser oder '-and), eine Quelle P für seismische Signale, einen Empfänger Q für die Signale und eine Grenzfläche A (beispielsweise der Meeresboden) erkennen, der einen ersten Reflektor bildet. Drei mögliche Wege oder Bahnen, wie sie die reflektierten seismischen Signale nehmen können, sind mit den Bahnen R₁. R₂ bzw. R_s angedeutet.1 and 2 show a surface S (water or sand), a source P for seismic signals, a receiver Q for the signals and an interface A (for example the sea floor) which forms a first reflector. Three possible paths or tracks how the reflected seismic signals can take are with tracks R ₁ . R ₂ and R _s indicated.

F i g. 2 zeigt außerdem einen zweiten Reflektor, der durch eine Grenzfläche B gebildet wird, wobei eine Anzahl der vielen möglichen Wege, denen die reflektierten Signale folgen können, mit den Bahnen A₄. /?5 bzw. Rb angedeutet ist.F i g. 2 also shows a second reflector, which is formed by a boundary surface B, wherein a number of the many possible paths which the reflected signals may follow, with the tracks A _fourth /? 5 or Rb is indicated.

In dem Diagramm la der Fig. 1 wird das Signal entsprechend der Bahn /?, auf seinem Weg von der Quelle P zum Empfänger Q nur einmal reflektiert und entsprechend als die primäre Reflexion oder Primärreflexion bezeichnet Entsprechend dem Diagramm Xb wird das Signal auf der Bahn Ri von der Quelic P zum Empfänger Q dreimal reflektiert, während es entsprechend Fig. Ic auf seinem der Bahn A3 entsprechenden Weg fünfmal reflektiert wird. Signale entsprechend den Bahnen Ri oder Ry, die also mehr als einmal reflektiert worden sind, werden als Mehrfachreflexionen erster bzw. zweiter Ordnung bezeichnet Es versteht sich, daß Mehrfachreflexionen anderer Ordnung auftreten können, wenn mehr Reflexionen als mit F i g. 1 gezeigt stattfinden. In the diagram la of FIG. 1, the signal corresponding to the path / ?, on its way from the source P to the receiver Q is reflected only once and is referred to as the primary reflection or primary reflection. According to diagram Xb , the signal is on the path Ri is reflected three times from the Quelic P to the receiver Q , while it is reflected five times on its path corresponding to the path A3, as shown in FIG. Signals corresponding to the paths Ri or Ry, which have thus been reflected more than once, are referred to as multiple reflections of the first or second order. It goes without saying that multiple reflections of a different order can occur if more reflections than with F i g. 1 shown take place.

Das im Diagramm 2a der Bahn Rt folgende Signal wird an der Grenzfläche A gebrochen, jedoch erfolgt an der Grenzfläche ßnur eine einzige Reflexion, so daß das der Bahn Λ folgende Signal eine Primärreflexion darstellt Das der Bahn A₅ des Diagramms 2b folgende Signal wird in der ersten Lage zweimal reflektiert, einmal davon an der Grenzfläche A, und einmal wird esThe following in the graph 2a of the web Rt signal is refracted at the interface A, however, takes place at the interface ßnur a single reflection, so that the path Λ following signal representing a primary reflection, the train A ₅ of the diagram 2b following signal is in the first layer is reflected twice, once at the interface A, and once it is

an der Grenzfläche H reflektiert. Entsprechend dem Diagramm 2c wird das Signal in der ersten Lage viermal und einmal an der Grenzfläche ö reflektiert. Es versteht sich wiederum, daß diese Diagramme nur einige der unendlichen Zahl möglicher Wege wiedergeben, denen die Signale folgen könnten.reflected at the interface H. According to diagram 2c, the signal in the first position is reflected four times and once at the interface δ. Again, it should be understood that these diagrams depict only some of the infinite number of possible paths that the signals could follow.

In Fig.3-6 ist in Abhängigkeit vom Abstand d(m) der einzelnen Beobachtungspunkte von der Signalquelle die ,> iit ' (s) der empfangenen Signale aufgetragen. Fig./ zeigt die quadrierte Einhüllende der Mehrfachreflexions-Schätzung entsprechend der ersten Signalspur der Aufzeichnung 6b in Abhängigkeit von der Zeit t. In FIG. 3-6, the 'iit' (s) of the received signals is plotted as a function of the distance d (m) of the individual observation points from the signal source. Fig./ shows the squared envelope of the multiple reflection estimation corresponding to the first signal trace of the recording 6b as a function of the time t.

Die mit F i g. I und 2 wiedergegebene seismische Situation ist insbesondere für marine Untersuchungen maßgeblich, wo der Wasserkige starke Mehrfachreflexionen zugeordnet sein könren. Backus (Veröff. A) berücksichtigte die Auswirkungen einer solchen Wasserlage im Rahmen einer iinesr Filteriheorie und schlug die Verwendung inverter Filter vor, um mehrfachThe with F i g. I and 2 reproduced seismic Situation is particular for marine investigations decisive where the watery strong multiple reflections can be assigned. Backus (publ. A) took into account the effects of such a water situation as part of an iinesr filter series and suggested the use of inverter filters to multiple times

schiedliche Typen sog. Backus Filter-Operatoren notwendig, um Mehrfachreflexi'inen der beiden verschiedenen Reflexionsarten zu beseitigen, wie sie mit F i g. 1 und 2 gezeigt sind. Keiner dieser Backus-Filter-Operatoren kann dabei vollständig zufriedenstellend arbeiten, wenn auf die zusammengesetzte seismische Auszeichnung eingewirkt wird, wie sie mit F i g. 5 wiedergegeben ist. Außerdem häng* die Gültigkeit dieser inversen Filteroperatoren von den Annahmen ab. daß die Wellen sich eben ausbreiten, die Lagen horizontal verlaufen und die seismischen Wellenfronten normal auf die reflektieren·. :n Oberflächen auftreffen. Diese Annahmen oder Voraussetzungen stehen im Widerspruch zu praktischen seismischen Untersuchungen.Different types of so-called Backus filter operators are necessary to avoid multiple reflections of the two different ones Eliminate types of reflection as shown in FIG. 1 and 2 are shown. None of these Backus filter operators can work completely satisfactorily when on the compound seismic markup is acted upon, as shown with FIG. 5 is reproduced. In addition, the validity of this inverse depends Filter operators depend on the assumptions. that the waves spread evenly, the layers run horizontally and the seismic wavefronts reflect normally on the ·. : n hit surfaces. These assumptions or Conditions contradict practical seismic investigations.

Die Auslegung von Backus-Filter-Operatoren verlangt die Kenntnis oder Schätzung der normalen seismischen Einfalls-Zweiweg-Fortpflanzungszeit in der Nachhall-Lage sowie aer Reflexionskocffizienten an den Grenzen dieser Lage. Andere Methoden der Ausgestaltung oder Auslegung von inversen Filtern, die diese Kenntnis nicht ausdrücklich erfordern und die Nachhall- oder Reflexionsproblemen, die mehr als eine einfache Lage betreffen, gleichzeitig gerecht werden können, sind beispielsweise von Robinson (Veröff. B) sowie Peacock und Treitel (Veröff. C) beschrieben worden. Diese weiteren Verfahren beruhen auf statistischen Kennwerten des Echoproblems, wie sie sich aus der Korrelationsanalyse aufgezeichneter seismischer Signalkurven herleiten lassen. Die Annahmen, daß die Wellen sich in einer Ebene fortpflanzen und die seismischen Wellenfronten normal auftreffen, sind aber auch bei diesem Verfahren implizit enthalten, da die Analyse davon ausgeht, daß das Nachhall- oder Echoproblem über erhebliche Zeitabschnitte der Aufzeichnungsdauer sich stationär verhält The design of Backus filter operators requires knowledge or estimation of the normal seismic incidence two-way propagation time in the reverberation layer as well as the reflection coefficients at the boundaries of that position. Other methods of designing or designing inverse filters that do not explicitly require this knowledge and that can deal with reverberation or reflection problems that affect more than a simple layer at the same time are, for example, by Robinson (publ. B) and Peacock and Treitel ( Publication C) has been described. These further methods are based on statistical characteristic values of the echo problem, as can be derived from the correlation analysis of recorded seismic signal curves. The assumptions that the waves propagate in a plane and that the seismic wave fronts strike normally are also implicit in this method, since the analysis assumes that the reverberation or echo problem remains stationary over considerable periods of the recording period

F i g. 1 und 2 zeigen, daß die bei Verschiebung des Empfängers gegenüber der seismischen QueUe für verschiedene Reflexionsstellen maßgeblichen Einfallswinkel von der Ordnung der Mehrfachreflexion abhängen. Aufgrund geometrischer Überlegungen ist es klar, daß aufeinanderfolgend höhere Ordnungen von Mehrfachreflexionen nicht mit einer genau regelmäßigen Periode auftreten, wie das der Fall wäre, wenn der Empfänger nicht gegenüber der Quelle versetzt würde. Daher ist das Zeitintervall zwischen aufeinanderfolgenden Mehrfachreflexions-Größen oder -Ereignissen in der idealisierten Wiedergabe einer seismischen Aufzeichnung entsprechend Fig.3 und 4 nicht konstant sondern es nimmt allgemein mit späterer Aufzeichnungszeit zu. Ebenso ist es klar, daß jede Änderung in der Stärke der mehrfach reflektierenden Lage oder der Neigung der Grenzen über einen Entfernungsabschnitt, , der die Reflexionsstellen enthält, zu einer weiteren Abweichung von einem konstanten Zeitintervall zwischen aufeinanderfolgenden Mehrfachreflexions-Größen oder -Ereignissen führt.
In einer praktischen seismischen Situation treten in F i g. 1 and 2 show that the angles of incidence which are decisive for different reflection points when the receiver is shifted with respect to the seismic source depend on the order of the multiple reflection. From geometrical considerations it is clear that successively higher orders of multiple reflections do not occur with an exactly regular period, as would be the case if the receiver were not offset from the source. Therefore, the time interval between successive multiple reflection quantities or events in the idealized reproduction of a seismic recording according to FIGS. 3 and 4 is not constant, but it generally increases with the later recording time. It is also clear that any change in the strength of the multi-reflective position or the inclination of the boundaries over a distance segment containing the reflection points leads to a further deviation from a constant time interval between successive multi-reflection quantities or events.
In a practical seismic situation, step in

m Verbindung mit dem Amplitudenverhalten aufeinanderfolgender Mehrfachreflexionen zusätzliche Komplikationen auf. Es ist allgemein bekannt, daß der Reflexionskoeffizient an einer seismischen Diskontinuität von den akustischen Impedanzen, die an der Grenzem connection with the amplitude behavior of successive Multiple reflections raise additional complications. It is well known that the Reflection coefficient at a seismic discontinuity from the acoustic impedances that are at the limit

,", auftreten, sowie von den Einfallswinkeln und der Übertragung seismischer Wellenfronten an der Grenze abhängt. Anhand Fig. I und 2 ergibt sich, daß diese Winkel von dem Abstand der Versetzung zwischen der seismischen Quelle und dem Empfänger, der Ordnung, ", occur, as well as from the angles of incidence and the Transmission of seismic wave fronts depends on the boundary. Based on Fig. I and 2 it can be seen that this Angle of the distance of the displacement between the seismic source and the receiver, the order

.'■'' Ul.! lfl\.IIM«V.III\.ll^Aliyir JUtVIV. Uli KLIgUIIg UV.I MkIIt* xionsflächcn abhängen. In einer praktischen Situation kann nicht davon ausgegangen werden, daß die akustischen Eigenschaften reflektierender Grenzen an jedem Punkt die gleichen sind. Diese Faktoren führen zu . '■''Ul.! lfl \ .IIM «V.III \ .ll ^ Aliyir JUtVIV. Uli KLIgUIIg UV.I MkIIt * xionsflächencn depend. In a practical situation it cannot be assumed that the acoustic properties of reflective boundaries are the same at every point. These factors lead to

.»> Unterschieden bei den Reflexionskoeffizienten, die für die seismische Welle an verschiedenen Reflexionsstellen tatsächlich maßgeblich sind. Ebenso ist es allgemein bekannt, daß seismische Wellenfrontflächen eine etwa kugelförmige Gestalt haben und daher infolge Diver-. »> Differences in the reflection coefficients used for the seismic wave at different reflection points are actually decisive. It is also general known that seismic wavefront surfaces have an approximately spherical shape and therefore due to diver-

in genz in ihrer Amplitude abfallen. Der Erfinder der vorliegenden Anmeldung, hat in Verbindung mit der vorerwähnten Veröffentlichung D gezeigt, daß der Amplitudenabfall infolge Divergenz durch Brechungseffekte beeinflußt wird und es im allgemeinen nichtdecrease in amplitude. The inventor of the present application, has shown in connection with the aforementioned publication D that the The decrease in amplitude due to divergence is influenced by refraction effects and it is generally not

ts möglich ist, eine einzige Korrektur für einen Amplitudenabfall infolge Divergenz anzugeben, die sich auf alle möglichen Reflexionsereignisse zu allen Aufzeichnungszeiten anwenden ließe.
Aus den vorstehenden Überlegungen läßt sichts is possible to specify a single correction for a drop in amplitude due to divergence, which can be applied to all possible reflection events at all recording times.
From the above considerations it can be

to ersehen, daß die Verfahren nach dem Stand der Technik zur Dämpfung oder Unterdrückung von Mehrfachreflexionen, die eine ebene Wellenausbreitung, ein normales Auftreffen und einen horizontalen Lagenverlauf voraussetzen oder inverse Filter-Operatoren aus einer to see that the methods according to the prior art for attenuating or suppressing multiple reflections, which require plane wave propagation, normal impingement and a horizontal position profile, or inverse filter operators from a

4ϊ statistischen Auswertung des Verhaltens seismischer Signalkurven ableiten, im günstigsten Fall nur teilweise zum Erfolg führea Es ist daher eine weitere Verfeinerung notwendig, um geringfügige, aber wichtige Änderungen in der Amplitude und der zeitlichen Festlegung von Mehrfachreflexions-Größen aufzunehmen, wie sie in einer praktischen Situation aufgezeichnet werden. 4ϊ Derive statistical evaluation of the behavior of seismic signal curves, in the best case only partially lead to success a Further refinement is therefore necessary in order to include minor but important changes in the amplitude and the time-dependent definition of multiple reflection quantities, as in a practical situation recorded .

Wie zuvor erwähnt besteht ein erster Schritt bei der Anwendung der vorliegenden Erfindung darin, eine anfängliche Schätzung des Mehrfachreflexions-SignaJs zu erhalten, das in der seismischen Signalkurve enthalten ist Bei der weiter ins einzelne gehenden Beschreibung der Erfindung wird davon ausgegangen, daß die zu verarbeitenden seismischen Signalkurven mit As previously mentioned, a first step in practicing the present invention is to obtain an initial estimate of the multiple reflection signal contained in the seismic waveform. In further describing the invention, it is assumed that the seismic signals to be processed Signal curves with

bo regelmäßigen diskreten Zeitintervallen ausgetastet und in digitaler Form aufgezeichnet worden sind, wie das derzeit üblich ist Ein bevorzugtes Verfahren, das Mehrfachreflexionssignal zu schätzen, arbeitet demnach mit einem Wiener-Voraussagefilter in analoger Weise, wie das von Peacock und Treitel in der oben erwähnten Veröffentlichung C beschrieben wird, wo der Voraussage-Abstand etwas kleiner als die kürzeste interessierende Mehrfachreflexionsdauer ist Wenn das Voraussage bo regular discrete time intervals have been blanked and recorded in digital form, as is currently customary. A preferred method of estimating the multiple reflection signal accordingly works with a Wiener prediction filter in an analogous manner to that of Peacock and Treitel in the above-mentioned publication C is described where the prediction distance is slightly less than the shortest multiple reflection duration of interest if that prediction

IOIO

Filter mit einer seismischen Signalkurve beaufschlagt und das Ausgangssigna! um eine Zeitspanne verzögert wird, die gleich der Voraussageentfernung ist, so entspricht das Ergebnis einer Schätzung des voraussagbaren Inhalts der Eingangs-Signalkurve, d. h. der Mehrfachreflexionsereignisse. Die Signalkurve enthält nach Faltung mittels des Voraussage-Filters keine primären Reflexionsgrößen, weil das Auftreten dieser Ereignisse im allgemeinen nicht voraussagbar ist.A seismic signal curve is applied to the filter and the output signal! delayed by a period of time which is equal to the prediction distance, the result corresponds to an estimate of the predictable Content of the input signal curve, d. H. of multiple reflection events. The signal curve contains after convolution by means of the prediction filter, no primary reflection quantities, because the occurrence of these Events are generally unpredictable.

Ebenso sind weitere Möglichkeiten für eine anfängliche Schätzung der Mehrfachreflexionssignale in Verbindung mit der Anwendung der vorliegenden Erfindung möglich. Beispielsweise stellt das Verfahren /ur Synthetisierung von Mehrfachreflexionen durch Retrokorrelation seismischer Signalkurven, wie es von Anstey und Newman in der britischen Patentschrift 11 95 204 bzw. der US-PS 34 96 529 beschrieben wird, ein direktes Verfahren zur Abschätzung dar. Das Verfahren voraussagender Deconvolution kann eine Schätzung bert-Transformation gebildet und die quadrierten Werte des ursprünglichen Signals und seiner Quadraturfunktion summiert werden.Likewise, there are other options for an initial Estimation of the multiple reflection signals in connection with the practice of the present invention possible. For example, the method / ur synthesizes multiple reflections by retrocorrelation seismic waveforms as described by Anstey and Newman in British Patent 11 95 204 or US Pat. No. 3,496,529 describes a direct method for estimation. The method predictive deconvolution can be an estimate formed bert transform and squared Values of the original signal and its quadrature function are summed.

Es sind weitere Verfahren zur Festlegung von Punkten auf der Zeitachse, wo die Mehrfachschätzung in Segmente unterteilt werden soll, möglich, wobei diese weiteren Verfahren mit unter die Erfindung fallen sollen. Die Einhüllende kann beispielsweise durch Bildung eines laufenden Mittelwertes des Signalspurbetrags angenähert werden, oder die quadrierte Einhüllende kann durch laufende Mittelwertbildung der Werte der quadrierten Signalspur angenähert werden. Wesentlich ist, daß Punkte auf der Zeitachse bestimmt werden müssen, an denen die Energie in der der Mehrfach-Schätzung entsprechenden Signalspur ein Minimum ist.There are other methods of establishing points on the time axis where the multiple estimate is to be divided into segments, possible, these further methods also falling under the invention should. The envelope can be created, for example, by forming a running mean value of the signal trace amount can be approximated, or the squared envelope can be calculated by averaging the values can be approximated to the squared signal trace. It is essential that points are determined on the time axis at which the energy in the signal trace corresponding to the multiple estimate is a minimum.

Ein Vergleich der Diagramme 6a und 6bzeigt, daß die Mehrzahl der im Verhältnis zueinander ausgerichteten Echos im Diagramm 6a ebenso im Diagramm 6b wiedergegeben ist, so daß diese somit als Mehrfachre-A comparison of diagrams 6a and 6b shows that the majority of the echoes aligned in relation to one another in diagram 6a are also reproduced in diagram 6b , so that they are thus represented as multiple re-

VUIl IVICIII ICtUIICIIVUIl IVICIII ICtUIICII

viii; tn.uiviii; tn.ui

bevorzugten Verfahren äquivalent ist, wenn die einer Deconvolution ausgesetzten Signalkurven von den Eingangskurven subtrahiert werden. Ebenso können Backus-Filter-Operatoren eine vernünftige Schätzung von Mehrfachreflexionen liefern, wenn der anfängliche Einheitsimpuls aus dem Operator ausgeschlossen wird und die verbleibenden Ausdrücke hinsichtlich ihrer Polarität invertiert werden. Allgemein lassen sich viele Mehrfach-Unterdrückungsverfahren darauf abstimmen, eine gewisse Schätzung von Mehrfachreflexionen zu liefern, indem einfach der Ausgangswert des Verfahrens vom Eingangswert subtrahiert wird. Die Zuverlässigkeit der so erhaltenen Vorausschätzung hängt vom Erfolg des Verfahrens ab.preferred method is equivalent if the deconvolution-exposed signal curves from the Input curves are subtracted. Likewise, Backus filter operators can give a reasonable estimate of multiple reflections if the initial unit pulse is excluded from the operator and invert the polarity of the remaining terms. In general, many can Multiple suppression methods adjust to a certain estimate of multiple reflections by simply subtracting the output value of the method from the input value. The reliability the forecast obtained in this way depends on the success of the procedure.

In Fig. 6 zeigt das Diagramm 6a eine seismische Original-Aufzeichnung, wie sie im Verlauf einer auf See durchgeführten Untersuchung in einer Wassertiefe von etwa 180 m gemacht wurde. Das Entfernungsintervall zwischen Empfängerstellen betrug 50 m, und die Quelle war gegenüber dem nächsten Empfänger um 300 m versetzt. Das Diagramm 6b zeigt die erste Abschätzung des Gehalts an Mehrfachreflexionen, sog. Multiplen, der gleichen Aufzeichnung, wie sie unter Verwendung eines Wiener-Voraussage-Filters erhalten wurde. Die Aufzeichnung 6c wurde durch Subtraktion der ersten Schätzung der Multiplen von der Originalaufzeichnung erhalten. Sie entspricht dem Ergebnis eines bekannten Voraussage-Dekonvolutionsverfahrens. Das Diagramm 6<y zeigt das Ergebnis, wie es durch das Verfahren nach der vorliegenden Erfindung erhalten wurde, wobei die erste Schätzung der Multiplen in der richtigen Weise für die Subtraktion von der Eingangsaufzeichnung ausgerichtet wurde.In FIG. 6, the diagram 6a shows an original seismic recording as it was made in the course of an investigation carried out at sea in a water depth of about 180 m. The distance interval between receiver locations was 50 m and the source was offset from the nearest receiver by 300 m. Diagram 6b shows the first estimate of the content of multiple reflections, so-called multiples, of the same recording as was obtained using a Wiener prediction filter. Record 6c was obtained by subtracting the first estimate of the multiples from the original record. It corresponds to the result of a known predictive deconvolution method. The diagram 6 <y shows the result as obtained by the method according to the present invention with the first estimate of the multiples properly aligned for subtraction from the input record.

Ds die erforderliche Anpassung nach Zeit und Amplitude sowohl für jede Signalkurve als auch für die einzelnen Multiplen unterschiedlich sind, ist es notwendig, auf der Zeitachse jeder Signalspur die Punkte zu bestimmen, an denen Änderungen der Anpassungsfaktoren am besten vorgenommen werden. Entsprechend dem Verfahren nach der vorliegenden Erfindung werden diese Punkte dorthin gelegt, wo der quadrierte Betrag der Einhüllenden der den geschätzten Multiplen entsprechenden Signalspur Minimalwerte annimmt F i g. 7 zeigt beispielhaft die quadrierte Hüllkurve, die der ersten Signalspur der Aufzeichnung 66 entspricht Die quadrierte Hüllkurve wird bei dieser bevorzugten Betriebsart bestimmt, indem die Quadransrfunktionen des geschätzten Mehrfachsignals durch diskrete HiI-zeigt jedoch, daß die Amplitude und die zeitliche Festlegung der geschätzten Multiplen nicht völlig mit den tatsächlichen Verhältnissen übereinstimmt. Fehler in der zeitlichen Festlegung werden besonders anhand von Größen deutlich, die zu Zeitpunkten von etwa 1,0 s oder weniger in der Signalspur auftreten, die der Aufzeichnung entspricht, die der Quelle am nächsten ist. Beispielsweise erscheint die Größe, die in dieser Spur bei etwa 0,8 s auftritt, in der Vorausschätzung annähernd 12 ms zu spät. Infolge dieser Fehler werden die frühen Mehrfachreflexionen im Diagramm 6c, das die von der Eingangsaufzeichnung subtrahierte geschätzte Mehrfachaufzeichnung wiedergibt, nicht nennenswert gedämpft.Ds the necessary adjustment according to time and Amplitude are different for each signal curve as well as for the individual multiples, it is necessary to on the time axis of each signal trace to determine the points at which changes in the adjustment factors best to be done. According to the method of the present invention these points are placed where the squared The amount of the envelope of the signal trace corresponding to the estimated multiples assumes minimum values F i g. 7 shows an example of the squared envelope curve that corresponds to the first signal trace of recording 66 In this preferred operating mode, the squared envelope is determined by the quadransr functions of the estimated multiple signal by discrete HiI- shows that the amplitude and the temporal Establishing the estimated multiples does not entirely correspond to the actual circumstances. failure in the time specification are particularly clear on the basis of variables that at points in time of around 1.0 s or fewer occur in the signal track that corresponds to the recording that is closest to the source. For example, the magnitude that occurs in this trace at around 0.8 s appears in the forecast approximately 12 ms too late. As a result of these errors the early multiple reflections in graph 6c, which is the estimated one subtracted from the input record Multiple recording reproduces, not significantly attenuated.

Um die Zeitkorrekturen zu bestimmen, der die Schätzung der Mehrfachen ausgesetzt werden muß, wird jede geschätzte Signalspur durch Bezugnahme auf die Signalkurven-Hüllfunktion in der zuvor beschriebenen Weise in Segmente unterteilt. Jedes Segment wird dann mit einem entsprechenden Teil der Eingangs-Signalkurve kreuzkorreliert, um die Zeitverschiebung zu bestimmen, die aufgebracht werden muß. um die beste Anpassung zwischen geschätzten Mehrfachen und tatsächlichen Mehrfachen zu erhalten. Die erforderliche Zeitverschiebung wird so gewählt, daß sich eine Maximierung der Kreuzkorrelationsfunktion ergibt. Der Maximalwert jeder Kreuzkorrelationsfunktion wird außerdem für die Verwendung bei der Bestimmung der Amplitudenmaßstabskorrektur gespeichert, wie sie ι für jedes Segment der Mehrfachschätzung erforderlich ist. Da der Bereich wahrscheinlicher Zeitfehler üblicherweise vorausbestimmt werden kann, lassen sich Begrenzungen vornehmen, um die Anzahl von Verzögerungen in der Kreuzkorrelationsfunktion zu beschränkea To determine the time corrections to which the multiple estimate must be subjected, becomes each estimated signal trace by referring to the signal curve envelope function in that previously described Way divided into segments. Each segment is then associated with a corresponding part of the input signal curve cross-correlated to determine the time shift that needs to be applied. to be the best Adjustment between estimated multiples and actual multiples. The required The time shift is chosen so that the cross-correlation function is maximized. The maximum value of each cross-correlation function is also used in the determination the amplitude scale correction is stored as required for each segment of the multiple estimate is. Since the range of probable timing errors can usually be predetermined, Set limits to limit the number of delays in the cross-correlation function a

Der für jedes Segment der Mehrfachschätzung erforderliche Amplitudenkorrekturfaktor ergibt sich annähernd ausThe amplitude correction factor required for each segment of the multiple estimation results approximately off

wobei Φ_ηΛ der zuvor bestimmte Kreuzkorrelations-Scheitelwert und Φ*_Λ der Nullverzögerucgs-Autokorrelationswert des entsprechenden Schätzungssegments ist. Um dies besser zu verstehen, sei berücksichtigt, daß das in einem bestimmten Segmentwhere Φ _{ηΛ is} the previously determined cross-correlation peak value and Φ * _{Λ is} the zero delay autocorrelation value of the corresponding estimation segment. To understand this better, consider that this is in a specific segment

enthaltene geschätzte Mehrfachsignal durch eine Folge von Austastwerten wiedergegeben wird:contained estimated multiple signal by a sequence is reproduced by blanking values:

m₀, Ot₁, m₂ ... m„ m ₀ , Ot ₁ , m ₂ ... m "

und daß das tatsächliche Mchrfachsignal, da₃ in einem entsprechenden Segment der Eingangssignalkurve enthalten ist, durch eine weitere Folge repräsentiert wird:and that the actual multiple signal, since _{3 is contained} in a corresponding segment of the input signal curve, is represented by a further sequence:

mn, nt\, mn, nt \,

... m„ ... m "

Die Eingangssignalkurve kann außerdem ein primäres Reflexionssignal enthalten, das durch eine weitere Folge repräsentiert wird:The input signal curve can also contain a primary reflection signal that is passed through another Sequence is represented:

O_n. 0\, Ü2 ... a„. O _n . 0 \, Ü2 ... a ".

Es sei weiter angenommen, daß die Schätzung der Mehrfachen sich jetzt in der richtigen zeitlichen ZunrHnnno befindet ^c<***h ahor KA7iinlirh jhrpr Amnlitiulp von den tatsächlichen Mehrfachen noch um einen MaßstabsL-ktor unterscheidet, so daßIt is further assumed that the estimate of the multiples now located in the right temporal ZunrHnnno ^c <*** h Ahor KA7iinlirh jhrpr Amnlitiulp different ctor MaßstabsL-from the actual multiple nor a so

km, =■■ m,km, = ■■ m,
und somitand thus

Φ™* ⁼ — Σ \a,m, + m,m,\ Φ ™ * ⁼ - Σ \ a, m, + m, m, \
ηη

= —- Σ [a,m, + m*]= —- Σ [a, m, + m *] knkn

- m: - m:

k^:nk ^: n

m,.m ,.

DaherTherefore istis frfr ηη ι ~^~ ^i j ι ~ ^ ~ ^ i j yy mimi \m \ m Φ,Φ,

(2)(2)

Die Beziehung (2) zeigt, daß der aus der Näherungsbeziehung (1) gewonnene Amplitudenkorrekturfaktor k genau wird, wenn kein primäres ReP.exionssignal in einem speziellen Segment der Eingangs-Signalkurve enthalten ist. Der Fehler in der Näherung, der sich aus dem Ausdruck ^anm in der Gleichung (2) ergibt, ist im allgemeinen klein, da er von der Zufälligkeits-Korrelation zwischen primären Reflexionssignalen und Mehrfachreflexionssignaicn abhängt, die sich aus früheren primären Größen ergeben. Das Vorzeichen des Fehler-Ausdrucks kann mit gleicher Wahrscheinlichkeit positiv oder negativ sein, und die Größe verliert an Bedeutung, wenn die Amplitude von Mehrfachsignalen hn Verhältnis zu primären Signalen zunimmt So wird das Verfahren in solchen Fällen effektiver, bei denen Mehrfachreflexionen einen starken Einfluß haben. Um eine Absicherung dagegen vorzunehmen, daß übermäßige Fehler auftreten, kann in einfacher Weise der Wert von k, der tatsächlich als Korrekturfaktor angewandt wird, begrenzt werden, so daß der wahrscheinliche Bereich der Mehrfachschätzungs-Fehler nicht überschritten wird. The relationship (2) shows that the amplitude correction factor k obtained from the approximate relationship (1) becomes accurate if no primary reflection signal is contained in a specific segment of the input signal curve. The error in the approximation that results from the expression ^ anm in equation (2) is generally small, since it depends on the randomness correlation between primary reflection signals and multiple reflection signals resulting from earlier primary quantities. The sign of the error term can be positive or negative with equal probability, and the size becomes less important as the amplitude of multiple signals increases in relation to primary signals.This makes the method more effective in cases where multiple reflections have a strong influence. In order to safeguard against the occurrence of excessive errors, the value of k which is actually used as a correction factor can be limited in a simple manner so that the probable range of the multiple estimate errors is not exceeded.

Zeit- und Amplitudenkorrekturen, wie sie oben bestimmt wurden, werden für jedes Segment der den geschätzten Mehrfachen entsprechenden Signalspuren vorgenommen, und diese Signalspuren werden anschließend von den Eingangs-Signalspuren subtrahiert, um eine Ausgangsaufzeichnung entsprechend Diagramm 6c/zu erzeugen. Ersichtlich sind in dieser Mehrfachreflexionssignale erheblich mehr gedämpft bzw. ganz unterdrückt.Time and amplitude corrections like those above are determined for each segment of the signal traces corresponding to the estimated multiples and these signal traces are then subtracted from the input signal traces to to generate an output record according to diagram 6c /. This shows multiple reflection signals considerably more attenuated or completely suppressed.

Wie zu ersehen, muß für eine perfekte Verwirk·', chung der vorliegenden Erfindung angestrebt werden, die Wellenform einer jeden speziellen Multiplen möglichst wiedergabegetreu zu schätzen und die Zeitverschiebungen und die Maßstabsfaktoren möglichst genau zu bestimmen. Wenn Übereinstimmung sichergestellt werden könnte, so wurden Mehrfa···'- Icxionssignale vollständig aus der seismischen Auf/-.),' nung ausgeschlossen. In der Pr.ixis reicht es aus, lediglich Mehrfachreflexionssignale zu dämpfen, so daß "rimäre Reflexionen nicht üb?r^Qch?*tt?t wprHen so daß nur eine nicht ganz so perfekte Verwirklichung benötigt wird.As can be seen, for a perfect implementation of the present invention, one must strive to estimate the waveform of each particular multiple as faithfully as possible and to determine the time shifts and the scaling factors as accurately as possible. If agreement could be ensured so r Multi ··· were '- Icxionssignale completely from the seismic On / -).' Planning impossible. In the Prixis it is sufficient to only attenuate multiple reflection signals so that "peripheral reflections do not sound over ^Q ch? * Tt? T so that only a less than perfect realization is required.

Mit Fig. 8 sind die verschiedenen Schritte schematisch veranschaulicht, wie sie bei einer bevorzugten Anwendung der Erfindung auftreten. Die einzelnen Schnitte bzw. die zu ihrer Durchführung dienenden Elemente und ihre Arbeitsweise sind mit den Blöcken 8Λ bis SG gezeigt. Der Einfachheit halber werden die Operationen in der Form digitaler Computerausdriicke beschrieben, jedoch können einer oder mehrere Schritte auch durch analoge oder sonstige Mittel ausgeführt werden, wie das in der vorstehenden Beschreibung angedeutet wurde.FIG. 8 schematically illustrates the various steps as they occur in a preferred application of the invention. The individual cuts or the elements used to carry them out and their mode of operation are shown with blocks 8Λ to SG . For convenience, the operations will be described in the form of digital computer expressions, but one or more steps may also be carried out by analog or other means as indicated in the foregoing description.

In F i g. 8 wird das Signal 51 der Eingangssignalkurve zu einem dem Block 8.4 entsprechenden Voraussageoperatorbild sowie mit dem Block SB einer Voraussageoperatoranwendung zugeführt, wobei der Block 8fl gleichzeitig auch vom Block SA gespeist wird. Der Ausgang des Blocks 85 gelangt zu einem eine Zeitverschiebung erzeugenden Block SC. In Fig. 8, the signal 51 of the input signal curve is fed to a prediction operator image corresponding to block 8.4 and, with block SB, to a prediction operator application, block 8fl also being fed from block SA at the same time. The output of block 85 goes to a block SC which generates a time shift.

Der Ausgang des Blocks SC ist die erste Schätzung des Mehrfachreflexionssignals und gelangt zu einem Block 8D, der einem Zeitgeber für die Segmente einer quadrierten Hüllkurve entspricht und für eine Unterteilung in Zeitintervalle sorgt. Das Ausgangssignal des Blocks SD gelangt zu einem Block SE, der zum Zeit- und Amplitudenvergleich dient. Der Block SE wird außerdem mit dem Eingangssignal Si sowie dem Ausgangssignal des Blocks 8Cgespeist. Ein Block 8Fsorgt für eine Korrektur der Schätzung, und sein Ausgang liefert die segmentartigen Schätzungen, die von den Segmenten des Signals S] mittels eines Blocks SG subtrahiert werden, der eine Summierfunktion ausübt und die verarbeitete Ausgangs»Signalkurve S₂ liefertThe output of the block SC is the first estimate of the multiple reflection signal and arrives at a block 8D, which corresponds to a timer for the segments of a squared envelope and provides for a division into time intervals. The output signal of the block SD arrives at a block SE, which is used to compare the time and amplitude. The block SE is also fed with the input signal Si and the output signal of the block 8C. A block 8F provides for a correction of the estimate, and its output supplies the segment-like estimates which are subtracted from the segments of the signal S] by means of a block SG which performs a summing function and supplies the processed output signal curve S ₂

In Verbindung mit dem ersten bevorzugten Verfahren löst Block SA die matrixartig aufgebaute Gleichung:In connection with the first preferred method, block SA solves the matrix-like equation:

r„-\ T_n-r "- \ T _n -

/ο/ ο ₌₌ r.-2r.-2 /ι/ ι II. ιι II. ιι II. II. roro

ga+\ga + \

ga + n-1 ga + n-1

worin γ, die Autokorrelation der Eingangs-Signalkurve und als Annäherung der Autokorrelation einer einzel-where γ, the autocorrelation of the input signal curve and as an approximation of the autocorrelation of a single

nen primären seismischen Welle zusammen mit einer zugehörigen Folge mehrfachreflektierte'. Wellen angesehen wird, gi entspricht der Kreuzkorrelation der gewünschten Ausgangswelle und der Eingangsweile, während f· der Voraussage-Operator ist, der durch Lösung der Gleichung (3) bestimmt werden muß. Gleichung (3) ist eine spezielle Version der sog. Wiener-Normalgleichungen. In diesem besonderen Fall ist der gewünschte Ausgang einfach eine zeitlich vorgerückte Version der Eingangswelle, und die Kreuzkorrelation g, ist das Äquivalent der Autokorrelation r_b die längs der Zeitachse um einen Betrag gleich dem Voraussage-Abstand α verschoben ist, d. h. g_: wird numerisch gleich /·,+*. Wie ersichtlich, stellt Gleichung (3) eine Reihe von π simultanen Gleichungen mit π Unbekannten dar, so daß sie sich leicht lösen läßt.A primary seismic wave along with an associated sequence of multiple reflections'. Waves is considered, gi corresponds to the cross-correlation of the desired output wave and the input wave, while f · is the prediction operator which must be determined by solving equation (3). Equation (3) is a special version of the so-called Wiener normal equations. In this particular case, the desired output is simply a temporally advanced version of the input wave, and the cross-correlation g, is the equivalent of the autocorrelation r _b shifted along the time axis by an amount equal to the prediction distance α, i.e. g _: becomes numerically the same / ·, + *. As can be seen, equation (3) represents a series of π simultaneous equations with π unknowns so that they can be easily solved.

Im Block 8J5 wird der Voraussage-Operator /, in Verbindung mit der Eingangs-Signalkurve x, gebracht, um so eine erste Schätzung des voraussagbaren Teils der Signalkurve zu liefern, d. h. das mehrfachreflektierte Signal, das sich in voraussagbarer Weise nach dem Auftreten jedes primären Reflexionssignals wiederholt Der damit verbundene Dekonvolutionsvorgang wird für diskrete ausgetastete Daten durch die nachstehende Beziehung bestimmt:In block 8J5, the prediction operator / is connected to the input signal curve x so as to provide a first estimate of the predictable part of the signal curve, ie the multiply reflected signal which is predictable after the occurrence of each primary reflection signal repeated The associated deconvolution process is determined for discrete blanked data by the following relationship:

(4)(4)

Realteil zu bildea Die Quadraturkomponente des analytischen Signals wird im Block SD bestimmt, indem die Mehrfachschätzung einer [Convolution mit dem diskreten Hilhert-Transformutions-Operator ausgesetzt wird:Real part to form a The quadrature component of the analytical signal is determined in block SD by subjecting the multiple estimate to a convolution with the discrete Hilhert transformation operator:

v, orin Xi der diskrete ausgetastete Operand, f, der Operator und y_T das Produkt der Konvolution sind.v, orin Xi is the discrete blanked operand, f is the operator and y _{T is} the product of the convolution.

Wenn der Operator /, im vorliegenden Fall ein Voraussage-Operator ist, der den Wert von x, zu einer bestimmten zukünftigen Zeit t+tx vorwegnehmen soll, so liegt die erste Schätzung von Mehrfach-Reflexionssignalen um eine Zeit gleich <x vor den tatsächlichen Signalen. Daher wird durch den Block SC eine diesem Wert entsprechende zeitliche Verzögerung eingeführt Am Ausgang des Blocks SC steht die Mehrfachschätzung nunmehr zum Vergleich mit tatsächlich aufgezeichneten Signalen zur Verfügung, um Korrekturen zu bestimmen, die zur Verbesserung der Schätzung notwendig sein können.If the operator /, in the present case a prediction operator, is intended to anticipate the value of x at a certain future time t + tx , then the first estimate of multiple reflection signals is before the actual signals by a time equal to <x . A time delay corresponding to this value is therefore introduced by the block SC . At the output of the block SC , the multiple estimation is now available for comparison with actually recorded signals in order to determine corrections which may be necessary to improve the estimation.

In der Praxis ist zu erwarten, daß Fehler in der Schätzung von Mehrfach-Reflexionssignalen für verschiedene Punkte der Zeitachse schwanken, so daß Korrekturen der Schätzung strikt zeitveränderlich sein sollten. In der vorliegenden bevorzugten Version wird eine Art Quasi-Zeitveränderung dadurch erzielt, daß die Schätzung der Mehrfachen in Segmente unterteilt und für jedes Segment gesondert Korrekturwerte bestimmt werden. Dementsprechend bestimmt der Block SD die quadrierte Hüllkurve vor der der Mehrfachschätzung entsprechenden Signalkurve. Die Hüllkurve ist als der Betrag (Modul) des analytischen Signals bestimmt, von dem die Mehrfachschätzung hergeleitet wird, um den , 2 Sin² T"
Λ, = — · In practice it is to be expected that errors in the estimation of multiple reflection signals will vary for different points on the time axis, so that corrections of the estimation should be strictly time-varying. In the present preferred version, a kind of quasi-time change is achieved in that the estimate of the multiple is divided into segments and correction values are determined separately for each segment. Accordingly, the block SD determines the squared envelope curve in front of the signal curve corresponding to the multiple estimate. The envelope is determined as the amount (module) of the analytical signal from which the multiple estimate is derived, by which "2 Sin ² T"
Λ, = - ·

ίΦΟίΦΟ

0,0,

Es wird dafür gesorgt, daß der Block 8 D den Operator h, durch Beaufschlagung mit einem Hamming-Fenster oder einer anderen geeigneten Funktion ein- '^ engt, um unerwünschte Abrundungseffekte auf einem Minimum zu halten. Die Quadraturkomponente des analytischen Signals ist dann:It is ensured that the block D 8 h the operator, by applying a Hamming window or other suitable function ^ mono- concentrated to hold unwanted rounding effects to a minimum. The quadrature component of the analytic signal is then:

z, = m,* ht, z, = m, * ht ,

worin das Stemzeichen die Konvolution andeutet und z, vollständig in der imaginären Ebene enthalten ist Die quadrierte Hüllkurve der Schätzung der Mehrfach-Reflexionssignale ist dann:where the stem symbol indicates the convolution and z, is completely contained in the imaginary plane.The squared envelope of the estimate of the multiple reflection signals is then:

Ein typisches Beispiel für die quadrierte Hülikurve isiA typical example of the squared envelope curve isi

jo mit Fig.7 wiedergegeben. Ein weiterer Zweck des Blocks SD besteht darin, die Zeitpunkte zu erfassen, zu denen wesentliche Minima m der quadrierten Hüllkurvenfunktion auftreten, und um die Schätzung der Mehrfach-Reflexionssignale zu diesen Zeitpunkten injo reproduced with Fig.7. Another purpose of the block SD is to record the times at which significant minima m of the squared envelope function occur and to estimate the multiple reflection signals at these times in

j-> Segmente zu unterteilen.j-> divide segments.

Der Block 8ffkorreliert Segmente der Schätzung der Mehrfach-Reflexionssignale fn, mit den entsprechender Zeitsegmenten der aufgezeichneten Signalkurve x_t Zeitkorrekturen, denen die einzelnen aufeinanderfolgenden Segmente der Schätzung ausgesetzt v/erder müssen, werden von den Zeitverschiebungen abgeleitet die zur Maximierung jeder Korrelation erfordernd sind. Der Block SE bestimmt außerdem den Nullver zögerungs-Autokorrelationswert jedes Segments dei Abschätzung, und er sorgt ferner für eine Berechnung des Verhältnisses der Kreuzkorrelations-Scheitc lwertc zu den Autokorrelations-Scheitelwerten, um den in dei Beziehung (1) der Beschreibung enthaltenen Maßstabs korrekturwert K zu liefern.The block 8ff correlates segments of the estimation of the multiple reflection signals fn with the corresponding time segments of the recorded signal curve x _t Time corrections to which the individual successive segments of the estimation must be subjected are derived from the time shifts which are required to maximize each correlation. The block SE also determines the zero delay autocorrelation value of each segment of the estimate, and also provides for a calculation of the ratio of the cross-correlation peaks to the autocorrelation peaks in order to give the correction value K contained in relationship (1) of the description deliver.

V) Korrekturen hinsidfitßch der zeitlichen Lage und dei Amplitude des geschätzten Mehrfach-Reflexionssignal· werden im Block SF vorgenommen, und die korrigiert« Schätzung wird schließlich von der Eingangs-Signalkur ve im Element SG subtrahiert V) Corrections with regard to the temporal position and the amplitude of the estimated multiple reflection signal are made in block SF , and the corrected estimate is finally subtracted from the input signal curve in element SG

Ilicr/ti 4 UIaII Zcichmir^cnIlicr / ti 4 UIaII Zcichmir ^ cn

Claims (1)

Patentansprüche:Patent claims: 1. Verfahren zur weitmöglichsten Beseitigung der Multiplen aus einer oder mehreren grafisch oder phonographisch reproduzierbaren, insbesondere digital abgetasteten und gespeicherten seismischen Datenspuren seismischer Wellen, die von einem Sender ausgehend mehrere vom Sender und voneinander einen Abstand aufweisenden Empfängern direkt und/oder nach mehrfacher Reflexion an horizontalen unterhalb der Oberflächen liegenden Schichtenerreichen, dadurch gekennzeichnet, daß jede einzelne der Datenspuren nach Reproduktion folgendermaßen behandelt wird:1. Procedure for the greatest possible elimination of the Multiples of one or more graphically or phonographically reproducible, especially digital sampled and stored seismic data traces of seismic waves generated by a Transmitter based on a plurality of receivers that are spaced apart from the transmitter and from one another directly and / or after multiple reflections on horizontal ones below the surfaces Reaching layers, characterized by that each of the data tracks is treated as follows after reproduction: