CN102236108B - Three-dimensional terrain correcting method for magnetic surface - Google Patents
Three-dimensional terrain correcting method for magnetic surface Download PDFInfo
- Publication number
- CN102236108B CN102236108B CN 201010169380 CN201010169380A CN102236108B CN 102236108 B CN102236108 B CN 102236108B CN 201010169380 CN201010169380 CN 201010169380 CN 201010169380 A CN201010169380 A CN 201010169380A CN 102236108 B CN102236108 B CN 102236108B
- Authority
- CN
- China
- Prior art keywords
- magnetic
- value
- anomaly
- node
- terrain
- 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.)
- Active
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 take the method for exploration of subsurface formations rock magnetism difference as the basis, utilizes magnetometer observational record changes of magnetic field, processes obtaining subsurface geologic structures through data.Existing petroleum industry standard (" ground magnetometer survey technical regulation (SY/T5771-2004) ") has provided method and the standard of the processing of magnetic exploration data, wherein not requiring to carry out the magnetic topographic correction processes, the distortion effects that landform causes do not have magnetic landscape topographic correction method, so that can not obtain proofreading and correct yet.
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 at present magnetic exploration songization quadratic method to process mountain area magnetic anomalies data; but these class methods belong to Qu Huaping; be about to relief surface and be converted to flatly list processing (LISP); 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 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 the Three-Dimensional Magnetic landform causes of proofreading and correct, can reflect the three-dimensional terrain correcting method for magnetic surface 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, earth's surface magnetic susceptibility distribution grid data.
Step 2) stratum, described earth's surface 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 with the calculation level height value, make up dimensional topography built-up pillar body Model;
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 consists of negative downward upright body, and the center of each main body is each landform node.
4) according to stratum, the earth's surface 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;
Step 4) the upright square body magnetization of described perpendicular magnetization is (1+Q) * k*T
0,
In the formula: Q is 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;
Step 5) described computer capacity is 30-50km.
Step 5) the upright square body magnetic anomaly of the spatial domain perpendicular magnetization computing formula in is:
Wherein: the calculation level coordinate is (x, y, z), the coordinate of landform 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 with the magnetic effect absolute value summation of each terrain mesh node in all computer capacities for calculation level, obtain the magnetic topographic correction value of this calculation level.
7) each is calculated magnetic topographic correction value and the step 1 of measuring point) the change utmost point magnetic anomaly data addition of the corresponding point of exploration acquisition, 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 processed 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 accuracy.
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 front magnetic anomaly for ground, and dotted line changes rear magnetic anomaly for ground) (isoline is apart from 20nT).
Embodiment
Before the invention process Three-Dimensional Magnetic topographic correction, 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 Magnetisability determination 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, earth's surface magnetic susceptibility distribution grid data (Fig. 2).
3) read in calculation level coordinate (x, y) and elevation (h), each node elevation data value of terrain mesh and calculation level height value are compared, 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 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) vertical magnetic anomaly that all landform altitude cylinders produce for measuring point (x, y, h) in the calculation level 40km of adjusting the distance:
Wherein: the calculation level coordinate is (x, y, z), the coordinate of landform 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 the Contrast on effect of seeing that the Three-Dimensional Magnetic topographic correction is forward and backward by Fig. 6.
Claims (6)
1. three-dimensional terrain correcting method for magnetic surface is characterized in that realizing by following steps:
1) magnetic exploration is carried out in the work area, obtain the measuring point magnetic survey, 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, earth's surface magnetic susceptibility distribution grid data;
3) according to the relative height relation of each node elevation data value of terrain mesh with the calculation level height value, make up dimensional topography built-up pillar body Model;
4) according to stratum, the earth's surface 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 with the magnetic effect absolute value summation of each terrain mesh node in all computer capacities for calculation level, obtain the magnetic topographic correction value of this 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. three-dimensional terrain correcting method for magnetic surface according to claim 1 is characterized in that step 2) magnetic susceptibility distribution grid data in stratum, described earth's surface is the panel data grid.
3. three-dimensional terrain correcting method for magnetic surface 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 designated as the positive upright cylinder that makes progress, the landform node grid that height value is lower than calculation level consists of negative downward upright cylinder, and the center of each cylinder is each landform node.
4. three-dimensional terrain correcting method for magnetic surface 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 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. three-dimensional terrain correcting method for magnetic surface according to claim 4 is characterized in that above-mentioned remanent magnetism sense magnetic is 0.1-2.0 than Q.
6. three-dimensional terrain correcting method for magnetic surface according to claim 1 is characterized in that the described computer capacity of step 5) is 30-50km.
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 CN102236108A (en) | 2011-11-09 |
| CN102236108B true 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) |
Families Citing this family (9)
| 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 |
| AU2016228027B2 (en) * | 2015-03-04 | 2018-11-22 | Institute of Mineral Resources, Chinese Academy of Geological Sciences | Method for automatically extracting tectonic framework of potential field |
| 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 |
| 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 | 中国石油天然气股份有限公司 | The method and apparatus for determining Basal ganglia hemorrhage distribution |
| CN109870723B (en) * | 2019-03-18 | 2020-06-23 | 云南航天工程物探检测股份有限公司 | High-power electrical sounding method and system based on mountainous area 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 (3)
| 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 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8055446B2 (en) * | 2006-10-11 | 2011-11-08 | Byerly Kent A | Methods of processing magnetotelluric signals |
-
2010
- 2010-05-06 CN CN 201010169380 patent/CN102236108B/en active Active
Patent Citations (3)
| 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 |
Non-Patent Citations (2)
| Title |
|---|
| 刘云祥.重磁弱异常处理研究与应用.《勘探地球物理进展》.2007,第30卷(第6期), |
| 重磁弱异常处理研究与应用;刘云祥;《勘探地球物理进展》;20071230;第30卷(第6期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102236108A (en) | 2011-11-09 |
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 | |
| CN101636669B (en) | Terrain correction systems | |
| CN102938095A (en) | Mining subsidence monitoring and early warning method based on multi-source data | |
| CN109031301A (en) | Alpine terrain deformation extracting method based on PSInSAR technology | |
| EP2235568A1 (en) | Geophysical data processing systems | |
| CN101950436A (en) | Method for manufacturing digital elevation model by using data of laser altimeter | |
| CN102901985A (en) | Depth domain layer speed correcting method suitable for undulating surface | |
| Oyeniyi et al. | Magnetic surveying as an aid to geological mapping: A case study from Obafemi Awolowo University Campus in Ile-Ife, Southwest Nigeria | |
| CN108427741B (en) | DEM relative error evaluation method based on large number of high-precision control points | |
| Buckley et al. | Combining aerial photogrammetry and terrestrial lidar for reservoir analog modeling | |
| CN110361788B (en) | Air-ground combined three-dimensional gravity data feature analysis and density inversion method | |
| 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 | |
| Whitelaw et al. | High-resolution gravity study of the Gray Fossil Site | |
| CN113552637A (en) | Collaborative three-dimensional inversion method for magnetic anomaly data in aviation-ground-well | |
| CN111399031B (en) | Method for acquiring and processing mountain land seismic data | |
| CN112596113A (en) | Method for identifying field source position based on intersection points of characteristic values of different gradients of gravity | |
| CN102236104B (en) | Three-dimensional density variable-velocity seismic reflection spatial correction method | |
| CN105093324B (en) | A kind of residual static corrections | |
| CN114398696A (en) | Alkaline granite type uranium deposit three-dimensional fine modeling method | |
| CN114674277A (en) | Method for monitoring surface subsidence of deep mining and evaluating land reclamation suitability by combining whole-field survey lines | |
| CN111190224B (en) | Dynamic sampling full-waveform inversion system and method based on three-dimensional seismic wave reverse illumination | |
| Prokhorenko et al. | Topographic correction of GPR profile based on odometer and inclinometer data | |
| Araffa et al. | Assessment of subsurface structural elements around greater Cairo by using integrated geophysical tools |
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 |