CN109752767A - A kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method - Google Patents

A kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method Download PDF

Info

Publication number
CN109752767A
CN109752767A CN201811517204.7A CN201811517204A CN109752767A CN 109752767 A CN109752767 A CN 109752767A CN 201811517204 A CN201811517204 A CN 201811517204A CN 109752767 A CN109752767 A CN 109752767A
Authority
CN
China
Prior art keywords
magnetic
field source
hydrothermal
grid
recognition method
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
Application number
CN201811517204.7A
Other languages
Chinese (zh)
Inventor
喻翔
杨龙泉
赵丹
吴国东
陈聪
周俊杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Research Institute of Uranium Geology
Original Assignee
Beijing Research Institute of Uranium Geology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Research Institute of Uranium Geology filed Critical Beijing Research Institute of Uranium Geology
Priority to CN201811517204.7A priority Critical patent/CN109752767A/en
Publication of CN109752767A publication Critical patent/CN109752767A/en
Pending legal-status Critical Current

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention belongs to hydrothermal-type uranium deposit magnetic data processing fields, and in particular to a kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method.The present invention is the following steps are included: step 1, earth's surface magnetic anomalies gridding processing, all measuring point accurate coordinates of RTK measurement, elevation information;Step 2 is handled for magnetic anomalies grid progress pole;Step 3 seeks total level directional derivative to magnetic anomaly normalizing polar net lattice;Step 4 is filtered the total level directional derivative grid that step 3 obtains using the peaceful smoothing filter of the Chinese;Step 5 compares nodes all in step 4 grid and its week boundary values;Step 6 projects to all peak values in step 5 in database in the obtained magnetic anomaly △ Tization polar net trrellis diagram of cloth rapid five, and field source boundary is then shown in figure.The present invention can effectively, truly Underground uranium source rock body or geophysics field source body boundary position form and range.

Description

A kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method
Technical field
The invention belongs to hydrothermal-type uranium deposit magnetic data processing fields, and in particular to a kind of hydrothermal-type uranium deposit earth object Manage field source Boundary Recognition method.
Background technique
Hydrothermal-type uranium deposit refers to hydrothermal solution containing mine under certain physical and chemical condition, in various advantageous constructions and rock mass In, the valuable mineral accumulation body that is formed by patterns of metallogeny such as filling, explanation and depositions.Including granite type and volcanic type U-ore Bed, it is China's key industry type.Uranium ore output position is generally on the inside of rock mass or outer contacting band.Therefore, uranium source rock The information such as body or geophysics field source body boundary position form, range size, planar distribution feature are particularly important, rock mass side Boundary and outer contacting band or great fracture have important ore-prospecting meaning to the control of rock mass inside, peripheral mineralization.
Due to being influenced by many factors such as oblique magnetization, the superposition of underground field source buried depth, scale, magnetic anomaly and noises, Directly the detection at field source edge is carried out based on earth's surface magnetic anomaly and identification is the difficult point to work at present.And in hydrothermal-type uranium deposit It is previous there is no specific method or not explicitly based on magnetic anomalies progress earth object in geophysical exploration appraisal The method flow for managing field source Boundary Recognition directly affects hydrothermal-type uranium deposit research and appraisal effect, therefore it is urgent to provide a kind of new The method on the identification geophysics field source boundary of type.
Summary of the invention
The technical problem to be solved in the present invention: the present invention provides a kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition Method, can effectively, truly Underground uranium source rock body or geophysics field source body boundary position form and range.
The technical solution adopted by the present invention:
A kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method, successively the following steps are included: step 1, earth's surface Magnetic anomalies gridding processing, RTK measure all measuring point accurate coordinates, elevation information;Step 2, for magnetic anomalies grid into The processing of rowization pole;Step 3 seeks total level directional derivative to magnetic anomaly normalizing polar net lattice;Step 4 utilizes the peaceful smoothing filter of the Chinese The total level directional derivative grid that step 3 obtains is filtered, the superficial part information of short arc, the series of the peaceful filter of the Chinese are suppressed It is set as 0;Step 5 compares nodes all in step 4 grid and its week boundary values, and peak detection rank is chosen to be 3 grades, Peak value is selected, then peak extractions all in grid are come out and are put into database;Step 6, will be all in database in step 5 Peak value project in the obtained magnetic anomaly △ Tization polar net trrellis diagram of cloth rapid five, field source boundary is then shown in figure.
The step 1 includes the following steps: step 1.1, is reconnoitring area's development High-accuracy ground magnetic survey, obtains magnetic survey Data, in order to reach better effect, magnetic survey overall accuracy should be higher than that 2nT, surveys area and is selected in the flat region of landform;Step 1.2, It reconnoitres area and carries out RTK measurement, obtain the accurate latitude and longitude coordinates of measuring point and altitude data;Step 1.3, the magnetic that step 1.1 is obtained Measured data carries out diurnal correction, the correction of normal field gradient, altitude correction and base correction, measuring point exception △ T value is obtained, to magnetic When abnormal data carries out diurnal correction, fundamental magnetic field of the earth parameter need to be searched, this is calculated by the humorous model of ball and reconnoitres area's background Field strength values, magnetic dip angle and magnetic declination;Step 1.4, magnetic anomaly △ T value resulting to step 1.3 carry out data smoothing processing, High-frequency Interference is eliminated, smooth magnetic anomaly curve is obtained and the magnetic anomaly △ T value after all smooth treatments is carried out at gridding Reason, is similarly formed 401 × 402 magnetic anomaly △ T grid file, the preferred minimum-curvature method of gridding difference approach.
The specific steps of the step 2 are as follows: the 1.4 obtained grid data progress pole magnetic anomaly △ T processing rapid to cloth, Magnetic dip angle needed for changing pole processing, magnetic declination are obtained by step 1.3, form 401 × 402 magnetic anomaly normalizing polar net lattice file, The preferred minimum-curvature method of gridding difference approach.
The specific steps of the step 3 are as follows: 2 obtained magnetic anomaly △ Tization rapid to cloth pole grid data seeks total level Directional derivative is sought total level directional derivative method and is calculated using formula (1):
THD indicates to calculate the magnetic anomaly total level directional derivative of point in formula, and Δ T is magnetic anomaly.
In the step 5, all boundary values are referred to: X-direction, Y-direction and two diagonal values;3 grades refer to: nodal value ratio 3 direction value height of surrounding.
In the step 6, geophysics field source boundary is drawn, puts on colour code to peak value, the symbol of drafting passes through long axis side It is shown to by edge trend, these symbols will track geology contact band boundary.
When projecting to peak value in the total level derivative figure that step 3 obtains, these symbols will track localized variation.
The beneficial technical effect of the present invention lies in: after utilizing technical solution of the present invention, to hydrothermal-type uranium deposit geophysics Field source boundary has carried out feature identification, effectively has rated the regularity of distribution of subsurface geophysical field source body;To reach quickly, Mesothermal gold deposits resource is accurately evaluated, there is important practical significance to Exploration.
Detailed description of the invention
Fig. 1 is a kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method flow diagram provided by the present invention.
Specific embodiment
Invention is further described in detail with reference to the accompanying drawings and examples.
As shown in Figure 1, a kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method provided by the invention, is successively wrapped Include following steps:
Step 1, earth's surface magnetic anomalies gridding processing, RTK measure all measuring point accurate coordinates, elevation information;
Step 1.1 is reconnoitring area's development High-accuracy ground magnetic survey, magnetic data is obtained, in order to reach better effect Fruit, magnetic survey overall accuracy should be higher than that 2nT, survey area and are selected in the flat region of landform;
Step 1.2 is reconnoitring area's development RTK measurement (real time kinematic survey system), obtains the accurate latitude and longitude coordinates of measuring point And altitude data;
Step 1.3, to step 1.1 obtain magnetic data carry out diurnal correction, normal field gradient correction, altitude correction and Base correction obtains measuring point exception △ T value, when carrying out diurnal correction to magnetic anomaly regular data, need to search fundamental magnetic field of the earth ginseng Number, calculates this by the humorous model of ball and reconnoitres area's background magnetic field intensity value, magnetic dip angle and magnetic declination;
Step 1.4, magnetic anomaly △ T value resulting to step 1.3 carry out data smoothing processing, eliminate High-frequency Interference, obtain Smooth magnetic anomaly curve simultaneously carries out gridding processing to the magnetic anomaly △ T value after all smooth treatments, it is similarly formed 401 × 402 magnetic anomaly △ T grid file, the preferred minimum-curvature method of gridding difference approach;
Step 2 is handled for magnetic anomalies grid progress pole;
The 1.4 obtained grid data progress pole magnetic anomaly △ T processing rapid to cloth, magnetic dip angle needed for changing pole processing, Magnetic declination is obtained by step 1.3, forms 401 × 402 magnetic anomaly normalizing polar net lattice file, and gridding difference approach is preferably minimum Curvature method;
Step 3 seeks total level directional derivative to magnetic anomaly normalizing polar net lattice;
2 obtained magnetic anomaly △ Tization rapid to cloth pole grid data seeks total level directional derivative, seeks total level direction Derivative method is calculated using formula (1):
THD indicates to calculate the magnetic anomaly total level directional derivative of point in formula, and Δ T is magnetic anomaly.
Step 4 is filtered the total level directional derivative grid that step 3 obtains using the peaceful smoothing filter of the Chinese, compacting The series of the superficial part information of short arc, the peaceful filter of the Chinese is set as 0;
Step 5, nodes all in step 4 grid and its week boundary values are compared (X-direction, Y-direction and two it is diagonal Value), peak detection rank is chosen to be 3 grades (i.e. nodal value is higher than 3 direction value of surrounding), selects peak value, then by institute in grid There is peak extraction to come out and is put into database;
All peak values in step 5 in database are projected to the obtained magnetic anomaly △ Tization pole of cloth rapid five by step 6 On grid chart, field source boundary is then shown in figure.Geophysics field source boundary is drawn, puts on colour code, the symbol of drafting to peak value Edge (geology contact relation) trend is shown by long axis direction, these symbols will track geology contact band boundary.When When peak value being projected in the total level derivative figure that step 3 obtains, these symbols will track localized variation.

