DE2325071A1

DE2325071A1 - PROCEDURE FOR DETERMINING INTERFERENCE CORRECTION SIZES FOR SEISMIC DATA AND ARRANGEMENT FOR PERFORMING THE PROCEDURE

Info

Publication number: DE2325071A1
Application number: DE19732325071
Authority: DE
Inventors: Charles D Dirks
Original assignee: Texas Instruments Inc
Current assignee: Texas Instruments Inc
Priority date: 1972-05-17
Filing date: 1973-05-17
Publication date: 1973-12-06
Also published as: GB1422657A; FR2184986A1; NL7305717A; JPS4966502A; AU5446673A

Description

Unser Zeichen; T 1384 Our sign; M 1384

TEXAS INSTRUMENTS INCORPORATED
13500 North Central Expressway
Dallas, Texas, V.St.A.TEXAS INSTRUMENTS INCORPORATED
13500 North Central Expressway
Dallas, Texas, V.St.A.

Verfahren zur Bestimmung von Störkorrekturgrößen für seismische Daten und Anordnung zur DurchführungProcedure for the determination of disturbance correction quantities for seismic data and arrangement for implementation

des Verfahrensof the procedure

Die Erfindung bezieht sich auf seismische Exploration und insbesondere auf die Entfernung von dicht bei der Oberfläche erzeugten Zeit- und Geschwindigkeitsverschiebungen bei seismischen Daten«The invention relates to seismic exploration and, more particularly, to removal from close to the surface generated time and speed shifts in seismic data «

Die seismische Exploration unter Ausnützung von Punkten gemeinsamer Tiefe (CDP-Exploration) wird derzeit bei der Suche nach Lagerstätten von Erdöl und anderen Mineralien häufig angewendet. Bei der CDP-Exploration werden Schuß«- punkte an aufeinanderfolgend voneinander entfernten Orten hintereinander zur Explosion gebracht, während die Verteilung von Geophonen entsprechend bewegt wird. Es werden seismische Spuren aufgezdchnet, die durch eine Tiefenschichtredündanz gekennzeichnet sind_s so daß sich besondere Vorteile bei der Verarbeitung zur Rauschminderung ergeben. Ein typischer seismischer CDP-Explorationsvorgang ist in der US-PS 2 732 906 beschrieben.; ■Seismic exploration using points of common depth (CDP exploration) is currently widely used in the search for deposits of crude oil and other minerals. In CDP exploration, bullet points are detonated one after the other at locations that are spaced apart from one another, while the distribution of geophones is moved accordingly. Seismic traces are counted which are characterized by a depth layer redundancy _s so that special advantages arise when processing for noise reduction. A typical CDP seismic exploration operation is described in U.S. Patent 2,732,906 .; ■

30 9849/0 469 "30 9849/0 469 "

Damit bei der CDP-Exploration·aussagekräftige Daten erhalten werden, müssen an den seismischen Aufzeichnungen Störkorrekturen durchgeführt werden, damit Fehler auf ' ein Minimum verringert werden,· die infolge von oberflächennahen Zuständen wie der Verwitterung u.dgl. auftreten. Zur Erzeugung von Störkorrekturen für CDP-Aufzeichnungen sind daher verschiedeneverfahren entwickelt worden, von denen die meisten Kreuzkorrelationsverfahren erfordern.So that during the CDP exploration · meaningful data are obtained Corrections must be made to the seismic records to ensure that errors occur. a minimum that occur as a result of conditions close to the surface such as weathering and the like. For generating noise corrections for CDP recordings therefore, various methods have been developed, most of which are cross-correlation methods require.

Ein weiteres bekanntes Verfahren zur Beseitigung von Störgrößen ist in der US-PS 3 539 982 beschrieben. Bei diesem Verfahren wird die Kreuzkorrelation zur Feststellung von Zeitverschiebungen jeder seismischen Spur angewendet, und die Zeitverschiebungen v/erden dann entsprechend gemeinsamen Schußpunkt- und Empfängerpositionen gruppiert und gemittelt. Dieses Verfahren hat sich zwar in der Praxis bewährt, doch erfordern Kreuzkorrelationsverfahren für ihre Wirksamkeit ganz allgemein eine relativ große Anzahl von Datenpunkten. Es ist jedoch häufig schwierig, in einer Zone mit schlechtem Rauschabstand eine genügende Anzahl von Datenpunkten zu finden, und die einzelnen Zeitverschiebungen werden oft zumindest teilweise von statistischem und kohärentem Rauschen beeinflußt.Another known method for eliminating disturbance variables is described in U.S. Patent 3,539,982. In this method, the cross-correlation is used to determine Time shifts of each seismic trace are applied, and the time shifts are then correspondingly common Shot point and receiver positions grouped and averaged. Although this method has proven itself in practice, it requires cross-correlation methods to be effective generally a relatively large number of data points. However, it is often difficult to be in a zone with bad S / N ratio to find a sufficient number of data points, and the individual time shifts are often influenced at least in part by statistical and coherent noise.

Ferner können mit dem in der US-PS 3 539 982 beschriebenen Verfahren nicht allgemein vorhandene Kenntnisse zuvor ausgewählter Störgrößen bei der Unterstützung der Voraussage der nächsten zu erzeugenden Störkorrekturgröße verwendet werden. Daß dies nicht möglich ist, ergibt sich daraus, daß aus den Korrelationen einzelne Zeitverschiebungen herausgegriffen werden, wobei die Zeitverschiebungen eineFurthermore, prior knowledge not generally available with the method described in US Pat. No. 3,539,982 selected disturbance variables are used to support the prediction of the next disturbance correction variable to be generated will. That this is not possible results from the fact that individual time shifts arise from the correlations be singled out, with the time shifts a

30 98 4.9/0 489-30 98 4.9 / 0 489-

Kombination von Schußpunkt- und Empfängerstörgrößen sind, so daß sie nicht allgemein irgendwelchen eingestellten Mustern folgendUe Unfähigkeit, bereits vorher vorhandene Kenntnisse anzuwenden_s macht es oft äußerst schwierig, Störgrößen festzustellen, die größer als die halbe Zeitperiode der seismischen Daten sind, wodurch die Peststellung niederfrequenter Störgrößen eingeschränkt wird.Combination of Schußpunkt- and Empfängerstörgrößen are so that they apply not generally any set patterns folgendUe inability previously existing knowledge _s often makes it extremely difficult to detect disturbances that are greater than half the period of seismic data whereby the Peststellung low frequency disturbances is restricted.

Ein weiterer Nachteil des in der US-PS 3 539 982 beschriebenen Verfahrens besteht darin₂ daß es nicht möglich ist, die Einflüsse der Störgrößen auf die Scheingeschwindigkeit der Daten zu beseitigen«, Da dieses bekannte Verfahren die Eingabe geschwindigkeitskorrigierter Daten erfordert, steht die Geschwindigkeitsinformation daher für eine Abfrage während der Störgrößenkorrekturen nicht zur Verfügung« Die oben erwähnte Unfähigkeit des bekannten-Verfahrens, niederfrequente Störungen festzusteilen, behindert ebenfalls die Geschwindigkeitskorrektur durch "das Verfahren.A further disadvantage of the method described in U.S. Patent No. 3,539,982 is ₂ that it is not possible, the effects of the disturbances on the apparent speed of data to eliminate "As this known method requires the input speed corrected data is available, the speed information, therefore, not available for an interrogation during the disturbance variable corrections. The above-mentioned inability of the known method to fix low-frequency disturbances also hinders the speed correction by the method.

Ein weiteres bekanntes Verfahren zur Störungsbeseitigung ist in einem Aufsatz mit dem Titel "Statistical Automatic Statics Analysis" von D.A.Di^vsher und anderen beschrieben, der in der Zeitschrift Geophysics, XXXVg Nr. 4, Seiten 574 bis 585 vom August 1970 veröffentlicht -worden ist. Das in diesem Aufsatz beschriebene Verfahren gleicht in gewisser Weise dem in der oben genannten US-PS 3 539 982 beschriebenen Verfahren, da Kreuzkorrelationen angewendet ,werden, um Zeitr verschiebungen zu erhalten, die dann zur Bestimmung von Störkorrekturgrößen verarbeitet werden. Das in dem genannten Aufsatz beschriebene Verfahren weist somit die gleichen Nachteile auf, wie oben bereits im Zusammenhang mit dem in der US-PS 3 539 982 beschriebenen Verfahren erwähnt worden sind. Außerdem wird das in dem Aufsatz beschriebene VerfahrenAnother known method of troubleshooting is described in an article entitled "Statistical Automatic Statics Analysis" by DADi ^v sher and others, published in August 1970 in the journal Geophysics, XXXVg No. 4, pages 574 to 585. The method described in this paper is somewhat similar to the method described in the aforementioned US Pat. No. 3,539,982 in that cross-correlations are used to obtain time shifts which are then processed to determine interference correction quantities. The process described in the cited article thus has the same disadvantages as have already been mentioned above in connection with the process described in US Pat. No. 3,539,982. In addition, the procedure described in the article

309849/0489309849/0489

■ 2325Q71■ 2325Q71

noch mehr als das in der genannten USA-Patentschrift genannte Verfahren von den Schwierigkeiten und schlechtem Rauschabstand beeinflußt, da keine gestapelten Spuren angewendet werden, so daß sich die daraus gewonnene Verbesserung des Rauschabstandes nicht ergibt. Ferner erfordern die beiden bekannten Verfahren allgemein vorherige Kenntnisse der Geschwindigkeit der Daten, und sie sind für Fehler auf diesem Gebiet empfindlich.even more than that mentioned in the aforementioned US patent Method affected by the difficulty and poor signal-to-noise ratio since no stacked tracks are used so that the improvement in the signal-to-noise ratio obtained therefrom does not result. Also require Both of the known methods generally have prior knowledge of the speed of the data, and they are sensitive to errors in this area.

Mit Hilfe der Erfindung wird ein Verfahren zur Bestimmung von Störkorrekturgrößen für seismische Daten mit Tiefenschichtredundanz geschaffen. Bei dem Verfahren werden aus einer digitalen Segmentdatei Hauptsegmente ausgewählt, die vorbestimmte Parameter erfüllen. In Abhängigkeit von den Hauptsegmenten werden aus den Sammelaufzeichnungen der seismischen Daten mehrere Wellenzüge (wavelets) herausgegriffen. Die Wellenzüge werden dann in ¥ellenzugkartenanordnungen so_gruppiert, daß die Einflüsse der Zeit, der Geschwindigkeit und der Geologie auf die Wellenzüge beseitigtWith the aid of the invention, a method for determining interference correction variables for seismic data with depth layer redundancy is provided created. In the process, main segments are selected from a digital segment file, meet predetermined parameters. Depending on the main segments, the collective recordings become several wave trains (wavelets) selected from the seismic data. The wave trains are then arranged in ellipsoidal train cards grouped in such a way that the influences of time, speed and geology on the wave trains are eliminated

