CN102236108A - Three-dimensional terrain correcting method for magnetic surface - Google Patents
Three-dimensional terrain correcting method for magnetic surface Download PDFInfo
- Publication number
- CN102236108A CN102236108A CN2010101693803A CN201010169380A CN102236108A CN 102236108 A CN102236108 A CN 102236108A CN 2010101693803 A CN2010101693803 A CN 2010101693803A CN 201010169380 A CN201010169380 A CN 201010169380A CN 102236108 A CN102236108 A CN 102236108A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- value
- anomaly
- node
- magnetic anomaly
- 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.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a three-dimensional terrain correcting method for a magnetic surface in geophysical exploration. The method comprises the following steps: acquiring a surface and stratum magnetisability distribution grid data; constructing a three-dimensional terrain combined cylinder model to form a vertical magnetic upright square cylinder; calculating a magnetic anomaly effect value of each magnetic terrain square cylinder for the vertical magnetic anomaly generated by all upright cylinders with corresponding magnetisability to a survey point; taking the absolute values of the magnetic anomaly effect values of node elevations of the survey point, and summing the absolute values; and adding a magnetic terrain correction value of the calculating point to data of a corresponding point of exploration to acquire corrected magnetic anomaly data. By the method, the magnetic anomaly distortion influenced by magnetic terrain is comprehensively corrected, the distorted magnetic anomaly value on a topographic relief position is raised and recovered, the reflection of magnetic anomaly to an underground magnetic target is more accurate, the magnetic anomaly becomes continuous, zonation is remarkably strengthened, and the mountain magnetic exploration information processing accuracy can be improved.
Description
Technical field
The invention belongs to geophysical exploration method, is the data processing method of a kind of ground magnetic exploration.
Background technology
Magnetic exploration is the method for exploration based on subsurface formations rock magnetism difference, utilizes magnetometer observational record changes of magnetic field, obtains subsurface geologic structures through data processing.Existing petroleum industry standard (" ground magnetometer survey technical regulation (SY/T5771-2004) ") has provided the method and the standard of magnetic exploration DATA PROCESSING, wherein not requiring to carry out the magnetic topographic correction handles, also do not have magnetic landscape topographic correction method, the distortion effects that makes landform cause can not obtain proofreading and correct.
Above-mentioned treatment technology still can be suitable in the low precision stage, but in the high precision magnetic exploration, the landform of mountain region magnetic survey particularly, and its harmful effect can not be ignored.Utilize magnetic exploration songization quadratic method to handle mountain area magnetic anomalies data at present; but these class methods belong to Qu Huaping; be about to relief surface and be converted to list processing (LISP) flatly; do not consider the magnetic effect of landform magnetic; do not calculate the measuring point magnetic anomaly distortion that the magnetic landform produces; the songization quadratic method is not a magnetic topographic correction method, can not get high-precision magnetic exploration effect.
Summary of the invention
The purpose of this invention is to provide a kind of magnetic anomaly distortion that three-dimensional magnetic landform causes of proofreading and correct, can reflect the magnetic landscape dimensional topography correcting method of the magnetic source regularity of distribution and numerical value.
The present invention realizes by following steps:
1) magnetic exploration is carried out in the work area, the measuring point magnetic survey of acquisition, change utmost point magnetic anomaly, work area and extend out district's landform altitude grid, geologic map and stratum magnetic statistics;
2) give the magnetic susceptibility value according to stratum magnetic susceptibility with corresponding stratum in the geologic map, obtain stratum, face of land magnetic susceptibility distribution grid data.
Step 2) stratum, described face of land magnetic susceptibility distribution grid data is the panel data grid.
3) according to the relative height relation of each node elevation data value of terrain mesh, make up dimensional topography combination type cylinder models with the calculation level height value;
The described structure three-dimensional array of step 3 type cylinder models is that the landform node grid that height value is higher than calculation level is the positive upright cylinder that makes progress, the landform node grid that height value is lower than calculation level constitutes negative downward upright body, and the center of each main body is each landform node.
4) according to stratum, the face of land magnetic susceptibility grid node value of each cylinder corresponding flat node location in the dimensional topography combination cylinder, give the cylinder magnetization and form the upright square body of perpendicular magnetization;
The upright square body magnetization of the described perpendicular magnetization of step 4) is (1+Q) * k*T
0,
In the formula: Q is a remanent magnetism sense magnetic ratio, surveys in the district and can be taken as constant, and k is the magnetic susceptibility value of corresponding node, T
0For surveying the terrestrial magnetic field mean value in the district, be constant,
Above-mentioned remanent magnetism sense magnetic is 0.1-2.0 than Q.
5) according to the position of calculation level, the upright cylinders with corresponding magnetic susceptibility all in the computer capacity for the vertical magnetic anomaly that measuring point produces, is calculated the magnetic anomaly effect value of each magnetic landform square body by the upright square body magnetic anomaly of spatial domain perpendicular magnetization computing formula;
The described computer capacity of step 5) is 30-50km.
The upright square body magnetic anomaly of spatial domain perpendicular magnetization computing formula in the step 5) is:
Wherein: the calculation level coordinate be (x, y, z), the coordinate of ground annular volume branch be (ξ,
ζ),
μ
0: permeability of vacuum.
6) will be higher than the magnetic anomaly effect value of node elevation of measuring point and the magnetic anomaly effect value that is lower than the measuring point elevation all takes absolute value, and, obtain the magnetic topographic correction value of this calculation level with the magnetic effect absolute value summation of each terrain mesh node in all computer capacities for calculation level.
7) each is calculated the change utmost point magnetic anomaly data addition of magnetic topographic correction value with the corresponding point of step 1) exploration acquisition of measuring point, obtain the magnetic anomaly data behind the magnetic topographic correction.
The present invention is comprehensively proofreaied and correct the magnetic anomaly distortion that is subjected to the magnetic influence of topography; field data is handled and is shown; behind the topographic correction; the magnetic anomaly value of the distortion at topographic relief place obtains raising and recovers; magnetic anomaly is more accurate to the reflection of underground magnetic target; magnetic anomaly becomes continuously, tapability obviously strengthens, and can improve mountain region magnetic exploration Data Processing precision.
Description of drawings
Fig. 1 is landform altitude figure (isoline is apart from 100m);
Fig. 2 is that (isoline is apart from 50*10 for the magnetic susceptibility distribution plan
-5SI);
Fig. 3 is magnetic topographic correction value figure (isoline is apart from 20nT);
Fig. 4 is the change utmost point MAM (isoline is apart from 100nT) before the magnetic topographic correction;
Fig. 5 is the magnetic anomaly data behind the magnetic topographic correction;
Fig. 6 is magnetic anomaly comparison diagram before and after correct in work area of the present invention (solid line changes preceding magnetic anomaly for ground, and dotted line changes the back magnetic anomaly for ground) (isoline is apart from 20nT).
Embodiment
Before the three-dimensional magnetic topographic correction of the invention process, the related data that needs to prepare comprises that exploration work area measuring point extends out the landform altitude grid data (Fig. 1) of computer capacity (40km), actual measurement magnetic force measuring point coordinate data (x, y, h), survey district's stratum magnetic susceptibility and measure statistical form and geologic map, change utmost point magnetic anomaly data, these data are data that conventional magnetic exploration all possesses.
Three-dimensional magnetic topographic correction method of work is implemented as follows:
1) reads in work area utmost point magnetic anomaly (Fig. 4), landform altitude grid (Fig. 1);
2) give the magnetic susceptibility value according to stratum, work area magnetic susceptibility with corresponding stratum in the geologic map, obtain stratum, face of land magnetic susceptibility distribution grid data (Fig. 2).
3) read in the calculation level coordinate (x y) and elevation (h), compares each node elevation data value of terrain mesh and calculation level height value, the node elevation difference that is higher than the calculation level elevation on the occasion of, the node elevation difference that is lower than the calculation level elevation is a negative value;
4) calculate the magnetization of respectively calculating cylinder corresponding flat node location: j=(1+Q) KT
0, T
0: local geomagnetic field intensity (getting 55400.0); The K-rock susceptibility; Q: the surplus sense of rock magnetism is than (getting 0.5).
5) adjust the distance in the calculation level 40km all landform altitude cylinder for measuring point (x, y, h) vertical magnetic anomaly of Chan Shenging:
Wherein: the calculation level coordinate be (x, y, z), the coordinate of ground annular volume branch be (ξ,
ζ),
μ
0: permeability of vacuum.
6) the magnetic anomaly effect value with each node elevation cylinder takes absolute value, and summation, obtains the magnetic topographic correction value of this calculation level.All measuring points are calculated, obtain the magnetic topographic correction value (Fig. 3) of each measuring point in work area.
7) each is calculated change utmost point magnetic anomaly (Fig. 4) addition of magnetic topographic correction value (Fig. 3) with the corresponding point of measuring point, obtain the magnetic anomaly data (Fig. 5) behind the magnetic topographic correction.
Can know by Fig. 6 and to see the forward and backward effect comparison of three-dimensional magnetic topographic correction.
Claims (7)
1. magnetic landscape dimensional topography correcting method is characterized in that realizing by following steps:
1) magnetic exploration is carried out in the work area, the measuring point magnetic survey of acquisition, change utmost point magnetic anomaly, work area and extend out district's landform altitude grid, geologic map and stratum magnetic statistics;
2) give the magnetic susceptibility value according to stratum magnetic susceptibility with corresponding stratum in the geologic map, obtain stratum, face of land magnetic susceptibility distribution grid data.
3) according to the relative height relation of each node elevation data value of terrain mesh, make up dimensional topography combination type cylinder models with the calculation level height value;
4) according to stratum, the face of land magnetic susceptibility grid node value of each cylinder corresponding flat node location in the dimensional topography combination cylinder, give the cylinder magnetization and form the upright square body of perpendicular magnetization;
5) according to the position of calculation level, the upright cylinders with corresponding magnetic susceptibility all in the computer capacity for the vertical magnetic anomaly that measuring point produces, is calculated the magnetic anomaly effect value of each magnetic landform square body by the upright square body magnetic anomaly of spatial domain perpendicular magnetization computing formula;
6) will be higher than the magnetic anomaly effect value of node elevation of measuring point and the magnetic anomaly effect value that is lower than the measuring point elevation all takes absolute value, and, obtain the magnetic topographic correction value of this calculation level with the magnetic effect absolute value summation of each terrain mesh node in all computer capacities for calculation level;
7) each is calculated the change utmost point magnetic anomaly data addition of magnetic topographic correction value with the corresponding point of step 1) exploration acquisition of measuring point, obtain the magnetic anomaly data behind the magnetic topographic correction.
2. magnetic landscape dimensional topography correcting method according to claim 1 is characterized in that step 2) magnetic susceptibility distribution grid data in stratum, the described face of land is the panel data grid.
3. magnetic landscape dimensional topography correcting method according to claim 1, it is characterized in that the described structure three-dimensional array of step 3) type cylinder models is that the landform node grid that height value is higher than calculation level is the positive upright cylinder that makes progress, the landform node grid that height value is lower than calculation level constitutes negative downward upright body, and the center of each main body is each landform node.
4. magnetic landscape dimensional topography correcting method according to claim 1 is characterized in that the upright square body magnetization of the described perpendicular magnetization of step 4) is (1+Q) * k*T
0,
In the formula: Q is a remanent magnetism sense magnetic ratio, surveys in the district and can be taken as constant, and k is the magnetic susceptibility value of corresponding node, T
0For surveying the terrestrial magnetic field mean value in the district, be constant,
5. according to claim 1 or 4 described magnetic landscape dimensional topography correcting methods, it is characterized in that above-mentioned remanent magnetism sense magnetic is 0.1-2.0 than Q.
6. magnetic landscape dimensional topography correcting method according to claim 1 is characterized in that the described computer capacity of step 5) is 30-50km.
7. magnetic landscape dimensional topography correcting method according to claim 1 is characterized in that the upright square body magnetic anomaly of the spatial domain perpendicular magnetization computing formula in the step 5) is:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010169380 CN102236108B (en) | 2010-05-06 | 2010-05-06 | Three-dimensional terrain correcting method for magnetic surface |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010169380 CN102236108B (en) | 2010-05-06 | 2010-05-06 | Three-dimensional terrain correcting method for magnetic surface |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102236108A true CN102236108A (en) | 2011-11-09 |
CN102236108B CN102236108B (en) | 2013-04-24 |
Family
ID=44886941
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010169380 Active CN102236108B (en) | 2010-05-06 | 2010-05-06 | Three-dimensional terrain correcting method for magnetic surface |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102236108B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064124A (en) * | 2013-01-31 | 2013-04-24 | 江苏大学 | Ratio method for correcting topographic influence in electromagnetic prospecting |
WO2016138874A1 (en) * | 2015-03-04 | 2016-09-09 | 中国地质科学院矿产资源研究所 | Method for automatically extracting tectonic framework of potential field |
CN106650192A (en) * | 2015-10-30 | 2017-05-10 | 核工业北京地质研究院 | Volcanic type uranium ore deposit magnetic interface retrieval method |
CN106855642A (en) * | 2016-11-21 | 2017-06-16 | 中国石油天然气股份有限公司 | Method and device for determining distribution of basement volcanic rocks |
CN108061922A (en) * | 2016-11-07 | 2018-05-22 | 北京自动化控制设备研究所 | A kind of distribution magnetic anomaly detection system dynamic magnetic compensation method |
CN108072910A (en) * | 2016-11-18 | 2018-05-25 | 北京自动化控制设备研究所 | A kind of distribution magnetic anomaly detection system environment magnetic compensation method |
CN109870723A (en) * | 2019-03-18 | 2019-06-11 | 云南航天工程物探检测股份有限公司 | High-power electric measuring method of depth and system based on alpine terrain correction |
CN110108951A (en) * | 2019-04-30 | 2019-08-09 | 中国人民解放军海军工程大学 | Realize the underwater electric scalar potential of the ship method that Qu Huazhi converts from the distant to the near |
CN114488344A (en) * | 2022-02-24 | 2022-05-13 | 核工业北京地质研究院 | Gravity-magnetic three-dimensional inversion method and system for hydrothermal uranium mineralization environment detection |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5502687A (en) * | 1993-07-01 | 1996-03-26 | Western Atlas International, Inc. | Method for datumizing seismic data by forward modeling |
CA2369566A1 (en) * | 1999-04-02 | 2000-10-12 | Alan Royce Huffman | A method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
US6615139B1 (en) * | 2002-03-28 | 2003-09-02 | Council Of Scientific & Industrial Research | Digitally implemented method for automatic optimization of gravity fields obtained from three-dimensional density interfaces using depth dependent density |
US20080091354A1 (en) * | 2006-10-11 | 2008-04-17 | Byerly Kent A | Methods of processing magnetotelluric signals |
-
2010
- 2010-05-06 CN CN 201010169380 patent/CN102236108B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5502687A (en) * | 1993-07-01 | 1996-03-26 | Western Atlas International, Inc. | Method for datumizing seismic data by forward modeling |
CA2369566A1 (en) * | 1999-04-02 | 2000-10-12 | Alan Royce Huffman | A method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
US6615139B1 (en) * | 2002-03-28 | 2003-09-02 | Council Of Scientific & Industrial Research | Digitally implemented method for automatic optimization of gravity fields obtained from three-dimensional density interfaces using depth dependent density |
US20080091354A1 (en) * | 2006-10-11 | 2008-04-17 | Byerly Kent A | Methods of processing magnetotelluric signals |
Non-Patent Citations (1)
Title |
---|
刘云祥: "重磁弱异常处理研究与应用", 《勘探地球物理进展》, vol. 30, no. 6, 30 December 2007 (2007-12-30) * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064124B (en) * | 2013-01-31 | 2015-09-16 | 江苏大学 | A kind of ratio approach correcting electromagnetic survey Considering Terrain Effect |
CN103064124A (en) * | 2013-01-31 | 2013-04-24 | 江苏大学 | Ratio method for correcting topographic influence in electromagnetic prospecting |
WO2016138874A1 (en) * | 2015-03-04 | 2016-09-09 | 中国地质科学院矿产资源研究所 | Method for automatically extracting tectonic framework of potential field |
RU2664488C1 (en) * | 2015-03-04 | 2018-08-17 | Инститьют Оф Минерал Рисорсиз, Чайниз Акедеми Оф Джиолоджикал Сайенсиз | Method of automatic generation of potential field data structure |
US10884161B2 (en) | 2015-03-04 | 2021-01-05 | Institute of Mineral Resources, Chinese Academy of Geological Sciences | Method for automatically extracting structural framework from potential field data |
CN106650192A (en) * | 2015-10-30 | 2017-05-10 | 核工业北京地质研究院 | Volcanic type uranium ore deposit magnetic interface retrieval method |
CN106650192B (en) * | 2015-10-30 | 2019-02-26 | 核工业北京地质研究院 | A kind of Volcanic Type Uranium Deposits magnetic interface inversion method |
CN108061922B (en) * | 2016-11-07 | 2019-06-11 | 北京自动化控制设备研究所 | A kind of distribution magnetic anomaly detection system dynamic magnetic compensation method |
CN108061922A (en) * | 2016-11-07 | 2018-05-22 | 北京自动化控制设备研究所 | A kind of distribution magnetic anomaly detection system dynamic magnetic compensation method |
CN108072910A (en) * | 2016-11-18 | 2018-05-25 | 北京自动化控制设备研究所 | A kind of distribution magnetic anomaly detection system environment magnetic compensation method |
CN108072910B (en) * | 2016-11-18 | 2019-06-07 | 北京自动化控制设备研究所 | A kind of distribution magnetic anomaly detection system environment magnetic compensation method |
CN106855642A (en) * | 2016-11-21 | 2017-06-16 | 中国石油天然气股份有限公司 | Method and device for determining distribution of basement volcanic rocks |
CN109870723A (en) * | 2019-03-18 | 2019-06-11 | 云南航天工程物探检测股份有限公司 | High-power electric measuring method of depth and system based on alpine terrain correction |
CN110108951A (en) * | 2019-04-30 | 2019-08-09 | 中国人民解放军海军工程大学 | Realize the underwater electric scalar potential of the ship method that Qu Huazhi converts from the distant to the near |
CN114488344A (en) * | 2022-02-24 | 2022-05-13 | 核工业北京地质研究院 | Gravity-magnetic three-dimensional inversion method and system for hydrothermal uranium mineralization environment detection |
CN114488344B (en) * | 2022-02-24 | 2024-07-12 | 核工业北京地质研究院 | Heavy magnetic three-dimensional inversion method and system for hydrothermal uranium mining environment detection |
Also Published As
Publication number | Publication date |
---|---|
CN102236108B (en) | 2013-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102236108B (en) | Three-dimensional terrain correcting method for magnetic surface | |
Tung et al. | Assessments of serious anthropogenic land subsidence in Yunlin County of central Taiwan from 1996 to 1999 by Persistent Scatterers InSAR | |
CN104111457B (en) | The inspection mutually of a kind of lift rail PSInSAR Ground Subsidence Monitoring result and time-series fusion method | |
CN101636669B (en) | Terrain correction systems | |
CN109031301A (en) | Alpine terrain deformation extracting method based on PSInSAR technology | |
Oyeniyi et al. | Magnetic surveying as an aid to geological mapping: A case study from Obafemi Awolowo University Campus in Ile-Ife, Southwest Nigeria | |
CN101950436A (en) | Method for manufacturing digital elevation model by using data of laser altimeter | |
CN102175213A (en) | GPS leveling height transfer method | |
CN103727920B (en) | The method of height differences by leveling is measured based on geoid's model | |
CN111191307A (en) | Earthwork virtual construction method based on BIM + GIS technology | |
CN102901985A (en) | Depth domain layer speed correcting method suitable for undulating surface | |
CN105137482A (en) | Method for calculating ancient gradient of sedimentary body | |
CN112666614B (en) | Debris flow source static reserve calculation method based on electrical prospecting and digital elevation model | |
CN110455367A (en) | The engineering reactor of joint unmanned plane and High Density Resistivity abandons scale of construction measurement method | |
CN101957193B (en) | Optimization method for sea island reef height transmission | |
CN108427741B (en) | DEM relative error evaluation method based on large number of high-precision control points | |
CN114674277B (en) | Deep mining earth surface subsidence monitoring method by combining full-field survey lines | |
CN113514900A (en) | Density constraint-based spherical coordinate system gravity and gravity gradient joint inversion method | |
CN102565303B (en) | Fast monitoring method for headward erosion rate of gully head | |
Buckley et al. | Combining aerial photogrammetry and terrestrial lidar for reservoir analog modeling | |
CN113552637A (en) | Collaborative three-dimensional inversion method for magnetic anomaly data in aviation-ground-well | |
Whitelaw et al. | High-resolution gravity study of the Gray Fossil Site | |
Sgattoni et al. | Combining single-station microtremor and gravity surveys for deep stratigraphic mapping | |
CN105093324B (en) | A kind of residual static corrections | |
CN114398696A (en) | Alkaline granite type uranium deposit three-dimensional fine modeling method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |