DE2631634A1

DE2631634A1 - Anordnung und verfahren zur bestimmung der position eines in der seismischen exploration eingesetzten seeschleppkabels

Info

Publication number: DE2631634A1
Application number: DE19762631634
Authority: DE
Inventors: Jun Robert Calvin Donald; Larry Wayne Moore; Jun Joseph Aaron Rice
Original assignee: Texas Instruments Inc
Current assignee: Texas Instruments Inc
Priority date: 1975-07-14
Filing date: 1976-07-14
Publication date: 1977-01-20
Also published as: FR2319911B1; NO155417B; DE2631634C2; JPS5238250A; AU1569176A; AU507411B2; NO762445L; FR2319911A1; US4068208A; NL187088C; NL7607784A; GB1553712A; CA1076702A; NL187088B; NO155417C

Description

Patentanwälte

Dipl.-Ing. Dipl.-Ohem. Dipl.-Ing.

E. Prinz - Dr. G. Hauser - G. Leiser

Ernsbergerstrasse 19 _

8 München 60 2631634

Unser^Zeichen; T 2052 9,JuIi 1976

Texas Instruments Incorporated
13500 North ^Central Expressway-Dallas, Texas, V.St.A,

Anordnung und Verfahren zur· Bestimmung der Position eines in der seismischen Exploration eingesetzten Seeschleppkabels

Die Erfindung bezieht sich auf die seismische Exploration auf See und insbesondere auf eine' Anordnung und ein Verfahren zum. Bestimmen der Position eines seismischen Seeschleppkabels«

In der seismischen Exploration auf See zieht ein Explorationsschiff ein seismisches Schleppkabel nach,in dem mehrere druckempfindliche Detektoren untergebracht sind, die üblicherweise als Hydrophone bezeichnet werden, Impulsquellen wie mit Ex~ plosivstoffen oder Luft arbeitende Impulsgeber, die ebenfalls vom Explorationsschiff geschleppt v/erden, v/erden so gezündet, daß auf das Wasser und schließlich auf die darunterliegende Erdkruste eine sich ausbreitende Energie übertragen wird. Da sich diese Energie in einer allgemein abwärts verlaufenden Richtung ausbreitet, v/erden Teile der Energie von Erdkrusten-Oberflächenunstetigkeiten reflektiert und schließlich von den Hydrophonen als Druckschwankungen empfangen, Schw/Ba

609883ΛΟ973

Die entsprechenden elektrischen Signale werden vom Schleppkabel zu einer Aufzeichnungsvorrichtung an Bord des Explorationsschiffs für die anschließende Verwendung bei der Interpretation der Struktur der Erdkruste übertragen. Damit solche Interpretationen aussagekräftig sind, ist es jedoch notwendig, den Ort des Explorationsschiffs zu der Zeit, an der die Daten erfaßt werden, mit gewisser Genauigkeit zu kennen. Zur Erzeugung dieser Information sind mehrere anspruchsvolle Navigationssysteme entwickelt worden. Eines dieser Systeme ist beispielsweise das von der Firma Texas Instruments Incorporated, Dallas, Texas hergestellte und für den Einsatz bei der seismischen Exploration bestimmte Navigationssystem GEONAV.

Solche Systeme können zwar zuverlässige Messungen über den Ort des Explorationsschiffs liefern, doch bleiben immer noch Mehrdeutigkeiten bei der Kenntnis der Lage der entfernten Teile des Schleppkabels bezüglich der vom Schiff befahrenen Bahn. Beispielsweise kommt es selten vor, daß das Schleppkabel genau längs der Bahn des Schiffs nachgezogen wird. Das Schleppkabel ist zwar an einem über der Viasseroberfläche befindlichen Punkt am Schiff befestigt, doch wird sein Hauptabschnitt, dessen Länge eine Meile oder mehr betragen kann, typischerweise in einer Nenntiefe unterhalb der Wasseroberfläche von mehreren Tief.enkontrollkörpern gehalten, die längs des Schleppkabels angeordnet sind. Ein solcher Tiefenkontrollkörper ist beispielsweise in der US-PS 3 372 666 beschrieben. Ein Grund dafür, daß das Schleppkabel nicht der Bahn des Schiffs nachfolgt, ist die Tatsache, daß die Querströmungsgeschwindigkeit in der Tiefe des Schleppkabels oft von der quer zur Bahn auftretenden Strömungsgeschwindigkeit abweicht, die auf den Rumpf des Schiffs selbst einwirkt. Weitere Faktoren, die zu diesem Problem

809883/097

beitragen, sind der Vorhaltewinkel des Schiffs und Strömungsgradienten längs der Bahn.

Bisher ist die Lage des Schleppkabels geschätzt worden, indem ein Radarreflektor oder Radartransponder an einer an der Oberfläche befindlichen Endboje verwendet wurde₅ die an dem vom Explorationsschiff entfernten Ende des Schleppkabels befestigt ist. Dies ermöglicht eine Abschätzung der Lage des entfernten Endes des Schleppkabels, jedoch sind die dabei auftretenden Ungenauigkeiten zu groß, als daß sie mit dem Auflösungsvermögen vereinbar wären-, das .mit modernen Auswertungsverfahren erreicht wird, die für die seismischen Daten selbst verfügbar sind. Außerdem wird bei dem Radarortungsverfahren für das Schleppkabel eine geradlinige Bahn Zwischen dem Heck des Schiffs und der Endboje angenommen.

Mit Hilfe der Erfindung sollen daher eine Anordnung und ein Verfahren zum Bestimmen des Orts eines seismischen Seeschleppkabels geschaffen werden. Dabei sollen Messungen angewendet werden, die zweckmässigerweise am seismischen Schleppkabel durchgeführt werden. Die. Kompliziertheit des Verfahrens und der Anordnung soll dabei auf die Genauigkeit der gewünschten Ergebnisse bezogen sein» Ferner soll mit Hilfe der Erfindung eine Vorrichtung geschaffen werden, mit der der Verlauf des dem Schleppschiff, benachbarten , Abschnitts des seismischen Schleppkabels bestimmt werden kann_o

In der anläßlich der zweiten jährlichen Off-Shore-Technologiekonferenz, Houston, Texas, vom 22. bis 24.April 1970 herausgegebenen Druckschrift Nr. 0TC1281 mit dem Titel "Three-dimensional Boundary Value Problems for Flexible Cables" von John W.Bedenbender, ist ein mathematisches Modell für die Gleichgewichtslage eines seismischen

609883/0973

Seekabels in Form von sechs Differentialgleichungen entwickelt worden. In der Druckschrift ist außerdem ein Yerfahren zum Lösen der Gruppe von Differentialgleichungen angegeben. Dieses Verfahren erfordert jedoch die Kenntnis einer großen Anzahl von Parametern^ von denen einige _s wenn überhaupt im Fall eines seismischen Schleppkabels nur unzureichend bekannt sind_o Bei Anwendung der .in der Druckschrift beschriebenen Yerfahren zur Untersuchung simulierter Fälle konnte beobachtet werden, daß die seitliche Verschiebung des Schleppkabels in enger Beziehung mit dem Gierungswinkel des Kabels in einem auf die Erde bezogenen Koordinatensystem (gemessen in der Horizontalebene) relativ zu der vom Schiff durchfahrenen Bahn steht.

Gemäß einem Merkmal der Erfindung wird daher eine besondere Vorrichtung zum Verbinden des seismischen Seeschleppkabels mit dem Explorationsschiff verwendete Der Körper der Verbindungsvorrichtung kann fest am Heck des Schiffsangebracht werden. Das Seeschleppkabel ist mit Hilfe eines Kreuzgelenks mit dem Körper der Verbindungsvorrichtung verbunden. In der bevorzugten Ausführungsform enthält die Vorrichtung zwei Drehmelder ₉ die mit dem Kreuzgelenk so zusammenwirken daß Signale geliefert werden, die den Gierungswinkel und den Längsneigungswinkel des an das Schiff angrenzenden Abschnitts des seismischen Schleppkabels repräsentieren. Diese Meßwerte des Gierungswinkels und des Längsneigungswinkels ergeben sich in einem auf dös Schiff bezogenen Koordinatensystem. Die Vorrichtung enthält außerdem Pendeleinheiten, die Signale liefern, die das Schlingern und die Längsneigung des Explorationsschiffs repräsentieren. Diese Signale werden mit den von den Drehmeldern gelieferten Signalen kombiniert, damit sich Messungen bezüglich der Gierung und der Längsneigung des Kabels in einem erdbezogenen Koordinatensystem ergeben. ¥ie später noch genauer erläutert wird, werden diese Messungen der

609883/0973

Gierung und der Längsneigung dazu verwendet, Modelle mit verschiedenen Graden der Kompliziertheit zur Darstellung der Lage des Schleppkabels zu entwickeln«

Bei einem dieser Modelle wird zusätzlich ein Magnetkompaß verwendet, der an einer vom Explorationsschiff entfernten Stelle am -Schleppkabel, angebracht ist. Das von dem Magnetkompaß gelieferte Signal repräsentiert nach einer Korrektur hinsichtlich der magnetischen Variation in dem Gebiet die seitliche Lage des Schleppkabels am Ort des Kompaß. Das Modell ist eine mit Hilfe von zwei Linien erhaltene Näherung der Lage des Schleppkabels. Die erste Linie, die durch das Heck des Explorationsschif£s verläuft, wird tangential zur gemessenen Gierung des Schleppkabels am Schiff konstruiert. Die zweite Linie wird parallel zur Korrigierten Kompaßanzeige konstruiert, und sie schneidet die erste Linie an einem Punkt, der mittels einer empirischen Beziehung abgeschätzt ist.

Weitere Modelle, die in manchen Fällen genauere Schätzwerte der Lage des Schleppkabels ergeben, werden später noch erläutert.

Die Erfindung wird nun an Hand der Zeichnung beispielshalber erläutert. Es zeigen:

Fig.1 eine Draufsicht auf ein Schiff, das ein seismisches Seeschleppkabel nachzieht,

Fig.2 die Gierungs- und Längsneigungswinkel des Schleppkabels,

Fig.3 eine Befestigungsvorrichtung zum "terbinden des Schleppkabels mit dem Schiff und zum Messen der Gierungs- und Längsneigungswinkel,

609 8 8 3/0973

Fig.4 die Verbindung des Schleppkabels mit der Befestigungsvorrichtung und die Anbringung eines Zugmessers im Schleppkabely

Fig.5 ein Blockschaltbild einer Anordnung zur Bestimmung der Lage eines Schleppkabels,

Fig.6 ein Flußdiagramm für die Erzeugung eines Positionsmodells in der bevorzugten Ausführungsform _s

Fig.7 das Positionsmodell der bevorzugten Ausführungsform und

Fig.8 ein Flußdiagramm für die Erzeugung einer anderen Ausführungsform eines Positionsmodells.