1 i 1 i

werden. Anschließend werden die Wellenzüge entsprechend der gemeinsamen"Schußpunktpositionen und der gemeinsamen Empfängerpositionen summiert. Es werden dann die Zeitintervalle bestimmt,-an denen ausgewählte Spitzen der summierten Wellenzüge auftreten. Im Anschluß daran, werden Schußpunkt- und Empfänger-Störkorrektursignale in Abhängigkeit von den festgestellten Zeitintervallen erzeugt.will. Then the wave trains are corresponding to the common "shot point positions" and the common receiver positions summed up. The time intervals at which selected peaks of the summed up are then determined Wave trains occur. Following this, the shot point and receiver clutter correction signals are dependent on the determined time intervals generated.

Ferner wird mit Hilfe der Erfindung ein Verfahren geschaffen, mit dem aus Sammelaufzeichnungen seismischer Daten für Punkte gemeinsamer Tiefe und aus einer seismischen digitalen Segmentdatei mit Schätzwerten für die Laufzeit und die Geschwindigkeit der seismischen Daten Störkorrekturgrößen bestimmt werden können. Aus der Segmentdatei werden mehrereFurthermore, with the help of the invention, a method is created with that from aggregate recordings of seismic data for points of common depth and from a seismic digital segment file with estimated values for the transit time and the speed of the seismic data can be determined. The segment file becomes several

309849/0469309849/0469

.5.5

Hauptsegmente ausgewählt, die vorbestimmte Parameter erfüllen. Aus den Sammelauf ze iclmungen werden dann mehrere Wellenzüge ausgewählt;, die allgemein, um die Laufzeiten und die Moveout-¥erte der Haupt segment e zentriert sind« (Unter "Moveout" wird die Zunahme der Laufzeit bedingt durch die Meßgeometrie verstanden) _β Die Wellenzüge werden in ¥ellenzugkartenanordnungen gruppiert, in denen die Wellenzüge bezüglich ihrer geschützten Laufzeiten längs. Koordinaten ausgerichtet werden₉ die den Punkt gemeinsamer Tiefe, gemeinsame- Schußpunktpositionen_s die gemeinsamen Empfängerpositionen und die gemeinsame Versetzung'angeben. Die ¥ellenzüge mit gemeinsamen Schußpunktpositionen werden dann summiert. Anschließend werden die Wellenzüge mit gemeinsamen Bmpf ängerpositionen summierte Dann werden die Zeitintervalle zwischen ausgewählten Spitzen der summierten Wellenzüge und der Mitte der summierten Wellenzüge bestimmt«, In Abhängigkeit von diesen Zeitintervallen werden dann Schußpunkt=· und Empfänger-Störkorrektursignale erzeugt.Main segments selected that meet predetermined parameters. Several wave trains are then selected from the collective records, which are generally centered around the transit times and the moveout values of the main segments ("Moveout" means the increase in the transit time due to the measurement geometry) _β The wave trains are grouped in ¥ ellenzugkarte arrangements in which the wave trains are longitudinal with respect to their protected transit times. Coordinates are aligned with the point ₉ of common depth point locations gemeinsame- weft _s the common receiver positions and the common Versetzung'angeben. The ¥ ellipses with common shot point positions are then summed up. Then the wave trains with common receiver positions are summed up. Then the time intervals between selected peaks of the summed wave trains and the middle of the summed wave trains are determined.

Gemäß einer¹ Weiterbildung der Erfindung wird nach der anfänglichen Bestimmung der SchuBpunkt-und Empfänger-Störkorrektursignale jeder der Wellenzüge in den Wellenzugkartenanordnungen gemäß den festgestellten, jedem Wellenzug zugeordneten Schußpunkt= und Empfängerstör·» großen verschoben«, Die Summierungs- und Feststellungssshritte des erfindungsgemäßen Verfahrens werden dann zur Erzeugung verbesserter Störkorrektursignale wiederholt durchlauf en..'Außer dem werden nach der Verschiebung der Wellenzüge in den Wellenzugkartenanordnungen und nach dem erwähnten wiederholten Durchlauf die Daten für den Punkt gemeinsamer Tiefe entsprechend den Schußpunkt- und Empfänger- Störgrößenkorrektursignalen verschoben, und die Schritte des erfindungsgemäßen Verfahrens werdenAccording to a ¹ further development of the invention, after the initial determination of the thrust point and receiver Störkorrektursignale each of the wave trains shifted in accordance with the identified, each wave train associated shot point = and Empfängerstör · "great in the Wellenzugkartenanordnungen" The sum and Feststellungssshritte be of the method according to the invention then repeatedly run through to generate improved interference correction signals .. 'In addition, after the shifting of the wave trains in the wave train card arrangements and after the mentioned repeated run, the data for the point of common depth are shifted according to the shot point and receiver interference variable correction signals, and the steps of the invention Procedure

309849/0469309849/0469

zur Erzeugung verbesserter Störkorrektursignale in einer ausreichenden Anzahl erneut durchlaufen. ■run through a sufficient number of times again to generate improved interference correction signals. ■

In einer weiteren Ausgestaltung der Erfindung' werden Spitzen der summierten Wellenzüge erfaßt, in dem Wellenzugsspitzen eliminiert werden, die unter einem vorgewählten Amplitudenverhältnis zur zuletzt festgestellten Spitze liegen. Die absolute Zeit jeder verbleibenden Spitze bezüglich der • Zeit des zuletzt festgestellten Wellenzugs wird dann bestimmt. Es wird die Spitze festgestellt, die der vorbestimmten absoluten Zeit am nächsten liegt, wenn nicht eine andere Spitze des Wellenzugs eine in einem vorbestimmten Maß größere Amplitude als die ausgewählte Spitze auf weist.Das Zeitintervall zwischen der ausgewählten Spitze und der Mitte des Wellenzugs wird dann zur Erzielung des gewünschten Störgrößenkorrektur-Schätzwerts bestimmt.In a further embodiment of the invention, tips are used of the summed wave trains are detected in which wave train peaks are eliminated which are below a preselected amplitude ratio to the last point determined. The absolute time of each remaining peak in relation to the • The time of the last determined wave train then becomes certainly. The peak is found which is closest to the predetermined absolute time, if not another peak of the wave train has an amplitude greater than the selected peak by a predetermined amount The time interval between the selected The peak and center of the wave train is then determined to obtain the desired disturbance correction estimate.

Ein Ausführungsbeispiel der Erfindung ist' in der Zeichnung dargestelli;. Darin zeigen:An embodiment of the invention is' in the drawing dargestelli ;. Show in it:

Fig.1 eine schematische Darstellung des. Auftretens von Störzeitverschiebungen bei der seismischen Exploration,1 shows a schematic representation of the occurrence of Disturbance time shifts in seismic exploration,

Fig.2 ein Flußdiagramm zur Ausführung des erfindungsgemäßen Störkorrekturverfahrens auf einem digitalen Computer,2 shows a flow chart for carrying out the inventive method Interference correction procedure on a digital computer,

Fig.3 ein Flußdiagramm des in Fig.2 gezeigten Segmentauswahlprozesses, Figure 3 is a flow chart of the segment selection process shown in Figure 2;

Fig.4 ein Flußdiagramm zur Veranschaulichung des in Fig. gezeigten Aufbaus der Wellenzugkartenanordnung,4 shows a flow chart to illustrate the process shown in FIG. the structure of the wave train card arrangement shown,

Fig.5 ein Diagramm der etwa auf einer hyperbolischen Kurve5 shows a diagram of the approximately on a hyperbolic curve

. liegenden Reflexionszeiten von von einem Tiefensehichtpunkt ausgehenden Peflexionswellenzügen auf einer. reflection times from a depth viewpoint outgoing reflection wave trains on a

309849/0469309849/0469

'Sammelaufzeichnung für einen Punkt gemeinsamer Tiefe,'' Collective recording for a point of common depth,

Fig.6 den Aufbau einer Wellenzugkarte für einen Punkt gemeinsamer Tiefe für den in Fig„5 dargestellten Wellenzug,Fig. 6 shows the structure of a wave train map for a point common depth for the one shown in FIG Wave train,

Fig.7 die hyperbolischen Kurven der Reflexionswellenzüge auf einer Sammelaufzeichnung für vier Punkte gemeinsamer Tiefe,Fig. 7 the hyperbolic curves of the reflection wave trains on a collective record for four points of common depth,

Fig.8 eine Wellenzugkarte für die in Fig.7 dargestellen Punkte gemeinsamer Tiefe,Fig.8 shows a wave train map for the illustrated in Fig.7 Points of common depth,

Fig«,9 ein Flußdi'agramm des in Fig.2 dargestellten Verfahrens "Wellenzugkarten-Smash",FIG. 9 shows a flow diagram of the method shown in FIG "Wave Train Card Smash",

Fig.10 ein Flußdiagramm des in Fig_e2 gezeigten Störgrößen erfassungsschritts und10 shows a flowchart of the interference detecting step shown in Figure 2 and _e

Fig. 11 bis 13 Diagramme verschiedener Wellenzüge zur Veranschaulichung des Spitzenerfassungsverfahrens von Fig.10. ' ■ . .,11 to 13 are diagrams of various wave trains for illustration of the peak detection method of Fig. 10. '■. .,

In Fig.1 ist in einem Diagramm eine seismische Operation für eine mit Punkten gemeinsamer Tiefe arbeitende Exploration (CDP-Sxploration) dargestellt. Eine Quelle 10 für seismischeIn Fig.1 is a diagram of a seismic operation for an exploration working with points of common depth (CDP exploration). A source 10 for seismic

Impulsej beispielsweise ein Dynamitschußpunkt, erzeugt seismische Wellen, die eine unregelmäßige Schicht mit niedriger Ausbreitungsgeschwindigkeit und dann Schichten mit höherer Ausbreitungsgeschwindigkeit durchlaufen, bis sie von einem reflektierenden Horizont 12 reflektiert werden. Die seismischen Reflexionen werden dann von > Geophonen 14a bis I4x empfangen»Impulse j, for example, a dynamite shot, creates seismic waves that traverse an irregular, slow-moving layer and then higher-speed layers until they are reflected off a reflective horizon 12. The seismic reflections are then received by > geophones 14a to I4x »

98 4 9/04698 4 9/046

- 3 - ■- 3 - ■

2325Q712325Q71

