Marine EM, the Past, the Present, and the Future. Steven Constable Scripps Institution of Oceanography

Similar documents
Ten years of marine CSEM for hydrocarbon exploration

Vertical and horizontal resolution considerations for a joint 3D CSEM and MT inversion

Special Section Marine Control-Source Electromagnetic Methods

Modeling of 3D MCSEM and Sensitivity Analysis

2011 SEG SEG San Antonio 2011 Annual Meeting 771. Summary. Method

Towed Streamer EM Integrated interpretation for accurate characterization of the sub-surface. PETEX, Tuesday 15th of November 2016

A Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods

Anisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements

Controlled-source electromagnetic imaging of the Middle America Trench o"shore Nicaragua

Towed Streamer EM data from Barents Sea, Norway

The feasibility of reservoir monitoring using time-lapse marine CSEM

D034 On the Effects of Anisotropy in Marine CSEM

Tu B3 15 Multi-physics Characterisation of Reservoir Prospects in the Hoop Area of the Barents Sea

Summary. Introduction

IPTC PP Challenges in Shallow Water CSEM Surveying: A Case History from Southeast Asia

Estimating vertical and horizontal resistivity of the overburden and the reservoir for the Alvheim Boa field. Folke Engelmark* and Johan Mattsson, PGS

G008 Advancing Marine Controlled Source Electromagnetics in the Santos Basin, Brazil

Repeatability and Detectability Requirements for 4D CSEM

Modeling and interpretation of CSEM data from Bressay, Bentley and Kraken area of East Shetland Platform, North Sea

Geophysical model response in a shale gas

National Geophysical Survey Science Opportunities

We A Multi-Measurement Integration Case Study from West Loppa Area in the Barents Sea

Seismic processing of numerical EM data John W. Neese* and Leon Thomsen, University of Houston

Using seismic guided EM inversion to explore a complex geological area: An application to the Kraken and Bressay heavy oil discoveries, North Sea

Downloaded 08/29/13 to Redistribution subject to SEG license or copyright; see Terms of Use at

PhD Thesis. Multiscale methods for CSEM data interpretation

Reducing Uncertainty through Multi-Measurement Integration: from Regional to Reservoir scale

X004 3D CSEM Inversion Strategy - An Example Offshore West of Greenland

Rock physics integration of CSEM and seismic data: a case study based on the Luva gas field.

Summary. Introduction

Improved Exploration, Appraisal and Production Monitoring with Multi-Transient EM Solutions

Summary. (a) (b) Introduction

Integrating seismic, CSEM, and well-log data for reservoir characterization

Downloaded 10/29/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

Identified a possible new offset location where the customer is currently exploring drill options.

Sub-basalt exploration in the Kutch-Saurashtra basin using EM

Downloaded 03/06/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

SEG Houston 2009 International Exposition and Annual Meeting

SEAFLOOR ELECTROMAGNETIC METHODS CONSORTIUM A research proposal submitted by

POTASH DRAGON CHILE GEOPHYSICAL SURVEY TRANSIENT ELECTROMAGNETIC (TEM) METHOD. LLAMARA and SOLIDA PROJECTS SALAR DE LLAMARA, IQUIQUE, REGION I, CHILE

The prediction of reservoir

SIO223A A Brief Introduction to Geophysical Modeling and Inverse Theory.

Electrical prospecting involves detection of surface effects produced by electrical current flow in the ground.

A Broadband marine CSEM demonstration survey to map the Uranus salt structure

SEG Houston 2009 International Exposition and Annual Meeting

Detection, Delineation and Characterization of Shallow Anomalies Using Dual Sensor Seismic and Towed Streamer EM data

Consider a Geonics EM31 with the TX and RX dipoles a distance s = 3 m apart.

Magnetotelluric (MT) Method

We E Estimation of Induction Noise in a Towed EM Streamer

An Introduction to Marine Controlled Source Electromagnetic Methods

Subsurface hydrocarbons detected by electromagnetic sounding

Geophysics Course Introduction to DC Resistivity

Pluto 1.5 2D ELASTIC MODEL FOR WAVEFIELD INVESTIGATIONS OF SUBSALT OBJECTIVES, DEEP WATER GULF OF MEXICO*

SPE DISTINGUISHED LECTURER SERIES is funded principally through a grant of the SPE FOUNDATION

RC 1.3. SEG/Houston 2005 Annual Meeting 1307

