A method for removing the effects of an airwave from marine electromagnetic data comprising providing an electromagnetic source and at least one receiver in the water; measuring the electromagnetic response at a first source-receiver separation; determining a scaled version of the airwave response at a source-receiver separation where the earth response is negligible and using the scaled airwave response to determine the earth response measured at the first separation. Using this method, an improved estimate of the earth's response can be achieved.

Patent
   8427165
Priority
Feb 12 2007
Filed
Feb 11 2008
Issued
Apr 23 2013
Expiry
Jan 06 2030

TERM.DISCL.
Extension
695 days
Assg.orig
Entity
Large
1
12
EXPIRED
10. A system for determining the earth response to an electromagnetic signal comprising:
an electromagnetic source;
at least one receiver in water, the at least one receiver being adapted to:
measure a first electromagnetic response at a first source-receiver separation; and
measure a second electromagnetic response at a second source-receiver separation where the earth response is negligible; and
a processor, the processor being adapted to:
determine from the second electromagnetic response the airwave response at the second source-receiver separation; and
use the determined airwave response and the measured first electromagnetic response to determine the earth response at the first source-receiver separation by:
finding an inverse filter f(t) of the determined airwave response â(t);
convolving the filter f(t) with the measured response at the first source-receiver separation to provide a result that includes an airwave impulse;
removing the airwave impulse from the result of the convolution; and
convolving the remainder with the determined airwave response â(t) to recover the earth response at the first source-receiver separation without the airwave;
wherein the earth response comprises a signal with a transmission path directly through the earth.
1. A method for determining the earth response to an electromagnetic signal using a source for transmitting the electromagnetic signal to at least one receiver in water, the method comprising:
measuring an first electromagnetic response at a first source-receiver seperation;
measuring a second electromagnetic response at a larger second source-receiver seperation where the earth response is negligible;
determining from the second electromagnetic response the airwave response at the second source-receiver separation; and
using the determined airwave response and the measured first electromagnetic response to determine the earth response at the first source-receiver separation,
wherein using the determined airwave response and the measured first electromagnetic response to determine the earth response involves:
finding an inverse filter f(t) of the determined airwave response â(t);
convolving the filter f(t) with the measured response at the first source-receiver seperation to provide a result that includes an airwave impulse;
removing the airwave impulse from the result of the convolution; and
convolving the remainder with the determined airwave response â(t) to recover the earth response at the first source-receiver separation without the airwave; wherein the earth response comprises a signal with a transmission path directly through the earth.
2. A method as claimed in claim 1 wherein the step of determining the earth response involves measuring or estimating the airwave response at the second source-receiver separation.
3. A method as claimed in claim 1 wherein determining the airwave response involves estimating its shape.
4. A method as claimed in claim 1 wherein the second source-receiver separation is approximately five or more times greater than the first source-receiver separation.
5. A method as claimed in claim 1 wherein both the first and the second source-receiver separations are at least five times greater than the depth of the water.
6. A method as claimed in claim 1 wherein the second source-receiver separation is selected so that the earth response measured at the receiver is substantially zero.
7. A method as claimed in claim 1 wherein the source is a current bi-pole or a magnetic loop.
8. A method as claimed in claim 1 the receiver is a bi-pole receiver or a magnetic loop.
9. A method as claimed in claim 1 comprising locating the source and receiver on or on the seabed.
11. A system as claimed in claim 10 that is adapted to determine the earth response by measuring or estimating the airwave response at the second source-receiver separation.
12. A system as claimed in claim 10 wherein the second source-receiver separation is approximately five or more times greater than the first source-receiver separation.
13. A system as claimed in claim 10 wherein both the first and second source-receiver separations are at least five times greater than the depth of the water.
14. A system as claimed in claim 10 wherein the second source-receiver separation is selected so that the earth response measured at the receiver is substantially zero.
15. A system as claimed in claim 10 wherein the source is a current bi-pole or a magnetic loop.
16. A system as claimed in claim 10 wherein the receiver is a bi-pole receiver or a magnetic loop.
17. A system as claimed in claim 10 wherein the source and receiver are located on or near the seabed.