Die seismische CDP-Exploration kann zwar mit einer Anzahl' von unterschiedlichen geometrischen Oberflächenanordnungen der Schußpunkte und der Geophone ausgeführt werden, doch wird typischerweise mit der üblichen Geometrie unter Verwendung von 24 Geophonen (off-end shot-spread geometry) gearbeitet, die für eine Sechsfach-Tiefenpunkt-Tiefenschicht» ■erfassung besonders brauchbar ist. ¥ie in der US-PS 2 732 906 beschrieben ist, werden bei CDP-Operationen an der Oberfläche Schußpunkte an aufeinanderfolgend voneinander entfernten Orten angebracht, während gleichzeitig die Geophonverteilung bewegt wird. Diese CDP-Operation erlaubt eine kontinuierliche Mehrfacherfassung einer Tiefenschicht zur Erzeugung von' Daten mit einer großen Tiefenschichtredundanz, die mit Hilfe des hier angegebenen Verfahrens vorteilhaft verarbeitet werden können.CDP seismic exploration can be carried out with a number of ' are carried out by different geometrical surface arrangements of the shot points and the geophones, but is typically done with the usual geometry using 24 geophones (off-end shot-spread geometry) worked for a six-fold depth point depth layer » ■ acquisition is particularly useful. ¥ ie in the US PS 2 732 906, in CDP operations on the surface, bullet points are made at successive points at a distance from each other while the geophone distribution is moved at the same time. This CDP operation allows a continuous multiple acquisition of a depth layer to generate 'data with a large Depth layer redundancy, which with the help of the given here Method can be processed advantageously.

Störkorrekturen sollen unregelmäßige Zeitverzögerungen beseitigen, die von dem Material zwischen der Oberfläche und einer beliebigen Bezugsebene verursacht v/erden. Die Schußpunkt-Störkorrektur umfaßt daher die Laufzeit der seismischen Welle vom Punkt a bis zum Punkt c, während die Empfänger-Störkorrektur die Laufzeit vom Punkt e bis zum Punkt g umfaßt. Die mit Hilfe des hier zu beschreibenden Verfahrens auszuführenden Störkorrekturen werden als unabhängig vom Austrittswinkel der Strahlenbahn angesehen. Diese Annahme ist zwar nicht genau richtig, doch ist sie im allgemeinen eine befriedigende Näherung , wobei wesentliche Ausnahmen allgemein nur dann auftreten, wenn die Basis der Verwitterungsschicht sehr unregelmäßig oder sehr tief ist.Disturbance corrections are intended to have irregular time delays eliminate grounds caused by the material between the surface and any reference plane. the Shot point interference correction therefore includes the running time of the seismic wave from point a to point c while the receiver interference correction the transit time from point e includes up to point g. The interference corrections to be carried out with the aid of the method to be described here are considered to be independent of the exit angle of the beam path. While this assumption is not accurate correct, but it is generally a satisfactory approximation, with major exceptions generally only then occur when the base of the weathering layer is very irregular or very deep.

Die Ausgangssignale der Geophone I4a bis 14x v/erden durch ein herkömmliches Aufzeichnungsgerät 16 aufgezeichnet, damit ein. Seismogramm mit Daten erzeugt wird, das insbesondere mit dem reflektierenden Horizont 12 und auch mitThe output signals of the Geophone I4a to 14x v / ground recorded by a conventional recorder 16 so that a. Seismogram is generated with data that in particular with the reflective horizon 12 and also with

■309849/0489■ 309849/0489

möglicherweise vorhandenen reflektierenden Horizonten in Beziehung steht. Vor der Anwendung der hier zu beschreibenden Verfahrensschritte werden dann auf die aufgezeichneten seismischen Spuren übliche Verarbeitungsschritte verwendet, bei denen in üblicher Weise a priori Bezugsstörgrößen oder andere Verarbeitungsverfahren wie eine Rücktransformation nach einer Faltung oder eine Wiedergewinnung der wahren Amplitudenwerte angewendet werden. Diese a priori Bezugsstörgrößen werden zur Korrektur der Daten hinsichtlich einer Bezugsebne 19 eingesetzt, indem Oberflächenerhebungen und geeignete oberflächennahe Geschwindigkeitsschätzwerte und Aufrechterhaltungszeiten bei Verwendung von Dynamitschußpunkten angewendet v/erden. Diese Bezugsstörgrößen werden üblicherweise bei der seismischen Explosion angewendet; wegen unregelmäßiger Änderungen der Dicke und der Geschwindigkeit · der Schicht mit niedriger Ausbreitungsgeschwindigkeit steLlen sie jedoch keine Lösung des Störproblems dar. Nach oben oder nach unten abweichende Fehler oder Reststörkorrekturen bezüglich der Bezugsebene bleiben icmer noch vorhanden, die mit Hilfe des hier zu beschreibenden Verfahrens beseitigt werden müssen, ehe aus dem Seismogramm aussagekräftige Daten erhalten werden können.possibly existing reflective horizons is related. Before applying the one to be described here Process steps are then used on the recorded seismic traces in the usual processing steps, in which in the usual way a priori Reference disturbance variables or other processing methods such as an inverse transformation after a convolution or a Recovery of the true amplitude values can be applied. These a priori reference disturbance variables are used for Correction of the data with regard to a reference plane 19 used by adding surface elevations and suitable near-surface speed estimates and maintenance times applied when using dynamite bullet points. These reference disturbance variables are usually applied to seismic explosion; due to irregular changes in thickness and speed however, they do not represent a solution to the disturbance problem in the layer with low propagation speed or downward deviating errors or residual interference corrections with regard to the reference plane, icmer still exist, which is eliminated with the help of the procedure to be described here must be before meaningful data can be obtained from the seismogram.

"Bei der Durchführung des hier zu beschreibenden Verfahrens wird das mit Hilfe des Aufzeichnungsgeräts Ϊ6 aufgezeichnete Seismogramm digitalisiert unizur Erzeugung einer Sammelaufzeichnung für Punkte gemeinsamer Tiefe (CDP-Samiiie !aufzeichnung) aller seismischen Spuren verarbeitet, die infolge der Aufzeichnungsgeometrie Reflexionen von dem-gleichen TLefenschichtreflexionspunkt enthalten. Außerdem wird das Seismogramm so verarbeitet, daß eine digitale Seginentdatei entsteht, in der Schätzwerte für die Laufzeit und die Geschwindigkeit der seismischen CDP-Daten enthält. Ein geeignetes- Verfahren zur Erzeugung der digitalen Segmentdatei"When performing the procedure to be described here is recorded using the recorder Ϊ6 Seismogram is digitized in order to generate a collective record for points of common depth (CDP-Samiiie! record) of all seismic traces processed, which are reflections of the same as a result of the recording geometry Tefe layer reflection point included. In addition, the seismogram is processed into a digital segment file which contains estimates of the transit time and velocity of the CDP seismic data. A suitable method for generating the digital segment file

30 9 84 9AH 68 " ■ ."30 9 84 9AH 68 "■."

ist das sogenannte "600 Processing Package", das derzeit von Geophysical Service, Inc.,einer Tochtergesellschaft der Firma Texas Instruments, Inc. vertreten wird. Das "600 Processing Package" ist in der Zeitschrift Geophysics, Dezember 1971, XXXVI, Seiten 1043 bis 1059 beschrieben. Die seismischen CDP-Daten können auch nach dem Verfahren "700 Processing Package" verarbeitet werden, das in der der USA-Patentanmeldung Serial Number 252 892 vom 12.Mai 1972 entsprechenden deutschen Patentanmeldung beschrieben ist.is the so-called "600 Processing Package" currently available from Geophysical Service, Inc., a subsidiary Texas Instruments, Inc. is represented. The "600 Processing Package" is in the journal Geophysics, December 1971, XXXVI, pages 1043-1059. The CDP seismic data can also be obtained according to the procedure "700 Processing Package" described in U.S. Patent Application Serial Number 252 892 dated May 12, 1972 corresponding German patent application is described.

Bei dem Verfahren "600 Processing Package" und "700 ProcessingPackage" werden CDP-Sammelaufzeichnungen so bearbeitet, daß digitale Segmentdateien erzeugt werden, die bei dem hier zu.beschreibenden Verfahren zur Erzeugung von Schußpunkt- und Empfänger-Störkorrekturen verwendet werden können.. Zur Durchführung des hier zu. beschreibenden Verfahrens kann zwar irgendein richtigprogrammierter digitaler Computer verwendet werden, doch wird vorzugsweise ein solcher digitaler Computer verwendet, wie er unter der Bezeichnung TIAC 870A von der Firma Texas Instruments Inc. Dallas, Texas hergestellt und vertrieben wird.In the case of the "600 Processing Package" and "700 ProcessingPackage "become CDP collective records processed in such a way that digital segment files are generated that are used in the process to be described here for generating shot point and receiver interference corrections can be used .. to carry out the here too. descriptive procedure can be any properly programmed digital computer are used, but it is preferred to use such a digital computer as manufactured and sold under the designation TIAC 870A by Texas Instruments Inc. Dallas, Texas will.

Fig.2 zeigt ein Funktionsflußdiagramm des hier zu beschreibenden Störkorrekturverfahrens. Die CDP-Sammelaufzeichnungen werden bei 20 in den Computer eingegeben und im Plattenspeicher gespeichert. In der bevorzugten Ausführungsform ist das Computermodell TIAC 870A so programmiert, daß es 2-bis 72-fach-CDP-Daten bearbeitet. Typischerweise werden im allgemeinen 6-bis 24-fach-Daten eingegeben. Jede der CDP-Spuren erstreckt sich seitlich von 1 bis 8 Sekunden, und sie enthält zwischen 250 und 2000 digitale Abtastpunkte, die bei typischen Fenstern mit einer Dauer von 4 Millisekunden genommen sind. Die CDP-SammelaufZeichnungsdaten sind gewöhnlich durch zeit-"2 shows a functional flow diagram of the here too descriptive interference correction procedure. The CDP collective records are entered into the computer at 20 and stored in disk memory. In the preferred In the embodiment, the TIAC 870A computer model is programmed to process 2 to 72-fold CDP data. Typically, 6 to 24 fold data are generally used entered. Each of the CDP lanes extends from 1 to 8 seconds laterally, and contains between 250 and 2000 digital sampling points that are taken in typical windows with a duration of 4 milliseconds. the CDP collective record data are usually defined by time "

309849/0469"309849/0469 "

abhängiges, bandbegrenzendes Filtern vorverarbeitet worden, das die Signal-Rausch-Bearbeitung nach dem hier zu beschreiebenden Verfahren unterstützt «Zusätzlich sind die CDP-Samm'elaufzeichnungdatteien durch eine Leistungsentzerrung vorverarbeitet worden, die dazu beiträgt, zu einem Schußpunkt-oder Einpfängerort gehörige Datenunregelmäsßigkeiten zu reduzieren. . 'dependent, band-limiting filtering has been preprocessed, which the signal-to-noise processing after that to be described here Procedure supported «In addition, the CDP collective recording files are through a power equalization has been preprocessed, which contributes to data irregularities associated with a shot point or receiver location to reduce. . '