Use of Seismic and EM Data for Exploration, Appraisal and Reservoir Characterization

Subsurface Geology and Resource Exploration

6298 Stress induced azimuthally anisotropic reservoir - AVO modeling

Predicting Gas Hydrates Using Prestack Seismic Data in Deepwater Gulf of Mexico (JIP Projects)

Geoelectricity. ieso 2010

Hydrocarbon Volumetric Analysis Using Seismic and Borehole Data over Umoru Field, Niger Delta-Nigeria

Seismoelectric Ground-flow DC-4500 Locator

KMS Technologies KJT Enterprises, Inc. Novel Marine Electromagnetics: from Deep into Shallow Water

MT Prospecting. Map Resistivity. Determine Formations. Determine Structure. Targeted Drilling

Gas hydrate-related sedimentary pore pressure changes offshore Angola

The Afar Rift Consortium

Seismic wavepropagation concepts applied to the interpretation of marine controlled-source electromagnetics

Downloaded 11/20/12 to Redistribution subject to SEG license or copyright; see Terms of Use at

3. Magnetic Methods / 62

Acoustic Anisotropy Measurements and Interpretation in Deviated Wells

Geophysics for Environmental and Geotechnical Applications

Determination of reservoir properties from the integration of CSEM and seismic data

H005 Pre-salt Depth Imaging of the Deepwater Santos Basin, Brazil

Reservoir properties inversion from AVO attributes

PART A: Short-answer questions (50%; each worth 2%)

A multigrid integral equation method for large-scale models with inhomogeneous backgrounds

The elastic properties such as velocity, density, impedance,

Marine CSEM of the Scarborough gas field, Part 2: 2D inversion

New EM technology offerings are growing quickly

Application of sensitivity analysis in DC resistivity monitoring of SAGD steam chambers

SEAM2: 3D NON-SEISMIC MODELING OF A COMPLEX MIDDLE EAST O&G PROSPECT

Electrical Impedance Tomography of a Seafloor Volcano R.V. Revelle shiptime proposal submitted by David Myer, Steven Constable, and Kerry Key

Integrated Geophysical Model for Suswa Geothermal Prospect using Resistivity, Seismics and Gravity Survey Data in Kenya

E : Ground-penetrating radar (GPR)

GY 402: Sedimentary Petrology

Some consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston

So I have a Seismic Image, But what is in that Image?

SEG Houston 2009 International Exposition and Annual Meeting SUMMARY

A11 Planning Time-lapse CSEM-surveys for Joint Seismic-EM Monitoring of Geological Carbon Dioxide Injection

Maxwell s Equations & Hertz Waves

Magnetotelluric Array and Magnetic Magnetotelluric Surveys

Pressure Prediction and Hazard Avoidance through Improved Seismic Imaging

Seismic Prospecting For Oil By C. Hewitt Dix READ ONLINE

Detecting fractures using time-lapse 3C-3D seismic data

Interpretation and Reservoir Properties Estimation Using Dual-Sensor Streamer Seismic Without the Use of Well

Introducing the Stratagem EH4. Electrical Conductivity Imaging System Hybrid-Source Magnetotellurics

Integrating rock physics and full elastic modeling for reservoir characterization Mosab Nasser and John B. Sinton*, Maersk Oil Houston Inc.

DEEP SAN ANDREAS FAULT BOUNDARY STRUCTURE FROM MARINE MAGNETOTELLURICS A

Electrical Methods. Resistivity Surveying

Site Characterization & Hydrogeophysics

Transcription:

Marine EM, the Past, the Present, and the Future Steven Constable Scripps Institution of Oceanography

Acknowledgements SEG and SEG Foundation Sponsored by Statoil Sponsored by Paradigm Scripps Institution of Oceanography, Seafloor Electromagnetic Methods Consortium

SEG Membership Benefits SEG Digital Library full text articles Technical Journals in Print and Online Join Today! Networking Opportunities Membership Discounts on! Continuing Education Training Courses! Publications (35% off list price)! Workshops and Meetings Membership materials are available today! Join Online at seg.org/about-seg/membership

Student Opportunities Sponsored Membership Student Chapter Programs SEG/Chevron Student Leadership Symposium SEG/ExxonMobil Student Education Program Challenge Bowl Scholarships Field Camp Grants Geoscientists Without Borders Student Expos & IGSCs Honorary and Distinguished Lecturers SEG Online Learn more at seg.org/education/students-early-career