Fig.1 zeigt schematisch ein "Explorationsschiff, das ein seismisches Seeschleppkabel nachzieht; die Darstellung ist dabei so ausgeführt, wie sie sich bei einer Ansicht der Meeresoberfläche direkt von oben ergibt. Das Schiff 10 bewegt sich längs einer Bahn, die durch die gestrichelte Linie 12 dargestellt ist. In Fig. 1 verläuft die Längsachse des Schiffs 10 zwar in der gleichen Richtung wie die Bahn 12, doch ist es im allgemeinen nicht der Fall, da das Schiff gewöhnlich einen Vorhaitwinkel bezüglich der Bahn einnimmt, längs der es sich bewegt. Die für das seismische Schleppkabel 14 angegebene Position kann als den Ort der Projektion des Schleppkabels auf die Meereeoberfläche angesehen werden. In regelraässigen Abständen längs der Bogenlänge des Schleppkabels sind mehrere Tiefenkontrollkörper 16 der oben erörterten Art dargestellt. Seismische Schleppkabel sind zwar allgemein so konstruiert, daß sie angenähert neutral schwimmfähig sind, doch gewährleisten die Tiefenkontrollkörper 16, daß die Abschnitte des Kabels,

609883/0973

an denen sie befestigt sind, in einer konstanten vorbestimmten Tiefe unterhalb der Meeresoberfläche gehalten werden. Es ist somit zu erkennen ₅ daß ein Abschnitt des Schleppkabels dicht beim Schiff 10 an einem über der Meeresoberfläche liegenden Punkt am Schiff befestigt istj, während der Abschnitt des Schleppkabels am ersten Tiefensteuerkörper 16 in der oben erwähnten vorbestimmten Tiefe gehalten wird. Am Ende des Schleppkabels 14«, d_oh_o an dem vom Schiff 10 am weitesten entfernt liegenden Ende_s ist am Schleppkabel eine Endboje 18 befestigt, die dazu dient_P das Ende des Schleppkabels 14 an der Meeresoberfläche zu halten_o

Die in Fig.1 angegebene Linie 20 ist in der Horizontalebene eines auf dfe Erde bez ogenen Koordinatensystems konstruiert, und sie verläuft tangential zur Projektion des Seekabels 10 auf die Meeresoberfläche. Der Tangentenpunkt ist der Punkt, an dem das Seekabel 14 am Schiff 10 befestigt ist. Die Linie 20 definiert mit der Längsachse des Schiffs einen Winkel Q-', der als der Gierungswinkel des Kabels am Schiff 10 bezeichnet wird. Bei der in Fig.1 dargestellten Lage des Schiffs 10 verläuft die Längsachse des Schiffs längs der Bahn 12.

Der Gierungswinkel θ und der Längsneigungswinkel 0 sind in Fig.2 dargestellt, in der auch ein Einheitsvektor 22 angegeben ist. Der Einheitsvektor 22 verläuft an dem Punkt, an dem das Schleppkabel am Schiff 10 befestigt ist, in der gleichen Richtung wie das Schleppkabel. Dieser Punkt ist durch den Ursprung des Koordinatensystems in Fig.2 gegeben. In dem in Fig.2 dargestellten, auf die Erde bezogenen Koordinatensystem fällt die x-Achse mit der Längsachse des Schiffs 10 zusammen, wobei für den Augenblick die Längsachse des Schiffs 10 so angenommen wird, als falle sie mit der Bahn 12 des Schiffs zusammen. Die

6 0 9883/0973

z-Achse, die ebenfalls in der Horizontalebene liegt, ist die Querachse des Schiffs während die y-Achse die vertikale Hochachse ist. Die Projektion des Einheitsvektors 22 auf die Horizontalebene ist durch den gestrichelten Vektor 24 dargestellt, während die Projektion des Einheitsvektors 22 auf die x~Achse durch den Abschnitt der x-Achse angegeben ist, der am Punkt 26 endet. Aus Fig.2 ist zu erkennen, daß der Gierungswinkel θ der Winkel zwischen der horizontalen Projektion des Schleppkabels am Schiff und der Schiffslängsachse ist. Der Längsneigungswinkel 0 ist der Winkel zwischen dem Einheitsvektor, der das Schleppkabel am Verbindungspunkt mit dem Schiff repräsentiert und der horizontalen Projektion dieses Einheitsvektors. -..'

Eine Vorrichtung 30 zum Messen de's Gierungswinkels θ .und des Längsneigungswinkels 0 am Heck des Schiffs ist in Fig.3 dargestellt. Die Vorrichtung enthält einen Träger 32, der mit Hilfe von Befestigungsbolzen 34· am Heck des Fahrzeugs 10 fest mit dem Deck verbunden ist. Mit Hilfe eines Kreuzgelenks 38 ist eine Halteplatte 36 am Träger befestigt. Ein Kabel 40 ist mit Hilfe einer Kabelrolle mit der;Halteplatte 36 verbunden. Wie aus Fig_e4a zu erkennen ist, ist das Kabel 40 an seinem anderen Ende an einem Wechselteil befestigt, das seinerseits das beim Schiff 10 liegende Ende des seismischen Schleppkabels 14 trägt. Da das Kreuzgelenk 38 zwei Freiheitsgrade zwischen der Halteplatte 36 und dem Träger 32 zuläßt, ergibt sich aus der vorhergehenden Beschreibung daß die Stirnfläche 44 der Halteplatte 36 stets senkrecht zum Einheitsvektor gehalten wird, der den Verlauf des an das Schiff 10 angrenzenden Abschnitts des seismischen Schleppkabels 14 repräsentiert. Es ist somit zu erkennen, daß die Ausrichtung der Halteplatte 36 bezüglich des Trägers 32 den Verlauf des seismischen

609883/0973

Schleppkabels 14 "bezüglich des Schiffs 10 direkt repräsentiert.

Eine Drehung der Halteplatte 36 bezüglich des Trägers 32 um eine parallel zur vertikalen Achse des Schiffs 10 verlaufende Achse wird von einem Drehmelder 46 festgestellt. Diese Winkeldrehung wird vom Kreuzgelenk 38 mit Hilfe einer biegsamen Welle 48 zum Drehmelder 46 übertragen. Das vom Drehmelder- 46 am elektrischen Kabel 50 gelieferte elektrische Signal repräsentiert daher den Gierungswinkel θ des Schleppkabels 14 bezüglich des Koordinatensystems des Schiffs.

In ähnlicher Weise wird eine- Drehung der Halteplatte 36 bezüglich des Trägers 32 um eine Querachse mit Hilfe einer zweiten (nicht dargestellten)biegsamen Welle zum Drehmelder 52 übertragen, der am elektrischen Kabel 54 ein Signal liefert, das den Längsneigungswinkel 0 des seismischen Schleppkabels 14 repräsentiert. Auch hier repräsentiert das Signal am Kabel 54 den Längsneigungswinkel 0 gemessen in einem auf das Schiff bezogenen Koordinatensystem. Die Drehmelder 46 und 52 können imHandel erhältliche Drehmelder der Firma Sperr/Marine Systems mit der Modellnummer 1802 sein. ' .

In der bevorzugten Ausführungsform ist erwünscht, den Gierungswinkel θ und den Längsneigungswinkel 0 in einem auf die Erde bezogenen Koordinatensystem auszudrücken«, Demnach ist am Träger 32 eine Pendeleinheit 56 starr befestigt. Die Pendeleinheit 56 enthält ein erstes Pendelpotentiometer, das in einer Dimension festgehalten ist, so daß es sich nur um die Längsachse des Schiffs 10 frei drehen kann. Es ist zu erkennen, daß dieses Potentiometer für Rollbewegungen des Schiffs 10 empfindlich ist.

609883/0 973

Die Pendeleinheit 56 enthält ein zweites Pendelpotentiometer «, das ebenfalls in einer Dimension festgehalten ist, so daß es sich nur um die Querachse des Schiffs 10 frei drehen kann. Dieses zweite Potentiometer ist somit für Stampfbewegungen des Schiffs 10 empfindlich. Die von diesen zwei Potentiometern gelieferten zwei elektrischen Signale sind an einem elektrischen Anschluß 58 zur weiteren Verarbeitung verfügbar. Eine Vorrichtung, die in der Pendeleinheit 56 als Pendelpotentiometer zufriedesteilend benutzt werden kann, ist von der Firma Humphrey Inc. San Diego,Kalifornien, als Modellnummer CP17-0601-1 im Handel erhältlich.

Aus Fig.2 ist zu erkennen, daß die x-, y~ und z-Komponenten des Einheitsvektors 22 die in der nachfolgend angegebenen Gleichung (1) angegebenen Werte haben, in der θ der Gierungswinkel des Kabels und 0 der Längsneigungswinkel des Kabels in einem auf die Erde bezogenen.Koordinatensystem sind. Die an den Kabeln 50 und 54 erscheinenden elektrischen Signale repräsentieren jedoch die Gierungsund Längsneigungswinkel des Kabels bezüglich des Koordinatensystems des Schiffs, dessen Längs-,Vertikal- und Querkoordinatenachsen mit 1, ν bzw. t bezeichnet sind. Die 1-, v-und t-Komponenten des Einheitsvektors sind durch die nachfolgend angegebene Gleichung (2) gegeben, in der θ und 0 die Gierungs- bzw. Langsneigungswinkel des Kabels bezüglich des Koordinatensystems des Schiffs sind; diese ¥inkel werden dabei durch die an den Kabeln 50 und 54 erscheinenden elektrischen Signale repräsentiert.

(D

cos φ cos sin 0_s
cos 0 sin

(2)

609883/0973

Die Einheitsvektoren U₂₂ und u stehen über die Trans-.formationsraatrix T miteinander in Beziehung wie die
Gleichung (3) zeigt.

U_s (3)

Die einzelnen Elemente der Transformationsmatrix/"^i?/sind in der Gleichung (4) angegeben.

cosP -sinP .0

sinP cosR cosP cosR 'ginR cosP

-sinR sinP cosP -sinR cos P cosR.

D 3D

(4)

wobei gilt: D= (1 - sin²R sin²P) ¹/² .

In der Gleichung^) sind die Größen R und P die Querbzw. Längsneigungswinkel des Schiffs, die von den Signalen repräsentiert werden, die am elektrischen Anschluß 58 erscheinen. Die an den elektrischen Anschlüssen 50, 54 und 58 erscheinenden Signale werden gemäß der Gleichung (3) so kombiniert, daß sich die x-, y- und z-Komponenten des Einheitsvektors 22 ergeben, wie in der Gleichung (1) angegeben ist«

Die einzelnen Komponenten des Kabelvektors ΰ*₂₂, die in der Gleichung (1) gegeben sind, werden so kombiniert, wie in der Gleichung (5) gezeigt ist, damit die Gierungsund Längsneigungswinkel'0 bzw. θ des Kabels erhalten
werden.

θ = tan"¹(u_z/u_x) ; 0=sin¹(Uy) (5)

603883/097 3

2831634

dabei sind: u die x-Komponente von U

u die y-Komponente von
u die z-Komponente von

Der Gierungswinkel θ ist bis hierher relativ zur Längsachse des Schiffs 10 gemessen. Der Gierungswinkel Q wird weiterhin, mit dem Steuerkurswinkel des Schiffs 10 kombiniert der von einem einen Teil des Navigationssystems des Schiffs bildenden Kreiselkompaß geliefert wird_? damit der Gierungswinkel θ bezüglich der Nordrichtung ausgedrückt wird» Es sei angenommen^ daß von nun an der berechnete Gierungswinkel θ so korrigiert worden ist, daß er-auf die Nordrichtung bezogen ist, _

