CA2468896C - Electromagnetic surveying for hydrocarbon reservoirs - Google Patents
Electromagnetic surveying for hydrocarbon reservoirs Download PDFInfo
- Publication number
- CA2468896C CA2468896C CA2468896A CA2468896A CA2468896C CA 2468896 C CA2468896 C CA 2468896C CA 2468896 A CA2468896 A CA 2468896A CA 2468896 A CA2468896 A CA 2468896A CA 2468896 C CA2468896 C CA 2468896C
- Authority
- CA
- Canada
- Prior art keywords
- survey
- electromagnetic
- source
- detector
- data sets
- 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
- 229930195733 hydrocarbon Natural products 0.000 title claims abstract description 91
- 150000002430 hydrocarbons Chemical class 0.000 title claims abstract description 90
- 239000004215 Carbon black (E152) Substances 0.000 title claims abstract description 89
- 238000000034 method Methods 0.000 claims abstract description 79
- 230000000694 effects Effects 0.000 claims abstract description 33
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 25
- 239000011435 rock Substances 0.000 claims abstract description 21
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 15
- 230000001939 inductive effect Effects 0.000 claims abstract description 14
- 230000005684 electric field Effects 0.000 claims description 28
- 230000006870 function Effects 0.000 claims description 22
- 238000010606 normalization Methods 0.000 claims description 12
- 230000005672 electromagnetic field Effects 0.000 claims description 9
- 238000004088 simulation Methods 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 4
- 230000001131 transforming effect Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 4
- 230000001419 dependent effect Effects 0.000 abstract description 3
- 239000013535 sea water Substances 0.000 description 13
- 230000001965 increasing effect Effects 0.000 description 11
- 238000013459 approach Methods 0.000 description 8
- 239000011148 porous material Substances 0.000 description 7
- 238000000926 separation method Methods 0.000 description 6
- 239000012530 fluid Substances 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 239000013049 sediment Substances 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 239000008186 active pharmaceutical agent Substances 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 238000011109 contamination Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000005358 geomagnetic field Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 101710083262 Ectin Proteins 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000007340 echolocation Effects 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005188 flotation Methods 0.000 description 1
- 239000013505 freshwater Substances 0.000 description 1
- 125000001183 hydrocarbyl group Chemical group 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
An electromagnetic survey method for surveying an area previously identified as potentially containing a subsea hydrocarbon reservoir, comprising obtaining first and second survey data sets with an electromagnetic source aligned end-on and broadside relative to the same or different receivers. The invention also relates to planning a survey using this method, and to analysis of survey data taken in combination allow the galvanic contribution to the signals collected at the receiver to be contrasted with the inductive effects, and the effects of signal attenuation, which are highly dependent on local properties of the rock formation, overlying water and air at the survey area. This is very important to the success of using electromagnetic surveying for identifying hydrocarbon reserves and distinguishing them from other classes of structure.
Description
TITLE OF THE INVENTION
ELECTROMAGNETIC SURVEYING FOR HYDROCARBON RESERVOIRS
BACKGROUND OF THE INVENTION
The invention relates to seafloor electromagnetic surveying for oil and other hydrocarbon reserves.
Determining the response of the sub-surface strata within the earth's crust to electromagnetic fields is a valuable tool in the field of geophysical research. The geological processes occurring in thermally, hydrothermally or magmatically active regions can be studied. In addition, electromagnetic sounding techniques can provide valuable insights into the nature, and particularly the likely hydrocarbon content, of subterranean reservoirs in the context of subterranean oil exploration and surveying.
Seismic techniques are often used during oil-exploration expeditions to identify the existence, location and extent of reservoirs in subterranean rock strata.
Whilst seismic surveying is able to identify such structures, the technique is often unable to distinguish between the different possible compositions of pore fluids within them, especially for pore fluids which have similar mechanical properties. In the field of oil exploration, it is necessary to determine whether a previously identified reservoir contains oil or just aqueous pore fluids. To do this, an exploratory well is drilled to determine the contents of the reservoir. However, this is an expensive process, and one which provides no guarantee of reward.
Whilst oil-filled and water-filled reservoirs are mechanically similar, they do possess significantly different electrical properties and these provide for the possibility of electromagnetic based discrimination testing. A known technique for electromagnetic probing of subterranean rock strata is the passive magneto-telluric (MT) method. The signal measured by a surface-based electromagnetic detector in response to electromagnetic (EM) fields generated naturally, such as within the earth's upper atmosphere, can provide details about the surrounding subterranean rock strata.
However, for deep-sea surveys, all but those MT signals with periods corresponding to several cycles per hour are screened from the seafloor by the highly conductive seawater. Whilst the long wavelength signals which do penetrate to the seafloor can be used for large scale undersea probing, they do not provide sufficient spatial resolution to examine the electrical properties of the typically relatively small scale subterranean reservoirs. Moreover, since MT surveying relies primarily on horizontally polarised EM fields, it is intrinsically insensitive to thin resistive layers.
Nonetheless, measurements of electrical resistivity beneath the seafloor have traditionally played a crucial role in hydrocarbon exploration and reservoir assessment and development. In industry, sub-seafloor resistivity data have generally been obtained almost exclusively by wire-line logging of wells. There are, though, clear advantages to developing non-invasive geophysical methods capable of providing such information. Although inevitably such methods would be unable to provide comparable vertical resolution to wireline logging, the vast saving in terms of avoiding the costs of drilling test wells into structures that do not contain economically recoverable amounts of hydrocarbon would represent a major economic advantage.
In research fields that are not of commercial interest, geophysical methods for mapping sub-seafloor resistivity variations by various forms of electromagnetic surveying have been under development for many years (e.g. Sinha et al., 1990;
Evans et al., 1994). WO 00/13046 and WO 00/57555 make proposals for finding hydrocarbon reservoirs using such electromagnetic surveying.
ELECTROMAGNETIC SURVEYING FOR HYDROCARBON RESERVOIRS
BACKGROUND OF THE INVENTION
The invention relates to seafloor electromagnetic surveying for oil and other hydrocarbon reserves.
Determining the response of the sub-surface strata within the earth's crust to electromagnetic fields is a valuable tool in the field of geophysical research. The geological processes occurring in thermally, hydrothermally or magmatically active regions can be studied. In addition, electromagnetic sounding techniques can provide valuable insights into the nature, and particularly the likely hydrocarbon content, of subterranean reservoirs in the context of subterranean oil exploration and surveying.
Seismic techniques are often used during oil-exploration expeditions to identify the existence, location and extent of reservoirs in subterranean rock strata.
Whilst seismic surveying is able to identify such structures, the technique is often unable to distinguish between the different possible compositions of pore fluids within them, especially for pore fluids which have similar mechanical properties. In the field of oil exploration, it is necessary to determine whether a previously identified reservoir contains oil or just aqueous pore fluids. To do this, an exploratory well is drilled to determine the contents of the reservoir. However, this is an expensive process, and one which provides no guarantee of reward.
Whilst oil-filled and water-filled reservoirs are mechanically similar, they do possess significantly different electrical properties and these provide for the possibility of electromagnetic based discrimination testing. A known technique for electromagnetic probing of subterranean rock strata is the passive magneto-telluric (MT) method. The signal measured by a surface-based electromagnetic detector in response to electromagnetic (EM) fields generated naturally, such as within the earth's upper atmosphere, can provide details about the surrounding subterranean rock strata.
However, for deep-sea surveys, all but those MT signals with periods corresponding to several cycles per hour are screened from the seafloor by the highly conductive seawater. Whilst the long wavelength signals which do penetrate to the seafloor can be used for large scale undersea probing, they do not provide sufficient spatial resolution to examine the electrical properties of the typically relatively small scale subterranean reservoirs. Moreover, since MT surveying relies primarily on horizontally polarised EM fields, it is intrinsically insensitive to thin resistive layers.
Nonetheless, measurements of electrical resistivity beneath the seafloor have traditionally played a crucial role in hydrocarbon exploration and reservoir assessment and development. In industry, sub-seafloor resistivity data have generally been obtained almost exclusively by wire-line logging of wells. There are, though, clear advantages to developing non-invasive geophysical methods capable of providing such information. Although inevitably such methods would be unable to provide comparable vertical resolution to wireline logging, the vast saving in terms of avoiding the costs of drilling test wells into structures that do not contain economically recoverable amounts of hydrocarbon would represent a major economic advantage.
In research fields that are not of commercial interest, geophysical methods for mapping sub-seafloor resistivity variations by various forms of electromagnetic surveying have been under development for many years (e.g. Sinha et al., 1990;
Evans et al., 1994). WO 00/13046 and WO 00/57555 make proposals for finding hydrocarbon reservoirs using such electromagnetic surveying.
SUMMARY OF THE INVENTION
The invention discloses a new approach for electromagnetic surveying to locate hydrocarbon layers. New source-detector geometries are used based upon an electromagnetic source.
According to a first aspect of the invention there is provided an electromagnetic survey method for surveying an area previously identified as potentially containing a subsea hydrocarbon reservoir, comprising: obtaining first and second survey data sets with an electromagnetic source aligned end-on relative to a first detector and aligned broadside relative to a second detector.
The terms source and detector are used interchangeably with transmitter and receiver respectively throughout this document.
The survey data from end-on and broadside alignments taken in combination allow the difference between galvanic and inductive contributions to the signals collected at the detector to be determined. Collecting survey data highlighting only the galvanic contribution in an end-on geometry is not reliable. As is demonstrated in the examples given below, it is generally impossible to differentiate between a rock formation containing a hydrocarbon reservoir and one which does not contain a hydrocarbon reservoir by studying the end-on survey data alone. The previously proposed electromagnetic survey methods for finding hydrocarbon reservoirs are thus believed to be at best highly unreliable.
The end-on survey data are sensitive to the presence of resistive hydrocarbon layers (exploiting largely galvanic effects and vertical components of induced current flow). By contrast, the broadside survey data are sensitive to the larger scale structure, but relatively insensitive to resistive hydrocarbon layers (exploiting the dominantly inductive effects). The reason why collection of survey data from the broadside geometry is essential for reliable electromagnetic surveying is that many features other than hydrocarbon reservoirs can affect the resistivity beneath the seafloor and the results of a survey. For example, resistivity often increases steadily with depth in submarine sedimentary basins, due to the progressive expulsion of pore fluids by rising overburden pressure. Such a resistivity profile will produce effectively the same response in the end-on survey data, as would presence of a hydrocarbon layer.
Comparative examples given below demonstrate this effect.
In the preferred implementation of the survey of the first aspect of the invention, the end-on and broadside alignments correspond to an arrangement of the electromagnetic source and the first and second detectors in which a right angle is formed between a first line leading from the first detector to the source and a second line leading from the second detector to the source, and wherein the source has its dipole axis aligned along the first line. However, in practice, an only approximate satisfaction of this condition will not greatly reduce the quality of the collected survey data. In any case, this ideal condition will not be satisfied exactly in practice, since the source is typically moved during surveying, being in the form of an antenna towed by an unmanned submarine craft. It would also be possible to obtain useable data if the above-mentioned right angle was changed to an angle away from 90 degrees, for example anywhere from 45-135 degrees may be satisfactory. How much the quality of the data deteriorates as the survey geometry moves away from right angles has not been studied, although this would be straightforward to do using the modelling techniques described herein.
With the first aspect of the invention, the first and second survey data sets are preferably obtained concurrently. This can be achieved during a single tow of the electromagnetic source. The data collected by the detectors can then be time-synchronised to the same absolute clock.
According to a second aspect of the invention there is provided an electromagnetic survey method for surveying an area previously identified as potentially containing a subsea hydrocarbon reservoir, comprising: obtaining first and second survey data sets with an electromagnetic source aligned end-on and broadside respectively relative to a first detector.
The invention discloses a new approach for electromagnetic surveying to locate hydrocarbon layers. New source-detector geometries are used based upon an electromagnetic source.
According to a first aspect of the invention there is provided an electromagnetic survey method for surveying an area previously identified as potentially containing a subsea hydrocarbon reservoir, comprising: obtaining first and second survey data sets with an electromagnetic source aligned end-on relative to a first detector and aligned broadside relative to a second detector.
The terms source and detector are used interchangeably with transmitter and receiver respectively throughout this document.
The survey data from end-on and broadside alignments taken in combination allow the difference between galvanic and inductive contributions to the signals collected at the detector to be determined. Collecting survey data highlighting only the galvanic contribution in an end-on geometry is not reliable. As is demonstrated in the examples given below, it is generally impossible to differentiate between a rock formation containing a hydrocarbon reservoir and one which does not contain a hydrocarbon reservoir by studying the end-on survey data alone. The previously proposed electromagnetic survey methods for finding hydrocarbon reservoirs are thus believed to be at best highly unreliable.
The end-on survey data are sensitive to the presence of resistive hydrocarbon layers (exploiting largely galvanic effects and vertical components of induced current flow). By contrast, the broadside survey data are sensitive to the larger scale structure, but relatively insensitive to resistive hydrocarbon layers (exploiting the dominantly inductive effects). The reason why collection of survey data from the broadside geometry is essential for reliable electromagnetic surveying is that many features other than hydrocarbon reservoirs can affect the resistivity beneath the seafloor and the results of a survey. For example, resistivity often increases steadily with depth in submarine sedimentary basins, due to the progressive expulsion of pore fluids by rising overburden pressure. Such a resistivity profile will produce effectively the same response in the end-on survey data, as would presence of a hydrocarbon layer.
Comparative examples given below demonstrate this effect.
In the preferred implementation of the survey of the first aspect of the invention, the end-on and broadside alignments correspond to an arrangement of the electromagnetic source and the first and second detectors in which a right angle is formed between a first line leading from the first detector to the source and a second line leading from the second detector to the source, and wherein the source has its dipole axis aligned along the first line. However, in practice, an only approximate satisfaction of this condition will not greatly reduce the quality of the collected survey data. In any case, this ideal condition will not be satisfied exactly in practice, since the source is typically moved during surveying, being in the form of an antenna towed by an unmanned submarine craft. It would also be possible to obtain useable data if the above-mentioned right angle was changed to an angle away from 90 degrees, for example anywhere from 45-135 degrees may be satisfactory. How much the quality of the data deteriorates as the survey geometry moves away from right angles has not been studied, although this would be straightforward to do using the modelling techniques described herein.
With the first aspect of the invention, the first and second survey data sets are preferably obtained concurrently. This can be achieved during a single tow of the electromagnetic source. The data collected by the detectors can then be time-synchronised to the same absolute clock.
According to a second aspect of the invention there is provided an electromagnetic survey method for surveying an area previously identified as potentially containing a subsea hydrocarbon reservoir, comprising: obtaining first and second survey data sets with an electromagnetic source aligned end-on and broadside respectively relative to a first detector.
The method of the second aspect of the invention thus differs from that of the first aspect in that a single detector can be used to collect both the end-on and broadside survey data. This can be done by towing the source twice, once in a direction along the line connecting the source to the detector, and again in a direction transverse thereto. Namely, the first and second survey data sets can be obtained consecutively. The relative alignment between the source and detector when the end-on and broadside survey data are collected can be varied in the same way as discussed above in relation to the first aspect of the invention. In other words, it is best if the dipole is aligned along a line connecting the source and detector when the end-on survey data is collected, and aligned perpendicular to that line when the broadside survey data is collected. However, deviations from that condition will occur and can be tolerated.
It will also be understood that the first survey data set can be obtained before or after second data survey set using the same electromagnetic source.
In principle, the first and second survey data sets could be obtained with separate sources carried by different transmitters, in which case the data sets could be obtained concurrently. However, in practice, it is likely that only a single source will be deployed and the first and second data sets will thus be obtained one after the other.
According to a third aspect of the invention there is provided a method of analysing results from an electromagnetic survey of an area potentially containing a subsea hydrocarbon reservoir, comprising: providing first and second survey data sets obtained from an electromagnetic source aligned respectively end-on and broadside relative to a detector; and combining the first and second survey data sets to obtain a results data set that represents a difference between the end-on and broadside alignments.
The method can be greatly improved by normalising each of the first and second survey data sets relative to respective first and second normalisation data sets or functions specific to the end-on and broadside alignments respectively, prior to the combining.
It will also be understood that the first survey data set can be obtained before or after second data survey set using the same electromagnetic source.
In principle, the first and second survey data sets could be obtained with separate sources carried by different transmitters, in which case the data sets could be obtained concurrently. However, in practice, it is likely that only a single source will be deployed and the first and second data sets will thus be obtained one after the other.
According to a third aspect of the invention there is provided a method of analysing results from an electromagnetic survey of an area potentially containing a subsea hydrocarbon reservoir, comprising: providing first and second survey data sets obtained from an electromagnetic source aligned respectively end-on and broadside relative to a detector; and combining the first and second survey data sets to obtain a results data set that represents a difference between the end-on and broadside alignments.
The method can be greatly improved by normalising each of the first and second survey data sets relative to respective first and second normalisation data sets or functions specific to the end-on and broadside alignments respectively, prior to the combining.
The first and second normalisation data sets or functions can be calculated from a rock formation model, or from the first and second survey data sets.
In a preferred embodiment, the first and second data sets each comprise radial and azimuthal components of electric field or magnetic field measured at the detector, and the method further comprises: transforming the radial and azimuthal components into at least one polarisation ellipse parameter, prior to the combining. The polarisation ellipse parameters) can be the amplitude and/or phase of the component of the electric field or magnetic field aligned along a major axis of the ellipse.
The method may advantageously further comprise: visually representing the results data set in a plot of at least two dimensions corresponding to the survey area.
The visual representation can be a two-dimensional (2D) plot in plan view, or a three-dimensional (3D) plot including depth, e.g. a perspective view of the survey area.
The plot may include markings of areas of equal or similar electromagnetic field strength. These markings could be contour lines, or colour or grey scale gradations with one colour or grey tone being used for a range of electromagnetic field strength values. Stepwise colour gradation is used in the preferred embodiment, with each colour representing a defined range of data values.
When the plot is of normalised survey data it is helpful if the plot also includes lines of equal absolute electromagnetic field strength, showing how the signal strength has decayed as one moves away from the source. The lines, which may appear as contours, can be labelled with either a relative or absolute decay value.
According to a fourth aspect of the invention, there is provided a method of planning an electromagnetic survey of an area identified as potentially containing a subsea hydrocarbon reservoir, comprising: creating a model of the area to be surveyed, including a rock formation containing a hydrocarbon reservoir and a body of water above the rock formation; setting values for water depth, depth below the seafloor of the hydrocarbon reservoir, and resistivity structure of the rock formation;
performing a simulation of an electromagnetic survey in the model of the survey area by calculating first and second survey data sets for an electromagnetic source aligned end-on and broadside relative to a detector; and combining the first and second survey data sets to obtain a results data set that represents a difference between the end-on and broadside alignments.
In typical use, the simulation will be repeated for a number of source frequencies and/or source-to-detector distances, and/or other parameters which can be varied during a survey. This iterative procedure can be used in order to select optimum surveying conditions in terms of source frequency and source-to-detector distance for probing the hydrocarbon reservoir. By optimum, it is not necessarily meant that the best survey conditions are found, but only that a set of survey conditions is found which will provide strong, unmistakable indications in the case that there is a hydrocarbon reservoir at the survey site. The iterative procedure may be purely under manual control. However, preferably, the simulator can allow the user the option of automatically optimising the survey conditions. The user can then switch between manual and automated iteration as desired.
The model should preferably include a body of air above the body of water, so that the simulation can take account of signal propagation paths including the body of air when calculating the first and second survey data sets. The propagation path through the air (the 'air wave') will in fact dominate for shallower water and longer distances between source and detector(s), as will be apparent from the examples described fiuther below. For deep water and shorter detector-source distances the effect is less important and may be omitted from the model.
The method preferably further comprises: normalising each of the first and second survey data sets relative to respective first and second normalisation data sets or functions specific to the end-on and broadside alignments respectively, prior to the combining. The first and second normalisation data sets or functions can be calculated from a rock formation model, for example.
The method preferably further comprises visually representing the results data set in a plot of at least two dimensions corresponding to the survey area. The visual representation can be a 2D plot in plan view, or a 3D plot including depth, e.g. a _$_ perspective view of the survey area. The other comments made above in relation to the plots of the third aspect of the invention also apply to the fourth aspect of the invention.
Another aspect of the invention relates to a computer program product bearing machine readable instructions for implementing the method of analysing results from an electromagnetic survey according to the third aspect of the invention.
Another aspect of the invention relates to a computer apparatus loaded with machine readable instructions for implementing the method of analysing results from an electromagnetic survey according to the third aspect of the invention.
Another aspect of the invention relates to a computer program product bearing machine readable instructions for implementing the method of planning an electromagnetic survey according to the fourth aspect of the invention.
Another aspect of the invention relates to a computer apparatus loaded with machine readable instructions for implementing the method of planning an electromagnetic survey according to the fourth aspect of the invention.
BRIEF DESCRIPTION OF THE DRAWINGS
For a better understanding of the invention and to show how the same may be carried into effect reference is now made by way of example to the accompanying drawings.
Figure 1 shows a schematic representation of an exploratory EM sounding survey.
Figure 2A is a graph showing the modelled electric field strength as a function of distance from, and along the axis of, the EM source. This is shown for several frequency components. The modelled sea depth is 750 m.
Figure 2B is a graph showing the modelled electric field strength as a function of distance from, and along the axis of, the EM source. This is shown for several frequency components. The modelled sea depth is 900 m.
Figure 2C is a graph showing the modelled electric field strength as a function of distance from, and along the axis of, the EM source. This is shown for several frequency components. The modelled sea depth is 1300 m.
Figure 3 is a schematic plan view defining a survey geometry coordinate system.
Figure 4A shows a schematic section view of a model subterranean strata configuration containing a hydrocarbon reservoir.
Figure 4B shows a schematic section view of a model subterranean strata configuration containing layers of increasing resistivity.
Figure 4C is a graph showing the electric field strength (Eo, E9o) as a function of distance from the source for two different detector positions for the subterranean strata shown in Figure 4A. Also shown are the normalising signals (No, N9o) expected from a uniform subterranean strata configuration.
Figure 4D is a graph showing the electric field strength (Eo, E9o) as a function of distance from the source for two different detector positions for the subterranean strata shown in Figure 4B. Also shown are the normalising signals (No, N9o) expected from a uniform subterranean strata configuration.
Figure 4E is a graph showing the electric field strengths (Eo, E9o) shown in Figure 4C normalised by the normalising field strengths (No, N9o) shown in the same.
Figure 4F is a graph showing the electric field strengths (Eo, E9o) shown in Figure 4D normalised by the normalising field strengths (No, N9o) shown in the same.
Figure 5 is a schematic plan view showing an arrangement of sixteen detectors on a section of sea floor above a subterranean reservoir, a particular source transmitter tow path is also shown.
Figure 6A shows a schematic section view of a model subterranean strata configuration containing a hydrocarbon reservoir.
Figure 6B is a plan view representing the normalised electric field strength as a function of location for the model subterranean strata configuration shown in Figure 6A for a first transmitter position.
Figure 6C is a plan view representing the normalised electric field strength as a function of location for the model subterranean strata configuration shown in Figure 6A for a second transmitter position.
DETAILED DESCRIPTION
Figure 1 of the accompanying drawing shows schematically a surface vessel 14 undertaking EM sounding of subterranean rock strata 8 within which a hydrocarbon reservoir 12 is located. The surface vessel 14 floats on the surface 2 of the sea 4. A deep-towed vehicle 18 is attached to the surface vessel 14 by an umbilical cable 16 which provides an electrical, optical and mechanical connection between the deep-towed vehicle 18 and the surface vessel 14. The deep-towed vehicle 18 is towed by the surface vessel 14 such that it remains consistently close to the seafloor 6. This is facilitated by an echo-location package 20 which relays information about the height of the deep-towed vehicle 18 above the seafloor 6 to the surface vessel 14. The deep-towed vehicle 18 receives electrical power from the ship's on-board power supply via the umbilical cable 16.
A cycloconverter unit 30 generates the chosen waveform to be supplied to an EM source in the form of a transmitting antenna 22 which is towed by the deep-towed vehicle 18. The transmitting antenna 22 broadcasts the EM signal into the sea 4, and this results in a component passing through the rock strata 8. One or more remote instrument packages 26 record the signal received by receiving antennae 24 in response to the transmitted EM signal. If the separation of the transmitting antenna 22 and the receiving antenna 24 is greater than a few hundred meters, the highly conductive seawater strongly attenuates the direct signal between them. The components of the EM signal that have travelled through the rock strata 8 and the reservoir 12 dominate the received signal and provide information about the electrical properties of these regions. At the end of the sounding experiment, a remotely operable release system allows the instrument package 26 to be detached from a ballast weight (not shown) so that an in-built flotation device 28 can carry the instrument package 26 to the surface 2 for recovery and retrieval of data for inversion analysis.
The transmitting antenna 22 emits signals that propagate outwards both into the overlying water column 4 and downwards into the seafloor 6 and the underlying strata 8, 12. In both cases, at practical frequencies for this method and given the typical resistivity of the media 4, 8, 12, propagation occurs by diffusion of electromagnetic fields. The rate of decay in amplitude and the phase shift of the signal are controlled both by geometric spreading and by skin depth effects. Because in general the underlying strata 8, 12 are more resistive than seawater 4, skin depths in the underlying strata 8, 12 are longer. As a result, electric fields measured at the seafloor 6 by a receiving antenna 24 at a suitable horizontal range axe dominated by the components of the source fields which have propagated downwards through the seafloor 6, along within the underlying strata 8, 12, and back up to the receiving antenna 24. Both the amplitude and the phase of the received signal depend on the resistivity structure of the underlying strata 8, 12 - and so, in principal, a survey consisting of many transmitter (source) and receiver. (detector) locations can provide a mufti-dimensional image, by geophysical inversion, of sub-seafloor resistivity.
The technique described here exploits the large resistivity contrast that exists between a hydrocarbon reservoir (typically tens of Sam or higher) and the over-and under-lying sediments (typically ~ 2 Sam or less). Such a contrast has a detectable influence on controlled source electromagnetic data collected at the seafloor 6 above the reservoir 12. The effect of the reservoir is most detectable in controlled source electromagnetic data at an appropriate frequency, and if the horizontal range from source 22 to receiver 24 is of the order of 2 to 5 times the depth of burial of the reservoir 12 in typical situations.
This following text describes specific geometric and data reduction approaches that allow the effect of a hydrocarbon reservoir on the outcome of a controlled source electromagnetic survey to be detected and analysed in practice.
Use of a mobile source and multiple fixed receivers In order to achieve a satisfactory survey outcome, it is essential to make controlled source electromagnetic measurements over a broad range of survey geometries - in other words, many transmitting locations and many receiving locations. The transmitter, i.e. the transmitting antenna, requires significant power to drive it, of the order tens of kilowatts or greater for signals detectable at ranges of several kilometres. It must therefore be connected by an umbilical cable 16 to the survey vessel 14 and this makes it relatively straightforward to make the transmitter mobile. It can then be towed in an appropriate survey pattern from the surface survey vessel 14. Since in many situations surveys of this kind are liable to take place over areas of the seafloor 6 where sensitive engineering installations exist or are planned, there are significant advantages to using a transmitter which does not physically come into contact with the seafloor 6. Provided that the separation between the transmitting antenna 22 and the seafloor 6 is small compared to a skin depth of the investigating field in seawater, the survey can still be completed satisfactorily. As a result, the most appropriate form of transmitter to use for this type of survey is one in which a neutrally buoyant horizontal electric dipole antenna 22 is towed behind the deep-towed vehicle 18 at a height of a few metres to a few tens of metres above the seafloor 6.
In the case of the receivers, there is an important advantage in using a static recording device. It is necessary to measure the alternating electric or magnetic field at the seafloor 6, resulting from the signal emitted by the transmitter. In typical applications, the signal-to-noise ratio of the received signal will be critical to the success and resolution of the survey, and so must be maximised. Moving the receiver inevitably generates noise, whether the signal measured is magnetic or electric field.
In the case of electric fields, any motion of the receiver through the conducting seawater medium 4 in the presence of earth's geomagnetic field will generate an electromotive force across the receiving antenna 24. Receiver movements will therefore map into spurious electric field signals in the recorded data. In the case of magnetic field recordings, there are also significant disadvantages to moving the receiver. Most importantly, if vector rather than scalar magnetometers are used (i.e.
measuring individual directional components of the magnetic field), any variation in the orientation of the receiving antenna 24 will again lead to significant spurious signals, since the magnetic detecting element will detect changes in the component of the geomagnetic field aligned with it. As a consequence of these two effects, any translational movement of an electric field sensor or rotational movement of a magnetic field sensor will result in contamination of the received signal by motionally induced noise.
For these reasons, it is desirable to carry out a controlled source electromagnetic survey to investigate or detect the presence of buried hydrocarbons using a combination of a mobile horizontal electric dipole source, equipped with a neutrally buoyant streamed antenna 22 and operated just above the seafloor 6;
and an array of electric and/or magnetic field sensing receivers 24 placed statically on the seafloor 6 as shown in Figure 1. The receiving instruments 26 can be recovered using acoustically actuated release mechanisms to separate them from their ballast weights, allowing them to return to the sea surface for recovery, following standard oceanographic and marine geophysical practice.
The effects of propagation through the atmosphere There is an important additional factor that is crucial to the success or otherwise of the technique described here. In shallow water depths, it is possible for signals from the transmitter to follow a propagation path upwards through the water column to the surface; horizontally through the air; and back down through the water column to the seafloor receiver. This 'air wave' component contains no information about the sub-seafloor resistivity. It tends to dominate the received signal in shallow water and at long source-to-receiver offsets. The effect of the air wave can be minimised by choosing appropriate transmission frequencies, and by targeting surveys on prospects in deep water and in which the target is at a relatively shallow depth below the seafloor.
Figures 2A, 2B and 2C show three cases of one-dimensional (1D) modelling showing signal amplitude as a function of source to receiver range and frequency.
Each of the figures corresponds to a different depth of seawater overlying the seafloor.
Figure 2A shows the modelled decay of electric field strength E as the distance R between the source (transmitter) and detector (receiver) increases. The detector is separated from the receiver along a line which runs parallel to, and passes through the transmitting dipole source. The model survey is undertaken above a semi-infinite subterranean strata of uniform resistivity 1 SZm, and in seawater of depth 750 m. The decay is shown for five different electromagnetic frequency components ranging from 0.5 Hz to 8 Hz.
Figure 2B represents the results of a similar modelling to that shown in Figure 2A, but for a seawater depth of 900 m.
Figure 2C represents the corresponding results for a seawater depth of 1300 m.
On all three sets of curves shown in Figures 2A, 2B and 2C, amplitude decreases rapidly with range, as expected. Additionally, amplitudes at a given range decrease with increasing frequency. This is because higher frequencies have shorter skin depths, and so experience increased attenuation. On each curve, the break in slope to a distinctly shallower gradient on the graph indicates the emergence of the air wave as the dominant signal at the receiver. For example, the curve corresponding to a signal component at 4 Hz shown in Figure 2B shows that the air wave becomes important in seawater of depth 900 m and using a 4 Hz source transmission frequency when the detector is more than about 3.5 km away from the source. It can be seen that this becomes an increasing problem in shallow water and at higher frequencies.
For the technique described here to work most efficiently, the signal at the receiver must be primarily due to propagation through the seafloor, and not due to the air wave. For instance, for this condition to be met at 5 to 6 km source-detector offset for the model subterranean strata described above, the frequency used for the survey should be no higher than 1 Hz in 1300 m water depth, or 0.5 Hz in 900 m water depth.
Modelling can thus be used to plan a survey in order to determine the maximum distance allowable between source and detectors for various frequencies.
Suitable detector deployment positions and source frequency ranges can thus be determined in advance of performing a survey.
The use of survey result parameters derived from polarisation ellipse analysis Use of a horizontal electric dipole antenna produces electromagnetic fields at the seafloor that can be measured using electrometers or scalar or vector magnetometers. Scalar magnetometers are not used widely in practice. Current practice is to measure two or more orthogonal components of either electric or magnetic field. Up to six channels of data can usefully be recorded by the receiving instruments, corresponding to three orthogonal directional components each of electric and magnetic field. For operational reasons of instrument complexity and data storage, in most cases a sub-set of these is recorded.
The commonest arrangement currently is to use an orthogonal pair of horizontal electric dipole receiving elements in each receiving instrument. In the horizontal plane, the transmitter generates both radial and azimuthal components of electric field at the receiver, and, since these have different amplitudes and phases in general, the electric field at the seafloor from a harmonic transmitter maps out a polarisation ellipse in the horizontal plane. The same applies in the case of horizontal magnetic field; and if vertical components are included, the planar ellipse becomes a polarisation ellipsoid.
A straightforward approach to analysing controlled source electromagnetic data is to resolve the observed fields into radial and azimuthal components.
However this suffers from contamination of field values by errors arising from small inaccuracies in source and receiver orientation and position. As an alternative, we therefore propose the novel approach of analysing the data from the recorded components in terms of polarisation ellipse properties - e.g. the amplitude and/or phase of the component of the signal aligned along the major axis of the ellipse. Since this value is much less dependent on the accuracy of the source and receiver positioning - and in particular removes the significant errors that arise from even small uncertainties in receiver orientation - use of polarisation ellipse values can lead to smaller uncertainties in measured field properties at the receiver. We show below as an example that use of the amplitude along the major axis of the polarisation ellipse for horizontal electric field can be used as a robust measurement parameter for detecting the presence of buried hydrocarbon layers. In principle, other polarisation ellipse parameters can be used in a similar way. Parameters that could readily be measured include amplitude or phase along the major axis; the horizontal polarisation ellipse of either electric or magnetic field; or in either case a polarisation ellipsoid including a vertical component.
Presentation of survey results using normalised values As seen in Figures 2A, 2B and 2C above, a feature of a controlled source electromagnetic survey over a sedimented seafloor is that the received field amplitude drops very rapidly with increasing range. The received signal properties also depend on the orientations of the source and receiving dipoles. In a typical survey, the signal amplitude is likely to vary by several orders of magnitude over the useful set of source-detector (transmitter-receiver) offsets; while the phase lag of the received signal also increases steadily with increasing offset. The effects of buried structure beneath the survey show up as deviations in the signal from that which would be expected if the sub-surface consisted only of a homogeneous, isotropic half space.
In order to visualise more clearly the effects of buried structure on controlled source electromagnetic signal properties, it is convenient to normalise the observed signals with respect to some reference model. A simple model to use would consist of a water layer of finite thickness above the survey of the true depth; and a homogeneous isotropic half space representing the seafloor below the survey.
If appropriate, and if the relevant a priori information is available, a more complex reference model may be used, although it is desirable to use always the simplest reference model that can broadly represent the large scale background properties of the seafloor.
The normalisation is carried out for amplitude by dividing the observed field by that calculated for the reference model, using the same source-detector geometry as for the observed data. In the case of phase, the normalisation is carried out by subtracting the phase lag calculated for the reference model from that of the observed data. We show below that use of normalised values calculated in this way can dramatically enhance the visual presentation of controlled source electromagnetic data from a survey over a hydrocarbon reservoir. Our illustrations use normalised amplitudes, but normalised phases could equally be used.
As an alternative, normalisation could be based on the survey data itself, for example using data collected adj acent to the target.
The varying~physical response of a survey de~endin~ on source orientation ~eometry Overall, the most crucial factor for the success or otherwise of the controlled source electromagnetic technique in practical applications related to hydrocarbon reservoirs is related to survey geometry.
Figure 3 is a schematic plan view from above illustrating a suitable co-ordinate system for describing the relative placement of a controlled source electromagnetic source 30 and an electromagnetic detector 32. The position of the detector 32 with respect to the source 30 is most suitably described in polar co-ordinates, with the centre of the source 30 providing the origin of the polar co-ordinate system.
The position of the detector 32 is defined by an azimuthal angle 8 and a distance R. The detector 32 and source 30 are considered to be co-planar. The angle 0 is measured clockwise from a line 33 passing through, and running parallel to, the source axis, as indicated in Figure 3 by the line marked 8=0°. A detector placed along the line 33, i.e.
such that is has an azimuthal angle 8 of 0°, is referred to as being in an end-on position. A detector with an azimuthal angle ~ of 90°, such that it lies on the line 34 marked A=90° in Figure 3, is referred to as being in a broadside position. The electric field strength at the detector is considered in terms of a radial component EP
and an orthogonal azimuthal component Ee as indicated in Figure 3.
Depending on the relative placements and orientations of sources and receivers, the physics of propagation of the signal through the seafloor - and so the net effect on the properties of the received signal - can be more or less sensitive to different classes of sub seafloor structure. As a simple illustration of this in the context of hydrocarbon surveys, we consider a case in which the sub-seafloor structure can be represented by a simple stack of horizontal layers (a useful first approximation for many situations in sedimentary basins).
Figure 4A shows in schematic vertical section an exemplary subterranean strata configuration. A section of seafloor 42 lies beneath 800 m of seawater 40. The strata beneath the seafloor 42 comprise a 1 km thick overburden layer 44, representing sediments above a hydrocarbon reservoir layer 46. This overburden layer 44 has low resistivity, primarily due to aqueous saturation of pore spaces. The middle layer 46 corresponds to a 100 m thick hydrocarbon reservoir and has a resistivity perhaps 100 times greater than the overburden 44. This is due to the presence of the non-conducting hydrocarbon occupying much of the pore spaces. Below the hydrocarbon reservoir layer 46 is a sub-reservoir region 48 which has low resistivity due to its similarity to the overburden layer 44 and extends downwards for an effectively infinite distance.
In the case of an electromagnetic survey, we can consider the differing behaviours of electric currents generated by the source resolved along horizontal and vertical directions. The relationships between the electric currents flowing in two adjacent regions of space are determined by both galvanic (i.e. direct transfer of charge) and inductive effects. Since charge is conserved, current leaving one volume of the sub-surface strata 44, 46, 48 and arriving in an adjacent volume along the direction of current flow are related to each other by a galvanic mechanism.
On the other hand, if two volumes are close together but separated from each other along a direction orthogonal to current flow, then the coupling between the currents flowing in the two volumes will be primarily inductive.
Applying this to our model of a sub-seafloor structure shown in Figure 4A, we can infer that the effect of the thin but resistive hydrocarbon reservoir layer 46 on the survey results will depend strongly on the direction of flow of the currents generated by the source. If the current at the base of the overburden layer 44 is dominantly horizontal, then coupling between this layer and the deeper layers 46, 48 will be dominated by inductive effects. Although little current will flow in the hydrocarbon reservoir layer 46, the inductive coupling between the overburden 44 and the underlying layer 48 will be only mildly affected. Currents in the two conducting layers 44, 48 (overburden and underlying layers) will therefore be similar to the case without a hydrocarbon reservoir layer 46. The results of such a survey will therefore be only weakly affected by the presence of the hydrocarbon reservoir layer 46.
In contrast, if our survey generates significant vertical components of electric current flow in the overburden layer 44, then galvanic effects along the current path will be strongly affected by the thin resistive hydrocarbon reservoir layer 46. The resistive hydrocarbon reservoir layer 46 will tend to block the current flow.
As a result, currents flowing in the deeper layer 48 will be much reduced; and the overall pattern of current flow in the overburden layer 44 will therefore be significantly altered. In such a survey situation, the presence of the hydrocarbon will significantly influence the outcome.
Hence it becomes of paramount importance in designing a practical survey approach for detecting buried hydrocarbon reservoirs to distinguish between source (transmitter) and detector (receiver) geometries in which coupling between layers is largely inductive between current sheets in a horizontal plane (in which case the survey has little sensitivity to the presence of the reservoir); and those in which a significant component of vertical current flow occurs, in which case blocking of the passage of this current flow by the reservoir leads to a survey which is strongly sensitive to the presence of the reservoir.
A specific approach to combining two contrasting source ~eometries when pros~ectin~ f, or hydrocarbon filled reservoirs By selecting a suitable survey geometry, it is possible to exploit the diverse properties of electromagnetic induction outlined above, by designing a survey in such a way that it provides data that are sensitive to the presence of thin resistive layers (exploiting largely galvanic effects and vertical components of induced current flow);
while simultaneously obtaining data that are more sensitive to the much larger scale structure, but relatively insensitive to a thin hydrocarbon layer (exploiting the dominantly inductive effects). The reason why the latter component is essential is that many features other than hydrocarbon reservoirs can affect the resistivity beneath the seafloor and the results of a survey. For example, resistivity often increases steadily with depth in submarine sedimentary basins, due to the progressive expulsion of pore fluids by rising overburden pressure.
Figure 4B shows in vertical section a highly schematic exemplary subterranean strata configuration which exhibits increasing resistivity with depth. A
section of seafloor 52 lies beneath 800 m of seawater 50. The strata beneath the seafloor comprise a series of layers of sediment of increasing resistivity. A first layer 54 has a uniform resistivity of 1 S2m and a thickness of 1 km. A second layer 56 has a uniform resistivity of 5 S2m and a thickness of 1 km. A third layer 58 has a uniform resistivity of 10 S2m and a thickness of 1 km. A fourth layer 60 has a uniform resistivity of 50 S2m and a thickness of 1 lmn. Beneath the fourth layer 60 is a fifth layer 62 which has a resistivity of 100 SZm and extends downwards for an infinite extent.
The increasing resistivity indicated in Figure 4B leads to longer skin depths, and hence to larger observed amplitudes and smaller observed phase lags for all detector placements. Although this does not exactly reproduce the same galvanic effects of a thin resistive layer, the overall effect on observed fields in the end-on geometry is likely to be very similar to (and difficult or impossible to distinguish from) the effect of a hydrocarbon reservoir. Designing a survey in such a way that different parts of the resulting data are more or less sensitive to the presence of certain key features is essential for removing ambiguity in the interpretation.
Recent studies of volcanic systems at mid-ocean ridges and of sediments beneath resistive basalt layers (MacGregor et al., 1998, 2001; MacGregor &
Sinha, 2000) have demonstrated that improved model resolution and reduced interpretational ambiguity can be achieved by using received electric field components from two distinct source-detector geometries. This approach is now applied to the case of surveying for thin hydrocarbon reservoirs.
Figures 4C, 4D, 4E and 4F show survey simulation results for the two subterranean strata configurations shown in Figures 4A and 4B. The method relies on collecting survey data using two distinct geometric modes. The first mode corresponds to an end-on geometry, in which, as described above, the receiver lies along the axis of the transmitting dipole. The source-detector azimuth B as defined in Figure 3 is 0°, and the field observed at the receiver is dominated by the radial component EP. The second mode corresponds to a broadside geometry, in which, as described above, the source-detector azimuth 8 is 90°, and the field observed at the receiver is dominated by the azimuthal component Ee.
For the case of hydrocarbon exploration, where the target is a thin resistive layer embedded between a more conductive underburden and overburden, such as that shown schematically in Figure 4A, the received signal depends on the azimuthal angle 0 of the detector. Because of the transmitter field pattern, the end-on, or EP, geometry results in a significant contribution to the observed field at the seafloor by the vertical component of current flow. The broadside, or Ee, geometry results in fields at the seafloor that are more dependent on the contribution of inductively coupled currents flowing in horizontal planes. As a result, the presence of a resistive hydrocarbon reservoir layer produces a significant increase at certain ranges in the end-on (EP) fields, while having virtually no effect on the broadside (Eo) fields. The result is 'splitting' of amplitudes between the two modes, and this splitting is characteristic of a structure in which resistivity first increases, in this case due to a hydrocarbon layer, and then decreases with depth.
Figure 4C is a graph showing the modelled amplitudes of electric field E along the major axes of the horizontal polarisation ellipses for two source-detector geometries and as a function of distance R between source and detector. These are determined for the model of sub-seafloor structure indicated in Figure 4A.
Four curves are shown. The curve labelled Eo shows the field strength as a function of distance for end-on geometries, the curve labelled Ego shows the same for broadside geometries.
The curves labelled No and N9o are the corresponding normalisation curves.
These are determined for the same geometries as for Eo and E9o, but represent the modelled amplitudes of electric field along the major axes of the horizontal polarisation ellipses where the subterranean strata comprise a homogeneous isotropic half space of resistivity 1 S2m, such as previously described.
Figure 4D is a graph showing the set of curves which correspond to those curves shown in Figure 4C, but for the subterranean strata configuration shown in Figure 4B.
Figure 4E is a graph showing the normalised values of the fields Eo and E9o shown in Figure 4C. These fields have been normalised by the fields No and N9o, also shown in Figure 4C, as discussed above.
Figure 4F is a graph showing the normalised values of the fields Eo and E9o shown in Figure 4D. These fields have been normalised by the fields No and N9o, also shown in Figure 4D, as discussed above.
It can be seen in Figure 4E that the presence of the~reservoir has relatively little effect on the broadside amplitudes (i.e. the curve labelled E9o/N9o is close to unity), but that between 4 and 6 km offset, it produces a substantial increase in the end-on amplitudes (i.e. the curve labelled Eo/No is strongly enhanced over this range). The hydrocarbon layer in this particular model leads to amplitude splitting between the two geometrical modes by about a factor of 10.
It can be seen in Figure 4F that the steadily increasing resistivity with depth model shown in Figure 4B strongly affects both geometries (i.e. both curves are strongly enhanced between 4 and 6 km offset) and no significant splitting between the two geometric modes occurs.
The end-on geometry data alone cannot distinguish between the two subterranean strata configurations shown in Figures 4A and 4B. These comparative examples conclusively demonstrate the fact that collection of survey data from both end-on and broadside geometries is needed, in order to reliably detect the hydrocarbon layer and eliminate "false-positive" detection of a rock formation containing no hydrocarbon reservoir, but merely having an increasing resistivity with depth, which is not unusual.
Lastly, the normalised curves in Figures 4E and 4F start off close to unity, at short ranges where the data are insensitive to the buried structure. The effects of the buried structure become greatest at between 5 and 6 km range. Beyond this, the normalised curves rapidly return towards a value of l, because at long ranges, air wave propagation (insensitive to sub-seafloor structure) begins to dominate, thus masking the effects of the sub surface. As a result, the presence of the target structure can best be observed in controlled source electromagnetic data of this type over a relatively narrow window of source-detector separation ranges.
An efficient electromagnetic survey may use a large number of detectors distributed over a target area and a signal transmitter which operates continuously as it is towed along an extended tow path. This follows from the relatively low cost of deploying detectors and the relatively high cost of deploying a transmitter.
Figure 5 is a schematic plan view showing an example layout of sixteen detectors 64 distributed across a section of seafloor 65 above subterranean reservoir 66. The reservoir 66 has a linear extent of several km and its boundary is indicated by a heavy line 67. In this example, the detectors 64 are uniformly distributed in a square-grid pattern so as to approximately cover the substantially square reservoir 66.
In performing a survey, a source (not shown) starts from the location marked 'A' in Figure 5 and is towed, whilst broadcasting continuously, along a path indicated by the broken line 68, the survey is completed when the source reaches the location marked 'B'. Data are continuously recorded by the detectors 64 throughout the towing process and the position of the source transmitter relative to the detector network is also logged.
During the towing process, each of the detectors 64 presents several different orientations to the source. For example, when the source is directly above the detector position D 1 and on the vertical section of the tow path, the detectors at positions D2 and D3 are at different ranges in an end-on position, the detectors at positions D4 and DS are at different ranges in a broadside position and the detector at position D6 is midway between. However, when the source later passes over the detector position D1 when on the horizontal section of the tow path, the detectors at positions D2 and D3 are now in a broadside position, and the detectors at position D4 and DS are in an end-on position. Thus, in the course of a survey, and in conjunction with the positional information of the source, data from the detectors 64 can be used to provide details of the signal transmission through the subterranean strata for a comprehensive range of distances and orientations between source and detector, each with varying galvanic and inductive contributions to the signal propagation. In this way a simple continuous towing of the source transmitter can provide a detailed survey which covers the extent of the subterranean reservoir 66.
Although the above example is based on a square detector grid, it will be understood that a wide variety of detector placements may be used, for example other high symmetry regular grids, such as triangular or rectangular, may be used.
In addition irregular grids may be used that have no high level of symmetry.
On the detection of the edges of a reservoir by this techniaue We showed above that our proposed arrangement of survey geometries can both detect a hydrocarbon reservoir, and distinguish it from a general increase in resistivity with depth. The modelling though was based only on 1-D
assumptions, i.e.
horizontal layers of infinite horizontal extent. By using higher dimensional modelling, we show below that the survey technique also works well in the case of reservoirs of finite extent; and indeed that it can be used as a reliable means of locating the edges of a reservoir structure.
Figure 6A shows a schematic section view of a model subterranean strata configuration containing a hydrocarbon reservoir layer of finite horizontal extent. A
section of seafloor 70 lies beneath a body of sea 72 which is 1.15 km deep and has a resistivity of 0.31 S2m. A planar hydrocarbon reservoir 74 is 0.15 km thick and its base is 1 km below the seafloor .70. The reservoir 74 has a resistivity of 100 Sam and extends infinitely out of the plane of Figure 6A, and semi-infinitely within the plane of Figure 6A. A uniform sedimentary structure 76 of resistivity 1 S2m surrounds the reservoir 74 and fills the remaining space beneath the seafloor 70. For the purposes of the model, two transmitter positions are considered. A first dipole transmitter 75 is aligned parallel to and directly above the edge 77 of the reservoir 74 on the seafloor 70. The dipole axis of the transmitter 75 is perpendicular to the plane of Figure 6A
and a cross is used to simply indicate its position rather than represent its structure. A
second transmitter 78 is also on the seafloor 70 and aligned parallel to the first transmitter 75. The second transmitter 78 is horizontally displaced from the first transmitter 75 by 1 km in a direction which positions it directly above the reservoir 74 as indicated in Figure 6A.
Figures 6B and 6C show the results of 2.5 D modelling (3-D source, 2-D
resistivity structure) of the exemplary subterranean strata configuration shown in Figure 6A for the first and second transmitter positions 75 and 78 respectively in order to investigate the effects of the finite areal extent of the resistive hydrocarbon reservoir 74. Figures 6B and 6C are plan views representing the normalised fields at the seafloor 70 as a function of position. The normalised field component shown in each case is the semi-major axis of the horizontal polarisation ellipse of electric field.
In Figure 6B the first transmitter 75 is at the origin of the local co-ordinate system and its dipole axis is parallel to the vertical axis in this system. In Figure 6C
the second transmitter 78 is at the origin of the local co-ordinate system and its dipole axis is parallel to the vertical axis in this system.
The normalised field amplitude is shown in a grey scale representation, while contour lines indicate the absolute value of amplitude. Figures 6B and 6C show that the reservoir structure causes substantial increases in the fields at positions close to end-on geometries, this is especially so when the source and detector are both over the reservoir, as in Figure 6C. If the source or detector are at or outside the edge of the reservoir, the amplitude increase is not seen. At detector locations close to broadside geometry, the increase in amplitude is again not seen, even if both source and detector are over the reservoir. Figures 6B and 6C illustrate that the amplitude splitting effect is observable over reservoirs of finite size, as well as over the 1-D
(layered) structures discussed previously.
The marked difference between Figures 6B and 6C shows that the observed fields are extremely sensitive to the location of the source (transmitter) with respect to the edge of the reservoir. Tn Figure 6B, the marked asymmetry of the fields about the y axis shows that - provided the source is over the reservoir - the observed fields are also highly sensitive to the location of the detector (receiver) with respect to the edge of the reservoir. Thus by careful siting of both source tow tracks and detectors relative to a suspected hydrocarbon bearing structure, the survey method described here can provide detailed information on the areal extent of the hydrocarbon, as well as on its existence or otherwise.
In addition to initially surveying for subterranean hydrocarbon reservoirs, the utility of the invention in detecting the edges of reservoirs makes it applicable to assessing changes in hydrocarbon reservoir content over a period of time. This is of _~g_ particular value, for example, where a hydrocarbon reservoir is being exploited. As water replaces the hydrocarbon which is drawn from the reservoir, regular follow up surveys can be used to determine the change in the distribution of the remaining hydrocarbon. Such surveys could therefore be used to characterise the evolution of the reservoir through time, and hence enable more efficient extraction and reservoir management.
Source and receiver locations - design considerations for real surveys The methodology described in this document can be readily applied to real surveys for hydrocarbon resources on continental margins. The required transmitter and receiver characteristics in terms of power, signal to noise ratios and operating parameters can be met by existing instrumental technology.
In order to plan a successful survey, the following factors are important:
An understanding of the need to collect observations over an appropriate set of source-detector separation ranges, and using both broadside and end-on surveying geometries;
A designed layout of the detector array and of source tow line locations and orientations optimised for the above;
A designed layout for the survey that also takes account of the need to locate the limits of the areal extent of any reservoir, based on the approach presented in Figures 6A, 6B and 6C;
Prior modelling to establish the range of transmission frequencies, and the set of source-detector separation ranges, that should be used, taking into consideration the expected resistivity structure of the seafloor, the depth of the taxget, the depth of the seafloor, and the influence of the air wave on the data collected.
Conclusions We have shown that the application of a specific set of methods to a controlled source electromagnetic survey of the seafloor can allow the technique to be successfully applied to the problem of detecting hydrocarbons beneath the seafloor.
Specifically, we have shown that:
The survey should be carried out using a combination of a mobile horizontal electric dipole source transmitter, equipped with a neutrally buoyant streamed antenna and operated just above the seafloor; and an array of electric and/or magnetic field sensing receivers placed statically on the seafloor.
The effects of propagation of the signal through the atmosphere (the 'air wave') are significant at high frequencies, at long source-detector separation distances, and in shallow water. This effect limits the applicability of the method by favouring detection of structures relatively shallowly situated beneath the seafloor, but in deep water; and limits the choice of frequencies and source-detector offsets for the survey.
Improvements in data analysis can be achieved by making use of field measurement parameters aligned along the major axis of the polarisation ellipse at the seafloor.
For purposes of data presentation and interpretation it is desirable to make use of field parameters that have been normalised with reference to an appropriate simplified model of the sub-seafloor structure.
In order to resolve the presence of any hydrocarbon, and to distinguish its effect on the data from other likely structures and reduce the ambiguity of interpretation, it is essential to collect survey data from both the end-on and broadside geometric configurations; and to analyse the data in terms of splitting between the two geometric modes.
2.5 D modelling shows that provided an appropriately designed array of receivers, and an appropriate set of transmitter tow tracks, are used, then the method can yield valuable information about the limits of the areal extent of any sub-seafloor hydrocarbon reservoir, as well as detecting its presence.
Finally it will be understood that the invention is equally applicable to surveying of freshwater, for example large lakes, so that references to seafloor, seawater etc. should not be regarded as limiting.
References:
1. Sinha, M.C., Patel, P.D., Unsworth, M.J., Owen, T.R.E. & MacCormack, M.R.G.
An active source electromagnetic sounding system for marine use. Mar. Geophys.
Res.,12, 1990, 59-68.
2. MacGregor, L.M., Constable, S.C. & Sinha, M.G. The RAMESSES experiment III: Controlled source electromagnetic sounding of the Reykjanes Ridge at 570 45' N. Geophysical Journal International, 135, 1998, 773-789.
3. MacGregor, L. M. & Sinha, M.G. Use of marine controlled source electromagnetic sounding for sub-basalt exploration. Geophysical Pf°ospecting, 48, 2000, 1091-1106.
4. MacGregor, L., Sinha, M. ~ Constable, S. Electrical resistivity structure of the Valu Fa Ridge, Lau Basin, from marine controlled source electromagnetic sounding. Geophys. J. Ir~t.,146, 217-236, 2001.
5. Evans, R. L., Sinha, M. C., Constable, S. C. & Unsworth, M. J. On the electrical nature of the axial melt zone at 130N on the East Pacific Rise. J. Geophys.
Res., 99, 1994, 577 - 588.
6. WO 00/13046 7. WO 01/57555
In a preferred embodiment, the first and second data sets each comprise radial and azimuthal components of electric field or magnetic field measured at the detector, and the method further comprises: transforming the radial and azimuthal components into at least one polarisation ellipse parameter, prior to the combining. The polarisation ellipse parameters) can be the amplitude and/or phase of the component of the electric field or magnetic field aligned along a major axis of the ellipse.
The method may advantageously further comprise: visually representing the results data set in a plot of at least two dimensions corresponding to the survey area.
The visual representation can be a two-dimensional (2D) plot in plan view, or a three-dimensional (3D) plot including depth, e.g. a perspective view of the survey area.
The plot may include markings of areas of equal or similar electromagnetic field strength. These markings could be contour lines, or colour or grey scale gradations with one colour or grey tone being used for a range of electromagnetic field strength values. Stepwise colour gradation is used in the preferred embodiment, with each colour representing a defined range of data values.
When the plot is of normalised survey data it is helpful if the plot also includes lines of equal absolute electromagnetic field strength, showing how the signal strength has decayed as one moves away from the source. The lines, which may appear as contours, can be labelled with either a relative or absolute decay value.
According to a fourth aspect of the invention, there is provided a method of planning an electromagnetic survey of an area identified as potentially containing a subsea hydrocarbon reservoir, comprising: creating a model of the area to be surveyed, including a rock formation containing a hydrocarbon reservoir and a body of water above the rock formation; setting values for water depth, depth below the seafloor of the hydrocarbon reservoir, and resistivity structure of the rock formation;
performing a simulation of an electromagnetic survey in the model of the survey area by calculating first and second survey data sets for an electromagnetic source aligned end-on and broadside relative to a detector; and combining the first and second survey data sets to obtain a results data set that represents a difference between the end-on and broadside alignments.
In typical use, the simulation will be repeated for a number of source frequencies and/or source-to-detector distances, and/or other parameters which can be varied during a survey. This iterative procedure can be used in order to select optimum surveying conditions in terms of source frequency and source-to-detector distance for probing the hydrocarbon reservoir. By optimum, it is not necessarily meant that the best survey conditions are found, but only that a set of survey conditions is found which will provide strong, unmistakable indications in the case that there is a hydrocarbon reservoir at the survey site. The iterative procedure may be purely under manual control. However, preferably, the simulator can allow the user the option of automatically optimising the survey conditions. The user can then switch between manual and automated iteration as desired.
The model should preferably include a body of air above the body of water, so that the simulation can take account of signal propagation paths including the body of air when calculating the first and second survey data sets. The propagation path through the air (the 'air wave') will in fact dominate for shallower water and longer distances between source and detector(s), as will be apparent from the examples described fiuther below. For deep water and shorter detector-source distances the effect is less important and may be omitted from the model.
The method preferably further comprises: normalising each of the first and second survey data sets relative to respective first and second normalisation data sets or functions specific to the end-on and broadside alignments respectively, prior to the combining. The first and second normalisation data sets or functions can be calculated from a rock formation model, for example.
The method preferably further comprises visually representing the results data set in a plot of at least two dimensions corresponding to the survey area. The visual representation can be a 2D plot in plan view, or a 3D plot including depth, e.g. a _$_ perspective view of the survey area. The other comments made above in relation to the plots of the third aspect of the invention also apply to the fourth aspect of the invention.
Another aspect of the invention relates to a computer program product bearing machine readable instructions for implementing the method of analysing results from an electromagnetic survey according to the third aspect of the invention.
Another aspect of the invention relates to a computer apparatus loaded with machine readable instructions for implementing the method of analysing results from an electromagnetic survey according to the third aspect of the invention.
Another aspect of the invention relates to a computer program product bearing machine readable instructions for implementing the method of planning an electromagnetic survey according to the fourth aspect of the invention.
Another aspect of the invention relates to a computer apparatus loaded with machine readable instructions for implementing the method of planning an electromagnetic survey according to the fourth aspect of the invention.
BRIEF DESCRIPTION OF THE DRAWINGS
For a better understanding of the invention and to show how the same may be carried into effect reference is now made by way of example to the accompanying drawings.
Figure 1 shows a schematic representation of an exploratory EM sounding survey.
Figure 2A is a graph showing the modelled electric field strength as a function of distance from, and along the axis of, the EM source. This is shown for several frequency components. The modelled sea depth is 750 m.
Figure 2B is a graph showing the modelled electric field strength as a function of distance from, and along the axis of, the EM source. This is shown for several frequency components. The modelled sea depth is 900 m.
Figure 2C is a graph showing the modelled electric field strength as a function of distance from, and along the axis of, the EM source. This is shown for several frequency components. The modelled sea depth is 1300 m.
Figure 3 is a schematic plan view defining a survey geometry coordinate system.
Figure 4A shows a schematic section view of a model subterranean strata configuration containing a hydrocarbon reservoir.
Figure 4B shows a schematic section view of a model subterranean strata configuration containing layers of increasing resistivity.
Figure 4C is a graph showing the electric field strength (Eo, E9o) as a function of distance from the source for two different detector positions for the subterranean strata shown in Figure 4A. Also shown are the normalising signals (No, N9o) expected from a uniform subterranean strata configuration.
Figure 4D is a graph showing the electric field strength (Eo, E9o) as a function of distance from the source for two different detector positions for the subterranean strata shown in Figure 4B. Also shown are the normalising signals (No, N9o) expected from a uniform subterranean strata configuration.
Figure 4E is a graph showing the electric field strengths (Eo, E9o) shown in Figure 4C normalised by the normalising field strengths (No, N9o) shown in the same.
Figure 4F is a graph showing the electric field strengths (Eo, E9o) shown in Figure 4D normalised by the normalising field strengths (No, N9o) shown in the same.
Figure 5 is a schematic plan view showing an arrangement of sixteen detectors on a section of sea floor above a subterranean reservoir, a particular source transmitter tow path is also shown.
Figure 6A shows a schematic section view of a model subterranean strata configuration containing a hydrocarbon reservoir.
Figure 6B is a plan view representing the normalised electric field strength as a function of location for the model subterranean strata configuration shown in Figure 6A for a first transmitter position.
Figure 6C is a plan view representing the normalised electric field strength as a function of location for the model subterranean strata configuration shown in Figure 6A for a second transmitter position.
DETAILED DESCRIPTION
Figure 1 of the accompanying drawing shows schematically a surface vessel 14 undertaking EM sounding of subterranean rock strata 8 within which a hydrocarbon reservoir 12 is located. The surface vessel 14 floats on the surface 2 of the sea 4. A deep-towed vehicle 18 is attached to the surface vessel 14 by an umbilical cable 16 which provides an electrical, optical and mechanical connection between the deep-towed vehicle 18 and the surface vessel 14. The deep-towed vehicle 18 is towed by the surface vessel 14 such that it remains consistently close to the seafloor 6. This is facilitated by an echo-location package 20 which relays information about the height of the deep-towed vehicle 18 above the seafloor 6 to the surface vessel 14. The deep-towed vehicle 18 receives electrical power from the ship's on-board power supply via the umbilical cable 16.
A cycloconverter unit 30 generates the chosen waveform to be supplied to an EM source in the form of a transmitting antenna 22 which is towed by the deep-towed vehicle 18. The transmitting antenna 22 broadcasts the EM signal into the sea 4, and this results in a component passing through the rock strata 8. One or more remote instrument packages 26 record the signal received by receiving antennae 24 in response to the transmitted EM signal. If the separation of the transmitting antenna 22 and the receiving antenna 24 is greater than a few hundred meters, the highly conductive seawater strongly attenuates the direct signal between them. The components of the EM signal that have travelled through the rock strata 8 and the reservoir 12 dominate the received signal and provide information about the electrical properties of these regions. At the end of the sounding experiment, a remotely operable release system allows the instrument package 26 to be detached from a ballast weight (not shown) so that an in-built flotation device 28 can carry the instrument package 26 to the surface 2 for recovery and retrieval of data for inversion analysis.
The transmitting antenna 22 emits signals that propagate outwards both into the overlying water column 4 and downwards into the seafloor 6 and the underlying strata 8, 12. In both cases, at practical frequencies for this method and given the typical resistivity of the media 4, 8, 12, propagation occurs by diffusion of electromagnetic fields. The rate of decay in amplitude and the phase shift of the signal are controlled both by geometric spreading and by skin depth effects. Because in general the underlying strata 8, 12 are more resistive than seawater 4, skin depths in the underlying strata 8, 12 are longer. As a result, electric fields measured at the seafloor 6 by a receiving antenna 24 at a suitable horizontal range axe dominated by the components of the source fields which have propagated downwards through the seafloor 6, along within the underlying strata 8, 12, and back up to the receiving antenna 24. Both the amplitude and the phase of the received signal depend on the resistivity structure of the underlying strata 8, 12 - and so, in principal, a survey consisting of many transmitter (source) and receiver. (detector) locations can provide a mufti-dimensional image, by geophysical inversion, of sub-seafloor resistivity.
The technique described here exploits the large resistivity contrast that exists between a hydrocarbon reservoir (typically tens of Sam or higher) and the over-and under-lying sediments (typically ~ 2 Sam or less). Such a contrast has a detectable influence on controlled source electromagnetic data collected at the seafloor 6 above the reservoir 12. The effect of the reservoir is most detectable in controlled source electromagnetic data at an appropriate frequency, and if the horizontal range from source 22 to receiver 24 is of the order of 2 to 5 times the depth of burial of the reservoir 12 in typical situations.
This following text describes specific geometric and data reduction approaches that allow the effect of a hydrocarbon reservoir on the outcome of a controlled source electromagnetic survey to be detected and analysed in practice.
Use of a mobile source and multiple fixed receivers In order to achieve a satisfactory survey outcome, it is essential to make controlled source electromagnetic measurements over a broad range of survey geometries - in other words, many transmitting locations and many receiving locations. The transmitter, i.e. the transmitting antenna, requires significant power to drive it, of the order tens of kilowatts or greater for signals detectable at ranges of several kilometres. It must therefore be connected by an umbilical cable 16 to the survey vessel 14 and this makes it relatively straightforward to make the transmitter mobile. It can then be towed in an appropriate survey pattern from the surface survey vessel 14. Since in many situations surveys of this kind are liable to take place over areas of the seafloor 6 where sensitive engineering installations exist or are planned, there are significant advantages to using a transmitter which does not physically come into contact with the seafloor 6. Provided that the separation between the transmitting antenna 22 and the seafloor 6 is small compared to a skin depth of the investigating field in seawater, the survey can still be completed satisfactorily. As a result, the most appropriate form of transmitter to use for this type of survey is one in which a neutrally buoyant horizontal electric dipole antenna 22 is towed behind the deep-towed vehicle 18 at a height of a few metres to a few tens of metres above the seafloor 6.
In the case of the receivers, there is an important advantage in using a static recording device. It is necessary to measure the alternating electric or magnetic field at the seafloor 6, resulting from the signal emitted by the transmitter. In typical applications, the signal-to-noise ratio of the received signal will be critical to the success and resolution of the survey, and so must be maximised. Moving the receiver inevitably generates noise, whether the signal measured is magnetic or electric field.
In the case of electric fields, any motion of the receiver through the conducting seawater medium 4 in the presence of earth's geomagnetic field will generate an electromotive force across the receiving antenna 24. Receiver movements will therefore map into spurious electric field signals in the recorded data. In the case of magnetic field recordings, there are also significant disadvantages to moving the receiver. Most importantly, if vector rather than scalar magnetometers are used (i.e.
measuring individual directional components of the magnetic field), any variation in the orientation of the receiving antenna 24 will again lead to significant spurious signals, since the magnetic detecting element will detect changes in the component of the geomagnetic field aligned with it. As a consequence of these two effects, any translational movement of an electric field sensor or rotational movement of a magnetic field sensor will result in contamination of the received signal by motionally induced noise.
For these reasons, it is desirable to carry out a controlled source electromagnetic survey to investigate or detect the presence of buried hydrocarbons using a combination of a mobile horizontal electric dipole source, equipped with a neutrally buoyant streamed antenna 22 and operated just above the seafloor 6;
and an array of electric and/or magnetic field sensing receivers 24 placed statically on the seafloor 6 as shown in Figure 1. The receiving instruments 26 can be recovered using acoustically actuated release mechanisms to separate them from their ballast weights, allowing them to return to the sea surface for recovery, following standard oceanographic and marine geophysical practice.
The effects of propagation through the atmosphere There is an important additional factor that is crucial to the success or otherwise of the technique described here. In shallow water depths, it is possible for signals from the transmitter to follow a propagation path upwards through the water column to the surface; horizontally through the air; and back down through the water column to the seafloor receiver. This 'air wave' component contains no information about the sub-seafloor resistivity. It tends to dominate the received signal in shallow water and at long source-to-receiver offsets. The effect of the air wave can be minimised by choosing appropriate transmission frequencies, and by targeting surveys on prospects in deep water and in which the target is at a relatively shallow depth below the seafloor.
Figures 2A, 2B and 2C show three cases of one-dimensional (1D) modelling showing signal amplitude as a function of source to receiver range and frequency.
Each of the figures corresponds to a different depth of seawater overlying the seafloor.
Figure 2A shows the modelled decay of electric field strength E as the distance R between the source (transmitter) and detector (receiver) increases. The detector is separated from the receiver along a line which runs parallel to, and passes through the transmitting dipole source. The model survey is undertaken above a semi-infinite subterranean strata of uniform resistivity 1 SZm, and in seawater of depth 750 m. The decay is shown for five different electromagnetic frequency components ranging from 0.5 Hz to 8 Hz.
Figure 2B represents the results of a similar modelling to that shown in Figure 2A, but for a seawater depth of 900 m.
Figure 2C represents the corresponding results for a seawater depth of 1300 m.
On all three sets of curves shown in Figures 2A, 2B and 2C, amplitude decreases rapidly with range, as expected. Additionally, amplitudes at a given range decrease with increasing frequency. This is because higher frequencies have shorter skin depths, and so experience increased attenuation. On each curve, the break in slope to a distinctly shallower gradient on the graph indicates the emergence of the air wave as the dominant signal at the receiver. For example, the curve corresponding to a signal component at 4 Hz shown in Figure 2B shows that the air wave becomes important in seawater of depth 900 m and using a 4 Hz source transmission frequency when the detector is more than about 3.5 km away from the source. It can be seen that this becomes an increasing problem in shallow water and at higher frequencies.
For the technique described here to work most efficiently, the signal at the receiver must be primarily due to propagation through the seafloor, and not due to the air wave. For instance, for this condition to be met at 5 to 6 km source-detector offset for the model subterranean strata described above, the frequency used for the survey should be no higher than 1 Hz in 1300 m water depth, or 0.5 Hz in 900 m water depth.
Modelling can thus be used to plan a survey in order to determine the maximum distance allowable between source and detectors for various frequencies.
Suitable detector deployment positions and source frequency ranges can thus be determined in advance of performing a survey.
The use of survey result parameters derived from polarisation ellipse analysis Use of a horizontal electric dipole antenna produces electromagnetic fields at the seafloor that can be measured using electrometers or scalar or vector magnetometers. Scalar magnetometers are not used widely in practice. Current practice is to measure two or more orthogonal components of either electric or magnetic field. Up to six channels of data can usefully be recorded by the receiving instruments, corresponding to three orthogonal directional components each of electric and magnetic field. For operational reasons of instrument complexity and data storage, in most cases a sub-set of these is recorded.
The commonest arrangement currently is to use an orthogonal pair of horizontal electric dipole receiving elements in each receiving instrument. In the horizontal plane, the transmitter generates both radial and azimuthal components of electric field at the receiver, and, since these have different amplitudes and phases in general, the electric field at the seafloor from a harmonic transmitter maps out a polarisation ellipse in the horizontal plane. The same applies in the case of horizontal magnetic field; and if vertical components are included, the planar ellipse becomes a polarisation ellipsoid.
A straightforward approach to analysing controlled source electromagnetic data is to resolve the observed fields into radial and azimuthal components.
However this suffers from contamination of field values by errors arising from small inaccuracies in source and receiver orientation and position. As an alternative, we therefore propose the novel approach of analysing the data from the recorded components in terms of polarisation ellipse properties - e.g. the amplitude and/or phase of the component of the signal aligned along the major axis of the ellipse. Since this value is much less dependent on the accuracy of the source and receiver positioning - and in particular removes the significant errors that arise from even small uncertainties in receiver orientation - use of polarisation ellipse values can lead to smaller uncertainties in measured field properties at the receiver. We show below as an example that use of the amplitude along the major axis of the polarisation ellipse for horizontal electric field can be used as a robust measurement parameter for detecting the presence of buried hydrocarbon layers. In principle, other polarisation ellipse parameters can be used in a similar way. Parameters that could readily be measured include amplitude or phase along the major axis; the horizontal polarisation ellipse of either electric or magnetic field; or in either case a polarisation ellipsoid including a vertical component.
Presentation of survey results using normalised values As seen in Figures 2A, 2B and 2C above, a feature of a controlled source electromagnetic survey over a sedimented seafloor is that the received field amplitude drops very rapidly with increasing range. The received signal properties also depend on the orientations of the source and receiving dipoles. In a typical survey, the signal amplitude is likely to vary by several orders of magnitude over the useful set of source-detector (transmitter-receiver) offsets; while the phase lag of the received signal also increases steadily with increasing offset. The effects of buried structure beneath the survey show up as deviations in the signal from that which would be expected if the sub-surface consisted only of a homogeneous, isotropic half space.
In order to visualise more clearly the effects of buried structure on controlled source electromagnetic signal properties, it is convenient to normalise the observed signals with respect to some reference model. A simple model to use would consist of a water layer of finite thickness above the survey of the true depth; and a homogeneous isotropic half space representing the seafloor below the survey.
If appropriate, and if the relevant a priori information is available, a more complex reference model may be used, although it is desirable to use always the simplest reference model that can broadly represent the large scale background properties of the seafloor.
The normalisation is carried out for amplitude by dividing the observed field by that calculated for the reference model, using the same source-detector geometry as for the observed data. In the case of phase, the normalisation is carried out by subtracting the phase lag calculated for the reference model from that of the observed data. We show below that use of normalised values calculated in this way can dramatically enhance the visual presentation of controlled source electromagnetic data from a survey over a hydrocarbon reservoir. Our illustrations use normalised amplitudes, but normalised phases could equally be used.
As an alternative, normalisation could be based on the survey data itself, for example using data collected adj acent to the target.
The varying~physical response of a survey de~endin~ on source orientation ~eometry Overall, the most crucial factor for the success or otherwise of the controlled source electromagnetic technique in practical applications related to hydrocarbon reservoirs is related to survey geometry.
Figure 3 is a schematic plan view from above illustrating a suitable co-ordinate system for describing the relative placement of a controlled source electromagnetic source 30 and an electromagnetic detector 32. The position of the detector 32 with respect to the source 30 is most suitably described in polar co-ordinates, with the centre of the source 30 providing the origin of the polar co-ordinate system.
The position of the detector 32 is defined by an azimuthal angle 8 and a distance R. The detector 32 and source 30 are considered to be co-planar. The angle 0 is measured clockwise from a line 33 passing through, and running parallel to, the source axis, as indicated in Figure 3 by the line marked 8=0°. A detector placed along the line 33, i.e.
such that is has an azimuthal angle 8 of 0°, is referred to as being in an end-on position. A detector with an azimuthal angle ~ of 90°, such that it lies on the line 34 marked A=90° in Figure 3, is referred to as being in a broadside position. The electric field strength at the detector is considered in terms of a radial component EP
and an orthogonal azimuthal component Ee as indicated in Figure 3.
Depending on the relative placements and orientations of sources and receivers, the physics of propagation of the signal through the seafloor - and so the net effect on the properties of the received signal - can be more or less sensitive to different classes of sub seafloor structure. As a simple illustration of this in the context of hydrocarbon surveys, we consider a case in which the sub-seafloor structure can be represented by a simple stack of horizontal layers (a useful first approximation for many situations in sedimentary basins).
Figure 4A shows in schematic vertical section an exemplary subterranean strata configuration. A section of seafloor 42 lies beneath 800 m of seawater 40. The strata beneath the seafloor 42 comprise a 1 km thick overburden layer 44, representing sediments above a hydrocarbon reservoir layer 46. This overburden layer 44 has low resistivity, primarily due to aqueous saturation of pore spaces. The middle layer 46 corresponds to a 100 m thick hydrocarbon reservoir and has a resistivity perhaps 100 times greater than the overburden 44. This is due to the presence of the non-conducting hydrocarbon occupying much of the pore spaces. Below the hydrocarbon reservoir layer 46 is a sub-reservoir region 48 which has low resistivity due to its similarity to the overburden layer 44 and extends downwards for an effectively infinite distance.
In the case of an electromagnetic survey, we can consider the differing behaviours of electric currents generated by the source resolved along horizontal and vertical directions. The relationships between the electric currents flowing in two adjacent regions of space are determined by both galvanic (i.e. direct transfer of charge) and inductive effects. Since charge is conserved, current leaving one volume of the sub-surface strata 44, 46, 48 and arriving in an adjacent volume along the direction of current flow are related to each other by a galvanic mechanism.
On the other hand, if two volumes are close together but separated from each other along a direction orthogonal to current flow, then the coupling between the currents flowing in the two volumes will be primarily inductive.
Applying this to our model of a sub-seafloor structure shown in Figure 4A, we can infer that the effect of the thin but resistive hydrocarbon reservoir layer 46 on the survey results will depend strongly on the direction of flow of the currents generated by the source. If the current at the base of the overburden layer 44 is dominantly horizontal, then coupling between this layer and the deeper layers 46, 48 will be dominated by inductive effects. Although little current will flow in the hydrocarbon reservoir layer 46, the inductive coupling between the overburden 44 and the underlying layer 48 will be only mildly affected. Currents in the two conducting layers 44, 48 (overburden and underlying layers) will therefore be similar to the case without a hydrocarbon reservoir layer 46. The results of such a survey will therefore be only weakly affected by the presence of the hydrocarbon reservoir layer 46.
In contrast, if our survey generates significant vertical components of electric current flow in the overburden layer 44, then galvanic effects along the current path will be strongly affected by the thin resistive hydrocarbon reservoir layer 46. The resistive hydrocarbon reservoir layer 46 will tend to block the current flow.
As a result, currents flowing in the deeper layer 48 will be much reduced; and the overall pattern of current flow in the overburden layer 44 will therefore be significantly altered. In such a survey situation, the presence of the hydrocarbon will significantly influence the outcome.
Hence it becomes of paramount importance in designing a practical survey approach for detecting buried hydrocarbon reservoirs to distinguish between source (transmitter) and detector (receiver) geometries in which coupling between layers is largely inductive between current sheets in a horizontal plane (in which case the survey has little sensitivity to the presence of the reservoir); and those in which a significant component of vertical current flow occurs, in which case blocking of the passage of this current flow by the reservoir leads to a survey which is strongly sensitive to the presence of the reservoir.
A specific approach to combining two contrasting source ~eometries when pros~ectin~ f, or hydrocarbon filled reservoirs By selecting a suitable survey geometry, it is possible to exploit the diverse properties of electromagnetic induction outlined above, by designing a survey in such a way that it provides data that are sensitive to the presence of thin resistive layers (exploiting largely galvanic effects and vertical components of induced current flow);
while simultaneously obtaining data that are more sensitive to the much larger scale structure, but relatively insensitive to a thin hydrocarbon layer (exploiting the dominantly inductive effects). The reason why the latter component is essential is that many features other than hydrocarbon reservoirs can affect the resistivity beneath the seafloor and the results of a survey. For example, resistivity often increases steadily with depth in submarine sedimentary basins, due to the progressive expulsion of pore fluids by rising overburden pressure.
Figure 4B shows in vertical section a highly schematic exemplary subterranean strata configuration which exhibits increasing resistivity with depth. A
section of seafloor 52 lies beneath 800 m of seawater 50. The strata beneath the seafloor comprise a series of layers of sediment of increasing resistivity. A first layer 54 has a uniform resistivity of 1 S2m and a thickness of 1 km. A second layer 56 has a uniform resistivity of 5 S2m and a thickness of 1 km. A third layer 58 has a uniform resistivity of 10 S2m and a thickness of 1 km. A fourth layer 60 has a uniform resistivity of 50 S2m and a thickness of 1 lmn. Beneath the fourth layer 60 is a fifth layer 62 which has a resistivity of 100 SZm and extends downwards for an infinite extent.
The increasing resistivity indicated in Figure 4B leads to longer skin depths, and hence to larger observed amplitudes and smaller observed phase lags for all detector placements. Although this does not exactly reproduce the same galvanic effects of a thin resistive layer, the overall effect on observed fields in the end-on geometry is likely to be very similar to (and difficult or impossible to distinguish from) the effect of a hydrocarbon reservoir. Designing a survey in such a way that different parts of the resulting data are more or less sensitive to the presence of certain key features is essential for removing ambiguity in the interpretation.
Recent studies of volcanic systems at mid-ocean ridges and of sediments beneath resistive basalt layers (MacGregor et al., 1998, 2001; MacGregor &
Sinha, 2000) have demonstrated that improved model resolution and reduced interpretational ambiguity can be achieved by using received electric field components from two distinct source-detector geometries. This approach is now applied to the case of surveying for thin hydrocarbon reservoirs.
Figures 4C, 4D, 4E and 4F show survey simulation results for the two subterranean strata configurations shown in Figures 4A and 4B. The method relies on collecting survey data using two distinct geometric modes. The first mode corresponds to an end-on geometry, in which, as described above, the receiver lies along the axis of the transmitting dipole. The source-detector azimuth B as defined in Figure 3 is 0°, and the field observed at the receiver is dominated by the radial component EP. The second mode corresponds to a broadside geometry, in which, as described above, the source-detector azimuth 8 is 90°, and the field observed at the receiver is dominated by the azimuthal component Ee.
For the case of hydrocarbon exploration, where the target is a thin resistive layer embedded between a more conductive underburden and overburden, such as that shown schematically in Figure 4A, the received signal depends on the azimuthal angle 0 of the detector. Because of the transmitter field pattern, the end-on, or EP, geometry results in a significant contribution to the observed field at the seafloor by the vertical component of current flow. The broadside, or Ee, geometry results in fields at the seafloor that are more dependent on the contribution of inductively coupled currents flowing in horizontal planes. As a result, the presence of a resistive hydrocarbon reservoir layer produces a significant increase at certain ranges in the end-on (EP) fields, while having virtually no effect on the broadside (Eo) fields. The result is 'splitting' of amplitudes between the two modes, and this splitting is characteristic of a structure in which resistivity first increases, in this case due to a hydrocarbon layer, and then decreases with depth.
Figure 4C is a graph showing the modelled amplitudes of electric field E along the major axes of the horizontal polarisation ellipses for two source-detector geometries and as a function of distance R between source and detector. These are determined for the model of sub-seafloor structure indicated in Figure 4A.
Four curves are shown. The curve labelled Eo shows the field strength as a function of distance for end-on geometries, the curve labelled Ego shows the same for broadside geometries.
The curves labelled No and N9o are the corresponding normalisation curves.
These are determined for the same geometries as for Eo and E9o, but represent the modelled amplitudes of electric field along the major axes of the horizontal polarisation ellipses where the subterranean strata comprise a homogeneous isotropic half space of resistivity 1 S2m, such as previously described.
Figure 4D is a graph showing the set of curves which correspond to those curves shown in Figure 4C, but for the subterranean strata configuration shown in Figure 4B.
Figure 4E is a graph showing the normalised values of the fields Eo and E9o shown in Figure 4C. These fields have been normalised by the fields No and N9o, also shown in Figure 4C, as discussed above.
Figure 4F is a graph showing the normalised values of the fields Eo and E9o shown in Figure 4D. These fields have been normalised by the fields No and N9o, also shown in Figure 4D, as discussed above.
It can be seen in Figure 4E that the presence of the~reservoir has relatively little effect on the broadside amplitudes (i.e. the curve labelled E9o/N9o is close to unity), but that between 4 and 6 km offset, it produces a substantial increase in the end-on amplitudes (i.e. the curve labelled Eo/No is strongly enhanced over this range). The hydrocarbon layer in this particular model leads to amplitude splitting between the two geometrical modes by about a factor of 10.
It can be seen in Figure 4F that the steadily increasing resistivity with depth model shown in Figure 4B strongly affects both geometries (i.e. both curves are strongly enhanced between 4 and 6 km offset) and no significant splitting between the two geometric modes occurs.
The end-on geometry data alone cannot distinguish between the two subterranean strata configurations shown in Figures 4A and 4B. These comparative examples conclusively demonstrate the fact that collection of survey data from both end-on and broadside geometries is needed, in order to reliably detect the hydrocarbon layer and eliminate "false-positive" detection of a rock formation containing no hydrocarbon reservoir, but merely having an increasing resistivity with depth, which is not unusual.
Lastly, the normalised curves in Figures 4E and 4F start off close to unity, at short ranges where the data are insensitive to the buried structure. The effects of the buried structure become greatest at between 5 and 6 km range. Beyond this, the normalised curves rapidly return towards a value of l, because at long ranges, air wave propagation (insensitive to sub-seafloor structure) begins to dominate, thus masking the effects of the sub surface. As a result, the presence of the target structure can best be observed in controlled source electromagnetic data of this type over a relatively narrow window of source-detector separation ranges.
An efficient electromagnetic survey may use a large number of detectors distributed over a target area and a signal transmitter which operates continuously as it is towed along an extended tow path. This follows from the relatively low cost of deploying detectors and the relatively high cost of deploying a transmitter.
Figure 5 is a schematic plan view showing an example layout of sixteen detectors 64 distributed across a section of seafloor 65 above subterranean reservoir 66. The reservoir 66 has a linear extent of several km and its boundary is indicated by a heavy line 67. In this example, the detectors 64 are uniformly distributed in a square-grid pattern so as to approximately cover the substantially square reservoir 66.
In performing a survey, a source (not shown) starts from the location marked 'A' in Figure 5 and is towed, whilst broadcasting continuously, along a path indicated by the broken line 68, the survey is completed when the source reaches the location marked 'B'. Data are continuously recorded by the detectors 64 throughout the towing process and the position of the source transmitter relative to the detector network is also logged.
During the towing process, each of the detectors 64 presents several different orientations to the source. For example, when the source is directly above the detector position D 1 and on the vertical section of the tow path, the detectors at positions D2 and D3 are at different ranges in an end-on position, the detectors at positions D4 and DS are at different ranges in a broadside position and the detector at position D6 is midway between. However, when the source later passes over the detector position D1 when on the horizontal section of the tow path, the detectors at positions D2 and D3 are now in a broadside position, and the detectors at position D4 and DS are in an end-on position. Thus, in the course of a survey, and in conjunction with the positional information of the source, data from the detectors 64 can be used to provide details of the signal transmission through the subterranean strata for a comprehensive range of distances and orientations between source and detector, each with varying galvanic and inductive contributions to the signal propagation. In this way a simple continuous towing of the source transmitter can provide a detailed survey which covers the extent of the subterranean reservoir 66.
Although the above example is based on a square detector grid, it will be understood that a wide variety of detector placements may be used, for example other high symmetry regular grids, such as triangular or rectangular, may be used.
In addition irregular grids may be used that have no high level of symmetry.
On the detection of the edges of a reservoir by this techniaue We showed above that our proposed arrangement of survey geometries can both detect a hydrocarbon reservoir, and distinguish it from a general increase in resistivity with depth. The modelling though was based only on 1-D
assumptions, i.e.
horizontal layers of infinite horizontal extent. By using higher dimensional modelling, we show below that the survey technique also works well in the case of reservoirs of finite extent; and indeed that it can be used as a reliable means of locating the edges of a reservoir structure.
Figure 6A shows a schematic section view of a model subterranean strata configuration containing a hydrocarbon reservoir layer of finite horizontal extent. A
section of seafloor 70 lies beneath a body of sea 72 which is 1.15 km deep and has a resistivity of 0.31 S2m. A planar hydrocarbon reservoir 74 is 0.15 km thick and its base is 1 km below the seafloor .70. The reservoir 74 has a resistivity of 100 Sam and extends infinitely out of the plane of Figure 6A, and semi-infinitely within the plane of Figure 6A. A uniform sedimentary structure 76 of resistivity 1 S2m surrounds the reservoir 74 and fills the remaining space beneath the seafloor 70. For the purposes of the model, two transmitter positions are considered. A first dipole transmitter 75 is aligned parallel to and directly above the edge 77 of the reservoir 74 on the seafloor 70. The dipole axis of the transmitter 75 is perpendicular to the plane of Figure 6A
and a cross is used to simply indicate its position rather than represent its structure. A
second transmitter 78 is also on the seafloor 70 and aligned parallel to the first transmitter 75. The second transmitter 78 is horizontally displaced from the first transmitter 75 by 1 km in a direction which positions it directly above the reservoir 74 as indicated in Figure 6A.
Figures 6B and 6C show the results of 2.5 D modelling (3-D source, 2-D
resistivity structure) of the exemplary subterranean strata configuration shown in Figure 6A for the first and second transmitter positions 75 and 78 respectively in order to investigate the effects of the finite areal extent of the resistive hydrocarbon reservoir 74. Figures 6B and 6C are plan views representing the normalised fields at the seafloor 70 as a function of position. The normalised field component shown in each case is the semi-major axis of the horizontal polarisation ellipse of electric field.
In Figure 6B the first transmitter 75 is at the origin of the local co-ordinate system and its dipole axis is parallel to the vertical axis in this system. In Figure 6C
the second transmitter 78 is at the origin of the local co-ordinate system and its dipole axis is parallel to the vertical axis in this system.
The normalised field amplitude is shown in a grey scale representation, while contour lines indicate the absolute value of amplitude. Figures 6B and 6C show that the reservoir structure causes substantial increases in the fields at positions close to end-on geometries, this is especially so when the source and detector are both over the reservoir, as in Figure 6C. If the source or detector are at or outside the edge of the reservoir, the amplitude increase is not seen. At detector locations close to broadside geometry, the increase in amplitude is again not seen, even if both source and detector are over the reservoir. Figures 6B and 6C illustrate that the amplitude splitting effect is observable over reservoirs of finite size, as well as over the 1-D
(layered) structures discussed previously.
The marked difference between Figures 6B and 6C shows that the observed fields are extremely sensitive to the location of the source (transmitter) with respect to the edge of the reservoir. Tn Figure 6B, the marked asymmetry of the fields about the y axis shows that - provided the source is over the reservoir - the observed fields are also highly sensitive to the location of the detector (receiver) with respect to the edge of the reservoir. Thus by careful siting of both source tow tracks and detectors relative to a suspected hydrocarbon bearing structure, the survey method described here can provide detailed information on the areal extent of the hydrocarbon, as well as on its existence or otherwise.
In addition to initially surveying for subterranean hydrocarbon reservoirs, the utility of the invention in detecting the edges of reservoirs makes it applicable to assessing changes in hydrocarbon reservoir content over a period of time. This is of _~g_ particular value, for example, where a hydrocarbon reservoir is being exploited. As water replaces the hydrocarbon which is drawn from the reservoir, regular follow up surveys can be used to determine the change in the distribution of the remaining hydrocarbon. Such surveys could therefore be used to characterise the evolution of the reservoir through time, and hence enable more efficient extraction and reservoir management.
Source and receiver locations - design considerations for real surveys The methodology described in this document can be readily applied to real surveys for hydrocarbon resources on continental margins. The required transmitter and receiver characteristics in terms of power, signal to noise ratios and operating parameters can be met by existing instrumental technology.
In order to plan a successful survey, the following factors are important:
An understanding of the need to collect observations over an appropriate set of source-detector separation ranges, and using both broadside and end-on surveying geometries;
A designed layout of the detector array and of source tow line locations and orientations optimised for the above;
A designed layout for the survey that also takes account of the need to locate the limits of the areal extent of any reservoir, based on the approach presented in Figures 6A, 6B and 6C;
Prior modelling to establish the range of transmission frequencies, and the set of source-detector separation ranges, that should be used, taking into consideration the expected resistivity structure of the seafloor, the depth of the taxget, the depth of the seafloor, and the influence of the air wave on the data collected.
Conclusions We have shown that the application of a specific set of methods to a controlled source electromagnetic survey of the seafloor can allow the technique to be successfully applied to the problem of detecting hydrocarbons beneath the seafloor.
Specifically, we have shown that:
The survey should be carried out using a combination of a mobile horizontal electric dipole source transmitter, equipped with a neutrally buoyant streamed antenna and operated just above the seafloor; and an array of electric and/or magnetic field sensing receivers placed statically on the seafloor.
The effects of propagation of the signal through the atmosphere (the 'air wave') are significant at high frequencies, at long source-detector separation distances, and in shallow water. This effect limits the applicability of the method by favouring detection of structures relatively shallowly situated beneath the seafloor, but in deep water; and limits the choice of frequencies and source-detector offsets for the survey.
Improvements in data analysis can be achieved by making use of field measurement parameters aligned along the major axis of the polarisation ellipse at the seafloor.
For purposes of data presentation and interpretation it is desirable to make use of field parameters that have been normalised with reference to an appropriate simplified model of the sub-seafloor structure.
In order to resolve the presence of any hydrocarbon, and to distinguish its effect on the data from other likely structures and reduce the ambiguity of interpretation, it is essential to collect survey data from both the end-on and broadside geometric configurations; and to analyse the data in terms of splitting between the two geometric modes.
2.5 D modelling shows that provided an appropriately designed array of receivers, and an appropriate set of transmitter tow tracks, are used, then the method can yield valuable information about the limits of the areal extent of any sub-seafloor hydrocarbon reservoir, as well as detecting its presence.
Finally it will be understood that the invention is equally applicable to surveying of freshwater, for example large lakes, so that references to seafloor, seawater etc. should not be regarded as limiting.
References:
1. Sinha, M.C., Patel, P.D., Unsworth, M.J., Owen, T.R.E. & MacCormack, M.R.G.
An active source electromagnetic sounding system for marine use. Mar. Geophys.
Res.,12, 1990, 59-68.
2. MacGregor, L.M., Constable, S.C. & Sinha, M.G. The RAMESSES experiment III: Controlled source electromagnetic sounding of the Reykjanes Ridge at 570 45' N. Geophysical Journal International, 135, 1998, 773-789.
3. MacGregor, L. M. & Sinha, M.G. Use of marine controlled source electromagnetic sounding for sub-basalt exploration. Geophysical Pf°ospecting, 48, 2000, 1091-1106.
4. MacGregor, L., Sinha, M. ~ Constable, S. Electrical resistivity structure of the Valu Fa Ridge, Lau Basin, from marine controlled source electromagnetic sounding. Geophys. J. Ir~t.,146, 217-236, 2001.
5. Evans, R. L., Sinha, M. C., Constable, S. C. & Unsworth, M. J. On the electrical nature of the axial melt zone at 130N on the East Pacific Rise. J. Geophys.
Res., 99, 1994, 577 - 588.
6. WO 00/13046 7. WO 01/57555
Claims (31)
1. An electromagnetic survey method for surveying an area previously identified as potentially containing a subsea hydrocarbon reservoir, comprising:
providing an electromagnetic source having a dipole axis and first and second detectors;
obtaining first and second survey data sets bar moving the electromagnetic source relative to each detector to collect data over a range of source-to-detector distances, wherein the first survey data set is obtained with the dipole axis of the electromagnetic source aligned end-on relative to the first detector so that the first survey data set is sensitive to resistive hydrocarbon layers exploiting largely galvanic effects, and wherein the second survey data set is obtained with the dipole axis of the electromagnetic source aligned broadside relative to the second detector so that the second survey data set is relatively insensitive to resistive hydrocarbon layers exploiting dominantly inductive effects.
providing an electromagnetic source having a dipole axis and first and second detectors;
obtaining first and second survey data sets bar moving the electromagnetic source relative to each detector to collect data over a range of source-to-detector distances, wherein the first survey data set is obtained with the dipole axis of the electromagnetic source aligned end-on relative to the first detector so that the first survey data set is sensitive to resistive hydrocarbon layers exploiting largely galvanic effects, and wherein the second survey data set is obtained with the dipole axis of the electromagnetic source aligned broadside relative to the second detector so that the second survey data set is relatively insensitive to resistive hydrocarbon layers exploiting dominantly inductive effects.
2. An electromagnetic survey method according to claim 1, wherein the first and second survey data sets are obtained concurrently.
3. An electromagnetic survey method according to claim 1 or claim 2, wherein a right angle is formed between a first line leading from the first detector to the electromagnetic source and a second line leading from the second detector to the electromagnetic source.
4. An electromagnetic survey method for surveying an area previously identified as potentially containing a subsea hydrocarbon reservoir, comprising:
providing an electromagnetic source having a dipole axis and a first detector;
obtaining first and second survey data sets by moving the electromagnetic source relative to the first detector to collect data over a range of source-to-detector distances and orientations, wherein the first survey data set is obtained with the dipole axis of the electromagnetic source aligned end-on relative to the first detector so that the first survey data set is sensitive to resistive hydrocarbon layers exploiting largely galvanic effects, and wherein the second survey data set is obtained with the dipole axis of the electromagnetic source aligned broadside relative to the first detector so that the second survey data set is relatively insensitive to resistive hydrocarbon layers exploiting dominantly inductive effects.
providing an electromagnetic source having a dipole axis and a first detector;
obtaining first and second survey data sets by moving the electromagnetic source relative to the first detector to collect data over a range of source-to-detector distances and orientations, wherein the first survey data set is obtained with the dipole axis of the electromagnetic source aligned end-on relative to the first detector so that the first survey data set is sensitive to resistive hydrocarbon layers exploiting largely galvanic effects, and wherein the second survey data set is obtained with the dipole axis of the electromagnetic source aligned broadside relative to the first detector so that the second survey data set is relatively insensitive to resistive hydrocarbon layers exploiting dominantly inductive effects.
5. An electromagnetic survey method according to claim 4, wherein the first and second survey data sets are obtained consecutively.
6. An electromagnetic survey method according to claim 4 or claim 5, wherein the electromagnetic source is moved in a direction along a line connecting the electromagnetic source to the detector to obtain other of the first and second survey data sets and again transverse thereto to obtain the other of the first and second survey data sets.
7. An electromagnetic survey method according to any one of claims 1 to 6, wherein the electromagnetic source is a mobile horizontal electric dipole source equipped with a streamed antenna which is towed above the seafloor to obtain the first and second survey data sets.
8. An electromagnetic survey method according to any one of claims 1 to 7, wherein each detector is static while the first and second survey data sets are obtained.
9. An electromagnetic survey method according to claim 8, wherein each detector is placed on the seafloor while the first and second survey data sets are obtained.
10. A method of analysing results from an electromagnetic survey of an area potentially containing a subsea hydrocarbon reservoir, comprising:
providing first and second survey data sets obtained using the method of any one of claims 1 to 9; and combining the first and second survey data sets to obtain a results data, set that represents a difference between the end on and broadside alignments as a function of the source-to-detector distances surveyed.
providing first and second survey data sets obtained using the method of any one of claims 1 to 9; and combining the first and second survey data sets to obtain a results data, set that represents a difference between the end on and broadside alignments as a function of the source-to-detector distances surveyed.
11. A method of analysing results from an electromagnetic survey according to claim 10, further comprising:
normalising each of the first and second survey data sets relative to respective first and second normalisation data sets or functions specific to the end-on and broadside alignments respectively, prior to the combining.
normalising each of the first and second survey data sets relative to respective first and second normalisation data sets or functions specific to the end-on and broadside alignments respectively, prior to the combining.
12. A method of analysing results from an electromagnetic survey according to claim 11, wherein the first and second normalisation data sets or functions are calculated from a rock formation model.
13. A method of analysing results from an electromagnetic survey according to claim 11, wherein the first and second normalisation data sets or functions are calculated from the first and second survey data sets.
14. A method of analysing results from an electromagnetic survey according to any one of claims 10 to 13, wherein the first and second data sets each comprise radial and azimuthal components of electric field or magnetic field measured at the detector, the method further comprising:
transforming the radial and azimuthal components into at least one polarisation ellipse parameter, prior to the combining.
transforming the radial and azimuthal components into at least one polarisation ellipse parameter, prior to the combining.
15. A method of analysing results from an electromagnetic survey according to claim 14, wherein the at least one polarisation ellipse parameter is the amplitude and/or phase of the component of the electric field or magnetic field aligned along a major axis of the ellipse.
16. A method of analysing results from an electromagnetic survey according to any one of claims 10 to 15, further comprising:
visually representing the results data set in a plot of at least two dimensions corresponding to the survey area.
visually representing the results data set in a plot of at least two dimensions corresponding to the survey area.
17. A method of analysing results from an electromagnetic survey according to claim 16, wherein the plot includes marking of areas of equal or similar electromagnetic held strength.
18. A method of analysing results from an electromagnetic survey according to claim 17, when appended to any one of claims 11 to 13, wherein the plot includes lines of equal absolute electromagnetic field strength.
19. A computer program product comprising a computer readable memory storing computer executable instructions thereon that when executed by a computer perform a method of analysing results from an electromagnetic survey according to any one of claims 10 to 18.
20. A computer readable memory having recorded thereon statements and instructions for execution by a computer to carry out the method of analysing results from an electromagnetic survey according to any one of claims 10 to 18.
21. A method of planning an electromagnetic survey of an area identified as potentially containing a subsea hydrocarbon reservoir, comprising:
creating a model of the area to be surveyed, including a rock formation containing a hydrocarbon reservoir and a body of water above the rock formation;
setting values for water depth, depth below the seafloor of the hydrocarbon reservoir, and resistivity structure of the rock formation;
performing a simulation of an electromagnetic survey in the model of the survey area by calculating first and second survey data sets by simulating an electromagnetic source having a dipole axis and first and second detectors and moving the electromagnetic source relative to each detector to collect data over a range of source-to-detector distances, wherein the first survey data set is obtained with the dipole axis of the electromagnetic source aligned end on relative to the first detector so that the first survey data set is sensitive to resistive hydrocarbon layers exploiting largely galvanic effects, and wherein the second survey data set is obtained with the dipole axis of the electromagnetic source aligned broadside relative to the second detector so that the second survey data set is relatively insensitive to resistive hydrocarbon layers exploiting dominantly inductive effects; and combining the first and second survey data sets to obtain a results data set that represents a difference between the end-on and broadside alignments as a function of the source-to-detector distances.
creating a model of the area to be surveyed, including a rock formation containing a hydrocarbon reservoir and a body of water above the rock formation;
setting values for water depth, depth below the seafloor of the hydrocarbon reservoir, and resistivity structure of the rock formation;
performing a simulation of an electromagnetic survey in the model of the survey area by calculating first and second survey data sets by simulating an electromagnetic source having a dipole axis and first and second detectors and moving the electromagnetic source relative to each detector to collect data over a range of source-to-detector distances, wherein the first survey data set is obtained with the dipole axis of the electromagnetic source aligned end on relative to the first detector so that the first survey data set is sensitive to resistive hydrocarbon layers exploiting largely galvanic effects, and wherein the second survey data set is obtained with the dipole axis of the electromagnetic source aligned broadside relative to the second detector so that the second survey data set is relatively insensitive to resistive hydrocarbon layers exploiting dominantly inductive effects; and combining the first and second survey data sets to obtain a results data set that represents a difference between the end-on and broadside alignments as a function of the source-to-detector distances.
22. A method of plug as electromagnetic survey of an area identified as potentially containing a subsea hydrocarbon reservoir, comprising:
creating a model of the area to be surveyed, including a rock formation containing a hydrocarbon reservoir and a body of water above the rock formation;
setting values for water depth, depth below the seafloor of the hydrocarbon reservoir, and resistivity structure of the rock formation;
performing a simulation of an electromagnetic survey in the model of the survey area by calculating first and second survey data sets by simulating an electromagnetic source having a dipole axis and a fist detector and moving the electromagnetic source relative to the first detector to collect data over a range of source-to-detector distances and orientations, wherein the first survey data set is obtained with the dipole axis of the electromagnetic source aligned end-on relative to the first detector so that the first survey data set is sensitive to resistive hydrocarbon layers exploiting largely galvanic effects, and wherein the second survey data set is obtained with the dipole axis of the electromagnetic source aligned broadside relative to the that detector sa that the second survey data set is relatively insensitive to resistive hydrocarbon layers exploiting dominantly inductive effects; and combining the first sad second survey data sets to obtain a results data set that represents a difference between the end-an and broadside alignments as a function of the source-to-detector distances.
creating a model of the area to be surveyed, including a rock formation containing a hydrocarbon reservoir and a body of water above the rock formation;
setting values for water depth, depth below the seafloor of the hydrocarbon reservoir, and resistivity structure of the rock formation;
performing a simulation of an electromagnetic survey in the model of the survey area by calculating first and second survey data sets by simulating an electromagnetic source having a dipole axis and a fist detector and moving the electromagnetic source relative to the first detector to collect data over a range of source-to-detector distances and orientations, wherein the first survey data set is obtained with the dipole axis of the electromagnetic source aligned end-on relative to the first detector so that the first survey data set is sensitive to resistive hydrocarbon layers exploiting largely galvanic effects, and wherein the second survey data set is obtained with the dipole axis of the electromagnetic source aligned broadside relative to the that detector sa that the second survey data set is relatively insensitive to resistive hydrocarbon layers exploiting dominantly inductive effects; and combining the first sad second survey data sets to obtain a results data set that represents a difference between the end-an and broadside alignments as a function of the source-to-detector distances.
23. A method of planning an electromagnetic survey according to claim 21 or 22, further comprising:
repeating the simulation for a number of source frequencies in order to select optimum surveying conditions in terms of source frequency and source-to-detector distances for probing the hydrocarbon reservoir.
repeating the simulation for a number of source frequencies in order to select optimum surveying conditions in terms of source frequency and source-to-detector distances for probing the hydrocarbon reservoir.
24. A method of planning an electromagnetic survey according to claim 21, 22 or 23, wherein tho modal includes a body of air above the body of water, and wherein the simulation takes account of signal propagation paths including the body of air when calculating the first and second survey data sets.
25. A method of planning an electromagnetic survey according to any one of claims 21 to 24, further comprising:
normalising each of the first and second survey data sets relative to respective first and second normalisation data seta or functions specific to the end-on and broadside alignments respectively, prior to the combining.
normalising each of the first and second survey data sets relative to respective first and second normalisation data seta or functions specific to the end-on and broadside alignments respectively, prior to the combining.
26. A method of planning an electromagnetic survey according to any one of claims 21 to 25, wherein the first and second data sets each comprise radial and azimuthal components of electric field or magnetic field measured at the detector, the method further comprising:
transforming the radial and azimuthal components into at least one polarisation ellipse parameter, prior to the combining.
transforming the radial and azimuthal components into at least one polarisation ellipse parameter, prior to the combining.
27. A method of planning an electromagnetic survey according to any one of claims 21 to 26, further comprising visually representing the results data sat in a plot of at least two dimensions corresponding to the survey area.
28. A method of planning an electromagnetic survey according to claim 27, wherein the plot includes markings of areas of equal or similar electromagnetic field.
29. A method of planning an electromagnetic survey according to claim 28, when appended to claim 25, wherein the plot includes lines of equal absolute electromagnetic field strength.
30. A computer program product comprising a computer readable memory storing computer executable instructions thereon that when executed by a computer perform the method of planning an electromagnetic survey according to any one of claims 21 to 29.
31. A computer readable memory having recorded thereon statements and instructions for execution by a computer to carry out the method of planning an electromagnetic survey according to any one of claims 21 to 29.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB0129407A GB2382875B (en) | 2001-12-07 | 2001-12-07 | Electromagnetic surveying for hydrocarbon reservoirs |
GB0129407.3 | 2001-12-07 | ||
PCT/GB2002/005355 WO2003048812A1 (en) | 2001-12-07 | 2002-11-28 | Electromagnetic surveying for hydrocarbon reservoirs |
Publications (2)
Publication Number | Publication Date |
---|---|
CA2468896A1 CA2468896A1 (en) | 2003-06-12 |
CA2468896C true CA2468896C (en) | 2013-02-19 |
Family
ID=9927254
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA2468896A Expired - Fee Related CA2468896C (en) | 2001-12-07 | 2002-11-28 | Electromagnetic surveying for hydrocarbon reservoirs |
Country Status (7)
Country | Link |
---|---|
US (2) | US7126338B2 (en) |
AU (1) | AU2002349125A1 (en) |
BR (2) | BRPI0214678B1 (en) |
CA (1) | CA2468896C (en) |
GB (1) | GB2382875B (en) |
NO (1) | NO325116B1 (en) |
WO (1) | WO2003048812A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20210341638A1 (en) * | 2019-09-05 | 2021-11-04 | Institute Of Geology And Geophysics Chinese Academy Of Sciences (Iggcas) | Method for Extracting IP Information from Grounded-wire Transient Electromagnetic Response |
Families Citing this family (123)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2382875B (en) * | 2001-12-07 | 2004-03-03 | Univ Southampton | Electromagnetic surveying for hydrocarbon reservoirs |
US6842006B2 (en) | 2002-06-27 | 2005-01-11 | Schlumberger Technology Corporation | Marine electromagnetic measurement system |
GB2395563B (en) | 2002-11-25 | 2004-12-01 | Activeem Ltd | Electromagnetic surveying for hydrocarbon reservoirs |
US7023213B2 (en) | 2002-12-10 | 2006-04-04 | Schlumberger Technology Corporation | Subsurface conductivity imaging systems and methods |
CN100339724C (en) | 2002-12-10 | 2007-09-26 | 加利福尼亚大学董事会 | System and method for hydrocarbon reservoir monitoring using controlled-source electromagnetic fields |
GB2402745B (en) | 2003-06-10 | 2005-08-24 | Activeem Ltd | Electromagnetic surveying for hydrocarbon reservoirs |
CA2526266C (en) * | 2003-06-26 | 2012-01-03 | Exxonmobil Upstream Research Company | Method for removing air wave effect from offshore frequency domain controlled-source electromagnetic data |
NO326506B1 (en) * | 2003-07-10 | 2008-12-15 | Norsk Hydro As | A marine geophysical collection system with a cable with seismic sources and receivers and electromagnetic sources and receivers |
GB2409900B (en) | 2004-01-09 | 2006-05-24 | Statoil Asa | Processing seismic data representing a physical system |
GB2411006B (en) * | 2004-02-16 | 2006-01-25 | Ohm Ltd | Electromagnetic surveying for hydrocarbon reservoirs |
GB2412741B (en) * | 2004-04-03 | 2009-02-25 | Statoil Asa | Electromagnetic data processing |
GB2413851B (en) | 2004-05-06 | 2006-08-09 | Ohm Ltd | Electromagnetic surveying for hydrocarbon reservoirs |
WO2006026361A1 (en) | 2004-08-25 | 2006-03-09 | The Regents Of The University Of California | Three-axis marine electric field sensor for seafloor electrical resistivity measurement |
JP3890344B2 (en) * | 2004-09-29 | 2007-03-07 | 株式会社日立ハイテクノロジーズ | Biomagnetic field measurement device |
GB2420855B (en) | 2004-12-02 | 2009-08-26 | Electromagnetic Geoservices As | Source for electromagnetic surveying |
NO323490B1 (en) * | 2005-01-10 | 2007-05-29 | Norsk Hydro As | Electromagnetic hydrocarbon exploration in shallow sea areas |
GB2422673B (en) * | 2005-02-01 | 2010-03-24 | Electromagnetic Geoservices As | Optimum signal for sea bed logging |
US20060186887A1 (en) * | 2005-02-22 | 2006-08-24 | Strack Kurt M | Method for identifying subsurface features from marine transient controlled source electromagnetic surveys |
GB2423370B (en) | 2005-02-22 | 2007-05-02 | Ohm Ltd | Electromagnetic surveying for resistive or conductive bodies |
US7333893B2 (en) * | 2005-03-07 | 2008-02-19 | Exxonmobil Upstream Research Company | Method for spatially interpreting electromagnetic data using multiple frequencies |
US7383132B2 (en) | 2005-03-07 | 2008-06-03 | Exxonmobil Upstream Research Co. | Method for identifying resistivity anomalies in electromagnetic survey data |
US7295013B2 (en) | 2005-04-11 | 2007-11-13 | Schlumberger Technology Corporation | Remotely operable measurement system and method employing same |
EA010950B1 (en) * | 2005-06-09 | 2008-12-30 | Эксонмобил Апстрим Рисерч Компани | Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys |
WO2006135568A2 (en) * | 2005-06-10 | 2006-12-21 | Exxonmobil Upstream Research Company | Method for controlled source electromagnetic reconnaissance surveying |
RU2005118534A (en) * | 2005-06-16 | 2006-12-27 | Евгений Дмитриевич Лисицын (RU) | METHOD FOR MARINE ELECTRIC EXPLORATION OF OIL AND GAS DEPOSITS AND THE EQUIPMENT COMPLEX FOR ITS IMPLEMENTATION OF "VESOTEM" |
WO2007018810A1 (en) * | 2005-07-22 | 2007-02-15 | Exxonmobil Upstream Research Company | Method for determining receiver orientations |
CA2625709A1 (en) | 2005-10-14 | 2007-04-26 | Exxonmobil Upstream Research Company | Method and apparatus for utilizing time division multiple waveform transmitting |
NO323889B3 (en) | 2005-11-03 | 2007-07-16 | Advanced Hydrocarbon Mapping As | Method for mapping hydrocarbon reservoirs and apparatus for use in carrying out the method |
GB0526303D0 (en) * | 2005-12-23 | 2006-02-01 | Wireless Fibre Systems Ltd | Transmission of underwater electromagnetic radiation through the seabed |
GB2434340B (en) * | 2006-01-20 | 2008-01-02 | Ohm Ltd | Underwater equipment recovery |
US7203599B1 (en) | 2006-01-30 | 2007-04-10 | Kjt Enterprises, Inc. | Method for acquiring transient electromagnetic survey data |
GB2434868B (en) | 2006-02-06 | 2010-05-12 | Statoil Asa | Method of conducting a seismic survey |
GB2435693A (en) | 2006-02-09 | 2007-09-05 | Electromagnetic Geoservices As | Seabed electromagnetic surveying |
NO326957B1 (en) * | 2006-02-13 | 2009-03-23 | Norsk Hydro As | Electromagnetic method on shallow water using controlled source |
US8165815B2 (en) | 2006-02-21 | 2012-04-24 | Exxonmobil Upstream Research Co. | Method for electromagnetic air-wave suppression by active cancellation and shielding |
GB0604829D0 (en) * | 2006-03-10 | 2006-04-19 | Mtem Ltd | Optimisation of mtem parameters |
US20070216416A1 (en) * | 2006-03-15 | 2007-09-20 | Baker Hughes Incorporated | Electromagnetic and Magnetostatic Shield To Perform Measurements Ahead of the Drill Bit |
US7328107B2 (en) * | 2006-04-28 | 2008-02-05 | Kjt Enterprises, Inc. | Integrated earth formation evaluation method using controlled source electromagnetic survey data and seismic data |
GB2438430B (en) * | 2006-05-22 | 2008-09-17 | Ohm Ltd | Electromagnetic surveying |
GB2439378B (en) | 2006-06-09 | 2011-03-16 | Electromagnetic Geoservices As | Instrument for measuring electromagnetic signals |
US7340348B2 (en) * | 2006-06-15 | 2008-03-04 | Kjt Enterprises, Inc. | Method for acquiring and interpreting seismoelectric and electroseismic data |
US7356411B1 (en) * | 2006-07-01 | 2008-04-08 | Kjt Enterprises, Inc. | Method for acquiring and interpreting transient electromagnetic measurements |
US7860655B2 (en) * | 2006-07-14 | 2010-12-28 | Westerngeco L.L.C. | Electromagnetically detecting thin resistive bodies in shallow water and terrestrial environments |
US7657391B2 (en) * | 2006-07-14 | 2010-02-02 | Westerngeco L.L.C. | Electromagnetically detecting thin resistive bodies in shallow water and terrestrial environments |
US8121789B2 (en) | 2006-07-25 | 2012-02-21 | Exxonmobil Upstream Research Co. | Method for correcting the phase of electromagnetic data |
US7574410B2 (en) * | 2006-08-22 | 2009-08-11 | Kjt Enterprises, Inc. | Fast 3D inversion of electromagnetic survey data using a trained neural network in the forward modeling branch |
AU2007288399B2 (en) * | 2006-08-24 | 2011-09-29 | Exxonmobil Upstream Research Company | Interpretation and mapping of electromagnetic survey data |
GB0616870D0 (en) * | 2006-08-25 | 2006-10-04 | Mtem Ltd | Improvements In Marine EM Exploration |
US8253418B2 (en) * | 2006-08-30 | 2012-08-28 | The Regents Of The University Of California | Method and system for detecting and mapping hydrocarbon reservoirs using electromagnetic fields |
CA2663662C (en) | 2006-09-13 | 2016-07-05 | Exxonmobil Upstream Research Company | Rapid inversion of electromagnetic reconnaissance survey data |
US7427862B2 (en) * | 2006-09-29 | 2008-09-23 | Baker Hughes Incorporated | Increasing the resolution of electromagnetic tools for resistivity evaluations in near borehole zones |
US7400977B2 (en) * | 2006-10-12 | 2008-07-15 | Schlumberger Technology Corporation | Computing values for surveying a subterranean structure based on measurements according to different electromagnetic survey techniques |
GB2442749B (en) | 2006-10-12 | 2010-05-19 | Electromagnetic Geoservices As | Positioning system |
ITBO20060811A1 (en) * | 2006-11-29 | 2008-05-30 | Pierburg Spa | OIL PUMP WITH VARIABLE DISPLACEMENT PALETTE. |
US9110183B2 (en) | 2006-12-06 | 2015-08-18 | Technoimaging, Llc | Systems and methods for remote electromagnetic exploration for mineral and energy resources using stationary long-range transmitters |
US8564296B2 (en) | 2006-12-06 | 2013-10-22 | Technoimaging, Llc. | Systems and methods for remote electromagnetic exploration for mineral and energy resources |
US7969152B2 (en) * | 2006-12-06 | 2011-06-28 | Technoimaging, Llc | Systems and methods for measuring sea-bed resistivity |
GB2445582A (en) | 2007-01-09 | 2008-07-16 | Statoil Asa | Method for analysing data from an electromagnetic survey |
GB0702661D0 (en) * | 2007-02-12 | 2007-03-21 | Mtem Ltd | Improvements in marine EM exploration |
AU2008244598B8 (en) | 2007-04-26 | 2012-10-11 | Exxonmobil Upstream Research Company | Method for electromagnetic survey design |
US7640110B2 (en) * | 2007-04-27 | 2009-12-29 | Schlumberger Technology Corporation | Pixel based inversion method for surface electromagnetic measurement |
US7746077B2 (en) * | 2007-04-30 | 2010-06-29 | Kjt Enterprises, Inc. | Method for measuring the magnetotelluric response to the earth's subsurface |
US8026723B2 (en) | 2007-04-30 | 2011-09-27 | Kjt Enterprises, Inc. | Multi-component marine electromagnetic signal acquisition method |
US7872477B2 (en) * | 2007-04-30 | 2011-01-18 | Kjt Enterprises, Inc. | Multi-component marine electromagnetic signal acquisition cable and system |
NZ581872A (en) | 2007-05-14 | 2012-07-27 | Ocean Floor Geophysics Inc | Underwater electric field electromagnetic prospecting system |
US8294466B2 (en) * | 2007-06-14 | 2012-10-23 | Exxonmobil Upstream Research Company | Scaled plots of electromagnetic data |
GB2450158B (en) | 2007-06-15 | 2009-06-03 | Ohm Ltd | Electromagnetic detector for marine surveying |
US8482287B2 (en) * | 2007-07-03 | 2013-07-09 | Shell Oil Company | Method of determining electrical anisotropy in a subsurface formation |
US7705599B2 (en) * | 2007-07-09 | 2010-04-27 | Kjt Enterprises, Inc. | Buoy-based marine electromagnetic signal acquisition system |
US7565245B2 (en) | 2007-09-20 | 2009-07-21 | Ohm Limited | Electromagnetic surveying |
RU2351958C1 (en) * | 2007-11-12 | 2009-04-10 | Екатерина Николаевна Рыхлинская | Method of sea geo-electro-survey with electrical focusing (versions) |
US7949470B2 (en) * | 2007-11-21 | 2011-05-24 | Westerngeco L.L.C. | Processing measurement data in a deep water application |
US8547783B2 (en) | 2007-12-12 | 2013-10-01 | Exxonmobil Upstream Research Company | Method and apparatus for evaluating submarine formations |
WO2009079355A1 (en) * | 2007-12-18 | 2009-06-25 | Schlumberger Canada Limited | System and method for improving surface electromagnetic surveys |
US8275592B2 (en) * | 2008-04-07 | 2012-09-25 | Westerngeco L.L.C. | Joint inversion of time domain controlled source electromagnetic (TD-CSEM) data and further data |
US20090265111A1 (en) * | 2008-04-16 | 2009-10-22 | Kjt Enterprises, Inc. | Signal processing method for marine electromagnetic signals |
US20110066379A1 (en) * | 2008-05-26 | 2011-03-17 | Mes Marius J | survey system for locating geophysical anomalies |
US9383475B2 (en) | 2008-06-09 | 2016-07-05 | Rock Solid Images, Inc. | Geophysical surveying |
US8063642B2 (en) * | 2008-06-11 | 2011-11-22 | Mtem Ltd | Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources |
CA2725922C (en) * | 2008-06-24 | 2016-08-30 | Exxonmobil Upstream Research Company | Method for electromagnetic prospecting waveform design |
US8239181B2 (en) | 2008-07-23 | 2012-08-07 | Exxonmobil Upstream Research Company | Inversion of CSEM data with measurement system signature suppression |
GB2462861B (en) * | 2008-08-22 | 2010-08-18 | Ohm Ltd | Electromagnetic surveying |
GB2466764B (en) | 2008-10-02 | 2013-03-27 | Electromagnetic Geoservices As | Method for enhanced subsurface electromagnetic sensitivity |
US8164340B2 (en) | 2008-10-23 | 2012-04-24 | Kjt Enterprises, Inc. | Method for determining electromagnetic survey sensor orientation |
RU2381531C1 (en) * | 2008-12-22 | 2010-02-10 | Екатерина Николаевна Рыхлинская | Method for marine geoelectrical exploration with electrical current focusing |
US8258791B2 (en) | 2009-01-27 | 2012-09-04 | Mtem Ltd. | Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources to impart electromagnetic signals into a subsurface formation and thereby determining a formation response to each signal |
WO2010105206A2 (en) * | 2009-03-13 | 2010-09-16 | Schlumberger Canada Limited | Electromagnetic surface-to-borehole look around systems and methods of monitoring in horizontal wells |
US20100235100A1 (en) * | 2009-03-16 | 2010-09-16 | Bruce Alan Hobbs | Method for determining resistivity anisotropy from earth electromagnetic responses |
NO330702B1 (en) * | 2009-03-20 | 2011-06-14 | Advanced Hydrocarbon Mapping As | Method and apparatus for electromagnetic mapping of subsea hydrocarbon deposits based on total magnetic field paints |
US9035657B2 (en) * | 2009-04-10 | 2015-05-19 | Schlumberger Technology Corporation | Electromagnetic logging between a cased borehole and surface |
US8554482B2 (en) * | 2009-05-05 | 2013-10-08 | Baker Hughes Incorporated | Monitoring reservoirs using array based controlled source electromagnetic methods |
US20110012601A1 (en) * | 2009-07-15 | 2011-01-20 | Bruce Alan Hobbs | Method for determining resistivity anisotropy from earth electromagnetic tansient step response and electromagnetic transient peak impulse response |
US8386221B2 (en) * | 2009-12-07 | 2013-02-26 | Nuovo Pignone S.P.A. | Method for subsea equipment subject to hydrogen induced stress cracking |
US20110144472A1 (en) * | 2009-12-11 | 2011-06-16 | Technoimaging, Llc | Methods of gravity and/or magnetic holographic imaging using vector and/or tensor data |
US9588250B2 (en) | 2010-04-14 | 2017-03-07 | Baker Hughes Incorporated | Three-coil system with short nonconductive inserts for transient MWD resistivity measurements |
US8498845B2 (en) | 2010-04-21 | 2013-07-30 | Exxonmobil Upstream Research Company | Method for geophysical imaging |
GB2481845B (en) | 2010-07-08 | 2014-04-30 | Electromagnetic Geoservices As | Low noise marine electric field sensor system |
CA2802722C (en) | 2010-07-27 | 2023-04-04 | Exxonmobil Upstream Research Company | Inverting geophysical data for geological parameters or lithology |
US8624969B2 (en) | 2010-08-02 | 2014-01-07 | Technoimaging, Llc. | Methods of electromagnetic migration imaging of geologic formation |
US9195783B2 (en) | 2010-08-16 | 2015-11-24 | Exxonmobil Upstream Research Company | Reducing the dimensionality of the joint inversion problem |
US8756017B2 (en) * | 2011-02-17 | 2014-06-17 | Yangtze University | Method for detecting formation resistivity outside of metal casing using time-domain electromagnetic pulse in well |
US9110195B2 (en) | 2011-04-14 | 2015-08-18 | Wen J. Whan | Electromagnetic and its combined surveying apparatus and method |
US9453929B2 (en) | 2011-06-02 | 2016-09-27 | Exxonmobil Upstream Research Company | Joint inversion with unknown lithology |
US9702995B2 (en) | 2011-06-17 | 2017-07-11 | Exxonmobil Upstream Research Company | Domain freezing in joint inversion |
US9322910B2 (en) | 2011-07-15 | 2016-04-26 | Technoimaging, Llc | Method of real time subsurface imaging using electromagnetic data acquired from moving platforms |
WO2013012470A1 (en) | 2011-07-21 | 2013-01-24 | Exxonmobil Upstream Research Company | Adaptive weighting of geophysical data types in joint inversion |
US9606256B2 (en) * | 2011-12-07 | 2017-03-28 | Pgs Geophysical As | Method and system of determining parameters associated with a hydrocarbon bearing formation beneath a sea bed |
US8587316B2 (en) | 2011-12-08 | 2013-11-19 | Pgs Geophysical As | Noise reduction systems and methods for a geophysical survey cable |
US8922214B2 (en) | 2011-12-27 | 2014-12-30 | Pgs Geophysical As | Electromagnetic geophysical survey systems and methods employing electric potential mapping |
US8898102B2 (en) | 2012-05-25 | 2014-11-25 | Pgs Geophysical As | Adaptive computational grid for inverse problems |
WO2014039851A1 (en) * | 2012-09-06 | 2014-03-13 | Mastinc. | System and method to monitor powerlines |
US10591638B2 (en) | 2013-03-06 | 2020-03-17 | Exxonmobil Upstream Research Company | Inversion of geophysical data on computer system having parallel processors |
US20140266215A1 (en) * | 2013-03-12 | 2014-09-18 | Pgs Geophysical As | Systems and methods for measuring water properties in electromagnetic marine surveys |
US9274241B2 (en) * | 2013-03-14 | 2016-03-01 | Pgs Geophysical As | Method and system for suppressing swell-induced electromagnetic noise |
US9846255B2 (en) | 2013-04-22 | 2017-12-19 | Exxonmobil Upstream Research Company | Reverse semi-airborne electromagnetic prospecting |
US9316758B2 (en) | 2013-05-29 | 2016-04-19 | Liquid Robotics Oil and Gas LLC | Earth surveying for improved drilling applications |
US9890632B2 (en) | 2014-06-20 | 2018-02-13 | Saudi Arabian Oil Company | Systems, computer medium and computer-implemented methods for logging using a logging tool with adjustable detectors |
CN104237970B (en) * | 2014-09-23 | 2017-07-07 | 中国石油天然气集团公司 | Electromagnetism of Earthquake joint exploration system and its data acquisition device and collecting method |
WO2016120814A2 (en) * | 2015-01-28 | 2016-08-04 | University Of The Witwatersrand, Johannesburg | Method of determining the location and depth of magnetic sources in the earth |
US10571592B2 (en) | 2015-08-31 | 2020-02-25 | Pgs Geophysical As | Direct resistivity determination |
NO342689B1 (en) * | 2016-05-30 | 2018-07-09 | Advanced Hydrocarbon Mapping As | Apparatus for orienting an electromagnetic field sensor, and related receiver unit and method |
US9863895B1 (en) | 2017-02-22 | 2018-01-09 | Saudi Arabian Oil Company | Systems and methods for monitoring casing cement integrity |
JP7148266B2 (en) * | 2018-04-26 | 2022-10-05 | 川崎重工業株式会社 | Work method using autonomous unmanned submersible |
Family Cites Families (87)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US1818331A (en) * | 1928-08-21 | 1931-08-11 | Radiore Company | Method for determining the character of ore bodies |
US2077707A (en) * | 1933-08-01 | 1937-04-20 | Melton Benjamin Starr | Electromagnetic prospecting method |
US2139460A (en) * | 1936-07-06 | 1938-12-06 | Geo Frequenta Corp | Means and method for geophysical prospecting |
US2268106A (en) * | 1939-05-13 | 1941-12-30 | Standard Oil Dev Co | Radio wave prospecting |
US2426918A (en) * | 1941-03-17 | 1947-09-02 | Engineering Res Corp | Method for electromagnetic-wave investigations of earth formations |
US2531088A (en) * | 1947-10-16 | 1950-11-21 | Standard Oil Dev Co | Electrical prospecting method |
US2766422A (en) * | 1949-07-02 | 1956-10-09 | Carbonetto Stefano | Methods and arrangements for detecting layers of hydrocarbons in the ground |
US2919397A (en) * | 1956-04-30 | 1959-12-29 | Lawrence W Morley | Method and apparatus for inductive prospecting |
US2953742A (en) * | 1957-09-04 | 1960-09-20 | Charles J Hughes | Geophysical prospecting apparatus |
US3052836A (en) * | 1957-12-24 | 1962-09-04 | Shell Oil Co | Method for marine electrical prospecting |
US3398356A (en) * | 1964-02-10 | 1968-08-20 | Westinghouse Electric Corp | Method utilizing a pair of subsurface antennas for determining the physical properties effecting radio energy propagation through earth |
GB1239953A (en) | 1967-06-06 | 1971-07-21 | Rech S Geol Et Minieres Bureau | Improvements in or relating to methods and apparatus for determining the electrical resistance of the sub-soil |
US4010413A (en) * | 1971-08-23 | 1977-03-01 | Geo-Nav, Inc. | Plural frequency geological exploration system and method with phase comparison |
US3806795A (en) * | 1972-01-03 | 1974-04-23 | Geophysical Survey Sys Inc | Geophysical surveying system employing electromagnetic impulses |
FR2288988A1 (en) * | 1974-07-30 | 1976-05-21 | Duroux Jean | METHOD AND APPARATUS FOR PROSPECTING AT SEA BY MEASURING ELECTROMAGNETIC FIELDS |
DE2550715C3 (en) * | 1975-11-12 | 1980-07-03 | Diether-Alfred 5305 Alfter- Oedekoven Schroeder | Circuit arrangement for geophysical prospecting out of a vehicle |
US4079309A (en) * | 1976-09-03 | 1978-03-14 | The United States Of America As Represented By The Secretary Of The Navy | Method for determining changes in earth resistivity by measuring phase difference between magnetic field components |
FR2390743A1 (en) * | 1977-05-09 | 1978-12-08 | Geophysique Cie Gle | ELECTROMAGNETIC PROSPECTION OF THE BASEMENT BY INDUCTION, ASSOCIATED WITH PROSPECTION BY ELECTRIC SURVEY |
US4258321A (en) * | 1978-03-09 | 1981-03-24 | Neale Jr Dory J | Radio geophysical surveying method and apparatus |
US4308499A (en) * | 1978-05-26 | 1981-12-29 | Kali Und Salz A.G. | Method utilizing electromagnetic wave pulses for determining the locations of boundary surfaces of underground mineral deposits |
US5025218A (en) * | 1979-04-23 | 1991-06-18 | The United States Of America As Represented By The Secretary Of The Navy | Pulsed field system for detecting the presence of a target in a subsurface environment |
US4218678A (en) * | 1979-05-11 | 1980-08-19 | Ensco, Inc. | Synthetic pulse radar including a microprocessor based controller |
MA18895A1 (en) * | 1979-07-09 | 1981-04-01 | Cie Generale De Geophysique Sa | METHOD AND DEVICE FOR GEOPHYSICAL PROSPECTION WITH TRANSIENT CURRENTS |
FR2497360A1 (en) * | 1980-12-31 | 1982-07-02 | Schlumberger Prospection | PHASE AND AMPLITUDE MEASUREMENT FOR A DIELECTRIC PROPERTIES LOGGING SYSTEM |
US4451789A (en) * | 1981-09-28 | 1984-05-29 | Nl Industries, Inc. | Logging tool and method for measuring resistivity of different radial zones at a common depth of measurement |
US4504833A (en) * | 1981-12-09 | 1985-03-12 | Xadar Corporation | Synthetic pulse radar system and method |
US4506225A (en) * | 1981-12-28 | 1985-03-19 | Barringer Research Limited | Method for remote measurement of anomalous complex variations of a predetermined electrical parameter in a target zone |
US4489276A (en) * | 1982-01-20 | 1984-12-18 | The United States Of America As Represented By The United States Department Of Energy | Dual-cone double-helical downhole logging device |
PL141895B1 (en) | 1983-03-03 | 1987-09-30 | Instytut Gornictwa Naftowego Gaz | Method of and system for direct prospecting of hydrocarbon accumulations |
US4594551A (en) * | 1983-03-31 | 1986-06-10 | Texaco Inc. | Method of deep penetration well logging using three receivers |
US4617518A (en) * | 1983-11-21 | 1986-10-14 | Exxon Production Research Co. | Method and apparatus for offshore electromagnetic sounding utilizing wavelength effects to determine optimum source and detector positions |
US4616184A (en) * | 1984-06-27 | 1986-10-07 | The United States Of America As Represented By The United States Department Of Energy | CSAMT method for determining depth and shape of a sub-surface conductive object |
US4652829A (en) * | 1984-12-28 | 1987-03-24 | Schlumberger Technology Corp. | Electromagnetic logging apparatus with button antennas for measuring the dielectric constant of formation surrounding a borehole |
DE3529466A1 (en) * | 1985-08-16 | 1987-04-09 | Pipeline Engineering Ges Fuer | METHOD FOR DETERMINING THE LIMITS OF UNDERGROUND NATURAL GAS DEPOSIT |
US5633590A (en) * | 1986-11-04 | 1997-05-27 | Paramagnetic Logging, Inc. | Formation resistivity measurements from within a cased well used to quantitatively determine the amount of oil and gas present |
US5570024A (en) * | 1986-11-04 | 1996-10-29 | Paramagnetic Logging, Inc. | Determining resistivity of a formation adjacent to a borehole having casing using multiple electrodes and with resistances being defined between the electrodes |
US4835474A (en) * | 1986-11-24 | 1989-05-30 | Southwest Research Institute | Method and apparatus for detecting subsurface anomalies |
GB8825435D0 (en) * | 1988-10-31 | 1988-12-29 | Cross T E | Detection of non metallic material |
US5025222A (en) * | 1989-07-07 | 1991-06-18 | Phase Dynamics, Inc. | System and method for monitoring substances and reactions |
US5066916A (en) * | 1990-01-10 | 1991-11-19 | Halliburton Logging Services, Inc. | Technique for separating electromagnetic refracted signals from reflected signals in down hole electromagnetic tools |
US5185578A (en) * | 1990-01-17 | 1993-02-09 | Stolar, Inc. | Method for detecting anomalous geological zones by transmitting electromagnetic energy between spaced drillholes using different frequency ranges |
US5877995A (en) * | 1991-05-06 | 1999-03-02 | Exxon Production Research Company | Geophysical prospecting |
US5192952A (en) * | 1991-06-11 | 1993-03-09 | Johler J Ralph | Method and apparatus for transmitting electromagnetic signals into the earth from a capacitor |
US5280284A (en) * | 1991-06-11 | 1994-01-18 | Johler J Ralph | Method of determining the electrical properties of the earth by processing electromagnetic signals propagated through the earth from a capacitor |
US5274386A (en) * | 1992-06-17 | 1993-12-28 | General Electric Co. | Reduced hardware antenna beamformer |
USH1490H (en) * | 1992-09-28 | 1995-09-05 | Exxon Production Research Company | Marine geophysical prospecting system |
US5486764A (en) * | 1993-01-15 | 1996-01-23 | Exxon Production Research Company | Method for determining subsurface electrical resistance using electroseismic measurements |
USH1524H (en) * | 1993-01-15 | 1996-04-02 | Exxon Production Research Company | Method for using electromagnetic grounded antennas as directional geophones |
US5357253A (en) * | 1993-04-02 | 1994-10-18 | Earth Sounding International | System and method for earth probing with deep subsurface penetration using low frequency electromagnetic signals |
US5373443A (en) * | 1993-10-06 | 1994-12-13 | The Regents, University Of California | Method for imaging with low frequency electromagnetic fields |
US6060885A (en) * | 1993-10-14 | 2000-05-09 | Western Atlas International, Inc. | Method and apparatus for determining the resistivity and conductivity of geological formations surrounding a borehole |
US5563513A (en) * | 1993-12-09 | 1996-10-08 | Stratasearch Corp. | Electromagnetic imaging device and method for delineating anomalous resistivity patterns associated with oil and gas traps |
US5400030A (en) * | 1994-02-09 | 1995-03-21 | Exxon Production Research Company | Detection and mapping of hydrocarbon reservoirs with radar waves |
US6060884A (en) * | 1994-03-14 | 2000-05-09 | Meyer, Jr.; Wallace Harold | Method and apparatus for measuring electromagnetic properties of materials in borehole environs and simultaneously determining the quality of the measurements |
US5811973A (en) * | 1994-03-14 | 1998-09-22 | Baker Hughes Incorporated | Determination of dielectric properties with propagation resistivity tools using both real and imaginary components of measurements |
US5892361A (en) * | 1994-03-14 | 1999-04-06 | Baker Hughes Incorporated | Use of raw amplitude and phase in propagation resistivity measurements to measure borehole environmental parameters |
NO314646B1 (en) * | 1994-08-15 | 2003-04-22 | Western Atlas Int Inc | Transient electromagnetic measuring tool and method for use in a well |
USH1561H (en) * | 1994-09-22 | 1996-07-02 | Exxon Production Research Company | Method and apparatus for detection of seismic and electromagnetic waves |
FR2729222A1 (en) * | 1995-01-10 | 1996-07-12 | Commissariat Energie Atomique | DETERMINATION OF THE POROSITY AND PERMEABILITY OF A GEOLOGICAL FORMATION FROM THE ELECTROFILTRATION PHENOMENON |
DE19518420C2 (en) * | 1995-05-19 | 1998-01-02 | Diether Alfred Schroeder | Circuit arrangement for use in a geophysical prospecting process |
US5553407A (en) * | 1995-06-19 | 1996-09-10 | Vermeer Manufacturing Company | Excavator data acquisition and control system and method of use |
US6023168A (en) * | 1995-08-21 | 2000-02-08 | Schlumberger Technology Corporation | Apparatus and method for measuring the resistivity of underground formations |
GB9521171D0 (en) * | 1995-10-17 | 1995-12-20 | Millar John W A | Detection method |
US5886526A (en) * | 1996-06-19 | 1999-03-23 | Schlumberger Technology Corporation | Apparatus and method for determining properties of anisotropic earth formations |
US5770945A (en) * | 1996-06-26 | 1998-06-23 | The Regents Of The University Of California | Seafloor magnetotelluric system and method for oil exploration |
US5777478A (en) * | 1996-08-27 | 1998-07-07 | John R. Jackson | Passive geophysical prospecting apparatus and method based upon detection of discontinuities associated with extremely low frequency electromagnetic fields |
US5825188A (en) * | 1996-11-27 | 1998-10-20 | Montgomery; Jerry R. | Method of mapping and monitoring groundwater and subsurface aqueous systems |
US5841280A (en) * | 1997-06-24 | 1998-11-24 | Western Atlas International, Inc. | Apparatus and method for combined acoustic and seismoelectric logging measurements |
US6188222B1 (en) * | 1997-09-19 | 2001-02-13 | Schlumberger Technology Corporation | Method and apparatus for measuring resistivity of an earth formation |
US6140819A (en) * | 1998-05-26 | 2000-10-31 | Heath Consultants, Inc. | Continuous-depth-indicating underground pipe and cable locator |
NO315725B1 (en) * | 1998-06-18 | 2003-10-13 | Norges Geotekniske Inst | Device for measuring and monitoring resistivity outside a well pipe in a petroleum reservoir |
US6236212B1 (en) * | 1998-06-22 | 2001-05-22 | The United States Of America As Represented By The Secretary Of The Interior | Induced polarization system using towed cable carrying transmitters and receivers for identifying minerals on the ocean floor |
US6188221B1 (en) * | 1998-08-07 | 2001-02-13 | Van De Kop Franz | Method and apparatus for transmitting electromagnetic waves and analyzing returns to locate underground fluid deposits |
GB9818875D0 (en) * | 1998-08-28 | 1998-10-21 | Norske Stats Oljeselskap | Method and apparatus for determining the nature of subterranean reservoirs |
WO2000013037A1 (en) * | 1998-08-31 | 2000-03-09 | Osaka Gas Co., Ltd. | Three-dimensional questing method, three-dimensional voxel data displaying method, and device therefor |
US6163155A (en) * | 1999-01-28 | 2000-12-19 | Dresser Industries, Inc. | Electromagnetic wave resistivity tool having a tilted antenna for determining the horizontal and vertical resistivities and relative dip angle in anisotropic earth formations |
US6184685B1 (en) * | 1999-02-22 | 2001-02-06 | Halliburton Energy Services, Inc. | Mulitiple spacing resistivity measurements with receiver arrays |
AU3622200A (en) * | 1999-03-12 | 2000-09-28 | Profile Technologies, Inc. | Dynamic electromagnetic methods for direct prospecting for oil |
US6294917B1 (en) * | 1999-09-13 | 2001-09-25 | Electromagnetic Instruments, Inc. | Electromagnetic induction method and apparatus for the measurement of the electrical resistivity of geologic formations surrounding boreholes cased with a conductive liner |
GB0002422D0 (en) * | 2000-02-02 | 2000-03-22 | Norske Stats Oljeselskap | Method and apparatus for determining the nature of subterranean reservoirs |
US7167526B2 (en) * | 2001-06-07 | 2007-01-23 | National Univ. Of Singapore | Wireless communication apparatus and method |
GB2383133A (en) * | 2001-08-07 | 2003-06-18 | Statoil Asa | Investigation of subterranean reservoirs |
GB2413188B (en) * | 2001-08-07 | 2006-01-11 | Electromagnetic Geoservices As | Method and apparatus for determining the nature of subterranean reservoirs |
GB2382875B (en) * | 2001-12-07 | 2004-03-03 | Univ Southampton | Electromagnetic surveying for hydrocarbon reservoirs |
US6842006B2 (en) * | 2002-06-27 | 2005-01-11 | Schlumberger Technology Corporation | Marine electromagnetic measurement system |
US6927728B2 (en) * | 2003-03-13 | 2005-08-09 | Motorola, Inc. | Method and apparatus for multi-antenna transmission |
US7385914B2 (en) * | 2003-10-08 | 2008-06-10 | Atheros Communications, Inc. | Apparatus and method of multiple antenna transmitter beamforming of high data rate wideband packetized wireless communication signals |
-
2001
- 2001-12-07 GB GB0129407A patent/GB2382875B/en not_active Expired - Lifetime
-
2002
- 2002-11-28 AU AU2002349125A patent/AU2002349125A1/en not_active Abandoned
- 2002-11-28 BR BRPI0214678A patent/BRPI0214678B1/en unknown
- 2002-11-28 BR BR0214678-9A patent/BR0214678A/en active IP Right Grant
- 2002-11-28 CA CA2468896A patent/CA2468896C/en not_active Expired - Fee Related
- 2002-11-28 WO PCT/GB2002/005355 patent/WO2003048812A1/en not_active Application Discontinuation
- 2002-11-28 US US10/497,807 patent/US7126338B2/en not_active Expired - Lifetime
-
2004
- 2004-07-06 NO NO20042852A patent/NO325116B1/en not_active IP Right Cessation
-
2006
- 2006-10-24 US US11/585,585 patent/US20070150201A1/en not_active Abandoned
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20210341638A1 (en) * | 2019-09-05 | 2021-11-04 | Institute Of Geology And Geophysics Chinese Academy Of Sciences (Iggcas) | Method for Extracting IP Information from Grounded-wire Transient Electromagnetic Response |
Also Published As
Publication number | Publication date |
---|---|
GB2382875A (en) | 2003-06-11 |
BR0214678A (en) | 2004-11-23 |
NO20042852L (en) | 2004-09-07 |
US20070150201A1 (en) | 2007-06-28 |
CA2468896A1 (en) | 2003-06-12 |
US7126338B2 (en) | 2006-10-24 |
BRPI0214678B1 (en) | 2020-04-14 |
NO325116B1 (en) | 2008-02-04 |
WO2003048812A1 (en) | 2003-06-12 |
GB0129407D0 (en) | 2002-01-30 |
GB2382875B (en) | 2004-03-03 |
AU2002349125A1 (en) | 2003-06-17 |
US20050077902A1 (en) | 2005-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2468896C (en) | Electromagnetic surveying for hydrocarbon reservoirs | |
US7679367B2 (en) | Electromagnetic surveying for hydrocarbon reservoirs | |
US7483792B2 (en) | Electromagnetic surveying for hydrocarbon reservoirs | |
CA2598468C (en) | Electromagnetic surveying for resistive or conductive bodies | |
AU2005215434B2 (en) | Electromagnetic surveying for resistive or conductive bodies | |
US7565245B2 (en) | Electromagnetic surveying |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request | ||
MKLA | Lapsed |
Effective date: 20201130 |