Claims (7)

1. a kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method, it is characterised in that: successively the following steps are included: step Suddenly (1), earth's surface magnetic anomalies gridding processing, RTK measure all measuring point accurate coordinates, elevation information;Step (2) is directed to magnetic The processing of power pixilated grid progress pole;Step (3) seeks total level directional derivative to magnetic anomaly normalizing polar net lattice;Step (4), benefit The total level directional derivative grid that step (3) obtains is filtered with the Chinese peaceful smoothing filter, suppresses the superficial part letter of short arc Breath, the series of the peaceful filter of the Chinese are set as 0;Step (5) compares all nodes and its week boundary values in step (4) grid, Peak detection rank is chosen to be 3 grades, selects peak value, then peak extractions all in grid are come out and are put into database;Step (6), all peak values in step (5) in database are projected in the obtained magnetic anomaly △ Tization polar net trrellis diagram of cloth rapid five, Field source boundary is then shown in figure.
2. hydrothermal-type uranium deposit geophysics field source Boundary Recognition method according to claim 1, it is characterised in that: described Step (1) includes the following steps: step (1.1), is reconnoitring area's development High-accuracy ground magnetic survey, obtains magnetic data, is Reach better effect, magnetic survey overall accuracy should be higher than that 2nT, surveys area and is selected in the flat region of landform;Step (1.2) is being reconnoitred Area carries out RTK measurement, obtains the accurate latitude and longitude coordinates of measuring point and altitude data;Step (1.3), the magnetic that step (1.1) are obtained Measured data carries out diurnal correction, the correction of normal field gradient, altitude correction and base correction, measuring point exception △ T value is obtained, to magnetic When abnormal data carries out diurnal correction, fundamental magnetic field of the earth parameter need to be searched, this is calculated by the humorous model of ball and reconnoitres area's background Field strength values, magnetic dip angle and magnetic declination;Step (1.4), magnetic anomaly △ T value resulting to step (1.3) carry out data smoothing High-frequency Interference is eliminated in processing, is obtained smooth magnetic anomaly curve and is carried out grid to the magnetic anomaly △ T value after all smooth treatments Change processing, is similarly formed 401 × 402 magnetic anomaly △ T grid file, the preferred minimum-curvature method of gridding difference approach.
3. hydrothermal-type uranium deposit geophysics field source Boundary Recognition method according to claim 2, it is characterised in that: described The specific steps of step (2) are as follows: pole processing is changed in (1.4) the grid data progress pole obtained magnetic anomaly △ T processing rapid to cloth Required magnetic dip angle, magnetic declination are obtained by step (1.3), form 401 × 402 magnetic anomaly normalizing polar net lattice file, gridding The preferred minimum-curvature method of difference approach.
4. hydrothermal-type uranium deposit geophysics field source Boundary Recognition method according to claim 3, it is characterised in that: described The specific steps of step (3) are as follows: (2) rapid to cloth obtained magnetic anomaly △ Tization pole grid data seeks total level directional derivative, Total level directional derivative method is sought to be calculated using formula (1):
THD indicates to calculate the magnetic anomaly total level directional derivative of point in formula, and Δ T is magnetic anomaly.
5. hydrothermal-type uranium deposit geophysics field source Boundary Recognition method according to claim 4, it is characterised in that: described In step (5), all boundary values are referred to: X-direction, Y-direction and two diagonal values;3 grades refer to: nodal value is than 3 sides of surrounding It is high to value.
6. hydrothermal-type uranium deposit geophysics field source Boundary Recognition method according to claim 5, it is characterised in that: described In step (6), geophysics field source boundary is drawn, puts on colour code to peak value, the symbol of drafting is walked edge by long axis direction To showing, these symbols will track geology contact band boundary.
7. hydrothermal-type uranium deposit geophysics field source Boundary Recognition method according to claim 6, it is characterised in that: as general When peak value is projected in the total level derivative figure that step (3) obtains, these symbols will track localized variation.
CN201811517204.7A 2018-12-12 2018-12-12 A kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method Pending CN109752767A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811517204.7A CN109752767A (en) 2018-12-12 2018-12-12 A kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811517204.7A CN109752767A (en) 2018-12-12 2018-12-12 A kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method