Die Befestigungsvorrichtung 30 arbeitet mit der in Pig_o4a dargestellten Haltevorrichtung so zusammen,, daß das Schleppkabel 14 mit dem Heck des Schiffs 10 verbunden wird. Nach Fig»4a verbindet das Kabel 40 die Kabelrolle 42 mit dem ¥echselteil 7O₀ Das Schleppkabel"14 ist am Punkt 72 starr am Wechselteil 70 befestigt. Am Punkt 72,, an dem das Seekabel 14 mit dem Wechselteil 70 zusammenwirkt₉ tritt aus dem Schleppkabel 40 ein Kabelbündel 74 aus_o Bei einer typischen Schleppkabelanordnung enthält das Kabelbündel 74 für jede längs des Schleppkabels 14 angebrachte Hydrophongruppe zwei elektrische Leiter„ Im vorliegenden Ausführungsbeispiel enthält das Kabelbündel 74 außerdem elektrische Leiter,, die zu Zugmessern und zu einem Magnetkompaß führen, die ebenfalls an Stellen längs des Schleppkabels 14 angebracht sind_o Das Kabelbündel 74 ist in bekannter Weise an eine an Bord des Schiffs 10 befindliche Signalaufzeichnungs- und Signalverarbeitungsanlage angeschlossen«,

609883/097 3

Im bevorzugten Ausführungsbeispiel ist an einer dicht beim Wechselteil 70 liegenden Stelle 74 im Schleppkabel 14 ein erster Zugmesser angebracht. Fig.4b zeigt eine Explosionsansicht des Inneren des Schleppkabels im Bereich des ersten Zugmessers. Ein geeigneter Zugmesser ist ein Wandler des Modells Nr. WML2-251-10K 3926, der im Handel von der Firma Transducers, Inc. Whittier, Kalifornien erhältlich ist. Nach Fig.4b hat der Zugmesser 76 einen allgemein sechseckigen Querschnitt. An jedem Ende des Zugmessers 76 befindet sich eine Gewindebohrung 78, mit deren Hilfe der Zugmesser 76 mit dem Schleppkabel verbunden ist. Zum Herstellen der externen elektrischen Verbindung mit der Wheatstone-Brücke des Zugmessers 76 sind vier Anschlußstifte 80 vorgesehen. Zum Verbinden eines Endes des Schleppkabels 14 mit demZugmesser 76 ist eine Halteplatte 82 vorgesehen. Die Halteplatte 82 enthält einen Gewindezapfen 84, der mit einer der Gewindebohrungen so zusammenwirkt, daß die Halteplatte 82 starr am Zugmesser 76 befestigt wird. Die drei aus Edelstahl bestehenden Zugglieder 86, die die Zugaufnahmeglieder in einem herkömmlichen seismischen Schleppkabel bilden, sind an den-Punkten 88 an der Halteplatte 82 befestigt. Teile der Halteplatte 82 sind wie bei 90 weggeschnitten, damit das Kabelbündel des Schleppkabels durch die Halteplatte 82, seitlich längs des Zugmessers 76 und weiter längs des Schleppkabels 14 zum Schiff 10 geführt werden kann. Weitere elektrische Leitungen, die an die Anschlußstifte 80 angeschlossen sind, werden ebenfalls ein Teil des zum Schiff 10 führenden Kabelbündels. Eine der Halteplatte 82 ähnliche zweite.Halteplatte ist zum Verbinden des vorderen Endes des Schleppkabels 14 mit dem vorderen Ende des Zugmessers 76 vorgesehen. Das gesamte seismische Schleppkabel einschließlich des den Zugmesser 76 enthaltenden Abschnitts ist in bekannter Weise von einem Mantel umgeben, der typischer-

609883/0973

~ 14 -

weise aus einem Material wie Polyvinylchlorid extrudiert ist.

Bei seismischen Seeschleppkabeln ist es allgemein üblich, in dem zwischen dem Schleppschiff und dem ersten Tiefensteuerkörper befindlichen Abschnitt des Schleppkabels einen Dehnungsabschnitt einzufügen. Der Dehnungsabschnitt unterscheidet sich vom herkömmlichen Schleppkabelabschnitt dadurch, daß die aus Edelstahl bestehenden Zugglieder durch Zugglieder aus Nylon ersetzt werden, die relativ elastisch sind_o Die Anwesenheit des Dehnungsabschnitts setzt die Wahrscheinlichkeit einer Zerstörung des Schleppkabel^ durch Stoßbeanspruchung herab, und er ergibt eine Isolierung der Hydrophone von mechanischen Geräuschen, die vom Schleppschiff hervorgerufen werden. Damit der genaue Bogenabstand zwischen jedem Punkt des Schleppkabels außerhalb des Dehnungsabschnitts und dem Schiff bekannt.ist, muß jedoch auch die prozentuale Dehnung des Dehnungsabschnitts bekannt sein.

Der Zugmesser 76 liefert ein Signal, das die in Längsrichtung auf das Schleppkabel an der Stelle einwirkende Zugbeanspruchung repräsentiert, an der das Schleppkabel mit.dem Heck des Schiffs 10 verbunden ist«, In der bevorzugten Ausführungsform-ist dicht bei dem dem Schiff 10 am nächsten liegenden Tiefensteuerkörper 16 ein zweiter ebensolcher Zugmesser im Schleppkabel 14 angebracht. Wie oben genauer erläutert wurde, können die von diesen zwei Zugmessern gelieferten Zugmeßergebnisse dazu verwendet werden, das Ausmaß der elastischen Dehnung des zwischen den zwei Zugmessern befindlichen Abschnitts des Schleppkabels 14 zu bestimmen.

Im Schleppkabel 14 ist ein Magnetkompaß an einer Stelle angebracht, die soweit vom Schiff 10 entfernt ist, daß sichergestellt ist, daß Verzerrungen des Brdmagnetfelds

609883/0973

durch das Feld des Schiffs 10 keine unannehmbaren Fehler der festgestellten magnetischen Nordrichtung herbeiführen. In der bevorzugten Ausführungsform befindet sich der Magnetkompaß beispielsweise an einer etwa in der Mitte der Länge des Schleppkabels 14 befindlichen Stelle. Ein geeigneter Magnetkompaß ist der von der Firma Digicourse Inc. erhältliche Steuerkursfühler Modell 314„ Dieser spezielle: Magnetkompaß hat die erwünschte Eigenschaft, daß er kardanisch aufgehängt ist, so daß er um volle 360° um die Längsachse des Schleppkabels 14 drehbar ist. Als Folge davon hat die längs des Seekabels 14 üblicherweise auftretende Verdrillung keine störende Auswirkung auf das richtige Arbeiten des Magnetkompasses. Zusätzliche elektrische Leiter im Kabelbündel übertragen das vom Magnetkompaß gelieferte elektrische Signal zu der an Bord des Schiffs 10 befindlichen Verarbeitungsanlageο '

Die gegenseitige Verbindung der Baueinheiten, die die bevorzugte Ausführungsform der Erfindung bilden_s ist in dem Blockschaltbild von Fig.5 dargestellt. Die zum Abtasten verschiedener Parameter in der Anordnung verwendeten Baueinheiten enthalten die Zugmesser 100 und 102, den Magnetkompaß 104, die Drehmelder 46 und 52 und die zwei Pendelpotentiometer in der Pendeleinheit 56. Die von den Zugmessern 100 und 102 gelieferten Analogsignale werden in Analog-Digital-Umsetzern 106 bzw«, 108 digitalisiert. In gleicher Weise werden auch die von den zwei Potentiometern der Pendeleinheit 56 gelieferten Analogsignale mittels der Analog-Digital-Umsetzer 116 und -118 digitalisiert, Die Analog-Digital-Umsetzer können beliebige im Handel erhältliche Umsetzer sein; ein dafür geeigneter Umsetzer ist beispielsweise das von der Firma Analog Devices vertriebene Modell Nr. AD2003 mit Digitalausgang.

609883/097 3

Der Magnetkompaß 104 ist an die Schnittstelleneinheit 110 angeschlossen«, In der bevorzugten Ausführungsform ist die Schnittstelleneinheit 110 die von der Firma Digicourse Inc. hergestellte Steuerkursfühler-Schnittstelleneinheit Modell 251. Die Schnittstelleneinheit 110 ist auch mit einem Block 120 verbunden, in dem zur Erzielung einer digitalen Anzeige der Kompaßablesung die vom Magnetkompaß 104 und von der Schnittstelleneinheit 110 gelieferte Impulsfolge gezählt wird. Die Einheit 120 dient auch dazu in ausgewählter Weise unter der Steuerung durch den Hauptprozessor 126 an den Magnetkompaß 104 und an die Schnittstelle ' 110 Versorgungsenergie anzulegen. Beim bevorzugten Betriebsverfahren wird die Versorgungsenergie des Magnetkompaß in den Zeitperioden abgeschaltet, in denen seismische Daten erfaßt werden, damit die' Möglichkeit einer Störverfälschung der seismischen Daten eliminiert wird. Der Magnetkompaß wird so eingeschaltet, daß unmittelbar vor der Auslösung einer seismischen Störung und der Aufzeichnung der resultierenden seismischen Daten ein Kompaßablesewert erhalten wird.

Die von den Drehmeldern 46 und 52 gelieferten Analogsignale werden in Analog-Digital-Wandler 112 bzw_e 114 digitaliaiert. Die Analog-Digital-Wandler 112 und 114 sind Synchro-Digital-Wandler, Modell M5OOOM36/3DP, die von der Firma AstroSystems, Inc. Lake Success, New York hergestellt und vertrieben werden.

Die gemäß^;den obigen Ausführungen erzeugten digitalen Signale werden im Hauptprozessor 126 so kombiniert, daß Modelle der Schleppkabelposition erhalten werden«, Der Hauptprozessor 126 kann ein von der Firma Texas Instruments Incorporated, Dallas, Texas hergestellter und vertriebener Computer mit der Modellbezeichnung 980A sei.n. Wie durch den Block 122

609883/0973

angegeben ist, erfolgt die Übermittlung von Datenwerten zum Hauptprozessor 126 über eine Eingabe/Ausgabe-Hauptprozessor-Schnitt stelle 122_e Das Vorsehen der Eingabe/Ausgabe-Schnittstelle zum Abfragen der verschiedenen digitalen Abtastwerte und zum Übermitteln der digitalen Abtastwerte zum Hauptprozessor 126 ist für den Fachmann ohne weiteres verständlich ₉ so daß hier keine weiteren Ausführungen erforderlich sind. Als Beispiel ist jedoch ein Abschnitt der Eingabe/Ausgabe-Schnittstelle 122 in Fig.5 ausdrücklich als Eingabe/Ausgabe-Datenmodul 124 bezeichnet worden. Dieser Abschnitt der Schnittstelle, der für die Zugmesser- und Magnetkompaß-Signale zuständig ist_? besteht aus einem Eingabe/Ausgabe-Datenmodul, Modell 961648-2 der Firma Texas Instruments Incorporated. Die Eingabe/Ausgabe-Hauptprozessor-Schnittstelle 122 erfaßt periodisch, d.h. einmal pro Sekunde, einen Datenabtastwert von jeder der Eingabeeinheiten, und sie gibt diese Datenabtastwerte zum Hauptprozessor 126 weiter, wo sie dazu verwendet werden, ein Modell der Schleppkabelposition zu erzeugen. Die verschiedenen Datenabtastwerte und das Modell werden in der Magnetbandeinheit 128 aufgezeichnet. Als Magnetbandeinheit 128 kann das Modell TI-979 der Firma Texas Instruments Incorporated verwendet werden.

