Item

Abstract

Knowledge of spatial and temporal distribution of fluids in the subsurface is crucialin a wide range of applications. During the production of crude oil typically high sa-line produced formation water is injected into the reservoir layer, aiming to push theoil towards production wells. While oil is commonly seen as an electrical insulator,the injected saline brines are characterised by low electrical resistivity. Thus, electro-magnetic (EM) methods and especially Controlled Source Electromagnetics (CSEM)attracted an increasing interest to monitor these resistivity changes inside the reser-voir over time.This thesis mainly reports on numerical aspects of modelling and inversion of landbased CSEM with particular focus towards hydrocarbon monitoring applications.Most of the presented developments were inspired by a superordinate research pro-ject including CSEM field surveys across an actively producing onshore oil field inNorthern Germany.In producing oil fields there exists a large number of steel-cased wells. Such existingoil field infrastructure and especially the presence of metal casings significantly altersthe propagation of EM fields in the subsurface. Their spatially unfavourable dimen-sions effectively prohibits a straightforward implementation into the modelling grid.Thus I developed a new modelling approach allowing consideration of such thin butvertically extended highly conductive structures including their mutual interaction.The developed methodology had been implemented into existing modelling and in-version codes. Using the new approach to investigate the influence of metal casingson CSEM data shows that they act as additional inductively coupled vertical electricdipole sources at depth and thereby increase resolution capabilities at depth. Thepresence of metal casings can thus be exploited by optimising the source receiverlayout in such a way that the strength of these additional vertical dipole sources ismaximised.An additional working package of the superordinate project was the measurement ofvertical electric fields in a shallow observation well. However, measurements of ver-tical electric fields requires long measurement dipoles to achieve satisfactory signal-to-noise ratios. Such extended dipoles span several modelling cells and are thereforein conflict with assumptions usually made for modelling, that receivers can be rep-resented as point dipoles. I therefore expanded the modelling and inversion codesto consider the physical receiver dimensions. The new algorithm implicitly considersimperfect alignment of the receiver with the corresponding field component. Withoutthe consideration of this effect inversion of vertical electric field measurements islikely to cause erroneous results.Finally I discuss different aspects of time-lapse inversion required to track changes influid saturation over time. The cascaded inversion scheme is applied to synthetic time-lapse data for a simplified oilfield undergoing brine flushing. The influence of variousinversion parameters in particular different regularisation techniques are examined.Surface based sources and receivers typically provide low sensitivity towards deeptargets in highly conductive backgrounds. Despite that using additional constraints,in particular a model weighting scheme together with energised steel casings allowedto track resistivity changes inside the reservoir based on synthetic time-lapse data.

Wissen über räumliche und zeitliche Verteilung von Fluiden im Untergrund is unerlässlich für eine Reihe von Anwendungen. Typischerweise wird während der Förderung von Rohöl salinares produziertes Formationswasser in die ölführende Formation injiziert um das Öl-Wasser-Gemisch in Richtung der Förderbohrungen zu spülen. Während Öl als elektrischer Isolator gilt, zeichnen sich die injizierten salinaren Fluide durch eine hohe elektrische Leitfähigkeit aus. Daher erfahren elektromagnetische Methoden und insbesondere Controlled Source Electromagentics (CSEM) zunehmendes Interesse diese Änderung des elektrischen Widerstands mit der Zeit zu überwachen. Diese Arbeit beschäftigt sich im wesentlichen mit numerischen Aspekten der Modellierung und Inversion von CSEM an Land mit speziellem Fokus auf der Überwachung der Kohlenwasserstoff Produktion. Die meisten der gezeigten Entwicklungen sind entwickelt im Zuge eines übergeordneten Forschungsprojektes inklusive CSEM Feldmessungen in einem produzierenden Ölfeld in Norddeutschland. Produzierende Ölfelder sind gekennzeichnet durch eine große Anzahl von Stahl verrohrten Bohrungen. Die Anwesenheit von Stahlinfrastruktur insbesondere von Stahlschutzrohren beeinflusst die Ausbreitung von elektromagnetischen Feldern im Untergrund. Deren unvorteilhafte Geometrie erlaubt keine direkte Berücksichtigung in dem Modellierungsgitter. Daher habe ich einen neuen Modellierungsanzatz entwickelt der es erlaubt solch dünne aber vertikal ausgedehnte hochgradig leitfähige Strukturen inklusive deren gegenseitige Wechselwirkung zu berücksichtigen. Die entwickelte Methode wurde in bestehende Modellierungs- und Inversionssoftware implementiert. Mithilfe dieses neuen Ansatzes konnte der Einfluss von Stahlverrohrungen auf CSEM Daten untersucht werden. Stahlverrohrungen wirken wie zusätzliche induktiv angeregte vertikale elektrische Dipolquellen im Untergrund und helfen daher die Auflösung in der Tiefe zu erhöhen. Die Anwesenheit von Stahlverrohrungen kann daher ausgenutzt werden in dem man die Quell-Empfänger-Geometrie in einer Art und Weise optimiert, die die Stärke dieser zusätzlichen vertikalen Dipolquellen maximiert. Ein weiteres Arbeitspaket des übergeordneten Forschungsprojektes bestand in der Messung von vertikalen elektrischen Feldern in flachen Beobachtungsbohrungen. Messungen des vertikalen elektrischen Feldes erfordert lange Messdipole um ein ausreichendes Signal-Rausch-Verhältnisses zu gewährleisten. Solch ausgedehnte Dipole überspannen mehrere Zellen des Modellierungsgitters und verletzen die übliche Annahme, wonach die Länge der Empfänger vernachlässigbar ist. Daher habe ich die bestehenden Modellierungs- und Inversionsprogramme erweitert um die physischen Dimensionen von elektrischen Feld Empfängern zu berücksichtigen. Der implementierte Algorithmus berücksichtigt implizit Abweichungen der Orientierung des Messdipols von der Richtung der zu messenden Feldkomponente. Ohne dieser Berücksichtigung führt eine Inversion von vertikalen elektrischen Feld Daten zu fehlerhaften Ergebnissen. Schließlich werden unterschiedliche Aspekte von time-lapse Inversion diskutiert, welche notwendig ist um Änderungen der Fluidzusammensetzung abzubilden und zu verfolgen. Eine kaskadiertes Inversionsschema wurde auf synthetische time-lapse Daten eines vereinfachten Ölfeldes angewendet. Untersucht wurde der Einfluss verschiedener Parameter insbesondere verschiedener Regularisierungstechniken. Sender und Empfänger an der Erdoberfläche sind typischerweise wenig sensitiv zu tiefen Strukturen in leitfähiger Umgebung. Anhand von synthetische Daten konnte gezeigt werden, dass das benutzen zusätzlicher Nebenbedingungen wie einer Modellgewichtung und dem ausnutzen von vorhandenen Stahlverrohrungen es dennoch erlaubt Änderungen innerhalb des Ölreservoirs zu lokalisieren.