CN103926626A - Height correcting method and system for aeromagnetic data - Google Patents
Height correcting method and system for aeromagnetic data Download PDFInfo
- Publication number
- CN103926626A CN103926626A CN201410182126.5A CN201410182126A CN103926626A CN 103926626 A CN103926626 A CN 103926626A CN 201410182126 A CN201410182126 A CN 201410182126A CN 103926626 A CN103926626 A CN 103926626A
- Authority
- CN
- China
- Prior art keywords
- height
- aeromagnetic data
- data
- aeromagnetic
- corrected
- 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.)
- Pending
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a height correcting method and system for aeromagnetic data. The height correcting method for the aeromagnetic data includes the steps of obtaining the flight height h of an aircraft and the aeromagnetic data T(h) measured when the aircraft is at the flight height h, presetting a to-be-corrected height h0 of the aeromagnetic data, and generating the corrected aeromagnetic data according to the flight height h, the aeromagnetic data T(h) and the to-be-corrected height h0. By the height correcting method and system for the aeromagnetic data, the stripped abnormal phenomenon caused by difference of intensities of the geomagnetic field can be avoided, false abnormities caused by the different flight heights can be eliminated, and a foundation is provided for further processing the aeromagnetic data.
Description
Technical field
The invention relates to geophysical exploration technology, especially about magnetic airborne surveys technology, particularly about a kind of Aeromagnetic data height antidote and system.
Background technology
Magnetic airborne surveys is one of important method of geophysical survey, is widely used in the aspects such as petroleum prospecting, ore exploration, Coal Exploration.Particularly in waters, the area of the topographic condition such as loess cutting, forest bog, ground magnetic survey cannot be carried out work, can only rely on magnetic airborne surveys to complete magnetic survey work.And at the Data processing of magnetic airborne surveys, the magnetic field data of airborne magnetometer collection changes with aircraft altitude, in the time that aircraft altitude is low, actual measurement magnetic field grow, when flying height is high, actual measurement magnetic field dies down.Therefore, the flying height difference between survey line, varies in size the geomagnetic field intensity observing, and causes the difference of magnetic anomaly between line-line, forms ribbon abnormal (as shown in Figure 1), cannot carry out post-processed and the explanation of observation data.So must carry out leveling processing to measurement data, eliminate because of the different spurious anomalies that produce of flying height.Conventional bearing calibration is mainly the method for proofreading and correct with line of cut leveling, but line of cut leveling correction method imperfection is difficult to proofread and correct local anomaly.
Summary of the invention
The invention provides a kind of Aeromagnetic data altitude correction method and system, abnormal to improve because of the vary in size ribbon that causes of geomagnetic field intensity, eliminate because of the different spurious anomalies that produce of flying height.
To achieve these goals, the embodiment of the present invention provides a kind of Aeromagnetic data altitude correction method, and described Aeromagnetic data height antidote comprises:
Obtain the flying height h of aircraft and the Aeromagnetic data T (h) that described aircraft is measured at this flying height h;
The height h to be corrected of default Aeromagnetic data
0;
According to described flying height h, Aeromagnetic data T (h) and height h to be corrected
0generate the Aeromagnetic data after proofreading and correct
Wherein, dT (h)/dh is magnetic field Vertical Derivative, unit: nT/m; Δ h=h
0-h, unit: m.
In one embodiment, described height h to be corrected
0for described aircraft is in the average flight altitude of the practical flight height of each measurement point.
In one embodiment, the height h to be corrected of described default Aeromagnetic data
0, comprising: carry out normal state statistics by the practical flight height to described each measurement point, obtain described average flight altitude.
To achieve these goals, the embodiment of the present invention provides a kind of Aeromagnetic data height correction system, and described Aeromagnetic data height correction system comprises:
Data capture unit, for obtaining the flying height h of aircraft and the Aeromagnetic data T (h) that described aircraft is measured at this flying height h;
Highly default unit, for the height h to be corrected of default Aeromagnetic data
0;
Data correction unit, for according to described flying height h, Aeromagnetic data T (h) and height h to be corrected
0generate the Aeromagnetic data after proofreading and correct
Wherein, dT (h)/dh is magnetic field Vertical Derivative, unit: nT/m; Δ h=h
0-h, unit: m.
In one embodiment, described height h to be corrected
0for described aircraft is in the average flight altitude of the practical flight height of each measurement point.
In one embodiment, the default unit of described height specifically for: carry out normal state statistics by the practical flight height to described each measurement point, obtain described average flight altitude.
The beneficial effect of the embodiment of the present invention is, the present invention has improved because of the geomagnetic field intensity ribbon abnormal occurrence causing that varies in size, and eliminates because of the different spurious anomalies that produce of flying height, for further Aeromagnetic data being processed and provided the foundation.
Brief description of the drawings
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, to the accompanying drawing of required use in embodiment or description of the Prior Art be briefly described below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, do not paying under the prerequisite of creative work, can also obtain according to these accompanying drawings other accompanying drawing.
Fig. 1 is actual measurement Aeromagnetic data planimetric map of the prior art;
Fig. 2 is the Aeromagnetic data altitude correction method process flow diagram of the embodiment of the present invention;
Fig. 3 is embodiment of the present invention flight path model flying height planimetric map;
Fig. 4 is the abnormal planimetric map of actual flight path that flying height changes;
Fig. 5 is the abnormal planimetric map of surface level apart from ground 70m;
Fig. 6 is the Vertical Derivative planimetric map being calculated by panel data;
Fig. 7 is 70 meters of height Vertical Derivative calculated value planimetric maps;
Fig. 8 is the H plane figure after the height correction of gross data Vertical Derivative method;
Fig. 9 is the H plane figure after the height correction of measured data Vertical Derivative method;
Figure 10 A is former survey curve and T1 Repeated Line Tl curve comparison schematic diagram before proofreading and correct;
Figure 10 B is former survey curve and T1 Repeated Line Tl curve comparison schematic diagram after proofreading and correct;
Figure 11 is the Aeromagnetic data height correction system architecture diagram of the embodiment of the present invention.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiment.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtaining under creative work prerequisite, belong to the scope of protection of the invention.
As shown in Figure 2, the embodiment of the present invention provides a kind of Aeromagnetic data altitude correction method, and described Aeromagnetic data height antidote comprises:
Step 201: obtain the flying height h of aircraft and the Aeromagnetic data T (h) that described aircraft is measured at this flying height h;
Step 202: the height h to be corrected of default Aeromagnetic data
0;
Step 203: according to described flying height h, Aeromagnetic data T (h) and height h to be corrected
0generate the Aeromagnetic data after proofreading and correct
Wherein, dT (h)/dh is magnetic field Vertical Derivative, unit: nT/m; Δ h=h0-h, unit: m.
In one embodiment, the height h to be corrected presetting in step 202
0for aircraft is in the average flight altitude of the practical flight height of each measurement point.
In one embodiment, when the concrete enforcement of step 202, can carry out normal state statistics by the practical flight height of each measurement point to aircraft measurement Aeromagnetic data, obtain average flight altitude.Normal state statistical formula is as follows:
(i=1,2,3 ..., N), h: average height; h
i: the height that i is ordered; N: total measuring point.
Height h0 to be corrected default in step 202 can not be also the average flight altitude of aircraft at the practical flight height of each measurement point, can be predefined other level altitude.
When the concrete enforcement of step 203, according to the first approximation of Taylor series, the magnetic survey relation of h0 and h both heights can be written as:
The beneficial effect of the embodiment of the present invention is, the Aeromagnetic data altitude correction method of the embodiment of the present invention has improved because of the geomagnetic field intensity ribbon abnormal occurrence causing that varies in size, eliminate because of the different spurious anomalies that produce of flying height, for further Aeromagnetic data being processed and provided the foundation.
Below in conjunction with concrete example, technique effect of the present invention is described:
1, numerical experimentation
The theoretical magnetic field of underground 5 magnetic spheres under simulation actual flight path (flight path of aircraft) condition of calculating, calculate height off the ground magnetic field of (cross track) when constant simultaneously, as normal data, Vertical Derivative method by the magnetic field under actual flight path condition is proofreaied and correct, with normal data comparison, evaluate and proofread and correct result.
1.1, modelling
5 magnetic spheres parameters of design are in table 1.
Table 1
The magnetic susceptibility of spheroid is 20000CGSM.Normally magnetic field 50000nT, drift angle: 0 degree; Inclination angle 67 is spent.
The simulation flight path that practical flight flight path using power glider in 2013 when the magnetic airborne surveys in many Golcondas area calculates as theory.
Survey line length: 14.325km, survey line number: 143, T1 Repeated Line Tl number: 3, line of cut: 10.
Flight path model as shown in Figure 3, average flight altitude 68.15m, standard deviation 14.25m.
1.2, just drilling result of calculation
According to the theoretical calculation formula of magnetic spheres, just drilling and calculating T (h) and T (70).
T (h): the Δ T of practical flight height is abnormal, as shown in Figure 4.
T (70): the Δ T of (ground is flat) 70m level height is abnormal apart from ground, as shown in Figure 5.Can be illustrated by Fig. 5, flight path height causes obvious ribbon abnormal, and T1 Repeated Line Tl checks that error (mean square deviation) is: 20.28nT.
1.3, the calculating of Vertical Derivative
With line of cut to measured data leveling after, according to potential field theory, by measured data calculate Vertical Derivative, result is as shown in FIG. 6 and 7.
1.4, Vertical Derivative height correction result
As shown in Figure 8, comparison diagram 8 and Fig. 3 can find Vertical Derivative height correction result, and after correction, ribbon extremely obviously improves, and through iterative testing statistics, mean square deviation becomes 3.11nT.
2, the correction of measured data
Fig. 1 is the actual measurement Aeromagnetic data planimetric map of power glider in 2013 boat magnetic system in many Golcondas of Heihe City area, Heilongjiang Province, as background technology is introduced, extremely obvious along survey line ribbon in Fig. 1, brings difficulty to the post-processed of data.
Fig. 9 is the height correction result that adopts Aeromagnetic data altitude correction method of the present invention (Vertical Derivative method), as seen from Figure 9, extremely substantially disappears along the ribbon of survey line.Figure 10 A and Figure 10 B are the former survey line of L508 line and the contrast situation schematic diagram of inspection line before and after proofreading and correct, and Figure 10 A is for proofreading and correct front former survey curve and T1 Repeated Line Tl curve comparison schematic diagram, and Figure 10 B is former survey curve and T1 Repeated Line Tl curve comparison schematic diagram after proofreading and correct.After proofreading and correct by Aeromagnetic data altitude correction method of the present invention, 2 flight is extremely identical, and precision is 4.0nT greatly.
As shown in figure 11, the embodiment of the present invention provides a kind of Aeromagnetic data height correction system, and described Aeromagnetic data height correction system comprises: data capture unit 1101, highly default unit 1102 and Data correction unit 1103.
Data capture unit 1101 can be for obtaining the flying height h of aircraft and the Aeromagnetic data T (h) that aircraft is measured at this flying height h.
Highly default unit 102 can be for the height h to be corrected of default Aeromagnetic data
0.In one embodiment, described height h to be corrected
0for aircraft is in the average flight altitude of the practical flight height of each measurement point, this highly default unit 102 can carry out normal state statistics by the practical flight height to each measurement point, obtains above-mentioned average flight altitude.
Data correction unit 103 can, for according to described flying height h, Aeromagnetic data T (h) and height h0 to be corrected, generate the Aeromagnetic data after proofreading and correct
Wherein, dT (h)/dh is magnetic field Vertical Derivative, unit: nT/m; Δ h=h0-h, unit: m.
The beneficial effect of the embodiment of the present invention is, the present invention has improved because of the geomagnetic field intensity ribbon abnormal occurrence causing that varies in size, and eliminates because of the different spurious anomalies that produce of flying height, for further Aeromagnetic data being processed and provided the foundation.
Those skilled in the art should understand, embodiments of the invention can be provided as method, system or computer program.Therefore, the present invention can adopt complete hardware implementation example, completely implement software example or the form in conjunction with the embodiment of software and hardware aspect.And the present invention can adopt the form at one or more upper computer programs of implementing of computer-usable storage medium (including but not limited to magnetic disk memory, CD-ROM, optical memory etc.) that wherein include computer usable program code.
The present invention is with reference to describing according to process flow diagram and/or the block scheme of the method for the embodiment of the present invention, equipment (system) and computer program.Should understand can be by the flow process in each flow process in computer program instructions realization flow figure and/or block scheme and/or square frame and process flow diagram and/or block scheme and/or the combination of square frame.Can provide these computer program instructions to the processor of multi-purpose computer, special purpose computer, Embedded Processor or other programmable data processing device to produce a machine, the instruction that makes to carry out by the processor of computing machine or other programmable data processing device produces the device for realizing the function of specifying at flow process of process flow diagram or multiple flow process and/or square frame of block scheme or multiple square frame.
These computer program instructions also can be stored in energy vectoring computer or the computer-readable memory of other programmable data processing device with ad hoc fashion work, the instruction that makes to be stored in this computer-readable memory produces the manufacture that comprises command device, and this command device is realized the function of specifying in flow process of process flow diagram or multiple flow process and/or square frame of block scheme or multiple square frame.
These computer program instructions also can be loaded in computing machine or other programmable data processing device, make to carry out sequence of operations step to produce computer implemented processing on computing machine or other programmable devices, thereby the instruction of carrying out is provided for realizing the step of the function of specifying in flow process of process flow diagram or multiple flow process and/or square frame of block scheme or multiple square frame on computing machine or other programmable devices.
In the present invention, applied specific embodiment principle of the present invention and embodiment are set forth, the explanation of above embodiment is just for helping to understand method of the present invention and core concept thereof; , for one of ordinary skill in the art, according to thought of the present invention, all will change in specific embodiments and applications, in sum, this description should not be construed as limitation of the present invention meanwhile.
Claims (6)
1. an Aeromagnetic data altitude correction method, is characterized in that, described Aeromagnetic data height antidote comprises:
Obtain the flying height h of aircraft and the Aeromagnetic data T (h) that described aircraft is measured at this flying height h;
The height h to be corrected of default Aeromagnetic data
0;
According to described flying height h, Aeromagnetic data T (h) and height h to be corrected
0generate the Aeromagnetic data after proofreading and correct
Wherein, dT (h)/dh is magnetic field Vertical Derivative, unit: nT/m; Δ h=h
0-h, unit: m.
2. Aeromagnetic data altitude correction method according to claim 1, is characterized in that, described height h to be corrected
0for described aircraft is in the average flight altitude of the practical flight height of each measurement point.
3. Aeromagnetic data altitude correction method according to claim 2, is characterized in that, the height h to be corrected of described default Aeromagnetic data
0, comprising: carry out normal state statistics by the practical flight height to described each measurement point, obtain described average flight altitude.
4. an Aeromagnetic data height correction system, is characterized in that, described Aeromagnetic data height correction system comprises:
Data capture unit, for obtaining the flying height h of aircraft and the Aeromagnetic data T (h) that described aircraft is measured at this flying height h;
Highly default unit, for the height h to be corrected of default Aeromagnetic data
0;
Data correction unit, for according to described flying height h, Aeromagnetic data T (h) and height h to be corrected
0generate the Aeromagnetic data after proofreading and correct
Wherein, dT (h)/dh is magnetic field Vertical Derivative, unit: nT/m; Δ h=h
0-h, unit: m.
5. Aeromagnetic data height correction system according to claim 4, is characterized in that, described height h to be corrected
0for described aircraft is in the average flight altitude of the practical flight height of each measurement point.
6. Aeromagnetic data height correction system according to claim 5, is characterized in that, the default unit of described height specifically for: carry out normal state statistics by the practical flight height to described each measurement point, obtain described average flight altitude.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410182126.5A CN103926626A (en) | 2014-04-30 | 2014-04-30 | Height correcting method and system for aeromagnetic data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410182126.5A CN103926626A (en) | 2014-04-30 | 2014-04-30 | Height correcting method and system for aeromagnetic data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103926626A true CN103926626A (en) | 2014-07-16 |
Family
ID=51144913
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410182126.5A Pending CN103926626A (en) | 2014-04-30 | 2014-04-30 | Height correcting method and system for aeromagnetic data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103926626A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104729463A (en) * | 2015-03-20 | 2015-06-24 | 广东小天才科技有限公司 | Height correction method and device |
| CN104809352A (en) * | 2015-05-11 | 2015-07-29 | 中国地质大学（北京） | Towed aeromagnetic full-tensor gradient data software compensation method based on forward calculation |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080125920A1 (en) * | 2004-10-08 | 2008-05-29 | Fugro Aiirborne Surveys Corporation | Unmanned Airborne Vehicle For Geophysical Surveying |
| CN102621584A (en) * | 2011-12-14 | 2012-08-01 | 吉林大学 | Aeromagnetic vector data processing method |
-
2014
- 2014-04-30 CN CN201410182126.5A patent/CN103926626A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080125920A1 (en) * | 2004-10-08 | 2008-05-29 | Fugro Aiirborne Surveys Corporation | Unmanned Airborne Vehicle For Geophysical Surveying |
| CN102621584A (en) * | 2011-12-14 | 2012-08-01 | 吉林大学 | Aeromagnetic vector data processing method |
Non-Patent Citations (6)
| Title |
|---|
| 常铮: "磁测资料预处理软件研制", 《中国优秀博硕士学位论文全文数据库（硕士）基础科学辑》 * |
| 徐国苍等: "高精度航磁测量数据各项修正方法", 《中国核科学技术进展报告（第三卷）》 * |
| 方迎尧等: "《DZ/T 0142-94. 航空磁测技术规范》", 27 January 1995, 中华人民共和国地质矿产部 * |
| 简燕等: "高精度磁测正常场改正和高度改正精度", 《物探与化探》 * |
| 袁桂琴等: "地球物理勘查技术与应用研究", 《地质学报》 * |
| 骆遥等: "航磁垂直梯度调整ΔT水平方法研究", 《地球物理学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104729463A (en) * | 2015-03-20 | 2015-06-24 | 广东小天才科技有限公司 | Height correction method and device |
| CN104809352A (en) * | 2015-05-11 | 2015-07-29 | 中国地质大学（北京） | Towed aeromagnetic full-tensor gradient data software compensation method based on forward calculation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103217688B (en) | Airborne laser radar point cloud adjustment computing method based on triangular irregular network | |
| CN103775071B (en) | Adopt the measuring method that coal and rock crack develops | |
| CN102561118B (en) | Quick fine tuning method for ballastless track of high-speed rail based on track deviation | |
| CN108871268B (en) | Tunnel under-excavation numerical calculation method based on laser point cloud | |
| CN105787281A (en) | Fusion inversion method and device for sea wave significant wave height field | |
| CN103207419B (en) | Three-dimensional measurement method for tunnel rock formation attitude | |
| CN102096944A (en) | Cloud recognition method of three-dimensional laser scanning point in geological body structural plane | |
| CN103884520A (en) | Method for manufacturing lunar surface topography simulator | |
| CN106291725A (en) | A kind of method of fast inversion underground geologic bodies locus | |
| CN110333543A (en) | Post non of low resistance body explanation and imaging method and system based on reflection coefficient analysis | |
| CN108427741B (en) | DEM relative error evaluation method based on large number of high-precision control points | |
| CN103256914A (en) | Method and system for calculating inundated areas of warp land dam based on DEM | |
| CN103245324B (en) | Method and system for elevation precision control and correction of island remote sensing mapping | |
| CN111551985A (en) | Underground karst cave detection method based on pile hammer shock and k nearest neighbor algorithm | |
| CN111538106A (en) | Unmanned archaeological exploration method suitable for ultra-shallow water area | |
| CN103926626A (en) | Height correcting method and system for aeromagnetic data | |
| CN104866682A (en) | Ground deformation-based method for inversion of tectonic stress field in shale gas exploration area | |
| CN104330830A (en) | Magnetic body top surface burial depth forecasting method and device | |
| CN108507576B (en) | Three-dimensional flight path planning method for aerial gravity measurement in middle and high mountainous areas | |
| CN108508479B (en) | Method for inverting three-dimensional gravity-magnetic data of open-ground well in cooperation with target position | |
| CN107255837A (en) | A kind of two dimensional terrain forward modeling and correcting method | |
| CN115423954A (en) | Numerical inversion method of nonlinear stress field | |
| CN110968656B (en) | Topographic map processing method and device | |
| CN107219560A (en) | Mine worked-out section deflection and stability assessment method based on gravity anomaly | |
| CN113189660A (en) | Array type land time-varying gravity and gradient field observation method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140716 |