The present application is the U.S. national stage application of International Application PCT/GB2008/000467, filed on Feb. 11, 2008, which international application was published on Aug. 21, 2008 as International Publication WO 2008/099151. The International Application claims priority of United Kingdom Patent Application GB 0702661.0 filed on Feb. 12, 2007.

The present invention relates to a method for removing the airwave from marine electromagnetic data. In particular, the invention relates to a technique for reducing the impact of the airwave in marine based multi-channel transient electromagnetic (MTEM) measurements.

Porous rocks are saturated with fluids. The fluids may be water, gas, or oil, or a mixture of all three. The flow of current in the earth is determined by the resistivities of such rocks, which are affected by the saturating fluids. For instance, brine-saturated porous rocks are much less resistive than the same rocks filled with hydrocarbons. By measuring the resistivity of geological formations, hydrocarbons can be detected. Hence, resistivity measurements can be made in an exploration phase to detect hydrocarbons prior to drilling.

Various techniques for measuring the resistivity of geological formations are known, for example time domain electromagnetic techniques, as described in WO 03/023452, the contents of which are incorporated herein by reference. Conventionally, time domain electromagnetic investigations use a transmitter and one or more receivers. The transmitter may be an electric source, that is, a grounded bipole, or a magnetic source, that is, a current in a wire loop or multi-loop. The receivers may be grounded bipoles for measuring potential differences, or wire loops or multi-loops or magnetometers for measuring magnetic fields and/or the time derivatives of magnetic fields. The transmitted signal is often formed by a step change in current in either an electric or magnetic source, but any transient signal may be used, including, for example, a pseudo-random binary sequence. Measurements can be taken on land or in an underwater environment.

FIG. 1 shows a view of a typical setup for transient electromagnetic marine surveying. This has a bi-pole current source with mid-point xs on or near the sea floor for transmitting a transient current between two electrodes. The time function of the current might be a simple step change in current or a more complicated signal such as a pseudo-random binary sequence. The response of the earth-water system is measured by a line of bi-pole receivers on or near the sea floor, each receiver characterised by its mid-point position xr and measuring the potential difference between a pair of electrodes. All the electrodes are essentially in the same vertical plane.

In use, the electromagnetic signal generated by the source of FIG. 1 can follow three transmission paths to the receiver electrodes, these being directly through the earth, directly through the water, and via the water through the air and back through the water again. The signal transmitted by this third path is known as the airwave. In deep water, the airwave has a negligible impact. In contrast in shallow water, the signal that is transmitted through the water is negligible, but the airwave can have a significant impact and so make interpretation of the data difficult.

According to one aspect of the present invention, there is provided a method for removing the effects of an airwave from marine electromagnetic data comprising providing a source and at least one receiver in the water; measuring the electromagnetic response at a first source-receiver separation; determining a scaled version of the airwave response at a source-receiver separation where the earth response is negligible and using the scaled airwave response to determine the earth response measured at the first separation. The step of determining the scaled airwave may involve measuring or estimating the airwave response.

For an impulsive source current, the airwave response decays approximately as the cube of the source-receiver separation at separations large compared with the water depth, whereas the response of the component travelling directly through the earth decays approximately as the fifth power of the source-receiver separation. At very large offsets, therefore, the airwave dominates the total response. Hence, by measuring the electromagnetic response at a large source-receiver separation where the earth response is negligible, a measure of the airwave at that source receiver separation can be determined. This is a scaled version of the airwave response at the first source-receiver measurement separation, and can be used to allow removal of the airwave from the measured data and consequently a more accurate measure of the earth response at the first source-receiver separation.

Using the scaled airwave response to determine the earth response may involve finding an inverse filter f(t) of the scaled airwave response â(t); convolving the inverse filter f(t) with the measured data x(t) to compress the airwave to an impulse; removing the airwave impulse from the result of the convolution with the air wave data and convolving the result with the scaled determined airwave â(t) to recover the earth response without the airwave.

Preferably, the source is a current bi-pole or a magnetic loop. Preferably, the receiver is a bi-pole receiver or a magnetic loop. Preferably, the source and receiver are located on or near the seabed.

The electromagnetic measurements are transient electromagnetic measurements; that is, the response to the input signal at the source must have reached essentially steady state before the measurement is complete.

