DE2323983A1 - PROCESS FOR GENERATING A DIGITAL SEGMENT FILE FROM SEISMIC COLLECTIVE RECORDS FOR POINTS OF COMMON DEPTH AND ARRANGEMENT FOR PERFORMING THE PROCESS - Google Patents

Description

Verfahren zur Erzeugung einer digitalen Segrcentdatei aus seismischen Sammelaufzeichnungen für Punkte gemeinsamer Tiefe und Anordnung zur Durchführung des VerfahrensMethod for generating a digital Segrcent file from seismic collective records for common points Depth and arrangement for carrying out the procedure

Die Erfindung bezieht sich auf ein Verfahren zum Interpretieren seismischer Aufschlußdaten und insbesondere auf ein Verfahren und eine Anordnung air Erzeugung digitaler Segmentdateien aus Sammelaufzeichnungen für Punkte gemeinsamer Tiefe.The invention relates to a method for interpreting seismic outcrop data and, more particularly, to a Method and arrangement air generation of digital segment files from collective records for points of common depth.

Seismische Explorationsverfahren unter Verwendung von Punkten gemeinsamer Tiefe (CDP-Explorationsverfahren) werden bei der Suche nach Erdöl-Lagerstätten allgemein angewendet. Bei der CDP-Exploration werden mehrere seismische Spuren aufgezeichnet, die durch eine Tiefenschichtredundanz gekennzeichnet sind-, was sich aus der Geometrie der gemeinsamen Tiefenschicht ergibt, bei der Gruppen von seismischen Signalen von identischenSeismic Exploration Techniques Using Points Of Common Depth (CDP Exploration Techniques) will be used in the The search for oil deposits is generally used. During the CDP exploration, several seismic traces are recorded, which are characterized by a depth layer redundancy - which results from the geometry of the common depth layer, at the groups of seismic signals of identical

-~j
Tiefenschichtpunkten reflektiert werden, jedoch zu unterschiedlichen Schußpunkt- und Empfängerorten gehören. Die CDP-Exploration ergibt infolge der Tiefenschichtredundanz besondere Vorteile bei der Störungsverminderung u.dgl. Eine typische CDP-Explöration ist in der US-PS 2 732 906 beschrieben. - ~ j
Depth layer points are reflected, but belong to different shot point and receiver locations. CDP exploration has particular advantages in reducing interference, etc. due to the redundancy of the depth layers, etc. A typical CDP exploration is described in US Pat. No. 2,732,906.

309847/0941309847/0941

OFUGlNAL INSPECTBDOFUGlNAL INSPECTBD

Bisher sind verschiedene Verfahren entwickelt worden, bei denen ein programmierter digitaler Computer zum Sichten und Herausgreifen von Tiefenschichten-Reflexionsinformationen aus CDP-Daten verwendet wird. Ein frühes, mit einem Computer arbeitendes Dichtungs- und Entnahmeverfahren arbeitete mit einer relativ einfachen eindimensionalen Festste llungs- und Spurverfolgungslogik zum Verfolgen von Amplitudenextremwerten längs eines gestapelten Aufzeichnungs- Schnitts, fast so, als würde ein Auswerter Horizonte mit einem Farbstift heraussuchen. In neuerer Zeit sind zweidimensionale Sichtungsalgorithmen entwickelt worden, die Amplitudenextremwerte entweder in den Zeit-Moveout-Bereichen oder in den Zeit-Dip-Bereichen feststellen, so daß eine Reflexionsinformation und eine verbesserte Sichtungswirkung erzielt v/erden.So far, various methods have been developed in which a programmed digital computer to Sifting and picking out deep layer reflection information from CDP data is used. An early The computerized sealing and extraction process used a relatively simple one-dimensional one Determination and tracking logic for tracking extreme amplitude values along a stacked recording slice, almost as if an evaluator were picking out horizons with a colored pencil. In more recent times are two-dimensional Sighting algorithms have been developed which have extreme amplitude values in either the time moveout ranges or in the time dip ranges so that reflection information and an improved sighting effect can be achieved.

Bei der Zeit-Moveout-Sichtungstechnik wird die CDP-Sammelaufzeichnung über einen Bereich von Moveout-Werten durch Zeitverschieben und Summieren längs hyperbolischer Bahnen abgetastet. Der resultierende Satz Λτοη summierten Ausgangsspuren , wobei für jeden unterschiedlichen Moveout-Wert eine Ausgangsspur vorhanden ist, bildet einen zweidimensional en Suchraum, in dem Amplitudenextremwerte liegen. Der Ort jedes Amplitudenextremwerts ergibt eine gemeinsame Schätzung der Laufzeit, der Amplitude und des Moveouts des Ereignisses. Die Extremwerte werden automatisch zwischen benachbarten Tiefenpunkten verglichen, und wenn die verglichenen Extremwerte innerhalb einer gegebenen Toleranzgrenze in den gewünschten Größen liegen, dann werden sie zur Bildung von Reflektorsegmenten aneinandergehängt. Das oben beschriebene Verfahren wird mit Hilfe des "600 Processing Package" ausgeführt,das von Geophysical Services Inc., einer Tochtergesellschaft der Texas Instruments Inc., Dallas, Texas, vertrieben wird. "600 Processing Package" ist in einer Broschüre mit dem Titel "600 -Processing Package Examples" , vom Juli 1970In the time moveout viewing technique, the CDP composite record is scanned over a range of moveout values by time shifting and summing along hyperbolic trajectories. The resulting set Λτοη summed output tracks, with an output track being available for each different moveout value, forms a two-dimensional search space in which extreme amplitude values lie. The location of each amplitude extreme gives a common estimate of the time of flight, amplitude, and moveout of the event. The extreme values are automatically compared between adjacent depth points, and if the compared extreme values are within a given tolerance limit in the desired sizes, then they are linked to one another to form reflector segments. The method described above is carried out using the "600 Processing Package" sold by Geophysical Services Inc., a subsidiary of Texas Instruments Inc., Dallas, Texas. "600 Processing Package" is in a booklet titled "600-Processing Package Examples", dated July 1970

309847/0941309847/0941

und außerdem in dem Aufsatz "Developments in Seismic Data Processing and Analysis ( 1968 bis 1970)" von William. A. Schneider in der Zeitschrift Geophysics, XXXVI, Nor. 6, Seiten 1043 bis 1073 vom Dezember 1971 beschrieben.and also in the article "Developments in Seismic Data Processing and Analysis (1968 to 1970) "by William. A. Schneider in the journal Geophysics, XXXVI, Nor. 6, Pages 1043-1073 of December 1971.

Die oben erwähnten zweidimensionalen Sichtungsverfahren haben zwar brauchbare Reflexionsinformationen ergeben, doch waren sie oft nicht dazu geeignet, Daten mit schlechter Qualität in befriedigender V/eise zu sichten. Solche bisher bekannte zweidimensionale Sichtungsverfahren ergaben somit oft Segmentanzeigen, in denen in erkennbarer Weise in Stör- und Rauschbereichen Daten fehlten. Daher ist ein Bedarf nach einem verbesserten Reflexionssegment-Sichtungsverfahren entstanden, bei dem auch eine exakte Sichtung durch Störungen hindurch und bei niedrigeren Rauschabständen ausgeführt v/erden kann, als es früher möglich war.While the two-dimensional sighting methods mentioned above have yielded useful reflection information, they were they are often not suited to sifting through poor quality data in a satisfactory manner. Such previously known Two-dimensional inspection processes thus often resulted in segment displays in which in a recognizable manner in disturbance and noise areas data were missing. Therefore, there is a need for an improved reflective segment sighting method arose, in which an exact sighting through disturbances and at lower signal-to-noise ratios is also carried out can than was previously possible.

Mit Hilfe der Erfindung wird ein Verfahren zur Erzeugung einer Segmentdatei aus digitalen Sammelaufzeichnungen seismischer Spuren mit Tiefenschichtredundanz geschaffen. Bei dem Verfahren werden die seismischen Spuren gestapelt, während unterschiedliche Moveout- und Dip-¥erte auf die seismischen Spuren angewendet werden. Innerhalb der gestapelten seismischen Spuren werden dann Amplitudenextremwerte festgestellt, damit Schätzungen der Laufzeit-, Amplituden-, Moveout- und , Dip-Werte des Ereignisses erzeugt werden. Jede der Schätzungen wird mit benachbarten Schätzungen verglichen, und es werden Reflektorsegmente gebildet, die sich zwischen benachbarten Schätzungen ähnlicher Werte erstrecken.With the aid of the invention, a method for generating a segment file from digital collective records becomes seismic Tracks created with depth layer redundancy. In the process, the seismic traces are stacked while different moveout and dip ¥ ¥ ¥ ¥ values can be applied to the seismic traces. Inside the stacked seismic Traces are then determined amplitude extremes, so that estimates of the transit time, amplitude, moveout and, Dip values of the event are generated. Each of the estimates is compared with neighboring estimates, and it reflector segments are formed that extend between adjacent estimates of similar values.

Gemäß einer Ausgestaltung der Erfindung wird eine Anordnung geschaffen, mit deren Hilfe eine Segmentdatei aus digitalen Sammelaufzeichnungen seismischer Spuren mit Tiefenschichtredundanz erzeug werden kann. Die Anordnung enthält Einrichtungen zum Programmieren eines digitalen ComputersAccording to one embodiment of the invention, an arrangement is created with the aid of which a segment file from digital Collective recordings of seismic traces with depth layer redundancy can be generated. The arrangement contains facilities for programming a digital computer

309847/0941309847/0941

derart, daß die seismischen Spuren gestapelt werden, während auf sie unterschiedliche Moveout- und Dip-Werte angewendet werden. Der Computer ist auch so programmiert, daß Amplitudenexi;remwerte innerhalb der gestapelten seismischen Spuren gesichtet werden, damit Schätzungen der Laufzeit-, Ampli-=- tuden-, Moveout-und Dip-Werte des Ereignisses erzeugt werden. Der Computer bestimmt dann Reflektorsegmente, die sich zwischen benachbarten Schätzungen mit ähnlichen Werten erstrecken.such that the seismic traces are stacked while having different moveout and dip values applied to them will. The computer is also programmed so that amplitude exits sighted within the stacked seismic traces to allow estimates of transit time, ampli - = - tuden, moveout and dip values of the event are generated. The computer then determines reflector segments that extend between neighboring estimates with similar values.

Gemäß einer Weiterbildung der Erfindung wird ein Verfahren zur Erzeugung einer Segmentdatei geschaffen, bei dem Sammelaufzeichnungen für Punkte gemeinsamer Tiefe während der Anwendung unterschiedlicher Moveout-und Dip-Werte zur Erzeugung eines dreidimensionalen Suchvolumens gestapelt worden sind. Das Verfahren umfaßt die sequentielle Anwendung eines vorbestimmten Testbereichs längs eines vorbestimmten Zeitkorridors des Suchvolumens. Während jeder Anwendung des Testbereichs wird der innerhalb des Testbereichs liegende höchste Amplitudenwert festgestellt. Für jeden der höchsten Amplitudenwerte werden dann die Zeit- und Moveout-Werte festgestellt. Dann wird der Dip-Wert jedes höchsten Amplitudenwerts festgestellt, und aus den die höchsten Amplitudenwerte definierenden Parametern wird eine Datei aufgebaut.According to a development of the invention, a method for generating a segment file is created in which Collective records for points of common depth while applying different moveout and dip values have been stacked to create a three-dimensional search volume. The method includes the sequential Application of a predetermined test area along a predetermined time corridor of the search volume. During each application of the test range, the highest amplitude value within the test range is determined. The time and moveout values are then determined for each of the highest amplitude values. Then the dip value of each highest amplitude value is determined, and out of those the highest A file is created for parameters defining amplitude values.

Gemäß einer weiteren Ausgestaltung der Erfindung wird ein Verfahren zur Erzeugung einer Segmentdatei aus einer digitalen Sammelaufzeichnung seismischer Spuren für Punkte gemeinsamer Tiefe geschaffen. Die digitale Sammelaufzeichnung wird für mehrere unterschiedliche Moveout-Werte zeitlich verschoben. Die zeitlich verschobene Sammelaufzeichnung wird dann für jeden der unterschied-According to a further embodiment of the invention, a method for generating a segment file is from created a collective digital record of seismic traces for points of common depth. The digital Collective recording is shifted in time for several different moveout values. The time-shifted Collective recording is then made for each of the different

309847/094 1309847/094 1

lichen Moveout-Werte summiert, so daß mehrere Moveout Abtastspursammelaufzeichnungen bei verschiedenen Tiefenpunkten erzeugt werden. Die Moveout-Abtastspursammelauf zeichnungen werden dann entsprechend unterschiedlichen Dip-Richtungen verschoben. Die Moveout-Abtastsammelspuraufzeichnungen v/erden dann für jede der unterschiedlichen Dip-Richtungen summiert, damit ein dreidimensionales Suchvolumen mit den Koordinaten Zeit, Moveout und Dip erzeugt wird. Innerhalb dieses dreidimensionalen Suchvolumens werden dann Amplitudenextremwerte festgestellt, damit Schätzungen für die Laufzeit-, Amplituden-, Moveout-und Dip-Werte des Ereignieses erzeugt werden. Die Schätzungen werden mit Schätzungen von Extremwerten an benachbarten Tiefenpunkten verglichen. Dann werden Reflektorsegmente gebildet, die sich zwischen benachbarten Extremwerten erstrecken, die innerhalb vorbestimmter Tbleranzgrenzen liegende Schätzungen aufweisen.lichen Moveout values are summed, so that several Moveout Collective tracking recordings at different depth points be generated. The moveout scan track accumulation records then become correspondingly different Dip directions shifted. The moveout scan heap records v / earths then summed up for each of the different dip directions, thus creating a three-dimensional search volume with the coordinates time, moveout and dip. Be within this three-dimensional search volume Then extreme amplitude values are determined to make estimates for the transit time, amplitude, moveout and dip values of the Events can be generated. The estimates are made with estimates of extreme values at neighboring depth points compared. Then reflector segments are formed that extend between adjacent extreme values that are within Have estimates lying within predetermined tolerance limits.

Ein Ausführungsbeispiel der Erfindung ist in der Zeichnung dargestellt.Darin zeigen:An embodiment of the invention is shown in the drawing.

Fig.1 eine schematische Darstellung eines typischen Explorationsverfahrens mit Punkten gemeinsamer Tiefe,1 is a schematic representation of a typical exploration process with points of common depth,

Fig.2 eine schematische Darstellung einer Sammelaufzeichnung für Punkte gemeinsamer Tiefe,2 shows a schematic representation of a collective recording for points of common depth,

Fig.3 ein Flußdiagramm der Hauptschritte des erfindungsgemäßen Verfahrens,Figure 3 is a flow chart of the main steps of the invention Procedure,

Fig.4 ein Funktionsflußdiagramm der zeitabhängigen Filterung und der Sperrentzerung des erfindungsgemäßen Verfahrens,4 shows a functional flow diagram of the time-dependent filtering and the blocking equalization of the method according to the invention,

309847/0941309847/0941

Fig.5 ein Funktionsflußdiagramm der bei der Erfindung angewendeten Parameterinitialisierung,Fig. 5 is a functional flow diagram of that used in the invention Parameter initialization,

Fig.6 eine schematische Darstellung der vom Benutzer eingegebenen Zeit-und Geschwindigkeitsgrenzen für den Betrieb des erfindungsgemäßen Verfahrens,6 shows a schematic representation of the entered by the user Time and speed limits for the operation of the method according to the invention,

Fig.7 eine schematische Darstellung der durch Berechnen des Moveout-Abtastmaßes bestimmten Ausgangsnadelimpulse,7 shows a schematic representation of the by calculating the Moveout scanning distance determined output needle impulses,

Fig.8a einen zeitabhängig gefilterten Wellenzug, der mit der Moveout-Bereichs-Impulsantwort gefaltet ist,Fig.8a shows a time-dependent filtered wave train, which with the moveout area impulse response is convoluted,

Fig.8b den zeitverschobenen, zeitabhängig gefilterten Wellenzug, der sich aus der Faltung des in Fig.8a dargestellten Wellenzugs ergibt,8b the time-shifted, time-dependent filtered Wave train, which results from the folding of the wave train shown in Fig. 8a,

Fig.9 einen zeitabhängig gefilterten Wellenzug mit Einheitsamplitude , 9 shows a time-dependent filtered wave train with a unit amplitude,

Fig.10 einen zeitabhängig gefilterten Wellenzug mit der Amplitude 0,55 als Folge einer Reihe von Faltungen,Fig. 10 a time-dependent filtered wave train with the Amplitude 0.55 as a result of a series of folds,

Fig.11 eine Reihe von im Abstand voneinander liegenden Nadelimpulsen, die um das Dip-Abtastmaß als Ergebnis der zu dem Wellenzug von Fig.10 führenden Faltungen voneinander getrennt sind,Fig. 11 a series of spaced needle pulses, the one by the dip sampling dimension as a result of the convolutions from one another leading to the wave train of FIG are separated

Fig. 12 eine schenatische Darstellung der Definition der - Geschwindigkeitsgrenzen,Fig. 12 is a schematic representation of the definition of - speed limits,

Fig. 13 eine Transformation der Grenzen von Fig.12 auf Moveout-Grenzen, FIG. 13 shows a transformation of the limits of FIG. 12 to moveout limits,

309847/0941309847/0941

Fig.i4 ein Funktionsflußdiagramm für die Berechnung der Moveout-Abtastgrenzwerttabelle,Fig.i4 a function flow diagram for the calculation of the Moveout sample limit table,

Fig.15 ein Diagramm zur Darstellung der Sichtungstechnik nach der Erfindung ohne Beteiligung der GesamtheitFig. 15 is a diagram showing the visualization technique according to the invention without involvement of the whole

aller Daten,
Fig.i6 eine schematische Darstellung des Vorgangs zum Aufbau und zur Interpolation eines zweidimensionalen Filters,all data,
Fig.i6 is a schematic representation of the process for building and interpolating a two-dimensional filter,

Fig.17 ein Funktionsflußdiagramm für die Ausführung des Moveout-Abtastaufbauverfahrens nach der Erfindung,Fig.17 is a functional flow diagram for carrying out the Moveout scanning construction method according to the invention,

Fig.18 eine typische Moveout-Abtastsammelaufzeichnungsspur,Fig. 18 shows a typical moveout scan collective recording track,

Fig.19 ein Funktionsflußdiagramm für die Durchführung des Dip-Abtastaufbauverfahrens nach der Erfindung,19 is a functional flow diagram for carrying out the Dip-scanning construction method according to the invention,

Fig.20 eine schematische Darstellung einer Gruppe von Moveout-Abtastungen bei den Tiefenpunkten j— ,Fig.20 is a schematic representation of a group of Moveout samples at depth points j—,