Die vom Hauptprozessor 126 in der bevorzugten Äusführungsform der Erfindung ausgeführten Operationen sind in dem Flußdiagramm von Fig.6 veranschaulicht. Die Verarbeitung, die zur Erzielung des Modells der Schleppkabelposition für jede Gruppe von Abtastwerten aus den verschiedenen Eingabeeinheiten angewendet wird, wird durch die Eintrittsanweisung 140 ausgelöst. Beim Schritt 142 erfaßt die Anordnung die Momentanwerte des Breitengrads, des Längengrads und des Azimutwinkels des Schiffs 10. Im bevorzugten Ausführungs-

609883/0973

- .18 -

beispiel arbeitet die Anordnung zur Bestimmung der Schleppkabelposition zusammen mit dem bereits erwähnten, für den Einsatz bei der seismischen Exploration auf See bestimmten Navigations- und Ortungssystem GEONAV. Insbesondere wird der Hauptprozessor 26 von der Anordnung zur Bestimmung der Schleppkabelposition und von dem NavigaUons-und Ortungssystem gemeinsam benutzt. Als Folge davon enthält der Hauptprozessor 26 jederzeit die neuesten Meßwerte über die geor graphische Länge, die geografische Breite und den Azimutwinkel des Schiffs, so daß diese Vierte wie beim Schritt stets für die Benutzung durch die Schleppkabelpositionsanordnung verfügbar sind. Wie später noch genauer erläutert wird, wird der Azimutwinkel des Schiffs dazu verwendet, die Berechnung des Schleppkabelmodells in einem auf die Erde bezogenen Koordinatensystem -zu ermöglichen. Die Erfindung ist jedoch nicht darauf beschränkt; falls der Azimutwinkel des Schiffs nicht verfügbar ist, kann sie zur Erzielung eines Schleppkabelmodells verwendet werden, das in einem auf das Schiff bezogenen Koordinatensystem ausgedrückt ist.

Beim Schritt 144 erfaßt der Hauptprozessor 126 den neuesten Datenabtastwert vom Magnetkompaß 144, der an einer Stelle längs des Schleppkabels 14 angebracht ist. Im Schritt 146 wird der unverarbeitete Kompaßablesewert zur Beseitigung von Ungenauigkeiten des Kompaß korrigiert. Der Kompaß 104 muß vor seiner Anbringung im Schleppkabel 14 getestet worden sein, damit eine empirische Kurve erhalten wird, die die wahre magnetische Ablesung abhängig von der unverarbeiteten Kompaßablesung ausdrückt. Diese empirische Kurve wird zuvor in den Hauptprozessor 126 in Form einer aus acht linearen Segmenten bestehenden Näherung eingegeben. Beim Schritt 146 wendet der Hauptprozessor 126 die unverarbeitete Kompaßablesung auf diese Kurve an, damit die wahre magnetische

609883/0973

Kompaßablesung erhalten wird. Beim Schritt 148 wird die magnetische Abweichung aus der Kompaßablesung entfernt, damit ein auf die wahre Nordrichtung bezogener Steuerkurs erhalten wird. Die magnetische Abweichung wird von einer Bedienungsperson in den Hauptprozessor 126 eingegeben; sie kann periodisch aktualisiert werden, so daß die vom Hauptprozessor verwendete magnetische Abweichung für das Gebiet paßt, in dem die Exploration ausgeführt wird.

Beim Schritt 149 wird der Azimutwinkel des Schiffs vom eingestellten Kompaßsteuerkurs subtrahiert, damit der Winkel Ö2erhalten wird, der die Richtung des Schleppkabels bezüglich des Koordinatensystems des Schiffs ausdrückt.

Beim Schritt 150 erfaßt der Prozessor von der Befestigungsvorrichtung 20 die aktuellsten Abtastwerte des Gierungswinkels und des Längsneigungswinkels des Schleppkabels sowie des Querneigungswinkels und des Längsneigungswinkels des Schiffs. Beim Schritt 152 werden die Sinuswerte und die Cosinuswerte dieser Winkel bestimmt, wobei diese trigonometrischen Funktionen gemäß den Gleichungen (1) bis (4) dazu verwendet werden, einen Schleppkabel-Einheitsvektor ausgedrückt in auf das Schiff bezogenen Koordinaten zu erhalten. Beim Schritt 154 wird die Größe D berechnet, die in der Gleichung (4) definiert ist. Beim Schritt 156 bestimmt der Prozessor den Schleppkabel-Einheitsvektor in einem auf das Schiff bezogenen Koordinatensystem, wie in der Gleichung (2) definiert ist. Beim Schritt 158 wird die durch die Gleichung (3) ausgedrückte Drehung ausgeführt. Im Schritt 160 werden dann entsprechend der Gleichung (5) die in einem erdbezogenen Koordinatensystem ausgedrückten Komponenten des Schleppkabel-Einheitsvektors dazu verwendet, den Gierungswinkel Θ1 und den Längsneigungswinkel 0 des Schleppkabels zu bestimmen.

60988 3/097 3

Beim Schritt 162 erfaßt der Hauptprozessor 126 die aktuellsten Abtastwerte von den Zugmessern 1 und 2. Beim Schritt 164 wird die Differenz zwischen den vom Zugmesser 1 und den vom Zugmesser 2 gemessenen Zugwerten dazu bestimmtj die prozentuale Dehnung des Dehnungsabschnitts des Kabels zu bestimmen. Auf empirische Weise ist bestimmt worden, daß die prozentuale Dehnung durch eine quadratische Funktion angenähert werden kann, wenn sie als Funktion der Zugdifferenz an den zwei Enden des Dehnungsabschnitts ausgedrückt wird. Die entsprechende quadratische Funktion wird zuvor in den Hauptprozessor 126 eingegeben und beim Schritt 164 zusammen mit der Differenz zwischen den zwei gemessenen Zugwerten zur Erzielung der prozentualen Dehnung verwendete Beim Schritt 166 wird diese prozentuale Dehnung mit der im Gleichgewichtszustand vorliegenden Länge des Dehnungsabschnitts multipliziert, damit die Dehnungslänge des -Schleppkabels erhalten wird. Beim Schritt 168 wird diese Dehnungslänge zu den Bogenabständen des Vorderendes und des Hinterendes des Schleppkabels addiert, d.h. den Längen des unbelasteten Schleppkabels vom Heck des Fahrzeugs zum Anfang bzw. zum Ende des aktiven Abschnitts des SchleppkabeJL's 14. In diesem Zusammenhang ist das Torderende des Schleppkabels diejenige Stelle, an der sich die am dichtesten beim Schiff 10 angebrachten Hydrophone befinden.

Die mit Hilfe von zwei"Linien durchgeführte Baherung, die im bevorzugten Ausführungsbeispiel der Erfindung das Modell der Schleppkabelposition bildet, ist in Fig.7 dargestellt. Die erste Linie 172 des Modells verläuft durch das Heck des Schiffs 10;.sie verläuft in einem Winkel ©1 zur x-Achse oder Längsachse des Schiffs 10. Die zweite Linie 174- verläuft in einem Winkel ©2 zur x-Achse. Der Punkt, an dem sich die Linien 172 und schneiden, hat die Koordinaten (XB, ZB) die beim Sehritt

609 88 3/0973

von Jig.6 durch die in der Gleichung (6) angegebene empirische Beziehung abgeschätzt werden-

XB = K(tan 02 - 1) / (tan Θ1 - tan Θ2);

ZB = XB(tan 01) (6)

4
wobei gilt: K = 2,586 (01/9,08)

Beim Schritt 176 werden die Polarkoordinaten des Anfangsund des Endes des Schleppkabels in einem auf die Nordrichtung bezogenen System unter Anwendung ihrer Bogenabstände, der Winkel 01 und 02 der Koordinaten XB und ZB sowie des Azimutwinkels des Schiffs berechnet„ Auch der Ort desKnickpunkts (XB, ZB) wird transformiert, damit dieser Knickpunkt im Polarkoordinatensystem ausgedrückt wird. Beim Schritt 178 werden die unverarbeiteten Daten, die Orte des Anfangs und des Endes des Kabels_s der Ort des Knickpunkts in Polarkoordinaten, sowie die geographische Breite, die geographische Länge und der Äzimutwinkel des Schiffs zur Magnetbandeinheit'128 ausgegeben.

■Beim Schritt 180 wird der nächste Abtastwert vom Magnetkompaß 104 zur Vorbereitung der nächsten Berechnung abgelesen. Nachdem dieser nächste Abtastwert vom Kompaß erfaßt worden ist, wird beim Schritt 182 der Befehl "Kompaßleerlauf" zum Datenmodul gesendet. Wie bereits erläutert wurde, versetzt dieser Befehl den Magnetkompaß während der Erfassung seismischer Daten in einen abgeschalteten Zustand. Am Austrittspunkt 184 ist die Erzeugung des Modells beendet. Die an Hand des Flußdiagramms von Pig.6 veranschaulichte Operationsfolge wird in periodischen Intervallen, typischerweise einmal pro Sekunde wiederholt. Auf diese Weise werden auf dem Magnetband die unverarbeiteten

603883/0973

Daten und ein periodisch aktualisiertes Modell der momentanen Schleppkabelposition aufgezeichnet«

Es werden nun ein Alternativverfahren und eineAlternativanordnung zur Bildung eines Modells der Schleppkabelposition unter Verwendung der Ausgangssignale der verschiedenen Fühler des in Fig.5 dargestellten Systems erörtert. In der erwähnten Druckschrift von John Bedenbender sind sechs Grleichgewichtsdifferentialgleichungen für ein seismisches Schleppkabel angegeben* Diese Gleichungen sind hier als Gleichungen(7) bis (12) wiedergegeben.

H + F - ¥ sin 0 = 0 (7)

Ϊ^;cos 0 H + H = 0 (8)

⁵ "TZ +G-W cos 0=0 (9)

= cos 0 cos θ (10)

= sin 0 - . (11)

= cos 0 sin θ (12)

In diesen Gleichungen sind x,_v y und ζ rechtwinklige Koordinaten,

T die Kabel spannung s die Bogenlänge längs des Schleppkabels ¥ das Gewicht' des Schleppkabels pro Längeneinheit im flüssigen Medium,

θ der Winkel zwischen der positiven x-Achse und der das Element ds enthaltenden Tertikaiebene,

609883/0973

0 der Winkel zwischen ds und der x-z-Ebene in der . Vertikalebene des Elements ds,

