imported>Invisigoth67: typo

2026-04-29T06:14:08Z

typo

Neue Seite

[[Datei:Seismic CMP NMO.png|mini|300px|CMP-Anordnung und Laufzeitunterschiede (Normal Move Out)]]

Das ''' Processing''' ist die Aufbereitung [[Seismik|seismischer Daten]] bei [[Reflexionsseismik|reflexionsseismischen]] Untersuchungen. Sie wird durchgeführt, um aus den Rohdaten eines seismischen Experiments ein sinnvolles geologisches Abbild des Untergrundes zu erhalten. Bei einem seismischen Experiment werden Zweiwegelaufzeiten (''Two Way Travel Time, TWT'') aufgezeichnet, also die Zeit von der Abgabe des Impulses an der seismischen Quelle hinunter zum [[Seismischer Reflektor|Reflektor]] und zurück an die Erdoberfläche zum [[Geophon]]. Die Daten müssen so bearbeitet werden, dass aus den Zeitmessungen Informationen zur Tiefe der Reflektoren werden, um erfolgreiche [[Bohrung]]en (z. B. bei der [[Kohlenwasserstoff]]-[[Exploration (Geologie)|Exploration]]) oder tiefenbezogene Erkenntnisse in der [[geologisch]]en und [[geophysik]]alischen Grundlagenforschung zu ermöglichen.