Preferably, the larger source-receiver separation is approximately five or more times greater than the shorter separation. Ideally, both source-receiver separations are at least five times greater than the depth of the water.

The larger source-receiver separation is selected so that the earth response measured at the receiver is substantially zero.

According to another aspect of the present invention, there is provided a computer program, preferably on a computer readable medium or data carrier, for removing the effects of an airwave from marine electromagnetic data as measured by at least one source and at least one receiver, the computer program having code or instructions for using an airwave response at a source-receiver separation where the earth response is negligible to allow removal of an airwave response in measurements at a shorter source-receiver separation, thereby to allow the earth response to be determined.

Preferably, the larger separation is approximately five or more times greater than the shorter separation.

Preferably, the larger separation is selected so that the earth response measured at the receiver is substantially zero.

Preferably, the program has code or instructions for subtracting the determined airwave from the measurement at the shorter separation.

According to yet another aspect of the invention, there is provided a system for determining the earth response to an electromagnetic signal using an electromagnetic source and at least one receiver in water, the system being adapted to measure the electromagnetic response at a first source-receiver separation; determine the airwave response or a scaled version thereof at a source-receiver separation where the earth response is negligible and use the determined airwave response to determine the earth response measured at the first separation.

The system may be adapted to use the airwave response to determine the earth response involves finding an inverse filter f(t) of the scaled airwave response â(t); convolve the filter f(t) with the measured response at the first source-receiver separation to provide a result that includes an airwave impulse; remove the airwave impulse from the result of the convolution and convolve the remainder with the determined airwave response â(t) to recover the earth response at the first source-receiver separation without the airwave.

A system may be adapted to determine the earth response by measuring or estimating the airwave response.

The larger separation may be approximately five or more times greater than the shorter separation. Both separations may be at least five times greater than the depth of the water.

The larger separation may be selected so that the earth response measured at the receiver is substantially zero.

The source may be a current bi-pole or a magnetic loop. The receiver may be a bi-pole receiver or a magnetic loop.

The source and receiver may be located on or near the seabed.

According to still another aspect of the present invention, there is provided a method of measuring an electromagnetic response, the method comprising positioning an electromagnetic source and at least one electromagnetic receiver underwater; transmitting an electromagnetic signal from the source to the receiver, wherein the transmitted signal comprises an earth component and an airwave component, and the receiver is positioned relative to the source so that the measured response is substantially due to the airwave.

In any aspect of the invention, the transmitted signal may be formed by a step change in current in either an electric or magnetic source, although any transient signal may be used, including, for example, a pseudo-random binary sequence.

Various aspects of the invention will now be described by way of example only and with reference to the accompanying drawings, of which:

FIG. 2 is a model of the system of FIG. 1;

FIG. 3 shows an impulse response at 2 km for the arrangement of FIG. 2 in which the earth resistivity is 1 ohm-m; FIG. 3(a) shows the response out to 9 s, and FIG. 3(b) shows it out to 1 s;

FIG. 4 shows the impulse response at 10 km offset for a 1 ohm-m half-space beneath 100 m of water, where the dipole source and receiver are on the sea floor;

FIG. 4 (a) shows the response out to 9 s, and FIG. 4 (b) shows it out to 1 s;

FIG. 5 shows the marine impulse response, the earth impulse response and the airwave for a 2 km offset measurement, and

FIGS. 6(a) and (b) show the result of removing the airwave from the impulse responses of FIG. 3.

The present invention provides a method for removing the effects of an airwave from shallow water marine electromagnetic data. The invention will be illustrated using the model of FIG. 2, which shows a seawater layer of depth d and resistivity ρw overlying a homogeneous isotropic halfspace. The dipole source is on the sea floor and the in-line electric field is calculated at a receiver position at some offset r from it on the sea floor. For the sake of example, it will be assumed that the water layer is 100 m deep and has a resistivity of 0.3 ohm-m, and the half space has a resistivity of 1 ohm-m.