ί, H und G hydrodynamische Kräfte pro Längeneinheit des Schleppkabels, die auf ds jeweils tangential zu ds, senkrecht zu ds sowie senkrecht zur Vertikalebene von 0 und senkrecht zu ds in der Vertikalebene von 0 wirken.

In diesen Gleichungen ist die Bogenlänge s längs des Schleppkabels die unabhängige Variable, während x, y, z, T, 0 und O die abhängigen YariaBlejL sind. In vielen lallen, so auch in der vorliegenden Anordnung zur Bestimmung der Schleppkabelposition, sind die Anfangswerte einiger abhängiger Variabler und die hydrodynamischen Kräfte nicht so ausreichend bekannt, daß eine bestimmte numerische Integration der Differentialgleichungen möglich wäre. Die erwähnte Druckschrift von Bedenbender gibt ein Lösungsverfahren an, das auf einige solcher Probleme anwendbar ist. Bei diesem Verfahren werden die unbekannten Anfangswerte aus der allgemeinen -Kenntnis der physikalischen Situation abgeschätzt. Unter Anwendung dieser Anfangswerte werden die Differentialgleichungen längs des Schleppkabels bis zu einer Stelle integriert, an der bestimmte abhängige Variable bekannt sind. Diese bekannten Werte der abhängigen Variablen (G-renzbedingungen) werden mit den aus der Integration erhaltenen Werten der Differentialgleichungen verglichen, und die Unterschiede "werden in einem Ite-:. rationsverfahren nach Newton-Raphson dazu verwendet, neue Schätzwerte für die unbekannten Anfangswerte zu erhalten. Dieses Verfahren wird iterativ unter Verwendung der neuen AnfangsSchätzwerte durchgeführt, bis die durch die integrierten Differentialgleichungen erhaltenen Werte innerhalb vorgewählter Abstände von den

609883/0973

bekannten Grenzbedingungen liegen.

Das erwähnte iterative.Lösungsverfahren ist für die Anwendung bei Eingangsvariablen erweitert worden, die in der Anordnung zur Bestimmung der Position eines seismischen Schleppkabels verfügbar sind. Die Integration der Differentialgleichungen beginnt an einem Punkt längs des Schleppkabels, beispielsweise an dem am dichtesten beim Schiff 10 liegenden Punkt, an dem sich einer der liefensteuerkörper 16 nach ]?ig.1 befindet. Da der Tiefensteuerkörper das Schleppkabel an diesem Punkt in einer vorbestimmten Tiefe hält, ist der Anfangswert der abhängigen Variablen y bekannt. Den abhängigen Positionsvariablen χ und ζ wird willkürlich der Wert 0 zugeordnet. Basierend auf den Kenntnissen der beim Schleppen eines seismischen Kabels allgemein auftretenden Bedingungen werden Schätzwerte für die abhängigen Variablen T, 0 und ©'bestimmt. Mit diesen Anfangsbedingungen wird dann begonnen,·die Differentialgleichungen numerisch längs des Sohleppkabels zum Schiff 10 hin zu integrieren, damit Lösungswerte für die abhängigen Variablen am Heck des Schiffs erhalten werden. Die Lösungswerte für die abhängigen •Variablen v_s T und 0 werden mit den bekannten Werten J₉ Ψ und "0 am Heck des Schiffs 10 verglichen; dabei sind f die bekannte Höhe des Schleppkabels am Schiff 10, ¥ die mit Hilfe des am Punkt 74 in Pig.4a angebrachten Zugmessers gemessene Kabelspannung und JZF der mit Hilfe der Befestigungsvorrichtung 30 gemessene Längsneigungswinkel des Schleppkabels. Im Iterationsverfahren nach Newton-Raphson wird dieser Vergleich dazu verwendet, neue Schätzwerte für die Anfangswerte von T, 0 und θ am ersten Tiefensteuerkörper zu erhalten. Dieses iterative Verfahren wird solange fortgesetzt, bis die

609833/0 9.7

Lösungswerte für y, T und 0 innerhalb vorgewählter Schwellenwerte der bekanngen Werte konvergieren.

Gewöhnlich sind die hydrodynamischen Kräfte F₉ H und G, die sich aus der Strömungsgeschwindigkeit des Wassers relativ zum Schleppkabel ergeben, nicht genau bekannt. Die Schleppkabelposition wird stark von der quer zur Schlepprichtung auftretenden Wasserströmungsgeschwindig-, keit beeinflußt, von der bekannt ist, daß sie sich manchmal abhängig von der Tiefe ändert. Bei der Durchführung der oben erwähnten Iteration mit ,Anfangsbedingungen am ersten Tiefensteuerkörper wird ein anfänglicher Querströmungsgeschwindigkeitgradient angenommen. Es wird eine äußere Iterationssch-leife nach Newton-Raphson durchgeführt, damit bessere Schätzwerte des Querst römungsgeschwindigkeitsgradientai erhalten werden. Die abhängige Variable, die in dieser äußern Hewton-Raphson-Schleife verweni et wird, ist der G-ierungs winkel θ am Heck des Schiffs -10. Nachdem die inneren und die äußeren Schleifen zur Erzielung befriedigender Schätzwerte sowie des QuerStrömungsgeschwindigkeitsgradienten konvergiert sind, bildet die abschliessende •Integration der Differentialgleichungen unter Verwendung dieser verbesserten Schätzwerte das Modell der Schleppkabelpositioi für den Abschnitt des Schleppkabels, -über den die Integration durchgeführt worden ist.

Unter Anwendung der neuen Schätzwerte der Anfangsbedingungen am ersten Tiefensteuerkörper als bekannte Werte für das nächste Segment des Schleppkabels, d.h. das zwischen dem ersten und dem zweiten Tiefensteuerkörper liegende Segment,kann das Verfahren zur Erzielung verbesserter Schätzwerte der Anfangsbedingungen am Ort

609883/0973

des zweiten Tiefensteuerkörpers wiederholt werden. Auf diese Weise kann die Lösung schrittweise längs des Schleppkabels angewendet werden, so daß schließlich ein Modell der Schleppkabelposition über die gesamte länge erhalten wird.

Wie noch erläutert wird, kann als Alternative eine hyperbolische Kurve an die Lösung für den zwischen dem Schiff 10 und dem ersten Tiefensteuerkörper liegenden Abschnitt des Schleppkabels angepaßt werden.

Eine PORTRAU-Liste zur Ausführung der iterativen Lösung ist als Anlage A beigefügt. Die in dem PORTRAl-Programm erscheinenden Variablen s-ind in der Tabelle I definiert. Die Figuren 8a bis 8i bildai ein Plußdiagramm,das dem PORTRAH-Programm entspricht.

In Pig.8a sind beim Schritt 200 verschiedene Eingangsgrößen aufgeführt. Diese Eingangsgrößen enthalten die Eingangsströmungsgeschwindigkeit und die Kabelzugkoeffizienten, die dazu dienen, die in den Differentialgleichungen erscheinenden hydrodynamischen Kräfte zu definieren. Verschiedene Parameter, die das beim Integrationsverfahren gewünschte Auflösungsvermögen definieren, werden ebenfalls festgelegt. Schließlich erfolgt die. Eingabe der Schätzwerte für die Anfangswerte der Kabelspannung, des Gierungswinkels und des Neigungswinkels am ersten Tiefensteuerkörper, sowie der bekannten Grenzwerte am Schiff 10 für die Kabelspannung, den G-ierungswinkel, den Neigungswinkel und die Höhe. Beim S ehritt wird eine erste Abschätzung des ersten Querströmungsgesehwindigkeitsgradienten VS1 durchgeführt, und die Variable YSLOPE wird gleich . VS1 gesetzt. Beim Schritt 202 wird

609883/097 3

ein Zähler IPLAG- für die äußere Iterationsschleife ausgelöst, und der Anfangswert der Höhe Y am ersten Tiefensteuerkörper wird auf die vorbestimmte Tiefe eingestellt, die von diesem Tiefensteuerkörper aufrechterhalten wird. Beim Schritt 203 wird der Zähler HT für die innere Iierationsschleife ausgelöst.Baim Schritt 204 wird ein Zähler JJ auf den Anfangswert 1 eingestellt. Die beim Schritt 205 angegebene bedingte Verzweigung bewirkt eine Verzweigung zum Punkt 505, wenn der Stand des Zählers JJ den Wert 1 hat. Beim Schritt 206 erfolgt die Zuordnung der ersten G-ruppe von Anfangswerten für die Kabelspannung, den Gierungswinkel und den Neigungswinkel am ersten Tiefensteuerkörper.

Unter Yerwendung dieser Anfangswerte werden beim Schritt die Differentialgleichungen vom Ort des ersten Tiefensteuerkörpers zum Heck des Fahrzeugs 10 längs des Schleppkabels numerisch integriert. Hierbei kann eines von verschiedenen Verfahren zur numerischen Integration angewendet werden; ein geeignetes Verfahren ist das Verfahren nach Adams-Moulton mit einem Beginn nach Runge-Kutta. Dieses bestimmte numerische Integrationsverfahren wird mit Hilfe eines Unterprogramms EEAM ausgeführt, für das die PCRTRAtf-Codierung in der Anlage A enthalten ist. " ' - . ·

Die Integration ergibt die Endwerte für die Variablen Y, T und 0, wie beim Schritt 208 angegeben ist. Beim Schritt 209 werden die Unterschiede zwischen diesen berechneten Bndwerten und den entsprechenden Grenzbedingungen bestimmt. Im Plußdiagramm erfolgt dann ein Sprung vom Schritt 210 zum Punkt 510 in Fig.8d, von dem aus beim Schritt 211 der Zähler JJ um den Wert 1 erhöht wird, so daß er dann den Zählerstand 2 hat. Als Folge davon ergibt die beim Schritt 212 gestellte Präge die.Antwort "Nein" und dasElußdJ^ramm kehrt wieder zum Punkt B von Pig.8a zurück.

60 9 8 8 37 0 9 7 3

In diesem Pall erfolgt "beim Schritt 205 im Plußdiagramm eine Verzweigung zum Punkt 506. Als Ergebnis entsprechen die "beim Schritt 216 gebildeten Fehlerausdrüeke einem geringfügig unterschiedlichen Anfangswert der Kabel-.Spannung als es bei den beim Schritt 209 berechneten !ehlerausdrücken der lall ist. Im Plußdiagramm erfolgt dann ein Sprung vom Schritt 217 zum Punkt 510, worauf die Variable JJ wieder erhöht wird und eine erneute Rückkehr zum Punkt B erfolgt, woran sich schließlich eine Verzweigung zum Punkt 507 anschließt. Die beim Schritt 218 eingestellten Anfangswerte unterscheiden sich von den Anfangswerten beim.Schritt 206 insofern, als der Anfangswert für den Gisrungswinkel geringfügig verändert ist, was beim Schritt 221 zu einer dritten Gruppe von I¹ ehlerausdrücken führt. Nachdem der Ablauf wieder zum Punkt 508 zurückgekehrt ist,wird in gleicher Weise der Anfangswert des Neigungswinkels beim Schritt geringfügig geändert, was zu einer vierten Gruppe, von Pehlerausdrücken beim Schritt 226 führt. Vom .Schritt aus erfolgt ein Sprung zum Punkt 510; in diesem Fall wird dar Zähler JJ beim Schritt 211 so erhöht_s daß er den Zählerstand 5 hat. Die beim Schritt 212 durchgeführte •Prüfung hat erneut kein bejahendes Ergebnis, so daß im ilußdiagramm eine Rückkehr zum Punkt B erfolgt, von wo aus beim Schritt 205 eine Verzweigung zum Punkt 509 erfolgt.

