Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
GEOPHYSICAL DATA PROCESSING SYSTEMS
Document Type and Number:
WIPO Patent Application WO/2011/098821
Kind Code:
A2
Abstract:
A system for measuring geophysical data from a potential field survey of a surveyed region of the earth to provide a representation of the geology of said surveyed region, the system comprising: a gradiometer for measuring gravity gradient data; and a gravimeter for measuring gravity data wherein the gradiometer and gravimeter are mounted on a common stabilised platform whereby the gradiometer and gravimeter are fixed relative to each other.

Inventors:
HOUGHTON PHILIP (GB)
LUMLEY JOHN MORRIS (GB)
DYER NEIL (GB)
Application Number:
PCT/GB2011/050251
Publication Date:
August 18, 2011
Filing Date:
February 11, 2011
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
ARKEX LTD (GB)
HOUGHTON PHILIP (GB)
LUMLEY JOHN MORRIS (GB)
DYER NEIL (GB)
International Classes:
G01V7/16
Domestic Patent References:
WO2007012895A22007-02-01
WO2007085875A22007-08-02
WO2008093139A22008-08-07
WO2009016348A22009-02-05
WO2008117081A22008-10-02
Foreign References:
US20090064778A12009-03-12
US4399694A1983-08-23
US3731537A1973-05-08
US5668315A1997-09-16
GB2008050041W2008-01-21
Other References:
LEE ET AL.: "Effects of data filtering on inversion of gravity gradient data", SEG 75TH ANNUAL EXPOSITION AND ANNUAL MEETING, 6 November 2005 (2005-11-06)
R.J. BLAKELY: "Potential Theory in Gravity and Magnetic Applications", 1995, CAMBRIDGE UNIVERSITY PRESS, pages: 185
JACOBY W.; SMILDE P.L: "Gravity Interpretation Fundamentals and Application of Gravity Inversion and Geological Interpretation", 2009, SPRINGER
"Numerical Recipes in C", 1997, PRESS W.H.
LUMLEY, A SUPERCONDUCTING GRAVITY GRADIOMETER TOOL FOR EXPLORATION, 2004
Attorney, Agent or Firm:
MARKS & CLERK LLP (Cambridge, Cambridgeshire CB2 1LA, GB)
Download PDF:
Claims:
CLAIMS:

1. A system for measuring geophysical data from a potential field survey of a surveyed region of the earth to provide a representation of the geology of said surveyed region, the system comprising:

a gradiometer for measuring gravity gradient data; and

a gravimeter for measuring gravitydata

wherein the gradiometer and gravimeter are mounted on a common stabilised platform whereby the gradiometer and gravimeter are fixed relative to each other.

2. A system according to claim 1, wherein the common stabilised platform forms part of a first gimbal assembly whereby the gradiometer and gravimeter are stabilised with respect to a first axis.

3. A system according to claim 2, comprising a second gimbal assembly for stabilisation with respect to a second axis which is orthogonal to the first axis and wherein the common stabilised platform is connected to the second gimbal assembly.

4. A system according to claim 3, comprising a third gimbal assembly for stabilisation with respect to a third axis which is orthogonal to both the first and second axes and wherein the common stabilised platform is connected to the third gimbal assembly.

5. A system according to any one of the preceding claims comprising a data acquisition module for acquiring data from the gradiometer and gravimeter and a GPS instrument for simultaneously obtaining the time and location for the data acquired from the gradiometer and gravimeter.

6. A system according to any one of the preceding comprising

means for filtering said measured gravity data and

a processor which is configured to:

input said measured gravity and gravity gradient data; and

generate a model of said surveyed region by predicting gravity data and gravity gradient data corresponding to said measured gravity and gravity gradient data;

filtering said predicted gravity data, wherein said filtering matches the filtering applied to the measured gravity data and

fitting said filtered predicted gravity data and said predicted gravity gradient data to said measured gravity data and gravity gradient data, respectively.

7. A system according to claim 6, comprising means for filtering said measure gravity gradient data whereby the measured gravity gradient data is filtered data and the generating step includes filtering said predicted gravity gradient data and fitting said filtered predicted data to said measured data.

8. A system according to claim 6 or claim 7, comprising a low pass filter for filtering the predicted data.

9. A system according to claim 8, wherein the low pass filter has a cutoff frequency of approximately 10-20 mHz.

10. A system according to any one of claims 6 to 9, wherein the processor is configured to filter the predicted data in the time domain.

11. A system according to claim 10, wherein the processor is configured to filter by applying a weight to each element in the predicted data set.

12. A system according to claim 11, wherein the processor is configured to calculate the weight for each element based on a time characteristic of the means for filtering and a time difference between measurement of a pair of elements.

13. A system according to any one of claims 1 to 5, comprising

means for filtering said measured gravity data and

a processor which is configured to:

input said measured gravity and gravity gradient data; and generate a model of said surveyed region by

