DE69818154T2 - Methode zur 3D Modellierung der Impedanz in einem heterogenen Volumen der Untergrundes - Google Patents
Methode zur 3D Modellierung der Impedanz in einem heterogenen Volumen der Untergrundes Download PDFInfo
- Publication number
- DE69818154T2 DE69818154T2 DE69818154T DE69818154T DE69818154T2 DE 69818154 T2 DE69818154 T2 DE 69818154T2 DE 69818154 T DE69818154 T DE 69818154T DE 69818154 T DE69818154 T DE 69818154T DE 69818154 T2 DE69818154 T2 DE 69818154T2
- Authority
- DE
- Germany
- Prior art keywords
- impedance
- model
- volumes
- priori
- waves
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 31
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 238000002847 impedance measurement Methods 0.000 claims description 4
- 239000013049 sediment Substances 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims 1
- 239000010409 thin film Substances 0.000 claims 1
- 230000006870 function Effects 0.000 description 19
- 230000008021 deposition Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000014509 gene expression Effects 0.000 description 4
- 238000011835 investigation Methods 0.000 description 3
- 230000000149 penetrating effect Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000004064 recycling Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 235000014676 Phragmites communis Nutrition 0.000 description 1
- JLQUFIHWVLZVTJ-UHFFFAOYSA-N carbosulfan Chemical compound CCCCN(CCCC)SN(C)C(=O)OC1=CC=CC2=C1OC(C)(C)C2 JLQUFIHWVLZVTJ-UHFFFAOYSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000010454 slate Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Description
- Die vorliegende Erfindung betrifft ein Verfahren zur 3D-Modellierung der Impedanz eines heterogenen Mediums oder Volumens, ausgehend von Wellen, die durch Diskontinuitäten des Mediums reflektiert sind, und zwar in Abhängigkeit von in dem Medium (Milieu) ausgesandten Wellen. Die verwendeten Wellen können elastische Wellen, ausgesandt durch eine seismische Quelle oder elektromagnetische Wellen, ausgesandt beispielsweise durch ein Radar etc., sein.
- Die vorliegende Erfindung findet Anwendungen insbesondere für Untersuchungen im Untergrund mit dem Ziel, beispielsweise Kohlenwasserstofflagerstätten aufzuspüren.
- Bekannt sind verschiedene Anläufe, um seismische Daten in akustische Impedanzen umzuwandeln, welche durch die Art der Parametrisierung des Impedanzmodells sich unterscheiden, und zwar durch die Definition der Zielfunktionen und durch die Möglichkeit, Information a priori einzuführen. Die Parametrisierung kann genau sein, die Impedanz ist dabei konstant oder linear entsprechend einer vertikalen Richtung in einer gewissen Anzahl von Schichten unbekannter Geometrie. Die unbekannten Parameter sind hier nicht nur die akustischen Impedanzen, sondern auch die Orte der Grenzflächen sowie ihre Anzahl.
- Ein Verfahren zur Modellierung eines physikalischen Parameters, wie der akustischen Impedanz eines heterogenen Mediums, ist beispielsweise beschrieben in der Patentschrift
EP 0 354 112 der Anmelderin. - Die Parametrisierung kann auch dichter sein, das Impedanzfeld wird durch die Verwendung eines 3D-Gitters diskretisiert. Die seismische Zielfunktion besteht in einer Norm eines bekannten Typs L1 beispielsweise oder L2 der Differenz oder Abweichung zwischen den synthetischen Daten und den beobachteten Daten. Die synthetischen Daten resultieren aus der Umwandlung eines gegebenen optimalen Wave lets, das auf Orte der Bohrlöcher kalibriert ist, wobei die Reihe der Reflektionskoeftizienten ausgehend von der Impedanzdiagraphie berechnet sind. Die in der geologischen Zielfunktion betrachtete Größe kann die Impedanz selbst oder gewisse ihrer Ableitungen, speziell in Richtung der Schichten, sein. Die an diese Größen gelegten Typen von Beanspruchungen oder Kräften können eine Norm vom Typ L1 oder das Quadrat einer Norm L2 der Abweichung zwischen dem optimalen Modell und dem a-priori-Modell sein, genauso gut aber wie Gleichheiten oder Ungleichheiten zwischen der gewünschten Größe und den bekannten numerischen Werten. Die am weitesten praktizierten Optimierungsmethoden sind die Methoden mit konjugiertem Gradient oder sog. Techniken des „simulierten Glühens".
- Aus Brac J. et al; „Inversion with A Priori Information: an Approach to Integrated Stratigraphic Interpretation"; in Sheriff R.E.Ed., Reservoir Geophysics, Soc. Expl. Geophys., Investigations in Geophysics 7, ist es bekannt, als geologische Zielfunktion das Quadrat der Norm L2 der Impedanzabweichung und der Abweichung des Impedanzgradienten nach Projektion auf die Hangrichtung zu wählen, gesehen in der Vertikalebene der seismischen Linien. Auf diese Weise besteht die Inversion eines 3D-Blocks aus mehreren 2D-Mehrspur-Inversionen, die an seismische Linie nach seismischer Linie angelegt werden, selbst wenn das a-priori-Modell in 3D ausgearbeitet wurde.
- Die Verlässlichkeit und die Genauigkeit der Speichermodelle hängen ab vom Integrationsgrad zwischen den seismischen und geologischen Daten. Ein realistisches Speichermodell soll, bei einem noch so frühen Stadium des Realisationsprozesses wie möglich, relativ exakte seismische Daten nach Summierung (post-stack) und zwar der wenig zahlreichen diagraphischen Daten des Bohrlochs seitlich, jedoch vertikal genaue sowie auch interpretative Daten, die von einer regionalen geologischen Studie geliefert werden, kombinieren. Die stratigraphische Inversion nach Summierung (post-stack) ist ein wichtiger Schritt zur Integration der Geowissenschaften.
- In der folgenden Beschreibung wird der allgemeine Ausdruck „Volumen aus Blättern" oder mit einer „Schieferung versehenes Volumen" verwendet, um das zu modellierende Medium zu definieren, und zwar in seinem geometrischen Sinn genommen. Es handelt sich um ein Volumen, das aus einer Gesamtheit von Blättern oder zusam menhängenden und getrennten Oberflächen gebildet ist, deren Vereinigung das Volumen bildet.
- Das Verfahren nach der Erfindung
- Das Verfahren nach der Erfindung hat zum Ziel, ein optimales Modell in drei Dimensionen (3D) der Impedanz zu erhalten, die eine Zone eines heterogenen Mediums bezüglich Wellen (elastische, elektromagnetische) bietet, wobei dieses Medium mehrere Volumina umfasst, die je aus einer Anordnung von Blättern oder Schichten gebildet sind (im Folgenden Blätter genannt), wobei die verschiedenen Blätter definiert werden, ausgehend von aufgezeichneten Daten entsprechend Wellen, die aus dem heterogenen Medium in Abhängigkeit von Wellen stammen, die hier hinein übertragen werden, bei dem man ein geometrisches 3D-Modell konstruiert, das die verschiedenen Blätter-Volumina umfasst, und man konstruiert ein a-priori-3D-Modell der Impedanz, ausgehend von dem geometrischen Modell des Mediums und einer Vielzahl von auf verschiedenen Tiefen genommenen Impedanzmessungen.
- In seiner allgemeinen Definition zeichnet sich das Verfahren dadurch aus, dass es umfasst:
-
- – die Bildung einer globalen Zielfunktion einschließlich einer Impedanzfunktion und einer Funktion, die relativ zu den aufgezeichneten Daten definiert wurde;
- – die Auswahl eines 2D-Kovarianzmodells längs der Blätter dieser Volumina, wodurch die Unsicherheiten am a-priori-Impedanzmodell modelliert werden, um diese Impedanzfunktion für die gesamte Zone zu definieren; und
- – die Bildung des optimalen Impedanzmodells durch eine stratigraphische globale Inversion der aufgezeichneten Daten, unter Verwendung des a-priori-Impedanzmodells unter Minimierung der globalen Zielfunktion.
- Das Verfahren lässt sich insbesondere anwenden, um ein optimales 3D-Modell der akustischen Impedanz einer unterirdischen Zone zu erhalten, wobei die verschiedenen Volumina hier Sedimenteinheiten sind und die Impedanzmessungen auf verschiedenen Tiefen wenigstens eines die Zone durchsetzenden Bohrlochs erhalten wurden und die Blätter der verschiedenen Volumina Isochrone der Abscheidung in jeder dieser Einheiten sind.
- Gemäß einer bevorzugten Verwirklichungsform wählt man ein exponentielles isotropes oder anisotropes Kovarianzmodell aus.
- Die Funktion bezüglich der aufgezeichneten Daten wird beispielsweise gebildet durch das Quadrat der Norm L2 der Abweichung zwischen den synthetischen Daten, die ausgehend vom a-priori-Modell gebildet wurden und den aufgezeichneten Daten, die Impedanzfunktion wird beispielsweise gebildet einerseits durch das Quadrat der Norm L2 der Impedanzabmessung und andererseits das Quadrat der Norm L2 des Impedanzgradienten nach Projektion des Gradienten auf die Tangentenebene an ein lokales Blatt der Volumina (Einfallebene beispielsweise).
- Man wählt ein exponentielles 2D-Kovarianzmodell, das isotrop oder anisotrop sein kann.
- Man formt das 3D-a-priori-Impedanzmodell beispielsweise durch einfache Interpolation oder durch Kriging längs der Abscheidungsisochronen der Werte der Impedanz, die in wenigstens einem das Medium durchsetzenden Bohrloch bekannt sind.
- Das Verfahren nach der Erfindung bietet eine flexible Näherung um die Impedanzen im Band der verwendeten Frequenzen, beaufschlagt durch die aufgezeichneten 3D-Daten, die Impedanzlogs (gemessen im Bohrloch beispielsweise) und die bekannten Informationen des Mediums (die geologischen Daten sowohl hinsichtlich stratigraphischer wie struktureller Art beispielsweise) abzuschätzen.
- Weitere Merkmale und Vorteile des Verfahrens nach der Erfindung werden sich beim Lesen der nachstehenden Beschreibung eines nicht als begrenzend anzusehenden Verwirklichungsbeispiels unter Bezug auf die beiliegenden Zeichnungen ergeben, in denen
- die
1A ,1B einen Impedanzquerschnitt i zeigen, der ein Bohrloch jeweils nach einer multiplen 2D-Inversion und einer 3D-Inversion schneidet; -
2 zeigt die Inkohärenzen zwischen benachbarten Linien eines Impedanzquerschnitts i, die nach einer Multi-2D-Inversion existieren können; -
3 zeigt gut, in wieweit die 3D-Inversion dazu beiträgt, die in2 sichtbaren Inkohärenzen zu schwächen; -
4 zeigt verschiedene Kurven (a), (b), (c) der Variation des Korrelationskoeffizienten zwischen den Reflektivitäten, die zwei benachbarten Impedanzspuren zugeordnet sind, jeweils den Fällen der2 und3 und dem Fall der nachstehenden5 entsprechend, wobei die geringen Korrelationen nach links in der Figur auf das Vorhandensein von Verwerfungen zurückzuführen sind; -
5 zeigt einen Horizontalschnitt an der Basis des betrachteten Blocks nach einer 3D-Inversion, wo die Unsicherheitsparameter zugunsten der seismischen Information gewählt wurden; - die
6A und6B zeigen gleichzeitig die invertierten Impedanzspuren und die Logs (oder Diagraphien) der Impedanz, wie sie in Bohrlöchern erhalten wurden; und -
7 zeigt ein Organigramm in Erläuterung des Verfahrens nach der Erfindung. - Das Verfahren umfasst, man hat dies in einem ersten Anlauf gesehen, die Konstruktion eines geometrischen 3D-Modells, das mehrere Blattvolumina umfasst, die im beschriebenen Ausführungsbeispiel Sedimenteinheiten sind, deren Blätter Abscheidungsisochrone sind. Man konstruiert dann ein a-priori-Modell der Impedanz, ausgehend von diesem geometrischen Modell und einer Gruppe von Impedanzmessungen, die beispielsweise längs eines oder mehrerer die explorierte Zone durchsetzender Bohrlöcher genommen wurden. Die Verwirklichung dieser beiden Stufen ist beispielsweise beschrieben in:
-
- – Joint EAEG/EAPG/EAGO Multidisciplinary Workshop on „Developing New Reservoirs in Europe", 1994, St Petersburg; Extended Abstracts Book.
- Eine wesentliche Stufe des Verfahrens nach der Erfindung besteht darin, eine globale Zielfunktion zu minimieren, die an der Gesamtheit des 3D-Modells vermittels einer Methode des konjugierten Gradienten berechnet wurde. Diese globale Zielfunktion wird gebildet durch zwei Ausdrücke. Der erste, bei dem es sich um den seismischen Ausdruck handelt, besteht aus dem Quadrat der Norm L2 der Abweichung zwischen den synthetischen Daten und den beobachteten Daten. Der zweite, bei dem es sich um den Impedanzausdruck handelt, wird gebildet einerseits durch das Quadrat der Norm L2 der Impedanzabweichung und andererseits durch das Quadrat der Norm L2 des Abweichungsgradienten der Impedanz nach Projektion auf die Einfallebene. Dieses stratigraphische 3D-Inversionsverfahren gibt in einem einzigen Durchgang einen optimalen Impedanzblock 3D, der gleichzeitig kompatibel mit dem seismischen 3D-Block und dem a-priori-3D-Modell ist.
- Vollständige stratigraphische 3D-Inversion
- Das vollständige 3D-Inversionsverfahren basiert auf zwei praktischen folgenden Hypothesen, die üblicherweise im Rahmen von Untersuchungen an Speichern verwendet werden, d. h.:
-
- – die Funktion, die an den Impedanzspuren auf dem Gebiet der Zeit die seismischen Spuren zuordnet, wird berechnet durch die Convolution eines gegebenen Wavelets mit den Reflexionskoeffizienten jeder Spur;
- – die geologische Struktur umfasst gemäßigte Einfallebenen der Verwertungen mit geringer Rückstrahlung und die geringen seitlichen Geschwindigkeitsveränderungen.
- – Im übrigen ist anzusetzen, dass das seismische Wavelet konstant über die Gesamtheit des betrachteten 3D-Blocks ist und dass das die seismischen Spuren verändernde Rauschen ein weißes Rauschen (Gaußsches Rauschen) oder eine weiße Störung ist, die von einer Spur zur anderen nicht korreliert ist.
- – Angenommen wird ebenfalls, dass die Unsicherheiten bezüglich des a-priori-Impedanzmodells modelliert werden vermittels einer exponentiellen 2D-Kovarianz. Zahlreiche Kovarianzmodelle können verwendet werden. Gewisse werden einfach durch eine Varianz und eine Korrelationslänge definiert, die Wahl eines exponentiellen Kovarianzoperators ist besonders interessant in dem Ausmaß, wo die seiner Umkehrung zugeordnete Norm bekannt und leicht zu berechnen ist.
- Die seismische Zielfunktion
-
- In dieser Gleichung σs = SRMS/p ist die Unsicherheit bezüglich der seismischen Amplitude. SRMS ist der wirkliche Wert der seismischen Daten und p das Verhältnis von Signal-Rauschen. R bedeutet die Spur des Reflexionskoeffizienten und W das Wavelet, das nach Kalibrierung der Bohrlochdaten und der seismischen Daten erhalten wurde. S obs / ijk ist die i-te Probe der j-ten Spur des k-ten Querschnitts oder der Linie des seismischen 3D-Blocks.
- Die Impedanzzielfunktion
- In jeder Sedimenteinheit wird angenommen, dass die Kovarianzfunktion 2D längs jeder Korrelationsoberfläche die folgende ist: wo δZ(x) die Impedanzabweichung bezogen auf das a-priori-Modell am Ort x und δZ(y) die entsprechende Abweichung bei y ist. σ2 ist die Varianz und λ die Korrelationslänge. Der Kovarianzoperator ist dann:
-
- Während die 2D-Inversion, wie man weiß, ein Kovarianzmodell ID längs der Korrelationslinien bedeutet, bedeutet die 3D-Inversion ein 2D-Kovarianzmodell längs der Korrelationsflächen (den Abscheidungsisochronen), abgeleitet aus der stratigraphischen und strukturellen Interpretation.
- Die Diskretisierung von ∥φ∥ 2 / c-1, ergibt die geologische Zielfunktion.
- Der Frequenzinhalt der invertierten Impedanzen
- Die im a-priori-Modell enthaltenen Frequenzen, die größer als die durch den seismischen Bandpass gelieferten sind, werden nicht beim Inversionsverfahren geregelt. Aus diesem Grunde werden sie im allgemeinen verworfen. Diese Niederfrequenzinformationen nach Inversion stammen nicht aus den seismischen Daten sondern aus dem a-priori-Impedanzmodell. In zahlreichen Fällen stammen diese Niederfrequenz informationen aus Bohrlöchern und haben eine wichtige Bedeutung auf dem Gebiet der Vorhersage quantitativer Eigenschaften der Speicher, ausgehend von invertierten Impedanzen. So könnten zusätzliche Niederfrequenzinformationen, wie seismsiche Geschwindigkeiten bei der Konstruktion des a-priori-Impedanzmodells Verwendung finden.
- Wahl der Unsicherheitsparameter
- Die Varianz σ2 misst die Variabilität der Impedanzabweichung bezogen auf das a-priori-Modell längs der Abscheidungsisochronen. Die Korrelationslänge λ regelt die seitliche Kontinuität dieser Impedanzabweichung längs der stratigraphischen Flächen. In Praxis kann man einen Wert für jede geologische Einheit entweder ausgehend von früheren geologischen Informationen oder ausgehend von variographischen Analysen von Werten der Impedanz in Bohrlöchern längs Abscheidungsisochronen schätzen. Ihre Werte können auch entsprechend dem Vertrauen gewählt werden, das der Interpretierende dem a-priori-Modell gibt (Anzahl von Borhlöchern, geologische an den Bohrlöchern bekannte Umgebung).
- Digitale Ergebnisse
- Die Wirksamkeit des Verfahrens nach der Erfindung wurde verglichen mit der, welche das bekannte 2D-Mehrspurinversionsverfahren gibt, das definiert ist in der EP 0 354 112 (oben genannt), basierend auf Daten eines seismischen Blocks 3D, der in der Nordsee gewonnen wurde. Dieser Block erstreckt sich über 8 km2 zwischen 2,0 und 2,8 s. Das Gitter bestand aus 35 Linien mit 130 Spuren und 201 Zeitproben, was etwa 900 000 Knoten entspricht. Die Knotennahmeschritte liegen bei 75 m zwischen den Linien, 25 m zwischen den Spuren und zeitlich 4 ms. Neun Diagraphien der Impedanz wurden zum Aufbau des a-priori-Modells gewählt.
- 3D-Inversion verglichen mit der Mehrspur-2D-Inversion
- Indem man das gleiche 3D-a-priori-Modell verwendet, realisiert man gleichzeitig Mehrspur-2D-Inversionen und 3D-Inversionen am Datensatz. Die
1A und1B bezüglich des das Bohrloch durchsetzenden Querschnitts (durchgezogene weiße Linie in2 ) und die2 und3 , die sich auf den Schnitt beziehen, der längs eines seismischen Horizonts gelegt wurde, zeigen die Verbesserung der Kontinuität längs der stratigraphischen Flächen, die durch die 3D-Inversion erreicht wurde. F repräsentiert hier eine Verwertung. Diese Beobachtung wird durch4 bestätigt, wo die Kurven (a), (b), (c) den Korrelationskoeffizient zwischen den Reflexionen darstellen, die den Spuren der optimalen Impedanz der Linien 8 (weiß gestrichelt in2 ) und 9 bei gleicher Abszisse x zugeordnet sind, und zwar jeweils nach einer Multi-2D-Inversion, einer nach einer 3D-Inversion und nach einer 3D-Inversion, wobei die seismischen Daten favorisiert sind. Für die beiden Inversionen (a) und (b) sind die Parameter die folgenden: σ2 = 30%·σ = 250 g·cm– 3·m.s.–1 und λ = 500 m, was geologische Informationen begünstigt. - Der Einfluss der Unsicherheitsparameter
- Eine andere 3D-Inversion wurde mit Unsicherheitsparametern realisiert, die zugunsten der seismischen Daten anstelle des a-priori-Impedanzmodells (σ2 = 30% σ = 500 g·cm– 3·m.s.–1 und λ = 100 m) gewählt wurden. Die seismischen Daten stellen eine wesentlich stärkere seitliche Variation als die des a-priori-Modells längs der Horizonte dar. Aus diesem Grunde ist wegen des gesteigerten Einflusses der seismischen Daten das in
5 dargelegte Prinzip „rigoroser" (seitlich verrauscht), als das der2 und3 . Mit Parametern, welche seismische Daten begünstigen, bleibt die seitliche Kontinuität der optimalen Impedanz gering. Mit die geologischen Informationen begünstigenden Parametern wird die seitliche Kontinuität nur längs der Querachse mit einer Mehrspur-2D-Inversion erhöht und wird gleichzeitig längs der Querachse und der Längsachse mit einer 3D-Inversion erhöht. - Für die Ergebnisse, die gleichzeitig mit der Mehrspur-2D-Inversion und der 3D-Inversion erhalten wurden, wurde ein Vergleich zwischen der Impedanzspur am Ort des Bohrlochs und der Impedanzspur, die von Bohrlochdiagraphien abgeleitet wurde, durchgeführt. Man kann beobachten, dass dieses Bohrloch verwendet wurde, um das a-priori-Modell zu konstruieren.
6 zeigt, dass die nach 3D-Inversion erhaltene Impedanzspur geringfügig besser als die ist, die nach der Mehrspur-2D-Inversion erhalten wurde. In dieser Figur zeigen die Symbole jeweils das Folgende: -
- – Tc, das Dach der Kreideformation
- – Bc, die Basis der Kreideformation
- – Bb, die Brentbasis
- – a), die Impedanzdiagraphie
- – b), die Impedanz nach Multi-2D-Inversion und
- – c), die Impedanz nach 3D-Inversion.
- Um die Inversionsstufe der aufgezeichneten Daten unter Verwendung des a-priori-Impedanzmodells zu realisieren, kann man eine bekannte Inversionssoftware benutzen und beispielsweise die INTERWELL ® der stratigraphischen 2D- oder 3D-Inversion, wie dies das Flussdiagramm der
7 zeigt. - Das Verfahren nach der Erfindung wurde in einer Anwendung auf die Modellierung der unterirdischen Impedanz beschrieben. Es ist klar, dass es sich auch anwenden lässt auf jeden anderen physikalischen Parameter an den aufgezeichneten Daten. Allgemeiner verlässt man nicht den Rahmen der Erfindung, indem man das Verfahren nach der Erfindung auf die Realisierung von Tomographien in sehr unterschiedlichen Gebieten, insbesondere im Gebiet der Medizinbildtechnik oder der zerstörungsfreien Prüfung anwendet.
Claims (8)
- Verfahren zum Erhalt eines optimalen dreidimensionalen (3D) Modells der Impedanz einer Zone eines bezüglich Wellen heterogenen Milieus oder Mediums, wobei dieses Milieu mehrere Volumina, je gebildet durch eine Gesamtheit dünner Schichten bzw. Blätter, umfasst, und die verschiedenen Blätter, ausgehend von aufgezeichneten Daten entsprechend Wellen definiert werden, die aus einem heterogenen Milieu in Abhängigkeit von hier hinein übertragenen Wellen stammen, bei dem man ein geometrisches 3D-Modell, das unterschiedliche Blatt- oder Dünnschichtvolumina umfasst, konstruiert, und man ein 3D-a-priori-Impedanzmodell, ausgehend vom geometrischen Modell des Milieus bzw. Mediums und von einer Vielzahl von auf unterschiedlichen Tiefen genommenen Impedanzmessungen konstruiert, dadurch gekennzeichnet, dass es umfasst: – die Bildung einer globalen Zielfunktion einschließlich einer Impedanzfunktion und einer bezüglich der aufgezeichneten Daten definierten Funktion; – die Auswahl eines 2D-Kovarianzmodells längs Blättern dieser Volumina, wodurch die Unsicherheiten hinsichtlich des a-priori-Impedanzmodells modelliert werden, um diese Impedanzfunktion für die gesamte Zone zu definieren; und – die Bildung des optimalen Impedanzmodells durch eine stratigraphische Globalinversion der aufgezeichneten Daten bei Verwendung des a-priori-Impedanzmodells, unter Minimierung der globalen Zielfunktion.
- Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass das heterogene Milieu eine unterirdische Zone ist, wobei die unterschiedlichen Volumina Sedimenteinheiten sind und die Messungen der akustischen Impedanz auf unterschiedlichen Tiefen wenigstens eines durch diese Zone geführten Bohr lochs erhalten werden, und die Blätter der verschiedenen Volumina Lagerungsisochrone jeder dieser Einheiten sind.
- Verfahren nach einem der Ansprüche 1 oder 2, dadurch gekennzeichnet, dass die im Milieu ausgesandten Wellen elastische Wellen sind.
- Verfahren nach einem der Ansprüche 1 oder 2, dadurch gekennzeichnet, dass die im Milieu ausgesandten Wellen elektromagnetische Wellen sind.
- Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass die Funktion bezüglich der aufgezeichneten Daten gebildet wird durch das Quadrat der L2-Norm der Abweichung zwischen synthetischen Daten, die ausgehend vom a-priori-Impedanzmodell gewonnen wurden sowie den aufgezeichneten Daten, und die Impedanzfunktion gebildet wird einerseits durch das Quadrat der L2-Norm der Impedanzabweichung und andererseits das Quadrat der L2-Norm des Gradienten der Impedanzabweichung nach Projektion dieses Gradienten auf die Tangentenebene an eine lokale Schicht bzw. ein lokales Blatt der Volumina.
- Verfahren nach Anspruch 5, dadurch gekennzeichnet, dass die Tangentenebene an ein lokales Blatt der Volumina eine Neigungsebene der Lagerung ist.
- Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass man ein exponentielles 2D-Kovarianzmodell, das isotrop oder anisotrop ist, wählt.
- Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass man das 3D-a-priori-Impedanzmodell durch einfache Interpolation oder durch sog. Kriging längs Lagerungsisochronen der bekannten Impedanzwerte in wenigstens einem das Milieu durchsetzenden Bohrloch bildet.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9708512 | 1997-07-04 | ||
FR9708512A FR2765692B1 (fr) | 1997-07-04 | 1997-07-04 | Methode pour modeliser en 3d l'impedance d'un milieu heterogene |
Publications (2)
Publication Number | Publication Date |
---|---|
DE69818154D1 DE69818154D1 (de) | 2003-10-23 |
DE69818154T2 true DE69818154T2 (de) | 2004-04-08 |
Family
ID=9508883
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE69818154T Expired - Fee Related DE69818154T2 (de) | 1997-07-04 | 1998-06-25 | Methode zur 3D Modellierung der Impedanz in einem heterogenen Volumen der Untergrundes |
Country Status (7)
Country | Link |
---|---|
US (1) | US6374201B1 (de) |
EP (1) | EP0889331B1 (de) |
CA (1) | CA2241293C (de) |
DE (1) | DE69818154T2 (de) |
DK (1) | DK0889331T3 (de) |
FR (1) | FR2765692B1 (de) |
NO (1) | NO319802B1 (de) |
Families Citing this family (70)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2792419B1 (fr) | 1999-04-16 | 2001-09-07 | Inst Francais Du Petrole | Methode pour obtenir un modele optimal d'une caracteristique physique dans un milieu heterogene, tel que le sous-sol |
ES2223156T3 (es) | 1999-06-03 | 2005-02-16 | Jason Geosystems B.V. | Metodo de determinacion de cambios espaciales en la estratigrafia de la estructura del subsuelo, litologia y contenido de fluidos y de reduccion del ruido sismico. |
US7412357B1 (en) | 1999-06-15 | 2008-08-12 | Philippi-Hagenbuch, Inc. | Process for three-dimensional modeling and design of off-highway dump bodies |
FR2798197B1 (fr) | 1999-09-02 | 2001-10-05 | Inst Francais Du Petrole | Methode pour former un modele d'une formation geologique, contraint par des donnees dynamiques et statiques |
EP2296013B1 (de) * | 1999-10-22 | 2016-03-30 | CGG Services (NL) B.V. | Verfahren zur Beurteilung der elastischen Parameter und der Zusammensetzungsparameter aus seismischen und echoakustischen Daten |
FR2800473B1 (fr) * | 1999-10-29 | 2001-11-30 | Inst Francais Du Petrole | Methode pour modeliser en 2d ou 3d un milieu heterogene tel que le sous-sol decrit par plusieurs parametres physiques |
US7415401B2 (en) * | 2000-08-31 | 2008-08-19 | Exxonmobil Upstream Research Company | Method for constructing 3-D geologic models by combining multiple frequency passbands |
US6853922B2 (en) * | 2001-07-20 | 2005-02-08 | Tracy Joseph Stark | System for information extraction from geologic time volumes |
US7257467B2 (en) * | 2001-11-02 | 2007-08-14 | Hagenbuch Leroy G | Method of estimating the volumetric carrying capacity of a truck body |
US7072768B1 (en) | 2003-05-02 | 2006-07-04 | Young Alan G | Method for laterally extrapolating soil property data using soil samples and seismic amplitude data within a seismic coverage area |
CN1898675A (zh) * | 2004-01-30 | 2007-01-17 | 埃克森美孚上游研究公司 | 储层模型建立方法 |
CN1898640A (zh) * | 2004-01-30 | 2007-01-17 | 埃克森美孚上游研究公司 | 储层评价方法 |
US20060041409A1 (en) * | 2004-08-20 | 2006-02-23 | Chevron U.S.A. Inc. | Method for making a reservoir facies model utilizing a training image and a geologically interpreted facies probability cube |
US20060047429A1 (en) * | 2004-08-24 | 2006-03-02 | Adams Steven L | Method of estimating geological formation depths by converting interpreted seismic horizons from the time domain to the depth domain |
US7719923B2 (en) * | 2006-06-05 | 2010-05-18 | Westerngeco L.L.C. | Low frequency model estimation |
US7901009B2 (en) * | 2006-09-14 | 2011-03-08 | Hagenbuch Leroy G | Severe application off-highway truck body |
EP2067112B1 (de) | 2006-09-28 | 2017-10-18 | Exxonmobil Upstream Research Company | Iterative umkehrung von daten aus simultanen geophysikalischen quellen |
EP2153246B1 (de) * | 2007-05-09 | 2015-09-16 | ExxonMobil Upstream Research Company | Inversion von seismischen 4d-daten |
US8547783B2 (en) * | 2007-12-12 | 2013-10-01 | Exxonmobil Upstream Research Company | Method and apparatus for evaluating submarine formations |
US8812282B2 (en) | 2008-03-21 | 2014-08-19 | Exxonmobil Upstream Research Company | Efficient method for inversion of geophysical data |
CA2726462C (en) * | 2008-08-11 | 2016-12-06 | Exxonmobil Upstream Research Company | Estimation of soil properties using waveforms of seismic surface waves |
US9146329B2 (en) | 2008-12-17 | 2015-09-29 | Exxonmobil Upstream Research Company | System and method for reconstruction of time-lapse data |
US8705317B2 (en) | 2008-12-17 | 2014-04-22 | Exxonmobil Upstream Research Company | Method for imaging of targeted reflectors |
WO2010077568A1 (en) | 2008-12-17 | 2010-07-08 | Exxonmobil Upstream Research Company | System and method for performing time-lapse monitor surveying using sparse monitor data |
US8451683B2 (en) * | 2009-04-03 | 2013-05-28 | Exxonmobil Upstream Research Company | Method for determining the fluid/pressure distribution of hydrocarbon reservoirs from 4D seismic data |
US8332154B2 (en) | 2009-06-02 | 2012-12-11 | Exxonmobil Upstream Research Company | Estimating reservoir properties from 4D seismic data |
US8537638B2 (en) * | 2010-02-10 | 2013-09-17 | Exxonmobil Upstream Research Company | Methods for subsurface parameter estimation in full wavefield inversion and reverse-time migration |
FR2957682B1 (fr) | 2010-03-16 | 2012-03-23 | Inst Francais Du Petrole | Procede de surveillance d'un site de stockage geologique de gaz par inversion stratigraphique de donnees sismiques |
US8223587B2 (en) * | 2010-03-29 | 2012-07-17 | Exxonmobil Upstream Research Company | Full wavefield inversion using time varying filters |
US8694299B2 (en) | 2010-05-07 | 2014-04-08 | Exxonmobil Upstream Research Company | Artifact reduction in iterative inversion of geophysical data |
US8756042B2 (en) | 2010-05-19 | 2014-06-17 | Exxonmobile Upstream Research Company | Method and system for checkpointing during simulations |
US8767508B2 (en) | 2010-08-18 | 2014-07-01 | Exxonmobil Upstream Research Company | Using seismic P and S arrivals to determine shallow velocity structure |
KR101908278B1 (ko) | 2010-09-27 | 2018-10-17 | 엑손모빌 업스트림 리서치 캄파니 | 풀 파동장 반전을 위한 실제적인 해결책으로서의 동시 소스 인코딩 및 소스 분리 |
US8437998B2 (en) | 2010-09-27 | 2013-05-07 | Exxonmobil Upstream Research Company | Hybrid method for full waveform inversion using simultaneous and sequential source method |
AU2011337143B2 (en) | 2010-12-01 | 2016-09-29 | Exxonmobil Upstream Research Company | Simultaneous source inversion for marine streamer data with cross-correlation objective function |
AU2012233133B2 (en) | 2011-03-30 | 2014-11-20 | Exxonmobil Upstream Research Company | Convergence rate of full wavefield inversion using spectral shaping |
CA2825469C (en) | 2011-03-31 | 2017-06-13 | Exxonmobil Upstream Research Company | Method of wavelet estimation and multiple prediction in full wavefield inversion |
EP2751710B1 (de) | 2011-09-02 | 2017-08-02 | Exxonmobil Upstream Research Company | Verwendung von projektion auf konvexe datensätze zur begrenzung von full-wavefield-inversion |
US9176930B2 (en) | 2011-11-29 | 2015-11-03 | Exxonmobil Upstream Research Company | Methods for approximating hessian times vector operation in full wavefield inversion |
EP2803043B1 (de) * | 2012-01-13 | 2018-04-04 | Geco Technology B.V. | Wellenforminversion auf basis von 3d-oberflächen |
AU2013230789B2 (en) | 2012-03-08 | 2016-02-11 | Exxonmobil Upstream Research Company | Orthogonal source and receiver encoding |
WO2014084945A1 (en) | 2012-11-28 | 2014-06-05 | Exxonmobil Upstream Resarch Company | Reflection seismic data q tomography |
CN105308479B (zh) | 2013-05-24 | 2017-09-26 | 埃克森美孚上游研究公司 | 通过与偏移距相关的弹性fwi的多参数反演 |
US10459117B2 (en) | 2013-06-03 | 2019-10-29 | Exxonmobil Upstream Research Company | Extended subspace method for cross-talk mitigation in multi-parameter inversion |
US9702998B2 (en) | 2013-07-08 | 2017-07-11 | Exxonmobil Upstream Research Company | Full-wavefield inversion of primaries and multiples in marine environment |
EP3351972A1 (de) | 2013-08-23 | 2018-07-25 | Exxonmobil Upstream Research Company | Iterative inversion feldcodierter seismischer daten an hand der erzeugung von pseudo-superquellendaten |
US10036818B2 (en) | 2013-09-06 | 2018-07-31 | Exxonmobil Upstream Research Company | Accelerating full wavefield inversion with nonstationary point-spread functions |
CN103869360B (zh) * | 2014-02-26 | 2016-08-17 | 中国石油天然气股份有限公司 | 逆冲推覆带储层波阻抗反演方法及装置 |
US9910189B2 (en) | 2014-04-09 | 2018-03-06 | Exxonmobil Upstream Research Company | Method for fast line search in frequency domain FWI |
CA2947847C (en) | 2014-05-09 | 2018-08-14 | Exxonmobil Upstream Research Company | Efficient line search methods for multi-parameter full wavefield inversion |
US10185046B2 (en) | 2014-06-09 | 2019-01-22 | Exxonmobil Upstream Research Company | Method for temporal dispersion correction for seismic simulation, RTM and FWI |
CA2947410A1 (en) | 2014-06-17 | 2015-12-30 | Exxonmobil Upstream Research Company | Fast viscoacoustic and viscoelastic full-wavefield inversion |
US10838092B2 (en) | 2014-07-24 | 2020-11-17 | Exxonmobil Upstream Research Company | Estimating multiple subsurface parameters by cascaded inversion of wavefield components |
US10422899B2 (en) | 2014-07-30 | 2019-09-24 | Exxonmobil Upstream Research Company | Harmonic encoding for FWI |
US20160061986A1 (en) * | 2014-08-27 | 2016-03-03 | Schlumberger Technology Corporation | Formation Property Characteristic Determination Methods |
US10386511B2 (en) | 2014-10-03 | 2019-08-20 | Exxonmobil Upstream Research Company | Seismic survey design using full wavefield inversion |
WO2016064462A1 (en) | 2014-10-20 | 2016-04-28 | Exxonmobil Upstream Research Company | Velocity tomography using property scans |
FR3028649B1 (fr) | 2014-11-18 | 2016-12-02 | Ifp Energies Now | Procede de construction d'un modele geologique |
US10359532B2 (en) | 2014-12-10 | 2019-07-23 | Schlumberger Technology Corporation | Methods to characterize formation properties |
EP3234659A1 (de) | 2014-12-18 | 2017-10-25 | Exxonmobil Upstream Research Company | Skalierbare planung von parallelen iterativen seismischen arbeiten |
US10520618B2 (en) | 2015-02-04 | 2019-12-31 | ExxohnMobil Upstream Research Company | Poynting vector minimal reflection boundary conditions |
SG11201704620WA (en) | 2015-02-13 | 2017-09-28 | Exxonmobil Upstream Res Co | Efficient and stable absorbing boundary condition in finite-difference calculations |
US10670750B2 (en) | 2015-02-17 | 2020-06-02 | Exxonmobil Upstream Research Company | Multistage full wavefield inversion process that generates a multiple free data set |
US11062061B2 (en) * | 2015-03-30 | 2021-07-13 | Volvo Construction Equipment Ab | System and method for determining the material loading condition of a bucket of a material moving machine |
WO2016195774A1 (en) | 2015-06-04 | 2016-12-08 | Exxonmobil Upstream Research Company | Method for generating multiple free seismic images |
US10838093B2 (en) | 2015-07-02 | 2020-11-17 | Exxonmobil Upstream Research Company | Krylov-space-based quasi-newton preconditioner for full-wavefield inversion |
RU2693495C1 (ru) | 2015-10-02 | 2019-07-03 | Эксонмобил Апстрим Рисерч Компани | Полная инверсия волнового поля с компенсацией показателя качества |
KR102021276B1 (ko) | 2015-10-15 | 2019-09-16 | 엑손모빌 업스트림 리서치 캄파니 | 진폭 보존을 갖는 fwi 모델 도메인 각도 스택들 |
US10768324B2 (en) | 2016-05-19 | 2020-09-08 | Exxonmobil Upstream Research Company | Method to predict pore pressure and seal integrity using full wavefield inversion |
CN107065034B (zh) * | 2017-06-07 | 2018-11-27 | 中海石油(中国)有限公司 | 一种基于高程差的河流相地层等时划分方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2535259A1 (de) * | 1975-08-07 | 1977-02-10 | Helmut Dipl Phys Blum | Verfahren und vorrichtung zur entdeckung geologischer inhomogenitaeten |
US4679174A (en) * | 1984-04-26 | 1987-07-07 | Western Geophysical Company Of America | Method for seismic lithologic modeling |
FR2635197B1 (fr) * | 1988-08-05 | 1990-09-28 | Inst Francais Du Petrole | Methode pour obtenir un modele representatif d'un milieu heterogene et notamment du sous-sol |
US5113192A (en) * | 1991-05-03 | 1992-05-12 | Conoco Inc. | Method for using seismic data acquisition technology for acquisition of ground penetrating radar data |
US5321613A (en) * | 1992-11-12 | 1994-06-14 | Coleman Research Corporation | Data fusion workstation |
US5583825A (en) * | 1994-09-02 | 1996-12-10 | Exxon Production Research Company | Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data |
FR2734069B1 (fr) * | 1995-05-12 | 1997-07-04 | Inst Francais Du Petrole | Methode pour predire, par une technique d'inversion, l'evolution de la production d'un gisement souterrain |
US5838634A (en) * | 1996-04-04 | 1998-11-17 | Exxon Production Research Company | Method of generating 3-D geologic models incorporating geologic and geophysical constraints |
US5798982A (en) * | 1996-04-29 | 1998-08-25 | The Trustees Of Columbia University In The City Of New York | Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models |
-
1997
- 1997-07-04 FR FR9708512A patent/FR2765692B1/fr not_active Expired - Lifetime
-
1998
- 1998-06-25 DK DK98401568T patent/DK0889331T3/da active
- 1998-06-25 EP EP98401568A patent/EP0889331B1/de not_active Expired - Lifetime
- 1998-06-25 DE DE69818154T patent/DE69818154T2/de not_active Expired - Fee Related
- 1998-07-02 US US09/108,974 patent/US6374201B1/en not_active Expired - Lifetime
- 1998-07-03 NO NO19983088A patent/NO319802B1/no not_active IP Right Cessation
- 1998-07-03 CA CA002241293A patent/CA2241293C/fr not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
CA2241293C (fr) | 2008-02-05 |
US6374201B1 (en) | 2002-04-16 |
CA2241293A1 (fr) | 1999-01-04 |
FR2765692A1 (fr) | 1999-01-08 |
DE69818154D1 (de) | 2003-10-23 |
NO319802B1 (no) | 2005-09-19 |
FR2765692B1 (fr) | 1999-09-10 |
EP0889331A1 (de) | 1999-01-07 |
DK0889331T3 (da) | 2003-10-20 |
EP0889331B1 (de) | 2003-09-17 |
NO983088D0 (no) | 1998-07-03 |
NO983088L (no) | 1999-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE69818154T2 (de) | Methode zur 3D Modellierung der Impedanz in einem heterogenen Volumen der Untergrundes | |
DE69415142T2 (de) | Auswertungsverfahren von seismischen spuren durch eine statistische kalibrierung um geologische eigenschaften abzuleiten | |
DE69519641T2 (de) | Methoden zur Erzeugung numerischer Modelle von seismischen Informationen | |
DE69806782T2 (de) | Verfahren zur verarbeitung von seismischen daten | |
DE19934520B4 (de) | Verfahren und Einrichtung zum Korrigieren von seismischen Daten | |
DE69418453T2 (de) | Verfahren zur bestimmung von der wasserbodenrefletivität in zweisenser-seismischen messungen | |
DE3885939T2 (de) | Verfahren zum Messen eines vertikalen seismischen Profils. | |
DE60212575T2 (de) | Verfahren zum spektralausgleich von seismischen daten mit offset | |
DE69415866T2 (de) | Verfahren zum Bestimmen von petrophysikalische Eigenschaften einer unterirdischen Schicht | |
DE69420593T2 (de) | Verfahren und Vorrichtung zum Erfassen von seismischen Signalen | |
DE69214251T2 (de) | Auswertung aeromagnetischer Daten mittels eines neutralen Netzwerks | |
DE69422188T2 (de) | Verfahren zur verarbeitung von datensignalen, die energie darstellen, welche durch ein medium fliesst | |
EP1421412A2 (de) | Verfahren zur bestimmung lokaler ähnlichkeit aus seismischen 3d-messdaten | |
DE3544602A1 (de) | Verfahren zur bildung einer oberflaechenkonsistenten korrektur fuer den einfluss von erdschichten geringer uebertragungsgeschwindigkeit bei der verarbeitung seismischer daten | |
DE4118772A1 (de) | Verfahren zur korrektur maritimer seismischer daten aufgrund von geschwindigkeitsunterschieden einer meeresstroemung | |
DE69804693T2 (de) | Verfahren und Anordnung zur Dämpfung von Wellenzügen in seismischen Daten | |
DE68913832T2 (de) | Methode zur Reduktion von Bohrstangemultiplen in Feldsignale. | |
DE68902950T2 (de) | Verfahren zum berechnen des organischen gehalts in sedimenten aus aufgezeichneten daten, welche durch bohrlochmessgeraete erzeugt wurden. | |
DE69409050T2 (de) | Verfahren zum Verarbeiten seismischer Daten einschliesslich Verlust-Kompensation | |
DE19933717C1 (de) | Verfahren zur seismischen Datenverarbeitung | |
DE60019778T2 (de) | Erzeugung von bildern unterirdischer formationen | |
DE68908775T2 (de) | Verfahren zur Bestimmung der Geschwindigkeit in einer Bohrstange. | |
DE69802547T2 (de) | Verfahren zur seismische daten verarbeitung | |
DE69003348T2 (de) | Verfahren für Meeresseismische Untersuchungen. | |
DE2735523A1 (de) | Detektor-aufstellung bzw. verfahren und vorrichtung zur seismischen untersuchung |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
8328 | Change in the person/name/address of the agent |
Representative=s name: VONNEMANN, KLOIBER & KOLLEGEN, 80796 MUENCHEN |
|
8364 | No opposition during term of opposition | ||
8339 | Ceased/non-payment of the annual fee |