Da der Operationsablauf an diesem Punkt die innere Grenzwertschleife zum ersten Mal durchläuft_s hat die beim Schritt 228 durchgeführte Prüfung ein verneinendes Ergebnis, und beim Schritt 229 wird eine fünfte Gruppe von Anfangswerten gebildet. Es ist zu erkennen, daß sich alle Werte der Anfangswertgruppe beim Schritt 229 geringfügig von den ursprünglich beim Schritt 206 eingegebenen Werten unterscheiden. Dies führt beim Schritt

609883/0973

2831634

"zu einer fünften Gruppe von Fehlerausdrücken und im Flußdiagramm erfolgt ein Sprung vom Schritt 233 zum Punkt 510. Der Zähler JJ wird nun so weitergescla ltet, daß beim Schritt 211 sein Zählerstand den Wert 6 hat, und die Antwort auf die "beim Schritt 212 gestellte Präge ist "Ja",so daß das Flußdiagramm nun. mit dem Schritt 234 fortfährt.

Die verschiedenen Fehlerausdrücke XLMAD (Iß),die sich aus den bei den Schritten 206 bis 233 durchgeführten Integrationen ergeben, werden bei dem oben erwähnten iterativen Verfahren nach ETewton-Raphson dazu verwendet, verbesserte AnfangsSchätzwerte TEE3_? THETA3 und PHI3 entsprechend der folgenden Matrixgleichung (13) zu erzeugen:

TEE2 - TEE 3 THETA2 - THETA3 PHI2 - PHI3

A(1,1) A(1,2) A(1,3) A(2,1) A(2,2) A(2,3) A(3,D A(3,2) A(3,3)

-1

XLMDA(I,5) XLMDA(2,5) XLMDA(3,5)

Die Matrixelemente 4. (I»J) werden bei den Schritten und 235 berechnet; es ist zu erkennen, daß es sich dabei um Näherungen der partiellen Ableitungen der Fehlerausdrücke bezüglich der drei Anfangswertgroßen TEE, THETA und PHI handelt. Die Matrixgleichung (13) wird bei den Schritten 234 bis 240 des Flußdiagramms ausgeführt, damit neue Anfangsschätzwerte TEE3, THETA3 und PHI3 erhalten werden.Die neue Gruppe der Anfangs Schätzwerte wird bei den Schritten 241 bis 244 benutzt, bei denen die Differentialgleichungen erneut integriert werden, damit entsprechend der neuen Gruppe von Anfangswerten eine

609883/0973

Gruppe von Fehlerausdrücken SAYE16, SAVE26 und SAVE36 erhalten wird. Diese Fehlerausdrücke, die Fehler des Y-Werts, des Werts der Kabelspannung und des Werts des Kabelneigungswinkels repräsentieren, werden bei den Schritten 245 bis 246 bezüglich vorgewählter oberer G-renzwerte EPSLHTI, EPSHT2, und EPSLN3 geprüft. Wenn einer dieser G-renzwerte von einem Fehlerausdruck überschritten wird, erfolgt im Flußdiagramm ein Übergang zum Schritt 248.

BeimSehritt 248 werden bei der Vorbereitung für den nächsten Durchlauf durch die innere Schleife die Abstände des Werts TEE3 von TEE1 und TEE2 bestimmt. Bei den Schritten 249 bis 251 wird der Wert TEE2 für den nächsten Durchlauf durch die innere Schleife auf den neu berechneten Wert TEE3 eingestellt. Der Wert TEE1 kann entsprechend dem Ergebnis der Prüfung beim Schritt 249 für den nächsten Durchlauf durch die Schleife modifiziert werden oder nicht. In ähnlicher Weise werden bei den Schritten 252 bis 259 die Anfangsbedingungen für THETA und PHI für den nächsten Durchlauf durch die innere Schleife in größtmöglicher Annäherung an die neu. berechneten Anfangswerte gewählt.

Beim Schritt 260 wird die Differenz zwischen den neuen Werten TEE2 und TEE1 bezüglich einem Minimum TMIET geprüft. Wenn diese Differenz äußerst klein wird, dann führt dies beim Schritt 234 zu einem unbegrenzten Ergebnis. .'Wenn die Antwort auf die beim Schritt 260 gestellte Frage "Ja" lautet, wird demnach beim Schritt 261 der neue Wert TEE2 so eingestellt, daß er wenigstens einem Minimumabstand von;2EE1 -entspricht. -

Ebenso werden bei den Schritten 262 bis 265 die Anfangswertgruppen für PHI und THETA geprüft und auf Minimaldifferenzen eingestellt, falls dies erforderlich ist.

609883/0973

Beim Schritt 266 wird der für die innere Grenzschleife vorgesehene Zähler IIT fortgeschaltet und beim Schritt 247 bezüglich eines Iterationsschleifenzählers HIIER geprüft. Wenn beim Schritt 267 die Iterationsgrenze nicht überschritten ist, erfolgt im Flußdiagramm eine Rückkehr vom Punkt A für den nächsten Durchlauf durch die innere Grenzschleife. Beim zweiten und den folgenden Durchläufen durch die innere Grenzschleife ergibt sich der einzige Unterschied gegenüber dem ersten Durchlauf durch diese Schleife aus der Tatsache, daß·die Integration und die Fehlerabschätzung der Schritte 229 bis 233 bereits beim letzten Durchlauf durch die Schleife ausgeführt worden ist. Dementsprechend erfolgt beim Schritt 228 bei diesen anschliessenden Durchlaufen im Flußdiagramm eine ,Verzweigung zum Punkt 5095_} wo diese bestimmten Fehlerausdrücke'.gleich den beim letzten Durchlauf durch die innere Schleife berechneten Vierten gesetzt wird. Das Verlassen der'inneren Grenzschleife erfolgt auf einem von zwei Wegen. Wenn letzlich keine Konvergenz erreicht wird, dann wird beim Schritt 267 die Iterationsgrenze überschritten, und das Flußdiagramm fährt beim Schritt 269 fort, wo das Programm verändert wird. In der normalen Betriebsart wird vor diesem Ereignis Konvergenz erreicht, und die Prüfungen bei den Schritten 245, 246'und 247 haben ein bejahendes Ergebnis, wodurch eine Verzweigung zum Punkt 600 verursacht wird. Es sei daran erinnert, daß der,bei allen Schätzungen bis zu dieser Stelle verwendete Wert von VSLOPE beim Schritt 201'gleich einem AnfangsSchätzwert VS1 eingestellt worden ist. Der Zähler TPLAG der äußeren Schleife wurde auch beim Schritt auf den Wert -1 eingestellt. Demnach hat die Prüfung beim Schritt 270 ein verneinendes Ergebnis, doch die beim Schritt 271 gestellte Frage erhält die Antwort "Ja", so daß im Flußdiagramm mit dem Schritt 272 fortgefahren wird.Hier wird

609883/0973

der Wert TSLOPE auf.einen neuen Schätzwert YS2 zur Vorbereitung auf den nächsten Durchlauf durch die äußere Schleife eingestellt. Der Fehlerwert THIeERI des Winkels Θ, der dem Querströmungsgeschwindigkeitsgradienten VS1 entspricht, wird ebenfalls beim Schritt 272 berechnet, und der Stand des Außenschleifenzahlers wird um 1 erhöht. Die Ablauffolge kehrt dann Vom Schritt 273 zum Punkt 1 zurück, damit die innere G-renzschleife für den neuen Wert von VSLOPE, d.h. für den Wert VS2,'ausgelöst wird.

Nachdem innerhalb der inneren Schleife Konvergenz erreicht worden ist, erfolgt im Flußdiagramm wieder eine Rückkehr zum Punkt 600. In diesem Pail hat der Zähler Ii¹LAG- den Stand 0, und die Antwort auf die beim Schritt 270 gestellte Präge ist "Ja". Es erfolgt die Fortsetzung am Punkt 3, und beim Schritt 274 wird der Fehler THIER2 des Winkels 0 entsprechend dem Querströmungsgeschwindigkeitsgradienten VS2 berechnet. Die zwei Schätzwerte des Geschwindigkeitsgradienten und die resultierenden Fehlerausdrücke werden beim Schritt 275 kombiniert, damit sich ein verbesserter Schätzwert des Querströmungsgeschwindigkeitsgradienten VS3 ergibt. Die Größe VSLOPE wird gleich VS3 gesetzt, und der Zähler wird beim Schritt 275 fortgeschaltet. Vom Schritt 276 erfolgt eine Rückkehr zum Punkt 1 für den nächsten Durchlauf durch die innere Grenzschleife.

ITachdem in der inneren Grenzschleife Konvergenz erzielt worden ist, erfolgt wieder eine Rückkehr zum Punkt 600. In diesem Fall hat der Stand des Zählers IFLAG jedoch den Wert 1, und keine der in den Schritten 270, 271 oder

277 durchgeführten Prüfungen hat ein bejahendes Ergebnis«,¹ Die Ablauffolge geht somit zum Punkt 4 über.-BeimSchritt

278 werden die Größen IHIER1 und THIER2 auf die Fehlerwerte eingestellt, die den GeschwLndigkeitsgradienten VS2

609883/0973

bzw. YS3 entspreChen.. Wenn sich, beim Schritt 279 der Fehlerausdruck THTER2 als geringer als ein vorbestimmter G-fenzwert EPD erweist, dann wird das Programm beim Schritt 280 geändert, da Konvergenz erreicht worden ist.

Wenn heim Schritt 279 keine Konvergenz festgestellt wird, dann werden bei den Schritten 281 bis 284 aus der Gruppe TB-I ,YS2 und YS3 zwei dem Wert YS3 am nächsten liegende Anfangswerte ausgewählt. Diese zwei Werte werden beim Schritt 285 in der Newton-Raphson-Iterationsformel dazu verwendet,. einen verbesserten AnfangsSchätzwert YS3 zu berechnen. Der Außenschleifenzähler wird beim Schritt fortgeschaltet, und es erfolgt eine Rückkehr vom Schritt 287 zum Punkt 1 für die Durchführung eines Durchlaufs durch die innere Schleife■unter Anwendung des neuerlich g eschätzt en Werts des Querströmungsgeschwindigkeitsgradienten. Die Außenschleifeniteration des Querströmungsgeschwindigkeitsgradienten wird solange fortgesetzt, bis entweder beim Schritt 279 Konvergenz festgestellt wird, oder bis ohne Erzielung der Konvergenz fünf Iterationen ■ durchgeführt worden sind, worauf das Programm beim Schritt 288 beendet wird. Im Normalfall, bei dem Konvergenz erreicht wird, sind die .für die Anfangsbedingungen der Schleppkabelspannung, des GierungswinkeIs O und des Neigungswinkels 0 bestimmten Werte sowie der für den beim letzten Durchlauf durch die Außenschleife bestimmte Wert des Quersträmungsgeschwindigkeitsgradienten die richtigen Werte für diese Größen. Die beim letzten Durchlauf durch die Außenschleife ausgeführte letzte Integration ist unter Anwendung dieser Größen durchgeführt worden, so daß demnach die Vierte der bei der Integration bestimmten Ortsgrößen x, y und ζ das lagemodell für das Schleppkabel definieren. Nachdem nun das Lösungsmodell für den zwischen dem Heck des Schiffs 10 und dem ersten Tiefensteuerkörper verlaufenden