predicting gravity data and gravity gradient data corresponding to said measured gravity and gravity gradient data;

weighting each element of said predicted data whereby predicted data corresponding to data which is filtered out of said measured gravity data has a lower weight than both predicted data which is not filtered out of said measured gravity data and said predicted gravity gradient data; and

fitting said weighted predicted data to said measured data to generate the model.

14. A system according to claim 13, wherein said weighting step generates the effect of a low pass filter on the predicted gravity data and leaves the predicted gravity gradient data unchanged.

15. A system according to any one of claims 6 to 14, wherein the processor is configured to generate the model using regularisation.

Description:
Geophysical Data Processing Systems FIELD OF THE INVENTION

This invention relates to methods, apparatus, and computer program code for processing geophysical data, more particularly gravity data and gravity gradiometer data from a potential field survey to provide a representation of the underlying geology of the surveyed region. Embodiments of the techniques we describe are particularly useful for processing data from airborne surveys, in particular gravity field surveys.

BACKGROUND TO THE INVENTION A potential field survey is performed by measuring potential field data which, for a gravity survey, may comprise one or more of gravimeter data (measuring gravity field) or gravity gradiometer data (measuring gravity field gradient), vector magnetometer data, true magnetic gradiometer data, and other types of data well-known to those skilled in the art. A common aim of a geophysical potential field survey is to search for signatures which potentially indicate valuable deposits such as minerals or hydrocarbons.

The goal of most potential field surveys is to deduce the sub-surface geology. It is well known to the geophysics industry that the high spatial frequency potential field signals are dominated by near surface anomalies. Further, in most situations, the terrain and the near surface geology also give rise to substantial low spatial frequency signals. In general, the low frequency signals from the near surface anomalies are extremely hard to distinguish from the low frequency signals coming from deeper geology. Accordingly, it is possible to postulate a geology having only a near surface mass distribution which matches the measured signal implying that there are no anomalies left in the measured signal coming from deep geologies. Fortunately most such mass distributions can be ruled out using a range of tools including geological knowledge, known typical density distributions, prior information, similar known structures elsewhere in the world etc.

It is known that gravity data is affected by high-frequency noise originating from movement of a platform on which the gravimeter is placed. Low-pass filters are often applied to remove such motion-related noise. However, the filtering does not discriminate between signal and noise and removes both from frequencies outside the pass band. If this is not taken into account, interpretation of the data from such filtered gravimeters produces source bodies which are deeper and wider than the true sources. This problem is discussed in "Effects of data filtering on inversion of gravity gradient data" by Lee et al, SEG 75 th Annual Exposition and Annual Meeting 6-11 November 2005, Houston Texas. The paper suggests incorporating filtering into the inversion process to improve results. US2009/0064778, US4399694, US3731537 and US5668315 describe various gravity gradiometers.

There remains a need, however, for improved techniques for processing geophysical data from such surveys in order to identify the underlying geology.

SUMMARY OF THE INVENTION

According to a first aspect of the invention there is provided a system for measuring geophysical data from a potential field survey of a surveyed region of the earth to provide a representation of the geology of said surveyed region, the system comprising: a gradiometer for measuring gravity gradient data; and

a gravimeter for measuring gravity data

wherein the gradiometer and gravimeter are mounted on a common stabilised platform whereby the gradiometer and gravimeter are fixed relative to each other.

In this way, the system uses two different instruments in tandem. The first instrument, the gravimeter delivers a better measure of the long wavelength component of the Earth's gravity. The second instrument, the gradiometer delivers a better measure of the short wavelength component of the Earth's gravity. For gradient data it is necessary for the accelerations to be known to a few parts in 10,000 but for the gravity data, the sensitivity is greater and needs to be to a few parts in 100,000. Accordingly, it is not possible to use the output from a gradiometer to measure gravity data with sufficient accuracy.

The common stabilised platform may form part of a first gimbal assembly whereby the gradiometer and gravimeter are stabilised with respect to a first axis. The system may further comprise a second gimbal assembly for stabilisation with respect to a second axis which is orthogonal to the first axis and wherein the common stabilised platform is connected to the second gimbal assembly. The system may further comprise a third gimbal assembly for stabilisation with respect to a third axis which is orthogonal to both the first and second axes and wherein the common stabilised platform is connected to the third gimbal assembly.

The system may comprise a data acquisition module for acquiring data from the gradiometer and gravimeter and a GPS instrument for simultaneously obtaining the time and location for the data acquired from the gradiometer and gravimeter.

The system may comprise means for filtering said measured gravity data and a processor for processing the data, the processor being configured to::

input said measured gravity and gravity gradient data; and

generate a model of said surveyed region by

predicting gravity data and gravity gradient data corresponding to said measured gravity and gravity gradient data;

filtering said predicted gravity data, wherein said filtering matches the filtering applied to the measured gravity data and

