WO2003052460A1 - Method of processing marine magnetic gradient data and exploration methods using that data - Google Patents
Method of processing marine magnetic gradient data and exploration methods using that data Download PDFInfo
- Publication number
- WO2003052460A1 WO2003052460A1 PCT/AU2002/001669 AU0201669W WO03052460A1 WO 2003052460 A1 WO2003052460 A1 WO 2003052460A1 AU 0201669 W AU0201669 W AU 0201669W WO 03052460 A1 WO03052460 A1 WO 03052460A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- ship
- data
- gradient
- bias
- sensors
- Prior art date
Links
- 230000005291 magnetic effect Effects 0.000 title claims abstract description 159
- 238000000034 method Methods 0.000 title claims abstract description 67
- 238000012545 processing Methods 0.000 title claims abstract description 30
- 101150080778 INPP5D gene Proteins 0.000 claims description 173
- 238000005553 drilling Methods 0.000 claims description 14
- 230000010354 integration Effects 0.000 claims description 10
- 230000007613 environmental effect Effects 0.000 claims description 9
- 230000035899 viability Effects 0.000 claims description 4
- 230000000694 effects Effects 0.000 description 6
- 238000000926 separation method Methods 0.000 description 5
- 238000003672 processing method Methods 0.000 description 4
- 238000010420 art technique Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 101100458361 Drosophila melanogaster SmydA-8 gene Proteins 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 150000004649 carbonic acid derivatives Chemical class 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 231100001261 hazardous Toxicity 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000005070 sampling 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/15—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for use during transport, e.g. by a person, vehicle or boat
-
- 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/15—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for use during transport, e.g. by a person, vehicle or boat
- G01V3/165—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for use during transport, e.g. by a person, vehicle or boat operating with magnetic or electric fields produced or modified by the object or by the detecting device
-
- 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/38—Processing data, e.g. for analysis, for interpretation, for correction
Definitions
- This invention relates to a method of processing marine magnetic gradient data and also to exploration methods using that data.
- seismic techniques provide an indication of various subsurface structures including porous sandstone or fractured carbonates which may contain oil or other valuable deposits.
- a magnetic gradient survey enables magnetic data relating to the survey region to be obtained and this magnetic data can be used to provide information relating to the nature of subsurface structures. If the magnetic survey is overlayed with the seismic survey, then structures which appear to be of interest from the point of view of the seismic survey, can be further considered in the light of the magnetic data so that a clearer indication can be formed as to whether the structure is a structure which may contain a deposit of interest such as an oil deposit, or whether the structure exhibits magnetic phenomena, such would be the case with a volcano. Thus, the location of drilling platforms can be decided with more precision to avoid subsurface structures which could be hazardous from the point of view of a drilling operation. However, conventional techniques for processing magnetic gradient data contain considerable distortion and anomalies due to unwanted magnetic effects including the ship bias referred to above .
- the conventional method of obtaining that data in relation to onshore exploration is to tow magnetometers behind an aircraft to obtain magnetic data in relation to the survey area.
- the magnetometers are towed behind a ship.
- two magnetic field sensors which are generally referred to as fish, are towed behind a ship.
- the magnetometers are connected to a tow line and the magnetometer which is closest to the ship is towed at a distance of some 300-600 metres behind the ship to avoid the effect of ship induced magnetic field. Furthermore, the magnetometers are separated by a distance of in excess of 100 metres. The reason for the length of the tow line and the separation of the magnetometers is to reduce ship bias and therefore provide data which is relatively free of that bias.
- the processing techniques used to acquire magnetic gradient data assume that after a given time period, the trailing magnetometer will be towed to a position which coincides with a previous position of the magnetometer which is closest to the ship. This assumption is made during the processing of the magnetic data.
- the tow line is extremely long and the distance of separation between the magnetometers is considerable, drift of the magnetometers due to sea currents and the like means it is unlikely that the trailing magnetometer will actually occupy the same position as the first magnetometer after that given time period.
- the ship is required to travel along predetermined survey lines and because of the drifting of the sensors, the sensors do not actually travel along the survey line.
- the magnetometers drift to the north of the survey line, the recorded gradient data between the two magnetometers will have larger distortion than if drifting towards the south of the survey line.
- the object of the invention is to provide an improved processing method which accounts for ship-induced bias and which therefore enables the magnetometers to be towed by a tow line which is much shorter than the conventional tow line and also be located closer together to reduce the amount of drift and also increase the likelihood that the trailing magnetometer will occupy the same position which was previously occupied by the leading magnetometer after a given time period.
- the invention may be said to reside in a method of processing marine magnetic data obtained by towing first and second spaced apart sensors behind a ship7 the method including: obtaining raw magnetic gradient data from the sensors; determining the trend of the gradient of the ship bias detected by the sensors; subtracting the trend from the raw magnetic gradient data to obtain corrected gradient data; and processing the corrected gradient data to provide a data output.
- the method determines the trend of the gradient of the ship bias, a more accurate estimation of the actual gradient bias introduced by the ship, having regard to any drifting of the sensors, is taken into account.
- the bias introduced by the ship can be removed from the data so that more accurate data is obtained.
- the present method eliminates the effect of the ship bias, the sensors can be towed much closer to the ship and at much closer separation distances, thereby reducing the amount of drift of the sensors and therefore resulting in the sensors being more likely to follow the actual survey line upon which the ship travels and also for the rear sensor to occupy the same position at a later time which was previously occupied by the front sensor.
- the quality of the output data is greatly improved and therefore provides data which gives a much more accurate indication of the subsurface magnetic characteristics of the survey area, which can be then used in combination with seismic data to make an assessment as to the viability of further exploration or a drilling, program.
- the method includes the step of determining an estimate of the gradient of the ship bias from the raw magnetic gradient data obtained by the sensors, and determining the trend of the gradient of the ship bias from that estimate of the gradient of the ship bias.
- the raw magnetic gradient data is calculated as the measured magnetic signal at the leading sensor less the measured magnetic signal at the trailing sensor, divided by the distance between the sensors .
- the estimate of the gradient of the ship bias may in fact simply be derived from the raw magnetic gradient data which is obtained from the sensors and which may include data relating to the gradient of the ship bias as well as data relating to other magnetic effects.
- the trend of the gradient of the ship bias is determined by applying the gradient of the ship bias data to a B-Spline filter with smoothness 0.8.
- the processing comprises the application of the corrected gradient data over sample intervals, integration of the corrected gradient data into total magnetic intensity data, and application of a low pass filter to the integrated total magnetic intensity data to provide the data output.
- the total magnetic intensity is obtained by integrating the corrected gradient data in accordance with the following equation:
- TMI ⁇ [G c ⁇ ,(t)] + M, ⁇ ,
- G c is the corrected gradient data obtained after subtraction of the trend of the gradient of the ship bias from the raw gradient data
- M s a is the total magnetic intensity value at a point representing the start of the survey or at a location where a survey line and a tie line cross each other.
- TMI TMI ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
- the gradient of the ship bias is calculated according to the following equation:
- M j (x) M e (x) + D(t ⁇ ) + M b (t x )
- M r (x) M e (x) + D(t 2 ) + M b (t 2 )
- M f is the magnetic field as measured by the front sensor and which consists of the environmental field M ⁇ (x), diurnal variation D(t ⁇ ), ship bias M b (t ⁇ ) caused by ship-induced field, sensor drifting and sensor perturbation, all at time ti and along track distance x, and at some later time t 2
- ⁇ l is the distance between the sensors
- D(t ) is the diurnal variations sensed by the rear sensor
- M b (t 2 ) is the ship bias at time t 2 detected by the rear sensor M r .
- the first and second sensors which are towed behind the ship are included in a group of three or more towed sensors.
- the number of sensors towed behind the ship comprises three sensors.
- data from any two of the sensors is used to provide the raw magnetic gradient data.
- the embodiment may obtain raw gradient data from all three sensors and determine the trend of the gradient of the ship bias detected by all three sensors.
- the invention may also be said to reside in a method of obtaining gradient data for an exploration program, the method including: towing first and second sensors behind a ship along predetermined survey lines; obtaining raw magnetic gradient data from the sensors; determining the trend of the gradient of the ship bias detected by the sensors; subtracting the trend from the raw magnetic gradient data to obtain corrected gradient data; and processing the corrected gradient data to provide a data output.
- the method includes the step of determining an estimate of the gradient of the ship bias from the raw magnetic gradient data obtained by the sensors, and determining the trend of the gradient of the ship bias from that estimate of the gradient of the ship bias.
- the raw magnetic gradient data is calculated as the measured magnetic signal at the leading sensor less the measured magnetic signal at the trailing sensor, divided by the distance between the sensors.
- the estimate of the gradient of the ship bias may in fact simply be derived from the raw magnetic gradient data which is obtained from the sensors and which may include data relating to the gradient of the ship bias as well as data relating to other magnetic effects.
- the trend of the gradient of the ship bias is determined by applying the gradient of the ship bias data to a B-Spline filter with smoothness 0.8.
- the further processing comprises the application of the corrected gradient data over sample intervals, integration of the corrected gradient data into total magnetic intensity data, and application of a low passed filter to the integrated total magnetic intensity data to provide the data output.
- the total magnetic intensity is obtained by integrating the corrected gradient data in accordance with the following equation:
- M st at is the total magnetic intensity value at a point representing the start of the survey or at a location where a survey line and a tie line cross each other.
- TM so calculated are then smoothed to remove residual noise using a low pass filter.
- the gradient of the ship bias is calculated according to the following equation:
- M j (x) M e (x) + D(t ⁇ ) + M b
- M f is the magnetic field as measured by the front sensor and which consists of the environmental field M e (x), diurnal variation D(t ⁇ ), ship bias M b (tj . ) caused by ship-induced field, sensor drifting and sensor perturbation, all at time ti and along track distance x, and at some later time t 2
- ⁇ l is the distance between the sensors
- D(t 2 ) is the diurnal variations sensed by the rear sensor
- M b (t 2 ) is the ship bias at time t 2 detected by the rear sensor M r .
- the first and second sensors which are towed behind the ship are included in a group of three or more towed sensors .
- the number of sensors towed behind the ship comprises three sensors.
- data from any two of the sensors is used to provide the raw magnetic gradient data.
- the embodiment may obtain raw gradient data from all three sensors and determine the trend of the gradient of the ship bias detected by all three sensors.
- the invention may also be said to reside in a method of subsurface exploration to determine viability of drilling in a marine environment by considering magnetic data relating to the environment, and wherein the magnetic data has been obtained by: obtaining raw magnetic gradient data from the sensors; determining the trend of the gradient of the ship bias detected by the sensors; subtracting the trend from the raw magnetic gradient data to obtain corrected gradient data; and processing the corrected gradient data to provide a data output.
- obtaining the magnetic data includes the step of determining the gradient of the ship bias from data obtained by the sensors, and determining the trend of the gradient of the ship bias from the gradient of the ship bias.
- the raw magnetic gradient data is calculated as the measured magnetic signal at the leading sensor less the measured magnetic signal at the trailing sensor, divided by the distance between the sensors.
- the estimate of the gradient of the ship bias may in fact simply be derived from the raw magnetic gradient data which is obtained from the sensors and which may include data relating to the gradient of the ship bias as well as data relating to other magnetic effects.
- the trend of the gradient of the ship bias is determined by applying the gradient of the ship bias data to a B-Spline filter with smoothness 0.8.
- the processing comprises the application of the corrected gradient data over sample intervals, integration of the corrected gradient data into total magnetic intensity data, and application of a low passed filter to - li the integrated total magnetic intensity data to provide the data output .
- the total magnetic intensity is obtained by integrating the corrected gradient data in accordance with the following equation:
- Msta is the total magnetic intensity value at a point representing the start of the survey or at a location where a survey line and a tie line cross each other.
- TMI TMI ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
- the gradient of the ship bias is calculated according to the following equation:
- M j (x) M e (x) + D( ) + M b (t ⁇ )
- M r (x) M e (x) + D(t 2 ) + M b (t 2 )
- M £ is the magnetic field as measured by the front sensor and which consists of the environmental field M ⁇ (x), diurnal variation D(t ⁇ ), ship bias M b (t ⁇ ) caused by ship-induced field, sensor drifting and sensor perturbation, all at time t x and along track distance x, and at some later time t 2/ ⁇ l is the distance between the sensors, D(t 2 ) is the diurnal variations sensed by the rear sensor, and
- M b (t 2 ) is the ship bias at time t 2 detected by the rear sensor M r .
- the first and second sensors which are towed behind the ship are included in a group of three or more towed sensors.
- the number of sensors towed behind the ship comprises three sensors.
- data from any two of the sensors is used to provide the raw magnetic gradient data.
- the embodiment may obtain raw gradient data from all three sensors and determine the trend of the gradient of the ship bias detected by all three sensors.
- the invention may still further be said to reside in a method of drilling for a deposit in a marine environment, including: determining the location of drilling from data which has been obtained and which indicates the possible existence of the deposit; and which location is also determined by magnetic data which has been obtained by towing magnetic sensors behind a ship, the magnetic data being processed by: obtaining raw magnetic gradient data from the sensors; determining the trend of the gradient of the ship bias detected by the sensors; subtracting the trend from the raw magnetic gradient data to obtain corrected gradient data; and processing the corrected gradient data to provide a data output .
- the method includes the step of determining an estimate of the gradient of the ship bias from the raw magnetic gradient data obtained by the sensors, and determining the trend of the gradient of the ship bias from that estimate of the gradient of the ship bias.
- the raw magnetic gradient data is calculated as the measured magnetic signal at the leading sensor less the measured magnetic signal at the trailing sensor, divided by the distance between the sensors.
- the estimate of the gradient of the ship bias may in fact simply be derived from the raw magnetic gradient data which is obtained from the sensors and which may include data relating to the gradient of the ship bias as well as data relating to other magnetic effects.
- the trend of the gradient of the ship bias is determined by applying the gradient of the ship bias data to a B-Spline filter with smoothness 0.8.
- the processing of the corrected gradient data comprises the application of the corrected gradient data over sample intervals, integration of the corrected gradient data into total magnetic intensity data, and application of a low passed filter to the integrated total magnetic intensity data to provide the data output.
- the total magnetic intensity is obtained by integrating the corrected gradient data in accordance with the following equation:
- G c is the corrected gradient data obtained after subtraction of the ship bias trend from the raw gradient data
- M st a is the total magnetic intensity value at a point representing the start of the survey or at a location where a survey line and a tie line cross each other.
- TMI TMI ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
- the gradient of the ship bias is calculated according to the following equation:
- M J (x) M e (x) + D( ) + M b (t l )
- M r (x) M e (x) + D(t 2 ) + M b (t 2 )
- M f is the magnetic field as measured by the front sensor and which consists of the environmental field Me(x), diurnal variation D(t ⁇ ), ship bias M b (t ⁇ ) caused by ship-induced field, sensor drifting and sensor perturbation, all at time ti and along track distance x, and at some later time t 2 ⁇ l is the distance between the sensors, D(t 2 ) is the diurnal variations sensed by the rear sensor, and
- M b (t 2 ) is the ship bias at time t 2 detected by the rear sensor M r .
- the first and second sensors which are towed behind the ship are included in a group of three or more towed sensors.
- the number of sensors towed behind the ship comprises three sensors .
- data from any two of the sensors is used to provide the raw magnetic gradient data.
- the embodiment may obtain raw gradient data from all three sensors and determine the trend of the gradient of the ship bias detected by all three sensors.
- Figure 1 is a view showing a marine magnetic data survey
- Figure 2 is a flow chart showing the preferred method of processing the marine magnetic gradient data according to the preferred embodiment of the invention
- Figure 2A is a graph showing the function of a low pass filter used in the preferred embodiment
- Figure 3A and Figure 3B are graphs showing data obtained according to the survey and as processed according to the preferred embodiment of the invention
- Figure 4 is a representation of a survey using data created according to a prior art technique
- Figure 5 is a diagram similar to Figure 4, but using data obtained according to the preferred embodiment of the present invention.
- Figure 6 is a view of a further embodiment of the invention.
- a marine survey to acquire magnetic data is illustrated.
- a ship 10 tows first magnetometer Ml and second magnetometer M2 by a tow line 12.
- the ship travels along survey line 14 and the line
- the front sensor measures magnetic field - M j consisting of the environment field M e (x) , the diurnal variation D( t ⁇ ) , ship bias M b ( t_ . ) caused by ship induced field, sensor drifting and sensor perturbation.
- M r (x) M e (x) + D(t 2 ) + M b (t 2 )
- M j (x) - M r (x) [ (t,) - D(t 2 )) + [M b (t ) - M b (t 2 ) ⁇
- Ship bias gradient (Step 2, Figure 2 ) can be expressed as :
- G is the raw magnetic gradient data.
- bias zig-zag will be around its mean values.
- the variable mean value of the bias along a survey line is regarded as bias trend " since the high frequencies of sensor perturbation only create random noise around the bias trend and after integration its effect can be eliminated from raw gradient data.
- TMI Total Magnetic Intensity
- Axi ( t) is the sampling distance along the survey line.
- M s a t is the TMI value at a point of survey start or at the location where the survey line and a tie line cross each other.
- the TMI value calculated in Step 6 is then smoothed using a low pass filter (Step 7) to remove any features which have a greater rate of change of TMI with distance than is expected in the particular survey area.
- a low pass filter An example of the action of this filter is provided in Figure 2A in which trace 50 is the smooth TMI curve and trace 51 is the TMI data prior to smoothing.
- the data output obtained at Step 7 may include line levelling and data griding to provide a final data output.
- Figures 3A and 3B represent an actual example of the preferred embodiment of the invention which relates to raw data obtained from a known region.
- the trace 20 represents the total ship bias.
- the line 21 represents the bias trend and it can be seen from the left hand side of the trace 20 to the right hand side of the trace 20, the bias trend 21 changes non- uniformly about the value of -0.08 on the graph of Figure 3A.
- a particular value for the bias trend can " therefore be obtained for various intervals and this value can be subtracted from the raw data in order to provide the corrected data.
- the trend of the gradient of the ship bias is preferably determined by calculating the gradient of the bias data 20 and then applying that gradient data to a B-Spline filter with smoothness 0.8 in order to provide the representation of the ship bias 21.
- the gradient of the bias data 20 is calculated simply by using the raw magnetic gradient data obtained from the sensors because when that data is applied to the B-Spline filter, all that remains is the gradient of the ship bias component.
- the ship bias 21 can then be subtracted from the calculated raw gradient.
- the trend of the gradient of the ship bias is a non-linear function which is represented by the line 21 in Figure 3A.
- the line is a measure of how the bias changes over time as the sensors are towed behind the ship.
- the trend is not constant or merely an average, but rather represents the fluctuation of the bias gradient and at some times, is above the value of 0.08 and sometimes below that value.
- the trend is determined by applying the gradient of the ship bias to a B-Spline filter with smoothness 0.8.
- the filter could have different smoothness depending on the region from which data is collected and the nature of the data which is collected.
- the function of the filter is to smooth the trace 20 so that some meaningful value of the bias trend at particular time periods can be obtained.
- the filter is effectively determining the peaks and troughs of the trace 20 and fitting a curved line between those peaks and troughs, which gives a measure of how the gradient of the ship bias is changing with time as the sensors are towed behind the ship.
- dashed line 30 represents total magnetic intensity data which is provided according to the prior art processing technique.
- the line 32 represents the integrated total magnetic intensity data from the bias corrected gradient data according to the preferred embodiment of the invention.
- the line 34 represents the diurnal variation at a station some 500 km away from the survey area and line 36 is the observed field data including the diurnal effect.
- the preferred embodiment of the present invention provides data which more accurately reflects the likelihood of magnetic structures which can tlien be used as a basis for determining the viability of a drilling operation in an exploration or recovery program.
- Figure 6 shows a second embodiment of the invention in which like reference numerals indicate like components to those described with reference to Figure 1.
- three magnetometers Mi, M 2 and M 3 are towed by the ship 10 on tow line 12.
- This embodiment therefore builds some redundancy into the system in that, should one of the sensors fails, two sensors are still available in order to provide the required gradient data, and therefore there is considerably less likelihood that a survey will be wasted should a ship complete the survey only to find that one of the magnetometers has not operated properly.
- the distance between the magnetometers Mi and M 2 and the distance between the magnetometers M 2 and M 3 is about 15 metres.
- the towing distance between the ship 10 and the first magnetometer i is preferably about 150 metres or less.
- This embodiment provides the further advantage that any group of two sensors can be used to provide the gradient data, and therefore the gradient data could be provided by the magnetometers M x and M 2 , the magnetometers M 2 and M 3 , or the magnetometers M x and M 3 .
- This embodiment provides the further advantage that all three magnetometers could be used to provide the data to enable the calculation of the raw gradient and the gradient of the ship bias. The use of three magnetometers may improve accuracy because of the increased amount of data which is collected.
- magnetometers could be used if desired, thereby increasing the amount of data collected and reducing the likelihood that a survey will be wasted because at least two of the magnetometers are not operated properly.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU2002366366A AU2002366366A1 (en) | 2001-12-18 | 2002-12-10 | Method of processing marine magnetic gradient data and exploration methods using that data |
CNB028301544A CN100399052C (zh) | 2001-12-18 | 2002-12-10 | 处理海洋磁场梯度数据的方法和使用该数据的勘探方法 |
US10/538,709 US7313495B2 (en) | 2002-12-10 | 2002-12-10 | Method of processing marine magnetic gradient data and exploration methods using that data |
GB0512665A GB2417086B (en) | 2001-12-18 | 2002-12-10 | Method of processing marine magnetic gradient data and exploration methods using that data |
NO20052790A NO20052790L (no) | 2001-12-18 | 2005-06-09 | Method of processing marine magnetic gradient data and exploration methods using that data |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AUPR9588 | 2001-12-18 | ||
AUPR9588A AUPR958801A0 (en) | 2001-12-18 | 2001-12-18 | Method of processing marine magnetic gradient data and exploration methods using that data |
AUPR9721 | 2001-12-21 | ||
AUPR9721A AUPR972101A0 (en) | 2001-12-21 | 2001-12-21 | Method of processing marine magnetic gradiant data and exploration methods using that data |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2003052460A1 true WO2003052460A1 (en) | 2003-06-26 |
Family
ID=25646861
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/AU2002/001669 WO2003052460A1 (en) | 2001-12-18 | 2002-12-10 | Method of processing marine magnetic gradient data and exploration methods using that data |
Country Status (5)
Country | Link |
---|---|
CN (1) | CN100399052C (zh) |
AU (1) | AU2002366366A1 (zh) |
GB (1) | GB2417086B (zh) |
NO (1) | NO20052790L (zh) |
WO (1) | WO2003052460A1 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102006059204A1 (de) * | 2006-12-13 | 2008-06-26 | Bundesanstalt für Geowissenschaften und Rohstoffe | Verfahren und Anlage zur Messung des Erdmagnetfeldes in seinen Komponenten auf offener See |
WO2009062236A1 (en) * | 2007-11-12 | 2009-05-22 | Commonwealth Scientific And Industrial Research Organisation | Method and apparatus for detecting marine deposits |
CN102338890A (zh) * | 2010-10-22 | 2012-02-01 | 中国石油天然气股份有限公司 | 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 |
WO2013138908A1 (en) * | 2012-03-23 | 2013-09-26 | Fugro Canada Corp. | System and method for geophysical surveying using electromagnetic fields and gradients |
CN104297336A (zh) * | 2014-10-22 | 2015-01-21 | 中国地质大学(武汉) | 一种基于埋地钢制管道磁异常提取及解释方法 |
CN114814963A (zh) * | 2022-03-03 | 2022-07-29 | 吉林大学 | 一种基于b样条抽道的时间域电磁数据甚低频噪声压制方法 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100434934C (zh) * | 2006-07-12 | 2008-11-19 | 杨辉 | 重磁延拓回返垂直导数目标优化处理方法 |
CN109839877B (zh) * | 2019-03-03 | 2021-07-27 | 上海卯瑞船舶设备有限公司 | 一种用于船舶远程管理的监控装置及其使用方法 |
CN114415073B (zh) * | 2022-03-29 | 2022-08-09 | 中国人民解放军海军工程大学 | 一种航磁矢量梯度仪误差模型的地面快速校准方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3875497A (en) * | 1973-12-26 | 1975-04-01 | Seatrek Ltd | Waterborne magnetic anomaly detection system and apparatus |
US5218300A (en) * | 1991-01-31 | 1993-06-08 | Sextant Avionique | Method and equipment for noise-reduction when detecting a target by means of a system of several detector elements |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2439035Y (zh) * | 2000-09-18 | 2001-07-11 | 中国地质大学(北京) | 采集海底大地电磁信号的组合电路集 |
-
2002
- 2002-12-10 GB GB0512665A patent/GB2417086B/en not_active Expired - Fee Related
- 2002-12-10 AU AU2002366366A patent/AU2002366366A1/en not_active Abandoned
- 2002-12-10 CN CNB028301544A patent/CN100399052C/zh not_active Expired - Fee Related
- 2002-12-10 WO PCT/AU2002/001669 patent/WO2003052460A1/en not_active Application Discontinuation
-
2005
- 2005-06-09 NO NO20052790A patent/NO20052790L/no not_active Application Discontinuation
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3875497A (en) * | 1973-12-26 | 1975-04-01 | Seatrek Ltd | Waterborne magnetic anomaly detection system and apparatus |
US5218300A (en) * | 1991-01-31 | 1993-06-08 | Sextant Avionique | Method and equipment for noise-reduction when detecting a target by means of a system of several detector elements |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102006059204A1 (de) * | 2006-12-13 | 2008-06-26 | Bundesanstalt für Geowissenschaften und Rohstoffe | Verfahren und Anlage zur Messung des Erdmagnetfeldes in seinen Komponenten auf offener See |
DE102006059204B4 (de) * | 2006-12-13 | 2008-10-30 | Bundesrepublik Deutschland, vertreten durch den Präsidenten der Bundesanstalt für Geowissenschaften und Rohstoffe | Verfahren und Anlage zur Messung des Erdmagnetfeldes auf offener See |
WO2009062236A1 (en) * | 2007-11-12 | 2009-05-22 | Commonwealth Scientific And Industrial Research Organisation | Method and apparatus for detecting marine deposits |
CN102338890A (zh) * | 2010-10-22 | 2012-02-01 | 中国石油天然气股份有限公司 | 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 |
WO2013138908A1 (en) * | 2012-03-23 | 2013-09-26 | Fugro Canada Corp. | System and method for geophysical surveying using electromagnetic fields and gradients |
CN104297336A (zh) * | 2014-10-22 | 2015-01-21 | 中国地质大学(武汉) | 一种基于埋地钢制管道磁异常提取及解释方法 |
CN104297336B (zh) * | 2014-10-22 | 2017-02-08 | 中国地质大学(武汉) | 一种基于埋地钢制管道磁异常提取及解释方法 |
CN114814963A (zh) * | 2022-03-03 | 2022-07-29 | 吉林大学 | 一种基于b样条抽道的时间域电磁数据甚低频噪声压制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1720465A (zh) | 2006-01-11 |
GB2417086B (en) | 2006-08-16 |
NO20052790L (no) | 2005-09-05 |
AU2002366366A1 (en) | 2003-06-30 |
GB0512665D0 (en) | 2005-07-27 |
GB2417086A (en) | 2006-02-15 |
NO20052790D0 (no) | 2005-06-09 |
CN100399052C (zh) | 2008-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cochran et al. | The Southeast Indian Ridge between 88 E and 118 E: Gravity anomalies and crustal accretion at intermediate spreading rates | |
US5986974A (en) | Method for geophysical processing and interpretation using seismic trace difference for analysis and display | |
US7065449B2 (en) | Method and system for evaluating geophysical survey data | |
CN102375155B (zh) | 使用在不同水深度的多个地震源来记录的地震数据的波场反虚反射 | |
US8615362B2 (en) | Near-surface geomorphological characterization based on remote sensing data | |
US20080294393A1 (en) | Near Surface Layer Modeling | |
WO1997033183A1 (en) | Method for geophysical processing and interpretation using instantaneous phase and its derivatives and their derivatives | |
WO2006020662A2 (en) | Method and system for processing geophysical survey data | |
EP2279436B1 (en) | System and technique to determine high order derivatives from seismic sensor data | |
US8717845B2 (en) | Quality-based steering methods and systems for 4D geophysical surveys | |
CN102016643A (zh) | 用于衰减双传感器地震拖缆中的低频噪声的方法 | |
US20130176822A1 (en) | Method for estimating elastic parameters by inverting 4d seismic measurements | |
US20100094556A1 (en) | Terrain correction systems | |
WO2003052460A1 (en) | Method of processing marine magnetic gradient data and exploration methods using that data | |
AU2014203490A1 (en) | Survey techniques using streamers at different depths | |
US20010047245A1 (en) | Method for imaging discontinuites in seismic data | |
CN108983300B (zh) | 一种隧道掘进机施工条件下的瞬变电磁隧道超前预报方法 | |
Bourillet et al. | Swath mapping system processing: Bathymetry and cartography | |
US7313495B2 (en) | Method of processing marine magnetic gradient data and exploration methods using that data | |
CN103543761B (zh) | 控制传感器拖缆的牵引速度的方法和系统 | |
Cowan et al. | Aeromagnetic gradiometers? a perspective | |
US10274635B2 (en) | Joint inversion of subsurface resistivity and noise parameters | |
Nielsen et al. | Integrating ground-penetrating radar and borehole data from a Wadden Sea barrier island | |
US7720607B2 (en) | Method for treating seismic cubes corresponding, for a common zone on the ground, to different source/receiver and/or angle of incidence offset values | |
RU2298815C2 (ru) | Способ обработки морских магнитных градиентных данных и способы поисково-разведочных работ с использованием этих данных |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AK | Designated states |
Kind code of ref document: A1 Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ OM PH PL PT RO RU SC SD SE SG SK SL TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A1 Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LU MC NL PT SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG |
|
DFPE | Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101) | ||
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
WPC | Withdrawal of priority claims after completion of the technical preparations for international publication | ||
WWE | Wipo information: entry into national phase |
Ref document number: 2002366366 Country of ref document: AU |
|
WWE | Wipo information: entry into national phase |
Ref document number: 0512665 Country of ref document: GB |
|
ENP | Entry into the national phase |
Ref document number: 2005121566 Country of ref document: RU Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 761/MUMNP/2005 Country of ref document: IN |
|
WWE | Wipo information: entry into national phase |
Ref document number: 20028301544 Country of ref document: CN |
|
122 | Ep: pct application non-entry in european phase | ||
WWE | Wipo information: entry into national phase |
Ref document number: 2006247857 Country of ref document: US Ref document number: 10538709 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: JP |
|
WWW | Wipo information: withdrawn in national office |
Ref document number: JP |
|
WWP | Wipo information: published in national office |
Ref document number: 10538709 Country of ref document: US |