609883/0973

Abschnitt des Schleppkabels erhalten worden ist, kann die Prozedur wiederholt werden, damit ein Model 1 für den Abschnitt des Schleppkabels zwischen dem ersten und dem zweiten Tiefensteuerkörper erhalten wird. Die aus dem ersten Modell am ersten liefensteuerkörper erhaltenen Parameterwerte werden als G-renzwerte für die Anwendung des Fewton-Raphson-Yerfahrens auf den Abschnitt des Schleppkabels zwischen dem ersten und dem zweiten Tiefensteuerkörper verwendet. ITachdem ein Modell für den zwischen dem ersten und "dem zweiten Tiefensteuerkörper liegenden Schleppkabelabschnitt entwickelt worden ist, wird es für die G-renzwerte am zweiten Tiefensteuerkörper für den zwischen dem zweiten und dem dritten Tiefensteuerkörper liegenden Schleppkabelabschnitt verwendet. Auf diese Weise schmtet die Lösung längs des Schleppkabels fort bis der entfernteste Abschnitt des Schleppkabels im Modell nachgebildet ist.

Ein weiteres "Verfahren zur Erzeugung eines Modells der' Schleppkabelposition besteht darin, daß die oben angegebenen Differentialgleichungen nur für den Abschnitt des Schleppkabels gelöst v/erden, der sich zwischen dem Schiff 10 und dem ersten Tiefensteuerkörper befindet. Es hat sich gezeigt, daß die Lage der Projektion des seismischen Schieppkabeis auf die Viasseroberfläche in jedem Zeitpunkt gut durch, eine Hyperbel dargestellt v/erden kann. Die x-und z-Koordinaten verschiedener Abtastpunkte längs des zwischen dem Schiff 10 und dem ersten Tiefensteuerkörper befindlichen Schleppkabelabschnitts werden aus dem Modell bestimmt, das durch Integrieren der Differentialgleichungen erhalten wird. Es wird dann eine Hyperbel bestimmt, die im Sinne des geringsten mittleren quadratischen Fehlers (LMSE) am besten zu den Orten dieser Abtastpunkte paßt. Die Hyperbel wird dann sowohl für den zwisehen dem Schiff 10 und dem ersten Tiefensteuerkörper befindlichenSchleppkabelabschnitt

6098 8 3/0973

2631534

als auch für die anderen Abschnitte des Schleppkabels als Modell für die Schleppkabelpositiai verwendet.

Die allgemeineForm der zu bestimmenden Hyperbel ist in der nachfolgenden Gleichung (H) angegeben, in der die Größen Z(I) und X(I) die bekannten Koordinaten des I-ten Ab&astpunkts längs des zwischen dem Schiff 10 und dem ersten Tiefensteuerkörper liegenden Schleppkabelabschnitts sind.

[Z(I) + d]² = 0(1) [x(l) + CC]. ² + 0(2) ' (14) wobei gilt: I = 1, 2 ..., Ή.

Die Größen C(1) und C(2) sind Parameter der Hyperbel .die bestimmt werden massen. Zusätzlich zu diesen zwei unbekannten Größen gilt auch, daß der Ursprungspunkt der Hyperbel, die im Sinne des kleinsten mittleren quadratischen Fehlers am besten zu den Datenabtastwerten paßt, nicht am Heck des Schiffs 10, d.h. am Ursprung des x-z-Koordinatensystems auftritt. Daher muß auch der Wert der Größe D, nämlich die z-Aehsenverschiebung der am besten passenden Hyperbel, und die Größe CC, nämlich die x-Achsenverschiebung der am besten passenden Hyperbel bestimmt werden. Wenn sie richtig bestimmt worden sind, definieren die Größen D und CC dann einen Ursprungsvektor, der sich vom Ursprung des x-y-Koordinatensystems zum Ursprung der am besten passenden Hyperbel erstreckt.

Für den Beginn der Bestimmung dieser Hyperbel werden Werte für die Größen D und CC angenommen.Dabei kann zur Vereinfachung die Gleichung (14) so umformuliert werden, wie in der nachfolgenden Gleichung (15) angegeben ist.

ZZ(I) = C(1) · XX(I) + C(2) (15)

609 383/0973

wobei gilt: ΖΖ(Ϊ)= (z(l) + D] ² XX(I)= [X(I) + CC) ²

Die Vierte von C(1) und C(2), die die Hyperbel mit dem kleinsten mittleren quadratischen Fehler für bekannte Abtastwerte definieren, werden dann entsprechend der Matrixgleichung (16) bestimmt.

σ =

wobei gilt:

(16)

F =

P¹ ist

XX(1)
XX(2)

Z =

zz(i)

ZZ(2)

ZZ(N)_-

■die Transpornierte von F. Sachdem die ¥erte dieser Koeffizienten bestimmt worden sind, kann die Gieichung (17) verwendet werden, um den mittleren quadratischen Fehler (MSE) wischen der Hyperbel und den bekannten Abtastwerten zu bestimmen.

MSE =

[ZZ(I) - C(1) · XX(I) - C(2)]

(17)

1=1

Diese Hyperbel ist jedoch nur dann die echte Hyperbel mit dem kleinsten mittleren quadratischen Fehler, wenn der Ursprungsvektor für die von den' angenommenen ¥erten für D und CC bestimmte Hyperbel der echte Ursprungsvektor, ist. Dies ist natürlich nur selten der Fall. Zum Suchen eines "besseren Ursprungsvektors werden Werte von D und CC variiert, und die eben beschriebene Prozedur

609883/09 7-3

- 57 -

wird wiederholt, so daß sich ein neuer, mittlerer quadratischer Fehler ergibt, wie er in der Gleichung (17) definiert ist. Dieser iterative Vorgang wird solange wiederholt, "bis Werte für D und CC gefunden worden sind, die zu einem mittleren quadratischen Fehler führen, der kleiner als ein vorgegebener oberer Grenzwert ist, oder bis eine vorgewählte Anzahl von Iterationen durchgeführt worden ist. -

Dieser Iterationsvorgang sei nun genauer betrachtet. Nachdem der mittlere quadratische Fehler für die Hyperbel mit dem ersten angenommenen Ursprungsvektor bestimmt worden ist, wird ein neuer beliebiger Wert für den Ursprungsvektor bestimmt. Dieser neue Ursprungsvektor hat die gieicheLänge wie der ursprünglich angenommene Ursprungsvektor, doch hat er eine beliebige Orientierung. Anschließend wird die Hyperbel mit dem kleinsten mittleren quadratischen Fehler für diesen neuen Ursprungsvektor bestimmt und ihr mittlerer quadratischer Fehler wird entsprechend der Gleichung (17) festgestellt. Die Erzeugung eines neuen beliebig orientierten Ursprungsvektors und die anschließende Bestimmung des entsprechenden mittleren quadratischen Fehlers werden wiederholt, bis ein Ursprungsvektor gefunden ist, der zu einem mittleren quadratischen Fehler führt, der kleiner als der sich bei zuvor verwendeten Ursprungsvektoren ergebende Fehler ist. Die Orientierung dieses neuen Ursprungsvektors wird für den betrachteten Zeitpunkt als die richtige Orientierung angenommen. Unter Anwendung dieser neuen Orientierung wird ein Suchprogramm nach dem goldenen Schnitt angewendet, um die Länge des Ursprungsvektor zu bestimmen, der zum kleinsten mittleren quadratischen Fehler für Ursprungsvektoren mit der neuei Orientierung führt.Nach

6098 8 3/0973

2637634

der Bestimmung dieser Länge wird die Orientierung des Vektors wieder "beliebig "bei der Suche nach einem Ursprungsvektor mit der neuen Länge variiert, der zu einem noch kleineren mittleren quadratischen Fehler führt. Wenn eine neue Orientierung tatsächlich einen kleineren mittleren quadratischen Fehler ergibt, dann wird erneut das Suchprogramm nach dem goldenen Schnitt dazu benutzt, die Länge mit dem kleinsten mittleren quadratischen Fehler für diese neue Richtung zu finden. Diese Prozedur wird abwechselnd zwischen einem Variieren der Länge und einem Variieren der Richtung des Ursprungsvektors fortgesetzt, -bis ein Ursprungsvektor gefunden wird, der einen mittleren quadratischen Fehler ergibt, der kleiner als ein vorgewählter oberer Grenzwert ist. Die Hyperbel, die diesen abschliessenden mittleren quadrat is chenFehler ergibt, ist das Modell der Schleppkabelposition.

Ein Satz von FORTRAU-Befehlen zur Durchführung der oben beschriebenen Prozedur ist als Anlage B beigefügt. Zusätzlich zu diesem Hauptprogramm enthält dieser Befehlssatz ein Funktionsprogramm FUlTC. Dieses Programm bestimmt die C-Eoeffizienten der Hyperbeln der kleinsten mittleren quadratischen Fehler, die durch die Gleichung (16) definiert sind, und es bestimmt außerdem den mittleren quadratischen Fehler für die Hyperbel, der durch die Gleichung (17) ausgedrückt wird. Die Liste enthält auch ein Unterprogramm EXPMK2, die eine 2· ·± 2-Matrix invertiert, wie es in der Gleichung (16) erforderlich ist. Der Befehlssatz stellt ferner ein Unterprogramm GSRCH, das das nach dem goldenen Schnitt arbeitende Suchprogramm zur Bestimmung der Länge mit dem kleinsten mittleren quadratischen Fehler für einen Ursprungsvektor mit

609883/0973

festliegender Orientierung ausführt. Schließlich enthält der Befehlssatz ein Funktionsprogramm YMllT1,das vom Unterprogramm G-SRCH aufgerufen wird,um zu bestimmen, welche Gruppe aus vier Variablen den Minimalwert hat.

Wie oben erläutert wurde, ergibt die Prozedur eine Hyperbel, die genau den Ort des Abschnitts des Schleppkabels repräsentiert, für den die Integration der Gleichgewicht sdifferentialgleichungen Ortswerte erzeugt bat. Da der Ort der weiter entfernt liegenden Abschnitte des Schleppkabels stark vom Ort -des.zwischen dem Schiff und den ersten Tief ens teuerkörper befindlichen Schleppkabelabschnitts beeinflußt wird, dient die Hyperbel auch dazu, den Ort der' entfernt liegenden Abschnitte des Schleppkabels zu repräsentieren..

Die Erfindung ist hier zwar im Zusammenhang mit -speziellen Ausführungsbeispielen beschrieben worden, doch ist für den Fachmann ohne, weiteres zu erkennen, daß im Rahmen der Erfindung auch weitere Abwandlungen möglich sind. ■ '

