CN104502987A - Refined processing method for measurement data of geophysical prospecting instantaneous change electromagnetic system - Google Patents
Refined processing method for measurement data of geophysical prospecting instantaneous change electromagnetic system Download PDFInfo
- Publication number
- CN104502987A CN104502987A CN201410749348.0A CN201410749348A CN104502987A CN 104502987 A CN104502987 A CN 104502987A CN 201410749348 A CN201410749348 A CN 201410749348A CN 104502987 A CN104502987 A CN 104502987A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- math
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 15
- 238000003672 processing method Methods 0.000 title claims abstract description 8
- 230000005540 biological transmission Effects 0.000 claims abstract description 16
- 230000001052 transient effect Effects 0.000 claims description 20
- 238000011835 investigation Methods 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 4
- 238000000034 method Methods 0.000 abstract description 15
- 101100087528 Mus musculus Rhoj gene Proteins 0.000 abstract 1
- 238000012545 processing Methods 0.000 description 19
- 238000004364 calculation method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000009933 burial Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000006698 induction Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000007670 refining Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a refined processing method for measurement data of a geophysical prospecting instantaneous change electromagnetic system. The method comprises the following steps: collecting a series of inducted secondary field values V(ti) of different delay time ti; obtaining the prospecting depth Hi, and the electrical resistivity rhoi, the time constant taui and the conductivity Si corresponding to the Hi as show in the specification, wherein i is a positive integer; the corresponding relationship between the meaning of each parameter and the unit is as follows: M is magnetic moment of a transmission wire frame; M=I/a/b; I is transmission current; a and b are edge lengths of the transmission wire frame; M is ampere/m<2>; I is ampere; a is m; b is m; q is dipole moment of a receiving coil; q=N/SR; N is turns per coil; SR is the area of the receiving coil; q is m<2>; N is unitless; SR is m<2>, mu0 is magnetic conductivity in vacuum; mu0=4phi*10<-7> henry/m; rho is ohm/m; tau is second; S is siemens; H is m. The refined processing method for the measurement data of the geophysical prospecting instantaneous change electromagnetic system disclosed by the invention is accurate in result and wide in application range.
Description
Technical Field
The invention relates to the technical field of geological exploration, in particular to a method for finely processing measured data of a geophysical prospecting transient electromagnetic system.
Background
With the rapid development of economic society, the demand of mineral resources is sharply increased, the reserve of the reserved resources of a plurality of metal minerals is rapidly attenuated or exhausted, the contradiction of sustainable development is increasingly highlighted, and thus, effective ways and methods for resolving the contradiction fall on the major breakthrough of geological prospecting. At the present stage, the ore searching situation is more and more severe, the ore at the outcrop surface and the ore bodies buried in the shallow part are gradually rare, the geological ore searching in the middle and deep part is necessary, but the difficulty of ore searching in the middle and deep part is multiplied.
The geophysical prospecting method for prospecting minerals in middle and deep geology is indispensable and has more and more important functions. Geophysical prospecting transient electromagnetic prospecting is one of the more available and effective methods at present, and analysis, research, processing and explanation of the prospecting data, namely the data inversion processing technology (solving the shape, scale, burial depth, attitude and the like of an ore body) is a difficult problem which is continuously researched and deeply explored by a plurality of geophysical prospecting workers for many years.
The current data processing method of the transient electromagnetic system has many problems, such as: the data processing process is complicated, only one component is processed at a time, the information is incomplete, and few or no parameters directly reflecting the conductivity of the underground geologic body are available, so that the investigation effect of the transient electromagnetic system is seriously influenced. Therefore, a refinement technique for geophysical prospecting data is needed to directly reflect the electrical parameters of geologic body characteristics, quantitatively determine the buried depth, morphology and attitude characteristics of an ore body, and practically improve the geological prospecting effect.
Disclosure of Invention
The invention aims to solve the technical problem of providing a method for processing the measurement data of the geophysical prospecting transient electromagnetic system, which has accurate and reliable results and wide application range.
Therefore, the technical scheme of the invention is as follows:
a refinement processing method for measurement data of a geophysical prospecting transient electromagnetic system collects a series of different delay times tiInduced secondary field value V (t)i) To obtain a depth of investigation HiAnd with HiCorresponding resistivityρiTime constant τiConductivity SiThe following were used:
wherein I is a positive integer, the meaning of each parameter and the corresponding relation between the meaning and the unit are that M is the magnetic moment of the transmission wire frame, M is I.a.b, I is the transmission current, a and b are the side length of the transmission wire frame, and M-ampere.m2I-ampere, a-meter, b-meter, q is receiving coil dipole moment, q is N · SRN is the number of turns of the coil, SRIs the area of the receiving coil, q-meter2N-unitless, SR-rice2，μ0Is magnetic permeability, mu, in vacuum0＝4π×10-7Henry/meter, rho-ohm meter, tau-second, S-Siemens, H-meter.
The invention provides a new technical thought and means, a field measuring instrument carries out data acquisition according to a large constant source loop device of a time domain transient electromagnetic system, and carries out inversion processing and graph compiling and drawing on the acquired data indoors by using a refined calculation formula, so that a conclusion according with objective geological actual conditions can be obtained, the factors such as the form, the scale, the burial depth, the shape and the like of an ore body are accurately determined, and the effect of prospecting and finding the ore is effectively improved.
Detailed Description
It should be noted that the embodiments and features of the embodiments of the present application may be combined with each other without conflict.
According to the method for finely processing the measurement data of the geophysical prospecting transient electromagnetic system, provided by the embodiment of the invention, a series of different delay times t are collectediInduced secondary field value V (t)i) To obtain a depth of investigation HiAnd with HiCorresponding resistivity ρiTime constant τiConductivity SiThe following were used:
wherein I is a positive integer, the meaning of each parameter and the corresponding relation between the meaning and the unit are that M is the magnetic moment of the transmission wire frame, M is I.a.b, I is the transmission current, a and b are the side length of the transmission wire frame, and M-ampere.m2I-ampere, a-meter, b-meter, q is receiving coil dipole moment, q is N · SRN is the number of turns of the coil, SRIs the area of the receiving coil, q-meter2N-unitless, SR-rice2，μ0Is magnetic permeability, mu, in vacuum0＝4π×10-7Henry/meter, rho-ohm meter, tau-second, S-Siemens, H-meter.
The method for refining the measurement data of the geophysical prospecting transient electromagnetic system comprises two parts, namely field transient electromagnetic system measurement and indoor refined inversion processing. The field transient electromagnetic system measurement comprises a generator, a transmitter, a transmitting loop, a receiver, an induction probe and the like which are hardware parts, wherein the generator is used for measuring a power supply in the working process; the transmitter is used for inputting transient current to the transmission loop to generate a primary magnetic field excitation signal; the receiver is connected with the induction probe and used for collecting secondary field response signals generated by the underground geologic body to obtain measurement data. The refined inversion processing part comprises a data input module, a data processing module and a data storage and display module, wherein the data input module completes data synthesis according to the measuring device and inputs data in a corresponding format; the data processing module completes the refined inversion processing calculation of the multi-parameter data; and the data storage and display module is used for storing data and displaying a processing result and a graph.
According to the invention, a new technical thought is provided, a field measuring instrument carries out data acquisition according to a large constant source loop device of a time domain transient electromagnetic system, and carries out inversion processing and graph compiling on the acquired data indoors by using a refined calculation formula, so that a visual conclusion which accords with the actual situation of objective geology can be obtained, the factors such as the form, the scale, the burial depth, the shape and the like of an ore body are accurately determined, and the effect of prospecting and finding the ore in geological exploration is effectively improved.
The geophysical prospecting transient electromagnetic system measurement data refining inversion processing technology adopted by the invention comprises the following steps:
firstly, according to the field geological condition and the working design requirement, a large fixed source loop device of a time domain transient electromagnetic system is arranged, and a series of different delay times t are collected under the condition of setting a transmitting wire frame magnetic moment M and a receiving coil effective dipole moment qiInduced secondary field value V (t)i). Wherein,Vx(ti)、Vy(ti)、Vz(ti) Is X and Y on the same measuring point,Y, Z measured in three directions.
Secondly, formula derivation is carried out on the basis of the theory research of the para-field, and multivariate physical quantities which can visually reflect the electrical characteristics and the depth of the underground geologic body are obtained, such as: resistivity, time constant, conductivity, depth of investigation.
And thirdly, adopting a calculation formula deduced in the second step to provide a refined processing technical idea, constructing a mathematical model and compiling inversion software to process field measurement data. Acquiring a series of different delay times tiInduced secondary field value V (t)i) To obtain a depth of investigation HiAnd with HiCorresponding resistivity ρiTime constant τiConductivity SiThe following were used:
wherein I is a positive integer, the meaning of each parameter and the corresponding relation between the meaning and the unit are that M is the magnetic moment of the transmission wire frame, M is I.a.b, I is the transmission current, a and b are the side length of the transmission wire frame, and M-ampere.m2I-ampere, a-meter, b-meter, q is receiving coil dipole moment, q is N · SRN is the number of turns of the coil, SRIs the area of the receiving coil, q-meter2N-unitless, SR-rice2，μ0Is magnetic permeability, mu, in vacuum0＝4π×10-7Henry/meter, rho-ohm meter, tau-second, S-Siemens, H-meter.
Namely, a group of resistivity, time constant, conductivity and detection depth value can be obtained by selecting an induction secondary field value with a delay time. It is obvious that i can take a series of positive integers corresponding to different delay times, which is common knowledge of those skilled in the art and will not be described again.
And fourthly, repeating the step III in sequence, finely processing a series of induction secondary field values with different delay times acquired from various field measuring points to finish an inversion process, namely achieving the purposes of effectively processing and correctly interpreting data, and drawing a depth-resistivity section diagram, a depth-time constant section diagram and a depth-conductivity section diagram on a measured section of the transient electromagnetic system by using the processed data. On the depth-resistivity, depth-time constant and depth-conductivity section diagram, each electrical parameter can effectively reflect the characteristics of the geologic body such as distribution, burial depth, form and occurrence, etc., and is very consistent with the real geologic condition, and the geologic prospecting effect is very outstanding.
The refined calculation processing method is a technology for deducing an inverse mathematical model through a mathematical formula and compiling software for method application on the basis of the theoretical analysis and research on the geophysical transient electromagnetic field. The software consists of a data input module, a processing module and a storage and display module. In order to meet the requirements of various field instruments and equipment, software is improved and perfected for many times, a series of transient electromagnetic system measurement data such as PROTEM, PEM, TerraTEM, V8, GDP-32 and the like can be processed, and the method is very suitable for field production exploration.
In summary, the present invention is not limited to the above-mentioned embodiments, and those skilled in the art can propose other embodiments within the technical teaching of the present invention, but such embodiments are included in the scope of the present invention.
Claims (1)
1. A fine processing method for measurement data of a geophysical prospecting transient electromagnetic system is characterized in that a series of different delay times t are collectediInduced secondary field value V (t)i) To obtain a depth of investigation HiAnd with HiCorresponding resistivity ρiTime constant τiConductivity SiThe following were used:
wherein I is a positive integer, the meaning of each parameter and the corresponding relation between the meaning and the unit are that M is the magnetic moment of the transmission wire frame, M is I.a.b, I is the transmission current, a and b are the side length of the transmission wire frame, and M-ampere.m2I-ampere, a-meter, b-meter, q is receiving coil dipole moment, q is N · SRN is the number of turns of the coil, SRIs the area of the receiving coil, q-meter2N-unitless, SR-rice2，μ0Is magnetic permeability, mu, in vacuum0＝4π×10-7Henry/meter, rho-ohm meter, tau-second, S-Siemens, H-meter.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410749348.0A CN104502987A (en) | 2014-12-09 | 2014-12-09 | Refined processing method for measurement data of geophysical prospecting instantaneous change electromagnetic system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410749348.0A CN104502987A (en) | 2014-12-09 | 2014-12-09 | Refined processing method for measurement data of geophysical prospecting instantaneous change electromagnetic system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104502987A true CN104502987A (en) | 2015-04-08 |
Family
ID=52944398
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410749348.0A Pending CN104502987A (en) | 2014-12-09 | 2014-12-09 | Refined processing method for measurement data of geophysical prospecting instantaneous change electromagnetic system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104502987A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104793268A (en) * | 2015-04-28 | 2015-07-22 | 中国科学院地质与地球物理研究所 | Transient electromagnetic detection blind depth measurement method and device |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080082269A1 (en) * | 2006-07-01 | 2008-04-03 | Kjt Enteprises, Inc. | Method for acquiring and interpreting transient electromagnetic measurements |
| CN102608665A (en) * | 2011-11-01 | 2012-07-25 | 蔡运胜 | Fine treatment technology of data information measured by geophysical prospecting time-domain transient electromagnetic system |
-
2014
- 2014-12-09 CN CN201410749348.0A patent/CN104502987A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080082269A1 (en) * | 2006-07-01 | 2008-04-03 | Kjt Enteprises, Inc. | Method for acquiring and interpreting transient electromagnetic measurements |
| CN102608665A (en) * | 2011-11-01 | 2012-07-25 | 蔡运胜 | Fine treatment technology of data information measured by geophysical prospecting time-domain transient electromagnetic system |
Non-Patent Citations (2)
| Title |
|---|
| 李蜀: "《堤坝安全鉴定新办法与现代化防汛抗洪技术指导全书》", 31 August 2003, 吉林电子出版社 * |
| 蔡运胜: "时间域瞬变电磁系统的应用软件开发与利用", 《物探与化探》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104793268A (en) * | 2015-04-28 | 2015-07-22 | 中国科学院地质与地球物理研究所 | Transient electromagnetic detection blind depth measurement method and device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102096113B (en) | Time Domain Earth-Space Electromagnetic Detection System and Calibration Method | |
| CN102183792B (en) | Artificial source frequency domain electromagnetic sounding device and measurement method | |
| CN102426393B (en) | Electric prospecting method | |
| CN103869371B (en) | Manual field source frequency domain full-gradient electromagnetic measuring method | |
| CN103389514B (en) | A kind of mt denoising device and mt denoising method | |
| CN101283292A (en) | Multi-transient DC resistivity measurements | |
| CN202330736U (en) | Electric-method exploration device | |
| CN103207412A (en) | Method for detecting solution leaching and groundwater pollution scope of acid in-situ leaching of uranium | |
| CN104020496A (en) | Ground controlled source magnetotelluric method based on axial collinear manner | |
| CN105204073B (en) | A Tensor Apparent Conductivity Measurement Method | |
| CN105824052A (en) | Underground target body transient electromagnetic multicomponent synthesis method and device | |
| Yang et al. | Cone-shaped source characteristics and inductance effect of transient electromagnetic method | |
| CN105301664B (en) | Artificial source's tensor electromagnetic exploration method that a kind of band is far referred to | |
| CN110488365A (en) | A kind of multipolarization telluric electromagnetic sounding method | |
| CN103852795A (en) | Method for extracting magnetic anomaly signals of underwater small targets | |
| CN104502986B (en) | Physical prospecting induced electrical sounding data chromatography processing method | |
| CN108761540B (en) | A kind of frequency domain natural electric field three-dimensional exploitation method | |
| CN102608665A (en) | Fine treatment technology of data information measured by geophysical prospecting time-domain transient electromagnetic system | |
| CN104502987A (en) | Refined processing method for measurement data of geophysical prospecting instantaneous change electromagnetic system | |
| Frazer | Geomagnetic deep sounding with arrays of magnetometers | |
| CN109901226B (en) | Controllable source tensor geoelectromagnetic system and control calculation method thereof | |
| CN104793268B (en) | A blind depth measurement method and device for transient electromagnetic detection | |
| Khalil et al. | Comparative study between filtering and inversion of VLF-EM profile data | |
| CN113504574B (en) | Electromagnetic data quasi-seismic profile display method and system based on normalization function | |
| Abdelaal et al. | Detection of Depth and Location of Subsurface Utilities Using Smart Phone Magnetic Application: A Case Study at Assiut University, Egypt |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150408 |