Die digitalen Segmentdateien werden bei 22 eingegeben oder/ berechnet und in weederauffindbarer Form im Plattenspeicher des Computers gespeichert. Wie oben erwähnt wurde, können die Segmentdateien entweder aus dem ⁿ6OO Processing Package" oder aus dem "700 Processing Package" erhalten werden, die von Geophysical Services, Inc. vertrieben werden. DieSegmentdateien enthalten geordnete Schätzwerte für die Ankunftzeit, die Amplitude und den Moveout-liert eines Ereignisses. Mit "Moveout¹¹ wird der Parameter bezeichnet, der die Zunahme der Laufzeit bedingt durch die Meßgeometrie ausdrückt. Bei 24 werden Hauptsegmente entsprechend voreingestellten Parametern ausgewählt, die vom Benutzer eingegeben werden. Dadurch wird ermöglicht, daß das System nur mit Daten höchster Qualität arbeitet, so daß die Endgenauigkeit der Systemausgabewerte optimiert wird. .The segment digital files are entered or / computed at 22 and stored in undetectable form in the computer's disk storage. As mentioned above, the segment files can be obtained from either the ⁿ 600 Processing Package or the 700 Processing Package sold by Geophysical Services, Inc. The segment files contain ordered estimates of time of arrival, amplitude, and moveout -lated an event. "Moveout ¹¹ " denotes the parameter which expresses the increase in the running time due to the measurement geometry. At 24, major segments are selected according to preset parameters entered by the user. This enables the system to only work with data of the highest quality, so that the final accuracy of the system output values is optimized. .

In Fig.3 ist die Auswahl der Hauptsegmente ausführlicher dargestellt. Die Segmentdatei wird bei 26 eingegeben. Bei 28 werden die Parameter der Segaente mit Geschwindigkeits-, Zeit- und Raurafenstern verglichen*, Die Segmente werden auch mit einem Minimallängenschwellwert verglichen. Die Fenster und der Schwellwert werden vom Benutzer eingegeben. Der Benutzer bestimmt das Geschwindigkeitsfenster durch Untersuchen einer typischen Geschwindigkeitsinformation für den bestimmten betroffenen Ort. Eine solche InformationIn Fig. 3 the selection of the main segments is more detailed shown. The segment file is entered at 26. At 28 the parameters of the segments with speed, Time and space windows compared *, The segments are also compared to a minimum length threshold. The window and the threshold value are entered by the user. The user determines the speed window by examining typical speed information for the specific affected location. Such information

301149/0468'301149/0468 '

kann aus Kurvenstreudiagrammen öder Geschwindigkeitsstapeldiagrammen, die vorher festgestellt worden sind, erhalten werden. Die Zeit- und Raunifenster werden durch die besonders interessierende geophysikalische Zone bestimmt. Ein Beispiel eines typischen Längenschwellwerts sind 24 Tiefenpunkte. Jedes der bei 26 eingegebenen Segmente wird mit den. eingegebenen Fenster- und Längenschwellwerten verglichen, und diejenigen Segmente, die nicht innerhalb der Fenster liegen oder die Bedingungen des minimalen Längenschwellwerts nicht erfüllen, werden bei 28 gelöscht.can be made from curve scatter diagrams or speed stack diagrams, which have been previously determined can be obtained. The time and space windows are special by the geophysical zone of interest determined. An example a typical length threshold is 24 depth points. Each of the segments entered at 26 is compared with the. compared entered window and length threshold values, and those segments that are not within the window or the minimum length threshold conditions fail, are deleted at 28.

Die Segmente, die innerhalb der Fenster liegen und die die Bedingung des minimalen Längenschwellwertes erfüllen, werden bei 30 entsprechend dem Produkt aus der Schwellwertlänge und der Mittelamplitude des Segments eingeordnet. Die Segmente werden dabei entsprechend dem Produkt aus der Länge und der Amplitude .hintereinandergereiht, wobei das Segment mit dem größten Produkt an erster Stelle angeordnet wird.The segments that lie within the window and that meet the condition of the minimum length threshold, are classified at 30 according to the product of the threshold length and the mean amplitude of the segment. The segments are lined up according to the product of the length and the amplitude. placing the segment with the largest product first.

Die Einordnung der Segmente bei 30 kann auch entsprechend den Segmentlangen oder auch entsprechend den Segmentamplitüden erfolgen. Mit diesem Einordnungsschritt wird eine Tabelle aufgebaut, in der die qualitativ hochwertigsten Segmente an der Spitze stehen, so daß gewährleistet ist, daß die weitere Verarbeitung unter Verwendung der besten zur Verfügung stehenden Daten erfolgt.The classification of the segments at 30 can also correspond to the segment lengths or also according to the segment amplitudes take place. With this classification step a Table in which the highest quality segments are at the top, so that it is guaranteed that further processing is carried out using the best available data.

Das oberste .dinge ordne to Hauptsegment wird bei 32 aus der geordneten Datei entnommen, und bei 34 wird eine Entscheidung gefällt,'ob die maximale Segmentzahl, die vom Computer behandelt werden kann, bearbeitet worden ist oder nicht. Diese maximale Zahl wird von der SpeicherkapazitätThe top .things sort to main segment is at 32 taken from the ordered file and at 34 a Decision made as to whether the maximum number of segments that can be handled by the computer has been processed or not. This maximum number is determined by the storage capacity

309849/CU69309849 / CU69

. - 13 -. - 13 -

23230712323071

des Computers bestimmt. Wenn die maximale Segmentzahl nicht gewählt -worden ist, wird das als nächstes eingeordnete Segment in der geordneten Datei bei 36 genommen.of the computer. If the maximum number of segments has not been chosen, the next one will be ranked Segment taken in the ordered file at 36.

Bei 38 wird entschieden, ob das zweite Segment eine örtliche Redundanzanforderurig erfüllt.Der örtliche Redundanztest soll sicherstellen, daß Segmente so ausgewählt werden, daß «ie im wesentlichen die gesamte seismische Trasse erfassen. Wenn also bei 38 entschieden wird, daß bereits ein Segment gewählt worden ist,das einen Teil der seismischen Trasse erfaßt, der von dem bei 38 gerade getesteten Segment ebenfalls erfaßt wird, dann erfüllt das gerade getestete Segment die örtliche Redundanzanforderung nicht, und es wird gelöscht. Der Benutzer kann auf Wunsch jedoch einen gewissen Prozentsatz an Segmentüberlappung zulassen. Wenn das Segment die örtliche Redundanzanforderung erfüllt, wird das Segment bei 40 ausgewählt, und bei 42 .wird entschieden, ob alle der eingeordneten Segmente in der Datei getestet worden sind oder nicht. Wenn dies nicht der Fall war, dann werden die Schritte 34 bis 40. wiederholt, bis ' alle eingeordneten Segmente oder eine maximale Anzahl von gewählten Segmenten getestet worden sind, worauf die ausgewählten eingeordneten Segmente ausgegeben werden. · -At 38 a decision is made as to whether the second segment meets a local redundancy requirement. The local Redundancy testing is intended to ensure that segments are selected so that essentially the entire Record seismic route. So if it is decided at 38 that a segment has already been selected, the records a part of the seismic route which is also covered by the segment just tested at 38, then the segment being tested does not meet the local redundancy requirement and it is deleted. Of the However, the user can allow a certain percentage of segment overlap if desired. When the segment meets the local redundancy requirement, this will be Segment selected at 40, and a decision is made at 42 whether all of the classified segments in the file are tested have been or not. If this was not the case, then steps 34 to 40 are repeated until ' all ranked segments or a maximum number of selected segments have been tested, whereupon the selected classified segments are output. -

Nach Fig.2 v/erden die ausgewählten eingeordneten Segmente bei 44 geglättet. Ein solcher Glättungsvorgang kann eine Mittelung der Segmentzeit bei der Versetzung O und des Moveout-Werts bei der Bezugsversetzung mit ähnlichen, benachbarten Werten sein; er kann aber auch ein herkömmlicher gedichteter Mittelungevorgang sein. Der Glattunesvorgang kann vom Benutzer gesteuert sein, undAccording to FIG. 2, the selected classified segments are grounded smoothed at 44. Such a smoothing process can be an averaging of the segment time at the offset O and des Moveout value at the reference offset with similar, neighboring values; but it can also be a conventional one be a sealed averaging process. The smooth tuning process can be controlled by the user, and

3098^9/0^693098 ^ 9/0 ^ 69

J 14 ^J 14 ^

er wird dazu verwendet, Hochfrequenzzacken und dergleichen herabzusetzen. Diß geglätteten Segmente . werden dann bei 46 zur Erzeugung von Wellenzugkarten (wavelet-Karten) verwendet.it is used to reduce high frequency spikes and the like. This smoothed segments. are then used at 46 to generate wave train maps (wavelet maps) are used.

Fig.4 zeigt-im einzelnen die Verwendung der Segmente bei der Erzeugung von Xiellenzugkartenanordnungen. Die ausgewählten Segmente aus der geordneten Datei werden bei 4β eingegeben, und das erste Segment der Datei wird bei 50 ausgewählt. Die Segmentzeit T₀ bei der Versetzung O und der Moveout-¥ert Δ T^, bei der Bezugsversetzung werden bei 52 für den ersten gemeinsamen Tiefenpunkt des Segments erfaßt. Diese Parameter werden bei 54. zur Berechnung, der Ankunftzeit T für jede Spur des gemeinsamen Tiefenpunkts verwendet. Die .Berechnungen der Ankunftzeit T_v erfolgt mit Hilfe der folgenden Formel:FIG. 4 shows in detail the use of the segments in the generation of multiple train card arrangements. The selected segments from the ordered file are entered at 4β and the first segment of the file is selected at 50. The segment time T ₀ at the offset O and the moveout ¥ ert Δ T ^, at the reference offset are detected at 52 for the first common depth point of the segment. These parameters are used at 54th to calculate the arrival time T for each trace of the common depth point. The arrival time T _v is calculated using the following formula:

1/21/2