Section/Associated Society Opportunities Host DL, HL, and DISC Programs Council Representation Annual Meeting Booth Discount Best Papers presented at SEG Annual Meetings Joint Conferences, Workshops, and Forums in partnership with SEG For more information, including a list of benefits, please visit: www.seg.org/resources/sections-societies

Marine EM, the Past, the Present, and the Future Steven Constable Scripps Institution of Oceanography

On Tuesday 17th August, 2004 the front page of the Wall Street Journal had a lead article about a new geophysical technique. Marine CSEM was suddenly in the news

On Tuesday 17th August, 2004 the front page of the Wall Street Journal had a lead article about a new geophysical technique. Marine CSEM was suddenly in the news except is wasn t new. Over 23 years earlier an oceanographer called Charles Cox had suggested marine CSEM for hydrocarbon exploration

On Tuesday 17th August, 2004 the front page of the Wall Street Journal had a lead article about a new geophysical technique. Marine CSEM was suddenly in the news except is wasn t new. Over 23 years earlier an oceanographer called Charles Cox had proposed using marine CSEM for hydrocarbon exploration to Exxon.

Outline: Basics of the EM methods Early history of marine EM Commercialization 10 things you need to know The future

What is CSEM? Controlled Source ElectroMagnetic sounding has been used since the 1930 s to map sub-surface geology through the proxy of electrical conductivity. On land it is mainly used for mining exploration (shallower, conductive ore bodies), but has been used for map geological structure for oil exploration. Transmitters Receivers Bz I grounded wire (E) By Bx I horizontal loop (B) Ex V Ey V

What is MT? CSEM is related to another EM technique called MagnetoTelluric (MT) sounding, which uses similar receivers to measure Earth s natural magnetic field variations and the induced electric currents. MT signals are either part of a useful complementary method or a source of noise for CSEM. Natural-source Magnetotelluric fields Transmitters Receivers Bz I grounded wire (E) By Bx I horizontal loop (B) Ex V Ey V

Both can be used in the marine environment: Marine CSEM more sensitive to resistive rocks Marine MT more sensitive to conductive rocks But both methods can detect contrasts Air (resistive) Natural-source Magnetotelluric fields CSEM Transmitter Seawater (very conductive) 100-300 m dipole, 100-1000 amps 25-100 m Oil, Gas (resistive) Seafloor (variable conductivity) Electric and magnetic field recorders Salt, carbonate, volcanics (resistive)

Both MT and CSEM sounding use electromagnetic induction, which describes what happens around a time-varying primary magnetic field: primary magnetic field

Faraday s Law says that a time varying (or moving) magnetic field will induce electric fields: E dl = d C dt ( is magnetic flux). primary magnetic field electric field

Faraday s Law says that a time varying (or moving) magnetic field will induce electric fields: C E dl = d dt ( is magnetic flux). primary magnetic field electric field and current Ohm s Law says that E will generate a current J in a conductor: J = E

Faraday s Law says that a time varying (or moving) magnetic field will induce electric fields: C E dl = d dt ( is magnetic flux). Ohm s Law says that E will generate a current J in a conductor: J = E Ampere s Law says that the current will generate a secondary magnetic field: primary magnetic field secondary magnetic field electric field and current C B dl = µi

The secondary field opposes the changes in the primary field. The consequence of this is that conductive rocks absorb variations in EM fields more than resistive rocks. Depth > Half-space

The secondary field opposes the changes in the primary field. The consequence of this is that conductive rocks absorb variations in EM fields more than resistive rocks. This absorption is exponential: E(z) =E o e z/z s The rate of absorption is given by the skin depth, which depends on rock resistivity and period: Depth > High resistivity, long periods = large skin depths, greater penetration. Low resistivity, short periods = small skin depths, greater attenuation. Half-space

CRIPPS INSTITUTION OF OCEANOGRAPHY UCSD Secondary fields are different in amplitude, phase, and direction than the primary magnetic and electric fields. Electric fields can also be distorted galvanically, as in DC resistivity methods. Secondary fields Primary fields Induced magnetic fields Galvanic distortion of electric currents