Fig.21 eine schematische Darstellung eines dreidimensionalen Suchvolumens, das von Mo^ceout-Abtastungen definiert ist, die in unterschiedlichen Dip-Richtungen nach der Erfindung gestapelt sind,Fig.21 is a schematic representation of a three-dimensional Search volume defined by Mo ^ ceout scans is stacked in different dip directions according to the invention,

Fig.22 ein Funktionsflußdiagramm für eine dreidimensionale Ereignissichtung nach der Erfindung,FIG. 22 is a functional flow diagram for a three-dimensional Event viewing according to the invention,

Fig.23 pine graphische Darstellung eines Testbereichs nach der Erfindung,Fig. 23 pine graphical representation of a test area according to the invention,

Fig.24 eine graphische Darstellung, in dem ein gewünschtes Maximum mit Hilfe eines zuvor ausgeführtenAmplitudensuchvorgangs liegen soll,Fig. 24 is a graph showing a desired Maximum should be with the help of a previously performed amplitude search process,

309847/0941309847/0941

Fig.25 eine graphische Darstellung von Datenpunkten und interpolierten Punkten, die mit Hilfe einer interpolierenden Filterung erzeugt werden,Fig. 25 is a graphical representation of data points and interpolated points that are generated with the help of interpolating filtering,

Fig.26 ein verkleinertes Interpolationsfeld,Fig. 26 shows a reduced interpolation field,

Fig.27 ein weiter verkleinertes Interpolationsfeld, "bei dem der absolute Wert eines Datenpunkts A größer als der absolute Wert eines Datenpunkts E ist,FIG. 27 a further reduced interpolation field, "in which the absolute value of a data point A is greater than the absolute value of a data point E,

Fig.28 ein weiteres verkleinertes Interpolationsfeld, bei dem der absolute Wert des . Datenpunkts A kleiner als der absolute Wert des Datenpunkts Ξ ist,28 shows a further reduced interpolation field in which the absolute value of the. Data point A is smaller than the absolute value of data point Ξ,

Fig.29 ein als Beispiel zu interpolierendes künstliches Datenfeld, 29 shows an artificial data field to be interpolated as an example,

Fig.30 eine Darstellung zur Veranschaulichung der Interpolation des in Fig.29 dargestellten Datenfeldes entsprechend der bei der Erfindung angewendeten Logik mit einem Logikbaum,30 shows a representation to illustrate the interpolation of the data field shown in Fig. 29 according to the logic used in the invention with a logic tree,

Fig.31 eine graphische Darstellung des Amplitudenextremwertsichtens in drei benachbarten Dip-Ebenen,Fig. 31 is a graphic representation of the amplitude extremes viewing in three neighboring dip levels,

Fig.32 eine graphische Darstellung der Interpolation der Amplituden- und Dip-Werte für.die Amplitudenextrem werte zwischen den in Fig.31 dargestellten benachbarten Dip-Ebenen,Fig. 32 is a graphic representation of the interpolation of the Amplitude and dip values for the extreme amplitude values between the neighboring ones shown in FIG Dip levels,

Fig.33 ein Funktionsflußdiagramm des Verfahrens nach der Erfindung zum Aufbau von Segmenten, undFIG. 33 shows a functional flow diagram of the method according to FIG Invention for building segments, and

Fig.34 ein Diagramm zur Veranschaulichung des Aufbaus eines Reflektorsegments zwischen zwei benachbarten Fenstern nach der Erfindung.Fig. 34 is a diagram to illustrate the structure of a Reflector segment between two adjacent windows according to the invention.

309847/09/4 1309847/09/4 1

Fig.1 zeigt schematisch ein mit gemeinsamen Tiefenpunkten (CDP) arbeitendes Explorationsverfahren. Eine genauere Beschreibung einer solchen CDP-Exploration findet sich in der US-PS 2 732 906. Schußpunkte 10, 12, 14 und 16 werden nacheinander längs der seismischen Trasse zur Explosion gebracht. Längs der seismischen Trasse sind Geophone 18 in gleichem Abstand voneinander angebracht, und sie werden gemeinsam mit der Explosion der Schußpunkte längs dieser Trasse bewegt. Die von den Geophonen 18 erzeugten Signale werden in einem herkömmlichen Aufzeichnungsgerät 20 verstärkt und aufgezeichnet.1 shows schematically an exploration method working with common depth points (CDP). A more accurate one A description of such a CDP exploration can be found in U.S. Patent 2,732,906. Shots 10, 12, 14 and 16 are detonated one after the other along the seismic route. Are along the seismic route Geophones 18 are placed at the same distance from each other, and they are shared with the explosion of the points of fire moved along this route. The signals generated by the geophones 18 are recorded in a conventional recording device 20 amplified and recorded.

Es ist zu erkennen, daß während des Betriebs des dargestellten Systems vom Schußpunkt 10 ausgehende Energie an einem Punkt 22 auf einem reflektierenden Horizont 24 reflektiert und beim Geophon 26 empfangen wird. Außerdem wird vom Schußpunkt 12 ausgehende Energie vom Punkt 22 reflektiert und vom Geophon 28 empfangen. Somit bildet der auf der vertikalen Linie 29 liegende Punkt 22 einen gemeinsamen Reflexionstiefenpunkt. Wenn die Schußpunkte und die Geophone längs der seismischen Trasse bewegt werden, werden somit viele solcher Punkte gemeinsamer Tiefe erzeugt. Die sich bei der Exploration mit Punkten gemeinsamer Tiefe ergebende Tiefenschichtredundanz ist für die anschließende digitale Verarbeitung von Vorteil, wie unten noch beschrieben wird. Im Aufzeichnungsgerät 20 wird ein herkömmliches Seismogramm mit zackenförmigen Aufzeichnungsspuren (wiggle trace seismogram) aufgezeichnet. Das Seismogramm gibt beim Zeitpunkt Null den Explosionszeitpunkt an, und es zeigt ferner die Signallaufzeit zu einem Reflexionspunkt und zur Erdoberfläche zurück in Längsrichtung der Aufzeichnung an. Bei bisher bekannten Verfahren wurden die einzelnen Aufzeichnungsspuren des Seismogramms digitalisiert und so gruppiert, It can be seen that energy emanating from the point of fire 10 during the operation of the system shown a point 22 on a reflective horizon 24 and received at the geophone 26. aside from that The energy emanating from the shot point 12 is reflected by the point 22 and received by the geophone 28. Thus forms the point 22 lying on the vertical line 29 has a common reflection depth point. When the shot points and the geophones are moved along the seismic route, thus many such points become more common Creates depth. The depth layer redundancy resulting from exploration with points of common depth is advantageous for subsequent digital processing, as will be described below. In the recorder 20 a conventional seismogram is recorded with jagged recording tracks (wiggle trace seismogram). The seismogram indicates the time of the explosion at time zero, and it also shows the signal propagation time to a reflection point and back to the surface of the earth in the longitudinal direction of the record. At so far known methods, the individual traces of the seismogram were digitized and grouped so

309847/0941309847/0941

daß Aufzeichnungsspuren von Punkten gemeinsamer Tiefe zusammen angeordnet sind. Diese Gruppierung wird als Sammelaufzeichnung für einen Punkt gemeinsamer Tiefe bezeichnet; eine typische Aufzeichnung nach der Anfangsverarbeitung ist in Fig.2 dargestellt. Fig.2 zeigt eine herkömmliche 24-fache Sammelaufzeichnung längs eines Zeitabschnitts von 6 Sekunden. Bei der Anwendung der hier beschriebenen Erfindung v/erden solche Sammelaufzeichnungen in digitalisierter Form auf Magnetband gespeichert und vor der Verarbeitung in den Plattenspeicher eines Computers eingegeben. that recording tracks of points of common depth are arranged together. This grouping is called Denotes collective record for a point of common depth; a typical record after initial processing is shown in Fig.2. Fig.2 shows a conventional one 24-fold collective recording over a period of time of 6 seconds. When using the invention described herein, such collective records are grounded stored in digitized form on magnetic tape and entered into the disk memory of a computer before processing.

Das hier beschriebene Verfahren kann mit jedem geeigneten, entsprechend programmierten digitalen Computer ausgeführt werden. In der bevorzugten Ausführung wird das Verfahren hinsichtlich einer Verarbeitung mit einen Computer 870A TIAC der Firma Texas Instruments , Inc. Dallas," Texas beschrieben. Der Computer 870A ist im Programmiererhandbuch für das TIAC Model 870A der Texas Instruments, Inc. 1968 beschrieben.The method described here can be carried out with any suitable, appropriately programmed digital computer will. In the preferred embodiment, the method is in terms of processing with a computer 870A TIAC of Texas Instruments, Inc. Dallas, "Texas. The 870A computer is in the Programmer's Guide for the TIAC Model 870A from Texas Instruments, Inc. in 1968.

Fig.3 zeigt die Hauptfunktionsschritte des hier zu beschreibenden Verfahrens. Die digitalen Samine lauf zeichnungen für Punkte gleicher Tiefe (CDP-Sammelaufzeichnungen) werden vom Magnetband bei 30 in den Computerspeicher eingegeben. Eingangsparameter werden bei 32 durch den Benutzer eingegeben. Die später im einzelnen genauer beschriebenen Parameter werden vom Benutzer entsprechend der vorher bereits vorhandenen Kenntnis der Daten und entsprechend den gewünschten Ergebnissen bestimmt. Die CDP-Sammelaufzeichnungen werden bei 34 konditioniert sowie einer zeitveraMerlichen Filterung und einer Spurentzerrung unterzogen. Die für die anschließende Verarbeitung benötigtenFig. 3 shows the main functional steps of the one to be described here Procedure. The digital samples are run for points of the same depth (CDP collective records) entered from the magnetic tape at 30 in the computer memory. Input parameters are entered at 32 by the user. The parameters described in more detail later are set by the user in accordance with the previously pre-existing knowledge of the data and determined according to the desired results. The CDP collective records are conditioned at 34 as well as time-dependent filtering and track equalization subjected. The ones required for subsequent processing

309847/0941309847/0941

Parameter werden "bei 35 initialisiert. Bei 36 wird jedes Zeitfenster (time gate) der CDP-Sammelaufzeichnungen durch Zeitverschiebung und Summierung längs hyperbolischer Bahnen über einen Moveout-Bereich abgetastet. Mit "moveout¹¹ wird die Zunahme der Laufzeit bedingt durch die Meßgeometrie bezeichnet. Die Gruppen der summierten Ausgangsspuren für jeweils einen unterschiedlichen Moveout-Wert bilden einen zweidimensionalen Snchraum, in dem Amplitudenextremwerte aufgefunden werden können. Die summierten Ausgangsspuren werden dann bei 38 über einen Bereich von Dip-Werten durch Zeitverschiebung und anschließendes Stapeln abgetastet. Mit "Dip" wird die Neigung eines reflektierenden Horizonts bezeichnet. Dadurch entsteht ein dreidimensionales Suchvolumen, das eine Gruppe von summierten Ausgangsspuren mit jeweils einem unterschiedlichen Moveout- und Dip-Parameter enthält. Dieses dreidimensionale Suchvolumen wird dann bei 40 einer dreidimensionalen Ereignissichtung unterzogen, so daß in dem Volumen Amplitudenextremwerte zur Erzielung eines simultanen Schätzwerts für die Laufzeit, die Amplitude, den Moveout-Wert und den Dip-Wert des Ereignisses festgestellt werden.Parameters are "initialized" at 35. At 36, each time window (time gate) of the CDP collective recordings is scanned by time shifting and summing along hyperbolic paths over a moveout area. "Moveout ¹¹ " denotes the increase in transit time due to the measurement geometry. The groups of the summed output tracks for a different moveout value each form a two-dimensional snapshot space in which extreme amplitude values can be found. The summed output traces are then sampled at 38 over a range of dip values by time shifting and then stacking. "Dip" is the slope of a reflective horizon. This creates a three-dimensional search volume that contains a group of summed output tracks, each with a different moveout and dip parameter. This three-dimensional search volume is then subjected to a three-dimensional event viewing at 40, so that extreme amplitude values are determined in the volume in order to obtain a simultaneous estimate for the transit time, the amplitude, the moveout value and the dip value of the event.

Die simultanen Schätzwerte werden bei 42 mit benachbarten Schätzwerten verglichen, und wenn die Schätzwerte in allen Dimensionen innerhalb von Toleranzwerten liegen, die der Benutzer angegeben hat, dann werden die Extremwerte einander logisch zugeordnet oder "aneinandergehängt", damit Reflektorsegmente gebildet werden. Die angesammelten Segmente werden dann aufgezeichnet und bei 44 mit einer konstanten Rauaabtastgeschwindigkelt versehen. Die Segmentdaten können dann bei 46 dazu benutzt werden, eine große Anzahl von wichtigen Anzeigen zu erzeugen. Beispiele solcher Anzeigen sind Gesamtsegmentdateien, Primärsegmentdateien, nichtprimäre Segmentdateien, Geschwindigkeitsstreudiagramme für Gesamtsegmentdateien oder für ausgewählte Tiefenpunkte,The simultaneous estimates are compared to neighboring estimates at 42, and if the estimates are in all Dimensions are within tolerance values specified by the user, then the extreme values will be logically assigned to one another or "linked" so that reflector segments are formed. The accumulated Segments are then recorded and given a constant rough scan rate at 44. The segment data can then be used at 46 to generate a large number of important displays. Examples of such Displays are total segment files, primary segment files, non-primary segment files, velocity scatter graphs for total segment files or for selected depth points,

309847/09309847/09

Segmentklassen, Primärsegmente für Spitzen oder Täler, Intervallgeschv/indigkeitsanzeigen u.dgl« Beispiele solcher Anzeigen finden sich in der Druckschrift mit dem Titel "600 Package Processing Examples" vom Juli 1970, die von Geophysical Service Inc., einer Tochtergesellschaft von Texas Instruments Inc., Dallas, Texas herausgegeben worden ist.Segment classes, primary segments for peaks or valleys, interval speed displays and the like. Examples of such displays can be found in the publication entitled "600 Package Processing Examples "dated July 1970, provided by Geophysical Service Inc., a subsidiary of Texas Instruments Inc., Dallas, Texas.

Datenkonditicnierun^i zeitabhängige Filterung und Sntζerrang. Data condition nierun ^ i time-dependent filtering and signal rank .

Die Figuren 4 bis 17 zeigen verschiedene Aspekte der Datenkonditionierungj der zeitabhängigen Filterung lind der Spurentzerrungj wie sie im Schritt 3-^ ausgeführt "werden, Nach Fig,4 werden die CDF-Sa^melaufzeichnungen bei 50 eingegeben, und sie sind S-, 12~, 24- oder mehrfachFigures 4 to 17 show various aspects of the Datenkonditionierungj the time-dependent filtering lind the Spurentzerrungj be as performed in step 3 ^ "As shown in FIG, 4, the CDF-Sa ^ melaufzeichnungen be entered at 50, and they are S-, ~ 12 , 24 or more times

XiC UJJJr^-1C-i-'vti i*3jl· ^^"Sg'i^i -^j. _ c-.'α .,!U.i i-Si _,--.' «·;■,,. να :;-_t.....,^,.,. i »:"XiC UJJJr ^-1 Ci-'vti i * 3jl · ^^ "Sg'i ^ i - ^ j. _ C -. 'Α.,! Ui i-Si _, -.'«·; ■ ,,. Να:; - _t ....., ^,.,. I»: "

filtert. Bei diesem Piltsrungs^rgans; xint :1a:"' Bsnutzsr vier Frecuenzpunkte ui:d dia ^Q^lsiBoi^i ΐ-ίΐ-ΐ^-ΐκη^δ sin. Das Filter wird dann mit der Sinsc.ngaspur s;er£l~-f--ti Insbesondere wird die Spur ^rulT-^pli^isr-ti s'xiinii^rt i:nd dann zeitverschoben, damit dia gemiDschtsn ^ejiiiltorten Ausgangspunkte erzeugt verderu Ti;ir,s derartige llgitala Filterung ist bereits entwickelt -worden, und ,=i3 ist zusammen mit der von Geophysical SerYicSj "nc, _} einer Tochtergesellschaft von Texas Instruments, Ins*, Dallas, Texas auf den Harkt gebrachten ^πβ00 Package" angewendet worden. Nach der Filterung wird die Spur bei 56 entzerrt. Der Benutzer gibt eine Startzeit (SQ) und eine Endzeit (EQ) des Entzerrungsfensters in Blöcken ein. Die Potenz jedes Zeitfensters ergibt sich dann axis der folgenden Gleichung :filters. With this Piltsrungs ^ rgans; xint: 1a: "'Bsnutzsr four frequency points ui: d dia ^ Q ^ lsiBoi ^ i ΐ-ίΐ-ΐ ^ -ΐκη ^ δ sin. The filter is then s; er £ l ~ -f-- ti In particular, the track ^ rulT- ^ pli ^ isr-ti s'xiinii ^ rt i: nd is then time-shifted so that the mixed starting points are generated i3 has been used in conjunction with the ^π β00 Package "brought to the market by Geophysical SerYicSj" nc, _} a subsidiary of Texas Instruments, Ins *, Dallas, Texas. After filtering, the track is equalized at 56. The user gives a start time . (SQ) and an end time (EQ) of the equalization window blocks in a the potency of each time slot is then given axis of the following equation:

309847/09^ 1309847/09 ^ 1

wobei gilt: EQ = Ende des Entzerrungsfensters in Blöcken SQ = Start des Entzerrungsfensters in Blöcken A, = Amplitude Jeder Spur.where: EQ = end of the rectification window in blocks SQ = start of the rectification window in blocks A, = amplitude of each track.

Die Spuren werden dann durch Multiplikation mit einem Faktor S^ entzerrt, der folgendermaßen definiert ist:The tracks are then equalized by multiplication with a factor S ^ which is defined as follows:

S-j ="V|r;J =1, 2, ..., NFOLD (2)S-j = "V | r; J = 1, 2, ..., NFOLD (2)

wobei X eine beliebige quadratische Konstante ist.where X is any quadratic constant.

Nachdem jede Spur bei 56 entzerrt worden ist, springt der Vorgang zum Herausgreifen einer weiteren einzelnen CDP-Spur bei 52 zurück, und dieser Vorgang wird ständig wiederholt. Die zeitabhängig gefilterten und entzerrten Spuren werden dann bei 58 im Plattenspeicher gespeichert.After each track has been equalized at 56, the jumps The process of picking out another single CDP track at 52, and this process is repeated over and over. The time-filtered and equalized tracks are then stored at 58 in disk storage.

Bei 60 wird die mittlere Potenz F folgendermaßen bestimmt: ₁ NFOLDAt 60, the mean power F is determined as follows: ₁ NFOLD

NFOLD i-rx ^P,NFOLD i-rx ^P ,

J- I JJ- I J