T— /"""* Φ j-f—\ ( A T¹ Λ. O AT¹ Φ ^ "X ίΛ\ T— / """* Φ jf— \ ( AT ¹ Λ. O AT ¹ Φ ^" X ίΛ \

^J v ^J v

wobei gilt: T_Q = Segmentzeit bei der Versetzung Null j 'X= Versetzung für die bestimmte Spur,where: T _Q = segment time at zero offset j 'X = offset for the specific track,

X = Bezugsversetzung entsprechend, dem Move-X = reference offset according to the move

out-Wert bei der Bezugsversetzung und /T= Moveout-Tiert bei der Bezugs Versetzung.out value at the reference offset and / T = Moveout-Tiert at the reference offset.

Nach der Berechunc von Τ_χ für jede Spur des gemeinsamen Tiefenpunkts wird ein auf Τ_χ ausgerichteter oder zentrierter Wellenzug (wavelet)bei 56 aus jeder Spur aus dem Speicher genommen und für die V/eiterverwendung in den Kernspeicher übertragen. In der bevorzugten Ausführungsform enthält an T/ellenzug 25 digitale Abtastwerte, und er überdeckt etwa 100 Millisekunden ,so daß die Erfassung von Ctörkorrekturen mit +50 Millisekunden ermöglicht wird.After Berechunc of _Τ χ for each track of the common depth point on an aligned or centered _Τ χ wave (wavelet) is taken at 56 from each track from the memory and transmitted using facilities in the core memory for the V /. In the preferred embodiment, T / ellenzug 25 digital samples, and it covers about 100 milliseconds, so that the detection is made possible by Ctörkorrekturen with +50 milliseconds.

3 0 9 8 k 9 / 0 k 6 93 0 9 8 k 9/0 k 6 9

Die "bei 56 entnommenen Wellenzüge sind somit um die Laufzeiten und die Moveout-Werte der "bei 2h ausgewählten Hauptsegmente zentriert. Diese Wellenzüge werden in der allgemei~ nen Form einer Wellenzugkarte längs einer Reihe für einen gemeinsamen Tiefenpunkt (CDP-Reihe) angeordnet.The wave trains taken at 56 are thus centered around the transit times and the moveout values of the main segments selected at 2h. These wave trains are arranged in the general form of a wave train map along a row for a common depth point (CDP row).

Die Figuren 5 "bis 8 tragen zur Erklärung der Bildung einer Wellenzugkarte bei. Eine ausführliche Beschreibung , von Wellenzugkarten findet sich in dem Aufsatz "Wavelt Maps: A new Analysis Tool for Reflection Seismograms" in der Zeitschrift Geophysics,XXXV, Nr. 3 ₉ Juni 1970, Seiten Ms 460.Figures 5 "to 8 help explain the formation of a wave train map. A detailed description of wave train maps can be found in the article" Wavelt Maps: A New Analysis Tool for Reflection Seismograms "in the journal Geophysics, XXXV, No. 3 ₉ June 1970, pages Ms 460.

Als Beispiel für die Bildung einer typischen Wellenzugrkarte wird ein typisches 6-fach-CDP-Aufzeichnungsverfahren verwendet. Bei einer solchen Aufzeichnungsgeometrie werden sechs Schußpunkte und 6 Empfänger nacheinander zu 6 verschiedenen Versetzungen bewegt, damit 6 Reflexionswellenzüge von einem einzigen Tiefenschichtpunkt aufgezeichnet v/erden. Die Reflexionszeiten der Wellenzüge beschreiben etwa eine hyperbolische Kurve 60 auf einer CDP-Sammelaufzeichnung, wie in Fig.5 dargestellt ist. Es ist zu erkennen, daß in Fig.5 der größte Teil der CDP-Spuren weggelassen ist, während für jede der 6 Spuren nur der Wellenzugabschnitt dargestellt ist.As an example of the formation of a typical wave access map becomes a typical 6-way CDP recording process used. With such a recording geometry six shot points and 6 receivers in succession to 6 different ones Dislocations moved so that 6 reflection wave trains are recorded from a single depth layer point. The reflection times of the wave trains describe, for example, a hyperbolic curve 60 on a CDP collective recording, as shown in Fig.5. It can be seen that most of the CDP tracks have been omitted in FIG. while only the wave train section is shown for each of the 6 tracks.

Zum Aufbau einer Wellenzugkarte für einen einzigen gemeinsamen Tiefenpunkt werden die- in Fig.5 dargestellten CDP-Spuren längs der geschätzten Hyperbel 60 abgetastet. Die einzelnen Wellenzüge werden dann unter ihren jeweiligen Versetzungszahlen horizontal ausgerichtet und dann um ihre geschätzten Reflexionszeiten zentriert,-wie in Fig.6 dargestellt ist. Durch den in Fig.6 dargestellten Ausrichtvorgang werden derTo build a wave train card for a single common At the bottom point, the CDP tracks shown in FIG. 5 are scanned along the estimated hyperbola 60. The single ones Wave trains are then aligned horizontally under their respective dislocation numbers and then around their estimated ones Reflection times centered, -as shown in Fig.6. By the alignment process shown in Fig.6, the

309849/0469309849/0469

• - 16 -• - 16 -

2S230712S23071

hyperbolische Normalmoveout sowie strukturelle Einflüsse beseitigt. Das heißt mit anderen Worten, daß durch die Darstellungsform der Wellenzugkarte alle Zeit-, Geschwindigkeits- und Geologieeinflüsse 'beseitigt werden, so daß nur der Einfluß von Störungen und andere durch Rauschen bedingte Fehler bleiben. Der in den Figuren 5 und 6 dargestellte Vorgang wird dann für die verbleibenden gemeinsamen Tiefenpunkte in der CDP-Sammelauf zeichnung wiederholt, und es kann für alle gemeinsamen Tiefenpunkte , die einen von dem Segment bestimmten Tiefenschichtrefiektor bilden, eine Wellenzugkarte aufgebaut werden.■hyperbolic normal moveout and structural influences eliminated. In other words, through the Form of representation of the wave train map all time, speed and geological influences' are eliminated so that only the influence of disturbances and other errors caused by noise remain. The one shown in Figures 5 and 6 The process is then repeated for the remaining common depth points in the CDP collective record. and it can be for all common depth points that form a depth layer reflector determined by the segment, a wave train map can be set up ■

Fig.7 zeigt beispielsweise den Verlauf von vier hyperbolischen Kurven, durch die Reflexionszeiten der Wellenzüge für vier gemeinsamen Tiefenpunkte, die in einer CDP-Sammelauf zeichnung enthalten sind. Diese hyperbolischen Kurven v/erden durch Verschieben und Summieren der CDP-Spuren längs anderer hyperbolischer Kurven und durch Untersuchen der gestapelten Kurven auf Maximalamplituden Scheitelwerte bestimmt, wie in dem Aufsatz " Estimation of Continuious Reflection Parameters From Seismic Data" von Stephen Merdler und anderen beschrieben ist, der auf der.38. Internationalen SEG-Jahrestagung, Denver, Colorado, 1968 präsentiert worden ist.7 shows, for example, the course of four hyperbolic ones Curves, through the reflection times of the wave trains for four common depth points, which are in a CDP collection drawing are included. These hyperbolic curves are grounded by shifting and summing the CDP tracks along other hyperbolic curves and by examining the stacked curves for maximum amplitude peaks determined, as described in the article "Estimation of Continuous Reflection Parameters From Seismic Data "by Stephen Merdler and others on der.38. International SEG Annual Meeting, Denver, Colorado, 1968 was presented.

Jeder längs der hyperbolischen Kurven von Fig.7 liegende Wellen-"¹ZUg wird dann horizontal längs seiner zugehörigen· Versetzungszahl ausgerichtet und auf die geschätzte Reflexionszeit zentriert, damit in der beschriebenen Weise die Einflüsse von Zeit, Geschwindigkeit und Geologie beseitigt werden. Diese Ausrichtung ist in Fig.8 dargestellt, die eine typische Wellenzugkarte zeigt. Wenn aus einer großen Anzahl von Punkten gemeinsamer Tiefe relativ große Wellenzugkarten gebildet werden, kann eine Karte ^{Each wave- "1} ZUg lying along the hyperbolic curves of FIG. 7 is then aligned horizontally along its associated offset number and centered on the estimated reflection time, so that the influences of time, speed and geology are eliminated in the manner described. This alignment is is shown in Fig. 8, which shows a typical wave train map: If relatively large wave train maps are formed from a large number of points of common depth, a map

3 09849/OA6 93 09849 / OA6 9

aufgebaut v/erden, die eine mehrdimensionale Darstellung von Wellenzügen in Abhängigkeit von vier unabhängigen Variablen darstellt. Beispielsweise können die Wellenzüge in horizontaler Richtung als Reflexionen ausgerichtet v/erden, die infolge eines Tiefenschichtpunkts auf dem Reflektor auftreten. Dann können Wellenzüge in vertikaler Richtung ausgerichtet werden, die der gemeinsamen Schußpunkt-Empfänger-Entfernung oder -Versetzung angehören. Wellenzüge, die längs einer ersten Diagonalen liegen, werden vom gleichen Schußpunkt erzeugt,während die längs der entgegengesetzten Diagonalen liegenden Wellenzüge am gleichen Empfängerort aufgezeichnet werden. In einer solchen"Karte sind die die Reflexionswellenzüge in dein definierten mehrdimensionalen Raum beeinflussenden Variablen der Punkt gemeinsamer Tiefe, die Versetzung, die Schußpunktposition und die.Empfängerposition.built v / earth that is a multi-dimensional representation of wave trains as a function of four independent variables. For example, the wave trains aligned in the horizontal direction as reflections, those due to a depth layer point on the reflector appear. Then wave trains can be aligned in the vertical direction, that of the common shooting point-receiver distance or relocation. Wave trains that lie along a first diagonal are taken from the same shot point is generated, while the one along the opposite one Wave trains lying diagonally at the same receiver location are recorded. In such a "card are the the reflection wave trains in your defined multidimensional Variables influencing space the point of common depth, the offset, the weft point position and the.recipient position.

Die bei 56 entnommenen Wellenzüge können nach Fig.4 bei in einer oben beschriebenen Wellenzugkarte angeordnet werden, in der die Wellenzüge längs Diagonalen liegen, die den gleichen Schußpunkt-und Empfängerpositionen.entsprechen. In der bevorzugten Ausführungsform des hier beschriebenen Verfahrens werden die bei 56 entnommenen Wellenzüge jedoch mit ihren entsprechenden Schußpunkt-und Zmpfangerzahlen gekennzeichnet und der Reihe nach längs der Reihen für die gemeinsamen Tiefenpunkte im Kernspeicher angeordnet. Die Kennzeichnung der Wellenzüge ermöglicht eine Flexiblität bei der Behandlung nicht gleichförmiger Sprengverfahren u. dgl. Für die hier vorliegende Beschreibung ist der Ausdruck Wellenzugkartenanärdnung so definiert, daß er entweder die mehrdimensionale Darstellung der Wellenzüge in Abhängigkeit von den vier oben erwähnten unabhängigen Variablen oder die Kennzeichnung der Wellenzüge mit den Schußpunkt- und Empfängerpositionen und das Anordnen in einer ReiheThe wave trains removed at 56 can be seen in FIG be arranged in a wave train map described above, in which the wave trains lie along diagonals that the correspond to the same shot point and receiver positions. In the preferred embodiment of the one described here However, the wave trains removed at 56 are processed with their corresponding shot point and receiver numbers marked and in turn along the rows for the common depth points arranged in the core memory. The marking of the wave trains enables flexibility in the treatment of non-uniform blasting processes and the like. For the purposes of the description herein, the term Wave train map definition so that it is either the multi-dimensional representation of the wave trains as a function of of the four independent variables mentioned above or the identification of the wave trains with the weft point and receiver positions and arranging in a row