== Preprocessing ==
Für die anschließende Bearbeitung werden die gemessenen Daten zunächst in [[Reflexionsseismik#Common-Midpoint-Technik und Normal Move Out|Common Midpoint-Gather]] sortiert.

Das daran anschließende '''Preprocessing''' umfasst die analoge [[Bandpassfilter]]ung des Signals, statische Korrekturen und die [[Dekonvolution]] des Signals.

=== Bandpassfilterung ===
Der Bandpassfilter besteht aus einem [[Hochpassfilter]] und einem [[Tiefpassfilter]]. Bei der [[Digitalisierung]] des analogen seismischen Signals kann es zu [[Alias-Effekt]]en kommen. Dazu kommt es, wenn das [[Nyquist-Shannon-Abtasttheorem|Abtasttheorem]] verletzt wurde und man eine zu niedrige Abtastfrequenz gewählt hat. Um dies zu verhindern setzt man einen analogen Tiefpassfilter vor der Abtastung des Signals ein. Außerdem wendet man einen Hochpassfilter an, um [[Rauschen (Physik)|Rauschen]], [[Netzbrummen]] und langperiodische Anteile von [[Seismische Welle#Oberflächenwellen|Oberflächenwellen]] zu unterdrücken.

=== Statische Korrektur ===
Durch die '''statische Korrektur''' wird das seismische Experiment auf ein Bezugsniveau reduziert, da die oberste Schicht des Untergrundes (Verwitterungsschicht) Höhenunterschiede und starke Geschwindigkeits[[Heterogenität (Naturwissenschaft)|inhomogenitäten]] aufweist. Zur Berechnung der statischen Korrektur benötigt man die [[Mächtigkeit (Geologie)|Mächtigkeit]] und die seismische [[Geschwindigkeit]] der Schicht zwischen dem Bezugsniveau und der Oberfläche. Diese lassen sich mittels [[Refraktionsseismik]] an den oberflächennahen Schichten oder einer [[Seismische Quelle|seismischen Quelle]] in einem [[Bohrloch]] ermitteln. In der marinen Seismik sind statische Korrekturen notwendig, um den Einfluss der Gezeiten sowie unterschiedliche Tiefen von Empfängern und Quelle auszugleichen.

== Geschwindigkeitsanalyse ==
{| class="wikitable float-right"
|-
!class="hintergrundfarbe6"| Modell
!class="hintergrundfarbe6"| Geschwindigkeitsinkrement
|-
| Horizontaler Reflektor
|style="text-align:center"| <math>v</math>
|-
| n horizontale Reflektoren
|style="text-align:center"| [[Effektivwert|<math>v_{RMS}</math>]]
|-
| Geneigter Reflektor
|style="text-align:center"| <math>\frac{v}{\cos\zeta}</math>
|-
| Beliebiges Modell
|style="text-align:center"| <math>\sqrt{\frac{1}{t_0\cdot\frac{\text{d}^2t}{\text{d}x^2}}}</math>
|}
Der Normal MoveOut (NMO) ist die Abweichung der Laufzeiten von der Einsatzzeit. Die Geschwindigkeitsanalyse muss vor der weiteren Bearbeitung der seismischen Daten erfolgen, um den Normal MoveOut zu bestimmen und zu korrigieren. Die [[Laufzeitkurve]] einer seismischen Welle ergibt sich aus
<math>t^2=t_0^2 + \frac{x^2}{v_{NMO}^2}</math>.
Die MoveOut-Geschwindigkeit <math>v_{NMO}</math> beinhaltet sowohl die Art der Welle als auch die Modellannahme (siehe Tabelle).

=== Geschwindigkeitsbestimmung ===
==== ''t''–Δ''t''-Methode ====
Die <math>t</math>-<math>\Delta t</math>-Methode basiert auf der [[Parabel (Mathematik)|parabolischen]] [[Näherung]]. Die zugrundeliegende Gleichung ist <math>2t_0\Delta t_n\approx\frac{x^2}{v_{NMO}^2}</math>. Nun trägt man <math>2t_0\Delta t_n</math> gegen <math>x^2</math> auf und kann aus der Steigung <math>\frac{1}{v_{NMO}^2}</math> die MoveOut-Geschwindigkeit berechnen.

==== ''t''<sup>2</sup>–''x''<sup>2</sup>-Methode ====
Die <math>t^2</math>-<math>x^2</math>-Methode ist genauer als die <math>t</math>-<math>\Delta t</math>-Methode, da sie auf der [[Hyperbel (Mathematik)|Hyperbelnäherung]] basiert. Man kann die Laufzeitkurve in ein <math>t^2</math>-<math>x^2</math>-Diagramm eintragen und die Geschwindigkeit bestimmen. Wie bei der <math>t</math>-<math>\Delta t</math>-Methode beträgt die Steigung der Gerade <math>\frac{1}{v_{NMO}^2}</math>.

==== Dynamische Korrektur ====
Die dynamische NMO-Korrektur reduziert die Laufzeitkurve offsetabhängig auf die Einsatzzeit. Für jedes CMP-Gather wird für jedes <math>t_0</math> eine NMO-Korrektur durchgeführt. Dies gibt eine grobe Abschätzung der seismischen Geschwindigkeiten im Untergrund, birgt allerdings Probleme. Die NMO-Korrektur ist eine nicht-lineare Streckung der Zeitachse entlang des Offsets. Entsprechend wird das Signal mit zunehmenden Offset gestreckt (NMO-Stretch) und die Frequenz wird gemindert.

Die Daten werden um den NMO-Stretch korrigiert, indem man einen linearen Filter anwendet, der die Anteile im Gather Null setzt. Die Daten sind nun bereit [[Rauschabstand#Mittelung|gestapelt]] zu werden.

==== Durchschnittsgeschwindigkeit ====
Um die Durchschnittsgeschwindigkeit zu ermitteln, versenkt man Geophone in ein Bohrloch und zündet eine seismische Quelle neben dem Bohrloch (''vertical seismic profiling''). Die Durchschnittsgeschwindigkeit bezeichnet den geometrisch kürzesten Weg von der Quelle zum Reflexionspunkt.

<math>v_a=\sqrt{\frac{\sum\limits_{k=1}^n z_k}{\sum\limits_{k=1}^n \frac{z_k}{v_k}}}</math>

==== Effektivgeschwindigkeit ====
Die Effektivgeschwindigkeit (''Root Mean Square'') bezeichnet den quadratischen Mittelwert der Geschwindigkeit.

<math>v_{RMS}=\sqrt{\frac{\sum\limits_{k=1}^n z_k\cdot v_k}{\sum\limits_{k=1}^n \frac{z_k}{v_k}}}</math>

==== Intervallgeschwindigkeit ====
Die Intervallgeschwindigkeit ergibt sich aus der Dix-Dürbaum-Krey-Formel:

<math> v^2_n=\frac{(v_{RMS}^n)^2\cdot t_0^n-(v_{RMS}^{n-1})^2\cdot t_0^{n-1}}{t_0^n-t_0^{n-1}}</math>

<math>v_n</math> ist die Geschwindigkeit der Schicht zwischen <math>t_0^n</math> und <math>t_0^{n-1}</math>. Für eine homogene Schicht ist <math>v_n</math> die Schichtgeschwindigkeit <math>v</math> und für einen geschichteten Bereich ist <math>v_n</math> die Effektivgeschwindigkeit <math>v_{RMS}</math>

==== Stapelgeschwindigkeit ====
Die Stapelgeschwindigkeit ergibt zusammen mit der Pseudoeinsatzzeit <math>t_{0s}</math> die Laufzeitkurve

<math>t_s=\sqrt{t^2_{0s}+\frac{x^2}{v^2_s}}</math>.

Diese Hyperbel nähert sich an die wahre Laufzeitkurve in der Art an, dass die Flächen der positiven und negativen Abweichung sich balancieren. Die Stapelgeschwindigkeit und die Pseudoeinsatzzeit variieren mit dem Offset der Laufzeitkurve.

=== Automatische Geschwindigkeitsanalyse ===
Die <math>t^2</math>-<math>x^2</math>-Methode muss von Hand ausgeführt werden, sie lässt sich weder automatisieren noch ist sie objektiv. Diese Ansprüche werden jedoch an moderne Processing-Algorithmen gestellt.

==== Constant Velocity Scan ====
Es werden konstante Geschwindigkeiten aus einem adäquaten Geschwindigkeitsbereich verwendet, um eine NMO-Korrektur anzuwenden. Anschließend wird für jede Reflexion die Geschwindigkeit ausgewählt, die die MoveOut-Zeiten auf die Einsatzzeit verschiebt. Dies Verfahren muss visuell ausgewertet werden und ist besonders bei schlechtem [[Rauschabstand]] sehr subjektiv.

==== Constant Velocity Stack ====
Beim Constant Velocity Stack oder Brute Stack werden wieder konstante Geschwindigkeiten aus einem Geschwindigkeitsbereich verwandt um die NMO-Korrektur auf das CMP-Gather anzuwenden. Anschließend werden die einzelnen Spuren des CMP aufsummiert. Durch destruktive Überlagerungen werden Laufzeitkurven, die mit der falschen NMO-Geschwindigkeit korrigiert werden, bei der Summation (Stapelung) ausgelöscht. Sind alle MoveOut-Zeiten in Phase mit der Einsatzzeit, ist die Amplitude der resultierenden Spur am größten und die gewählte Geschwindigkeit optimal.

Dieses Verfahren muss trotzdem noch visuell ausgewertet werden und unterliegt einer gewissen Subjektivität.

==== Geschwindigkeitsspektrum ====
[[Datei:Seismic Velocity Analysis.png|mini|Geschwindigkeitsanalyse mit [[Seismic Unix]]]]
Die Berechnung von Geschwindigkeitsspektren automatisiert die Auswahl der Geschwindigkeitsfunktion <math>v(t_0)</math>.

Für jede Einsatzzeit wird für eine Geschwindigkeit <math>v</math> die NMO-Korrektur angewandt. Daraufhin wird der Semblance-Koeffizient <math>S</math> für ein Zeitfenster um <math>t_0</math> berechnet.

<math>S=\frac{\sum\limits_{i=1}^M\left(\sum\limits_{j=1}^N a_{ij}\right)^2}{N\sum\limits_{i=1}^M\sum\limits_{j=1}^Na^2_{ij}}</math>

Hierbei bezeichnet N die Anzahl der Spuren im CMP-Gather und M die Anzahl der diskreten Werte in der j-ten Spur.
Der Semblance-Koeffizient normiert die Energie der eingehenden Spur auf die Energie aller Spuren, weswegen er nur Werte zwischen 1 und 0 annehmen kann. Diese Berechnung wird für alle <math>t_0</math> and alle <math>v_{NMO}</math> wiederholt und in ein <math>(v,t_0)</math>-Diagramm eingetragen. Man erhält die Funktion <math>S(v,t_0)</math> für das Geschwindigkeitsspektrum. Nun kann automatisiert für jede Reflexion der höchste Semblance-Koeffizient ausgewählt werden, der die optimale Stapelgeschwindigkeit charakterisiert.

Der Semblance-Koeffizient ist ein Maß für die [[Kohärenz (Physik)|Kohärenz]] und wird wegen der hohen Auflösung in Zeit und Geschwindigkeit verwendet.

=== Ergebnis der Geschwindigkeitsanalyse ===
Aus der Geschwindigkeitsanalyse und NMO-Korrektur der CMP-Gather mit anschließender Stapelung erhält man eine Zeitsektion aus Lotzeiten, auch Zero-Offset-Sektion oder Stapelsektion genannt. Diese simulieren einen vertikalen Strahleinfall von der Quelle zum Reflektor. Das Datenvolumen wird um den Faktor des Überdeckungsgrades des CMP-Gathers reduziert.

== Migration ==
[[Datei:SimpleSyncline.jpg|mini|Nichtmigrierte synthetische Zero-Offset-Sektion mit Triplikation.]]

'''Migration''' (von lateinischen: ''{{lang|la|migratio}}'', Wanderung, Übersiedlung) in der [[Seismik]] bezeichnet Verfahren zur Erstellung eines Untergrundbildes aus dem gemessenen [[Wellenfeld]] an der Oberfläche.

Ziel der Migration ist es, aus den [[Reflexion (Physik)|Reflexionen]] im Untergrund ein geologisches Abbild des Untergrundes mit den korrekten Neigungen, Längen und Positionen der Reflektoren zu erstellen.

=== Motivation ===
CMP-Gather können Triplikationen und andere kreuzende Schichten beinhalten, die geologisch nicht vorkommen. Zudem werden diese Zeiten in Tiefen migriert um exakte Bohrungen zu ermöglichen.

=== Migrationsarten ===
Man kann verschiedene Kriterien zur Unterscheidung der Migration heranziehen.

Zum einen wird zwischen der Zeitmigration in den <math>(x,t)</math>-Raum und der Tiefenmigration in den <math>(x,z)</math>-Raum unterschieden. Die Tiefenmigration reagiert dabei empfindlicher auf Fehler im Geschwindigkeitsmodell des Untergrundes, liefert im Gegensatz zur Zeitmigration dafür ein authentisches Abbild des Untergrundes. Zum anderen kann man zwischen der Prestack- und der Poststack-Migration unterscheiden. Die Poststack-Migration wird nach der Geschwindigkeitsanalyse und Stapelung gemacht. Die Prestack-Migration findet vor der Geschwindigkeitsanalyse und der Stapelung statt, um auch Amplituden der eingehenden Signale zu berücksichtigen.

Durch speziellere Prestack-Migrationen können Aussagen über die [[Reflexionsfaktor|Reflexionskoeffizienten]] im Untergrund getroffen werden.

== Siehe auch ==
* [[Laufzeitdiagramm]]
* [[Refraktionsseismik]]

== Literatur ==
* O. Yilmaz: ''Seismic data processing.'' Society of Exploration Geophysicists, Tulsa 2001, ISBN 0-931830-40-0.
* Dirk Gajewski: [http://www.geophysics.zmaw.de/fileadmin/documents/teaching/Skripten/Skript.pdf ''Vorlesungsskript: Angewandte Geophysik II.''] (PDF; 32,5 MB) Hamburg 2010, abgerufen am 26. Februar 2010.
* Claudia Vanelle: ''Vorlesungsskript: Migration.'' Hamburg 2009; [http://www.geophysics.zmaw.de/fileadmin/documents/StudiumLehre/scripts/MIG-01.pdf Kapitel 1] (PDF) [http://www.geophysics.zmaw.de/fileadmin/documents/StudiumLehre/scripts/MIG-02.pdf Kapitel 2] (PDF) [http://www.geophysics.zmaw.de/fileadmin/documents/StudiumLehre/scripts/MIG-03a.pdf Kapitel 3] (PDF); abgerufen am 26. Februar 2010.

[[Kategorie:Geophysik]]
[[Kategorie:Seismik]]

Processing (Seismik) - Versionsgeschichte

imported>Invisigoth67: typo