Dieser V/ert P" wird dann für die Verwendung in der anschließenden maßstäblichen Änderung jeder Spur während der Moveout-Abtasterzeugung bei 36 (Fig.3) gespeichert. V/enn es erwünscht ist, kann jede Spur aber auch auf den Uert P" zu diesem Zeitpunkt entzerrt werden. Während der Moveout-Abtasterzeugung bei 36 wird ein korrigierender Faktor S'rf folgendermaßen gebildet:This V / ert P "will then be used in the subsequent scale change of each track during moveout scan generation is stored at 36 (Figure 3). V / enn it is desired, but each track can also be on the Uert P " be equalized at this point. During the moveout scan generation at 36, a corrective factor is applied S'rf formed as follows:

2, ..., NFOLD (4)2, ..., NFOLD (4)

Jede der Spuren wird mit dem korrigierenden Faktor S¹ . multipliziert, damit jede Spur hinsichtlich des Unterschieds zwischen X und P korrigiert wird.Each of the tracks is given the correcting factor S ¹ . multiplied so that each track is corrected for the difference between X and P.

309847/094 1309847/094 1

Der korrigierende Faktor S¹^ wiTd im Computer-Trommelspeicher bei 62 zusammen mit zusätzlichen, vom Tiefenpunkt abhängigen Informationen gespeichert, die bei 64 vom 25. Kanal des die CDP-Sammelaufzeichnungen enthaltenden Magnetbands eingegeben werden. Solche vom Tiefenpunkt abhängige Informationen enthalten beispielsweise die Wassertiefe, die Empfängerabstände, Störungen (statics) und dgl. The corrective factor S ¹ ^ wiTd is stored in the computer drum memory at 62 along with additional depth point dependent information entered at 64 from the 25th channel of the magnetic tape containing the CDP composite recordings. Such information, which is dependent on the depth point, contains, for example, the water depth, the receiver distances, disturbances (statics) and the like.

Parameterinitialisierung tParameter initialization t

Fig.5 zeigt ein Funktionsflußdiagramm für die Initialisierung von Zeitfensterparametern, die für die anschließende Verarbeitung benötigt werden« Der Benutzer gibt Eingangsparameter wie Start- und Endzeiten hei 70 ein, und die gefilterten und entzerrten CDP-Spuren werden bei 72 aus dem Speicher eingegeben. Die Bezugszeit T_Q der Moveout-Kurve wird bei 64 nach der folgenden Gleichung berechnet:5 shows a functional flow diagram for the initialization of time window parameters required for subsequent processing. The user enters input parameters such as start and end times at 70, and the filtered and equalized CDP tracks are entered from memory at 72. The reference time T _{Q of} the moveout curve is calculated at 64 according to the following equation:

τ - T(J-H) + TCi) ,^ τ - T (JH) + TCi) , ^

wobei gilt: T(i) = vom Benutzer eingegebenewhere: T (i) = entered by the user

Startzeit,Start time,

T(i+1) = vom Benutzer eingegebene Endzeit.T (i + 1) = user entered end time.

Die maximalen und minimalen Moveout-Werte, die bei der Suche angewendet werden sollen, werden bei 76 berechnet. Der Benutzer gibt einen durch Knicke und Schenkel nach Fig.6 bestimmten Geschwindigkeitsbereich ein. Unter Anwendung des vom Benutzer eingegebenen Zeitfensters und definiert durch den Bereich zwischen T^ und Tg werden die minimalen und maximalen Geschwindigkeiten als V^ und V₂ The maximum and minimum moveout values to be used in the search are calculated at 76. The user enters a speed range determined by kinks and legs as shown in FIG. Using the time window entered by the user and defined by the range between T ^ and Tg, the minimum and maximum speeds are expressed as V ^ and V ₂

309847/0941309847/0941

bestimmt. Zur Bestimmung der Minimal-und Maximalgeschwindigkeit aus den in Fig.6 angegebenen Grenzen werden die Indices der Geschwindigkeitswerte bestimmt, die die Startzeit (j) und die Endzeit (k) des Zeitfensters begrenzen (k > j). Zur Bestimmung der zu suchenden Minimal- und MaximalgeEchwindigkeiten werden die folgenden Gleichungen angewendet:certainly. To determine the minimum and maximum speed from the limits indicated in FIG. 6, the indices of the speed values are determined which limit the start time (j) and the end time (k) of the time window (k> j). The following equations are used to determine the minimum and maximum speeds to search for applied:

(6a)(6a)

W ^VH<^k> ⁺ < )) (TCi+D- TVCk)) (6b)W ^VH < ^k > ⁺ <)) (TCi + D- TVCk)) (6b)

wobei gilt: VL(j) = niedrige Geschwindigkeit der Schenkel, where: VL (j) = low speed of the legs,

mit j = 1, 2, ... , Zahl der Schenkel + 1, with j = 1, 2, ..., number of legs + 1,

T(i) = Startzeit eines Zeitfensters, TV(i) = Zeit des Geschwindigkeitsioiicks, VH(k) = hohe Geschv/indigkeit des SchenkelsT (i) = start time of a time window, TV (i) = time of the speed limit, VH (k) = high speed of the thigh

mit k = 1, ... Zahl der Schenkel + 1 , T(i41)= Endzeit des Zeitfensters.with k = 1, ... number of legs + 1 , T (i41) = end time of the time window.

Aus den minimalen und maximalen Geschwindigkeitswerten werden die für den Suchvorgang benötigten minimalen und maximalen Mouveout-Werte aus der folgenden Gleichung berechnet: The minimum and maximum mouveout values required for the search process are calculated from the following equation from the minimum and maximum speed values:

wobei gilt: T = Bezugszeit,where: T = reference time,

p = Bezugsversetzung,p = reference offset,

, = Geschwindigkeit des bezeichneten Punkts, = Speed of the designated point

309847/0941309847/0941

- 4b- - 4b-

Bei 78 wird das Zeitabtastmaß (At) berechnet. Der Scheitelpunkt jedes zeitabhängig gefilterten Tfellenzugs (wavelet) v/ird bestimmt, und anschließend wird der Abstand vom Scheitelpunkt des Wellenzugs zu einem Punkt des Uellenzugs bestimmt, der einem vorgewälilten Prozentsatz des Scheitelwerts entspricht. In einer bevorzugten Ausführungsform kann dieser vorbestimmte Prozentsatz des Scheitelwerts im Bereich von 0,8 bis 0,55 liegen. Der Abstand zwischen dem Scheitelpunkt und dem vorbestimmten Punkt entspricht dem zeitlichen Abtastmaß» Wenn der zeitabhängig gefilterte Wellenzug in Schritten zu jeweils einem Zehntel von 25 /um erzeugt v/ird, dann wird At zum nächstliegenden 1,25 yum-Schritt bestimmt. In der bevorzugten Ausführungsform ist dem Abtastmaß At eine Beschränkung in der Hinsicht auferlegt, daß der Bereich zwischen IOC und 250 um liegt, da dieser Zeitbereich die am besten durchführbaren Bandbreitenbedingungen überdeckt. At 78 the time sample amount (At) is calculated. The vertex of each time-dependently filtered wave train (wavelet) is determined, and then the distance from the vertex of the wave train to a point on the wave train is determined which corresponds to a preselected percentage of the peak value. In a preferred embodiment, this predetermined percentage of the peak can range from 0.8 to 0.55. The distance between the vertex and the predetermined point corresponds to the temporal sampling measure. If the time-dependent filtered wave train is generated in steps of one tenth of 25 / µm, then At is determined to the nearest 1.25 yum step. In the preferred embodiment the Abtastmaß At is imposed a limitation in the sense that the region between the IOC and 250 is to, as this time period covers the feasible best bandwidth conditions.

Das Mcveout-Abtastmaß Am wird bei 60 dadurcli bestimmt, daß zunächst das Feld der Versetzungsabstände aus der oben erwähnten Information im Kanal 25 entnommen wird. Aus dem Feld der Versetzungsabstände wird der Ort der Moveout-Bereichs-lmpulsantwort M^ zum nächstliegenden Zehntel von 25 lom entsprechend der folgenden Gleichung berechnet;The Mcveout sampling dimension Am is determined at 60 dadurcli, that first the field of the dislocation distances from the above-mentioned information in channel 25 is taken. The location of the Moveout area impulse response M ^ to nearest Tenths of 25 lom according to the following equation calculated;

309847/0941309847/0941

wobei gilt: i = 1, 2, ..., max (24, NFOLD),where: i = 1, 2, ..., max (24, NFOLD),

Tq = Bezugszeit (Mittelzeit des Zeitfensters), I!· = Absoluter Moveout-Wert (MittelmoveoutwertTq = reference time (mean time of the time window), I! = Absolute moveout value (mean moveout value

des Zeitfenster),of the time window),

X₁ = Abstand zum Versetzungswert i, Am· ss Fehler des absoluten Moveout-Werts (2,5 At), = Bezugsversetzung.X ₁ = distance to offset value i, Am · ss error of the absolute moveout value (2.5 At), = reference offset.

Die Berechnung mit Hilfe der Gleichung (8) ergibt eine Reihe von Nadeliapulsen, die in Fig.7 schematisch dargestellt sind. Die Faltung von M = ^M. mit dem zuvor bei 34 (Fig.3) bestimmten zeitabhängig gefilterten V/ellenzug, der bei Fig.8a schematisch dargestellt ist, ergibt eine Verschiebung des Scheitelwerts des Uellenzugs von T nach T = At¹ , wie in Fig.8b dargestellt ist. Das festgestellte Abtastmaß At wird zur Bestimmung des Moveout-Abtastmaßes mit Hilfe der folgenden Gleichung verwendet:The calculation with the aid of equation (8) results in a series of needle pulses, which are shown schematically in FIG. The convolution of M = ^ M. with the time-dependent filtered V / ellenzug previously determined at 34 (FIG. 3), which is shown schematically in FIG. 8a, results in a shift of the peak value of the Uellenzug from T to T = At ¹ , as shown in FIG. 8b. The determined sampling dimension At is used to determine the moveout sampling dimension with the aid of the following equation:

AtAm' _ 2,5 (At)² TT« - At¹ AtAm ' _ 2.5 (At) ² TT «- At ¹

Diese Bestimmung des Moveout-Abtastmaßes ergibt eine optimale Anpassung des Abtastmaßes an Δ t und an das zeitabhängige Filter. Venn weniger als die Datengesamtheit verarbeitet wird, ist die oben erwähnte Näherung des Moveout-Abtastmaßes zu fein, und somit wird dasMoveout-Abtastmaß in diesem Fall verdoppelt.This determination of the moveout sampling rate results in an optimal one Adaptation of the sampling dimension to Δ t and to the time-dependent filter. Venn processed less than the total data becomes, the above-mentioned approximation of the moveout sample rate is too fine, and thus the moveout sample rate becomes in this case doubled.

Das Dip-Abtastmaß Ad wird bei 82 berechnet. Sin SMSH ist als räumliche Aufsummierung oder als Moveout-Abtastungen von Daten für Punkte gleicher Tiefe definiert. Für einen N über 1 SMASH ist der Dip-Antv/ortoperator eine Folge von N gleich weit voneinander entfernten Nadelimpulsen , die einemThe dip sample amount Ad is calculated at 82. Sin SMSH is as spatial summation or as moveout samples defined by data for points of equal depth. For an N over 1 SMASH, the dip-answer operator is a sequence of N needle pulses equally spaced from each other,

309847/0941309847/0941

zusammengefassten (smashed ) Punkt gleicher Tiefe entsprechen. Diese Folge von Fadelimpulsen wird mit dem zuvor bestimmten, zeitabhängig gefilterten Wellenzug von Fig.9 gefaltet, der eine Einheitsamplitude aufweist. Wenn der korrekte Dip-Wert.geliefert wird, ist der Operator ein Nadelimpuls mit Einheitsamplitude, und der Abstand zwischen den N-Nadelimpulsen ist Null. Für einen Fehler von 2,5 wm pro Dip-Wert liegen die N Nadelimpulse jeweils 2,5yum voneinander entfernt, und sie haben eine Amplitude von
1/N. Die Dip-Antwortoperatoren werden auf diese Weise
mit dem in Fig.9 dargestellten, zeitabhängig gefilterten Wellenzug gefaltet, wobei jeder nachfolgende Djp-Operator einen um 2,5 Aim größeren Abstand zwischen den Nadelimpulsen als der vorhergehende Operator aufweist. Jede Faltungsserie ergibt eine Abnahme ^edes anschließenden Scheitelamplitudenwerts des Wellens.ugs· Fig.10 zeigt
die Herabsetzung der Amplitude des Wellenzugs auf 0,55 seines Einheitswerts. An diese Stelle ist das Dip-Abtastmaß Δ d der Abstand zwischen den Nadelimpulsen
des Dip-Antwortoperators, wie in Fig.11 dargestellt ist. Die Berechnungen des Dip-Abtastmaßes werden für jedes
aufeinanderfolgende Zeitfenster durchgeführt.combined (smashed) point correspond to the same depth. This sequence of thread pulses is convoluted with the previously determined, time-dependent filtered wave train of FIG. 9, which has a unit amplitude. If the correct Dip value is supplied, the operator is a needle pulse with unit amplitude and the spacing between the N needle pulses is zero. For an error of 2.5 µm per dip value, the N needle pulses are each 2.5 µm apart and have an amplitude of
1 / N. The dip answer operators are done this way
with the time-dependent filtered wave train shown in FIG. 9, each subsequent Djp operator having a distance between the needle pulses which is 2.5 Aim greater than that of the preceding operator. Each series of convolution results in a decrease in the subsequent peak amplitude value of the wave
reducing the amplitude of the wave train to 0.55 of its unit value. At this point, the dip scanning dimension Δ d is the distance between the needle pulses
of the dip response operator as shown in Figure 11. The dip sweep calculations are made for each
successive time windows carried out.

Bei 84 wird die Dip-Überlappung T_QD entsprechend der
folgenden Gleichung berechnet:At 84, the dip overlap T _QD becomes corresponding to FIG
calculated using the following equation:

wobei gilt: SMASH = Zahl der wahrend jedes Zeitfensters zuwhere: SMASH = number of times during each time window

betrachtenden Tiefenpunkte (vom Benutzer eingegeben),viewing depth points (entered by the user),

^Dmin ^{= minima}l^er Dip-Wert (vom Benutzer eingegeben) , ^D min ^{= minima} l ^he dip value (entered by the user),

^Dmax ^{= maxima}-l^er Dip-Wert (vom Benutzer eingegeben) . ^D max ^{= maxima} -l ^he dip value (entered by the user).

309847/0941309847/0941

In Fig.5 wird bei 86 eine Tabelle mit den Moveout-Abtastgrenzwerten berechnet und gespeichert. Fig.12 und 13 zeigen die Transformation der vom Benutzer definierten Geschwindigkeitsgrenzen in der Zeit-Moveout-Darsteilung. Fig.12 zeigt eine Geschwindigkeitsgrenze ähnlich der zuvor im Zusammenhang mit Fig.6 beschriebenen Grenze.An jedem in Fig.12 dargestellten Zeitpunkt existiert ein entsprechender Moveout-Wert, der mit Hilfe der Gleichung 7 berechnet werden kann. Diese Moveout-Werte für ein einziges Zeitfenster sind in Fig.13 aufgetragen. In Fig.13 sind auch Moveout-Spuren 90 dargestellt, die umÄm voneinander getrennt sind.Zur Erzielung der exakten Moveout-Werte der Schnittpunkte der Moveout-Spuren mit den in Fig.13 gezeigten Moveout-Grenzen muß eine lineare Interpolation durchgeführt v/erden. Diese Interpolationen werden in der nachfolgend beschriebenen Weise berechnet und sie werden bei 88 (Fig.5) in die modifizierte Moveout-Abtasttabelle eingefügt.In FIG. 5, a table with the moveout sampling limit values is shown at 86 calculated and saved. 12 and 13 show the transformation of the speed limits defined by the user in the time moveout representation. Fig.12 shows a Speed limit similar to the limit previously described in connection with Fig. 6, at each shown in Fig. 12 At the time, there is a corresponding moveout value that can be calculated using equation 7. These Moveout values for a single time window are shown in Fig. 13 applied. 13 also shows moveout tracks 90, which are separated from each other by Äm. To achieve the exact Moveout values of the intersection points of the moveout tracks with the moveout limits shown in FIG. 13 must be linear interpolation carried out v / earth. These interpolations are calculated in the manner described below, and they are performed at 88 (Fig. 5) inserted into the modified moveout scan table.

Fig.14 zeigt ein Funktionsflußdiagramia für die Berechnung der Moveout-Abtastgrenztabelle bei 86. Beim Schritt 100 wird bestimmt, welche Knickpunkte der vom Benutzer eingegebenen Geschwindigkeitsübersicht bestimmt werden sollen. Für die beim Schritt 100 dargestellte Gleichung hat der äußere Schleifenzählindex i die Vierte 2, 3, ... N. Das Abtastmaß At wird für > 0,75 X_R auf 100 Millisekunden und für < 0,75 X_R auf 10 Millisekunden eingestellt. Beim Schritt 102 werden die Geschwindigkeitsgrenzindices festgestellt, die die gewünschte Zeit begrenzen. Beim Schritt 103 wird entschieden, ob die vom Benutzer eingegebenen Geschwindigkeitsgrenzen größer oder kleiner als Null sind, damit die Eingangsgeschwindigkeitsmöglichkeiten geprüft werden. Wenn die vom Benutzer eingegebenen Geschwindigkeitsgrenzen kleiner als Null sind, wird beim Schritt 104 entschieden, ob Yy- kleiner als Null ist oder nicht. Wenn V_K 14 shows a functional flow diagram for the computation of the moveout scan limit table at 86. At step 100 it is determined which breakpoints of the speed overview entered by the user are to be determined. For the equation shown in step 100, the outer loop count index i has the fourth 2, 3, ... N. The sampling dimension At is set _{to 100 milliseconds for> 0.75 X R} and to 10 milliseconds for <0.75 X _R. At step 102, the speed limit indices that limit the desired time are determined. In step 103 a decision is made as to whether the speed limits entered by the user are greater or less than zero, so that the input speed possibilities are checked. If the speed limits entered by the user are less than zero, a decision is made in step 104 as to whether or not Yy- is less than zero. If V _K

309847/094 1309847/094 1

23239B323239B3

nicht kleiner als Null ist, dann wird "beim Schritt 106 nach der Gleichung 7 der Moveout-Wert M^ für (V^, T^) berechnet. Beim Schritt 108 v/ird der Zählindex um 1 erhöht. Es wird beim Schritt 108 festgestellt, daß v/egen der unendlichen Geschwindigkeit der Moveout-Wert gegen Null geht, so daß die Schleife beim Schritt 110 beendet werden kann, da die nachfolgenden Zeiten konstant sind. Beim Schritt 112 wird festgestellt, daß nur zwei Eintragungen in der Tabelle vorhanden sind und daß beide Eintragungen die gleiche Zeit aufweisen. Somit ist die Startzeit korrekt, und keine zusätzlichelnterpolation ist erforderlich, so daß die Schleife beim Schritt 110 beendet wird.is not less than zero, then the moveout value M ^ for (V ^, T ^) is calculated in step 106 according to equation 7. In step 108 the counting index is increased by 1. In step 108, it is determined that v / ee the infinite speed the moveout value tends to zero, so that the loop can be ended since the subsequent times are constant at step 110. At step 112 it is determined that there are only two entries in the table and that both entries have the same time, so the start time is correct and no additional interpolation is required, so the loop at step 110 is exited.