30 984 9/046930 984 9/0469

für einen gemeinsamen Tiefenpunkt umfaßt. Die gruppierten. Wellenzüge »werden zum Plattenspeicher des Computers übertragen , in dem sie in wiederauffindbarer Form gespeichert werden. In der bevorzugten Ausführungsform enthält jeder Wellenzug 30 digitale Wörter mit 5 Wörtern als Kennzeichnungsvorsatz und 25 Wörtern, die Abtastpunkte darstellen. for a common depth point. The grouped. Wave trains »are transferred to the computer's hard disk, where they are stored in a retrievable form will. In the preferred embodiment, each waveform contains 30 digital words with 5 words as prefixes and 25 words representing sample points.

Bei 62 von Fig.4 wird entschieden, ob von dem gerade bearbeiteten Segment weiterePunkte gemeinsamer Tiefe erfaßt werden oder nicht. Wenn dies nicht der Fall ist, wird die Wellenzugkartenanordnung bei 64 ausgegeben. Wenn weitere Punkte gemeinsamer Tiefe von dem gerade bearbeiteten Segment erfaßt werden, dann werden die Werte T und AT für den nächsten gemeinsamen Tiefenpunkt bei 66 erfaßt, und die Schritte 54 bis 62 werden erneut durchlaufen. Nach der Ausgabe der Wellenzugkartenanordnung bei 64 wird bei 68 entschieden, ob noch weitere Segmente vorhanden sind. Sind keine weiteren Segmente vorhanden, dann fährt das Programm mit der in Fig.2 dargestellten ¹¹ Smash"-Wahlmöglichkeit fort. Wenn die Entscheidung ein positives Ergebnis hat, werden bei 50 zusätzliche Segmente ausgewählt, und die nachfolgenden Schritte des Verfahrens zur Erzeugung der Wellenzugkarte werden erneut durchlaufen.At 62 of FIG. 4, a decision is made as to whether or not further points of common depth are to be detected from the segment which has just been machined. If not, the wave train card assembly is issued at 64. If further points of common depth are detected from the segment being machined, then the values T and AT for the next common depth point are detected at 66 and steps 54 through 62 are repeated. After the output of the wave train card arrangement at 64, a decision is made at 68 as to whether there are any further segments. If there are no more segments, the program continues with the ¹¹ smash "option shown in FIG run through again.

Nach dem Aufbau der Wellenzugkartenanordnung wie bei 46 ist nach Fig.. 2 bei 70 eine Wahlmöglichkeit hinsichtlich einer Verarbeitung der Wellenzüge nach einem Verfahren vorgesehen, das mit "Wellenzugkarten-Smash" bezeichnet ist. Wenn der Benutzer von dieser Möglichkeit Gebrauch machen will, v/erden die Wellenzugdaten bei 72 entsprechend dem in Fig.9 angegebenen Prozeß verarbeitet.After building the wave train card assembly as at 46 is according to Fig. 2 at 70 an option with respect to a processing of the wave trains is provided according to a method called "wave train card smash" is. If the user wishes to make use of this possibility, the wave train data are correspondingly grounded at 72 processed by the process indicated in FIG.

9/04899/0489

Nach Fig.9 werden die Wellenzugkartenanordnungen bei 74 eingegeben, und die erste Wellenzugkartenanordnung wird bei 76 zur Verarbeitung ausgewählt. Bei 78 wird entschieden, ob das beteiligte Segment ein Wellental war oder nicht. Liegt ein Wellental vor, dann wird die Polarität der ¥ellenzüge in der Wellenzugkartenanordnung umgekehrt, damit die Wellenzüge in Scheitelwertanordnungen geändert werden. Bei 82 wird eine kumulative Wellenzugkarte aufgebaut, indem die Wellenzüge in jeder Wellenzugkartenanordnung gestapelt werden, die von der gleichen seismischen Trasse kommen, jedoch auf unterschiedliche Zeiten ausgerichtet sind,' Bei 84 wird entschieden,ob jede der Wellenzugkartenanordnungen verwendet worden ist, .oder nicht, Wenn dies nicht der Fall war, dann werden die Schritte 76 bis 82 erneut durchlaufen, bis eine vollständige kumulative Wellenzugkarte durch Stapeln aller Wellenzüge innerhalb aller eingegebenen Wellenzugkartenanordnungen aufgebaut worden ist.According to Figure 9, the wave train card assemblies at 74 is entered and the first wave train card assembly is selected at 76 for processing. At 78 will decided whether the segment involved had a trough was or not. If there is a wave trough, then the polarity of the ellipses in the wave train map arrangement vice versa, so that the wave trains in peak value arrangements be changed. At 82 a cumulative wave train map is built by dividing the wave trains in each Wave train map arrangement are stacked that come from the same seismic route, but on different Times are aligned, 'At 84 it is decided whether each of the wave train card assemblies has been used, . or not, If this was not the case, then steps 76 to 82 are run through again until a complete Cumulative wave train map by stacking all wave trains within all entered wave train map arrangements has been built.

Da jeder Wellenzug, der bei 82 an einem bestimmten Punkt gestapelt wird, von der gleichen seismischen Trasse kommt, hat jeder Wellenzug die gleiche störende statische Zeit-; verschiebung. Das Verfahren "Wellenzugkarten-smash" ergibt somit eine Wellenzugkartenanordnung mit den gleichen Zeitverschiebungen wie die einzelnen Wellenzugkartenanordnungen, die gestapelt worden sind, wobei jedoch wegen des Stapelvorgangs ein Vorteil hinsichtlich des RauschabStandes vorliegt. Die sich ergebende kumulative Wellenzugkarte, die bei 86 ausgegeben wird, enthält eine einzige Reihe von gemittelten WelLenzügen in einer Folge für einen gemeinsamen Tiefenpunkt. Wenn von der Wahlmöglichkeit der Anwendung des Verfahrens «We llenzugkarten-Snash". Gebrauch gemacht wird, dann wird bei den nachfolgenden Schritten des Prozesses nur eine einzige kumulative Wellenzugkartenanordnung und nicht dieBecause each wave train that at 82 at some point is stacked, comes from the same seismic route, each wave train has the same disturbing static time; shift. The procedure "wave train map smash" results thus a wave train card arrangement with the same time shifts like the individual wave train card assemblies that have been stacked, but because of the stacking process an advantage in terms of the noise ratio is present. The resulting cumulative wave train map output at 86 contains a single series of averaged wave trains in a sequence for a common Depth point. If use is made of the option of using the "Wave train ticket snash" procedure, then, in subsequent steps of the process, only a single cumulative wave train card arrangement will be used and not the

30.9 84 9/048930.9 84 9/0489

-20- ■ , -20- ■,

größere Anzahl der Wellenzuskartenanordhungen verwendet. Die Anwendung des Verfahrens "Wellenzugkarten-smash" . ergibt somit Zeiteinsparungen, während der nachfolgenden Verarbeitungsschritte sowie einen verbesserten Rauschabstand, ι ·larger number of shaft additional card arrangements used. The application of the "wave train card smash" procedure. thus results in time savings during the subsequent processing steps as well as an improved signal-to-noise ratio, ι ·

Wenn von der Möglichkeit der Anwendung des Verfahrens "Welienzugkarten-smash" kein Gebrauch gemacht worden ist, wird nach Fig.2 bei 88 die erste Wellenzugkartenanordnung ausgewählt. Bei90 werden alle Wellenzüge mit gemeinsamen Empfängerorten summiert oder gestapelt, damit eine Gruppe von gemittelten Wellenzügen für einen gemeinsamen Empfängerort erzeugt wird, die zur Erzeugung von Eiopf anger Störgrößen verarbeitet v/erden kann. Alle Wellenzüge innerhalb der ausgewählten Kartenanordnung oder innerhalb der sich beim Verfahren "Wellenzugkarten-smash" ergebenden Wellenzugkarte, die gemeinsame Schußpunktorte aufweisen, werden dann bis 92 summiert oder gestapelt. Der beim Schritt 92 ausgeführte Summierungsschritt ergibt eine Gruppe von gemittelten¹ Wellenzügen für einen gemeinsamen Schußpunkt, de zur Erzeugung von Schußpunktstörgrößen verarbeitet werden können. Die EmpfangerStörgrößen werden bei 9h erfaßt, während die Schußpunktstörgrößen bei 96 erfaßt werden. Zum Erfassen der EmpfangerStörgrößen und der Schußpnnktstörgrößen wird der gleiche Prozeß angewendet; er ist im einzelnen in Fig. 1-0 dargestellt.If the possibility of using the "wave train card smash" method has not been used, the first wave train card arrangement is selected at 88 according to FIG. At 90, all wave trains with common receiver locations are summed or stacked to generate a group of averaged wave trains for a common receiver location that can be processed to generate small disturbances. All wave trains within the selected card arrangement or within the wave train map resulting from the "wave train card smash" method which have common weft point locations are then added up to 92 or stacked. The summing step carried out in step 92 results in a group of averaged ¹ wave trains for a common shot point, which can be processed to generate shot point disturbances. The receiver perturbations are detected at 9h , while the shot point perturbations are detected at 96 . The same process is used to detect the receiver disturbance and the shot point disturbance; it is shown in detail in Fig. 1-0.

Nach Fig.10 wsrden bei 100 entweder die gemittelten Wellenzuggruppen für einen gemeinsamen Schußpunkt (ACS) oder die gemittelten Wellenzuggruppen für einen gemeinsamen Empfänger (ACR) eingegeben. Die vorherige Störgröße (PS) und die vorherige Amplitude (PA) werden bei 102 auf Null gesetzt. Bei 104 wird der erste zu erfassende Wellenzug ■ ausgewählt, und bei 106 werden alle Spitzenwerte des Wellen-According to Fig. 10, at 100 either the averaged Wave train groups for a common weft point (ACS) or the averaged wave train groups for a common one Recipient (ACR) entered. The previous disturbance (PS) and the previous amplitude (PA) become zero at 102 set. At 104 the first wave train ■ to be detected is selected, and at 106 all peak values of the wave-

30984B/0.A8930984B / 0.A89

zugs festgestellt, die den Anforderungen hinsichtlich der maximalen Störgröße, der maximalen relativen Störgröße und der minimalen Amplitude erfüllen.. Die Anforderungen hinsichtlich der maximalen Störgröße , der maximalen relativen Störgröße und der minimalen Amplitude werden zur Beseitigung unbedeutender Spitzenwerte vom Benutzer eingegeben. Bei 108 wird entschieden, ob Spitzenwerte festgestellt -worden sind, die die minimalen Eingangsparameteranforderungen erfüllen. Wenn solche Spitzenwerte festgestellt worden sind, wird bei 110 unter Verwendung der folgenden Formel der richtige Spitzenwert erfaßt:Zugs found that the requirements regarding the maximum disturbance variable, the maximum relative disturbance variable and the minimum amplitude. The requirements with regard to the maximum disturbance variable, the maximum relative disturbance variable and the minimum amplitude are to be eliminated insignificant peaks entered by the user. At 108 it is decided whether peak values have been determined, that meet the minimum input parameter requirements. If such peaks have been detected, will at 110, the correct peak value is detected using the following formula:

TS. TS. W. = A₁ ( 1 --^ _+i>3JW) ) (2)TS. TS. W. = A ₁ (1 - ^ _{+ i} > 3JW)) (2)

wobei gilt: W. = ein bewerteter Wertfaktorwhere: W. = an evaluated value factor

A. = Amplitude des untersuchten Spitzenwertes, PSP = eine vom Benutzer eingegebene Spitzenwert-A. = amplitude of the examined peak value, PSP = a peak value entered by the user

auswahlparameter-Konstante (normalerweise 5) ■ .· undselection parameter constant (usually 5) ■. · And

TS₁ = Absolutwert der vorherigen Störgröße abzüglich der Erscheingungszeit des untersuchten Spitzenwerte.TS ₁ = absolute value of the previous disturbance minus the time of occurrence of the examined peak value.

Bei Anwendung der Gleichung (2) -werden alle Spitzen, deren Amplitude kleiner als die Amplitude der zuvor^: erfaßten Spitze ist, gelöscht. Bei 1.10 wird für jede bei 106 festgestellte Spitze ein Wert für den Wertfaktor W. berechnet. Die Spitze mit dem höchsten Wert des Wertfaktors W. wird als gültige Spitze ausgewählt, damit die Störkorrektur für den bestimmten gerade, untersuchten Wellenzug bestimmt wird.In applying the equation (2) all peaks whose amplitude smaller than the amplitude of the previously ^-Be: detected peak is deleted. At 1.10, a value for the value factor W is calculated for each peak found at 106. The peak with the highest value of the value factor W. is selected as the valid peak, so that the interference correction is determined for the particular straight, investigated wave train.

309849/0489309849/0489

Das Zeitintervall zwischen der ausgewählten Spitze jedes Wellenzugs und der NuHachse des Wellenzugs wird dann bei 112 als Störgröße für diesen Wellenzug festgelegt, Die auf diese Weise bestimmte Störkorrekturgröße wird als die vorherige Störgröße (PS) eingestellt^fumd im Speicher abgespeichert. Die Amplitude der erfaßten Spitze wird ebenfalls als vorherige Amplitude (PA) abgespeichert.Wenn die' bei 108 zu fällende Entscheidung ein negatives Ergebnis hat, dann wird die vorherige Störgröße (PS) auf die Hälfte des bei 114 festgestellten Werts eingestellt, und die Störgröße wird auf Null gesetzt. Dieses Wiedereinstellen der vorherigen Störgröße trägt zur Feststellung gültiger Spitzen in aufeinanderfolgenden Wellenzügen bei. Bei 116 wird entschieden, ob noch zusätzliche Wellenzüge abgefragt werden sollen oder nicht. Sind noch abzufragende Wellenzüge vorhanden, dann werden die Schritte 104 bis 114 erneut durchlaufen. Wenn alle Wellenzüge erfaßt worden sind, wird das Unterprogramm zu den nachfolgenden Schritten von Fig.2 verlassen.The time interval between the selected peak of each wave train and the NuHachse of the wave is then determined at 112 as a disturbance variable for this wavetrain The thus determined Störkorrekturgröße is set as the previous disturbance variable (PS) ^f umd stored in memory. The amplitude of the detected peak is also stored as the previous amplitude (PA). If the decision to be made at 108 has a negative result, then the previous disturbance variable (PS) is set to half the value determined at 114 and the disturbance variable becomes set to zero. This resetting of the previous disturbance contributes to the determination of valid peaks in successive wave trains. At 116 it is decided whether or not additional wave trains are to be interrogated. If there are still wave trains to be interrogated, steps 104 to 114 are run through again. When all wave trains have been recorded, the subroutine for the following steps in FIG. 2 is exited.

Die Figuren 11 bis 13 zeigen drei Beispiele der Spitzenerfassung und der StörgrößenbeStimmung gemäß dem in Zusammenhang mit Fig. 10beschriebenen Verfahren. Fig.11 zeigt einen Wellenzug mit vier Spitzen P^ bis P^. Eine Berechnung nach der Gleichung(2)mit PSP= 5 ergibt die folgenden Viertfaktoren: W₁ = 6,8, W₂ = 92, \U = 27,6 und W^ = 1,8. Bs wird daher die Spitze P₂ ausgewählt, und für die nachfolgenden Berechnungen bezüglich der in den Figuren 12 und 13 dargestellten Wellenzüge .-v/erden die vorherige Störgröße (PS) auf -5 und die vorherige Amplitude (PA) auf 100 eingestellt.FIGS. 11 to 13 show three examples of peak detection and disturbance variable determination in accordance with the method described in connection with FIG. Fig.11 shows a wave train with four peaks P ^ to P ^. A calculation according to equation (2) with PSP = 5 results in the following fourth factors: W ₁ = 6.8, W ₂ = 92, \ U = 27.6 and W ^ = 1.8. The peak P ₂ is therefore selected and the previous disturbance variable (PS) is set to -5 and the previous amplitude (PA) is set to 100 for the subsequent calculations with respect to the wave trains shown in FIGS.

3Q9 849/Q4693Q9 849 / Q469

Fig. 12 zeigt vier Spitzen P,- bis Pg. Wenn die Werte der Spitzen P5 bis Pg in die Gleichung(2) eingesetzt werden, werden die Spitzen P,- und Pg gelöscht, weil ihre Amplituden kleiner als die Hälfte der vorherigen Amplitude von 1Op sind. Für die Spitzen Pg und P₇ ergeben sich die folgenden W^-Werte: Wg = 110,4 und W₇ = 45,6. -Fig. 12 shows four peaks P, - to Pg. If the values of the peaks P5 to Pg are substituted into equation (2), the peaks P, - and Pg are deleted because their amplitudes are less than half the previous amplitude of 1Op are. The following W ^ values result for the peaks Pg and P ₇ : Wg = 110.4 and W ₇ = 45.6. -

Somit ist zu erkennen, daß der Wertfaktor W5 am. größten ist, so daß die Spitze Pg ausgewählt wird; dabei wird eine Störgröße von -10 als vorherige Störgröße eingestellt. Die vorherige Amplitude wird auf 120 eingestellt.It can thus be seen that the value factor W5 is greatest is so that the peak Pg is selected; thereby will a disturbance variable of -10 is set as the previous disturbance variable. The previous amplitude is set to 120.

Fig.13 zeigt einen Wellenzug mit drei Spitzen Pq bis P^^. Die Spitze Pg wird gelöscht, da ihre Amplitude kleiner als die Hälfte der vorherigen Amplitude von 120 ist. Das Einsetzen der Werte für die Spitzen P₁₀ und P₁₁ in die Gleichung (2) ergibt die folgenden Wertfaktoren: '■ W₁₀ ~ 71, 2 und W₁₁ = 42_e . _x .Fig. 13 shows a wave train with three peaks Pq to P ^^. The peak Pg is deleted because its amplitude is less than half the previous amplitude of 120. The insertion of the values for the peaks P ₁₀ and P ₁₁ in the equation (2) yields the following value factors: '■ W ₁₀ ~ 71, 2, and W ₁₁ = 42 _e. _x .

Der Faktor W₁₀ ist am größten, so daß der Scheitelpunkt P₁ ausgewählt wird, und die vorherige Störgröße vlrd auf eingestellt, und die vorherige Amplitude wird auf 100 eingestellt. Fig,13 zeigt, daß mit dem hier beschriebenen Verfahren nicht notwendigerweise der Wellenzug mit dem größten Scheitelwert ausgewählt wird. Die in den in Fig.10 dargestellten Srfassungsvorgang verwendete Trendfolgelogik macht von der Zeitverschiebung zuvor ausgewählter Störgrößen als dem wahrscheinlichen Ort der derzeitigen Zeit- \ verschiebung Gebrauch, Diese Technik ermöglicht die zuverlässige Fassung „ von Störgrößen, die größer als die Hälfte der Hauptperiode.des Wellenzugs ist, und sie erlaubtdie Feststellung niederfrequenter Störgrößen.The factor W ₁₀ is the largest, so the vertex P _{1 is} selected, and the previous disturbance vlrd is set to, and the previous amplitude is set to 100. 13 shows that the method described here does not necessarily select the wave train with the largest peak value. The trend following logic used in the illustrated in Figure 10 Srfassungsvorgang makes use of the time difference previously selected disruptions as the probable location of the current time \ shift use, this technology allows the reliable version "of disturbances that are greater than half of the wave train Hauptperiode.des and it allows the detection of low-frequency disturbances.

3093-43/04*93093-43 / 04 * 9

Nach Fig.2 wird eine vorhandene Stqrgrößendatei nach der Erfassung der Empfänger- und Schußpunktgrößen bei 120 auf den neuesten Stand gebracht, und bei 122 wird entschieden, ob der unten zu beschreibende innere Schleifenvorgang vollendet worden ist oder nicht. Wenn der innere Schleifendurchlauf noch nicht vollendet worden ist, dann wird die bei 88 bestimmte Wellenzugkartenanordnung bei 124 entsprechend den zuletzt festgestellten Schußpunkt-Empfänger-Störgrößen korrigiert. Die störungskorrigierte Wellenzugkartenanordnung wird dann den-Verfahrensschritten 90 bis 120 ausgesetzt. Der Benutzer gibt die Zahl ein, die_ festlegt, wie oft die Innenschleife in der obengenannten Weise erneut durchlaufen werden soll. Wenn die vorgeschriebene Zahl von Innenschleifendurchläufen erreicht worden ist, hat der Entscheidungsschritt bei 122 ein positives Ergebnis und bei 126 wird entschieden, ob die Störgrößen aller existierenden Wellenzugkartenanordnungen geschätzt worden sind oder nicht.According to Figure 2, an existing Stqr size file after the Record the receiver and shot point sizes at 120 is brought up to date and a decision is made at 122 whether the inner looping process to be described below is completed has been or not. If the inner loop iteration has not yet been completed, then that determined at 88 is determined Corrected wave train card arrangement at 124 according to the last detected shot point receiver disturbance variables. The perturbation corrected wave train map assembly is then exposed to process steps 90 to 120. The user Enters the number that_ determines how often the inner loop should be run through again in the above-mentioned manner. When the prescribed number of internal loop passes has been reached, the decision step at 122 has a positive result and a decision is made at 126 whether the Disturbances of all existing wave train map arrangements have been estimated or not.

In der Praxis hat sich ein, viermaliger Durchlauf der Innen-' schleife als gut herausgestellt. Ein solcher Durchlauf ergibt eine wesentliche Verbesserung der Störgrößenschätzwerte für jede Wellenzugkartenanordnung.In practice, a 'four-time pass of the inside' loop turned out to be good. Such a run results in a substantial improvement in the disturbance variable estimates for any wave train card arrangement.

Wenn bei 126 festgestellt worden ist, daß die Störgrößen aus allen Wellenzugkartenanordnungen geschätzt worden sind, dann werden die Störgrößen aus jeder einzelnen Wellenzugkartenanordnung bei 130 gemittelt.Diese Mittelung ist eine bewertete Mittelung, bei der die Zahl der Wellenzüge, die in den ge-'mittelten Wellenzug eingeht, und die Amplitude des gemittelten Wellenzugs berücksichtigt werden. Bei 132 wird entschieden, ob eine nachfolgend zu definierenden "Außenschleife" vollendet worden ist oder nicht. Ein negatives Ergebnis der bei 132 gefällten Entscheidung führt dazu, daß die bei 130 festgestellte Störgröße auf alle eingegebenen CDP-Sammelaufzeichnungen bei 13^ angewendet wird.-Die CDP-Sammelaufzeichnungen werden dann mit " 600 oderIf it is determined at 126 that the disturbance is off all wave train card arrangements have been estimated, then the disturbance variables from each individual wave train card arrangement are estimated averaged at 130. This averaging is a weighted one Averaging, in which the number of wave trains in the determined Wave train is received, and the amplitude of the averaged wave train are taken into account. At 132 it is decided whether or not an "outer loop" to be defined below has been completed. A negative result of the decision made at 132 leads to that the disturbance determined at 130 is applied to all entered CDP collective records at 13 ^ CDP collective records are then marked with "600 or

3Q9 84 9/04S93Q9 84 9 / 04S9

- .25 -- .25 -

700 Processing Package ". bei 22 erneut bearbeitet, damit für die v/eitere Verarbeitung eine störkorrigierte Segmentdatei entsteht. Die Schritte 22 bis 130_werden in einer vom Benutzer vorbestimmten Zahl erneut durchlaufen.700 Processing Package ". Processed again at 22 so that an error-corrected segment file is created for further processing. Steps 22 to 130_ are in a run through the number predetermined by the user again.

Es hat sich gezeigt, daß ein je nach der Größe der Störung und der Anwesenheit von Niederfrequenzstörungen zwei - bis viermaliger erneuter Durchlauf der Außenschleife eine verbesserte Ankunftszeit und ein verbessertes Geschv/indigkeitsbild der Daten ergibt. Der wiederholte Durchlauf der Außenschleife ermöglicht auch, daß mit dein hier beschriebenen Verfahren sehr große Störungen dadurch korrigiert werden, daß ein Teil der großen Störungen während jeder AußensdiLeife elminiert wird. Nachdem die gewünschte Zahl der Außenschleifen vollendet worden'ist, wird bei 136 eine Störgrößendatei für die Verwendung bei den nachfolgenden Verarbeitungsschritten ausgegeben. . ■It has been shown that depending on the size of the disturbance and the presence of low frequency disturbances Two to four times running through the outer loop an improved arrival time and an improved The speed of the data. The repeated passage of the outer loop also enables that with your method described here very large Disturbances are corrected by eliminating some of the major disturbances during each outward bout will. After the desired number of outer loops has been completed, a disturbance file is created at 136 for use in the following processing steps issued. . ■

Somit ist gezeigt worden, daß mit dem hier beschriebenen Verfahren mit extremer Genauigkeit automatische Störgrößenkorrekturen an digitalen seismischen Daten für Punkte gemeinsamer Tiefe ausgeführt werden können. Das hier beschriebene Verfahren _kann vorteilhaft mit einem • richtig ,programmierten digitalen Computer ausgeführt werden, und es hat sich gezeigt, daß die Fehlerkorrekturen genauer als mit bisherigen Verfahren erfolgen, die die Kreuzkorrelation von Daten erfordern.'Thus it has been shown that with that described herein Process with extreme accuracy, automatic correction of disturbance variables can be performed on digital seismic data for points of common depth. That The method described here can advantageously be used with a • Properly, programmed digital computers are running and bug fixes have been shown more accurate than previous methods that require data to be cross-correlated. '

Pat e ntansnrüche Patent claims

3Q98497CH693Q98497CH69

Claims

Claims - -

1J Procedure for the determination of disturbance correction quantities for seismic Data with depth layer redundancy, characterized in that several wave trains from the seismic data (wavelets) that the wave trains are grouped in wave train map arrangements in which the Influences of time, speed and geology on the wave trains are eliminated that those on common shooting point positions related! Wave trains summed up that the wave trains related to common receiver positions are summed so that time intervals are determined at which selected peaks of the summed wave trains occur and that shot point and receiver clutter correction signals are generated in response to these time intervals.

2. The method according to claim 1, characterized in that

- Wave trains taken from the same seismic recording track before the summation of the wave trains for formation stacked on a single wave train card.

3. The method according to claim 1, characterized in that after the determination of the time intervals of each of the wave trains in the wave train card arrangements according to the shot point and receiver disturbance variables assigned to each wave train is postponed and that for education improved Disturbance correction variables, the summation and determination steps are repeatedly run through.

4.Verfahren according to claim 1, characterized in that according to the determination of the time intervals, the seismic data corresponding to the shot point and receiver interference correction variables moved \ erden and that the steps of claim 1 repeated be run through.

309849/048

ij. method according to. Claim 1, characterized in that the Wave trains selected depending on the main segments that will be obtained from a digital segment file can be selected from seismic data.

6. The method according to claim 5, characterized in that the main segments according to a predefined speed range to be selected. '

7 - .. The method according to claim 5, characterized in that the Main segments can be selected according to a predefined time range.

8. The method according to claim 5, characterized in that the

Main segments can be selected in accordance with a pre-defined • room area. ",

9 ^ The method according to claim 5, characterized in that the Main segments can be selected according to a predefined length threshold.

10. The method according to claim 5, characterized in that

the main segments smoothed before capturing the wave trains will.

11. Method for determining seismic disturbance correction quantities from collective records of seismic data for points common depth and from a digital segment file with estimates of the transit time and the speed of the seismic data, characterized in that several main segments are selected from the segment file, the Adhere to predetermined parameters that several wave trains are selected from the collective recordings, the

3-09849/04893-09849 / 0489

generally on the runtimes and the moveout values of the main segments are centered that the wave trains are put together in wave train card arrangements in which the Wave trains according to their estimated transit times are grouped along coordinates that correspond to the common Depth point and a common offset value correspond to the wave trains with common weft point positions be summed up that wave trains with common receiver positions that time intervals between selected peak values of the summed wave trains are summed and the center of the summed wave trains can be determined and that in dependence on these time intervals, shot point and receiver interference correction signals are generated will.

12. The method according to claim 11, c & characterized in that before the wave train summation wave trains that are of the same seismic trace are stacked to form a single wave train map, which is then summed.

13. The method according to claim 11, characterized in that after the determination of the time intervals of each of the wave trains in the wave train maps according to those assigned to each wave train Shot point and receiver interference is shifted, and that the summing and determining steps to achieve improved disturbance correction variables are run through repeatedly.

14. The method according to claim 11, characterized in that after the determination of the time intervals, the collective recordings according to the shot point and receiver interference correction variables be shifted, and that the steps of claim Ί1 be run through repeatedly.

309849/0469

β Method according to claim Hj, characterized in that the Main segments smoothed before picking out the wave trains ■ become. . ■ _

16. The method according to claim 11, characterized in that

the main segments are selected in accordance with a predetermined speed range j, a time range, a space range and a length threshold _β

17. The method according to claim 11, characterized in that at the determination step wave trains are eliminated whose Vertices are below a predetermined amplitude ratio with respect to the most recently determined wave train vertex, that the absolute time of each remaining Vertex with respect to the time of the last determined wave train vertex it is determined that the the determined absolute time is selected ρ if not another 'vertex has an amplitude that is greater than a predetermined one Size and that the time interval between the selected vertex and <fer center of the wave train is determined *

18. Method for determining seismic disturbance correction variables ' ^: from collective records of seismic data for points of common depth and from a digital segment file of seismic data in which wave trains with common shot point and receiver positions have been averaged, characterized in that the absolute .time of each vertex' one Wave train with respect to a reference time of the last determined wave train apex

• it is determined that the vertex of the wave train closest to the reference time is selected, if not another Γ vertex of the wave train has an amplitude which is greater than a predetermined value, and ■ · that the time interval between the selected vertex and the center of the wave train to determine the

309849/0469

. Disturbance correction variable is determined *. ■

19. Arrangement for the determination of disturbance correction quantities, for seismic data with spatial redundancy, characterized by storage facilities for retrievable storage digital collective recordings of seismic data and digital segment files with estimated values for the transit times and the velocities of the seismic data, selectors to select several single blocks from the Collective recording »where the wave trains determined Main segments correspond within the segment files, which have predetermined parameters, means for • Grouping the wave trains in time train cards so that the Influences of time, speed and geology eliminated are, sunraising devices for totaling.the.up wave trains related to common shooting point positions and wave trains related to common receiver positions " , Means for determining the time intervals between vertices of the summed wave trains and the midpoints of the summed wave trains and means for Creation of shot point and receiver interference correction. Signals that are proportional to the time intervals.

20. Arrangement according to claim 19, characterized by devices for stacking wave trains which are taken from the same seismic track _s so that a single wave train map is formed before the summation of the wave trains »

21. Arrangement according to claim 19, characterized by devices to move each of the wave trains in the wave train maps according to the shot point and receiver disturbance variables assigned to each wave train and means for carrying it out of repeating the summation and determination steps to achieve improved disturbance correction values.

309849/0489

22. Arrangement according to claim 19 _S, characterized by devices for shifting the seismic data in accordance with the determined visual point and receiver gate correction variables and devices for repeated actuation of the devices of claim 19 «

23. The arrangement according to claim 19, characterized _s that the predetermined parameters used to Aüsxfahl the main segments contain a selected speed range. -..- ■.

24 »Arrangement according to claim 19 _S characterized in that the. for the selection of the main segments, predetermined parameters used contain a selected time range.

25. The arrangement according to claim 19 _S, characterized in that the predetermined parameters used to select the main segments contain a selected Raumbereicii.

26. The arrangement according to claim 19 _S, characterized in that the predetermined parameters used to select the main segments contain a selected length threshold value.

27. Arrangement mch claim 1 ^: 9, characterized by devices for smoothing the main segments before the wave trains are picked out.

28. Arrangement according to claim 19 * characterized by devices to classify the main segments according to the segment length multiplied by the segment mean amplitude.

308849 / 0AS9