MT and CSEM are methods to measure electrical conductivity. Conductivity varies over 5 orders of magnitude in common Earth materials, and provides the ideal mechanism for studying fluids, water, and geology, including hydrocarbons. Conductivity, S/m 10 1 0.1 0.01 10-3 10-4 water 20 C 0 C seawater fresh water sediments mantle marine lower mantle land 1400 C 1200 C peridotite 1100 C basalt, salt, carbonates basaltic melt 10% 3% 1% partial melt ice, hydrate hydrocarbons 0.1 1 10 Resistivity, Ωm 100 1,000 10,000

The Past

Interestingly, Cagniard proposed adaption to the marine environment in the 1953 paper that first presented the MT method. GEOPHYSICS, VOL 18, NO. 3 (JULY 1953), P. 605-635.

Charles (Chip) Cox and Jean Filloux of Scripps Institution of Oceanography made the first marine MT measurements in 1965, in 4300 m water offshore California. Charles Cox Tucson mag Jean Filloux seafloor mag Tucson mag seafloor E Cox, Filloux, and Larsen, 1968

Theory pre-dates practice, but even so marine CSEM is a young field. The earliest marine CSEM work was carried out by the British and US navies. This 1968 paper out of the US Navy Underwater Sound Lab appears to be the first proposal for marine CSEM as we now know it.

The same Chip Cox (1922-2015) developed practical deepwater CSEM in the early 1980 s, partly under funding by DARPA to improve understanding of seafloor communications, and ONR for submarine noise studies.

In 1979 Chip carried out an experiment in nearly 3,000m water with transmissions of 80 amps on an 800 m antenna. Frequencies of 0.25-2.25 Hz were detected 19 km away. His target was a mid-ocean ridge. Tow cable Transmitter Antenna Receiver Electrodes Young and Cox, 1981 Cox et al, 1981

Remarkably, within 2 years, in March 1981, Chip proposed the CSEM method for oil exploration to Exxon:

In the proposal, Chip explicitly described the direct detection of a hydrocarbon reservoir. Fields measured 2500 m from a 10 4 Am CSEM transmitter

I checked Chip s calculations and they are basically correct. However, the proposal was declined - Chip was too far ahead of the times. 10 2 Amplitude (nv/m) 10 1 10 0 10 1 10 1 10 0 10 1 10 2 Hz Fields measured 2500 m from a 10 4 Am CSEM transmitter

Meanwhile, in the late 1980 s Martin Sinha of Cambridge (later Southampton) developed a UK CSEM capability based on Scripps, but with an flown transmitter capable of working over mid-ocean ridges. (From Sinha et al, 1990)

1 km 2 km thin resistive layer overhang Current industry funding for marine EM started in 1995 in order to develop MT imaging of salt in the Gulf of Mexico. Clearly, the price of oil was not a factor (adjusted for inflation, 1998 was the lowest price in history). 4 km 6 km SW NW Resistive (4 Ohm-m) Conductive (0.3 Ohm-m) What drove interest in such an expensive, nonseismic method?

m 1k overhang thin resistive layer m km 2k 6 km 4 SW Resistive (4 Ohm-m) NW Conductive (0.3 Ohm-m) Current industry funding for marine EM started in 1995 in order to develop MT imaging of salt in the Gulf of Mexico. Clearly, the price of oil was not a factor (adjusted for inflation, 1998 was the lowest price in history). What drove interest in such an expensive, nonseismic method? First, the seismic method was having great difficultly imaging beneath salt

1 km 2 km thin resistive layer overhang Current industry funding for marine EM started in 1995 in order to develop MT imaging of salt in the Gulf of Mexico. Clearly, the price of oil was not a factor (adjusted for inflation, 1998 was the lowest price in history). 4 km 6 km SW NW Resistive (4 Ohm-m) Conductive (0.3 Ohm-m) What drove interest in such an expensive, nonseismic method? But probably the biggest factor was the development of the tension leg platform, which enabled production of oil and gas from deepwater prospects but deepwater drilling is very expensive.

180 160 Wells drilled Oil price, nominal $ Gulf of Mexico, wells >1,000 m water depth and oil price Wells; Oil price 140 120 100 80 60 First Scripps industry funding starts First Scripps industry funding stops Current industry funding starts 40 20 1970 1975 1980 1985 1990 1995 2000 2005 2010 2015 Year First oilfield CSEM surveys, Angola Commercial marine MT was offered starting in 1996. The first oilfield CSEM tests were done in late 2000 - early 2002, by Statoil and ExxonMobil. They both used Scripps receivers and the Southampton transmitter. Again, the driver was the high cost of drilling combined with limitations in the seismic method.