FIG. 3 shows an impulse response at 2 km for the arrangement of FIG. 2. FIG. 3(a) shows the response out to 9 s, and FIG. 3(b) shows it out to 1 s. A sharp initial peak occurs at 0.008 s and has an amplitude of 1.82E-10 V/m/s. This is the airwave. This has a finite rise time to the peak and a long tail. There is a second much broader peak at about 0.45 s with an amplitude of 2.4E-11 V/m/s. This is caused by the earth response. The objective is to separate the airwave and the earth impulse response.

At a given offset and after deconvolution for the system response, the impulse response data, which are causal, can be expressed as:
x(t)=a(t)+g(t), t≧0,
x(t)=0, t<0.  (1)
in which a(t) is the air wave and g(t) is the earth impulse response, with g(0)=0. Suitable deconvolution techniques are described for example by Wright, D., et al in the article Hydrocarbon detection and monitoring with a multichannel transient electromagnetic (MTEM) survey: The Leading Edge, 21, 852-864.

The airwave is expected to be minimum-phase. By this it is meant that without the water the airwave is an impulse and the propagation effects in the water must be factors involving decaying exponentials, which are minimum-phase.

An estimate â(t) of the air wave is assumed such that
(t)≈a(t).  (2).

This is an estimate of the shape, but the exact airwave amplitude A is not known. Next a filter f(t) has to be found, such that
f(t)*â(t)=δ(t),  (3)
in which the asterisk denotes convolution and δ(t) is the impulse function:
δ(t)=0, t≠0
−∞δ(t)dt=1.  (4)

As mentioned above a(t) and hence â(t) are expected to be minimum-phase. Therefore the inverse f(t) is causal, realizable and minimum-phase; that is, it exists only for t≧0, and it has finite amplitude.

Convolving the filter f(t) with the data yields
f(t)*x(t)≈f(t)*(t)+f(t)*g(t)=(t)+f(t)*g(t).  (5)

Now the airwave is removed by considering only t>0, so that:
f(t)*x(t)=f(t)*g(t), t>0,  (6)

This result can then be convolved with â(t):
â(t)*f(t)*g(t)=δ(t)*g(t)=g(t).  (7)

This allows the earth impulse as measured to be recovered, but without the airwave. This is what is needed, and it turns out that the unknown amplitude factor A is irrelevant. Hence, by merely having an estimate or measure of the shape of the airwave, the earth impulse can be determined.

Because of the presence of noise, an exact inverse of â(t) cannot be used. First equation (1) is re-written to include the noise:
x(t)=a(t)+g(t)+n(t)  (8)

Then equation (2) is modified to define a new filter f1(t), such that
f1(t)*â(t)=d(t−τ),  (9)
in which d(t) is a band-limited impulse whose integral is unity, for example a Gaussian, of the form d(t)=αexp(−πα2t2) with α a frequency close to the highest frequency of interest in the data, and τ is a time delay introduced to make f1(t) causal. The integral of d(t) is 1. Convolving this filter with the data yields
f1(t)*x(t)=f1(t)*(t)+f1(t)*g(t)+f1(t)*n(t)=Ad(t−τ)+f1(t)*g(t)+f1(t)*n(t)=y(t).  (10)

The value of α is chosen to be as large as possible but not so large that the high frequency components of the noise are amplified too much. Then
d(t)≈0 for values of |t|>3/α.  (11)

The function y(t) is set to zero for |t−τ|≦3/α to yield z(t). That is,
z(t)≈f1(t)*g(t)+f1(t)*n(t), t−τ>3/α.  (12)

The approximate inverse of f1(t) is â(t), as defined in equation (9).

Then â(t) is convolved with z(t) to yield
â(t)*z(t)≈â(t)*f1(t)*g(t)+â(t)*f1(t)*n(t)≈d(t−τ)*g(t)+d(t−τ)*n(t)≈g(t−τ)+n(t−τ).  (13)

Apart from the known time delay τ this is what is desired, plus noise. The time delay can easily be removed.

The shape of the airwave can be estimated by looking at an offset that is large compared with this offset and where the attenuation of the earth response is relatively high. It is suspected that the shape of the airwave is invariant with offset for a 1-D earth—beyond a certain offset—and the hypothesis is that an almost pure airwave could be measured, if the offset is large enough and the time window restricted. This gives an estimate of the shape of the airwave, which according to the method described above would allow the effects of the airwave at the measurement source-receiver separations to be removed.

