DE2460040A1

DE2460040A1 - Strata spacing seismic vel detmn for sloping inter-strata interface - in automatic process involving corrected reflections from at least three profile lines

Info

Publication number: DE2460040A1
Application number: DE19742460040
Authority: DE
Inventors: William Harlold Ruehle
Original assignee: Mobil Oil Corp
Current assignee: ExxonMobil Oil Corp
Priority date: 1974-12-19
Filing date: 1974-12-19
Publication date: 1976-06-24

Abstract

Prodn. of a recording of the actual strata spacing velocity (i.e. the velocity of the seismic waves to and from the underground reflectors) dimensionless groups for the ground in the presence of a sloping stratum, using an automatic technique in which (a) seismic reflexions are obtained from a series of >=3 seismic profile lines crossing in a common surface point, (b) the reflexion times for each reflexion, which are all common to the seismic profile lines are ascertained, (c) each reflection is corrected with respect to the inclination for those common to all profile lines in order to obtain the root mean square (r.m.s.) velocities from which a continuous and smoothed curve of r.m.s. velocity is produced, (d) an arbitrarily selected sloping Snellius Law - Layer model is fitted to the smoothed r.m.s. velocity curve so as to obtain the actual strata spacing velocities. Used in oil prospecting etc.

Description

Verfahren zum Bestimmen der Intervallgeschwindigkeit bei Gegenwart einer schrägliegenden Einbettung Die Erfindung bezieht sich auf die Bestimmung der echten Schichtintervallgeschwindigkeit aus einer Vielzahl Gruppen seismischer Aufzeichnungsspuren, insbesondere auf die Bestimmung der echten Schichtintervallgeschwindigkeit bei Gegenwart einer schrägliegenden Einbettung.Method for determining the interval speed in the presence an inclined embedding The invention relates to the determination of the real shift interval velocity from a large number of groups of seismic recording tracks, in particular on the determination of the real shift interval speed in the present an inclined embedding.

Bei der seismischen Exploration ist die genaue Geschwindigkeitsbestimmung für die Verarbeitung und Interpretation seismischer Daten sehr wichtig. Die Bestimmung der Geschwindigkeitscharakteristik ist beschrieben in "Seismic Velocities From Subsurface Measurements", C. H. Dix, GEOPHYSICS, Vol. 20, Seiten 68-86, 1955.In seismic exploration, the exact speed determination is very important for the processing and interpretation of seismic data. The determination the speed characteristic is described in "Seismic Velocities From Subsurface Measurements ", C. H. Dix, GEOPHYSICS, Vol. 20, pp. 68-86, 1955.

Eine verbreitet benutzte Feldoperationsmethode zum Erhalt seismi scher Aufzeichnungsspuren zwecks Geschwindigkeitsbestimmungen ist die Punkt-gemeinsamer-Tiefe-Methode (kurz: PGT-Methode) zum Erhalt einer Mehrfachüberdeckung unterirdischer Reflexionspunkte. Die Ankunftszeit von Aufzeichnungsspur zu Aufzeichnungsspur Jeder Reflexion in einer Gruppe von Seismogrammen bezüglich eines Punktes gemeinsamer Tiefe (kurz: PGT-Seismogrammgruppe) ändert sich nach einer hyperbolischen Funktion, die allgemein als Normalaustritt (normal moveout) bezeichnet wird. Durch eine Korrektur jeder Reflexion bezüglich ihres Normal austrittes kann eine jede der PGT-Aufzeichnungsspuren überlagert werden, um das Reflexionssignal zu verstärken. PGT-Überlagerung ist ebenso bei schrägliegenden wie bei flachverlaufenden unterirdischen Einbettungen anwendbar. Die Beziehung zwischen der scheinbaren Geschwindigkeit, die zur Korrektur der CDP-Daten für Normalaustritt vor der Überlagerung benutzt wird, und der echten Schichtintervallgeschwindigkeit, mit der die seismischen Wellen zu und von den unterirdischen Reflektoren laufen, ist beschrieben in dem Artikel "Apparent Velocity From Dipping Interface Reflections", Franklyn K. Levin, Velocity Symposium, Houston, Texas, 1969.A widely used field operation method for obtaining seismic The point-shared-depth method is used to record traces for the purpose of determining speed (short: PGT method) to obtain a multiple coverage of underground reflection points. The arrival time from recording track to recording track of each reflection in one Group of seismograms relating to a point of common depth (short: PGT seismogram group) changes according to a hyperbolic function, which is generally called normal exit (normal moveout) is called. By correcting any reflection regarding of their normal exit, each of the PGT recording tracks can be superimposed, to amplify the reflection signal. PGT overlay is the same with inclined ones as applicable for shallow underground embeddings. The relationship between the apparent speed used to correct the CDP data for normal exit is used before the overlay, and the real slice interval velocity, with which the seismic waves travel to and from the underground reflectors, is described in the article "Apparent Velocity From Dipping Interface Reflections", Franklyn K. Levin, Velocity Symposium, Houston, Texas, 1969.

Entsprechend der Erfindung werden seismische CDP-Reflexionen erhalten aus einer Vielzahl von wenigstens drei seismischen Profillinien, d.ie sich in einem gemeinsamen Oberflächenpunkt schneiden. Die mittlere Reflexionszeit für jene Reflexionen, die sämtlichen Profillinien gemeinsam sind, wird bestimmt.In accordance with the invention, CDP seismic reflections are obtained from a plurality of at least three seismic profile lines, i.e., which are in one intersect common surface point. The mean reflection time for those reflections the all profile lines are common is determined.

Die mittlere scheinbare Geschwindigkeit und die mittlere Schräglage für åede-Reflexion längs der Profillinien wird.The mean apparent speed and the mean lean angle for åede reflection along the profile lines.

