AU2002366366A1 - 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
- AU2002366366A1 AU2002366366A1 AU2002366366A AU2002366366A AU2002366366A1 AU 2002366366 A1 AU2002366366 A1 AU 2002366366A1 AU 2002366366 A AU2002366366 A AU 2002366366A AU 2002366366 A AU2002366366 A AU 2002366366A AU 2002366366 A1 AU2002366366 A1 AU 2002366366A1
- Authority
- AU
- Australia
- Prior art keywords
- ship
- data
- gradient
- bias
- sensors
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 230000005291 magnetic effect Effects 0.000 title claims description 156
- 238000000034 method Methods 0.000 title claims description 65
- 238000012545 processing Methods 0.000 title claims description 28
- 101150080778 INPP5D gene Proteins 0.000 claims description 174
- 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
- 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
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)
Description
WO 03/052460 PCT/AU02/01669 METHOD OF PROCESSING MARINE MAGNETIC GRADIENT DATA AND EXPLORATION METHODS USING THAT DATA Field of the Invention 5 This invention relates to a method of processing marine magnetic gradient data and also to exploration methods using that data. Background Art 10 Exploration of the earth to locate oil, minerals or other valuable deposits, uses seismic techniques in order to provide information about the subsurface structure of the earth so that an assessment can be made as to whether a particular deposit may or may not be present. In 15 particular, seismic techniques provide an indication of various subsurface structures including porous sandstone or fractured carbonates which may contain oil or other valuable deposits. 20 Whilst this technique does provide an indication of subsurface structures, it is possible a particular structure which may appear to be of interest, in fact represents a significant hazard to drilling. For example, it is possible that subsurface volcanos exhibit very 25 similar characteristics to anticlines which may contain oil. If a drilling operation is initiated and it is found that the anticline is in fact a volcano, then a considerable financial loss results because of the cost of establishing the drilling platform. This is particularly 30 the case in relation to marine exploration because of the significantly higher drilling cost. In order to provide more information concerning the nature of subsurface structures, magnetic data of the exploration 35 area is obtained.
WO 03/052460 PCT/AU02/01669 - 2 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 5 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 10 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 15 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. 20 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. In the case of marine exploration, the 25 magnetometers are towed behind a ship. The usefulness of marine magnetic data which is obtained in this manner has been limited by the quality of the data which is obtained. One of the major problems with marine 30 magnetic data acquisition is interference or so-called ship bias which is created by the magnetic field induced by the ship which tows the magnetometers. In conventional marine magnetic data acquisition, two 35 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 WO 03/052460 PCT/AU02/01669 - 3 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 5 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. However, because the towing line is so long and the 10 distance of separation so great, the degree of drift of the magnetometers in the ocean as the magnetometers are towed behind the ship is considerable. Furthermore, the processing techniques used to acquire magnetic gradient data assume that after a given time period, the trailing 15 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. However, because the tow line is extremely long and the distance of separation 20 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. 25 Furthermore, in the acquisition of the magnetic data, 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. In 30 the southern hemisphere, if 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. 35 Summary of the Invention The object of the invention is to provide an improved processing method which accounts for ship-induced bias and WO 03/052460 PCT/AU02/01669 -4 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 5 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 10 processing marine magnetic data obtained by towing first and second spaced apart sensors behind a ship, the method including: obtaining raw magnetic gradient data from the sensors; 15 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 20 a data output. Because 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 25 drifting of the sensors, is taken into account. By subtracting the trend of the gradient bias from the raw gradient data, the bias introduced by the ship can be removed from the data so that more accurate data is obtained. Because the present method eliminates the 30 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 35 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. Thus, the quality of the WO 03/052460 PCT/AU02/01669 -5 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 5 make an assessment as to the viability of further exploration or a drilling. program. Preferably the method includes the step of determining an estimate of the gradient of the ship bias from the raw 10 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. Preferably the raw magnetic gradient data is calculated as 15 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 20 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. 25 Preferably 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. Preferably the processing comprises the application of the 30 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. 35 Preferably the total magnetic intensity is obtained by integrating the corrected gradient data in accordance with WO 03/052460 PCT/AU02/01669 - 6 the following equation: N TMI= f [GcAx 1 (t)] + Maa i=1 wherein G. is the corrected gradient data 5 obtained after subtraction of the trend of the gradient of the ship bias from the raw gradient data, and Metat 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 10 other. The values of TMI so calculated are then smoothed to remove residual noise using a low pass filter. 15 Preferably the gradient of the ship bias is calculated according to the following equation: Mf(x)-M,(x) [D(t,)-D(t2)]+[Mb(t) -M b(t 2 )] Gbias = Al Al 20 where: M(x) = Me(x) + D(ti) + Mb(tl) Mr(x) = Me(x) + D(t 2 ) + Mb(t2) 25 where Mf 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 Mb(tl) caused by ship-induced field, sensor drifting and sensor 30 perturbation, all at time tj and along track distance x, and at some later time t 2 , Al is the distance between the sensors, D(t 2 ) is the diurnal variations sensed by the rear sensor, and Mb(t2) is the ship bias at time t 2 detected by the 35 rear sensor Mr.
WO 03/052460 PCT/AU02/01669 - 7 In one embodiment of the invention the first and second sensors which are towed behind the ship are included in a group of three or more towed sensors. 5 In this embodiment, preferably the number of sensors towed behind the ship comprises three sensors. According to this embodiment data from any two of the sensors is used to provide the raw magnetic gradient data. 10 However, 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. 15 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; 20 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 25 gradient data to obtain corrected gradient data; and processing the corrected gradient data to provide a data output. Preferably the method includes the step of determining an 30 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. 35 Preferably 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, WO 03/052460 PCT/AU02/01669 - 8 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 5 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. Preferably the trend of the gradient of the ship bias is 10 determined by applying the gradient of the ship bias data to a B-Spline filter with smoothness 0.8. Preferably.the further processing comprises the application of the corrected gradient data over sample 15 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. 20 Preferably the total magnetic intensity is obtained by integrating the corrected gradient data in accordance with the following equation: N TMI= , [GcAxi(t)] + Ms, i=1 25 wherein G is the corrected gradient data obtained after subtraction of the ship bias trend from the raw gradient data, and Metat is the total magnetic intensity value at a point representing the start of the survey or at a 30 location where a survey line and a tie line cross each other. The values of TMI so calculated are then smoothed to remove residual noise using a low pass filter. 35 Preferably the gradient of the ship bias is calculated WO 03/052460 PCT/AU02/01669 -9 according to the following equation: Mf(x)-M,(x) [D(t,)-D(t2)]+[M,(t,)-Mb2(t) ] Gbias - Al Al 5 where: Mx) = Me(x) + D(t) + Mb(t1) Mr(x) = Me(x) + D(t2) + M(t2) 10 where Mf is the magnetic field as measured by the front sensor and which consists of the environmental field Nl(x), diurnal variation D(t 1 ), ship bias Mb(t) caused by ship-induced field, sensor drifting and sensor 15 perturbation, all at time tj and along track distance x, and at some later time t 2 , Al is the distance between the sensors, D(t 2 ) is the diurnal variations sensed by the rear sensor, and Mb(t 2 ) is the ship bias at time t 2 detected by the 20 rear sensor Mr. In one embodiment of the invention the first and second sensors which are towed behind the ship are included in a group of three or more towed sensors. 25 In this embodiment, preferably the number of sensors towed behind the ship comprises three sensors. According to this embodiment data from any two of the 30 sensors is used to provide the raw magnetic gradient data. However, 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. 35 The invention may also be said to reside in a method of subsurface exploration to determine viability of drilling WO 03/052460 PCT/AU02/01669 - 10 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 5 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 10 processing the corrected gradient data to provide a data output. Preferably.obtaining the magnetic data includes the step of determining the gradient of the ship bias from data 15 obtained by the sensors, and determining the trend of the gradient of the ship bias from the gradient of the ship bias. Preferably the raw magnetic gradient data is calculated as 20 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 25 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. 30 Preferably 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. Preferably the processing comprises the application of the 35 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 WO 03/052460 PCT/AU02/01669 - 11 the integrated total magnetic intensity data to provide the data output. Preferably the total magnetic intensity is obtained by 5 integrating the corrected gradient data in accordance with the following equation: N TMI= [GcAxi(t)] + M i=1 wherein G, is the corrected gradient data 10 obtained after subtraction of the ship bias teend from the raw gradient data, and Mstat 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 15 other. The values of TMI so calculated are then smoothed to remove residual noise using a low pass filter. 20 Preferably the gradient of the ship bias is calculated according to the following equation: MI(x)-M r(x) [D(tj)-D(t 2 )]+[Mb(tl)-Mb(t 2 )] Gias = A Al Al 25 where: M(x) = Me(x) + D(t) + Mb(t1) Mr(x) = Me(x) + D(t2) + Mb(t2) 30 where Mf is the magnetic field as measured by the front sensor and which consists of the environmental field
M
0 (x), diurnal variation D(tj), ship bias Mb(tl) caused by ship-induced field, sensor drifting and sensor 35 perturbation, all at time tj and along track distance x, and at some later time t 2 , Al is the distance between the WO 03/052460 PCT/AU02/01669 - 12 sensors, D(t 2 ) is the diurnal variations sensed by the rear sensor, and Mb(t2) is the ship bias at time t 2 detected by the rear sensor Mr. 5 In one embodiment of the invention the first and second sensors which are towed behind the ship are included in a group of three or more towed sensors. 10 In this embodiment, preferably the number of sensors towed behind the ship comprises three sensors. According to this embodiment data from any two of the sensors is used to provide the raw magnetic gradient data. 15 However, 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. 20 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 25 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 30 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; 35 and processing the corrected gradient data to provide a data output.
WO 03/052460 PCT/AU02/01669 - 13 Preferably 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 5 determining the trend of the gradient of the ship bias from that estimate of the gradient of the ship bias. Preferably the raw magnetic gradient data is calculated as the measured magnetic signal at the leading sensor less 10 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 15 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. Preferably the trend of the gradient of the ship bias is 20 determined by applying the gradient of the ship bias data to a B-Spline filter with smoothness 0.8. Preferably the processing of the corrected gradient data comprises the application of the corrected gradient data 25 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. 30 Preferably the total magnetic intensity is obtained by integrating the corrected gradient data in accordance with the following equation: N TMI = 1 [GAx(t)] + M,,,, 35 wherein Go is the corrected gradient data obtained after subtraction of the ship bias trend from the WO 03/052460 PCT/AU02/01669 - 14 raw gradient data, and MNtat 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 5 other. The values of TMI so calculated are then smoothed to remove residual noise using a low pass filter. 10 Preferably the gradient of the ship bias is calculated according to the following equation: M= (x)-Mr(X)_ [D(tI)-D(t2)]+[Mb(t1)-Mb(t 2 )] Al Al 15 where: Mj(x) = Me(x) + D(tl) + Mb(tl) Mr(x) = Me(x) + D(t 2 ) + Mb(t2) 20 where Mf is the magnetic field as measured by the front sensor and which consists of the environmental field Me(X), diurnal variation D(t 1 ), ship bias Mb(tl) caused by ship-induced field, sensor drifting and sensor 25 perturbation, all at time tj and along track distance x, and at some later time t 2 , Al is the distance between the sensors, D(t 2 ) is the diurnal variations sensed by the rear sensor, and Mbl(t 2 ) is the ship bias at time t 2 detected by the 30 rear sensor Mr. In one embodiment of the invention the first and second sensors which are towed behind the ship are included in a group of three or more towed sensors. 35 In this embodiment, preferably the number of sensors towed behind the ship comprises three sensors.
WO 03/052460 PCT/AU02/01669 - 15 According to this embodiment data from any two of the sensors is used to provide the raw magnetic gradient data. 5 However, 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. Brief Description of the Drawings 10 A preferred embodiment of the invention will be described, by way of example, with reference to the accompanying drawings in which: Figure 1 is a view showing a marine magnetic data survey; 15 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; 20 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; 25 Figure 5 is a diagram similar to Figure 4, but using data obtained according to the preferred embodiment of the present invention; and Figure 6 is a view of a further embodiment of the invention. 30 Description of the Preferred Embodiment With reference to Figure 1, a marine survey to acquire magnetic data is illustrated. A ship 10 tows first magnetometer M1 and second magnetometer M2 by a tow line 35 12. The ship travels along survey line 14 and the line 16, which are dotted, represent the ship-induced magnetic field.
WO 03/052460 PCT/AU02/01669 - 16 In Figure 1, it is assumed the ship is in the southern hemisphere and is heading north-northwest, and the sensors are drifting towards the northeast. In this case, the 5 recorded gradient data between sensors M1 and M2 will have a larger distortion than if the sensors were drifting towards the southwest. Ideal conditions would assume that the sensors M1 and M2 10 are towed directly behind the ship on the survey line and that therefore, after a given time period, the sensor M2 would occupy the same position as was previously occupied by the sensor M1. However, due to drifting of the sensors because of ocean currents and the like, this ideal 15 situation in general cannot be achieved. Nevertheless, the processing methods used to process magnetic data do assume that the sensor M2 will occupy the same position as the sensor M1, previously occupied, but at a later time. The processing method according to the preferred 20 embodiment enables ship bias to be accounted for much more accurately than any prior art techniques and therefore enables the sensors to be towed closer to the ship and to be separated by a smaller distance. This therefore provides much greater accuracy in the measurements because 25 the amount of drift is not as great due to the shorter tow line and also because the separation between the sensors is smaller, the likelihood of the rear sensor occupying the same position as the front sensor at a later period in time is greatly increased. 30 At any time tj and along-track distance x, the front sensor measures magnetic field - Mf consisting of the environment field M(x), the diurnal variation D(tj), ship bias 14(tl) caused by ship induced field, sensor drifting and sensor 35 perturbation. M,(x) = Me(x) + D(tj) + Mb(tl) WO 03/052460 PCT/AU02/01669 - 17 At some later time t 2 , the rear sensor makes a measurement at the same spatial point: .5 Mr(x) = M(x) + D(t 2 ) + Mb(t 2 ) The difference between the two measurements at the same location is: 10 Mx) - Mr(x) = [D(t) - D(t 2 )] + [Mb() - Mb(t2)] Ship bias gradient (Step 2, Figure 2) can be expressed as: Gbi. = M (x) - M (x) = [D(t, ) - D(t 2 )] + [b ( 1 ) - MI (t 2 )] Al Al 15 It should be mentioned that the calculated ship bias gradient is still affected by diurnal variation. In practice we notice that the sea current drifting is a kind of long wave-length variation so ship bias should be 20 coincident with such current drifting. Step 2 in Figure 2, which is the calculation of the gradient of the ship bias, is actually an estimate of the gradient of the ship bias which is obtained from the raw data collected by the sensors M1 and M2. The data collected by the sensors will 25 include raw magnetic gradient data which contains many signals including environmental magnetic signals, diurnal signals and the ship bias, as well as instrument bias and drift. The trend of the gradient of the ship bias is obtained from this estimate by applying the estimate of 30 the gradient of the ship bias to a B-Spline filter, as will be described in more detail hereinafter. Therefore the trend of the gradient of the ship bias (Gtrend) (Step 3, Figure 2) is used for the correction along a track line. 35 The gradient correction of the ship bias (Gc)can be WO 03/052460 PCT/AU02/01669 - 18 expressed as follows: Gc = G - Gtre,,d (Step 4, Figure 2) 5 In this equation, G is the raw magnetic gradient data. As said before such drifting can cause large bias effects. Consequently, when the sensor drifting and perturbation take place simultaneously the bias zig-zag will be around 10 its mean values. The variable mean value of the bias along a survey line is regarded as bias trend'gince 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. 15 After removal of ship bias effects TMI (Total Magnetic Intensity) can be calculated (Step 6) through integration of magnetic gradient data: 20 TMI= i ' [GcAxi(t)] + Mstat Here, Axi(t) is the sampling distance along the survey line. Mta is the TMI value at a point of survey start or at the location where the survey line and a tie line cross 25 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 30 expected in the particular survey area. 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. 35 The data output obtained at Step 7 may include line levelling and data griding to provide a final data output.
WO 03/052460 PCT/AU02/01669 - 19 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. 5 In Figure 3A, 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 10 uniformly about the value of -0.08 on the graph of Figure 3A. A particular value for the bias trend caf-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 15 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. 20 In the preferred embodiment of the invention, 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 25 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. 30 The line is a measure of how the bias changes over time as the sensors are towed behind the ship. As is clearly seen from line 21, 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 35 sometimes below that value. As noted above, the trend is determined by applying the gradient of the ship bias to a B-Spline filter with smoothness 0.8. However, in other WO 03/052460 PCT/AU02/01669 - 20 embodiments, the filter could have different smoothness depending on the region from which data is collected and the nature of the data which is collected. In general, the function of the filter is to smooth the trace 20 so 5 that some meaningful value of the bias trend at particular time periods can be obtained. Thus, 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 10 ship bias is changing with time as the sensors are towed behind the ship. In the example of data obtained from the known region, dashed line 30 represents total magnetic intensity data 15 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 20 at a station some 500 km away from the survey area and line 36 is the observed field data including the diurnal effect. It can be seen from the conventional processing technique, 25 which gives line 30, false anomalies are provided which do not show up on the line 32, which is produced according to the present invention. In general, the results of the magnetic survey are 3,0 produced in a colour chart. Black and white representations of a colour chart are shown as Figure 4 in relation to conventional processing of data obtained from a known region and Figure 5 by a processing method according to the preferred embodiment of this invention. 35 The characteristics of the Igneous basement of the known region are well known and it can be seen that data WO 03/052460 PCT/AU02/01669 - 21 processed according to the present invention produces a much clearer indication of actual magnetic structures than prior art technique which include considerable interference and data which may mislead analysers into 5 considering that magnetic structures exist which are not actually there, or that magnetic structures do not exist when a magnetic structure actually is present. Thus, the preferred embodiment of the present invention 10 provides data which more accurately reflects the likelihood of magnetic structures which can tlen be used as a basis for determining the viability of a drilling operation in an exploration or recovery program. 15 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. In this embodiment, three magnetometers M 1 , M 2 and M 3 are towed by the ship 10 on tow line 12. This embodiment therefore 20 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 25 that one of the magnetometers has not operated properly. In this embodiment, preferably the distance between the magnetometers M, and M2, and the distance between the magnetometers M 2 and M3 is about 15 metres. The towing distance between the ship 10 and the first magnetometer M, 30 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 35 the magnetometers MI and M2, the magnetometers M 2 and M3, or the magnetometers M, and M 3 . This embodiment provides the further advantage that all three magnetometers could be WO 03/052460 PCT/AU02/01669 - 22 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. 5 Obviously, more than three 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 10 properly. Since modifications within the spirit and scope of the invention may readily be effected by persons skilled within the art, it is to be understood that this invention 15 is not limited to the particular embodiment described by way of example hereinabove.
Claims (36)
1. A method of processing marine magnetic data obtained by towing first and second spaced apart sensors 5 behind a ship, 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; 10 subtracting the trend from the raw magnetic gradient data to obtain corrected gradient dafa; and processing the corrected gradient data to provide a data output. 15
2. The method of claim 1 including 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 20 bias.
3. The method of claim 1 wherein the raw magnetic gradient data is calculated as the measured magnetic signal at the leading sensor less the measured magnetic 25 signal at the trailing sensor, divided by the distance between the sensors.
4. The method of claim 1 wherein the trend of the gradient of the ship bias is determined by applying the 30 gradient of the ship bias data to a B-Spline filter with smoothness 0.8.
5. The method of claim 1 wherein the processing comprises the application of the corrected gradient data 35 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 WO 03/052460 PCT/AU02/01669 - 24 magnetic intensity data to provide the data output.
6. The method of claim 5 wherein the total magnetic intensity is obtained by integrating the corrected 5 gradient data in accordance with the following equation: N TMI= [GcAx,(t)] + Mstat i=1 wherein Gc is the corrected gradient data obtained after subtraction of the trend of the gradient of 10 the ship bias from the raw gradient data, and MBtat 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. 15
7. The method of claim 1 wherein the gradient of the ship bias is calculated according to the following equation: MI (x)-M,(x)_ [D(t ) - D(t 2 )] +[Mh(tI)- Mb(t2)] 20 Gbias Al Al where: M/x) = Me(x) + D(ti) + Mb(tI) 25 M,(x) = Me(x) + D(t2) + Mb(t2) where Mf is the magnetic field as measured by the front sensor and which consists of the environmental field 30 M4.(x), diurnal variation D(tj), ship bias Mb(tl) caused by ship-induced field, sensor drifting and sensor perturbation, all at time tj and along track distance x, and at some later time t 2 , Al is the distance between the sensors, D(t 2 ) is the diurnal variations sensed by the rear 35 sensor, and Mb(t2) is the ship bias at time t 2 detected by the WO 03/052460 PCT/AU02/01669 - 25 rear sensor Mr.
8. The method of claim 1 wherein the first and second sensors which are towed behind the ship are 5 included in a group of three or more towed sensors.
9. The method of claim 8 wherein the number of sensors towed behind the ship comprises three sensors.
10 10. 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 15 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 20 processing the corrected gradient data to provide a data output.
11. The method of claim 10 wherein the method includes the step of determining an estimate of the 25 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. 30
12. The method of claim 10 wherein 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. 35
13. The method of claim 10 wherein the trend of the gradient of the ship bias is determined by applying the WO 03/052460 PCT/AU02/01669 - 26 gradient of the ship bias data to a B-Spline filter with smoothness 0.8.
14. The method of claim 10 wherein the further 5 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 10 data output.
15. The method of claim 14 wherein the total magnetic intensity is obtained by integrating the corrected gradient data in accordance with the following equation: 15 N TMI= [ [GAx;(t)] + MIa i=1 wherein G. is the corrected gradient data obtained after subtraction of the ship bias trend from the raw gradient data, and 20 Mstat 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. 25
16. The method of claim 10 wherein the gradient of the ship bias is calculated according to the following equation: Mf (x) - M, (x) [D(t,) -D(t2)]+ [MI, (t,)- M (tz) ] Gbias Al Al 30 where: Mj(x) = Me(x) + D(t ) + Mb(tI) 3 5 Mr(x) = Me(x) + D(t) + M,(t2) WO 03/052460 PCT/AU02/01669 - 27 where Mf is the magnetic field as measured by the front sensor and which consists of the environmental field M.(x), diurnal variation D(tj), ship bias Mb(tl) caused by ship-induced field, sensor drifting and sensor 5 perturbation, all at time tj and along track distance x, and at some later time t 2 , Al is the distance between the sensors, D(t 2 ) is the diurnal variations sensed by the rear sensor, and Mbl(t 2 ) is the ship bias at time t 2 detected by the 10 rear sensor Mr.
17. The method of claim 10 wherein the first and second sensors which are towed behind the ship are included in a group of three or more towed sensors. 15
18. The method of claim 17 wherein the number of sensors towed behind the ship comprises three sensors.
19. A method of subsurface exploration to determine 20 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 a plurality of sensors; 25 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 30 a data output.
20. The method of claim 19 wherein the magnetic data includes the step of determining the gradient of the ship bias from data obtained by the sensors, and determining 35 the trend of the gradient of the ship bias from the gradient of the ship bias. WO 03/052460 PCT/AU02/01669 - 28
21. The method of claim 19 wherein 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 5 between the sensors.
22. The method of claim 19 wherein 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 10 smoothness 0.8.
23. The method of claim 19 wherein the processing comprises the application of the corrected gradient data over sample intervals, integration of the corrected 15 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.
24. The method of claim 23 wherein the total magnetic 20 intensity is obtained by integrating the corrected gradient data in accordance with the following equation: N TMI= [ [GcAxi(t)] + M,, i=1 wherein Go is the corrected gradient data 25 obtained after subtraction of the ship bias trend from the raw gradient data, and Mtat 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 30 other.
25. The method of claim 19 wherein the gradient of the ship bias is calculated according to the following equation: 35 WO 03/052460 PCT/AU02/01669 - 29 Mf(x)-M,(x) [D(ti)-D(t 2 )]+[Mb(t)- Mb(t 2 )] Gbias = = Al Al where: 5 M/x) = Me(x) + D(ti) + Mb(tl) Mr(x) = Me(x) + D(t 2 ) + Mb(t2) where Mf is the magnetic field as measured by the 10 front sensor and which consists of the environmental field M.(x), diurnal variation D(t 1 ), ship bias Mb(tj caused by ship-induced field, sensor drifting and sensor perturbation, all at time tj and along track distance x, and at some later time t 2 , Al is the distance between the 15 sensors, D(t 2 ) is the diurnal variations sensed by the rear sensor, and Mb(t2) is the ship bias at time t 2 detected by the rear sensor M,. 20
26. The method of claim 19 wherein the number of sensors towed behind the ship comprises three sensors.
27. A method of drilling for a deposit in a marine environment, including: 25 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 30 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; 35 subtracting the trend from the raw magnetic gradient data to obtain corrected gradient data; and WO 03/052460 PCT/AU02/01669 - 30 processing the corrected gradient data to provide a data output.
28. The method of claim 27 wherein the method 5 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. 10
29. The method of claim 27 wherein the faw 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 15 between the sensors.
30. The method of claim 27 wherein 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 20 smoothness 0.8.
31. The method of claim 27 wherein the processing of the corrected gradient data comprises the application of the corrected gradient data over sample intervals, 25 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. 30
32. The method of claim 31 wherein the total magnetic intensity is obtained by integrating the corrected gradient data in accordance with the following equation: N TMI = I [Godxi(t)] + M,,, i=j 35 wherein Go is the corrected gradient data obtained after subtraction of the ship bias trend from the WO 03/052460 PCT/AU02/01669 - 31 raw gradient data, and Mtat 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 5 other.
33. The method of claim 27 wherein the gradient of the ship bias is calculated according to the following equation: 10 M f(x) - Ml(x) [D(t, ) - D(t 2 b b 2 )] bias Al Al where: 15 M(x) = Me(x) + D(t 1 ) + Mb(tI) Mr(x) = Me(x) + D(t 2 ) + Mb(t2) where Mf is the magnetic field as measured by the 20 front sensor and which consists of the environmental field NMe(x), diurnal variation D(ta), ship bias Mb(tl) caused by ship-induced field, sensor drifting and sensor perturbation, all at time tj and along track distance x, and at some later time t 2 , Al is the distance between the 25 sensors, D(t 2 ) is the diurnal variations sensed by the rear sensor, and Mb(t 2 ) is the ship bias at time t 2 detected by the rear sensor Mr. 30
34. The method of claim 27 wherein the first and second sensors which are towed behind the ship are included in a group of three or more towed sensors.
35. The method of claim 34 wherein the number of 35 sensors towed behind the ship comprises three sensors.
36. The method of claim 35 wherein data from any two WO 03/052460 PCT/AU02/01669 - 32 of the sensors is used to provide the raw magnetic gradient data.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| AUPR9588A AUPR958801A0 (en) | 2001-12-18 | 2001-12-18 | Method of processing marine magnetic gradient data and exploration methods using that data |
| AUPR9721A AUPR972101A0 (en) | 2001-12-21 | 2001-12-21 | Method of processing marine magnetic gradiant data and exploration methods using that data |
| PCT/AU2002/001669 WO2003052460A1 (en) | 2001-12-18 | 2002-12-10 | Method of processing marine magnetic gradient data and exploration methods using that data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| AU2002366366A1 true AU2002366366A1 (en) | 2003-06-30 |
Family
ID=25646861
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| AU2002366366A Abandoned AU2002366366A1 (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 (en) |
| AU (1) | AU2002366366A1 (en) |
| GB (1) | GB2417086B (en) |
| NO (1) | NO20052790L (en) |
| WO (1) | WO2003052460A1 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114415073A (en) * | 2022-03-29 | 2022-04-29 | 中国人民解放军海军工程大学 | Ground quick calibration method and system for aeromagnetic vector gradiometer error model |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100434934C (en) * | 2006-07-12 | 2008-11-19 | 杨辉 | Optimization processing technology for heavy magnetism by using continuation returning and vertical derivation technology |
| DE102006059204B4 (en) * | 2006-12-13 | 2008-10-30 | Bundesrepublik Deutschland, vertreten durch den Präsidenten der Bundesanstalt für Geowissenschaften und Rohstoffe | Method and installation for measuring the geomagnetic field on the open sea |
| WO2009062236A1 (en) * | 2007-11-12 | 2009-05-22 | Commonwealth Scientific And Industrial Research Organisation | Method and apparatus for detecting marine deposits |
| CN102338890B (en) * | 2010-10-22 | 2013-04-24 | 中国石油天然气股份有限公司 | Round window band-pass amplitude preservation filtering data processing method in geophysical exploration |
| WO2013138908A1 (en) * | 2012-03-23 | 2013-09-26 | Fugro Canada Corp. | System and method for geophysical surveying using electromagnetic fields and gradients |
| CN104297336B (en) * | 2014-10-22 | 2017-02-08 | 中国地质大学（武汉） | Buried-steel-pipeline-based magnetic anomaly extraction and interpretation method |
| CN109839877B (en) * | 2019-03-03 | 2021-07-27 | 上海卯瑞船舶设备有限公司 | Monitoring device for ship remote management and use method thereof |
| CN114814963B (en) * | 2022-03-03 | 2023-03-21 | 吉林大学 | B-spline channel extraction-based time domain electromagnetic data very low frequency noise suppression method |
Family Cites Families (3)
| 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 |
| FR2672395B1 (en) * | 1991-01-31 | 1993-04-30 | Sextant Avionique | METHOD AND DEVICE FOR REDUCING THE EFFECTS OF INTERFERENCE NOISE ON THE DETECTION OF A TARGET BY A SYSTEM COMPRISING A PLURALITY OF ELEMENTAL SENSORS. |
| CN2439035Y (en) * | 2000-09-18 | 2001-07-11 | 中国地质大学（北京） | Combined circuits integration for collecting subsea and earth electromagnetic signals |
-
2002
- 2002-12-10 GB GB0512665A patent/GB2417086B/en not_active Expired - Fee Related
- 2002-12-10 WO PCT/AU2002/001669 patent/WO2003052460A1/en not_active Application Discontinuation
- 2002-12-10 CN CNB028301544A patent/CN100399052C/en not_active Expired - Fee Related
- 2002-12-10 AU AU2002366366A patent/AU2002366366A1/en not_active Abandoned
-
2005
- 2005-06-09 NO NO20052790A patent/NO20052790L/en not_active Application Discontinuation
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114415073A (en) * | 2022-03-29 | 2022-04-29 | 中国人民解放军海军工程大学 | Ground quick calibration method and system for aeromagnetic vector gradiometer error model |
| CN114415073B (en) * | 2022-03-29 | 2022-08-09 | 中国人民解放军海军工程大学 | Ground quick calibration method and system for aeromagnetic vector gradiometer error model |
Also Published As
| Publication number | Publication date |
|---|---|
| GB0512665D0 (en) | 2005-07-27 |
| GB2417086A (en) | 2006-02-15 |
| NO20052790L (en) | 2005-09-05 |
| NO20052790D0 (en) | 2005-06-09 |
| CN100399052C (en) | 2008-07-02 |
| CN1720465A (en) | 2006-01-11 |
| WO2003052460A1 (en) | 2003-06-26 |
| GB2417086B (en) | 2006-08-16 |
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 | |
| US5724309A (en) | Method for geophysical processing and interpretation using instantaneous phase and its derivatives and their derivatives | |
| Fedi et al. | A method to estimate the total magnetization direction from a distortion analysis of magnetic anomalies 1 | |
| US20100094556A1 (en) | Terrain correction systems | |
| GB2497226A (en) | Method for estimating elastic parameters by inverting 4D seismic measurements | |
| AU2002366366A1 (en) | Method of processing marine magnetic gradient data and exploration methods using that data | |
| CN108983300B (en) | Transient electromagnetic tunnel advanced prediction method under tunnel boring machine construction condition | |
| Nelson | Leveling total-field aeromagnetic data with measured horizontal gradients | |
| CA2110433C (en) | Magnetic cartography process and apparatus | |
| US7313495B2 (en) | Method of processing marine magnetic gradient data and exploration methods using that data | |
| CN103543761B (en) | Control the method and system of the hauling speed of sensor towing cable | |
| EP0966693A1 (en) | A method of and apparatus for determining the quality of seismic data | |
| Eichkitz et al. | Enhanced coherence attribute imaging by structurally oriented filtering | |
| Reford | Gradient enhancement of the total magnetic field | |
| CN107038700A (en) | Soil moisture spatial distribution characteristic quantitative analysis method based on GSSIM | |
| Nielsen et al. | Integrating ground-penetrating radar and borehole data from a Wadden Sea barrier island | |
| RU2298815C2 (en) | Method of treatment of sea magnetic gradient data and methods of exploration on base of the data | |
| US7830746B2 (en) | Method for treating seismic data corresponding to acquisitions obtained for a common zone by means of underwater seismic receivers and by means of surface receivers | |
| US20070055465A1 (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 | |
| CN113552632A (en) | Seismic cross-correlation signal pickup method and system based on wavelet domain Kalman filtering | |
| Collet et al. | Spot satellite data analysis for bathymetric mapping | |
| Shon et al. | Simple data processing procedures for seismic section noise reduction | |
| CN112925020B (en) | Parametric array type shallow profile intelligent energy compensation method | |
| O’Connell | A heuristic method of removing micro-pulsations from airborne magnetic data | |
| Nobes* et al. | Enhancing form and structure: Complex attributes as aids for ground penetrating radar interpretation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MK4 | Application lapsed section 142(2)(d) - no continuation fee paid for the application |