FIG. 4 shows the impulse response at 10 km offset for a 1 ohm-m half space beneath 100 m of water, where the dipole source and receiver are on the sea floor. FIG. 4 (a) shows the response out to 9 s, and FIG. 4 (b) shows it out to 1 s. The sharp initial peak occurs at 0.008 s and has an amplitude of 1.46E-12 V/m/s. There is no discernible second peak. It should be noted that the peak of the airwave occurs at exactly the same time as for the 2 km impulse response.

FIG. 5 shows the marine impulse response, the earth impulse response and the airwave for a 2 km offset measurement. The earth impulse response was obtained according to the scheme described above; the airwave was obtained by subtracting the earth impulse response from the marine impulse response. A measure of the quality of this result is to compare the recovered impulse response with the analytic function: the impulse response of a 1 ohm-m half space at an offset of 2 km. FIG. 6 shows the result of removing the airwave from the impulse responses shown in FIG. 3. This is very similar to the analytic function for the impulse response of a 1 ohm-m half space. In particular the time and amplitude of the peak of the two functions are in close agreement.

The present invention provides a very simple and effective technique for removing the airwave from measured MTEM data. This can be summarised as follows: estimate the shape of the airwave using an airwave at a large source-receiver separation and find a filter f(t) based on the estimated airwave. This filter is then convolved with the measured data x(t). The result of this is very messy data, but there is a very sharp peak at time t=τ. This is the airwave. The airwave peak is then sliced out. The remaining data is returned to its original state by convolving it with the scaled airwave a(t). This operation undoes the effect of convolving with the filter f(t). The result is the original data, i.e. the earth response, without the airwave.

The data captured may be processed in accordance with the invention using any suitable hardware or software.

A skilled person will appreciate that variations of the disclosed arrangements are possible without departing from the invention. Accordingly the above description of the specific embodiment is made by way of example only and not for the purposes of limitation. It will be clear to the skilled person that minor modifications may be made without significant changes to the operation described.

Ziolkowski, Antoni Marjan

Patent Priority Assignee Title
10878142, Oct 26 2015 PGS Geophysical AS Bipole source modeling
Patent Priority Assignee Title
7982465, Feb 13 2006 MULTIFIELD GEOPHYSICS AS Electromagnetic method on shallow water using a controlled source
8099239, Feb 16 2004 PGS Geophysical AS Electromagnetic surveying for resistive or conductive bodies
8264227, Aug 25 2006 PGS EM Limited Marine EM exploration
20040239330,
20070061078,
20080150538,
20090278541,
20120011130,
WO3023452,
WO3048812,
WO2005010560,
WO2007094676,
///
Executed onAssignorAssigneeConveyanceFrameReelDoc
Feb 11 2008PGS EM Limited(assignment on the face of the patent)
Jan 08 2010ZIOLKOWSKL, ANTONI MARJANMTEM, LT5ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0237970944 pdf
Aug 21 2012MTEM LIMITEDPGS EM LimitedCHANGE OF NAME SEE DOCUMENT FOR DETAILS 0291030273 pdf
Date Maintenance Fee Events
Oct 24 2016M1551: Payment of Maintenance Fee, 4th Year, Large Entity.
Dec 14 2020REM: Maintenance Fee Reminder Mailed.
May 31 2021EXP: Patent Expired for Failure to Pay Maintenance Fees.


Date Maintenance Schedule
Apr 23 20164 years fee payment window open
Oct 23 20166 months grace period start (w surcharge)
Apr 23 2017patent expiry (for year 4)
Apr 23 20192 years to revive unintentionally abandoned end. (for year 4)
Apr 23 20208 years fee payment window open
Oct 23 20206 months grace period start (w surcharge)
Apr 23 2021patent expiry (for year 8)
Apr 23 20232 years to revive unintentionally abandoned end. (for year 8)
Apr 23 202412 years fee payment window open
Oct 23 20246 months grace period start (w surcharge)
Apr 23 2025patent expiry (for year 12)
Apr 23 20272 years to revive unintentionally abandoned end. (for year 12)