gleichfalls bestimmt. Die mittlere Schräglage und der Azimuth jeder Profillinie werden dann zur Bestimmung der echten Deklination (strike) und Schräglage jener Reflexionen verwendet, die allen Profillinien gemeinsam sind. Jene mittlere scheinbare Geschwindigkeit wird dann bezüglich der echten Deklination und Schräglage korrigiert, um eine VMittelwertquadrat- Geschwindigkeit (RMS velocity) für jede Reflexion längs jeder Profillinie zu bestimmen. Eine durchschnittliche telwertquadrat- Geschwindigkeit wird dann für jene Reflexionen bestimmt, die allen Profillinien gemeinsam sind. Die durchschnittlichen gMittelwertquadrat-Geschwindigkeiten werden dann über den durchschnittlichen Reflexionszeiten für jene Reflexionen aufgetragen, die allen Profillinien gemeinsam sind. Das Diagramm wird kontinuierlich gemacht, indem jeder (Mittelwertquadrat-Geschwindigkeitspunkt durch eine gerade Linie verbunden wird. Die kontinuierliche {Mittelwertquadrat-Geschwindigkeitsktrve wird dann mit Hilfe eines Glättungsfilter-Operator geglättet. Die echte Intervallgeschwindigkeit wird dann durch einen Iterationsprozeß bestimmt, der ein willkürliches, schrägliegendes Schichtenmodell entsprechend dem Snellius-Gesetz an die geglättete Kurve der Mittelwertquadrat-Geschwindigkeiten anpaßt.also determined. The mean skew and azimuth of each profile line are then used to determine the true declination (strike) and skew of those reflections that are common to all profile lines. That mean apparent speed is then corrected for true declination and bank angle by one V mean square Determine the velocity (RMS velocity) for each reflection along each profile line. An average telvalue square Speed is then determined for those reflections that are common to all profile lines. The average g mean square velocities are then plotted over the average reflection times for those reflections that are common to all profile lines. The chart is made continuously by everyone (Mean square speed point connected by a straight line. The continuous {Mean square velocity curve is then smoothed using a smoothing filter operator. The real interval speed is then determined by an iterative process, which an arbitrary, inclined layer model according to Snellius law to the smoothed curve of the Mean square velocities adapts.

Die Erfindung ist in den Ansprechen gekennzeichnet und in der Zeichnung im einzelnen beschrieben; es zeigen: Fig. 1 eine Feldmethode zum Erhalt seismischer PGT-Aufzeichnungsspuren, Fig. 2 die Ankunftszeiten d.er Reflexionen an den seismischen Aufzeichnungsspuren, wie diese auf die in Fig. 1 dargestellte Weise erhalten worden sind, Fig. 3 die Beziehung einer schrägliegenden Grenzfläche zur Erdoberfläche, Fig. 4 ein Überkreuzungsmuster seismischer Profillinien auf der Erdoberfläche, Fig. 5A und B Flußdiagramme für den Verfahrensablauf, Fig. 6A, 6B, 7A, 7B und 7C Matrizen zum Speichern der seismischen Daten, wie diese entsprechend. den Schritten nach Fig. 5A und B verarbeitet worden sind, Fig. 8 Vektor-Diagramm von Deklinations- oder Neigungsgrößen, wie diese mit d.en Schritten nach Fig. 5A und B erhalten werden, Fig. 9 die Richtung der Deklination der schrägliegenden Grenzfläche bezüglich der Profillinien in Fig. 4 Fig. 10 ein Diagramm der geschätzten quadratischen Geschwindigkeitsmittelwerte, wie diese erhalten worden sind, entsprechend den Schritten nach Fig. 5A und B, Fig. 11 einen Glättungsfilter-Operator zur Verwendung beim Glätten der in Fig. 10 aufgetragenen Kurve, Fig. 12 ein Diagramm von Neigung und Deklination, wie diese entsprechend den Schritten nach Fig. 5A und 5B erhalten worden sind, Fig. 13 ein Diagramm der vMittelwertquadrat-Geschwindigkeit, der neigungskorrigierten Intervallgeschwindigkeit und der mittleren Geschwindigkeit, wie diese entsprechend den Schritten nach Fig. 5A und B erhalten worden sind, und Fig. 14A, 14B und 14C Flußdiagramme zum Durchführen der Schritte nach Fig. 5A und 5B.The invention is characterized in the responses and described in detail in the drawing; 1 shows a field method for obtaining seismic PGT recording tracks, FIG. 2 shows the arrival times of the reflections on the seismic recording tracks as obtained in the manner shown in FIG Interface to the earth's surface, Fig. 4 shows a crossover pattern of seismic profile lines on the earth's surface, Figs. 5A and B are flow charts for the process sequence, Figs. 6A, 6B, 7A, 7B and 7C are matrices for storing the seismic data, as these correspond. 5A and B have been processed; with respect to the profile lines in FIG. 4, FIG. 10 is a diagram of the estimated root mean square velocity values, as they have been obtained, corresponding to the steps according to FIGS. 5A and B, FIG. 11 shows a smoothing filter operator for use in smoothing those plotted in FIG Curve, FIG. 12 is a diagram of inclination and declination as obtained in accordance with the steps according to FIGS v mean square velocity, the incline corrected interval speed and the mean speed obtained in accordance with the steps of Figs. 5A and B; and Figs. 14A, 14B and 14C are flow charts for performing the steps of Figs. 5A and 5B.

In Fig. 1 ist ein Verfahren zum Erhalt einer Gruppe seismischer Aufzeichnungsspuren an einem Punkt gemeinsamer Tiefe (kurz: PGT-Gruppe) dargestellt. Seismische Energie, die an jedem einer Vielzahl im Abstand längs einer seismischen Profillinie voneinander gelegenen Detonations- oder Schußpunkten S1-Sn erzeugt wird, wird von einer unterirdischen Grenzfläche an einem Punkt gemeinsamer Tiefe reflektiert und an einer Vielzahl längs der seismischen Profillinie im Abstand voneinander gelegenen Empfangsstellen R1-Rn zum Erzeugen einer PGT-Gruppe seismischer Aufzeichnungsspuren empfangen. Eine derartige PGT-Gruppe von Reflexionen Z1-Zn ist in Fig. 2 dargestellt. Die Reflexionssignale Z1-Zn, die von der PGT-Grenzfläche empfangen werd.en, sind als längs der Kurve W liegend angegeben. Die Kurve W stellt die AnkunStszeit von Aufzeichnungsspur zu Aufzeichnungsspur für jedes der Reflexionssignale in der Seismogramm-Schar dar. Diese Kurve ist durch die hyperbolische Funktion bestimmt: TX2 = T02 + X2/Va2 (1) Hierin bedeuten: Tx die Reflexionszeit auf einer bestimmten seismischen Aufzeichnungsspur T0 die Zeit jener Reflexion auf einer idealisierten seismischen Aufzeichnungsspur mit einer Reflexionsstelle direkt unterhalb einer Schußstelle X die horizontale Entfernung zwischen der Schußstelle und dem Empfänger, der die betreffende seismische Aufzeichnungsspur erzeugt, und Va die mittlere oder scheinbare akustische Geschwindigkeitskenngröße der Schicht, durch die die seismische Energie läuft.Referring to Fig. 1, there is a method of obtaining a group of seismic recording tracks shown at a point of common depth (short: PGT group). Seismic energy, those at each of a plurality spaced from one another along a seismic profile line located detonation or shot points S1-Sn is generated by an underground Boundary surface reflected at one point of common depth and longitudinally at a plurality the seismic profile line at a distance from each other receiving points R1-Rn to create a PGT group of seismic logging tracks. Such a one PGT group of reflections Z1-Zn is shown in FIG. The reflection signals Z1-Zn received from the PGT interface are shown along the curve W. indicated lying down. The curve W represents the arrival time of the recording track Represents the recording track for each of the reflection signals in the seismogram family. This curve is determined by the hyperbolic function: TX2 = T02 + X2 / Va2 (1) Herein mean: Tx the reflection time on a certain seismic Recording track T0 the time of that reflection on an idealized seismic Recording track with a point of reflection directly below a point of fire X is the horizontal distance between the point of fire and the recipient who received the relevant seismic trace is generated, and Va is the mean or apparent acoustic speed parameter of the layer through which the seismic energy runs.

Die Zeitverschiebung von Aufzeichnungsspur zu Aufzeichnungsspur d.er Reflexionen von der PGT-Grenzfläche ist allgemein bekannt als Normalaustritt und. ist spezifiziert durch die Funktion: Bezüglich einer weiteren Beschreibung der Bestimmung des Normalaustrittes wird verwiesen auf SEISMIC PROSPECTING FOR OIL von C. Hewitt Dix, 1952, Section 8.2.3, Seiten 134-137, Die scheinbare Geschwindigkeit Va, die zur Korrektur für den Normalaustritt benutzt wird, ist jedoch nicht dieselbe wie die echte Schichtintervallgeschwindigkeit Vi. Die Beziehung zwischen der scheinbaren Geschwindigkeit Va und der Intervallgeschwindigkeit Vi ist bestimmt durch die Richtung der seismischen Profillinie gegenüber der Neigung der schrägliegenden unterirdischen Grenzfläche. Entsprechend Fig. 3 ist der Neigungswinkel zwischen dem Lot auf d.er schrägverlaufenden Grenzfläche und dem Lot auf der Erdoberfläche durch den Winkel ffi dargestellt. Der Deklinationswinkel zwischen der seismischen Profillinie und der Neigungslinie längs der Erdoberfläche ist durch den Winkel @ dargestellt. Der Ausdruck für die echte Schichtintervallgeschwindigkeit kann deshalb angegeben werden durch die Funktion: Bezüglich der Ableitung von Gleichung (3) sei verwiesen auf die Arbeit "Apparent Velocity From Dipping Interface Reflections" von Franklyn K. Levin, Velocity Symposium, Houston, Texas, 1969.The time shift from record track to record track of the reflections from the PGT interface is commonly known as normal leakage and. is specified by the function: For a further description of the determination of the normal outlet, reference is made to SEISMIC PROSPECTING FOR OIL by C. Hewitt Dix, 1952, Section 8.2.3, pages 134-137. However, the apparent velocity Va, which is used to correct for the normal outlet, is not the same as the real slice interval velocity Vi. The relationship between the apparent velocity Va and the intermittent velocity Vi is determined by the direction of the seismic profile line versus the slope of the sloping subterranean interface. According to FIG. 3, the angle of inclination between the perpendicular on the sloping interface and the perpendicular on the earth's surface is represented by the angle ffi. The angle of declination between the seismic profile line and the slope line along the earth's surface is represented by the angle @. The expression for the real shift interval speed can therefore be given by the function: Regarding the derivation of equation (3), reference is made to the work "Apparent Velocity From Dipping Interface Reflections" by Franklyn K. Levin, Velocity Symposium, Houston, Texas, 1969.

Insoweit die Richtung der Schräglage unbekannt ist, wird wohl kaum die Empfänger-Linie mit dem Azimuth d.er Schräglage bzw.Insofar as the direction of the incline is unknown, it will hardly be the receiver line with the azimuth of the lean or

Neigung zusammenfallen. Entsprechend der Erfindung ist es möglich, Daten zu erhalten, die für die echte Schichtintervallgeschwindigkeit in Gebieten repräsentativ sind, die schrägliegende Formationen haben und bei denen der Deklinationswinkel unbekannt ist, indem Gleichung (3) nach diesen Unbekannten aufgelöst wird.Slope coincide. According to the invention it is possible Obtain data relevant for the real shift interval velocity in areas are representative that have sloping formations and where the angle of declination is unknown by solving equation (3) for these unknowns.

Da drei Unbekannte vorhanden sind, ist eine Kreuzungsstelle von drei oder mehr seismischen Profillinien erforderlich. Entsprechend einer Ausführungsform der Erfindung werden eine Vielzahl seismischer Profillinien auf der Erdoberfläche in einem Netzmuster entsprechend Fig. 4 ausgelegt. Das d.urch diese Linien erzeugte Netzmuster liefert eine Vielzahl Überkreuzungspunkte dreier Profillinien. Beispielsweise haben d.ie seismischen Profillinien A, B, C einen Kreuzungspunkt a. Seismische Reflexionen, die von diesen drei Profillinien A, B und C erhalten werden, werden zur Lösung von Gleichung (3) für den Kreuzungspunkt a benutzt. In ähnlicher Weise können seismische Reflexionen von den drei Profillinien, die jeden der übrigen Kreuzungspunkte des Netzmusters bilden, zur Lösung von Gleichung (3) nach diesen Kreuzungspunkten verwendet werden.Since there are three unknowns, there is an intersection of three or more seismic profile lines required. Corresponding one Embodiment of the invention are a plurality of seismic profile lines the earth's surface laid out in a network pattern according to FIG. The d.by The mesh pattern produced by these lines provides a multiplicity of crossing points of three Profile lines. For example, the seismic profile lines A, B, C have a point of intersection a. Seismic reflections obtained from these three profile lines A, B and C. are used to solve equation (3) for the intersection point a. In Similarly, seismic reflections can be seen from the three profile lines that define each of the remaining intersection points of the mesh pattern, to solve equation (3) can be used after these crossing points.

Fig. 5A und 5B sind Elußdiagramme einer auf einem Allzweckrechner entsprechend der Erfindung durchgeführten Methode zur Bestimmung der echten Schichtintervallgeschwindlgkeit Vi, des Neigungswinkels ld und Deklinationswinkel Q für jeden Kreuzungspunkt dreier sich schneidender seismischer Profillinien. Der Kreuzungspunkt a seismischer Profillinien A, B und C wird zu Erläuterungszwecken benutzt. Der erste Schritt ist, wie bei 20 im Flußdiagramm angegeben, aus d.en seismischen PGT-Gruppen für die drei Profillinien A, B und C die Zeit, die Entfernung, die scheinbare Geschwindigkeit und die Amplitude jeder seismischen Reflexion auf jeder seismischen PGT-Gruppe zu bestimmen. Eine solche Bestimmung kann ausgeführt werden durch fortlaufende Geschwindigkeitsbestimmungsmethoden, wie diese im einzelnen beschrieben ist in der US-PS 3,651,451.Figures 5A and 5B are flow charts of one on a general purpose computer Method carried out according to the invention for determining the real shift interval speed Vi, the angle of inclination ld and the angle of declination Q for each intersection of three intersecting seismic profile lines. The intersection of a seismic profile lines A, B, and C are used for explanatory purposes. The first step is like at 20 given in the flow chart, from the seismic PGT groups for the three profile lines A, B and C are the time, the distance, the apparent speed and the amplitude determine each seismic reflection on each seismic PGT group. One such determination can be carried out by continuous speed determination methods, as described in detail in U.S. Patent 3,651,451.

Der nachsta Schritt im Verfahren ist, wie bei 21 angegeben, die Daten des Schritts 2Q in drei getrennten Speichereinheiten zu speichern, von denen Je einer für die drei seismischen Profillinien A, B und C vorgesehen ist. dede Speichereinheit umfaßt zwei getrennte Matrizen die erste Matrix ist dreidimensional und hat die Zeit-Daten in einem Abschnitt, die Daten der scheinbare Geschwindigkeit in einem zweiten Abschnitt und die Amplitude in einem dritten Abschnitt gespeichert. Ein Beispiel für eine solche Matrix für die seismische Profillinie A ist in Fig. 6A dargestellt. Lediglich zu Erläuterungszwecken ist eine 6 x 6'x 3 - Matrix dargestellt, bei der nur die Zeit-Dateneingänge spezifiziert sind. Beispielsweise ist der Zeiteingang T11 die Zeit der Reflexion Z1 der P(1T-(iruppe 1 der Profillinie A.The next step in the method, as indicated at 21, is the data of step 2Q in three separate storage units, each of which one is provided for the three seismic profile lines A, B and C. dedicated storage unit comprises two separate matrices the first matrix is three-dimensional and has the Time data in one section, apparent speed data in one second section and the amplitude stored in a third section. A An example of such a matrix for seismic profile line A is shown in FIG. 6A shown. A 6 x 6'x 3 matrix is shown for explanatory purposes only, in which only the time data inputs are specified. For example is the time input T11 the time of the reflection Z1 of the P (1T- (i group 1 of the profile line A.

Die zweite Matrix ist eindimensional zur Speicherung der Entrernungsdaten und ist in Fig. 6B dargestellt. Beispielsweise ist die Entfernung D1 die Entfernung in 0,3 m zwischen der PGiT-Gruppe 1 und der PG'T-(1ruppe 2 für eine Reflexion der Profillinie A.The second matrix is one-dimensional for storing the evacuation data and is shown in Figure 6B. For example, the distance D1 is the distance 0.3 m between the PGiT group 1 and the PG'T- (1 group 2 for a reflection of the Profile line A.

Der nächste Schritt im Verfahren ist wie bei 22 angegeben, die Bestimmung der mittleren Zeit T, der durchschnittlichen scheinbaren Geschwindigkeit Va und der mittleren Neigung > für Jede Reflexion auf einer Profillinie. Zuerst werden die mittlere Zeit und die mittlere scheinbare G"eschwindigkeit Jeder Reflexion bestimmt, entsprechend den folgenden Gleichungen Hierin bedeuten N= die Anzahl der PGT-Gruppen-Zeiteingänge für eine gegebene Reflexion j= Reflexion und.The next step in the method, as indicated at 22, is to determine the mean time T, mean apparent velocity Va, and mean slope> for each reflection on a profile line. First, the mean time and the mean apparent velocity of each reflection are determined according to the following equations Here, N = the number of PGT group time inputs for a given reflection j = reflection and.

k= PGT-Gruppe Zeit für die Beispielsweise ist die mittlere/Reflexion Z1 der PGT-Gruppen 1 bis 6 auf Linie A in Fig. 6A T(zl)(A) T11 + T12 + T13 + T14 +- T15 + T16 (6) 6 In ähnlicher Weise ist die mittlere scheinbare Geschwindigkeit für die Reflexion Zi der PGT-Gruppen 1 bis 6 auf Linie A Va(Z1)(A) = V11 +V12 + V13 + V14 + V15 + V16 Die mittlere Neigung längs Jeder Reflexion ist dann bestimmt durch den folgenden Ausdruck: hierin bedeuten Dn = Entfernung zwischen benachbarten PGT-Gruppen für eine Reflexion längs einer gegebenen Profillinie.k = PGT group time for the For example, the mean / reflection Z1 of PGT groups 1 to 6 on line A in Fig. 6A is T (zl) (A) T11 + T12 + T13 + T14 + - T15 + T16 (6 ) 6 Similarly, the mean apparent velocity for the reflection Zi of PGT groups 1 to 6 on line A is Va (Z1) (A) = V11 + V12 + V13 + V14 + V15 + V16 The mean slope along each reflection is then determined by the following expression: herein Dn = distance between adjacent PGT groups for a reflection along a given profile line.

Beispielsweise ist die durchschnittliche Neigung > für die Reflexion Z1 der PGT-Gruppen 1 bis 6 von Linie A Vor dem-8bergang auf die nachfolgenden Schritte des Verfahrens wird eine Bestimmung hinsichtlich der Annehmbarkeit jeder Reflexion als eine gültige Reflexion für d.ie weitere Verarbeitung durchgeführt. Diese Bestimmung ist ein Prozeß mit zwei Schritten.For example, the average slope> for the reflection Z1 of PGT groups 1 through 6 of line A. Before proceeding to the subsequent steps of the method, a determination is made as to the acceptability of each reflection as a valid reflection for further processing. This determination is a two step process.

Zuerst wird die mittlere Amplitude A entsprechend der folgenden Gleichung bestimmt Beispielsweise ist die mittlere Amplitude der Reflexion Z1 der PGT-Gruppen 1 bis 6 auf Linie A gegeben durch A +A + A (Z1)(A) = A11 + 12 + A1 +-A14 + A15 + A16 (11) Die Reflexion beispielsweise Z1 auf Linie A wird dann als akzeptabel für die weitere Verarbeitung betrachtet, wenn die nachstehenden Bedingungen erfüllt sind W(Zi)(A) > Q1 (12) N > Q2 (13) Hierin ist Q1 ein Gütefaktor, der d.ie mittlere Mindestamplitude darstellt, die für eine bei der weiteren Verarbeitung zu betrachtende Reflexion akzeptabel ist, und. Q2 ist ein Gütefaktor, der die Mindestanzahl von PGT-Zeiteingängen darstellt, die für eine bei der weiteren Verarbeitung zu berücksichtigenden Reflexion akzeptabel ist. Sollte die mittlere Amplitude oder die Anzahl der PGT-Gruppen-Zeiteingänge nicht Q1 bzw. Q2 überschreiten, dann werden die errechneten Werte für die mittlere Zeit T, die mittlere scheinbare Geschwindigkeit Va und die mittlere Neigung a entsprechend Verfahrensschritt 22 nicht für die weitere Verarbeitung benutzt.First, the mean amplitude A is determined according to the following equation For example, the mean amplitude of the reflection Z1 of PGT groups 1 to 6 on line A is given by A + A + A (Z1) (A) = A11 + 12 + A1 + -A14 + A15 + A16 (11) The reflection for example Z1 on line A is considered acceptable for further processing if the following conditions are met W (Zi) (A)> Q1 (12) N> Q2 (13) Here, Q1 is a figure of merit, which is the minimum mean amplitude represents that is acceptable for a reflection to be considered in further processing, and. Q2 is a quality factor that represents the minimum number of PGT time inputs that is acceptable for a reflection to be taken into account in further processing. If the mean amplitude or the number of PGT group time inputs does not exceed Q1 or Q2, then the calculated values for the mean time T, the mean apparent speed Va and the mean inclination a according to method step 22 are not used for further processing .

Der nächste Schritt im Verfahren ist, wie bei 23 angegeben, d.ie Kontinuität der Reflexion zu bestimmen, d. h. ob dieselbe Reflexion auf allen drei Profillinien festgestellt worden ist. Es wird Reflexionsstelle beispielsweise auf d.er Linie A ausgewählt.As indicated at 23, the next step in the process is continuity to determine the reflection, d. H. whether the same reflection on all three profile lines has been established. It becomes a reflection point, for example on the line A selected.

Als zweites werden dann jede der Linien B und C daraufhin geprüft, ob dieselbe oder gemeinsame Reflexion zu finden ist. Dieser Schritt kann am besten anhand der Fig. 7A, 7B und 7C verstanden werden. Die mittlere Zeit für jede der Reflexionen Z1 bis Z6 für Linie A, wie diese nach Schritt 2? bestimmt ist, wird in der für Linie A in Fig. 7A bezeichneten Spalte gespeichert. Von der Reflexion Z1, die zu der mittleren Zeit TZ1(A) auf Linie A auftritt, kann während der Untersuchung der Linien B und C gefunden werden daß sie zu der mittleren Zeit Tz6(B) auf Linie B auftritt und zur mittleren Zeit TZ3(C) auf Linie C. D.h., die erste Reflexion Z1 auf Linie A, die sechste Reflexion Z6 auf Linie B und die dritte Reflexion Z3 auf Linie C treten sämtlich zur etwa selben mittleren Zeit auf. Diese beiden Zeiten Tz6(B) und TZ3(C) werden dann in der Zeitmatrix in den Spalten B und C in Fig. 7A gespeichert. Diese Uberprüfungsprozedur ermittelt jene Reflexionen auf Linien B und C, die in der Zeit einer Reflexion auf Linie A gemäß folgenden Beziehungen entsprechen: | T(A) - T(B) | # K1 und (14) |T(A) - T(C) | # K1 (15) hierin ist K1 eine Konstante.Second, each of the lines B and C are then checked to see whether the same or common reflection can be found. This step can best be understood with reference to Figures 7A, 7B and 7C. The mean time for each of the reflections Z1 through Z6 for line A, like this one after step 2? is determined is stored in the column designated for line A in Figure 7A. The reflection Z1 that occurs at the mean time TZ1 (A) on line A can be found during the examination of lines B and C that it occurs at the mean time Tz6 (B) on line B and at the mean time TZ3 ( C) on line CDh, the first reflection Z1 on line A, the sixth reflection Z6 on line B, and the third reflection Z3 on line C all occur at approximately the same mean time. These two times Tz6 (B) and TZ3 (C) are then stored in the time matrix in columns B and C in FIG. 7A. This verification procedure determines those reflections on lines B and C which correspond in time to a reflection on line A according to the following relationships: | T (A) - T (B) | # K1 and (14) | T (A) - T (C) | # K1 (15) here, K1 is a constant.

Beispielsweise kann es sein> daß die gesuchte gemeinsame Reflexion nicht genau zum selben Zeitpunkt auf jeder der drei Linien auftritt. Deshalb werden die Reflexionen auf Linien B und C, die zeitlich von der Reflexion auf Linie A sich nicht mehr als um eine Konstante, beispielsweise um 10 Millisekunden, voneinander tmterscheid.en, ausgewählt und in die Zeitmatrix nach Fig. 7A gegeben. Falls keine Zeitreflexion sowohl auf Linie B als auch auf Linie C vorhanden sind, die mit der gegebenen Konstante verträglich sind, dann wird keine Zeitreflexion in die Zeitmatrix für diesen Punkt gegeben. Dieses ist durch die Nullen in der zweiten Zeile von oben in jeder Spalte d.er Fig. 7A dargestellt.For example, it can be the common reflection you are looking for does not occur at exactly the same time on each of the three lines. Therefore be the reflections on lines B and C that are timed by the reflection on line A. no more than a constant, e.g. 10 milliseconds, from each other tmterscheid.en, selected and given in the time matrix of Fig. 7A. If not There are time reflections on both line B and line C, which coincide with the given constant compatible then there will be no time reflection given in the time matrix for this point. This is indicated by the zeros in the second Row from the top in each column of Fig. 7A.

Man sieht deshalb, daß jede Zeile, die gemeinsame mittlere Zeitreflexions-Eingänge hat, ein kontinuierliches Ereignis angibt, während jede Zeile, die keine Eingänge enthält, kein kontinuierliches Ereignis angibt.It can therefore be seen that each row has the common mean time reflection inputs has, indicating a continuous event, while each line has no inputs does not indicate a continuous event.

Nach Vervollständigung der Zeitmatrix nach Fig. 7A werden die Geschwindigkeitsmatrix der Fig. 7B und die Neigungsmatrix der Fig. 7C vervollständigt durch Eingeben der geeigneten mittleren Geschwindigkeiten und mittleren Neigungen entsprechend d.er mittleren Reflexionszeiteingänge für Linie A der Zeitmatrix der Fig. 7A. Beispielsweise sind die Geschwindigkeitseingänge für das kontinuierliche Ereignis zu den Zeiten TZ1(A)> TZ6(B) und TZ3() die Werte VZl(A)> Vz6 (B) und Vz3(c). In ähnlicher Weise sind die Neigungs-Eingänge #Z1(A), #Z6(B) und #Z3(C).After completing the time matrix of Figure 7A, the speed matrix of Fig. 7B and the slope matrix of Fig. 7C completed by entering the suitable mean speeds and mean inclinations according to d.er mean reflection time inputs for line A of the time matrix of Figure 7A. For example are the speed inputs for the continuous event at the times TZ1 (A)> TZ6 (B) and TZ3 () the values VZl (A)> Vz6 (B) and Vz3 (c). In a similar way Way are the incline inputs # Z1 (A), # Z6 (B) and # Z3 (C).

Gleichermaßen sind keine Eingänge in der zweiten Zeile von oben in Fig. 7B und 7C für nichtkontinuierliche Ereignisse vorhanden.Likewise, there are no inputs in the second row from the top in Figures 7B and 7C are for discontinuous events.

Der nächste Schritt im Verfahren ist, wie bei Schritt 24 angegeben, die Bestimmung des Mittelwertes der Reflexionszeiten für jene Reflexionen, die den drei Linien A, B und C gemeinsam sind. Dieses wird ausgerührt durch Ausmitteln der Reflexionszeiten auf jeder Zeile in Fig. 7A entsprechend folgender Beziehung: TCR = TA + TB + TC/3 (16) Beispielsweise ist der Mittelwert der Reflexionszeiten für die Reflexionen auf Linien B und C die der Reflexionszeit Tzl(A) auf Linie A gegeben durch TORi = Z1(A) Z6(B) + TZ3(C) (17) 3 An dieser Stelle ist es noch nicht möglich, eine echte Intervallgeschwindigkeit oder eine echte Neigung zu bestimmten, da die Werte in der Geschwindigkeits- und der Neigungsmatrix nach Fig. 7B und 7C Vektorkomponenten und nicht Vektorbeträge sind wie in Fig. 7A. Der nächste Schritt im Verfahren ist, wie bei 25 im Flußd.iagramm angegeben, die Verwendung der mittleren Neigung längs jeder Reflexion auf jeder Profillinie, wie dieses in der Neigungsmatrix 7C gespeichert ist, zusammen mit dem tatsächlichen Azimuth jeder Profillinie bezogen auf die echte magnetische Nordrichtung, um die echten Neigungs- und Deklinationsgrößen für jede Reflexion zu bestimmen. Diese Berechnungen der echten Deklination und Neigung kann, wie in SEISMIC PRO-SPECTING FOR OIL, von C. Hewitt Dix, 1952, Kapitel 9, Seiten 163-174, angegeben, für die drei möglichen Paarungen der drei Profillinien ausgeführt werden; d. h. für die Linien A und B, die Linien B und C und die Linien C und A. Aus diesen Rechnungen können drei Vektoren AB, BC und CA aufgetragen werden, wie dieses in Fig. 8 dargestellt ist.The next step in the procedure is, as indicated in step 24, the determination of the mean value of the reflection times for those reflections which the three lines A, B and C are common. This is carried out by averaging the Reflection times on each line in Fig. 7A correspond to the following relationship: TCR = TA + TB + TC / 3 (16) For example, is the mean value of the reflection times for the reflections on lines B and C that of the reflection time Tzl (A) on line A. given by TORi = Z1 (A) Z6 (B) + TZ3 (C) (17) 3 At this point it is not yet possible to determine a real interval speed or a real inclination, since the values in the speed and incline matrix of Figures 7B and 7C Vector components and not vector magnitudes are as in Figure 7A. The next step in the procedure, as indicated at 25 in the flow diagram, the use of the middle one Slope along each reflection on each profile line, like this one in the slope matrix 7C, along with the actual azimuth of each profile line to true magnetic north to get true pitch and declination quantities to be determined for each reflection. These calculations of real declination and inclination can, as in SEISMIC PRO-SPECTING FOR OIL, by C. Hewitt Dix, 1952, Chapter 9, pages 163-174, indicated, carried out for the three possible pairings of the three profile lines will; d. H. for lines A and B, lines B and C and lines C and A. From these calculations three vectors AB, BC and CA can be plotted, such as this is shown in FIG.

Um als gültiger Vektor betrachtet zu werden, muß der Vektor in ein vorgeschriebenes Verträglichkeitsgebiet fallen, wie dieses durch den schraffierten Bereich in Fig. 8 dargestellt ist und durch die nachstehenden Beziehungen gegeben ist: | #AB - #BC | # K2 (18) | #AB - #AC | # K2 (19) worin K2 eine Konstante, beispielsweise 15° für die Deklination ist, und | #AB - #BC | # K3 (20) | #AB - #AC | # K3 (21) worin K3 eine Konstante, beispielsweise 2°, für die Neigung ist.In order to be considered a valid vector, the vector must fall within a prescribed tolerance zone, as shown by the hatched area in Fig. 8 and given by the following relationships: | #AB - #BC | # K2 (18) | #AB - #AC | # K2 (19) where K2 is a constant, for example 15 ° for the declination, and | #AB - #BC | # K3 (20) | #AB - #AC | # K3 (21) where K3 is a constant, e.g. 2 °, for the slope.

Alle Deklinations- oder Neigungsvektoren, die nicht in das Verträglichkeitsgebiet fallen, werd.en verworfen. Die Deklinations-und Neigungsbestimmungen werden für gemeinsame Reflexionen gemittelt, ebenso wie dieses für die Zeitbestimmungen im Verfahrensschritt 24 erfolgte, um entsprechend den folgenden Beziehungen einen echten Deklinationswert QCR und einen echten Neigungswert #CR zu erhalten.All declination or inclination vectors that are not in the tolerance zone fall are discarded. The declination and inclination determinations are made for common reflections averaged, as well as this for the time determinations in Method step 24 was carried out to obtain a real one according to the following relationships To get a declination value QCR and a real slope value #CR.

#CR = #AB + #AC + #BC/3 (22) #CR = #AB + #AC + #BC/3 (23) Die vorstehenden Gleichungen für OCR und CR stellen die echte Neigung und Deklination der schrägverlaufenden Einbettung dar. #CR = #AB + #AC + # BC / 3 (22) #CR = #AB + #AC + # BC / 3 (23) the The equations above for OCR and CR represent the real slope and declination the inclined embedding.

Jedoch wird, wie erwähnt, eine Empfängerlinie kaum längs der Linie der echten Neigung verlaufen. Es ist deshalb notwendig, die mittleren scheinbaren Geschwindigkeiten Va, wie diese in Schritt 22 bestimmt worden sind, bezüglich der Differenz in der Beziehung zwischen der echten Neigungslinie und der Empfängerlinie auszukorrigieren. Die Korrektion, wenn im Verfahrensschritt 26 entsprechend der Gleichung (3) durchgeführt, liefert eine vMittelwertquadrat-Geschwindigkeit VRMS.However, as mentioned, a receiver line will hardly be along the line of true slope. It is therefore necessary to correct the mean apparent velocities Va as determined in step 22 for the difference in the relationship between the real incline line and the receiver line. The correction, if carried out in method step 26 in accordance with equation (3), provides a v mean square velocity VRMS.

Im Verfahrensschritt 26 wird die mittlere scheinbare Geschwindigkeit Va für jede Reflexion auf jeder Profillinie korrigiert bezüglich der echten Neigung und. der echten Deklination für jede Profillinie, um eine / Mittelwertquadrat-Geschwindigkeit entsprechend den folgenden Ausdrücken für jede der drei Profillinien zu erhalten. In method step 26, the mean apparent speed Va for each reflection on each profile line is corrected with respect to the real inclination and. the true declination for each profile line to one / Mean square speed according to the following expressions for each of the three profile lines.

Es sei daran erinnert, daß nur ein Neigungswinkel CR' d. <1. h. der CR Winkel zwischen der schrägverlaufenden Grenzfläche und der Erdoberfläche, vorhanden ist. Weiterhin sei daran erinnert, daß die Schnittlinie zwischen der Ebene der schrägverlaufend.en Grenzfläche und der Ebene der Erdoberfläche bezüglich einer bekannten oder gegebenen Richtung auf der Erdoberfläche beschrieben werden muß. In diesem speziellen Fall ist es die Richtung der Deklination gegenüber jeder der drei Profillinien A, B und C, die zur Korrektur der mittleren scheinbaren Geschwindigkeiten bezüglich der Neigung für jede der drei Profillinien erforderlich ist. Die Deklinationswinkel #A, #B und #C sind in Fig. 9 dargestellt.Recall that only one inclination angle CR 'd. <1. H. the CR angle between the sloping interface and the earth's surface, is available. It should also be remembered that the line of intersection between the plane the sloping interface and the plane of the earth's surface with respect to a known or given direction on the earth's surface must become. In this particular case it is the direction opposite to the declination each of the three profile lines A, B and C, which are used to correct the mean apparent Velocities related to the slope required for each of the three profile lines is. The declination angles #A, #B and #C are shown in FIG.

Entsprechend Gleichung (24) ist beispielsweise die gMittelwertquadrat- Geschwindigkeit für die Reflexion Z1 auf Linie A, korrigiert bezüglich der Neigung, Die Verträglichkeit der drei cttelwertquadrat-Geschwindigkeiten, die jeder gemeinsamen Reflexion (Fig. 7B) entsprechen, ist d.ann aus folgenden Beziehungen bestimmt: | VRMS(A) - VRMS(B) | # K4 (28) | VRMS(A) - VRMS(C) | # K4 (29) worin K4 eine Konstante ist.According to equation (24), for example, is g mean square Speed for the reflection Z1 on line A, corrected for the inclination, The compatibility of the three cttelwertquadrat-Velocities, which correspond to each common reflection (Fig. 7B) is then determined from the following relationships: | VRMS (A) - VRMS (B) | # K4 (28) | VRMS (A) - VRMS (C) | # K4 (29) where K4 is a constant.

Wenn der Absolutwert der Differenz zwischen VRMS(A) und VRMS(B) für eine gemeinsame Reflexion gleich oder kleiner als eine Konstante, beispielsweise 155 m/sec (500 Fuß/sec), dann werden die beiden Geschwindigkeiten als verträglich angesehen. In ähnlicher Weise werd.en VRMS(A) und VRMS(c) als verträglich betrachtet, wenn der Absolutwert ihrer Differenz für eine gemeinsame Reflexion kleiner als die Konstante von 155 m/sec ist.If the absolute value of the difference between VRMS (A) and VRMS (B) for a common reflection equal to or less than a constant, for example 155 m / sec (500 feet / sec) then the both speeds regarded as compatible. Similarly, VRMS (A) and VRMS (c) are found to be compatible considered when the absolute value of their difference for a common reflection is smaller than is the constant of 155 m / sec.

Beispielsweise würde für die gemeinsame Reflexion CR1 in der oberen Zeile in Fig. 7 B die Verträglichkeitsbestimmung die folgende sein: | VRMS(Z1)(A) - VRMS(Z6)(B) | # 155 m/sec, und (30) | VRMS(Z1)(A) - VRMS(Z3)(C) | # 155 m/sec (31) Alle verträglichen Geschwindigkeiten für jede gemeinsame Reflexion werden dann gemittelt, um die neigungskorrigierte Bestimmung der Mittelwertquadrat-Geschwind.igkeit zu erhalten. Beispielsweise wenn alle drei Profillinien-Geschwindigkeiten für eine gemeinsame Reflexion verträglich sind, dann gilt VRMS = VRMS(A) + VRMS(B) + VRMS(C)/3 (32) Beispielsweise ist die gMittelwertquadrat-Geschwindigkeit für die gemeinsame Reflexion CR1 in der oberen Zeile von Fig. 7B gegeben durch VRMS(CR1) = VRMS(Z1)(A) + VRMS(Z6)(B) + VRMS(Z3) (C)/3 Der nächste Schritt im Verfahren ist, wie bei 27 angegeben, die Erzeugung einer kontinuierlichen Kurve der errechneten g Mittelwertquadrat-Geschwindigkeiten. For example, for the common reflection CR1 in the top line in Fig. 7B, the compatibility determination would be as follows: | VRMS (Z1) (A) - VRMS (Z6) (B) | # 155 m / sec, and (30) | VRMS (Z1) (A) - VRMS (Z3) (C) | # 155 m / sec (31) All compatible velocities for each common reflection are then averaged to provide the tilt-corrected determination of the Mean square speed to obtain. For example, if all three profile line velocities are compatible for a common reflection, then VRMS = VRMS (A) + VRMS (B) + VRMS (C) / 3 (32), for example gMean square speed for the common reflection CR1 on the top row of Figure 7B given by VRMS (CR1) = VRMS (Z1) (A) + VRMS (Z6) (B) + VRMS (Z3) (C) / 3 The next step in the process is, as indicated at 27, the generation of a continuous curve of the calculated g mean square velocities.

Dieses wird bewerkstelligt durch Ziehen einer geraden Linie zwischen einem Zeitdiagramm der errechneten vMittelwertquadrat-Geschwindigkeit für jeden gemeinsamen Reflexionspunkt, wie dieses in Fig. 10 dargestellt ist.This is done by drawing a straight line between a timing diagram of the calculated v mean square velocity for each common reflection point as shown in FIG.

Der nächste Schritt im Prozeß, wie bei 28 angegeben, ist die Glättung der durch Schritt 27 erzeugten kontinuierlichen Kurve, um alle Unstetigkeiten in der Steigung zu beseitigen. Dieses wird bewerkstelligt durch Abtasten der kontinuierlichen Kurve zu periodischen Zeitintervallen und durch Anwenden eines Glättungsfilter-Operators auf d.ie abgetasteten Daten entsprechend der folgenden Beziehung: Cs C *h (25) Hierin bedeuten 0 die geglättete Kurve 0c die kontinuierliche Kurve, die zu glätten ist, und h der Glättungsfilter-Operator.The next step in the process, as indicated at 28, is smoothing of the continuous curve generated by step 27 to remove all discontinuities in to eliminate the slope. This is done by scanning the continuous Curve at periodic time intervals and by applying a smoothing filter operator on d. the sampled data according to the following relationship: Cs C * h (25) Herein 0 means the smoothed curve 0c the continuous curve to be smoothed, and h is the smoothing filter operator.

Die Glättungsoperation wird durch sukzessives Abwärtsverschieben des Glättungsfilter-Operators längs der kontinuierlichen Kurve entsprechend dem folgenden Integral ausgeführt: Hierin ist T das Schiebezeitintervall.The smoothing operation is carried out by successively shifting the smoothing filter operator downwards along the continuous curve according to the following integral: Here, T is the shift time interval.

Der Glättungsfilter-Operator, d.er bei dem vorliegenden Verfahren bevorzugt benutzt wird, ist der in Fig. 11 dargestellte Verschiebungsoperator. Das Verschiebungszeitintervall T ist das Zeitintervall über dem es gewünscht ist, die kontinuierliche Kurve zu glätten. Die Länge des Schiebeintervalls T kann variieren. Jedoch beträgt ein besonders geeignetes Intervall zum Glätten der kontinuierlichen Kurve 100 Millisekunden. Eine geeignete Abtastperiode für die Abtastpunkte längs der kontinuierlichen Kurve, auf die der Glättungsoperator angewandt werden soll, ist eine Periode von 4 Millisekunden; der Glättungsoperator wirkt daher auf 25 Abtastpunkte innerhalb seines Glättungszeitintervalls ein.The smoothing filter operator, i.e. the one in the present method is preferably used is the displacement operator shown in FIG. That Shift time interval T is the time interval over which it is desired that smooth continuous curve. The length of the shift interval T can vary. However, a particularly suitable interval for smoothing the continuous Curve 100 milliseconds. A suitable sampling period for the sampling points along the continuous curve to which the smoothing operator is to be applied, is a period of 4 milliseconds; the smoothing operator therefore acts on 25 sampling points within its smoothing time interval.

Der nächste Schritt im Verfahren, wie bei 29 angegeben, bestimmt die echte Intervallgeschwindigkeit Vi in einem Iterationsprozeß, der ein willkürlich schrägverlaurendes Snellius-Gesetz-Schichtmodell an die Daten anpaßt, die durch die in Schritt 28 erhaltene Glättungskurve der / Mittelwertquadrat-Geschwindig- keiten dargestellt sind. Eine solche Operation kann ausgeführt werden durch die Intervallgeschwindigkeits-Bestimmungsmethode, STRATV, d.ie im einzelnen in der US-PS 3,611,278 beschrieben ist.The next step in the method, as indicated at 29, determines the true interval velocity Vi in an iterative process which fits an arbitrarily sloping Snellius law layer model to the data represented by the smoothing curve obtained in step 28 / Mean square speed are shown. Such an operation can be carried out by the Interval Rate Determination Method, STRATV, that is described in detail in U.S. Patent 3,611,278.

Der nächste Verfahrensschritt ist, wie bei 30 angegeben, die Bestimmung der mittleren Geschwindigkeit V(t). Dieses wird bewerkstelligt durch Integrieren der in Schritt 29 bestimmten Intervallgeschwindigkeit Vi über der Zeit Der letzte Schritt im Prozeß ist, wie bei 31 angegeben, die Darstellung der Daten, wie dieses in Fig. 12 und 13 angegeben ist. Fig. 12 stellt ein Neigungs- und Deklinationsdiagramm dar, wobei d.ie Verschiebung der Pluszeichen die Größe der Neigung darstellt und wobei die Richtung der Neigung gegenüber Norden gegeben ist durch die den Pluszeichen beigegebenen Linien.The next method step is, as indicated at 30, the determination of the mean speed V (t). This is accomplished by integrating the interval velocity Vi determined in step 29 over time The final step in the process, as indicated at 31, is to display the data as indicated in FIGS. Fig. 12 shows a slope and declination diagram, where the displacement of the plus signs represents the magnitude of the slope and where the direction of the slope to north is given by the lines attached to the plus signs.

Drei Geschwindigkeiten sind in Fig. 13 dargestellt; Die ttelwertquadrat-Geschwindigkeit die die neigungskorrigierte Intervallgeschwindigkeit Vi und die mittlere Geschwindigkeit V(t).Three speeds are shown in Figure 13; the ttel value square speed the incline-corrected interval speed Vi and the mean speed V (t).

Aus dem Vorstehenden ist ersichtlich, daß das vorliegende Verfahren unter Verwendung verschiedener allgemein bekannter Rechenapparatetypen durchgeführt werden kann. Der Prozeß ist insbesondere geeignet zur Durchführung mit einem Allzweckrechner. Während das Verfahren durch zahlreiche Programme realisiert werden kann, zeigen Fig. 14 A - C ein Flußdiagramm eines Rechnerprogramms in FORTRAN-Sprache, zum Instruieren eines Allzweckrechners, um den Prozeß durchzuführen. Eine Rechenanlage, die speziell für die vorliegenden Zwecke geeignet ist, ist der Control Data Corporation Rechner Model 6600 und enthält die folgenden Eingangs/Ausgangs-Komponenten: Steuerrechner, 65K-Speicher 6602 Bildschirmgerät 0681 Datenkanal-Converter 405 Kartenleser 3447 Kartenleser-Steuergerät 501 Zeilendrucker 3256 Zeilendrucker-Steuergerät Ein entsprechender Plotter, der sich in Verbindung mit dieser Rechenanlage eignet, ist das Calcomp Plotter Model 763.From the foregoing it can be seen that the present process performed using various well known types of computing apparatus can be. The process is particularly suitable for being carried out on a general purpose computer. While the process can be realized through numerous programs, show 14 A-C is a flow chart of a computer program in FORTRAN language, for instruction a general purpose computer to perform the process. A computing system that specifically suitable for the present purposes is the Control Data Corporation calculator Model 6600 and contains the following input / output components: Tax calculator, 65K memory 6602 display unit 0681 data channel converter 405 card reader 3447 Card reader control unit 501 Line printer 3256 Line printer control unit A corresponding Plotter that is suitable in connection with this computer system is the Calcomp Plotter model 763.

Die meisten FORTRAN-Programme umfassen Schleifenbildung, die Wiederholung von Gruppen von Angaben mit sich ändernden Parametern. Die Steuerangabe, die in den Fig. 14A - 14C betroffen ist, ist die "DO"-Angabe. Für ein detailliertes Verständnis der Anwendung von FORTRAN-Angaben sei verwiesen auf "Introduction to FORTRAN' von S. C. Plumb, McGraw-Hill Book Company, New York, New York (1964).Most FORTRAN programs involve looping, repetition of groups of information with changing parameters. The tax information included in 14A-14C is the "DO" indication. For a detailed understanding the use of FORTRAN information is referred to "Introduction to FORTRAN" by S. C. Plumb, McGraw-Hill Book Company, New York, New York (1964).

Entsprechend Fig. 14A - 14C gibt die READ-Angabe 1002 die Eingangsdaten ein, wie diese bei Block 50 angegeben sind. Die Eingangsparameter sind wie folgt: AL(I) = seismische Profillinien-Verlauf DRIFT = zulässiger Reflexionszeit-Gleichstandfehler NSETS = Anzahl kreuzender Linien NC = Reflexionsgütefaktor DS = kürzeste zulässige Reflexionsstelle DL = längste zulässige Reflexionsstelle Bei 51, DO LOOP, DO 5: sequentielle Eingabe von PGT-Daten von 1 bis NSETS.14A-14C, READ indication 1002 provides the input data as indicated at block 50. The input parameters are as follows: AL (I) = seismic profile line course DRIFT = permissible reflection time tie error NSETS = number of crossing lines NC = reflection quality factor DS = shortest permissible Reflection point DL = longest permitted reflection point At 51, DO LOOP, DO 5: sequential entry of PGT data from 1 to NSETS.

Bei 52, READ-Angabe 1000 Eingabe der Daten über Zeit (T), Entfernung (D), mittlere Entfernung (H), Amplitude (A) und Geschwindigkeit (V) für die Anzahl Reflexionen von 1 bis I-te Linie.At 52, READ indication 1000 input of the data over time (T), distance (D), mean distance (H), amplitude (A) and speed (V) for the number Reflections from 1 to I-th line.

Bei 53, DO LOOP, DO 500: sortieren der Reflexionen entsprechend.At 53, DO LOOP, DO 500: sort the reflections accordingly.

dem Bewertungsfaktor N(J,I) und eliminieren jenerReflexionen von der-weiteren Verarbeitung, deren Mindestzahl von PGT-Zeiteingängen kleiner sind als der Gütefaktor NC.the weighting factor N (J, I) and eliminate those reflections from the other Processing whose minimum number of PGT time inputs is smaller than the quality factor NC.

Bei 54, DO LOOP, DO 230: sortieren der Reflexionen nach der mittleren Entfernung zum Kreuzpunkt der drei seismischen Linien und eliminieren jener Reflexionen von d.er weiteren Verarbeitung, die größer als DL aber kleiner als DS sind. Die Steuerung wird dann zurückübertragen auf Block 51 durch die 5 CONTINUE-Angabe im Block 55, und es wird eine andere PGT-Datengruppe zur Verarbeitung eingegeben. Nachdem alle die PGT-Gruppen durch die Blöcke 51, 52, 53 und 54 verarbeitet worden sind, wird die Steuerung auf Block 56 übertragen.At 54, DO LOOP, DO 230: sort the reflections according to the middle one Distance to the intersection of the three seismic lines and eliminate those reflections of the further processing, which are larger than DL but smaller than DS. the Control is then transferred back to block 51 by the 5 CONTINUE indication in the Block 55 and another PGT data group is entered for processing. After this all of the PGT groups have been processed by blocks 51, 52, 53 and 54, control is transferred to block 56.

Bei 56, DO LOOP, DO 1010: Speicherung der Zeit, Geschwindigkeit und Neigung für alle Reflexionen auf der ersten Linie in der Zeit-, Geschwindigkeits- bzw. Neigungsmatrix.At 56, DO LOOP, DO 1010: Storage of time, speed and Inclination for all reflections on the first line in time, speed or slope matrix.

Bei 57, DO LOOP, DO 1021: Eingabe von Zeit, Geschwindigkeit und Neigung für alle Reflexionen auf der zweiten und übrigen der NSETS der seismischen Linien. Dann DO LOOP, DO 1020: suchen unter allen Reflexionen nach jenen, die in der Zeit von T(J,O) um nicht mehr als die Konstante DIFT differieren und Eingabe dieser Reflexionen in die Zeit-, Geschwindigkeits-und Neigungsmatrizen als gemeinsame Reflexionen, um d.adurch ein kontinuierliches Ereignis anzuzeigen.At 57, DO LOOP, DO 1021: Enter time, speed and incline for all reflections on the second and remaining of the NSETS of the seismic lines. Then DO LOOP, DO 1020: look among all reflections for those who are in time differ from T (J, O) by no more than the constant DIFT and input these reflections into the time, speed and incline matrices as common reflections, to thereby indicate a continuous event.

Bei 58, DO LOOP, DO 1040: Eingabe aller gemeinsamer Reflexionen von der ersten Zeile der Zeitmatrix.At 58, DO LOOP, DO 1040: Enter all common reflections from the first row of the time matrix.

Bei 59, DO LOOP, DO 1041: berechnen der normierten Neigungskomponenten für jede Linie.At 59, DO LOOP, DO 1041: calculate the normalized slope components for each line.

Bei 60, DO LOOP, DO 1042: berechnen der Neigungsvektoren für Jede Linie.At 60, DO LOOP, DO 1042: calculate the slope vectors for each Line.

Bei 61, DO LOOP, DO 1043: vergleichen der Neigungsvektoren für jedes Profilpaar für jede Reflexion, abweisen jener Neigungsvektoren, die nicht innerhalb eines vorgeschriebenen Verträglichkeitsgebietes liegen.At 61, DO LOOP, DO 1043: compare the slope vectors for each Profile pair for each reflection, reject those slope vectors that are not within a prescribed tolerance area.

Bei 62, DO LOOP, DO 1044: Berechnung der Intervallgeschwindigkeiten unter Verwendung der Deklination- und Neigung-Grenzflächen mit Hilfe entweder OPTION I, wenn der KEY gleich oder grösser als 1 ist, oder mit Hilfe von OPTION I; wenn der KEY kleiner als 1 ist. Die Parameter für OPTION I sind. die folgenden: V = SQRT (1-SIN (a)²COS(ß)²) # VL, VL = gemessene Überlagerungsgeschwindigkeit für jede seismische Linie a = Größe der Neigung p = Winkel zwischen seismischer Linie und Neigungsvektor, ( ) = Durchschnitt über der Anzahl Überkreuzungslinien VI = Intervallgeschwindigkeit zwischen I- und. (I+1)-Reflexion, und a = Differenzoperator für I- und (I+1)-Werte Die Parameter für OPTION II sind wie folgt: VI = Intervallgeschwindigkeit Z = Tiefe zu jeder Grenzfläche Y = Grenzflächenneigung und # = Grenzflächenneigungs-Azimuth.At 62, DO LOOP, DO 1044: Calculation of the interval speeds using the declination and slope interfaces using either OPTION I, if the KEY is equal to or greater than 1, or with the help of OPTION I; if the KEY smaller than 1 is. The parameters for OPTION I are. the the following: V = SQRT (1-SIN (a) ²COS (ß) ²) # VL, VL = measured superimposition speed for each seismic line a = magnitude of the slope p = angle between seismic Line and slope vector, () = average over the number of crossing lines VI = interval speed between I- and. (I + 1) reflection, and a = difference operator for I and (I + 1) values The parameters for OPTION II are as follows: VI = interval speed Z = depth to each interface, Y = interface slope and # = interface slope azimuth.

STRATV ist die Intervallgeschwindigkeit smethod.e, die im einzelnen in der US-PS 3,611>278 beschrieben ist.STRATV is the interval speed smethod.e, which in detail is described in U.S. Patent 3,611> 278.

Bei 63 werden die Deklinations-, Neigungs- und Geschwindigkeitsdaten gegen die Zeit aufgetragen. Die Steuerung wird zurückgegeben auf Block 58 mit Hilfe der 1040 CONTINUE-Angabe des Blocks 64, und eine weitere Zeile gemeinsamer Reflexionen von der Zeitmatrix werden durch die DO 1040-Angabe eingegeben.At 63 the declination, incline and speed data plotted against time. Control is returned to block 58 with help the 1040 CONTINUE indication of block 64, and another line of common reflections from the time matrix are entered using the DO 1040 specification.

Claims

1. A method for producing a record of the earth's true interval velocity characteristics in the presence of an inclined embedding using an automatic system characterized by (a) obtaining seismic reflections from a plurality of at least three seismic profile lines which intersect at a common point on the surface, ( b) averaging the reflection times for those reflections common to all of the seismic profile lines; (c) correcting those reflections for inclination common to all profile lines to mean value {Average-

square speeds

to obtain, (d) generating a continuous curve of the \ / Mean square velocities

(e) smoothing the continuous curve of the f mean square

Velocities and (f) fitting an arbitrarily chosen inclined Snellius law layer model to the smoothed curve of the f / Mean square velocities,

in order to get real interval speeds. 2. The method according to claim 1, characterized in that the y mean square velocities

are determined by: (a) averaging the velocities for each reflection along each profile line to obtain an average apparent velocity) (b) averaging the slopes for each reflection along each profile line to obtain an average inclination, (c) calculating the true declination and real slope for those reflections common to all profile lines from the mean slope for each reflection along each profile line and from the azimuth of each profile line, (d) correcting each of the mean apparent velocities for real declination and real slope for Receiving a / Mean square

Velocity for each reflection along each profile line and (e) averaging the 2 mean square velocities

for those reflections that are common to all profile lines.

3. The method according to claim 2, characterized in that the real Interval velocity, the real declination, and the real incline over it Time scale to be plotted.

4. The method according to claim 1, characterized in that the generation of a continuous curve of the Y 'mean square' -

Velocities is done by (a) applying the Z mean square velocities

about the average reflection times for those reflections that are common to all profile lines and (b) connecting the / Mean square speed

score with a straight line below each other.

5. Process according to claim 1, characterized in that the smoothing of the continuous curve of the JMZ / Mean square

Velocities are accomplished by (a) sampling at periodic time intervals the data represented by the continuous curve and (b) applying a smoothing filter operator to the sampled data during a period of time for which the continuous curve is to be smoothed.

6. The method according to claim 1, characterized in that the real Interval speed for those reflections are integrated over time, the are common to all profile lines in order to maintain the average speed.

Blank page