The seismic method has the problem that small gas saturations produce big velocity changes. Seawater/gas in 50% porous sandstone 1.6 Archie s with fizz-gas P-wave velocity, km/s 1.2 0.8 Velocity Fizz-gas 0.4 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Gas saturation

However, electrical resistivity does not change until the gas saturation gets large. This means that marine CSEM can be used to assess targets prior to drilling. Seawater/gas in 50% porous sandstone 10 4 1.6 Archie s with fizz-gas P-wave velocity, km/s 1.2 0.8 Velocity Fizz-gas Least resistivity contrast easily detected with EM Archie s 10 3 10 2 10 1 Resistivity, Ωm, log scale 0.4 0 0.1 HS lower bound 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 100 Gas saturation

Oil and gas reservoirs are too thin to be sensed with MT - now we need to use a man-made source of energy. Again, EM energy propagates best through the more resistive rocks. In this movie, the electric fields are absorbed more rapidly in seawater: Air Sea Sediment movie Kerry Key, SIO.

CRIPPS INSTITUTION OF OCEANOGRAPHY UCSD If we add a thin resistive layer to the model (an oil or gas reservoir?), then energy propagates even better through this layer than the surrounding sediments. Then we add a transmitter and receivers to actually collect these data. Air Sea Oil layer Sediment movie Kerry Key, SIO.

The effect can be up to an order of magnitude for a large oil or gas reservoir. This is what I call the canonical model, and is similar to the target chosen by Statoil for the first test of the marine CSEM method, and very similar to the model that Chip used in his proposal. In-line electric field, V/m/(Am) 10-8 10-10 10-12 10-14 E-field amplitudes larger over hydrocarbon reservoirs no oil, gas 1000 m, 0.3 Ωm 1000 m, 1 Ωm 100 m, 100 Ωm oil, gas 0 2 4 6 8 10 Source-receiver offset, km

Statoil s tests were made in November 2000 over the Girasoll oil field off Angola: 35 30 tow A 10-10 Distance North, (km) 25 20 Radial Electric Field, (V/m)/(source dipole moment) 10-11 10-12 10-13 10-14 Infinite Reservoir Oil Reservoir No Reservoir Terminating Reservoir 15 10 5 V approximate extent of reservoir 0 5 10 15 20 21 Distance East (km) 10-15 -6000-4000 -2000 0 2000 4000 6000 N Source-Receiver Offset, m Constable and Srnka, 2007 http://www.oilpro.com/

Similar studies were carried out by ExxonMobil in Jan 2002, also off Angola. These studies used 30 new instruments designed and built by Scripps for ExxonMobil, operated by AOA Geophysics. Distance (km) 9 6 3 0 0 receiver 21 Salt 3 seismic O/W contact 6 Well B 9 Distance (km) Over a known discovery Well A tow 2 12 15 Normalized field (V/(Am 2) ) 10-10 10-11 10-12 10-13 10-14 10-15 Prior to drilling noise level reservoir limits 0 3 6 9 Distance along line (km) Inline E Field, 0.25 Hz Receiver 21, Tow Line 2 2 min. Stacked Data With Reservoir Backgound 12 15 Distance North (km) 15 12 9 6 3 Salt Oligocene tow 3 Salt Miocene well receiver 21 Radial Electric Field, (V/(Am 2 )) 10-10 10-11 10-12 10-13 10-14 noise level Rotated Inline E Field Receiver 21, Tow Line 3 0.25 Hz Data Mio + Olig HC Filled Oligocene HC Filled Brine Filled 0 0 10-15 0 3 6 9 12 15 3 6 9 12 15 Distance (km) Distance East (km) Constable and Srnka, 2007

Then things got interesting... By the end of 2002 three contractors were formed to offer marine MT/CSEM as a commercial product. At least 3 main contractors offering marine EM data collection and interpretation services Several special-purpose ships have been built since then. Much based on early academic science and technology.

Since then more companies have joined the party

The Present 10 things you should know about marine CSEM that are not likely to change

1. EM does not map oil or gas, or even rock type. We cannot tell the difference between an evaporite, a basalt, or an oil sand based on resistivity alone. Resistivity of Rocks and Minerals (Data from Telford et al.) basalts/lavas limestones wet OIL SANDS clays sandstones, quartzites dry pyrrhotite graphite seawater evaporites groundwater intrusive igneous rocks wet dry quartz -5-4 -3-2 -1 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 Log 10 (Resistivity, Ωm) Thus, to manage false positives you need additional data - seismics, geology, magnetics, gravity, Integration with other data sets is perhaps the only safe way to use marine EM data.