Beim Schritt 114 wird die Geschwindigkeit V. bei der Zeit T^ interpoliert. Die Interpolation beim Schritt 114 wird ausgeführt, da. es erwünscht ist, den Endpunkt der Moveout-Spur festzustellen, bei dem die Spur zwischen zwei der festgestellten Moveout-Punkte fällt. Beim Schritt 114 werden Daten erzeugt, damit in einer ansdiießend zu beschreibenden ¥eise eine Interpolation ermöglicht wird. Beim Schritt 116 wird der Moveout-Wert l-L für bestimmte Geschwindigkeitsund Zeitwerte (V^, T.) berechnet. Beim Schritt 118 wird entschieden, ob noch weitere Geschwindigkeitswerte zu berechnen sind. Wenn weitere Geschwindigkeitswerte zu berechnen sind, werden der Schritt 100 und die nachfolgenden Schritte erneut ausgeführt. Bei dem erneuten Durchlauf wird zum Zeitwert vor dem Durchlauf das Abtastmaß At addiert. "Wenn beim Schritt 118 entschieden wird, daß keine weiteren Geschwindigkeitswerte zu berechnen sind, werden beim Schritt 120 der Wert M_n auf den minimalenIn step 114, the speed V is interpolated at the time T ^. The interpolation at step 114 is carried out because. it is desirable to determine the end point of the moveout trace at which the trace falls between two of the detected moveout points. At step 114, data is generated to enable interpolation in a manner to be described below. In step 116 the moveout value IL is calculated for certain speed and time values (V ^, T.). In step 118 it is decided whether further speed values are to be calculated. If further speed values are to be calculated, step 100 and the subsequent steps are carried out again. In the next run, the sampling dimension At is added to the time value before the run. "If it is decided in step 118 that no further speed values are to be calculated, the value M _{n is set} to the minimum in step 120

30 98 4 7 /0930 98 4 7/09

Moveout-Wert und der Wert t auf den Wert t_max gesetzt, worauf die Schleife beim Schritt 110 beendet wird. Die mit Hilfe des in Fig.14 dargestellten Vorgangs bestimmten Moveout-Abtastwerte und Zeitwerte v/erden bei 86 (Fig.5) in die Moveout-Abtastgrenzwerttabelle eingegeben.Moveout value and the value t are set to the value t _max , whereupon the loop is ended in step 110. The moveout sampling values and time values determined with the aid of the process shown in FIG. 14 are entered into the moveout sampling limit value table at 86 (FIG. 5).

Die Gesamtzahl der benötigten Moveout-Abtastungen N_ wird aus der folgenden Gleichung berechnet:The total number of moveout samples required N_ is calculated from the following equation:

N_s = 2S + 1 (11)N _s = 2S + 1 (11)

wobei gilt:where:

1515th

H = (12b)H = (12b)

Aus den obigen Gleichungen kann der Wert der minimalen und maximalen Moveout-Abtastungen folgendermaßen bestimmt werden:From the equations above, the value of the minimum and maximum moveout samples can be determined as follows:

Minimaler Abtastwert = M - ^SAm<T''^I _mi_n (13a) Maximaler Abtastwert = Vi + S 4 m >M_max (13b)Minimum sample value = M - ^{SAm <T} '' ^I _m i _n (13a) Maximum sample value = Vi + S 4 m> M _max (13b)

Die Berechnung der obigen Werte ergibt eine Tabelle mit minimalen und maximalen Abtastwerten zusammen mit der Gesamtzahl der erforderlichen Abtastungen.Calculating the above values results in a table of minimum and maximum samples along with the total the required scans.

Es ist erwünscht, eine Dip-Überlappung in die bei 86 bestimmte Tabelle einzubauen, so daß keine Daten infolge der anschließenden Anwendung von Filtern verlorengehen. Die Moveout-Abtasttabelle wird daher bei 86 so modifiziert, daß sie eine Dip-Überlappung enthält. Überdies ist eine Interpolation erforderlich, damit der exakte Moveout von Spuren erzielt wird.It is desirable to have a dip overlap in the one at 86 to include a specific table so that no data is lost as a result of the subsequent application of filters. The moveout scan table is therefore modified at 86 to include dip overlap. Moreover, there is one Interpolation required in order to achieve the exact moveout of tracks.

309847/094 1309847/094 1