Publications (1)

Publication Number Publication Date
CN109752767A true CN109752767A (en) 2019-05-14

Family

ID=66403728

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811517204.7A Pending CN109752767A (en) 2018-12-12 2018-12-12 A kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method

Country Status (1)

Country Link
CN (1) CN109752767A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113933260A (en) * 2021-09-15 2022-01-14 核工业北京地质研究院 Identification method of hydrothermal uranium deposit fluid activity center

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070162235A1 (en) * 2005-08-25 2007-07-12 Schlumberger Technology Corporation Interpreting well test measurements
CN104965232A (en) * 2015-06-04 2015-10-07 中国地质科学院矿产资源研究所 Automatic extraction method of magnetic structure grillwork in low latitude region
CN106291725A (en) * 2015-05-13 2017-01-04 核工业北京地质研究院 A kind of method of fast inversion underground geologic bodies locus
CN106650192A (en) * 2015-10-30 2017-05-10 核工业北京地质研究院 Volcanic type uranium ore deposit magnetic interface retrieval method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070162235A1 (en) * 2005-08-25 2007-07-12 Schlumberger Technology Corporation Interpreting well test measurements
CN106291725A (en) * 2015-05-13 2017-01-04 核工业北京地质研究院 A kind of method of fast inversion underground geologic bodies locus
CN104965232A (en) * 2015-06-04 2015-10-07 中国地质科学院矿产资源研究所 Automatic extraction method of magnetic structure grillwork in low latitude region
CN106650192A (en) * 2015-10-30 2017-05-10 核工业北京地质研究院 Volcanic type uranium ore deposit magnetic interface retrieval method

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
RICHARD J. BLAKELY ET AL.: "Approximating edges of source bodies from magnetic or gravity anomalies", 《GEOPHYSICS》 *
余钦范等: "《青藏高原磁力测量与地壳结构研究》", 31 July 1993, 地质出版社 *
喻翔等: "诸广山岩体重磁异常特征与铀成矿关系", 《世界核地质科学》 *
朱黎明等: "基于Canny理论的磁场源边界检测", 《工程地球物理学报》 *
杨斯涵: "重磁位场分离及边界识别方法研究", 《中国博士学位论文全文数据库 基础科学辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113933260A (en) * 2021-09-15 2022-01-14 核工业北京地质研究院 Identification method of hydrothermal uranium deposit fluid activity center
CN113933260B (en) * 2021-09-15 2024-06-11 核工业北京地质研究院 Identification method for fluid activity center of hydrothermal uranium deposit