2. The resolution of EM induction is between wave propagation and potential fields: High frequency (megahertz) Radar Wave equation: Resolution ~ wavelength 2 E = µ E t + µ 2 E t 2 Seismics r 2 u = @u @t + 1 @ 2 u c 2 @t 2 Mid frequency (0.001-1000 Hz) Inductive EM Diffusion equation: Resolution ~ size/depth 2 E = µ E t Zero frequency DC Resistivity Laplace equation: Resolution ~ bounds only 2 E =0 Gravity/ Magnetism r 2 U =0 µ = electrical conductivity ~ = magnetic permeability ~ 3 10 6 S/m 10 4 10 6 = electric permittivity ~ 10 9 10 11 F/m H/m

3. Frequency is too low for wave propagation in rocks, but does constrain length scale through skin depth : z s At a single frequency 2 E = µ E becomes t 2 E = i µ E =2 f ( ) which for uniform fields has solutions of the form where we have defined a skin depth E(z) =E o e z/z s+i( t z/z s ) z s = 2/ µ For every skin depth the fields decay by 1/e (~37%) and phase lags by one radian (~57 ). In 1 Ωm sediments: frequency skin depth 0.1Hz 1,600 m 0.25 Hz 1,000 m 1.0 Hz 500 m 5.0 Hz 225 m

Although you need a low enough frequency to reach your target, sensitivity can increase with frequency (Eradial, 3,000 m, phase): 7830 7825 7820 degrees 3000m, E rad 0.25Hz 40 50 degrees 140 3000m, E rad 0.75Hz 150 Northing, km 7815 7810 7805 7800 7795 60 0.25 Hz 70 0.75 Hz 80 160 170 180 7790 7785 N = 623 rms = 4 90 N = 642 rms = 3 190 7830 7825 3000m, E rad 1.75Hz degrees 280 degrees 3000m, E rad 3.25Hz 410 7820 7815 290 420 Northing, km 7810 7805 7800 300 1.75 Hz 310 3.25 Hz 430 440 7795 320 450 7790 N = 602 rms = 5 7785 710 720 730 740 750 760 Easting, km 330 N = 348 rms = 8 710 720 730 740 750 760 Easting, km 460 (But your signal size gets smaller.) These are data over the Scarborough gas field offshore NW Australia, collected by Scripps with support from BHP Billiton.

Multiple frequencies Single frequency - 0.25 Hz Single frequency - 0.75 Hz Single frequency - 1.75 Hz 1 y transmitter Constable, Orange, and Key, 2015 Depth (km) 1.5 2 2.5 3 4. Good inversion technology is key to getting the most out of EM data, and inversions of multiple frequencies are much better than single frequency inversions. 3.5 4 Truth 1 Hz 0.1, 1Hz 0.1, 0.3, 1, 3, 10 Hz 4.5 10 1 10 0 10 1 10 2 Resistivity (ohm m) Key, 2009

With a single frequency and no phase data, a shallow resistor (say, hydrate) has an almost identical response to a deep hydrocarbon reservoir. This led to actual drilling errors. Reservoir Ey Shallow slab Receiver Reservoir 0.33 ohm-m Shallow slab 0.33 ohm-m 1 ohm-m 1 ohm-m Courtesy Arnold Orange