fitting said filtered predicted gravity data and said predicted gravity gradient data to said measured gravity data and gravity gradient data, respectively. Both sets of data may be input to the system from the same measuring instrument or may be measured separately and input for subsequent processing together. The second set of measured potential field data may also be filtered data and the generating step may include filtering said second set of predicted data and fitting said second set of filtered predicted data to said second set of filtered measured data. If the second set of data is filtered, the filtering is preferably different to that applied to the first set of data. The method may further comprise input additional sets of measured data. In this case, the model may be generated by filtering the predicted data if the corresponding measured data is filtered data and fitting the additional sets of predicted data (whether filtered or unfiltered) to the corresponding additional sets of measured data.

Incorporating the filtering into the model generation process allows the solution to correctly focus on the relevant parts of each of the data sets' spectra. In other words, the different sets of data are combined in an optimal way. Thus, information that may have been lost in one set of data (e.g. due to noise levels at certain frequencies) is made up by another set whose measurements are unfiltered or are filtered over a different bandwidth.

The second set of measured potential field data may be gravity gradient data and/or said first set of filtered measured potential field data may be gravity data. As explained above, typically gravity data is heavily filtered. By contrast, gravity gradient data is essentially unfiltered. By using the above method, gravity and gravity gradient data may be combined in a way that correctly accounts for the filtering applied to the measured data. As a result, better results over a wider bandwidth should be obtainable. Furthermore, the above method allows the gravity and gravity gradient measurements to be made simultaneously and subsequently processed together.

The method may comprise filtering the predicted data (i.e. the first set of data alone or both sets of data if the second set of measured data is also filtered) in the time domain. In this way, noisy parts of the spectrum may be removed to minimise corruption of the model generation process. For measured data which is noisy at high frequencies, e.g. airborne gravity data, the filtering may use a low pass filter, e.g. with a cut-off frequency of 3.5mHz. Alternatively, different types of filter may be used to match the filtering applied to the measured data.

The predicted data may be filtered by applying a weight to each element in the data set.

Alternatively, the system may comprise means for filtering said measured gravity data and a processor for processing the data, the processor being configured to:

input said measured gravity and gravity gradient data; and

generate a model of said surveyed region by

predicting gravity data and gravity gradient data corresponding to said measured gravity and gravity gradient data;

weighting each element of said predicted data whereby predicted data corresponding to data which is filtered out of said measured gravity data has a lower weight than both predicted data which is not filtered out of said measured gravity data and said predicted gravity gradient data; and

fitting said weighted predicted data to said measured data to generate the model.

The weighting step may generate the effect of a low pass filter on the first set of predicted data and may leave the second set of predicted data unchanged, e.g. if said second set of measured potential field data is gravity gradient data and if said first set of filtered measured potential field data is gravity data. Alternatively, the weighting may generate the effect of a filter, preferably different types of filter, on both sets of predicted data.

The weight for each element may be calculated based on a time characteristic of the filter and a time difference between measurement of a pair of elements in the data set. In other words, each measurement may be weighted by frequency. The model generating step may further comprise regularisation which may be of any standard form, e.g. Tikhonov and gradient regularisation. The model generating step may comprise any known fitting step, e.g. equivalent source minimisation, in which there is minimisation of the error between the measured data and the predicted data generated by the equivalent source model. The minimisation may use standard optimisation algorithms such as the conjugate gradient method or Monte Carlo methods.

In embodiments, in particular where the field comprises a gravity field, the analysis to generate a model comprises identifying spatial features by identifying one or more of maxima, minima and points/lines of inflection in the filtered or residual potential field data. Thus, for example, with the gravity gradient tensor the on-diagonal components G xx , G yy , and G zz , (which are differential signals) are interpreted to determine inflection points or changes in slope since these generally correspond to geologically significant features of the surveyed region such as an interface between two different