Der Interpolations-Vorgang der Moveout-Abtasttabelle beginnt mit der Bestimmung derjenigen Moveout-Sintragungen, die den Abtast-Moveout begrenzen. Der Abtast-Moveoutbereich für T_ß ist M + (S-1) Am-^I - (S-2) Am . Der Abtast-Moveoutbereich für T^ ist Έ + (S-2)Äm->M -(S-1)Z\m, Es sei ein Abtast-Moveout von Mg angenommen, bei dem giltj M + (S-DAm < M_s < H - (S-2) 4m. Die Moveouteintragungen werden dann in der zeitkonstant abgetasteten Tabellengrenze Mn derart abgetastet, daß gilt:The moveout scan table interpolation process begins by determining those moveout entries that limit the scan moveout. The scan moveout range for T _ß is M + (S-1) Am- ^ I - (S-2) Am. The scan moveout range for T ^ is Έ + (S-2) Äm-> M - (S-1) Z \ m, Assume a scan moveout of Mg, in which j M + (S-DAm <M _s <H - (S-2) 4m. The moveout entries are then scanned in the table boundary Mn, which is scanned at a constant time, in such a way that the following applies:

Die Startzeit für die Abtastung mit dem V/ert Mg ist folgendermaßen definiert:The start time for sampling with the V / ert Mg is defined as follows:

(1 - 1) At ₊ () At (14)(1 - 1) At ₊ () At (14)

309847/094 1309847/094 1

Die Moveout-Äbtasttabelle ist in der folgenden Tabelle
angegeben:The moveout scan table is in the following table
specified:

Tabelle ITable I.

AbtastindexSample index AbtastwertSample S AmS Am AmAt the StartzeitStart time EndzeitEnd time 11 M +M + (s-D(s-D AmAt the φ ( O\ φ
B^ ' OD φ ( O \ φ
B ^ 'OD T_E(3) ₊ T _E (3) ₊ 22 H +H + (S-2)(S-2) AmAt the ^XB^K^^} ~ ¹OD ^X B ^K ^ ^} ~ ¹ OD T_E(3) ₊ T _E (3) ₊ 33 M +M + (S-3)(S-3) ¹B^ ^; " OD ¹ B ^ ^; "OD I_E(4) ♦I _E (4) ♦ 44th M +M + ΓΠ / "X \ ΓΠ
¹B^^' ^X0DΓΠ / "X \ ΓΠ
¹ B ^^ ^'X 0D 1.(5) ♦1. (5) ♦

^Ns- ^N s- 33 H -H - 33 H -H - 11 M -M - N^N ^ M -M -

(S-3) Am T_B(N_s-4)-T_0D (S-2) Am T_B(N_s-3)-T_0D (S-DAm T_B(N_s-2)-T_0D (S-3) Am T _B (N _s -4) -T _0D (S-2) Am T _B (N _s -3) -T _0D (S-DAm T _B (N _s -2) -T _0D

T_E(N_S-1)₊T_OD T _E (N _S -1) ₊ T _OD

t_e(n_s-D₊t_0D t _e (n _s -D ₊ t _0D

T^Ne-D+TT ^ Ne-D + T

Entsprechende Punkte des Diagramms von Fig.13 sind in
der Tabelle I zu erkennen. Die in der Tabelle enthaltenen Werte enthalten die Dip-Überlappung um T_QD, damit verhindert wird, daß während der Anwendung von Filtern Daten verlorengehen, wie oben erwähnt wurde.Corresponding points in the diagram of FIG. 13 are in
Table I can be seen. The values in the table include the dip overlap around T _QD in order to prevent data from being lost during the application of filters, as mentioned above.

Für die Tabelle I ist die vom Benutzer eingegebene
Geschwindigkeit-Grenzspezifikation folgendermaßen :For table I is that entered by the user
Speed limit specification as follows:

v_L(D, V₁₁(D, T(D > οv _L (D, V ₁₁ (D, T (D> ο

V_L(2), V_H(2), T(2) V_L(n),V _L (2), V _H (2), T (2) V _L (n),

rj(rj (

T(n) < AufzeichnungslängeT (n) <record length

309847/0941309847/0941

Für die Tabelle I gelten folgende Grenzbedingungen: Grenzbedingungen The following boundary conditions apply to Table I: Boundary conditions

1. V_L(i+i) > V_L(i) und V_H(i+i) >. V_H(i)1. V _L (i + i)> V _L (i) and V _H (i + i)>. V _H (i)

2. Nur Vjr(i) kann negativ sein.2. Only Vjr (i) can be negative.

3. Wenn Vrr(i) negativ ist, dann müssen alle folgenden Werte V_H(i +1), V_H(i+2) ...V_H(n) ebenfalls kleiner als Null sein.3. If Vrr (i) is negative, then all following values V _H (i +1), V _H (i + 2) ... V _H (n) must also be less than zero.

Wenn es erforderlich ist, wird beim Schritt 122 (Fig.5) eine Modifikation ohne Beteiligung der Datengesamtheit durchgeführt. Dies -wird zur Einsparung von Verarbeitungszeit ausgeführt, daIf necessary, at step 122 (Fig. 5) a Modification carried out without involvement of the data population. This is done to save processing time because

sich gezeigt hat, daß es nicht notwendig ist,Teile der Datengesamtheit ebeno- stark wie die Gesamtheit der Daten abzutasten. Fig.15 zeigt die Anwendung der Teilverarbeitungstechnik.In dem Diagramm von Fig.15 sind die Moveout-Werte über der Zeit aufgetragen, wobei gilt : T. =0,3 ^ρ/^πο* T₂ = 0,5 X_R/^VCO' T₃ = 0,75 X_R/V_C0. Die Steigung der T.-Linie beträgt At/4Am, die Steigung der Tp-Linie beträgt At/2 Am, und die Steigung der T^-Linie beträgt / Am.It has been shown that it is not necessary to scan parts of the data as a whole as the whole of the data. Fig. 15 shows the application of the partial processing technique. In the diagram of Fig. 15 the moveout values are plotted against time, where: T. = 0.3 ^ ρ / ^ πο * T ₂ = 0.5 X _R / ^V CO 'T ₃ = 0.75 X _R / V _C0 . The slope of the T. line is At / 4Am, the slope of the Tp line is At / 2 Am, and the slope of the T ^ line is / Am.

Es ist zu erkennen, daß für im wesentlichen volle Faltdaten die Daten mit dem Moveout-Abtastschritt Am abgetastet werden. Für die Moveout-Werte während des Zeitbereichs zwischen Tp und T, ist der Moveout-Abtastschritt verdoppelt. Ebenso ist der Moveout-Abtastschritt für jeden darauffolgenden Zeitbereich von T. bis T₂ und von T₁ bis 0 verdoppelt.Die obige TeildatenmodifikationIt can be seen that for substantially full folded data, the data are scanned with the moveout scanning step Am. For the moveout values during the time range between Tp and T, the moveout sampling step is doubled. Likewise, the moveout sampling step is doubled for each subsequent time range from T. to T ₂ and from T ₁ to 0. The above partial data modification

309847/0941309847/0941

wird nur dann durchgeführt, wenn der vom Benutzer eingegebene Wert T₁ kleiner oder gleich T^₁ oder 0,75 der Gesamtheitszeit ist, was anzeigt, daß ein zu großes Abtastmaß verwendet wird.is only performed if the user entered value T _{1 is} less than or equal to T ^ ₁ or 0.75 of the total time, indicating that too large a sample size is being used.

Beim Schritt 130 von Fig.5 werden die Zeitgrenzen eines Punkts gemeinsamer Tiefe aus der folgenden Gleichung berechnet:At step 130 of Figure 5, the time limits of a point common depth is calculated from the following equation:

minmin

maxMax

(15a)(15a)

= T(i+1) + T_nn + (M + SAm) < Aufzeichnungs- (15b)= T (i + 1) + T _nn + (M + SAm) <record- (15b)

längelength

wobei gilt: M + SAm = M„,_QV(aus der Tabelle I),where: M + SAm = M ", _QV (from Table I),

IQcLjCIQcLjC

^T0D ⁼ Zeitbreite der Dip-Überlappung. ^T 0D ^{= time} width of the dip overlap.

Beim Schritt 133 werden zweidimensionale Filter für die spätere Verwendung bei dem Verfahren berechnet und gespeichert. Bei der Auslegung der zweidimensionalen Filter wird zunächst eine Matrix aus Koordinatenpositionen gebildet.At step 133, two-dimensional filters are created for later Use in the process calculated and stored. When designing the two-dimensional filter, a Matrix formed from coordinate positions.

Es sei angenommen, daß X der gegebene Eingangsdatenvektor ist,Assume that X is the given input data vector,

^X1 ^X 1

wobei gilt: Xwhere: X

X₂ X ₂

^X3 ^X4 ^X 3 ^X 4

XrrXrr

(16)(16)

X» —X »-

₁, M₁) i = 1, 2, ... 5 ₁ , M ₁ ) i = 1, 2, ... 5

Die Werte t. und M₁ sind die Koordinaten von χ in der Zeit-Mouveout-Ebene. The values t. and M ₁ are the coordinates of χ in the time mouveout plane.

309847/094 1309847/094 1

Fig. 16 dient der Veranschaulichung des Aufbaus des zweidimensionalen Filters sowie der Verdeutlichung der anschliessenden Anwendung dieses Filters bei den späteren Verarbeitungsschritten. Fig.i6 zeigt die Koordinatenzeit und Moveout und es sind fünf Eingangsdatenpunkte 1, 2, 3, 4 und 5 angegeben, die mit stark ausgezogenen Punkten dargestellt16 is used to illustrate the structure of the two-dimensional Filters as well as the clarification of the subsequent application of this filter in the later Processing steps. Fig.i6 shows the coordinate time and moveout and there are five input data points 1, 2, 3, 4 and 5 indicated, those shown with strongly drawn out dots

Aus den gewünschten Interpolationspunkten wird dann eine Matrix gebildet.A matrix is then formed from the desired interpolation points.

Y sei der gewünschte Ausgangsvektor.Let Y be the desired output vector.

^y7 ^y 7

^r8 ^r 8

(17)(17)

y-17y-17

Die Werte t. und M. sind die Koordinaten von y in der Zeit-Mo ve out-Ebene .The values t. and M. are the coordinates of y in time Mo ve out level.

Nun muß eine Filtermatrix A (eine Gruppe, von Filtern) gefunden werden, so daß gilt:Now a filter matrix A (a group of filters) has to be found so that the following applies:

-Y = A¹X-Y = A ¹ X

(18)(18)

309847/0941309847/0941

wobei gilt A = where A =

^a1,6 ^a1,7 ^a1,8 ^a1,18 ^a 1.6 ^a 1.7 ^a 1.8 ^a 1.18

^a2,6 ^a2,7 ^a2,8 ^a2,18 ^a 2.6 ^a 2.7 ^a 2.8 ^a 2.18

^a3,6 ^a3,7 ^a3,8 * ^a3,18 ^a 3.6 ^a 3.7 ^a 3.8 * ^a 3.18

^a4,6 ^a4,7 ^a4,8 · ^a4,18 ^a 4.6 ^a 4.7 ^a 4.8 ^a 4.18

^a5,6 ^a5,7 ^a5,8 ^a5,18 ^a 5.6 ^a 5.7 ^a 5.8 ^a 5.18

Jede Spalte der Matrix 4 ist ein Fünf-Punkt-Interpolationsfilter. Each column of the matrix 4 is a five point interpolation filter.

Unter Anwendung des Kriteriums des minimalen mittleren Fehlerquadrats j[ boxcar spectrum model) und unter Berücksichtigung der Zeit-Moveout-Schräglage, kann die Filtermatrix A folgendermaßen hergeleitet werden: Using the minimum mean criterion Error square j [boxcar spectrum model) and taking into account the time moveout inclination, the filter matrix A can be derived as follows:

A = ΙΓ¹ΦA = ΙΓ ¹ Φ

wobei
gilt: R =whereby
the following applies: R =

^r1,1 ^r1,2 ^Γ1,3 ^Γ1,4 ^Γ1,5 ^r 1.1 ^r 1.2 ^Γ 1.3 ^Γ 1.4 ^Γ 1.5

^r2,1 ^Γ2,2 ^Γ2,3 ^Γ2,4 ^Γ2,5 ^r 2.1 ^Γ 2.2 ^Γ 2.3 ^Γ 2.4 ^Γ 2.5

^Γ^ 1 ^Γ^ ? ^Γ^ ^ ^Γ"3 4 ^*^ 5 ^Γ ^ 1 ^Γ ^? ^Γ ^ ^ ^Γ "3 4 ^ * ^ 5

^Γ4,1 ^4,2 ^r4,3 ^r4,4 ^r4,5 ^Γ 4.1 ^ 4.2 ^r 4.3 ^r 4.4 ^r 4.5

^r5,V ^r5,2 ^r5,3 ^r5,4 ^r5,5 ^r 5, V ^r 5.2 ^r 5.3 ^r 5.4 ^r 5.5

Die zur Entwicklung der Matrix g verwendete Gleichung lautet folgendermaßen:The equation used to develop the matrix g is as follows:

SimSim

0.4π0.4π

+ tane+ tane

Sin<Sin <

^ -t r ^ -t r

0.4ΤΓ0.4ΤΓ

- tane (t^- tane (t ^

- taue- thaw

t_i j Jt _i j J

309847/09309847/09

wobei gilt: ϊwhere: ϊ

At, a mAt, a m

— ι» d_f ..., 5»- ι » d _f ..., 5»

— Λ Ο ¹^- Λ Ο ¹ ^

= constant= constant

= Abtastintervall (in ms) längs der Zeit- und der= Sampling interval (in ms) along the time and the

Moveout-Achse und
θ = Schrägwinkel (gegeben tanG).Moveout axis and
θ = skew angle (given tanG).

Weiterhin gilt:The following also applies:

*1,6 .*1,7 ^Φ1,18* 1.6. * 1.7 ^Φ 1.18

^Φ2,6 ^Φ2,7 ····.· ·₂,ΐ8 ^Φ 2.6 ^Φ 2.7 ····. · · ₂ , ΐ8

^Φ3,6 ^Φ3,7 ^Φ3,18 ^Φ 3.6 ^Φ 3.7 ^Φ 3.18

^Φ4,6 ^Φ4,7 ^Φ4,18 ^Φ 4.6 ^Φ 4.7 ^Φ 4.18

^Φ5,6 ^Φ5,7 ^Φ5,18 ^{Φ Φ} 5.6 5.7 ^Φ 5.18

(21)(21)

Φ kann aus folgender Gleichung entwickelt werden: Φ can be developed from the following equation:

SinSin

0.4π0.4π

^t₁) + tane^ t ₁ ) + tane

-t_ij +-t _i j +

0.4U0.4U

(l+tan6)At(l + tan6) At

C22)C22)

mit i = 1, 2, ..., 5
Ö = 6, 7, ..., 18with i = 1, 2, ..., 5
Ö = 6, 7, ..., 18

undand

309847/0941309847/0941

Als Beispiel für die Bildung von Filtern nach der Erfindung sei angenommen, daß die relative Lage von X und Y in der Zeit-Moveout-Ebene wie in Fig.16 gegeben ist, und daß gilt:As an example of the formation of filters according to the invention, it is assumed that the relative position of X and Y in the time moveout plane as given in Fig. 16 is, and that the following applies:

tan© = 0,4
At = 4(ms)
Am = 10(ms)tan © = 0.4
At = 4 (ms)
Am = 10 (ms)

Zur Bestimmung von r^ ■* und Φ^^ To determine r ^ ■ * and Φ ^^

aus Fig.16 ergeben sich die folgenden Zeit-Moveout-Koordinaten vonthe following time moveout coordinates result from FIG

V ^X3> ^yi6'V ^X 3> ^y i6 '

ZeitTime MoveoutMoveout ^Xl ^X l AtAt -Am-At the X₃ X ₃ 00 00 ruru AtAt Am
"TAt the
"T Y₁₆ Y ₁₆ At
"T At
"T AmAt the

Die relativen Positionen vonThe relative positions of

und (y^gand (y ^ g

sind;are;

^t3~^tl ^t 3 ~ ^t l

^m3"^ml ^{= m} 3 " ^m l ⁼

= m = 10(ms)= m = 10 (ms)

. _ -At At _ -3At ^t16~^t14 ~2 T 3. _ -At At _ -3At ^t 16 ~ ^t 14 ~ 2 T 3

Am Am Am _ ^m16"^m14 ⁼ "T 4" " "T "Am Am Am _ ^m 16 " ^m 14 ⁼ " T 4 """T"

309847/0941309847/0941

Aus Gleichung 20 ergibt sich somit :From equation 20 we get:

Sin / [-4 + 0.4 · 10]Sin / [-4 + 0.4 · 10]

^rl,3 L-4 +0.4 · IUJ ^r l, 3 L-4 +0.4 IUJ

Sin / [10-0.4 · 10]Sin / [10-0.4 · 10]

„ 1(1+0.4)·10"1 (1 + 0.4) · 10

^x UU-0.4 · -4J ^x UU-0.4 · -4J

ir · SJn(0.4tr) ir SJn (0.4tr)

0.116 * ΓΤ750.116 * ΓΤ75

Durch Einsetzen in die Gleichung 22 erhält man:Substituting it into equation 22, we get:

f \f \

Sin / [-3+0.4 · 2.5]^V Sin / [-3 + 0.4 * 2.5] ^V

^Φ14,16 ⁼ t-3 + 0.4 · Z.bJ ^Φ 14.16 ⁼ t-3 + 0.4 Z.bJ

Sin I [2.5 - 0.4 · -31]Sin I [2.5 - 0.4 · -31]

1(1+0.4) 10 1 (1 + 0.4) 10

Sin (0.172t) · Sin (0.128ττ) Z 3T7 Sin (0.172t) Sin (0.128ττ) Z 3T7

309847/0941309847/0941

Die beim Schritt 132 (Fig.5) berechneten zweidimensionalen Filter werden für die spätere Vervrendung zum Interpolieren von Amplituden während der dreidimensionalen Ereignissichtung gespeichert.The two-dimensional ones calculated in step 132 (FIG. 5) Filters are for later use to interpolate amplitudes during three-dimensional event viewing saved.

Nach Fig.5 wird die Zahl der Dip-Ebenen beim Schritt 134 aus der folgenden Gleichung berechnet:According to Figure 5, the number of dip levels in step 134 calculated from the following equation:

D^f - D' j
ND₁P ₌ -JSS2 2iS (23)D ^f - D ' j
ND ₁ P ₌ -JSS2 2iS (23)

AdAd

wobei gilt:where:

^Dimin^{= D}min ~ 2~
max⁼ max ⁺ 2Γ" ^Di min ^{= D} min ~ 2 ~
max ⁼ max ⁺ 2Γ "

Der Benutzer hat die minimalen und maximalen Dip-Parameter eingegeben. Der minimale Dip-Parameter wird um die Hälfte von Ad reduziert, damit gewährleistet wird, daß der gewünschte Benutzerbereich wegen Datenverlusten an den Bereichsrändern infolge von Interpolationsschritten völlig erfaßt wird. Mit der Beendigung des Schritts 134 ist auch die Parameterinitialisierung beendet.The user has entered the minimum and maximum dip parameters. The minimum dip parameter is cut in half of Ad is reduced, so that it is guaranteed that the desired user area due to data loss to the Area edges is completely detected as a result of interpolation steps. With the completion of step 134 is also the parameter initialization is finished.

Move out-AbtastaufbauMove out scan setup

Fig.17 zeigt ein Funktionsflußdiagramm für den Schritt 36 zum Aufbau der Moveout-Abtastung , der oben in Zusammenhang mit Fig.3 beßchrieben wurde. Die CDP-Sammelaufzeichnungen werden beim Schritt 150 eingegeben, und beim Schritt 152 wird eine Spur für einen gemeinsamen Tiefenpunkt gelesen. Die Spur wird beim Schritt 154 in Schritten von At abgetastet, wenn sie nicht bereits zuvor mit dem gewünschten Wert von At abgetastet worden ist. Das Wiederabtasten erfolgt mit herkömmlichen Interpolationsfiltcrn. Beim17 shows a functional flow diagram for step 36 to set up the moveout scan, related to the above was described with Fig.3. The CDP collective records are entered at step 150 and at step 152 a track for a common depth point is read. The track is scanned in steps of At at step 154, if it has not already been sampled beforehand with the desired value of At. The resampling takes place with conventional interpolation filters. At the

309847/0941309847/0941

Schritt 156 wird die Spur des gemeinsamen Tiefenpunkts um M₁ und um irgendeine gewünschte Störkorrekturgröße verschoben. Solche Störkorrekturen können durch irgendein herkömmliches digitales Störkorrekturverfahren bestimmt werden. Der Wert I>L ergibt sich aus der folgenden Gleichung:Step 156 shifts the track of the common depth point by M ₁ and any desired clutter correction amount. Such clutter corrections can be determined by any conventional digital clutter correction technique. The value I> L results from the following equation:

M. = 1/To² + Ρ- /Tn Am)² + 2(n Am)To J- To (24) ¹ γ ^Ä24M. = 1 / To ² + Ρ- / Tn Am) ² + 2 (n Am) To J- To (24) ¹ γ ^Ä 24

wobei gilt:where:

ηΔιη = Moveout-Wert bei J ηΔιη = moveout value at J

Am = Moveout-AbtastmaßAm = moveout sampling dimension

An = eine ganze Zahl , die systematisch um Eins erhöht wird, bis ^_3x überschritten ist.An = an integer that is systematically increased by one until ^ _{3x is} exceeded.

Die verschobene Spur wird dann beim Schritt 158 innerhalb des kumulativen Feldes summiert,das für die bestimmte Moveout-Abtastung aufgebaut ist. Beim Schritt 160 wird entschieden, ob weitere Moveout-Abtastspuren erzeugt werden sollen oder nicht. Wenn weitere Spuren erzeugt werden sollen, wird die nächste Moveout-VerSchiebung durch eine Erhöhung beim Schritt 162 berechnet, und die Spur für den Punkt gemeinsamer Tiefe wird beim Schritt 156 um diese Moveout-VerSchiebung verschoben. Die verschobene Spur wird dann beim Schritt 158 innerhalb des kumulativen Feldes für die bestimmte Moveout-Verschiebung summiert. Die Schleife mit den Schritten 156 bis 162 wird erneut durchlaufen, bis die Spur des Punkts gemeinsamer Tiefe, die bearbeitet wird, um jeder der Moveout-VerSchiebungen verschoben worden ist.The shifted track then becomes within at step 158 of the cumulative field established for the particular moveout sample. At step 160 decided whether to generate more moveout scan tracks should be or not. If more tracks are to be created, the next moveout shift is carried out an increment is calculated at step 162 and the track for the point of common depth is increased by that at step 156 Moveout shifted. The moved track becomes then summed at step 158 within the cumulative field for the particular moveout offset. the The loop with steps 156 to 162 is repeated until the trace of the point of common depth, the is processed to each of the moveout shifts has been postponed.

Wenn die Spur des Punkt gemeinsamer Tiefe um jede der Moveout-VerSchiebungen verschoben und innerhalb jedes einerWhen the trail of the point of common depth is shifted by each of the moveout shifts and within each one

309847/0941309847/0941

Moveout-Verschiebung entsprechenden Feldes summiert worden ist, wird "beim Schritt 164 entschieden, ob weitere CDP-Sammelspuren für die Verarbeitung übrigbleiben. Wenn dies der Fall ist, wird beim Schritt 152 eine weitere Spur eines Punkts gemeinsamer Tiefe gelesen , und die oben beschriebenen Schritte werden erneut durchlaufen. Wenn Sede der CDP-Sammelspuren bearbeitet worden ist, ist ein summiertes kumulatives ~ Feld für ^ede Moveout-Verschiebung aufgebaut worden.Moveout shift has been summed, it is determined at step 164 whether there are more CDP accumulation tracks left for processing. If so, another track of a common depth point is read in step 152 and the steps described above When each of the CDP collection tracks has been processed, a summed cumulative field has been built for each moveout.

Diese Spuren werden beim Schritt 166 entsprechend der Zahl der CDP-Spurbeiträge, die für die Information jedes kumulativen Feldes zur Verfügung stehen, maßstäblich verändert (recovery scaled). Eine solche maßstäbliche Veränderung mit einem Teilerfaktor ist erforderlich, da einige der Spuren wegen der Grenzflanke (cutoff ramp) des Benutzers nicht so stark wie andere Spuren zu der Summe beitragen. In einem Puffer wird ein Zählerstand der beitragenden Spuren, die in dem Feld gespeichert sind, festgehalten, und die Amplitude der Durchschnittsspuren wird entsprechend der Zahl der verwendeten Spuren gemäß dem folgenden Teilerfaktor reduziert:These traces are made in step 166 according to FIG Number of CDP track contributions that are required for the information of each cumulative field are available, changed to scale (recovery scaled). Such a scale Modification with a division factor is necessary because some of the traces are due to the cutoff ramp of the user do not contribute to the total as much as other tracks. A counter reading is stored in a buffer of the contributing traces stored in the field and the amplitude of the average traces is reduced according to the number of tracks used according to the following division factor:

scale _
^fact0r ~scale _
^fact0r ~

wobei gilt: i = Zahl der zu einem bestimmten Abschnittwhere: i = number of to a certain section

der Moveout-Spur beitragenden Spuren,tracks contributing to the moveout track,

NFOLD = maximale Anzahl der beitragenden Spuren. Die gemittelten Moveout-Abtastspuren werden beim Schritt 168 zum Plattenspeicher ausgegeben.NFOLD = maximum number of contributing tracks. The averaged moveout scan traces are at step 168 output to disk storage.

Zur weiteren Erklärung des in Fig.17 beschriebenen Moveout-Abtastaufbauverfahrens wird nun auf die Figuren 2 bisTo further explain the moveout scan setup method described in FIG will now refer to Figures 2 to

309847/0941309847/0941

Bezug genommen.Wie in Fig.2 dargestellt ist,werden die Spuren für Punkte gemeinsamer Tiefe längs vier hyperbolischer Moveout-Kurven A, B, C und D summiert. Jede der Moveout-Kurven A bis D sind durch das Moveout-Abtastfeld in Millisekunden pro Abtastung (dm) voneinander getrennt. Wie oben erwähnt wurde, wird jede der in Fig.2 dargestellten Spuren zunächst entsprechend dem Moveout-Wert, der von der Moveout-Kurve A beschrieben ist,zeitlich verschoben. Die verechobenen Spuren werden dann summiert. Anschließend v/erden die in Fig.2 dargestellten Spuren für jede der Moveout-Abtastspuren B bis D verschoben und summiert. Fig. 18 zeigt eine typische resultierende Moveout-Abtastspursammelaufzeichnung für einen Tiefenpunkt j.As shown in Figure 2, the Tracks for points of common depth along four hyperbolic moveout curves A, B, C, and D are summed. Each of the Moveout curves A through D are through the moveout scan field separated from each other in milliseconds per scan (dm). As mentioned above, each of the in Fig.2 The tracks shown initially in time according to the moveout value described by the moveout curve A. postponed. The displaced traces are then summed up. Subsequently, those shown in FIG. 2 are grounded Tracks for each of the moveout scan tracks B through D shifted and summed. Fig. 18 shows a typical resultant Moveout scan track bulk record for a depth point j.

Dip-AbtastaufbauDip scanning structure

Nachdem der Moveout-Abtastaufbau beendet worden ist, wird der in Fig.3 beim Schritt 33 angegebene Dip-Abtastaufbau ausgeführt. Fig.19 zeigt ein Funktionsflußdiagramm der Schritte des Dip-Abtastaufbaus. Kurz gesagt werden die Moveout-Abtastspuren für jeden Tiefenpunkt innerhalb eines bestimmten Raumfensters systematisch um verschiedene Dip-Werte verschoben, dann summiert und maßstäblich verändert, damit ein Zahlenfeld erzeugt wird, in dem jeder Amplitudenwert eine Funktion der Zeit, des Moveouts und des Dips ist. Die auf diese Weise erzeugte Zeit-Moveout-Fläche wird als Dip-Abtastung bezeichnet.After the moveout scan setup has been completed, the dip-scanning structure indicated in FIG. 3 in step 33 executed. Fig. 19 shows a functional flow diagram of the steps of the dip-scan construction. In short, the Moveout scan tracks for each depth point within a given space window systematically around different ones Dip values shifted, then summed and scaled to create a number field in which each amplitude value is a function of time, moveout and dip. The time moveout area created in this way is called dip scanning.

Wie aus Fig. 19 zu erkennen ist, werden die entsprechend dem inFig.17 dargestellten Vorgang berechneten Moveout-Abtastspuren beim Schritt 200 aus dem Speicher eingegeben. Beim Schritt 202 wird, eine einzelne Moveout-Abtastspur eingelesen, und der Dip-Wert für die Spur wird anfänglich beim SchrittAs can be seen from FIG. 19, the moveout scanning tracks calculated in accordance with the process shown in FIG. 17 are input from the memory in step 200. At step 202 , a single moveout scan track is read and the dip value for the track is initially set at step

309847/0941309847/0941

auf einen minimalen Dip-Wert eingestellt. Beim Schritt wird die Spur um einen Dip-Wert verschoben, der gleich dem minimalen Dip-Wert plus dem Dip-Abtastmaß (Ad) ist. Die verschobene Spur wird innerhalb des kumulativen Feldes summiert. Beim Schritt 210 wird entschieden, ob noch weitere Dip-Abtastspuren zur Erhöhung übrig sind. Wenn dies der Fall ist, wird die Moveout-Abtastspur erneut um das Dip-Abtastmaß Δd zeitlich verschoben und innerhalb eines anderen kumulativen Feldes beim Schritt summiert.set to a minimum dip value. When stepping, the track is shifted by a dip value that is equal to the minimum dip value plus the dip sampling dimension (Ad). The track moved will be within the cumulative Field totaled. At step 210 it is decided whether there are still more dip-scan tracks left to increase. If this is the case, the moveout scanning track is again shifted in time by the dip scanning amount Δd and summed up within another cumulative field at the step.

Die Schritte 206 bis 210 werden v/iederholt, bis der Dip-Bereich erschöpft ist. Beim Schritt 212 wird eine Entscheidung dahingehend getroffen, ob noch weitere Tiefenpunkte der Hoveout-Abtastspuren für die Verarbeitung übrig sind. Wenn noch solche Tiefenpunkte übrig sind, werden die Schritte 202 bis 212 erneut durchlaufen. Wenn alle Tiefenpunkte bearbeitet worden sind, werden die durchschnittlichen kumulativen Feuer beim Schritt 214 durch Multiplikation mit einem Faktor von 1/SIiASH , der gleich der Zahl der Beiträge für jedes kumulative Feld ist, maßstäblich verändert. Die maßstäblich veränderten Durchschnittspuren werden dann beim Schritt 216 zum Plattenspeicher ausgegeben.Beim Schritt 218 wird entschieden, ob noch weitere Moveout-Abtastspuren zu verarbeiten sind. Wenn noch weitere zu verarbeitende Abtastspuren Y>rhanden sind, werden die Schritte 202 bis 216 erneut durchlaufen, bis alle Moveout-Abtastspuren dip-mässig abgetastet worden sind, woraufhin das Dip-Abtastprogramm beim Schritt 220 endet.Steps 206 to 210 are repeated until the dip area is exhausted. At step 212 a decision is made as to whether there are any further depth points of the Hoveout scan tracks are left for processing. If there are still such depth points left, the Go through steps 202 through 212 again. When all depth points have been processed, the average cumulative fire at step 214 by multiplying by a factor of 1 / SIiASH, which is equal to the number of Contributions for each cumulative field are scaled. The scaled average tracks are then output to disk storage at step 216. At step 218, a decision is made as to whether there are any more Moveout scanning tracks are to be processed. If still there are further scanning tracks Y> to be processed, steps 202 to 216 are run through again until all moveout scanning tracks have been scanned dip-wise, whereupon the dip-scanning program at step 220 ends.

Die Figuren 20 und 21 tragen zum Verständnis des Dip-Abtastaufbauverfahrens bei. Fig.20 zeigt schematischFigures 20 and 21 aid in understanding the dip-scan design process at. Fig.20 shows schematically

309847/0941309847/0941

die Eingangsmoveout-Abtastspuren, die beim Schritt 200 (Fig.19) in das Programm eingegeben werden, Die Moveout-Abtastspuren bestehen für verschiedene Tiefenpunkte aus Moveout-Abtastspuren, die den in Fig.18 dargestellten Spuren gleichen. In dem dargestellten Ausführungsbeispiel von Fig.20 werden Moveout-Abtastspuren bei fünf Tiefenpunkten aus dem Computerspeicher für eine Verarbeitung nach dem in Fig.17 dargestellten Verfahren eingegeben. Die Moveout-Abtastspuren von Fig.20 werden dann entsprechend vorgewählten Dip-Werten verschoben und dann summiert oder gestapelt. Beispielsweise sei angenommen, daß die Moveout-Abtastspuren von Fig.20 zunächst für einen Dip-¥ert von Null verschoben und dann gestapelt werden, damit eine Dip-Abtastung 230 von Fig.21 entsteht. the input moveout scanning tracks which are entered into the program in step 200 (FIG. 19). The moveout scanning tracks consist of moveout scanning tracks for different depth points which are similar to the tracks shown in FIG. In the exemplary embodiment shown in FIG. 20, moveout scanning tracks at five depth points are input from the computer memory for processing according to the method shown in FIG. The moveout scanning tracks of Fig. 20 are then shifted according to preselected dip values and then summed or stacked. For example, it is assumed that the moveout scanning tracks of FIG. 20 are first shifted for a dip ¥ ert from zero and then stacked so that a dip scanning 230 of FIG. 21 is produced.

In allen Fällen wird die mittlere Moveout-Abtastspur 232 beim Tiefenpunkt j nicht in Dip-Richtung verschoben, jedoch werden die Tiefoapunkte zu beiden Seiten des Tiefenpunkts j um unterschiedliche Zeitintervalle verschoben. Damit beispielsweise die Moveout-Abtastungen von Fig.20 um Ad verschoben werden, wird der Tiefenpunkt j+1 zeitlich um ein Intervall Ad nach rechts verschoben, so daß er aufder gestrichelten Linie 234 liegt. Die Moveout-Abtastung beim Tiefenpunkt j+2 wird um 2&d verschoben. Die Moveout-Abtastung j-1 wird um das Intervall Ad nach links verschoben, während die Moveout-Abtastung am Tiefenpunkt j-2 um ein Intervall von 2Ad nach links verschoben wird. Die fünf Hoveout-Abtastungen v/erden dann summiert, damit sich eine Dip.-Abtastung 236 ergibt. In gleicher ¥eise kann die Moveout-Abtastung von Fig.20 in Ausrichtung auf eine gestrichelte Linie 240 verschoben und dann zur !Erzeugung der Dip-Abtastung 242 summiert werdexi.Die Moveout-Abtastungen von Fig.20 werden dann um zusätzliche Bip-Größen verschoben und summiert, damit die weiteren Dip-AbtastungenIn all cases, the middle moveout trace 232 is not shifted in the dip direction at depth point j, however, the deep oa points become on both sides of the deep point j shifted by different time intervals. For example, so that the moveout scans of FIG are shifted by Ad, the depth point becomes j + 1 in time shifted to the right by an interval Ad so that it lies on the dashed line 234. The moveout scan at the depth point j + 2 it is shifted by 2 & d. The moveout scan j-1 moves to the left by the interval Ad shifted while the moveout sample is shifted left by an interval of 2Ad at depth point j-2. the five hoveout samples are then summed to make up a dip. scan 236 results. In the same way, the moveout scanning of Fig. 20 in alignment with a dashed line 240 and then summed to generate the dip sample 242. The moveout samples 20 are then shifted by additional bip sizes and summed up so that the further dip samples

309847/0941309847/0941

von Fig.21 erzeugt werden. Die resultierenden Dip-Abtastdaten, die in Fig.21 dargestellt sind, werden im Plattenspeicher für die Verwendung in der anschließend zu beschreibenden dreidimensionalen Ereignissichtung gespeichert.of Fig.21 can be generated. The resulting dip scan data, which are shown in Fig.21, are in the disk storage for use in the subsequently to be written three-dimensional event viewing saved.

Dreidimensionale EreignisSichtungThree-dimensional event sighting

Die Figuren 22 bis 31 dienen der Beschreibung des Schritts 40 (Fig.3) zur dreidimensionalen Ereignissichtung. Kurz gesagt werden bei diesem Stadium des Verfahrens die mit Hilfe der in Fig.19 dargestellten Technik erzeugten Abtastflächen auf dreidimensionale Ainplitudenextremwerte abgesucht. Die Lage der Extremwerte im Zeit-Moveout-Dip-Bereich wird so interpretiert, daß der räumliche Ort einer reflektierenden Grenzfläche , eines vom Verfahren erzeugten künstlichen Gebildes oder einer nichtppimären Energiekonzentration innerhalb des zu analysierenden Raurar fensters definiert wird.FIGS. 22 to 31 serve to describe step 40 (FIG. 3) for three-dimensional event viewing. Short What is said at this stage of the process is the scanning surfaces generated with the aid of the technique shown in FIG searched for three-dimensional amplitude extreme values. The position of the extreme values in the time moveout dip area is interpreted to mean that the spatial location of a reflective interface, one of the process generated artificial structure or a non-propimeric energy concentration within the Raurar to be analyzed window is defined.

Ein zweidimensionaler Amplitudenextremwert wird immer dann erzeugt, wenn jeder Versetzungsposition auf einer gegebenen CDP-Sainraelaufzeichnung zugeordneten Wellenzüge (wavelets) längs einer der angewendeten Moveout-Hyperbeln liegen und sich phasenrichtig summieren. Ferner wird ein zweidimensionaler Amplitudenextremwert auch dann erzeugt, wenn der i-ten Hovout-Abtastspur aus N aufeinanderfolgenden Tiefenpunkten zugeordnete Wellenzüge auf einer bestimmten Dip-Neigung in einer Richtung liegen und sich phasenrichtig summieren. Jedesmal wenn, ein zweidimensionaler Extremwert sowohl hinsichtlich Moveout als Dip gleichzeitig auftreten, dann ist ein dreidimensionaler Extremwert definiert. Es hat sich gezeigt,daß diese dreidimensionalen Extremwerte mit vernünftiger Gewißheit festgestellt v/erden können, indem nur etwa 20 aufeinanderfolgende logische Tests in einemA two-dimensional amplitude extreme value is always generated when wave trains (wavelets) assigned to each displacement position on a given CDP system recording lie along one of the moveout hyperbolas used and add up in phase. Furthermore, a two-dimensional amplitude extreme value is also generated when the wave trains assigned to the i-th Hovout scanning track of N successive depth points lie on a specific dip inclination in one direction and add up in the correct phase. Each time, a two-dimensional extremum, both in terms moveout as a dip occur simultaneously, then a three-dimensional extreme value is defined. It has been found that these three-dimensional extreme values can be determined with reasonable certainty by only about 20 consecutive logical tests at one time

309847/0941309847/0941

Feld von 21 benachbarten Amplitudenwerten ausgeführt werden. Alle Tests müssen ein positives Ergebnis haben, und ein Ausfall an irgendeiner Stufe der Folge beendet eine weitere Prüfung dieses bestimmten Amplitudenfeldes.Field of 21 neighboring amplitude values. All tests must have a positive result, and failure at any stage in the sequence terminates further testing of that particular amplitude field.

Fig.22 zeigt ein Funktionsflußdiagramm der Schritte des hier beschriebenen Verfahrens zur dreidimensionalen Ereignissichtung. Anfänglich werden drei Moveout-Abtastspurenbeim Schritt 250 in jeder der Dip-Ebenen A und B gelesen. Beim Schritt 252 v/erden drei Moveout-Spuren in einer dritten Dip-Ebene aus dem Speicher gelesen. Während der Initialisierung dieses Programms ist die dritte Dip-Ebene die Dip-Ebene C. Während der anschliessenden wiederholten Durchläufe des Verfahrens werden die übrigen Dip-Ebenen nacheinander beim Schritt 252 gelesen.FIG. 22 shows a functional flow diagram of the steps of the method described here for three-dimensional Event viewing. Initially, there are three moveout scan tracks Read step 250 in each of dip levels A and B. At step 252, three moveout lanes ground read from memory in a third dip level. During the initialization of this program, the third dip level the dip level C. During the subsequent repeated runs of the process, the read remaining dip levels one after the other at step 252.

Beim Schritt 254 wird eine Prüfung in einem 3x3-Zeitkorridor ausgelöst. Der Zeitkorridor umfaßt drei Dip-Ebenen, drei Moveout-Abtastspuren (beginnend mit der Unterseite der Abtastspuren) und einen Zeitkorridor, der sich in Längsrichtung der Spuren erstreckt. Am Schritt 260 wird der Prüfvorgang auf Amplitudenextremwerte längs der Spur längs der Mittelspur innerhalb des Zeitkorridors eingeleitet. Wenn ein möglicher Amplitudenextremwert auf der Spur gefunden ist, wird die Prüfung innerhalb eines in Fig.23 dargestellten vordefinierten Prüfbereichs eingeleitet. Der Prüfbereich ist um den möglichen Amplitudenextremwert A zentriert. Der Prüfbereich enthält 27 Punkte, jedoch werden die Endpunkte der Ecken 256 und 258 jeder Moveout-Zeitebene nicht untersucht, da die Diagonale mit diesen Punkten hinsichtlich ihrer Neigung der Moveout-Neigung entgegengesetzt liegt. Da angenommen wird,daß die Ereignisachse derAt step 254 a test is made in a 3x3 time corridor triggered. The time corridor comprises three dip levels, three moveout scanning tracks (starting with the Underside of the scanning tracks) and a time corridor that extends in the longitudinal direction of the tracks. At the Step 260 begins the testing of amplitude extremes along the track along the center track within the Time corridor initiated. If a possible extreme amplitude value is found on the track, the test is carried out initiated within a predefined test area shown in FIG. The test area is around the possible extreme amplitude value A centered. The test area contains 27 points, but the end points are corners 256 and 258 of each moveout time plane do not investigated, since the diagonal with these points opposes the inclination of the moveout inclination lies. Since it is assumed that the event axis is the

309847/0941309847/0941

Moveout-Neigung nicht entgegengesetzt liegt, werden die Punkte 256 und 258 zur Einsparung von Suchzeit nicht untersucht.Moveout inclination is not opposite, the Points 256 and 258 not examined to save search time.

Wie in Fig.23 dargestellt ist, repräsentieren 21 schwarze Punkte die räumlichen Orte der Ainplitudenwerte. Ein möglicher Extremwert A muß größer oder gleich dem Wert jedes der 20 umliegenden Amplitudenwerte sein, damit er als ein dreidimensionaler Extremwert klassifiziert wird. Die Amplitude des Punkts A wird somit beim Schritt 262 im Hinblick auf die anderen 20 Werte überprüft, und beim Schritt 264 wird festgestellt, ob der ausgewählte "Scheitelpunkt" die Amplitudentests durchläuft. Wenn er die Tests nicht durchläuft, wird beim Schritt 260 ein anderer Extremwert oder Scheitelpunkt ausgewählt, und die Schritte 262 und 264 werden wieder durchlaufen. Wenn beim Schritt 264 ein Amplitudenextremwert festgestellt wird, werden beim Schritt 262 zweidimensionale Interpolationsfilter aufgegriffen und angewendet, die zuvor beim Schritt 132 von Fig.5 angegeben worden sind. Die zweidimensionalen Interpolationsfilter werden zum interpolieren des exakten Amplituden- und Dip-Werts des Extremswerts angewendet.As shown in Fig. 23, 21 represent black Points the spatial locations of the amplitude values. A possible one Extreme value A must be greater than or equal to the value of each of the 20 surrounding amplitude values in order for it to be considered a three-dimensional Extreme value is classified. The amplitude of point A is thus determined at step 262 with respect to the other 20 values are checked and a step 264 determines if the selected "vertex" is the Amplitude tests. If it does not pass the tests, at step 260 another extreme value or Vertex is selected and steps 262 and 264 are repeated. If at step 264 an amplitude extreme value is determined, two-dimensional interpolation filters are picked up and applied in step 262, which were previously indicated in step 132 of FIG. The two-dimensional interpolation filters are used to interpolate the exact amplitude and dip value of the extreme value.

Der Bereich, in dem der Extremwert A von Fig.23 sitzt, kann durch vorbestimmte Punkte definiert v/erden. Die Interpolationsfilter werden verwendet, damit der Ort der interpolierten Punkte präzise festgestellt wird. Es ist erwünscht, den genauen Ort mit der geringsten Anzahl von Interpolationsvorgängen festzustellen; somit wird zur Begrenzung der Zahl der Interpolationsvorgänge eine Logik mit einem logischen Baum (treeing logic) angewendet.The area in which the extreme value A of FIG. 23 sits can be defined by predetermined points. The interpolation filters are used to precisely determine the location of the interpolated points. It is desirable to determine the exact location with the fewest number of interpolations; thus, a logic with a logic tree (treeing logic) is used to limit the number of interpolation processes.

Zur Erklärung des Interpolationsverfahrens wird auf die Figuren 24 bis 30 Bezug genommen. Fig.24 zeigt einen Bereich,Reference is made to FIGS. 24 to 30 to explain the interpolation method. Fig. 24 shows an area

309847/0941309847/0941

UO -UO -

in dem mit Hilfe des oben beschriebenen Amplitudensuchvorgangs der gewünschte Extremwert festgestellt worden ist. In dem definierten Bereich sind A, B, C, D und E Datenpunkte, die innerhalb des Bereichs von +Am und + At liegen.in that using the amplitude search process described above the desired extreme value has been determined. In the defined area are A, B, C, D and E. Data points that are within the range of + Am and + At lie.

Damit der Datenpunkt C der Ort eines Maximums ist, muß sein Wert größer als der der vier anderen Datenpunkte sein. Durch Interpolationsfilter , unter Verwendung der vorbestimmten zweidimensionalen Filter können nun bei dem Zeitinterpolationsabtastmaß in Millisekunden (ΔΔΐ = Δt/4) und dem'Moveout-Interpolationsabtastmaß in Millisekunden pro Abtastung (Δΐη/Δ4) Punkte gefunden werden. Auf diese Weise ergibt sich das in Fig.25 dargestellte Feld in dem A bis E Datenpunkte sind, während die Kreise interpolierte Punkte sind.In order for data point C to be the location of a maximum, its value must be greater than that of the four other data points be. Through interpolation filters, using the predetermined two-dimensional filter, you can now at the time interpolation sample in milliseconds (ΔΔΐ = Δt / 4) and the moveout interpolation sample in milliseconds per sample (Δΐη / Δ4) points found will. This results in the field shown in FIG. 25 in which A to E are data points, while the circles are interpolated points.

Durch Arbeiten mit künstlichen Datenfeldern ist festgestellt worden, daß der Extremwert nur bei bestimmten Punkten der in Fig.25 dargestellten Interpolations« punkte sitzen kann und die Bedingung, daß C einen größeren Wert hat als A, B, D und E noch erfüllt. Das in Fjg.25 dargestellte Originalinterpolationsfeld kann somit auf das in Fig.26 dargestellte Feld reduziert werden, das 23 Punkte enthält.By working with artificial data fields it has been found that the extreme value is only for certain Points of the interpolation points shown in Fig. 25 and the condition that C has a has a greater value than A, B, D and E still fulfilled. The original interpolation field shown in Fig. 25 can thus be reduced to the field shown in Fig. 26, which contains 23 points.

Es ist zu erkennen, daß das in Fig.26 dargestellte Feld hinsichtlich des Ursprungs symmetrisch ist, und eine weitere Reduzierung der Zahl der erforderlichen Interpolationsvorgänge kann dann ausgeführt werden. Durch einen Vergleich von. A mit E kann festgestellt werden, in welcher Hälfte des Datenfeldes der Extremwert liegt. Wenn somitIt can be seen that the field shown in FIG is symmetrical with respect to the origin, and a further reduction in the number of interpolations required can then be executed. By comparing. A with E can be determined in which Half of the data field is the extreme value. If so

309847/0941309847/0941

|Al größer als |E| ist, kann das zu interpolierende Feld nach Fig.27 so definiert werden, daß es nur 14 Punkte enthält. Wenn andererseits | A| kleiner als IEj ist, dann kann das zu interpolierende Feld so wiedergegeben werden, wie es in Fig.28 dargestellt ist.| Al greater than | E | can be the field to be interpolated according to Fig.27 so that it contains only 14 points. On the other hand, if | A | is less than IEj then the field to be interpolated can be reproduced as shown in Fig. 28.

Zur Veranschaulichung des hier zu beschreibenden Verfahrens sei angenommen, daß ein künstliches Datenfeld interpoliert werden soll, das in Fig.29 dargestellt ist. Das in Fig.29 dargestellte Datenpunktfeld enthält fünf Datenpunkte 891» 720, 987, 734 und 733. Da 891 größer als 733 ist, wird entsprechend dem in Fig.27 dargestellten leid nur die obere Feldhälfte interpoliert. Somit enthält das gesamte interessierende Feld 13 Punkte, wie in Fig,29 dargestellt ist.To illustrate the method to be described here, it is assumed that an artificial data field is interpolating which is shown in Fig. 29. The one in Fig. 29 The data point field shown contains five data points 891 » 720, 987, 734 and 733. Since 891 is greater than 733, only the upper one is used, as shown in FIG Field half interpolated. Thus, as shown in FIG. 29, the entire field of interest contains 13 points.

Zur Bestimmung des im Feld von Fig.29 vorhandenen Maximums auf kürzestem Weg wird eine Logik mit logischem Baum angewendet. Die für das vorliegende Beispiel angewendete Logik läßt sich folgendermaßen zusammenfassen:To determine the maximum present in the field of Fig. 29 a logic with a logical tree is applied in the shortest possible way. The logic used for this example can be summarized as follows:

Tabelle IITable II

1. 987 - Datenpunkt C- wird als das laufende Zentrum bezeichnet;1. 987 - data point C- is used as the running center designated;

2. die Position 11 wird für den Wert 968 interpoliert;2. position 11 is interpolated for the value 968;

3. 987>9β8, somit bleibt 987 das laufende Zentrum;3. 987> 9β8, so 987 remains the current center;

4. die Position 7 wird für den Wert 1000 interpoliert;4. Position 7 is interpolated for the value 1000;

5. 987<1000 , somit wird 1000 zum neuen laufenden Zentrum;5. 987 <1000, so 1000 becomes the new running center;

309847/0941309847/0941

6. die Position 8 wird für den Wert 962 interpoliert;6. Position 8 is interpolated for the value 962;

7. 1000>982, somit bleibt 1000 das laufende Zentrum;7. 1000> 982, so 1000 remains the current center;

8. die position 6 wird für den Wert 982 interpoliert;8. position 6 is interpolated for the value 982;

9. 1000 >982;9. 1000> 982;

10. die Position 2 wird für den Wert 987 interpoliert;10. Position 2 is interpolated for the value 987;

11. 1p00>982;11. 1p00> 982;

12. die Position 3 wird für den Wert 968 interpoliert;12. Position 3 is interpolated for the value 968;

13. 1000 >968;13. 1000> 968;

14. die Position 4 wird für den Viert 914 interpoliert;14. position 4 is interpolated for fourth 914;

15. 100O914;15. 100O914;

Ende-1000 ist der Scheitelwert.End-1000 is the peak value.

Somit ist zu erkennen, daß nicht alle in Fig.29 dargestellten Datenpunkte interpoliert werden müssen. Fig.30 gibt an, daß nur sieben der möglichen 13 Punkte zur Interpolation erforderlich waren; es ist die in der Tabelle II definierte Logik angegeben. Da bei dem vorliegenden zveidimsnsinnalen Baumschema nur maximal 13 Punkte interpoliert v/erden müssen, sind nur 13 Filter erforderlich, von denen jedes für eine Gesamtzahl von 65 Orten eine Länge von 5 Wörtern haben muß.It can thus be seen that not all shown in FIG Data points need to be interpolated. Fig. 30 indicates that only seven of the possible 13 points for Interpolation were required; the logic defined in Table II is given. Since with the present zveidimsnsinnalen Tree schema only need to interpolate a maximum of 13 points, only 13 filters are required, of which each must be 5 words long for a total of 65 locations.

309847/0941309847/0941

Die oben beschriebene Filterung wird durch herkömmliche Faltungsverfahren ausgeführt, wobei jeder Filter- und Datenpunkt zur Erzeugung von Ausgangspunkten multipliziert und dann summiert wird. Durch Bestimmung des Maximalwerts sind nun die Zeit- und Moveout-¥erte für das Ereignis festgelegt worden. Nun muß der Dip-Wert für den Extremwert gefunden werden.The filtering described above is carried out by conventional convolution techniques, each filter and Data point is multiplied and then summed to create starting points. By determining the maximum value the time and moveout values for the event have now been set. Now the dip value must be for the extreme value being found.

Beim Schritt 270 von Fig.22 werden die gleichen Filter auf die äußenliegenden Dip-Ebenen angewendet. Fig.31 zeigt die drei Moveout-Zeit-Ebenen, die durch die Größe Dip im Abstand voneinander liegen. Nachdem der Maximumamplitudenpunkt a (k) mit Hilfe des zuvor beschriebenen Verfahrens in der mittleren Dip-Ebene festgestellt worden ist, werden in den zwei benachbarten Moveout-Zeit-Ebenen zusätzliche Maximumamplitudenpunkte a(k-1) und a(k+1) entsprechend dem gleichen Filterverfahren festgestellt. Diese interpolierten drei Punkte sind in Fig.31 dargestellt. Es ergibt sich, daß nur a(k-1) und a(k+1) wegen der Symmetrie der Amplitudenextremwerte im Dip-Abtastraum bstimmt werden müssen. Die drei in Fig.31 dargestellten Punkte werden dann als Funktion der Größe Dip und der Amplitude gemäß Fig.32 dargestellt. Beim Schritt 272 (Fig.22) wird eine hyperbolische Kurve durch die drei Punkte a(k-1), a(k) und a(k+1) gelegt, wie durch die Kurve 274 zweiter Ordnung in Fig.32 angegeben ist. Die Amplitude A der Kurve 274 wird dann beim Schritt 276 von Fig.22 als der Amplitudenv/ert des Ereignisses festgestellt.Der Abstand zwischen dem Punkt A und dem Dip-Wert von a(k) wird als der halbe Dip-Wert des Ereignisses festgestellt.At step 270 of Figure 22 the same filters applied to the outer dip levels. Fig.31 shows the three moveout time levels indicated by the size dip im Distance from each other. After the maximum amplitude point a (k) with the aid of the previously described method in the middle dip level has been determined, additional moveout time levels are added in the two adjacent moveout levels Maximum amplitude points a (k-1) and a (k + 1) corresponding to the same filtering process found. These interpolated three points are shown in Fig. 31. It turns out that only a (k-1) and a (k + 1) because of the symmetry of the extreme amplitude values must be tuned in the dip scanning room. the three points shown in FIG. 31 are then shown as a function of the quantity Dip and the amplitude according to FIG. 32. In step 272 (Fig. 22) a hyperbolic curve is laid through the three points a (k-1), a (k) and a (k + 1), as indicated by the second order curve 274 in Figure 32 is. The amplitude A of curve 274 is then used in step 276 of FIG. 22 as the amplitude value of the event The distance between point A and the dip value of a (k) is considered to be half the dip value of the event established.

Beim Schritt 280 von Fig.22 wird entschieden, ob alle möglichen Extremwerte längs der Spur innerhalb des definierten Korridors geprüft worden sind oder nicht. Wenn dies nicht der Fall war, werden die Schritte 260 bis 276 erneutAt step 280 of FIG. 22 it is decided whether all possible extreme values along the lane within the defined corridor have been checked or not. If this was not the case, steps 260 to 276 are repeated

309847/0941309847/0941

durchlaufen. Wenn alle möglichen Extremwerte längs der Spur geprüft worden sind, wird "beim Schritt 282 festgestellt, ob alle Dip-Ebenen verarbeitet worden sind oder nicht.Wenn dies nicht der Fall war, werden zusätzliche Moveout-Spuren in einer weiteren Dip-Ebene aus dem Speicher beim Schritt 252 gelesen,und die darauffolgenden Schritte werden erneut durchlaufen. Wenn alle Dip~Ebenen verarbeitet worden sind, wird beim Schritt 284 entschieden, ob noch zusätzliche Moveout-Abtastspuren zu verarbeiten sind. Sind noch solche Spuren vorhanden, dann wird die zuletzt verarbeitete Spur beim Schritt 286 umSins erhöht, und drei zusätzliche Moveout-Abtastspuren in zwei Dip-Ebenen werden aus dem Speicher gelesen.run through. When all possible extreme values along the track have been checked, it is determined at step 282, whether or not all dip levels have been processed. If not, additional moveout lanes are added is read from memory in another dip level at step 252 and subsequent steps are repeated run through. When all dip levels have been processed, a decision is made at step 284 whether there are any additional Moveout scanning tracks are to be processed. If there are still such tracks, the last track processed will be used incremented by sins at step 286, and three additional moveout scan tracks read from the memory in two dip levels.

Wenn al3ß Moveout-Abtastspuren verarbeitet worden sind, dann v/erden die Extremereignisses beim Schritt 290 zur Einsparung von ComputerSpeicherplatz zeitlich geordnet, falls es erforderlich ist. Beim Schritt 292 werden die Extremwertereignisse entsprechend einer zeitlich veränderlichen Amplitude sortiert. Bei diesem Sortiervorgang werden die zeitlich geordneten Extremwerte in 32-Ereignisblöcke aufgeteilt. Die gewünschte Anzahl der Ereignisse (Picks) ist festgelegt, und eine eingestellte Zahl von Ereignissen wird aus jeder Ereignisdatei entfernt. Auf diese Weise wird eine eingestellte Menge der Ereignisse mit den höchsten Amplitudenwerten pro Block festgehalten und beim Schritt 294 in den Plattenspeicher geschrieben, worauf das Unterprogramm zur Ereignissichtung beendet ist.When all moveout scan tracks have been processed, then the extreme events are ordered in time at step 290 to save computer memory space, if necessary. In step 292, the extreme value events are changed over time according to a Sorted amplitude. During this sorting process, the time-ordered extreme values are divided into 32 event blocks divided up. The desired number of events (picks) is specified, and a set number of events is removed from each event file. This way, a set amount of events will be the highest Amplitude values per block are recorded and written to disk memory at step 294, whereupon the subroutine for the event viewing has ended.

SegmentaufbauSegment structure

Der Segmentaufbauschritt 42 von Fig.3 wird im Zusammenhang mit den Figuren 33 und 34 erläutert. Nach Fig.33 werden die Ereignisauswahldateien aus dem Speicher 300 beim SchrittThe segment building step 42 of FIG. 3 is described in connection with with FIGS. 33 and 34 explained. According to Fig. 33 the Event selection files from memory 300 at step

309847/0941309847/0941

in Puffer mit zwei oder mehr zusammenhängenden Auswahldateien i, i+1 und dergleichen geladen. Beim Schritt 304 wird das laufende Ereignis A gespeichert, das von der Datei position i ausgeht. Beim Schritt 306 wird das laufende Ereignis Z gespeichert, das von der Dateiposition i+1 ausgeht. Beim Schritt 308 v/erden die zugehörigen logischen Tests hinsichtlich Zeit, Dip, Amplitude und Hoveout angewendet.in buffer with two or more contiguous selection files i, i + 1 and the like loaded. At step 304, the current event A is saved from the file position i. In step 306, the current event Z is stored, which is from the file position i + 1 goes out. At step 308, the associated logic tests for time, dip, amplitude, and hoveout are grounded applied.

Diese logischen Tests umfassen die folgenden Schwellwertbereiche :These logical tests include the following threshold ranges:

Zeit: T_D- At < T < T + at Time: T _D - At <T <T + at

* — q. — ρ ττπτ
wobei gilt T_p = T_A + ■■ * * - q. - ρ ττπτ
where T _p = T _A + ■■ *

Dip: D. -Ad < D „< D„ +Δ dDip: D. -Ad <D "<D" + Δ d

A — Δ— J^ A - Δ - J ^

Amplitude: A. * 1(T^a < Λ < Α_Λ * 10^a Move out: M_A - Am< M₂ <M_A + AmAmplitude: A. * 1 (T ^a <Λ <Α _Λ * 10 ^a Move out: M _A - Am < M ₂ <M _A + Am

Ein einziges Nichteinhalten der oben definierten Grenzwerte beendet ein Segment.A single non-compliance with the limit values defined above ends a segment.

Fig.34 zeigt graphisch die Anwendung des Verlängerungsverfahrens nach der Erfindung.Der Punkt A gibt die Ereigniszeit bezogen auf das Zentrum des ersten Fensters MASH an. Der Punkt Z gibt die Ereigniszeit bezogen auf den Mittelpunkt des zweiten Fensters an. Der Punkt P repräsentiert die Voraussagezeit in der Mitte zwischen den Fenstermittelpunkten , während der Punkt 0. die Ereigniszeit vom Punkt Z zwischen den Fenstermittelpunkten darstellt. Ein Fenster mit der Breite + AtFig. 34 graphically shows the application of the extension method according to the invention. Point A indicates the event time based on the center of the first window MASH. The point Z indicates the event time in relation to the center point of the second window. The point P represents the prediction time in the middle between the window centers, while point 0 is the event time from point Z between the window centers represents. A window with the width + At

309847/0941309847/0941

ist zu beiden Seiten des Punkts P gebildet. Da der Punkt Q innerhalb des Zeitfensters liegt, hält er die oben definierten Zeitschranken ein. Wenn der Punkt P die übrigen Grenzwerte hinsichtlich Amplitude, Moveout und Dip einhält, dann ist eine Verlängerung definiert, und das Segment wird zum nächsten Fenster fortgesetzt.is formed on both sides of the point P. Since the point Q lies within the time window, it holds the ones defined above Time limits. If point P complies with the other limit values with regard to amplitude, moveout and dip, then is an extension is defined and the segment continues to the next window.

Beim Schritt 320 von Fig.33 wird eine Entscheidung hinsichtlich des NichtbeStehens der logischen Tests ausgeführt. Wenn einNichtbestehen nicht auftritt, dann wird beim Schritt 322 die Verlängerung gespeichert. In manchen Fällen erfüllen mehrere Ereignisse in dem zweiten Gatter die oben erwähnten Verlängerungsgrenzwerte, Wenn dies eintritt, wird ein auf die Größen Zeit und Dip bezogener Voraussagefehler formuliert. Das Ereignis mit dem minimalen Fehler wird dann ausgewählt. Der Vorwärtsvoraussagefehler ist folgendermaßen definiert:At step 320 of Fig. 33, a decision is made as to of failing the logical tests. if a failure does not occur, then at step 322 the renewal saved. In some cases several events in the second gate satisfy those mentioned above Extension limits, when this occurs, a prediction error related to the time and dip quantities is formulated. The event with the minimum error is then selected. The forward prediction error is defined as follows:

= T. + IDEC*D. - T_v (26)= T. + IDEC * D. - T _v (26)

wobei gilt: T„ = Voraussagezeit vor der Mitte des nächstenwhere: T "= prediction time before the middle of the next

Dip-Abtastfensters,Dip scanning window,

T₂ = VQraussagezeit hinter der Mitte des letzten Dip-Abtastfensters,T ₂ = prediction time after the middle of the last dip sampling window,

T_A = Zeit des derzeitigen Ereignisses,T _A = time of the current event,

IDEC =s das Aufwärtsmaß von einem SMASH-Fenster zum nächsten, wobei das Aufwärtsmaß die Zahl der nichtverwendeten CDP-Spuren ist.IDEC = s the upward dimension from a SMASH window to the next, where the upward measure is the number of unused CDP lanes.

Der Rückwärtsvoraussagefehler ist folgendermaßen definiert? ⁶B = ^TB * ^TA = ^TZ * ^IDEC*^DZ - ^TAThe backward prediction error is defined as follows? ⁶ B = ^T B * ^T A = ^T Z * ^IDEC * ^D Z - ^T A

309847/0941309847/0941

_fwobei gilt: T_ß = Zeit des letzten Ereignisses. _f where: T _ß = time of the last event.

Der gemeinsame Voraussagefehler ist folgendermaßen definiert:The common prediction error is defined as follows:

⁵F ^{+ 6}B ⁵ F ^{+ 6} B

= δ/ + _δ2 (28)= δ / + _δ 2 (28)

Für die Segmentverlängerung -wird das Ereignis mit dem niedrigsten gemeinsamen Voraussagefehler ausgewählt.For the segment lengthening, the event with the lowest common prediction error selected.

Beim Schritt 324 wird entschieden, ob der Verlängerungspuffer voll ist oder nicht. Wenn er voll ist, wird er beim Schritt 326 zur Zwischenspeicherung in die Verlängerungspuffer 328 entleert. Wenn einer der logischen Tests nicht bestanden wird, wird der Verlängerungspuffer beim Schritt direkt geleert. Wenn der Verlängerungspuffer nicht voll ist, wird beim Schritt 330 entschieden, ob bei i noch weitere Ereignisse vorliegen. Wenn dies der Fall ist, werden die Schritte 304 bis 330 erneut durchlaufen. Wenn keine weiteren Ereignisse vorliegen, wird beim Schritt 332 entschieden, ob noch weitere Ereignisdateien zur Bearbeitung vorhanden sind.Sind noch weitere Dateien vorhanden, dann werden die Puffer mit zwei oder mehr zusätzlichen Ereignisdateien beim Schritt 302 geladen, und die darauffolgenden Schritte werden erneut durchlaufen.At step 324 a decision is made as to whether the extension buffer is full or not. If it is full, it is emptied into extension buffers 328 for intermediate storage at step 326. If one of the logical tests doesn't is passed, the extension buffer is emptied directly at step. If the extension buffer is not full, a decision is made in step 330 as to whether there are still further events at i. If so, the Go through steps 304 through 330 again. If there are no further events, a decision is made in step 332, Whether there are any further event files to be processed. If there are further files, the Buffers loaded with two or more additional event files at step 302, and subsequent steps are run through again.

Wenn alle Ereignisdateien bearbeitet worden sind, wird mit den in den Verlängerungspuffern gespeicherten Daten beim Schritt 340 ein Sortier- und Mischvorgang durchgeführt, und die gemischten Daten werden dann in den Verlange rung spuf fern für die sortierten und gemischten Daten abgespeichert. Beim Schritt 344 werden die Segmente dann formiert, wobei aus der Zahl der Eintragungen für jedes Segment und dergleichen ein Zählerstand gebildet wird. Beim Schritt 346 wird eine Arbeitsdatei erzeugt, und beimWhen all event files have been processed, the data stored in the extension buffers is used a sort and merge operation is performed at step 340, and the merged data is then requested tion is stored remotely for the sorted and mixed data. Then at step 344 the segments are formed, whereby a counter reading is formed from the number of entries for each segment and the like. At the A work file is created in step 346, and at

309847/0941309847/0941

Schritt 348 werden Summendateien vorbereitet. Die Arbeitsdatei enthält Segmente geordnet im Zeit- und Tiefenpunktraum sowie zugehörige Informationen für den Leser. Die Summendateien enthalten über einen bestimmten Zeit-Raum-Bereich verdichtete Segmentinformationen sowie verschiedene, von den Grobsegmentereignissen hergeleitete statistische Daten.Step 348 prepares sum files. The working file contains segments ordered in time and depth point space and related information for the reader. The totals files contain over a certain time-space range Condensed segment information as well as various derived from the coarse segment events statistical data.

Die Summendateien werden so bearbeitet, daß eine' konstante räumliche Abtastung erzielt wird. Vor der Bearbeitung der Summendateien werden Verlängerungen für verschiedene Zeitfenster gebildet, die mit unterschiedliehen räumlichen Abtastungen abgetastet werden.Beim Schritt"348 werden auf diese Weise Operationen ausgeführt, damit die erforderliche konstante räumliche Abtastung durchgeführt werden kann. Die Summendateien werden schließlich als Segmentdateien ausgegeben, die für viele oben erwähnte Anwendungsfälle von Nutzen sind. Die Segmentdateien sind hinsichtlich der Startzeit und des Starttiefenpunkts doppelt geordnet.The totals files are processed in such a way that a 'constant spatial scanning is achieved. Before the processing of the totals files, extensions are made for different time windows which are scanned with different spatial scans. At step "348, operations are performed in this way in order to perform the required constant spatial sampling. Finally, the totals files are output as segment files, which are used by many of the above-mentioned applications Benefits are. The segment files are ordered twice with regard to the start time and the start depth point.

Die Segmentdateien enthalten Segmente mit auf jeder Tiefenpunktbasis abgetasteten Zeit-Amplituden-Moveout-Werten. Alle Arbeitsdatei-Zeitfenstergrenzwerte und die sich ergebenden räumlichen Dezimierungsraten sind in dieser Segmentdatei offensichtlich.The segment files contain segments with on each depth point basis sampled time-amplitude moveout values. All working file time window limits and the resulting spatial decimation rates are in this Segment file obviously.

Die Segmentdateien enthalten in der bevorzugten Ausführung N Aufzeichnungen, mit jeweils einer maximalen Zahl von 64 512 Wörtern. Der Kopfteil der Segmentdateien ist anschließend wiedergegeben:The segment files contain in the preferred embodiment N records, each with a maximum number of 64,512 words. The header of the segment files is then reproduced:

3G984T/Ö9413G984T / Ö941

Wortword Symbolsymbol 11 NSEG :NSEG: CVlCVl IiWDS. :IiWDS. : 33 CW :CW: 44th NWF :NWF:

Tabelle IIITable III Bezeichnungdescription

Zahl der Segmente in der Datei Zahl der 18 Bit-Wörter in der Datei einschließlich des Kopfteils Segmentdatei Name CW = 1111 Zahl der zur Bildung dieser Segmentdatei vereinigten ArbeitsdateienNumber of segments in the file Number of 18-bit words in the file including the header Segment file name CW = 1111 Number of work files combined to form this segment file

Die Segmenteintragung ist eine organisierte Zahlenkette mit insgesamt NWDS-6 Wörtern die auf NSEG Segmente aufgeteilt sind. Die Reihenfolge in der Datei wird vom Segmentstarttiefenpunkt oder vom Segmentstartraumfenster bezüglich früherer Segmente festgelegt. Somit gehen Segmente, die beim Tiefenpunkt i eingeleitet v/erden, jenen Segmenten voraus, die beim Tiefenpunkt i+1 beginnen. Ferner geht das bei i zur Zeit T angegebene Segment j dem bei i zur Zeit T+ ΔΤ angegebenen Segment k voraus.The segment entry is an organized chain of numbers with a total of NWDS-6 words on NSEG segments are divided. The order in the file is determined by the segment starting depth point or the segment starting space window with respect to previous segments. Thus, segments that are initiated at depth point i go to those segments ahead, which begin at the depth point i + 1. Furthermore, the segment j indicated at i at time T corresponds to that at i at time T + ΔΤ specified segment k ahead.

Eine typische Segmenteintragung ist in der anschließenden Tabelle IV wiedergegeben.A typical segment entry is in the following Table IV reproduced.

309847/0941309847/0941

Tabelle IVTable IV

Wortword Symbolsymbol ηη SGN0SGN0 n+1n + 1 SSSSSS n+2n + 2 NSPNSP n+3n + 3 IT(1)IT (1) n+4n + 4 IA(1)IA (1) n+5n + 5 IM(1 )IM (1) n+6n + 6 IT(2)IT (2) n+7n + 7 IA(2)IA (2) n+8n + 8 IM(2).IM (2).

Beschreibungdescription

Segmentzahl; für das Anfangssegment der Datei ist SGN0 =1;Segment number; for the starting segment of the File is SGN0 = 1;

Zahl des Starttiefenpunkts; Zahl der aufeinanderfolgenden, zur Bildung von SGN0 aneinandergefügten Segmentpunkten; Zeit des Ereignisses in ms bei SSS Amplitude des Ereignisses bei SSS Moveout des Ereignisses in ms bei SSS Zeit bei SSS + 1Number of starting depth point; Number of consecutive, to form of SGN0 joined segment points; Time of the event in ms for SSS Amplitude of the event with SSS Moveout of the event in ms with SSS Time at SSS + 1

Amplitude bei SSS+1Amplitude at SSS + 1

Moveout bei SSS+1Moveout at SSS + 1

n+3*NDP ' n+3*NDP+1 n+3*NDP+2n + 3 * NDP 'n + 3 * NDP + 1 n + 3 * NDP + 2

IT(NDP) Zeit bei SSS+NDP-1 IA(NDP) Amplitude bei SSS+NDP-1 IM(NDP) Moveout bei SSS+NDP-1IT (NDP) time for SSS + NDP-1 IA (NDP) amplitude for SSS + NDP-1 IM (NDP) moveout with SSS + NDP-1

309847/0941309847/0941

In der oben angegebenen Tabelle IV ist der Dip-Wert nicht angegeben, da er leicht aus zwei oder mehr aufeinanderfolgenden Zeiteintragungen erhalten werden kann.In Table IV above, the dip value is not because it can easily be obtained from two or more consecutive time entries.

Somit ist zu erkennen, daß .mit Hilfe des hier beschriebenen Verfahrens eine Sammelaufzeichnung für Punkte gemeinsamer Tiefe (CDP-Sammelaufzeichnung) über einen Moveout-Bereich durch zeitliche Verschiebung und Summierung längs hyperbolischer Bahnen abgetastet werden kann. Die summierten Ausgangsspuren werden dann über einen ausgewählten Dip-Bereich abgetastet und in dem definierten Volumen werden Amplitudenextremwerte festgestellt, damit eine gleichzeitige Schätzung der Laufzeit, der Amplitude, des Moveout-Werts und des Dip-Werts eines Ereignisses erzielt werden kann, die mit benachbarten Schätzwerten verglichen und zur Bildung von Reflektorsegmenten zusammengefügt werden kann. las hier beschriebene Verfahren ergibt Reflektorsegmentdaten, die bei nachfolgenden Verarbeitungsschritten äußerst nützlich sind.Thus it can be seen that .with the help of the here described Process a collective record for points of common depth (CDP collective record) over a moveout area can be scanned along hyperbolic orbits by time shifting and summing. The summed output tracks are then scanned over a selected dip range and amplitude extremes are in the defined volume determined so that a simultaneous estimate of the transit time, the amplitude, the moveout value and the dip value of an event can be obtained, which can be compared with neighboring estimates and to form reflector segments can be put together. The method described here yields reflector segment data for the following Processing steps are extremely useful.

Beispielsweise können die gelieferten Segmentdateien in einem Störungskorrekturverfahren verwendet werden, wie es in der der USA-Patentanmeldung Serial Number 254 127 vom 17.Mai 1972 entsprechenden deutschen Patentanmeldung beschrieben ist. Ferner können die mit Hilfe des hier beschriebenen Verfahrens gebildeten Segmentdateien in vorteihafter Weise inf einem Aufbereitungsverfahren angewendet werden, das in einem Dialog zwischen programmiertem Computer und Mensch ausgeführt wird, bei denen Tiefenschichtreflexionsflachen oder Segmente in mehrdimensionalen Zusammenhängen klassifiziert und sortiert werden können. Beschreibungen solcher Dialogsysteme für die Anwendung der hier beschriebenen Segmentdateien finden sich in der Patentanmeldung P 22 64 035.3 sowie in der Patentanmeldung P 22 64 082.0.For example, the segment files supplied can be used in a disturbance correction method, as described in the U.S. Patent Application Serial Number 254 127 dated May 17, 1972 corresponding German patent application is described. Furthermore, with the help of the here described Process-formed segment files are used in an advantageous manner in a processing method that is in a dialogue between programmed computer and human is carried out, in which deep-layer reflection surfaces or segments can be classified and sorted in multi-dimensional contexts. Descriptions of such dialogue systems for the application of the segment files described here in the patent application P 22 64 035.3 and in the patent application P 22 64 082.0.

PatentansprücheClaims

309847/0941309847/0941

Claims (22)

PatentansprücheClaims M.1 Verfahren zur Erzeugung einer Segmentdatei, bei dem Sammel-M.1 Procedure for generating a segment file, in which collective ^S aufzeichnungen für Punkte gemeinsamer Tiefe während der Anwendung verschiedener Moveout-¥erte zur Erzeugung mehrerer Moveout-Abtastspur-Sammelaufzeichnungen bei verschiedenen Tiefenpunkten gestapelt worden sind, dadurch gekennzeichnet, daß die Moveout-Abtastspur-Sammelaufzeichnungen gestapelt werden, während zur Erzeugung eines dreidimensionalen Suchvolumens verschiedene Dip-Richtungen auf sie angewendet werden, und daß Reflektorsegmente bestimmt werden, die sich zwischen benachbarten Amplitudenextremvrerten erstrecken, die innerhalb des dreidimensionalen Suchvolumens festgestellt worden sind. ^ S records for points of common depth have been stacked during the application of different Moveout- ¥ erte to generate multiple Moveout-scanning track collective records at different depth points, characterized in that the Moveout-scanning track collective records are stacked, while different to generate a three-dimensional search volume Dip directions are applied to them and that reflector segments are determined which extend between adjacent amplitude extremes detected within the three-dimensional search volume. 2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daß die Moveout-Abtastspur-Sammelaufzeichnungen über einen Dip-Bereich um den zentralen Tiefenpunkt verschoben werden, damit ein dreidimensiona3.es Suchvolumen mit den Koordinaten Zeit, Moveout und Dip um den zentralen Tiefenpunkt zentriert erzeugt wird.2. The method according to claim 1, characterized in that the moveout scanning track collective recordings over a dip area be shifted around the central depth point, so that a three-dimensional search volume with the coordinates Time, Moveout and Dip are generated centered around the central depth point. 3. Verfahren nach Anspruch 2, dadurch gekennzeichnet,daß bei der Bestimmung der Reflektorsegmente mögliche Amplitudenextremwerte innerhalb eines vordefinierten Zeitkorridors des Suchvolumens festgestellt werden, daß zur Feststellunggültiger Amplitudenextremwerte ein Prüfvorgang innerhalb eines vorbestimmten Prüfbereichs um jeden möglichen Extremwert ausgeführt wird, und daß die Zeit-, Moveout- und Dip-Werte jedes gültigen Amplitudenextremwerts bestimmt werden.3. The method according to claim 2, characterized in that at the determination of the reflector segments possible extreme amplitude values within a predefined time corridor of the search volume that a test process within a predetermined test range is carried out around every possible extreme value, and that the time, moveout and dip values of each valid amplitude extreme value can be determined. 4. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daß bei der Bestimmung der Reflektorsegmente Schätzwerte der Laufzeit, der Amplitude, des Moveout-Werts und des Dip-Werts für die Ainplitudenextremwerte festgestellt werden,4. The method according to claim 1, characterized in that When determining the reflector segments, estimate the transit time, the amplitude, the moveout value and the dip value for the extreme amplitude values are determined, 309847/0941309847/0941 daß jeder Schätzwert mit den Schätzwerten von Extremwerten benachbarter Tiefenpunkte verglichen wird und daß Reflektorsegmente bestimmt werden, die sich zwischen benachbarten Extremwerten erstrecken, deren Schätzwerte innerhalb vorbestimmter Toleranzgrenzen liegen.that each estimated value is compared with the estimated values of extreme values of neighboring depth points and that reflector segments which extend between neighboring extreme values whose estimated values are within predetermined tolerance limits. 5. Verfahren nach Anspruch 3» dadurch gekennzeichnet, daß der vordefinierte Zeitkorridor drei von drei Dip-Ebenen begrenzte Moveout-Abtastspuren enthält, die sich über die Länge der Spuren der Sammelaufzeichnung erstrecken.5. The method according to claim 3 »characterized in that the predefined time corridor contains three moveout scanning tracks, delimited by three dip levels, which extend over the Extend the length of the tracks of the collective recording. 6. Verfahren nach Anspruch 5, dadurch gekennzeichnet, daß die möglichen Amplitudenextremwerte längs der Mitte der drei Moveout-Abtastspuren festgestellt werden.6. The method according to claim 5, characterized in that the possible extreme amplitude values along the center of the three moveout scan traces are detected. 7. Verfahren nach Anspruch 3, dadurch gekennzeichnet, daß der vordefinierte Testbereich von dem Rechteck begrenzt ist, das von acht Datenpunkten auf jeder der drei im Dip-Abstand voneinander entfernten Moveout-Zeit -Ebenen definiert ist.7. The method according to claim 3, characterized in that the predefined test area is delimited by the rectangle, which is defined by eight data points on each of the three moveout time planes, which are separated from one another by the dip distance. 8. Verfahrennach Anspruch 7» dadurch gekennzeichnet, daß ein gültiger Amplitudenextremwert die größte Amplitude innerhalb des Testbereichs aufweist.8. The method according to claim 7 »characterized in that a valid amplitude extreme value has the greatest amplitude within of the test area. 9. Verfahren nach Anspruch 3, dadurch gekennzeichnet, daß bei der Bestimmung der Zeit-Moveout-Werte für jeden gültigen Amplitudenextremwert eine Interpolation auf die genaue Lage des Extremwerts durch Anwendung von zweidimensionalen digitalen Filtern ausgeführt wird.9. The method according to claim 3, characterized in that when determining the time moveout values for each valid Amplitude extreme value an interpolation to the exact position of the extreme value by using two-dimensional digital filtering is performed. 10. Verfahren nach Anspruch 9» dadurch gekennzeichnet, daß der Extremwert in benachbarten Dip-Abtastebenen gemäß der Amplitude dargestellt wird, und daß als der Dip-Wert des10. The method according to claim 9 »characterized in that the Extreme value is represented in neighboring dip-scanning planes according to the amplitude, and that as the dip-value of the 309847/0941309847/0941 Extremwerts eine parabolische Kurve angepaßt wird.Extreme value a parabolic curve is fitted. 11. Verfahren nach Anspruch 3» dadurch gekennzeichnet, daß Reflektorsegmente bestimmt werden, die sich zwischen benachbarten Amplitudenextremwertei] mit Parametern erstrecken, die innerhalb vordefinierter Toleranzen liegen.11. The method according to claim 3 »characterized in that reflector segments are determined which are between adjacent ones Amplitude extremes i] with parameters that are within predefined tolerances. 12. Verfahren nach Anspruch 11, dadurch gekennzeichnet, daß bei jedem Extremwert längs seines geschätzten Dip-Werts Verlängerungen bis zu den Grenzen seines zugeordneten Fensters erzeugt werden, daß die Parameter an den Endpunkten der Verlängerungen in benachbarten Fenstern verglichen und daß die benachbarten Verlängerungen miteinander verbunden werden, wenn die Parameter innerhalb vorbestimmter Toleranzen liegen.12. The method according to claim 11, characterized in that at each extreme value along its estimated dip value Extensions up to the limits of its associated window are generated that the parameters at the end points of the extensions in neighboring windows and that the neighboring extensions are compared with each other be connected if the parameters are within predetermined tolerances. 13. Verfahren nach Anspruch 12, dadurch gekennzeichnet, daß die Parameter Laufzeit-, Dip-, Amplituden- und Moveout-Bereiche sind.13. The method according to claim 12, characterized in that the parameters transit time, dip, amplitude and moveout ranges are. 14.Verfahren nach Anspruch 3, dadurch gekennzeichnet, daß bei der Bestimmung der Zeit-und Moveout-Vferte für jeden gültigen Amplitudenextremwert eine Interpolation auf die genaue Position des Extremswerts durch Anwendung zweidimensionaler digitaler Filter durchgeführt wird.14.Verfahren according to claim 3, characterized in that in determining the time and moveout values for each valid amplitude extreme value an interpolation to the exact position of the extreme value by applying two-dimensional digital filter is performed. 15.Verfahren nach Anspruch 14, dadurch gekennzeichnet, daß die Extremwerte in benachbarten Dip-Abtastebenen entsprechend der Amplitude aufgezeichnet werden, daß durch die Ainplitudenv/erte der Extremwerte in benachbarten Dip-Abtastebenen eine parabolische Kurve gelegt wird und daß der Scheitelpunkt der parabolischen Kurve als der Dip-Wert des Extremwerts bezeichnet wird.15.Verfahren according to claim 14, characterized in that the extreme values are recorded in adjacent dip-scanning planes according to the amplitude that is determined by the amplitude values of the extreme values in adjacent dip-scan planes a parabolic curve is laid and that the vertex of the parabolic curve is referred to as the extreme dip value. 309847/0941309847/0941 16. Verfahr en nach Anspruch 14, dadurch gekennzeichnet, daß zum Suchen der interpolierten Position des Extremwerts eine Logik mit einem logischen Baum angewendet wird.16. Verfahr s according to claim 14, characterized in that logic with a logical tree is applied to search for the interpolated position of the extreme value. 17- Anordnung zur Erzeugung einer Segmentdatei aus einer digitalen Sammelaufzeichnung seismischer Spuren mit Tiefenschichtredundanz unter Anwendung des Verfahrens nach einem der vorhergehenden Ansprüche, gekennzeichnet durch Einrichtungen zum Programmieren eines digitalen Computers zum Stapeln der seismischen Spuren während der Anwendung unterschiedlicher Zeitverschiebungen und Dip-Werte auf die seismischen Spuren, Einrichtungen zum Programmieren des Computers zum Herausgreifen von Anrplitudenextremverten aus den gestapelten seismischen Spuren zur Erzeugung von Schätzungen für die Laufzeit-, Amplituden-, Moveout- und Dip-Verte eines Ereignisses, und Einrichtungen zum Programmieren des Computers zur Bestimmung von Reflektorsegmenten, die sich zwischen benachbarten Schätzungen von gleichen Werten erstrecken.17- Arrangement for generating a segment file from a digital collective record with seismic traces Depth layer redundancy using the method according to one of the preceding claims, characterized by devices for programming a digital Computers to stack the seismic traces while applying different time shifts and dip values on the seismic traces, facilities for programming the computer to pick out extreme amplitudes from the stacked seismic traces to generate estimates for the transit time, amplitude, Moveout and Dip values of an event, and facilities to program the computer to determine reflector segments that are between adjacent estimates of the same values. 18. Anordnung nach Anspruch 17, gekennzeichnet durch Einrichtungen zur zeitabhängigen Filterung und Entzerrung der seismischen Spuren vor der Stapelung.18. Arrangement according to claim 17, characterized by devices for time-dependent filtering and equalization of the seismic traces before stacking. 19. Anordnung nach Anspruch 17, dadurch gekennzeichnet, daß die Programmierungseinrichtungen zum Herausgreifen von Amplitudenextreisverten derart ausgestaltet sind, daß sie Aiuplitudenextremwerte innerhalb eines vordefinierten Zeitkorridors eines dreidimensionalen Suchvolumens feststellen, das durch das Stapeln definiert ist.19. The arrangement according to claim 17, characterized in that the programming devices for picking out amplitude extra verts are designed in such a way that they achieve extreme values determine within a predefined time corridor of a three-dimensional search volume that is defined by the stacking. 309847/094 1309847/094 1 20. Anordnung nach Anspruch 19, gekennzeichnet durch Einrichtungen zum Anwenden zweidimensionaler Filter zum Interpolieren der genauen Position der Extremwerte innerhalb des Zeitkorridors.20. Arrangement according to claim 19, characterized by devices to apply two-dimensional filters to interpolate the exact position of the extreme values within the time corridor. 21. Anordnung nach Anspruch 20, gekennzeichnet durch Einrichtungen zum -Anwenden der Filter auf "benachbarte Dip-Abtastebenen zur Entwicklung zusätzlicher Dip^ Orte der Extremwerte, Einrichtungen zum Darstellen der Dip-Orte entsprechend der Amplitude, Einrichtungen zum Legen einer Parabel durch die Amplitudenwerte und Einrichtungen zum Bestimmen des Scheitelpunkts der Parabel als Dip-Wert für den Extremwert.21. Arrangement according to claim 20, characterized by means for applying the filter to "adjacent Dip scan planes to develop additional Dip ^ Locations of the extreme values, devices for displaying the dip locations according to the amplitude, devices for laying a parabola through the amplitude values and means for determining the vertex of the Parabola as a dip value for the extreme value. 22. Anordnung nach Anspruch 17, gekennzeichnet durch Einrichtungen zur' Bildung einer geordneten Liste von Parametern der Reflektorsegmente und Magnetspeichereinrichtungen zum wiederauffindbaren Speichern der geordneten Liste.22. Arrangement according to claim 17, characterized by devices to create an ordered list of parameters the reflector segments and magnetic storage devices for the retrievable storage of the ordered List. '309847/094 1'309847/094 1