Similar Documents

Publication Publication Date Title
Ravat Analysis of the Euler method and its applicability in environmental magnetic investigations
CN103824133B (en) A kind of granite type U-ore field prospective area Comprehensive Prediction Method
CN104597519B (en) A kind of deep gold mine staged method of prospecting
CN102798898B (en) Three-dimensional inversion method for nonlinear conjugate gradient of magnetotelluric field
CN103913781B (en) Combined prospection method for predicting uranium ore body burial depth
CN106650192B (en) A kind of Volcanic Type Uranium Deposits magnetic interface inversion method
Kalscheuer et al. Two-dimensional magnetotelluric modelling of ore deposits: Improvements in model constraints by inclusion of borehole measurements
CN109782360A (en) A kind of low-resistance area of coverage deep sand body detection method
CN105589108A (en) Rapid three-dimensional inversion method for transient electromagnetism based on different constraint conditions
CN106405664B (en) A kind of magnetic anomaly normalizing pole method
CN104727813B (en) The porosity measurement method on one introduces a collection storage symbiotic type stratum
CN103675944A (en) Comprehensive geophysical method for sedimentary basin uranium resource exploration
CN109902315A (en) A method of delineation Hidden Granite rock mass deep boundary
Yan et al. Concealed faults and intrusions identification based on multiscale edge detection and 3D inversion of gravity and magnetic data: A case study in Qiongheba area, Xinjiang, Northwest China
CN108008456B (en) A method of delineation mesothermal gold deposits deep three-dimensional emphasis U metallogeny Favourable Target Areas
CN109298464A (en) Tight sandstone reservoir Diagenetic Facies Logging Identification Method and device
Wack et al. Cenozoic magnetostratigraphy and magnetic properties of the southern Issyk-Kul basin, Kyrgyzstan
CN109839670A (en) A kind of hydrothermal-type uranium deposit substrate interface inversion method
CN109752767A (en) A kind of hydrothermal-type uranium deposit geophysics field source Boundary Recognition method
CN106289241A (en) Utilize positioning correction method and the device of magnetic signature
Obiora et al. Interpretation of aeromagnetic data of Idah area in north central Nigeria using combined methods
Kaikkonen et al. Audiomagnetotelluric measurements across the Lake Ladoga—Bothnian Bay Zone in Central Finland
CN117115360B (en) Method and system for searching favorable ore forming part of granite type tin-niobium-tantalum-lithium ore deposit
Fusheng et al. Structural setting of the Zoujiashan-Julong’an region, Xiangshan volcanic basin, China, interpreted from modern CSAMT data
CN110058318A (en) A kind of one-dimensional square bearing calibration of the natural electric field based on measuring electrode coordinate

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20190514

RJ01 Rejection of invention patent application after publication