AU2016102088A4 - Method and system for measuring total earth’s magnetic field for ocean magnetic survey through elimination of geomagnetic disturbance and recording medium therefor - Google Patents
Method and system for measuring total earth’s magnetic field for ocean magnetic survey through elimination of geomagnetic disturbance and recording medium therefor Download PDFInfo
- Publication number
- AU2016102088A4 AU2016102088A4 AU2016102088A AU2016102088A AU2016102088A4 AU 2016102088 A4 AU2016102088 A4 AU 2016102088A4 AU 2016102088 A AU2016102088 A AU 2016102088A AU 2016102088 A AU2016102088 A AU 2016102088A AU 2016102088 A4 AU2016102088 A4 AU 2016102088A4
- Authority
- AU
- Australia
- Prior art keywords
- probe
- magnetic field
- earth
- magnetic
- measured
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired
Links
- 230000005291 magnetic effect Effects 0.000 title claims abstract description 176
- 238000000034 method Methods 0.000 title claims description 42
- 230000008030 elimination Effects 0.000 title description 4
- 238000003379 elimination reaction Methods 0.000 title description 4
- 239000000523 sample Substances 0.000 claims abstract description 86
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 239000011159 matrix material Substances 0.000 claims description 38
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 238000000691 measurement method Methods 0.000 abstract 2
- 238000012937 correction Methods 0.000 description 8
- 239000002184 metal Substances 0.000 description 7
- 229910052751 metal Inorganic materials 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000005389 magnetism Effects 0.000 description 2
- 239000007769 metal material Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 210000003462 vein Anatomy 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 229910001338 liquidmetal Inorganic materials 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 150000002739 metals Chemical class 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
- 230000005386 paleomagnetism Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Provided are ocean magnetic survey measurement method and system. The ocean magnetic survey measurement method includes a measuring operation for 5 measuring an earth's magnetic field value using a magnetometer according to a heading angle of a probe, which is measured in clockwise direction on a basis of the true north, while the probe moves along a circular path based on a center point in an exploration target area with the magnetometer taken in tow, a data input operation for receiving, by a calculation unit, data for the heading angle of the probe and data 10 for the earth's magnetic field value measured at the heading angle, a calculating operation for setting, by the calculation unit containing a calculation program therein, a magnetic disturbance value induced by the probe, which is included in the measured earth's magnetic field value, and calculating an actual earth's magnetic field value by subtracting the magnetic disturbance value from the measured earth's 15 magnetic field value.
Description
METHOD AND SYSTEM FOR MEASURING TOTAL EARTH’S MAGNETIC FIELD FOR OCEAN MAGNETIC SURVEY THROUGH ELIMINATION OF GEOMAGNETIC DISTURBANCE AND RECORDING
MEDIUM THEREFOR
BACKGROUND
[0001] The present disclosure relates to an ocean magnetic survey for exploring metal natural resources, large-scale geological structures, and pipelines in the sea floor, and more particularly, to a method and system for accurately obtaining a total earth’s magnetic field value by correcting distortion, caused by a probe constructed with a metal, in value of magnetic survey measurement at the time of undertaking the ocean magnetic survey, and a computer-readable recording medium with which the method is executable in a computer.
[0002] FIG. 1 is a schematic diagram of a view that an ocean magnetic survey is undertaken.
[0003] Referring to FIG. 1, typically, an ocean magnetic survey is a method for measuring a total earth’s magnetic field with a marine magnetometer taken in tow at the stem of a probe.
[0004] The earth’s magnetic field is established with about 98% of internal factors caused by liquid metal movement(referring to FIG.2) in the earth core and remaining about 2% of external factors by the sun. Unlike the earth’s gravitational field, the earth's magnetic field continuously varies and such a variation includes a diurnal variation of every day, a secular variation of several hundred or thousand years, and a sudden variation of magnetic storm.
[0005] As illustrated in FIG. 3, the earth’s magnetic field is represented with a vector having the magnitude and the direction, which is represented with three elements of earth’s magnetic field including the total magnetic field F, which is the magnitude of the earth’ s magnetic field, a declination D, which is an angle between a horizontal component H of the earth’s magnetic field and the true north, and an inclination I that is an angle between H and F. A unit of the earth’s magnetic field is nano-Tesla (nT), and values thereof at the equator and poles are respectively about 30,000 nT and 65,000 nT.
[0006] [Table 1] Application fields of magnetic survey (use of geophysical exploration, Jisu Kim et. al)
Positioning
Pipes, cables, and metal objects
Embedded military supplies (shells, bombs, etc.)
Buried metal drums having contaminated or harmful wastes therein Shaft or horizontal adit of collapsed mine Mapping
Archaeological remains Hidden dyke of igneous rock Vein containing metals
Magnetically contrasted geological boundary (stratum etc.) of rock floor Large-scale geological structures [0007] As shown in Table 1, the magnetic survey has a wide application range from small-scale exploration for searching a pipeline or cable near the surface of the earth to large-scale exploration such as exploration of petroleum resources. As the magnetic survey is frequently used, researches on magnetism of rock or mineral become active. In particular, researches on and applications of residual magnetism of a rock are developed to paleomagnetism, which occupies an important position in geophysics since 1950s, and to be presented as a quantitative evidence for proving main theories of modem geology such as a continent drift theory or a seafloor spreading theory. In addition, in view of industry, they are particularly used for investigating a subsurface structure such as bed rock at the time of exploring petroleum resources and are usefully used together with gravity survey and elastic wave exploration data.
[0008] In other words, while the total earth’s magnetic field of an exploration target area is measured through the ocean magnetic survey, sudden variation of the earth’s magnetic field is determined and elements such as metal cables or pipelines, or metal veins, which influence on the magnetic field, may be explored. Accordingly, it is very important to accurately measure a value of total earth’s magnetic field in the ocean magnetic survey. However, since the value measured with marine magnetometer includes a value of magnetic field induced by a probe (manufactured from a ferromagnetic steel material), it is difficult to obtain a pure value of the total earth’s magnetic field.
SUMMARY
[0009] The present disclosure provides a method, system, and computer recording medium for measuring a total earth’s magnetic field for an ocean magnetic survey through elimination of geomagnetic disturbance, which is capable of accurately measuring the total earth’s magnetic field in an exploration target area by eliminating influence of a magnetic field induced by a probe manufactured from a metal material at the time of ocean magnetic survey.
[0010] In accordance with an exemplary embodiment of the present invention, a method for measuring a total earth's magnetic field for an ocean magnetic survey, including : a measuring operation for measuring an earth’s magnetic field value using a magnetometer according to a heading angle of a probe, which is measured in clockwise direction on a basis of the true north, while the probe moves along a circular path based on a center point in an exploration target area with the magnetometer taken in tow; a data input operation for receiving, by a calculation unit, data for the heading angle of the probe and data for the earth’s magnetic field value measured at the heading angle; a calculating operation for setting, by the calculation unit containing a calculation program therein, a magnetic disturbance value induced by the probe, which is included in the measured earth’s magnetic field value ; and calculating an actual earth's magnetic field value by subtracting the magnetic disturbance value from the measured earth’s magnetic field value.
[0011] According to an embodiment, in the calculating operation, the magnetic disturbance value may be calculated by the following Equation (1),
•••(I) where Θ denotes the heading angle of the probe, which is measured from the true north in clockwise direction, Co, Ci, and C2, denote constants according to magnetic characteristics of the probe and distances from the exploration target area and the magnetometer to the probe, and the calculation program may determine the constants of Equation (1).
[0012] According to another embodiment, the calculation program may determine the constants Co, Ci, and C2 of Equation (1) of the magnetic disturbance value so as to minimize differences between variations of the earth’s magnetic field values continuously measured while the probe moves and changes the heading angle on a basis of an earth’s magnetic field value measured at an initial point (the true north direction of the heading angle of 0° of the probe) and each magnetic disturbance value induced by the probe which moves and changes the heading angle.
BRIEF DESCRIPTION OF THE DRAWINGS
[0013] Exemplary embodiments can be understood in more detail from the following description taken in conjunction with the accompanying drawings, in which: [0014] FIG. 1 is a schematic diagram for explaining an ocean magnetic survey; [0015] FIG. 2 is a view for explaining constitution of the earth’s magnetic field; [0016] FIG. 3 is a view for explaining three elements of the earth’s magnetic field; [0017] FIG. 4 is a graph representing a total earth’s magnetic field measured according to the heading of a probe and a best fit curve of the total earth’s magnetic field; [0018] FIG. 5 is a schematic flowchart of a method for measuring the total earth’s magnetic field according to the present disclosure; [0019] FIG. 6 is a view for explaining a circular survey; [0020] FIG. 7 is a graph representing a concept of the least squares method; [0021] FIG. 8 is a schematic diagram representing that a Joides Resolution probe vessel undertakes an ocean survey, and [0022] FIG. 9 is a graph representing magnetic field values measured by a magnetometer in an experiment, and a best fit curve formed by using the magnetic field values and an operation program.
DETAILED DESCRIPTION OF EMBODIMENTS
[0023] As described above, it is difficult to determine an accurate value of the total earth’s magnetic field in the exploring target area due to an influence of the magnetic field induced by the probe in the ocean magnetic survey. Thus in 1961, Bullard and Mason proposed a method (i.e. the following Equation (1)) for eliminating an influence of a probe from magnetic measurement data, and the present disclosure is based on the Bullard and Mason’s method.
......(1) where Fq denotes a value measured by a marine magnetometer, F denotes a value of a total earth’s magnetic field, Θ denotes the heading of the probe, which is measured from the north in a clockwise direction, and Co, Ci, C2, Si, and S2 denote constants according to magnetic characteristics of the probe, and distances from an exploration target area and the magnetometer to the probe.
[0024] Here, for a probe having a bilateral symmetry, an influence of sine values is very smaller than that of cosine values, and Si and S2 may be set as 0 (Bullard and Mason, 1961). Accordingly, Equation (1) may be arranged to the following Equation (2).
.......(2) [0025] The underlined part of Equation (2) is caused by a magnetic field induced by the probe, and is an element for distorting a measured value Fq of the earth’s magnetic field. Accordingly, once values of Co, Ci, and C2, are known, a pure total earth’s magnetic field F in an exploring target area may be obtained by eliminating the underlined elements from the measured magnetic field Fq.
[0026] ‘Bullard and Mason’ proposed a method for eliminating the probe’s influence.
[0027] FIG. 4 is a graph representing a total earth’s magnetic field measured according to the heading of the probe and a best fit curve of the total earth’s magnetic field.
[0028] Referring to FIG. 4, the heading direction (a horizontal axis) of the probe is represented as 0 to 360°. In other words, the probe moves along a path of a circle drawn from a specific point and returned thereto. In addition, the earth’s magnetic values are measured at each point (a heading angle of the probe) on the path and are represented as points on a vertical axis of the graph. In addition, ‘Bullard and Mason’ suggested that the constants Co, Ci, and C2 may be determined by obtaining a best fit curve to the measured values. However, since ‘Bullard and Mason’ did not propose a method for obtaining the best fit curve to the measured values, there is a limitation for using the method.
[0029] In the present disclosure, with respect to the measured values, a matrix and a least squares method are introduced to obtain the constants Co, Ci, and C2, and through this, the total earth’s magnetic field may be accurately measured.
[0030] Hereinafter, a description will be provided in detail about a method for measuring a total earth’s magnetic field in an ocean magnetic survey through eliminating geomagnetic disturbance (hereinafter, total earth’s magnetic field measuring method, total earth’s magnetic field measuring system) according to an embodiment of the present disclosure. The total earth’s magnetic field measuring method according to the present disclosure may be realized with a time series method, and may be more concretely realized by a computer readable program. Accordingly, the present disclosure is performed by executing a program in a computer in which the program is installed or in a computer using a recording medium on which the program is recorded.
[0031] Firstly, the total earth’s magnetic field measuring system according to the present disclosure includes a probe, a marine magnetometer, and a calculation unit with a computer program contained therein. As described above, the probe is manufactured from a metal material, and the marine magnetometer is disposed in the sea in a state of being connected to the probe through a wire line. According to movement of the probe, the marine magnetometer is taken in tow and measures a magnetic field value of each point. In addition, magnetic field data obtained by the marine magnetometer is transmitted to the calculation unit by an electrical connection means, and the calculation means calculates a pure value of total earth’s magnetic field in an exploration target area by removing an earth’s magnetic field distortion element caused by the probe through a calculation algorithm described later.
[0032] Using the foregoing system, the total earth’s magnetic field measuring method will be described.
[0033] FIG. 5 is a schematic flowchart of the total earth’s measuring method according to the present disclosure, and FIG. 6 is a view for explaining a circular survey.
[0034] Referring to FIGS. 5 and 6, the total earth’s measuring method according to the present disclosure firstly starts from measuring an earth’s magnetic field value at each point (or at each heading angle of the probe), while the probe draws a circle in the exploration target area. Regardless of a start direction, the probe draws a circle and measures the magnetic value. The measured data is included in the wire line or is transmitted to the calculation unit through a wired or wireless communication line separated from the wire line.
[0035] For the measured data, the calculation unit calculates the constants Co, Ci, and C2 for determining the influence of the magnetic field induced by the probe by applying a calculation algorithm according to the present disclosure.
[0036] In addition, the pure value F of the total earth’s magnetic field in the exploration area may be obtained by removing the influence by the probe from the total earth’s magnetic field value Fq that is obtained by substituting the constants to Equation (2).
[0037] Hereinafter, the core equation of the calculation program, which is developed by researchers of the present disclosure, is expressed as the following Equation (3) by subtracting a matrix Y from multiplication of matrixes A and C and introducing the least squares method (parts expressed with min and square) and 2-Norm (parts expressed with a subscript of 2 and double absolute value bars).
...........(3) [0038] From Equation (3), a matrix C for the magnetic disturbance constants Co, Ci, and C2 may be obtained like the following Equation (4) (detailed procedure will be described later).
..............(4) where a superscript T in Equation (4) denotes a transpose matrix of each matrix. For example, A is a transpose matrix of a matrix A.
[0039] Firstly, concepts of the least squares method and 2-Norm introduced to the calculation program are well-known in Algebra, and therefore they are briefly explained herein. Referring to the graph of FIG. 7, the measured values (red points) and a best fit curve (blue line) to the measured values are expressed with a function y. In addition, errors between actually measured values and function values are represented. The best fit curve to the measured values may be obtained by using the least squares method. Additionally, the researchers combine a mathematical tool of Norm to the least squares method. Norm is for efficiently analyzing an error in numerical solution of an equation.
[0040] In other words, when Χ=(χι, X2, X3 ···· xn), 1-Norm and 2-Norm are defined as the following.
[0041] 1-Norm is a sum of absolute values of each value, and 2-Norm is to put a sum of squares of each value to the power of 1/2. The present calculation program uses 2-Norm.
[0042] Equation (9), which is the core of the calculation program employed herein, is derived by defining a difference between the measured magnetic field value and an actual magnetic field value, namely, the value of magnetic field induced by the probe (herein a magnetic disturbance value). In other words, through applying the binominal theorem to Equation (2), the pure value F of magnetic field in the exploration area is subtracted from the actually measured magnetic value Fq to be expressed as the following.
[0043] In addition, since values are generated according to the heading direction Θ of the probe, a matrix is introduced to express the values as the following.
a * ε = γ (here, when n<3, where n is the number of data measured through the circular survey, C may not exist) [0044] In other words, multiplication A»C of matrixes A and C is a combined expression of a magnetic disturbance value at each point of the circular survey.
[0045] In addition, when a total earth's magnetic field value measured at 0° (start point) of the heading direction of the probe is set as a reference, Y may indicate a variation of the total earth’s magnetic field value represented according to the probe’s heading direction. In other words, since the probe does not move when the heading thereof is at 0°, the actually measured value Fq by the magnetometer may be considered as a pure value that is not influenced by the magnetic disturbance value caused by the probe. In addition, as the probe changes the heading thereof and moves, a variation value Fq-Fo of the actually measured earth’s magnetic field may be considered as the magnetic disturbance value yn occurred by the probe.
[0046] Accordingly, a matrix function AC-Y for which the matrix Y is subtracted from the multiplication AC of the matrixes A and C is set. For the value of AC-Y, the least squares method and 2-Norm are introduced to set [minllAC-YII 2 ] like Equation (3), and the matrix C for minimizing the value of Equation (3) is obtained.
[0047] Obtaining a solution of the matrix C for minimizing Equation (3) means that the best fit curve to the actually measured values may be obtained.
[0048] Hereinafter, a procedure for obtaining the solution of matrix C will be described.
[0049] From the matrix of Equation (3), the 2-Norm and squares may be expressed as the following. minllAC-YII22 = min(AC-Y)T(AC-Y) = min [(AC)T(AC)-(AC)TY-YT(AC)+YTY] = min (CTATAC-2CTATY+YTY)... Intermediate equation (1) [0050] In the intermediate equation (1), the superscript T denotes a transpose matrix.
[0051] Accordingly, in the intermediate equation (1), CTATAC-2CTATY+YTY may be changed to C2A2-2CAY+Y2. This is because the transpose matrix is just to transpose rows and columns of a matrix and an absolute value thereof is identical.
[0052] In order to obtain the matrix C, the following intermediate equation (2) is derived from the intermediate equation (1) through partial differentiation.
......intermediate equation (2) [0053] In other words, partially differentiating C2A2-2CAY+Y2 with respect λ to C becomes 2CA -2AY. However, since a sequence change is not allowed in the matrix operation, the intermediate equation (2) is arranged in that order.
[0054] Accordingly, the following Equation (4) may be derived by solving the intermediate equation (2) of 2A AC-2A Y=0 with respect to C. Equation (4) is like the following.
..........(4) -1 τ where (ΑΤΑ)" is an inverse matrix of A A, Y denotes a variation A of the total magnetic field measured along the heading of the probe, when the total magnetic field value measured at 0° of the probe’s heading is set as a reference.
[0055] Accordingly, by solving Equation (4), the solution (Co, Ci, C2) of the matrix C may be obtained. The calculation unit, which is an element of the present disclosure, contains Equation (4), and when data is input, calculation for obtaining the solution of matrix C is performed.
[0056] When the solution of the matrix is obtained, since the magnetic disturbance value (Co + Qcos + C2COS2 ) induced by the probe may be known, the pure value of total magnetic value in the exploration target area may be calculated by subtracting the magnetic disturbance value from the actually measured total magnetic field value at each heading of the probe.
[0057] In other words, in the present disclosure, it is set (AC-Y=0) that there is no influence of disturbance magnetic field induced by the probe by selecting the total magnetic field value measured at a initial start time ( a heading angle 0° of the probe) of the probe as an actual magnetic value. Thereafter, a solution for the matrix C is obtained for minimizing a difference between total magnetic field values Y (Y values actually mean a variation from the start time, since it is set that AC-Y=0 at the start time) measured while the probe moves. Once the solution of the matrix C is obtained, since the disturbance magnetic field value induced by the probe is capable of being calculated, the actual total magnetic field value may be accurately calculated by subtracting the magnetic disturbance value from the actually measured total magnetic value.
[0058] Hereinafter, an example of a procedure for mathematically obtaining the matrix C of Equation (4) will be described. In order to avoid mathematical complexity, the procedure is simplified with values measured at four angles i, 2, 3» and 4. The operation is performed as the following four steps of STEP 1 to 4. STEP1: ATA (here, AT is a transpose matrix of A)
Aij = value that horizontal components of AT are respectively multiplied by vertical components of A and summed STEP2: obtain (ATA)'* (i.e. inverse matrix of ATA)
When
, (ATA)_1 may be expressed as the following by the inverse matrix theorem.
STEP3: obtain ATY
STEP4: substitute (ATA)_1 ·ΑΤΥ
At this point, matrix and matrix
[0059] As described above, the calculation unit may receive data for the heading angle Θ of the probe and data for magnetic values measured at each angle, solve the above-described matrixes to obtain the magnetic disturbance constants C(Co, Ci, C2) and may obtain the magnetic disturbance value induced by the probe. In addition, an actual magnetic field value in the exploration target area may be accurately determined by subtracting the magnetic disturbance value from the actually measured magnetic value.
[0060] The researchers of the present disclosure performed experiments in order to investigate effects of the total magnetic field measuring method and calculation program. FIG. 8 is a schematic diagram representing that a Joides Resolution probe vessel undertakes an ocean survey.
[0061] As illustrated in FIGS. 8, a circular survey was performed in an area of “Ori Massif’ by using a probe vessel of “Joides Resolution”. In the experiment, the magnetic field value measured by the magnetometer and a best fit curve formed by the present calculation program were represented in FIG. 9.
[0062] Through the above-described experiment, a magnetic survey was performed by using the vessel of “Joides Resolution” in area of “Ori Massif’ and magnetic disturbance constants C(Co, Ci, C2) were obtained as the following by using the calculation program, and the disturbance magnetic field value Fh by the probe vessel was obtained.
[0063] Fh = Co + Cicos +C2COS2 = 40.91 - 45.89cos +3.66cos2 [0064] Referring to FIG. 9, as a result of deriving a best fit curve by determining a constant value C with the calculation program, it is shown that the actually measured value (points on the graph) and the best fit curve are almost matched.
[0065] The following Table 2 and Table 3 show results of cross-over error comparison before and after correction of the magnetic disturbance caused by the probe with the present calculation program in the above-described experiment.
[0066] [Table 2] Cross-over error before correction of magnetic disturbance «coord ycoord XOE Ref Line RefCnt RefVai CmpLine CmpCnt CmpVai WAveVai 158 804*778 36.1067347 6« EXP324-L5T 87 41730 8 EXP324-L5T 159 41724 1 41727 5 158.8058571 36.1085933 5.7 EXP324-L5T 88 41731.2 EXP324-L5T 160 41725.5 41728.4 158.8825014 36.0945155 48.2 EXP324-L5T 117 41817.7 EXP324-L5T 182 41769.5 41733.6
Before Correct: Num ofXOPs = 3, Mean XOE = 20-2, RMS XOE == 20.2, Max XOE = 48.2, Min XOE = 5.7 [0067] Before the correction of magnetic disturbance, it shows an average error of 20.2 nT, a maximum error of 48.2 nT, and a minimum error of 5.7 nT.
[0068] [Table 3] Cross-over error after correction of magnetic disturbance xcoord ycoord XOE RefLine RefCnt ReiVa! CmpLine CmpCnt CmpVai WAveVai 158 8044778 381067347 6.5 EXP324-L5T 87 41726.9 EXP324-LST 159 41720A 41722.7 158.8058571 361085933 121 EXP324-L5T 88 41728.0 EXP324-LST 160 41715.9 41721.9 158.8825014 36 0945155 -0.1 EXP324-LST 117 41727.5 EXP324-L5T 182 41727.5 41727.5
Before Correct : Num of XOPs - 3, Mean XOE - 6.2, RMS XOE - 6,2, Max XOE - 12.1, Min XOE - -Q.1 [0069] After the correction of magnetic disturbance, it shows an average error of 6.2 nT, a maximum error of 12.1 nT, and a minimum error of -0.1 nT. In other words, as a result of performing the correction with the present calculation program, it may be seen that an error range is remarkably reduced. This means that when the circular survey is performed and the calculation program is executed according to the present disclosure, an influence of the disturbance magnetic field induced by the probe is minimized in an ocean magnetic survey to enable a reliable magnetic survey.
[0070] In particular, the present disclosure is not applied only to a specific magnetometer, but is useable to all magnetometers regardless of types thereof. Furthermore, even an area where the exploration is already performed, when there are magnetic field measurement values according to a heading angle of the probe, correction therefor may be possible through post-processing data by executing the present calculation program.
[0071] According to the method and system for measuring a total earth’s magnetic field in the ocean magnetic survey, when a magnetic survey is undertaken by using a vessel in an exploration target area, the value of total earth’s magnetic field in the exploration target area may be accurately measured by removing an influence of the magnetic field induced by the vessel.
[0072] In addition, post-processing correction for the magnetic survey data is possible by using already measured data, even when the exploration is completed.
[0073] In addition, the system and method according to the present disclosure may be applied to all the magnetometers used for magnetic survey regardless of a type thereof.
[0074] Although the method and system for measuring total earth’s magnetic field in ocean magnetic survey through elimination of geomagnetic disturbance and Recording medium therefor have been described with reference to the specific embodiments, they are not limited thereto. Therefore, it will be readily understood by those skilled in the art that various modifications and changes can be made thereto without departing from the spirit and scope of the present invention defined by the appended claims.
Claims (6)
- WHAT IS CLAIMED IS:1. A method for measuring a total earth's magnetic field for an ocean magnetic survey, the method comprising: a measuring operation for measuring an earth’s magnetic field value using a magnetometer according to a heading angle of a probe, which is measured in clockwise direction on a basis of the true north, while the probe moves along a circular path based on a center point in an exploration target area with the magnetometer taken in tow; a data input operation for receiving, by a calculation unit, data for the heading angle of the probe and data for the earth’s magnetic field value measured at the heading angle; a calculating operation for setting, by the calculation unit containing a calculation program therein, a magnetic disturbance value induced by the probe, which is comprised in the measured earth’s magnetic field value ; and calculating an actual earth's magnetic field value by subtracting the magnetic disturbance value from the measured earth’s magnetic field value.
- 2. The method of claim 1, wherein in the calculating operation, the magnetic disturbance value is calculated by the following Equation (1),•••(I) where Θ denotes the heading angle of the probe, which is measured from the true north in clockwise direction, Co, Ci, and C2, denote constants according to magnetic characteristics of the probe and distances from the exploration target area and the magnetometer to the probe, and the calculation program determines the constants of Equation (1).
- 3. The method of claim 2, wherein the calculation program determines the constants Co, Ci, and C2 of Equation (1) so as to minimize differences between variations of the earth’s magnetic field values continuously measured while the probe moves and changes the heading angle on a basis of an earth’s magnetic field value measured at an initial point (the true north direction of the heading angle of 0° of the probe) and each magnetic disturbance value induced by the probe which moves and changes the heading angle.
- 4. The method of claim 2, wherein the calculation program determines the constants Co, Ci, and C2 by calculating the following matrix of Equation (2),... (2) where a matrix A, expressed as the following, denotes a matrix having the heading angle of the probe as a variable, a matrix C denotes a matrix for the constant values of Equation (1), Y denotes variations of the earth’s magnetic field values continuously measured while the probe moves and changes the heading angle on a basis of the earth’s magnetic field value actually measured by the probe or an earth’s magnetic field value measured at an initial time (the heading angle of 0° of the probe ), a superscript T denotes a transpose matrix, and a superscript -1 denotes an inverse matrix.(here, when n<3, where n is the number of data measured through the circular survey, C may not exist)
- 5. The method of claim 4, where Equation (2) is derived from the following Equation (3) which introduces mathematical tools of “least squares method" and “2-Norm”.-(3)
- 6. A system for measuring a total earth's magnetic field for an ocean magnetic survey, the system comprising: a probe operated in water; a magnetometer which is connected to the probe through a wire line and is taken in tow by the probe and which measures an earth’s magnetic field in an exploration target area; and a calculation unit receiving, from the magnetometer, data for the measured earth’s magnetic field and a heading angle of the probe and containing therein a computer readable calculation program for calculating a magnetic disturbance value induced by the probe, wherein the calculation program performs the calculating operation in any one of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| AU2016102088A AU2016102088A4 (en) | 2015-05-20 | 2016-12-07 | Method and system for measuring total earth’s magnetic field for ocean magnetic survey through elimination of geomagnetic disturbance and recording medium therefor |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR10-2015-0070256 | 2015-05-20 | ||
| PCT/KR2015/013641 WO2016186283A1 (en) | 2015-05-20 | 2015-12-14 | Method and system for measuring earth's magnetic field via marine magnetic force exploration for eliminating magnetic field disturbance, and recording medium therefor |
| AU2016102088A AU2016102088A4 (en) | 2015-05-20 | 2016-12-07 | Method and system for measuring total earth’s magnetic field for ocean magnetic survey through elimination of geomagnetic disturbance and recording medium therefor |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/KR2015/013641 Division WO2016186283A1 (en) | 2015-05-20 | 2015-12-14 | Method and system for measuring earth's magnetic field via marine magnetic force exploration for eliminating magnetic field disturbance, and recording medium therefor |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| AU2016102088A4 true AU2016102088A4 (en) | 2017-01-12 |
Family
ID=57724821
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| AU2016102088A Expired AU2016102088A4 (en) | 2015-05-20 | 2016-12-07 | Method and system for measuring total earth’s magnetic field for ocean magnetic survey through elimination of geomagnetic disturbance and recording medium therefor |
Country Status (1)
| Country | Link |
|---|---|
| AU (1) | AU2016102088A4 (en) |
-
2016
- 2016-12-07 AU AU2016102088A patent/AU2016102088A4/en not_active Expired
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20160341846A1 (en) | Method and system for measuring total earth's magnetic field for ocean magnetic survey through elimination of geomagnetic disturbance and recording medium therefor | |
| Munschy et al. | Magnetic mapping for the detection and characterization of UXO: Use of multi-sensor fluxgate 3-axis magnetometers and methods of interpretation | |
| Fedi | DEXP: A fast method to determine the depth and the structural index of potential fields sources | |
| Lelièvre et al. | A 3D total magnetization inversion applicable when significant, complicated remanence is present | |
| Hinze et al. | Gravity and magnetic exploration: Principles, practices, and applications | |
| Schmidt et al. | The magnetic gradient tensor: Its properties and uses in source characterization | |
| Lelièvre et al. | Joint inversion of seismic traveltimes and gravity data on unstructured grids with application to mineral exploration | |
| Pilkington et al. | Mitigating remanent magnetization effects in magnetic data using the normalized source strength | |
| Lu et al. | Paleomagnetism and anisotropy of magnetic susceptibility of the Tertiary Janggalsay section (southeast Tarim basin): Implications for Miocene tectonic evolution of the Altyn Tagh Range | |
| Kabirzadeh et al. | Dynamic error analysis of measurement while drilling using variable geomagnetic in-field referencing | |
| Liu et al. | Distortion magnetic field compensation of geomagnetic vector measurement system using a 3-D Helmholtz coil | |
| Liu et al. | Extracting induced and remanent magnetizations from magnetic data modeling | |
| CN105388533A (en) | Onshore correction method of magnetic interference of magnetometer installed in submarine vehicle | |
| Deng et al. | Suitability of magnetometry to detect clandestine buried firearms from a controlled field site and numerical modeling | |
| Virgil et al. | Reorientation of three‐component borehole magnetic data | |
| Schiffler et al. | Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data | |
| Liu et al. | A new method for distortion magnetic field compensation of a geomagnetic vector measurement system | |
| AU2016102088A4 (en) | Method and system for measuring total earth’s magnetic field for ocean magnetic survey through elimination of geomagnetic disturbance and recording medium therefor | |
| Pan et al. | Magnetic anomaly inversion using magnetic dipole reconstruction based on the pipeline section segmentation method | |
| Beggan et al. | Estimating external magnetic field differences at high geomagnetic latitudes from a single station | |
| Xu et al. | Simulation Analysis of Magnetic Gradient Full‐Tensor Measurement System | |
| Ribeiro | Recovering the total magnetization direction using a novel stable reduction to the pole filter (RTP) | |
| Vaish et al. | Interpretation of magnetic anomaly data over east Basuria region using an enhanced local wavenumber (ELW) technique | |
| Narkhov et al. | The perspective of full-gradient magnetometry system application in searching for magnetic objects | |
| Flanagan et al. | Improvements in magnetic depth estimation: application of depth and width extent nomographs to standard depth estimation techniques |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FGI | Letters patent sealed or granted (innovation patent) | ||
| MK22 | Patent ceased section 143a(d), or expired - non payment of renewal fee or expiry |