6098 8 3/0973

Claims

Patentansprüche

1. Anordnung zur Bestimmung des Verlaufs eines Abschnitts eines seismischen Seeschleppkabels, gekennzeichnet durch

(a) eine Vorrichtung zum Koppeln eines Endes des Schlepp- - kabeis mit einem- Punkt an einem Schleppschiff und

(b) eine Vorrichtung zur Erzeugung eines die S-fcellung des Schleppkabelabschnitts repräsentierenden Signals-p

2. Anordnung nach Anspruch 1, dadurch gekennzeichnet, daß der Schleppkabelabschnitt an dem beim Schleppschiff liegenden Ende des Schleppkabels liegt und daß das Signal die Stellung des Abschnitts bezüglich des Schleppschiffs repräsentiert.

5. Anordnung nach Anspruch 1, dadurch gekennzeichnet, daß der Abschnitt vom Schleppschiff entfernt liegt und daß das Signal, das die Stellung dieses Abschnitts repräsentiert auf ein Erdkoordinatensystem bezogen ist.

4. Anordnung zur Bestimmung des Verlaufs eines Abschnitts eines an ein Schleppschiff angrenzenden Abschnitts eines bei der seismischen Meeresexploration von einem Schleppschiff gezogenen Schleppkabels, gekennzeichnet durch

(a) eine, erste Winkelanzeigevorrichtung zur Erzeugung - ■ eines ersten Signals, das den G-ierungswinkel des Schleppkabels bezüglich eines Koordinatensystems des.Schiffs repräsentiert und

(b) eiiiß zweite Winkelanzeigevorrichtung zur Erzeugung eines zweiten Signals, das den Längsneigungswinkel

609883/0973

des Schleppkabels bezüglich dieses Koordinatensystems repräsentiert.

5. Anordnung nach Anspruch 4» dadurch gekennzeichnet, daß die Winkelanzeigevorrichtungen jeweils aus einer Drehmeldereinheit bestehen.

6. Anordnung nach Anspruch 4> gekennzeichnet durch

(a) eine erste Torrichtung zur Erzeugung eines dritten Signals, das das Rollen des Schleppschiffs in einem auf die Erde bezogenen Koordinatensystem repräsentiert,

(b) eine zweite Vorrichtung zur Erzeugung eines vierten Signals, das die längsneigung des Schleppschiffs in dem auf die Erde bezogenen Koordinatensystem repräsentiert und

(c) eine Signalverarbeitungsvorrichtung, die abhängig von dem ersten, dem zweiten, dem dritten und dem vierten Signal Signale erzeugt, die den G-ierungswinkel und den Längsneigungswinkel des Schleppkabels in dem auf die Erde bezogenen Koordinatensystem repräsentieren.

7. Anordnung nach Anspruch 6, dadurch gekennzeichnet, daß die erste Vorrichtung und die zweite Vorrichtung jeweils von einerPendeleinheit gebildet sind.

8. Anordnung zur Bestimmung der Position eines Schleppkabels, das bei der seismischen Meeresexploration von einem Schleppschiff gezogen wird, gekennzeichnet durch

609883/0973

2631834

(a) eine Vorrichtung zur Erzeugung eines Signals, das den Verlauf eines Abschnitts des Schleppkabels repräsentiert, und

(b) eine Signalverarbextungsanordnung, die abhängig

von dem Signal ein Modell der Position des Schleppkabels erzeugt.

9· Anordnung nach Anspruch 8, dadurch gekennzeichnet, daß das Signal den Verlauf des Abschnitts des Schleppkabels repräsentiert, der an das Schiff angrenzt.

10. Anordnung nach Anspruch 9, dadurch gekennzeichnet, daß das Signal eine G-ierungskomponente und eine Längsneigungs— komponente hat und daß der G-ierungswinkel des Modells bei dem Schiff im wesentlichen gleich der G-ierungskomponente i st.

11. Anordnung zur Bestimmung der Position eines seismischen Schleppkabels, gekennzeichnet 'durch-

(a) eine Vorrichtung zur Erzeugung eines Signals, die den Verlauf eines vom Schleppschiff entfernt liegenden Abschnitts des Schleppkabels repräsentiert und

(b) eine Signalverarbeitungsanordnung, die abhängig von dem

1
Signal ein Modell der Position des Schleppkabels erzeugt.

12. Anordnung nach Anspruch 11, dadurch gekennzeichnet, daß die Signalerzeugungsvorrichtung von einem Kompaß gebildet ist, der ein den seitlichen Verlauf des Schleppkabelabschnitts in einem auf die Erde bezogenen Koordinatensystem repräsentierendes

609883/0973

Signal erzeugt, und daß das Modell einen Abschnitt enthält, der im wesentlichen tangential zu dem seitlichen Verlauf liegt.

13. Anordnung zur Bestimmung der Position eines seismischen Seeschleppkabels, gekennzeichnet durch

(a) eine Vorrichtung zur Erzeugung eines ersten Signals, das den V erlauf eines ersten Abschnitts d?s Schleppkabel repräsentiert, '

(b) eine Vorrichtung zur Erzeugung eines zweiten Signals, das den Verlauf eines zweiten Abschnitts des Schleppkabels repräsentiert und

(c) eine Anordnung, die abhängig von dem ersten Signal und Ton dem zweiten Signal ein Modell der Position des Schleppkabels erzeugt.

14. Anordnung nach Anspruch 13, dadurch gekennzeichnet, daß der erste Abschnitt an das' Schleppschiff angrenzt, daß das erste Signal sowohl eine Gierungskomponente als auch eine Längsneigungskomponente aufweist und daß der Gierungswinkel des Modells bei dem Schleppschiff im wesentlichen gleich der Gierungskomponente ist.

15. Anordnung nach Anspruch 14, dadurch gekennzeichnet, daß die Vorrichtung zur Erzeugung des zweiten Signals von einem Kompaß gebildet ist, der den seitlichen V erlauf des zweiten Abschnitts anzeigt, und daß das Model 1 ein Segment enthält, äis im-wesentlichen tangential zu dem seitlichen Verlauf liegt.

16. Anordnung nach Anspruch 15, dadurch gekennzeichnet, daß wenigstens zwei Zugmesser vorgesehen sind, die Signale erzeugen, die die Kabelspannung an wenigstens zwei Punkten

• 609883/0973

längs des Schleppkabels repräsentieren , und daß Vorrichtungen vorgesehen sind, die abhängig von den Zugmessern die Drehung in wenigstens einem Abschnitt des Schleppkabels bestimmen.

17. Anordnung nach Anspruch 13, dadurch gekennzeichnet, daß das Modell im wesentlichen in Echtzeit erzeugt wird.

18. Anordnung zur Bestimmung der Position eines seismischen Seeschleppkabels,gekennzeichnet durch

(a.) eine Vorrichtung zur Erzeugung eines Signals, das den Verlauf eines Abschnitts des Schleppkabels repräsentiert,

(b) eine Vorrichtung, die abhängig von dem Signal ein Modell der Schleppkabelposition erzeugt, wobei wenigstens ein Abschnitt des Modells der Lösung von Differentialgleichungen entspricht, die den Gleichgewichtsverlauf des Schleppkabels beschreiben, und

(c) daß diese Lösung wenigstens eine von dem Signal gelieferte Schleppkabelgrenzbedingung erfüllt.

19· Anordnung nach Anspruch 18, dadurch gekennzeichnet, daß der Schleppkabelabschnitt an das Schleppschiff angrenzt, daß das Signal eine Gierungskomponente enthält und daß die Grenzbedingung diese Gierungskomponente ist.

20. Anordnung nach Anspruch 19» dadurch gekennzeichnet, daß eine Vorrichtung zur Erzeugung eines zweiten Signals vorgesehen ist, das die vom Schleppkabel am Schiff ausgeübte Zugspannung repräsentiert, und daß die Zugspannung eine zweite Grenzbedingung bildet, die von dieser

83/0973

Lösung erfüllt wird.

21. Anordnung nach Anspruch 20, dadurch gekennzeichnet, daß die Lösung außerdem dritte und vierte G-renzbedingungen erfüllt, die von einer Längsneigungskomponente des Signals und von der Tiefe an der Stelle gebildet sind, an der das Schleppkabel bezüglich der gemessenen Tiefe einer anderen Stelle des Schleppkabels in das Wasser eintritt.

22. Anordnung zur Erzeugung eines die Position eines seismischen Seeschleppkabels repräsentierenden Modells, gekennzeichnet durch

(a) eine Vorrichtung zur Erzeugung eines Signals, das den Verlauf eines Abschnitts des Schleppkabels repräsentiert,

(b) eine Signalverarbeitungsanordnung, die abhängig von dem Signal ein Zwisehenmod'ell erzeugt, das die Position eines Abschnitts des Schleppkabels darstellt, wobei dieses Zwischenmodell mit der Lösung von Gleichgewichtsdifferentialgleichungen für das Schleppkabel übereinstimmt, und -

(c) eine zweite Signalverarbeitunganordnung zur Erzeugung eines die Position des Schleppkabels darstellenden

• . Modells, das in. Porm einer konischen Punktion vorliegt, und das ferner im wesentlichen dem Zwischenmodell in dem Abschnitt des Schleppkabels entspricht, das vom Zwischenmodell dargestellt ist.

6 0988 3/0973

2531634

23. Anordnung nach Anspruch 22, dadurch gekennzeichnet, daß die konische Punktion eine Hyperbel ist.

24. Anordnung zum Bestimmen der Position eines seismischen Seeschleppkabels,gekennzeichnet durch

(a) eine Vorrichtung zuar Bestimmung der Position

an wenigstens zwei Punkten des Schleppkabels und

(b) eine Vorrichtung zur Erzeugung eines Modells der Schleppkabelposition, das eine Anpassung an die zwei Punkte im Sinne des kleinsten mittleren quadratischen Fehlers darstellt.

25· Anordnung nach Anspruch 24, dadurch gekennzeichnet, daß · das Modell eine Hyperbel ist.

26. Verfahren zum Bestimmen der Position eines seismischen Seeschleppkabels, dadurch gekennzeichnet,

(a) daß der Verlauf eines Abschnitts des Schleppkabels bestimmt wird, wobei dieser Verlauf eine seitliche Komponente" in einem auf die Erde bezogenen Koordinatensystem enthält und

(b) daß ein Modell der Position des Schleppkabels erzeugt wird, wobei ein Teil des Modells, das den Abschnitt repräsentiert, eine seitliche Komponente aufweist, die mit der seitlichen Komponente des Schleppkabel-"—-Verlaufs übereinstimmt.

27. Verfahren nach Anspruch 26, dadurch gekennzeichnet, daß der Schleppkabelabschnitt an das Schleppschiff angrenzt und daß der seitliche Verlauf eines zweiten Abschnitts des Schleppkabels bestimmt wird, wobei der Teil des Modells, der dem zweiten Abschnitt entspricht, einen seitlichen Verlauf aufweist, der mit dem seitlichen Verlauf des zweiten Abschnitts des Schleppkabels übereinstimmt.

609883/097?