types/densities of rock. For off-diagonal components, in particular G zx and G zy (which emphasise symmetries in the x- and y- direction respectively, maxima and/or minima are preferably identified; for G xy points are preferably identified by locating pairs of dipoles since these tend to identify corners of a subterranean body. Off-diagonal elements G Z i (where i is x or y) tend to emphasise symmetries in the i-direction.

Optionally the co-ordinates system may be rotated about one or more axis to potentially identify further geologically useful information. In embodiments the co-ordinate axis may be selected to maximise the apparent useful geological information. Similarly magnetic data may be processed to identify maximum/minimum inflection points/lines.

Referring again to the above-mentioned maxima, minima and lines of inflection, multiple sets of spatial features may be identified for the filtered or residual potential field data, for example from different vector or tensor components of a surveyed gravity field, magnetic field or gravity gradient field and/or from survey data other than from a potential field survey. Other survey data which may be combined with the data obtained from a potential field survey includes (but is not limited to): topographic information, for example determined by lidar, spectral or more preferably hyperspectral imagery, gas saturation data, chemical analysis data (from soil sampling), and other soil survey data. The potential field survey may be conducted from a moving platform such as an aircraft.

The invention also provides a method of prospecting for oil or a mineral from the earth, the method including conducting a potential field survey according to an aspect or embodiment of the invention as described above to generate a representation of the underlying geology of the surveyed region, and then using this representation to locate the desired oil or mineral. A further aspect of the invention also provides oil or mineral extracted using this technique.

The invention further provides processor control code to implement the above-described methods, in particular on a data carrier such as a disk, CD- or DVD-ROM, programmed memory such as read-only memory (Firmware), or on a data carrier such as an optical or electrical signal carrier. Code (and/or data) to implement embodiments of the invention may comprise source, object or executable code in a conventional programming language (interpreted or compiled) such as C, or assembly code, code for setting up or controlling an ASIC (Application Specific Integrated Circuit) or FPGA (Field Programmable Gate Array), or code for a hardware description language such as Verilog (Trade Mark) or VHDL (Very high speed integrated circuit Hardware Description Language). As the skilled person will appreciate such code and/or data may be distributed between a plurality of coupled components in communication with one another, for example distributed across a network.

The invention further provides a data processing system configured to implement embodiments of the above-described methods, to determine one or more parameters relating to physical properties of the Earth's interior from processed geophysical data. Such a data processing system may comprise: data memory for storing measured potential field data and plot data for representing the underlying geology of the surveyed region, program memory storing processor control code as described above; and a processor coupled to said data memory and to said program memory to load and implement said control code.

BRIEF DESCRIPTION OF THE DRAWINGS These and other aspects of the invention will now be further described, by way of example only, with reference to the accompanying figures in which: Figure 1 shows an aircraft with flight survey data, and an example of a data processing system configured to implement an embodiment of a method according to the invention;

Figure 2a shows the signal to noise ratio (SNR) along a survey line (60m/s) flying over a wide bandwidth geological model for both gravity gradient Gzz and gravity gz;

Figure 2b shows the power spectral density (PSD) impulse response (from a point source) as a function of spatial frequency for both gravity gradient Gzz and gravity gz;

Figure 3 is a flowchart of the methodology used to process measured data;

Figure 4a is an illustration of a synthetic model having dimensions 150 x 150 km used to test the methodology of Figure 3;

Figure 4b shows the gravity measurements to be expected from the model of Figure 4a;

Figure 5 a shows the error distribution in the predicted gravity when calculated from an equivalent source model using only gravity gradient simulated measurements;

Figure 5b shows the error distribution in the predicted gravity when calculated using the methodology of Figure 3;

Figure 6a is a cross-section construction of a gradiometer along one axis;

Figure 6b is a circuit diagram for the gradiometer of Figure 6a, and

Figure 7 is a schematic showing the mounting of a gradiometer and gravity instrument on a stabilised platform in accordance with an aspect of the present invention. DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

Figure 1 shows an example of an aircraft 10 for conducting a potential field survey to obtain data for processing in accordance with a method as described below. The aircraft 10 comprises an inertial platform 12 on which is mounted a gravity gradiometer 14 and a gravimeter 15 which both provide potential field survey data to a data collection system 16. Alternatively, the functionality of the gravimeter and gravity gradiometer may be incorporated into a single measuring instrument which is configured to measure both gravity and gravity gradiometry.

The inertial platform 14 is fitted with an inertial measurement unit (IMU) 18 which also provides data to data collection system 16 typically comprising attitude data (for example, pitch, roll and yaw data), angular rate and angular acceleration data, and aircraft acceleration data. The aircraft is also equipped with a differential GPS system 20 and a LIDAR system 22 or similar to provide data on the height of the aircraft and of the underlying terrain. The aircraft 10 may also be equipped with other instrumentation 24 such as a magnetometer, TDEM system and/or hyperspectral imaging system, again feeding into the data collection system. The data collection system 16 also has an input from general aircraft instrumentation 26 which may comprise, for example, an altimeter, air and/or ground speed data and the like. The data collection system 16 may provide some initial data pre-processing, for example to correct the LIDAR data for aircraft motion and/or to combine data from the IMU 18 and DGPS 20. The data collection system 16 may be provided with a communications link 16a and/or non- volatile storage 16b to enable the collected potential field and position data to be stored for later processing. A network interface (not shown) may also be provided.

Data processing to generate map data for the potential field survey is generally (but not necessarily) carried out offline, sometimes in a different country to that where the survey data was collected. As illustrated a data processing system 50 comprises a processor 52 coupled to code and data memory 54, an input/output system 56 (for example comprising interfaces for a network and/or storage media and/or other communications), and to a user interface 58 for example comprising a keyboard and/or mouse. The code and/or data stored in memory 54 may be provided on a removable storage medium 60. In operation the data includes data collected from the potential field survey and the code comprises code to process this data to generate map data, for example in accordance with the procedure shown in Figure 3, described below.

Consider an airborne potential field survey such as a gravity survey, flown on a grid pattern defined by orthogonal sets of parallel lines (flight paths) on a two-dimensional surface which is draped over the underlying terrain. When looking for underlying anomalies the nearby mass has a dominating effect and to provide an accurate representation of deep features a good representation of surface features is desirable so as to be able to perform terrain correction by subtracting-off particularly the higher frequencies (which dominate the power spectrum). A signal with wavelength λ falls off with height z as exp(-kz) where k = 2π/λ (from which it can be seen that longer wavelengths are less attenuated) and the wavelength scale corresponds to a signature expected given a target' s size and depth.

For gravity, the relevant potential is the gravity scalar potential, Φ(Γ) , defined as where r , p{r'), G are respectively, the position of measurement of the gravity field, the mass density at location r' , and the gravitational constant. The gravitational acceleration, which is how a gravitational field is experienced, is the spatial derivative of the scalar potential. Gravity is a vector in that it has directionality. It is represented by three components with respect to any chosen Cartesian coordinate system as:

Each of these three components varies in each of the three directions and the nine quantities so generated form the Gravity gradient tensor:

Although there appear to be nine components of the gravity gradient tensor, there are only five independent components. Firstly, the tensor is symmetric as the order of differentiation of a scalar quantity does not matter (implying that G xy = G yx ). Secondly, outside of the source, the sum of the diagonal terms equals zero (Laplace's equation). The ability to measure several spatially independent gravity components has obvious advantages over conventional gravity measurements, which only recovers the vertical component (g z )

There is a relationship between the depth (and shape) of a buried object and the wavelength (and amplitude) of the detected signal. In general, a measured quantity - say a component of the gravity vector or of the gravity gradient tensor will be a summation of the form shown below. Here we use gg as notation for the measured quantity, for example G zz . = ∑ m mass— element measure mass— element )

all— masses

In the above equation F is called a Greens function (see for example, R.J. Blakely, "Potential Theory in Gravity and Magnetic Applications", Cambridge University Press, 1995, at page 185, incorporated by reference) and r^ s - e i e m e n t defines the location of the mass element (for example the centre of mass or some other defined point).

The functions F are standard functions, essentially, the influence a source (mass element) of unity mass or density and defined shape would have at the relevant (measurement) point. The source may be a point source, sphere or ellipsoid but, in practice is more often a prism, which may be irregular. For example, if the presence of a particular geological layer or, say, geological anomaly, e.g. a kimberlite pipe, is suspected a shape can be defined to take account of this. A number of textbooks list Greens functions for simple shapes; functions for more complex source geometries can be found in the literature. Also the source influence superposes so that if a complex shape can be discretised into a plurality of simpler shapes then the Greens functions for the discrete shapes can be added together. This in principle allows numerical values for the Greens function of any arbitrary shape to be determined, although in practice relatively simple shapes are generally preferable. By way of example, the Green's function F for a rectangular prism (Blakely, ibid, at page 187), has 8 terms each of which corresponds to a vertex of the prism.

Airborne gravity gradiometry data is known to be of much higher resolution than conventional gravity data. This arises firstly because the power in the gradient signal is concentrated at higher spatial frequencies and secondly since the gradiometer, being a differencing instrument, is less sensitive to aircraft motion. The latter negates the need for heavy filtering of the data to suppress noise. Airborne gravimeters, intrinsically sensitive to aircraft motion, depend on filtering to reduce the noise to an acceptable level. Ultimately, the resolution of processed airborne gravity data is limited by the accuracy in which GPS derived acceleration corrections can be applied.

Figure 2a shows the signal to noise ratio (SNR) of actual survey measurements for both gravity data gz and gravity gradiometry data Gzz. The gravity gradiometry data Gzz is deduced from Full Tensor Gradiometer (FTG) measurements resulting in a white noise level of 8 E / VHz. The gravity data gz is deduced from airborne gravity measurements limited by GPS derived acceleration errors. The greyed out area highlights where the SNR is less than one and the signal is essentially lost in the noise. At higher

frequencies (i.e. above lOmHz) and thus at shorter wavelengths the gravimeter SNR deteriorates rapidly rendering the signals in this data with wavelengths less than 4 km unrecoverable. By contrast, the SNR for the gravity gradiometry does not deteriorate rapidly until after lOOmHz. Current gradiometer technology in commercial use (the FTG), measuring gradients with increased power in the high frequencies, allows the recovery of much higher resolution data down to 350 m. Nevertheless as shown in Figure 2b gravity data is also useful. Figure 2b shows that, for point source buried at a depth of 1km, the signal power spectral density (PSD) continues to increase with wavelength for gravity data whereas the PSD dramatically drops off for gravity gradient data from a spatial frequency below 10 - " 3 m - " 1.

The example data of Figure 2a shows that the Gzz SNR appears to go flat at low frequency for this geological signal. This maintenance of power is a result of the geological distribution of mass which tends to be correlated with a reddish spectrum therefore boosting the low frequency signal. Nevertheless, for frequencies less than 3.5 mHz (equivalent to roughly 18 km), the gravity measurements start to offer greater SNR than the gravity gradient measurements.

To make optimum use of both sets of measurements, one should combine the long wavelength signals in the gravity measurements with the high bandwidth signals in the gravity gradient data. This will benefit regional style surveys ( > 100 km in extent) most noticeably where there is more bandwidth over which the gravity data is superior to the gravity gradient. For smaller surveys, the SNR within the gravity gradient measurements alone is sufficient to accurately image the entire bandwidth. Referring now to Figure 3 this shows an example of a procedure for implementing on a data processor which may, in embodiments, comprise a general purpose computer system, for processing data from a flight survey in accordance with the previously described techniques. Thus, at step S200 the procedure inputs the measured potential field data, e.g. from the gravity gradiometer and the gravimeter or from a combined instrument, together with other data, e.g. associated 3D position data. Optionally at step S200a, some pre-processing may be applied, for example to remove outliers and/or to select ranges of the data to be processed.

Because of the high level of acceleration noise that impinges on airborne gravimeters, the instruments tend to have an intrinsically low bandwidth to prevent saturation.

Further, due to the poor accuracy in which GPS derived acceleration corrections can be applied at high frequencies, the recorded data tends to be heavily time domain filtered ( > 50 seconds). When combining gravity measurements with essentially unfiltered gravity gradient measurements, proper account needs to be taken of this filtering otherwise the gravity measurements will downgrade the high-resolution information in the gradients. At step S202, the data is combined using an equivalent source inversion that is augmented with information regarding the measurement filtering to generate an equivalent source model. In this scheme, the standard equivalent source minimisation is augmented with a function that filters the predictions from the equivalent source model in the time domain to match the filtering and bandwidth of the measurements,

Minimise { II Filter[Predictions(p)] - Measurements ll 2 + regularisation(p) } (1)

where Filter is a time domain filter matching the filtering and bandwidth of the measurements, p is the density distribution of the equivalent source model which is to be solved to minimise the above functional and the fit measure is the standard least squares L 2 norm. Standard algorithmic methods can be applied to adjust the model parameters (i.e. the density distribution p) to achieve a minimum.

For the case of combining gravity and gravity gradient measurements, the filter ensures that the attenuated and distorted high frequencies in the gravity data do not compromise the more high-resolution gravity gradient data.

The equation also incorporates optional regularisation. With the correct choice of this regularisation and/or constraints, the resulting inversion can actually deconvolve the effect of the filtering and recreate bandwidth that was lost in the measurements.

After a successful inversion, i.e. the generation of an optimal equivalent source model, at step S204, the model can be used to forward calculate gravity and gravity gradient components either back to the original measurement positions or at a series of new locations such as level grids for example. Such forward calculations proceed without the introduction of the filter matrices therefore providing the full bandwidth combination of all the measurement data. The forward calculated data may then be processed using any known technique, e.g. those taught in the applicant's earlier applications WO2007/012895, WO2007/085875, WO2008/093139, WO2009016348, WO2008117081 and PCT/GB2008/050041 which are herein incorporated by reference.

More details of the filter are provided for the case in which gravity and gravity gradient data is combined. The filter may be adapted as required to suit other types of data. The filter operates as a time domain convolution specifically for the pre-filtered gravity measurements. Representing (1) using matrix algebra results in

Minimise + 2 {Rp - r) (2).

Here the matrices have been partitioned vertically to show how the gravity, m g and gravity gradient, moo measurements are handled differently. A G represents the superposition matrix of Green's function integrals that when multiplied by the source density vector p generates the set of gravity predictions from the model. In a similar way, A GG p gives the gravity gradient predictions. F is the filter convolution matrix that replicates the filtering already performed on the gravity measurements m g . This commonly takes the form of a low pass filter, i.e. a filter that passes low-frequency signals but attenuates signals with frequencies higher than the cutoff frequency. Each row of F contains the weights that average together the predicted measurements from the source model. The gravity gradient model predictions in (2) are not filtered in this case and accordingly, the identity matrix / leaves them unchanged. Thus, the inversion uses only low frequency gravity data and full range gravity gradient data.

If parts of the gravity gradient data spectrum are known to suffer from extraneous noise, then a filter could be designed to remove this noise from the measurements and subsequently be built into the inversion as another filter matrix in place of /. The noisy part of this spectrum would not then compromise the inversion and the missing bandwidth would be made up by the gravity data.

A set of weights generating the effect of an exponential low pass filter with

characteristic time T c could be formed by the following

where ti j is the time between measurement points i and j and n is the filter order. For a first order filter (n = 1), above the cut-off frequency, the signal amplitude is reduced by a factor e = 2.718 every time the frequency doubles. A second order frequency attenuates higher frequencies more steeply.

Regularization is standard in these inversion problems, depicted in (2) by the matrix R and the vector r, can take many forms and aims to stabilise and control the solution.

Common choices include Tikhonov regularization where λ controls the relative amounts in which the solution fits the data and a priori information specified by R and r (see for example Gravity Interpretation Fundamentals and Application of Gravity Inversion and Geological Interpretation Jacoby W., Smilde P.L, 2009, Springer, Berlin). Other forms of regularisation can attempt to smoothen the density distribution in space by making the matrix R a gradient operator.

The gravity and gravity gradient data is inverted into the same model described by the vector of parameters p. The minimisation (2) can be achieved using standard

optimisation algorithms such as the conjugate gradient method or Monte Carlo methods (see for example: Numerical Recipes in C, 2 nd Edition, Press W.H., et. al. 1997, Cambridge University Press.) Figure 4a shows a synthetic model having dimensions 150 x 150 km which was constructed to illustrate the advantages of combining gravity and gravity gradient data for a large survey. The model includes an approximately 200 km feature in the data which peaks along the y = x diagonal. Forward calculations from this model to a set of survey lines spaced at 5 km were performed to simulate the gravity and gravity gradient signals that would be measured during a survey. Noise appropriate to the instrumentation and data correction errors were added and subsequent filtering applied to get as close as possible to realistic measurement data. The expected gravity measurements as predicted by these calculations are shown in Figure 4b.

Figure 5a shows the results of using only the simulated gravity gradient measurements in an equivalent source inversion. The resulting error distribution in the predicted gravity distribution is shown in Figure 5 a. Error plots of the predicted gravity gradient signals (not shown) reveal only small uncorrelated error distributions. The correlated nature of the gravity error highlights the deficiencies of accurate long wavelength signal in the gravity gradient measurements resulting in poor imaging of the 200 km feature peaking along the y = x diagonal.

Figure 5b shows the results of incorporating the simulated gravity measurements into the inversion and shows the resulting gravity error distribution. As shown, the error in the gravity becomes essentially uncorrelated and lower in amplitude showing the benefit of combining the data with the filter augmented inversion as described above. The predictions of the gravity gradient fields after incorporating the gravity data remain essentially unchanged showing that the high-resolution part of the inversion is not distorted by the filtered gravity measurements.

In summary, the signal to noise ratio provided by gravity gradiometer surveys is adequate to provide an accurate complete picture of the geological signal for moderately sized surveys supporting wavelengths up to approximately 200 km. For larger regional style surveys, the lack of gradient signal power at low frequencies starts to compromise the recovery of the long wavelength geological features. To recover this missing part of the spectrum, equivalent source inversions combining gravity gradient and gravity measurements are used in an optimum way to yield superior results giving an accurate signal over a wide bandwidth. To correctly invert measurement sets that have undergone significant time domain filtering, the equivalent source inversion is augmented with information about the filter. In this way the inversion fits the measurements without distorting the actual inversion model. Such a model can then be used to predict high resolution signals without the affects of filtering.

The gravity and gravity gradient data sets can derive from separate survey systems acquired at different times or simultaneously from within the same system. In the latter case, the gravimeter and gradiometer may be mounted on the same stabilised platform (as shown in figure 1). Alternatively, gravity data can also be provided by a suitable output channel of a gradiometer.

Figure 6a shows an accelerometer module for a gradiometer comprising a pair of proof masses 32 with coils mounted on coil supports 34 acting on surfaces of the proof masses 32. Niobium sensor modules 36 monitor movement of the proof masses. The accelerometer is mounted in a titanium mounting cube 30. Between two and six accelerometer modules may be used in a gradiometer. Figure 6a is taken from Lumley et. al., A superconducting gravity gradiometer tool for exploration (2004) ASEG-PESA Airborne Gravity Workshop, Sydney which describes the operation of the gradiometer in more detail and is hereincorporated by reference.

As set out in the paper, it is necessary to consider what happens when the pancake coil- proof mass system is placed within a closed superconducting loop to understand the operation of the gradiometer. Understanding the operation is further complicated since there are two masses and two loops. Various schematic circuits are illustrated in the paper, with the preferred differential mode circuit shown in Figure 6. The gradiometer described in the paper is designed to utilise the illustrated differential mode circuit to measure the gravity gradient signal. As explained in more detail below, it is possible to adapt such a gradiometer which has a principle form of operation revolving around differencing two or more acceleration measurements taken over a defined baseline to provide a common mode output which is a measure of linear acceleration. For example, in the case of an inline Gzz gradiometer that measures the vertical acceleration experienced at two points separated vertically, taking the difference of the two measurements provides for the gravity gradient whereas the common mode sum corresponds to the vertical acceleration. After conventional gravity corrections, processing, reduction and filtering, this common mode measurement provides the vertical gravity component. The differencing and summing stages are intrinsic to the gradiometer and are most effectively performed within suitably designed circuits that link the signals from each accelerometer module.

Figure 6b of the present application illustrates one such suitably designed circuit which may be used in the gradiometer described in the paper to provide both gravity gradient signals and gravity signals. In the schematic of Figure 6b there are two proof masses mi and m 2 , each interacting with a pair of coils; (C la , Cn,) and (C 2a , C 2b ). A first inductor 38 is connected to a first sensor (not shown) to measure the current across the inductor 38. A second inductor 40 is coupled to a second sensor (not shown) to measure the current in another adjacent loop. As explained in more detail below, the first sensor detects the gravity gradient and is operating in differential model; the second sensor detects the gravity signal and is operating in common model.

Both masses are supported by suspension structures which provide mass-on-spring type behaviour where the displacements of the proof masses are proportional to the total acceleration they experience. The proof masses and circuits are superconducting whereby persistent currents / / and h can be stored and provide sensitivity to the proof mass displacements. For example, if proof mass mi moves upwards towards coil Cn,, the inductance of coil Cn, will decrease and the inductance of coil C la will increase. By applying the rules of flux conservation in superconducting loops, (e.g. as described in Lumley et. al. (2004)), such a change in coil inductances will cause a redistribution of currents resulting in a signal current i \ commensurable with the proof mass

displacement.

If the gradiometer is subjected to a change in gravity, this causes displacements of both proof masses in the same direction. Flux quantisation for each superconducting loop requires both Ii and I 2 to increase. Since these two currents flow in opposite directions through the inductor 38, the changes tend to cancel one another and thus there is no current flowing through the central arm of the circuit. Thus a sensor coupled to this inductor measures no current ¾ and does not see the common mode gravity signal.

However, the current flowing through the outer loop is the average of Ii and I 2 and results in currents i \ and i 2 . These currents are measurable by providing inductors 41 and 42 which both couple to another loop which measures a common mode response current i c which is a sum of i \ and i 2 .

In the presence of a gravity gradient, the gravitational acceleration at proof mass mi is different from that at proof mass m 2. If the differential gravity causes the top proof mass to move upwards and the bottom mass to move downwards, Ii decreases and I 2 increases. Thus, currents ij and i 2 have opposite sense resulting in a net current flowing through the central arm of the circuit which may be termed a differential mode response current i d .

By monitoring currents i c and ¾ one therefore has simultaneous measures of gravity and gravity gradient changes from within the same sensor.

Figure 7 shows a gravity instrument 50 and a gravity gradiometer 52 mounted on a stabilised platform. The platform is supported on a support frame 54 which comprises a base 56, four pillars 58 mounted generally perpendicular to the base, one at each corner and an outer frame 60 supported on the pillars 58. As detailed above, gravity has three components with respect to any chosen Cartesian coordinate system which may be the three mutually perpendicular axes (X, Y, Z) shown in Figure 7. These axes are sometimes referred to as the yaw, pitch and roll axes. Alternatively, a different set of axes may be defined.

For accurate measurements of both gravity and gravity gradient, the measuring instruments need to be stabilised relative to the axis (or axes) along which measurements are taken. Figure 7 shows a pair of instruments that are stabilised with respect to all three axes so that the instruments are effectively mounted on a stable platform which remains fixed in inertial space. However, this level of stabilisation may not be required for all applications. Stabilisation in each axis is provided by a gimbal assembly comprising a gimbal driven by a respective torque motor. The torque motor controls the orientation of the respective gimbal. As shown in Figure 7, the X-axis torque motor 62 is mounted to the outer frame and connected to the X-axis gimbal 64. The X-axis gimbal 64 fits within the outer frame and has a similar shape to the outer frame; both are generally square. The Y-axis torque motor 66 is mounted to the X-axis gimbal 64 and connected to the Y-axis gimbal 68. The Y-axis gimbal fits within the X-axis gimbal and is an annular ring. The Z-axis gimbal assembly 70 is mounted within the Y-axis gimbal. The gravity instrument 50 and gravity gradiometer 52 are mounted on the Z-axis gimbal assembly so that they are stabilised at least with respect to this single axis. Optionally, the instruments may be stabilised with respect to all axes by activation of the appropriate gimbal and torque motor.

The instruments are fixed relative to each other and are thus acquire data at the same time from the same point in space. Information on this location and time may be captured by a GPS 74 or similar instrument. All the instruments have data outputs which feed to a data acquisition module 72 for processing as described above. If the gravity and gravity gradient measurements are not acquired from the same location at the same time, corrections are required before the data can be processed together. However, these corrections are not required with the system shown in Figure 7.

No doubt many other effective alternatives will occur to the skilled person. It will be understood that the invention is not limited to the described embodiments and encompasses modifications apparent to those skilled in the art lying within the spirit and scope of the claims appended hereto.