Northing, km Northing, km 7830 7825 7820 7815 7810 7805 7800 7795 7790 7785 7830 7825 7820 7815 7810 7805 7800 7795 7790 o 0 Azimuth Field purely radial (in-line geometry) Azimuth 4000m, E rad 0.75Hz, log 10 (V/Am 2 ) N = 690 rms = 10% 13.25 4000m, E az 0.75Hz N = 297 rms = 9% 7785 710 720 730 740 750 760 Easting, km log 10 (V/Am 2 ) Receiver 200 5. Transmitter fields have a dipole geometry. The radial (inline) and azimuthal 210 (broadside) fields behave differently. The greater vertical component 220 of the radial fields creates greater sensitivity to thin, 230 sub-horizontal resistors (reservoirs). You can use 12.8 this to tell thick from thin 240 resistors. The modes also 12.85 behave differently to anisotropy. 12.9 12.8 4000m, E rad 0.75Hz 12.85 12.9 12.95 13 13.05 13.1 13.15 13.2 13.2 13.25 13.3 13.35 13.4 13.45 13.5 13.55 13.6 13.65 N = 690 rms = 4 Radial field Er Eø Azimuthal field N = 297 rms = 6 7820 7810 710 720 730 740 750 760 Easting, km 7805 Northing, km degrees 250 At intermediate degrees azimuths the 100 4000m, E az 0.75Hz data are a trigonometric mixture of the 110 radial and azimuthal fields. 130 7830 7825 140 150 7815 12.95 13.05 13.1 13.15 13.2 120 N = 690 rms = 10% 13.25 0.75 Hz amplitude, 4km Tx-Rx offset 4000m, E az 0.75Hz 13 log 10 (V/Am 2 ) 13.2 13.25 13.3 13.35 13.4 13.45 N = Transmitter o 90 Azimuth Field purely azimuthal (broadside geometry) 7800 7795 7790 N = 297 rms = 9% 7785 710 720 730 740 750 760 Easting, km 13.5 13.55 13.6 13.65 710 N =

Inverting either inline Bx or Ey provides best resolution to reservoir targets. Note Ez: don t confuse sensitivity with resolution. x Crossline Tx x Inline Tx x Vertical Tx E B Bx B B E E Ez Ey y y y 1 1.5 x transmitter 1 1.5 y transmitter 1 1.5 z transmitter Depth (km) 2 2.5 3 3.5 4 Truth Ex By Bz Ex,By,Bz Depth (km) 2 2.5 3 3.5 4 Truth Bx Ey Ez Bx,Ey,Ez Depth (km) 2 2.5 3 3.5 4 Truth Bx Ey Ez Bx,Ey,Ez 4.5 10 1 10 0 10 1 10 2 Resistivity (ohm m) 4.5 10 1 10 0 10 1 10 2 Resistivity (ohm m) 4.5 10 1 10 0 10 1 10 2 Resistivity (ohm m) 1D Inversions of synthetic 0.1 and 1.0 Hz data with 1% noise added: Key, 2009

CRIPPS INSTITUTION OF OCEANOGRAPHY UCSD 6. CSEM works best in deep water, where the air wave is small. But, it works OK in shallow water because the air wave and target response are coupled. Don t try to remove the air wave. 'Air wave' Seawater propagation host rock propagation target signal Flood From: Orange, Key and Constable, Geophysics, 2009 Orange, Key, and Constable, 2009

7. Depth sensitivity in CSEM is largely determined by geometry. Sensitivity is always greatest right below the Tx and Rx. Depth sensitivity is about half Rx-Tx spacing, and so is sensitivity to off-line structure. Plan 0.25 Hz and perturbation of 1 Ωm halfspace Section inline (m) Constable, 2010

Targets have to be bigger than their depth of burial to be detected (note that this has nothing to do with transmitter power!). -8 Radial mode fields (1 Hz) (3D calculation courtesy C. Weiss) -9 log10(e field magnitude, V/m/(Am)) -10-11 -12-13 -14-15 1 km 2 km 3 km 1000 m, 1 Ωm 4 km 5 km diameter 100 m, 100 Ωm 1D layer -16-17 -18 1D halfspace 1000 2000 3000 4000 5000 6000 7000 8000 Range, m Constable and Weiss, 2006

8. Anisotropy is important. Aligned crystals and inter-bedding can produce anisotropy on microscopic to macroscopic scales. CSEM cannot tell micro- from macro-anisotropy until the layers get very thick. NR O Brien, J. Sed. Res. 1987 http://www.sleepingdogstudios.com/

To a good approximation (10% or so), the anisotropic response follows the vertical resistivity for the radial mode, and the horizontal resistivity for the azimuthal mode. Inverting the modes separately using isotropic models may give reasonable results, but joint inversion of both modes will fail without anisotropy. 1 Ωm 0.5 Ωm Red - magnetic fields Blue - electric fields Constable, 2010 Anisotropic Horizontal Vertical -8-8 Amplitude (log 10 V/m/(Am), T/(Am)) -9-10 -11-12 -13-14 -15-16 Radial Bx: Anisotropy Ey: Anisotropy Bx: HS_half Ey: HS_half Bx: HS_one Ey: HS_one Amplitude (log 10 V/m/(Am), T/(Am)) -9-10 -11-12 -13-14 -15-16 Azimuthal By: Anisotropy Ex: Anisotropy By: HS_half Ex: HS_half By: HS_one Ex: HS_one -17 RADIAL: Frequency: 0.25 Hz -18 0 1 2 3 4 5 6 7 8 9 10 Horizontal Range (km) -17 AZIMUTHAL: Frequency: 0.25 Hz -18 0 1 2 3 4 5 6 7 8 9 10 Horizontal Range (km)

9. Seawater conductivity matters. Getting it wrong can introduce spurious structure into CSEM models (red = correct, black = 1-layer approximation). Expensive mistakes have been made... 0 0.5 1 Truth Free sea, 1 layer, rms 3.6 Fixed sea, 20 layers 1.5 Depth (km) 2 2.5 3 3.5 4 4.5 10 1 10 0 10 1 10 2 Resistivit y (ohm m) Key, 2009

10. Like DC resistivity, CSEM is mainly sensitive to the resistivity-thickness product of thin resistive layers. This is called T-equivalence, where T = transverse resistance (i.e. resistivity times thickness). 10 1 I V 100 m, 1 Ωm 10 m, 100 Ωm I V 100 m, 1 Ωm 1 m, 1000 Ωm Apparent Resistivity 10 0 10 1 10 2 10 3 10 4 Electrode Spacing If resistivity is proportional to net pay, this may not be much of a problem.

The Future

Marine EM was almost certainly over-sold between 2002 and 2007 or so. WesternGeco pulled out of the marine EM business in 2010 and OHM was absorbed into EMGS in 2012. EM contractor share price The marine EM market was recently worth around $200m/year, or about 5% of the marine seismic market. Marine EM is a capital-intensive business and there is a persistent fear of patent lawsuits, so it is difficult for smaller, innovative companies to enter the market. On the other hand, there probably isn t enough profit for the big companies to care. But there is some progress.

PGS surface towed streamer system is providing good data in water depths up to 400 m or so, and allows data collection co-incident with 2D seismics. PGS website

Petromarker has a novel vertical-vertical system which can operate in deep water, but has limited lateral capabilities. Petromarker website

On the other hand, there is only one contractor offering the standard approach that works so well in deep water. EMGS website However, the custom vessel, large (3D) survey model has made CSEM a relatively expensive commodity.

Does it work? Jonny Hesthammer thinks so. False positives have been drilled (mostly lithological resistors), but false negatives are very rare. That is, without a CSEM signature the likelihood of a commercial discovery is very low - but industry just keeps drilling and providing more data for Jonny s plot - without more confidence in the method it is difficult to realize its value. Strong EM response (NAR>15%) Weak EM response (NAR<15%) Increasing EM signature 4.5 350% 300% 4 250% 3.5 200% 3 150% 2.5 100% 2 50% 1.5 01-50% 0.5 Ivory 32 mmboe Goliat 238 mmboe Havis 250 mmboe Wisting 132 mmboe Skrugard 286 mmboe Norvarg 189 mmboe Skavl 35 mmboe Hanssen 18-56 mmboe Storbarden - Dry Salina 38 mmboe Eik - Dry Mercury 6-12 mmboe Atlantis Minor gas Apollo - Dry* Phoenix Dry* Breiflabb Heilo - Dry 1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47 49 51 53 55 57 59 61 63 65 67 69 71 73 75 77 79 81 83 85 Mostly discoveries (76%) Mostly dry wells (68%) (2.9% commercial discovery rate) Discovery - Exploration Discovery - Calibration Discovery Unknown commerciality Discovery Non-commercial Dry well - Exploration Dry well - Calibration Wisting Alt. - Dry Darwin - Dry Byrkje - Dry Bønna - Dry Recent CSEM wells NAR=15% Hesthammer et al., First Break, 2012; Hesthammer pers. comm. 2015

Much of the CSEM data is collected by the majors (ExxonMobil, Shell, Chevron, etc.) and the large NOCs (Petrobras, Pemex, etc.). These companies keep their success rate statistics fairly confidential. But they are still collecting data I think the future lies in smaller, portable, ship-of-opportunity systems (like the ones academics must use), expanding the market to smaller client companies. This may take the passing of a few years and a rise in oil prices, but I think we can count on both of these. Thank you!

Connecting the world of applied geophysics