CN114002747A - Geologic body resistivity detection method and system for ground transmission and air reception - Google Patents

Geologic body resistivity detection method and system for ground transmission and air reception Download PDF

Info

Publication number
CN114002747A
CN114002747A CN202111349047.5A CN202111349047A CN114002747A CN 114002747 A CN114002747 A CN 114002747A CN 202111349047 A CN202111349047 A CN 202111349047A CN 114002747 A CN114002747 A CN 114002747A
Authority
CN
China
Prior art keywords
current
step response
magnetic
emission source
resistivity
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
Application number
CN202111349047.5A
Other languages
Chinese (zh)
Other versions
CN114002747B (en
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.)
Southwest University of Science and Technology
Original Assignee
Southwest University of Science and Technology
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 Southwest University of Science and Technology filed Critical Southwest University of Science and Technology
Priority to CN202111349047.5A priority Critical patent/CN114002747B/en
Publication of CN114002747A publication Critical patent/CN114002747A/en
Application granted granted Critical
Publication of CN114002747B publication Critical patent/CN114002747B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/083Controlled source electromagnetic [CSEM] surveying
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • G01R27/20Measuring earth resistance; Measuring contact resistance, e.g. of earth connections, e.g. plates
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/083Controlled source electromagnetic [CSEM] surveying
    • G01V2003/084Sources
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Abstract

The invention discloses a geologic body resistivity detection method and a system for ground emission and aerial reception, which comprises the steps of firstly, acquiring a magnetic total field data time sequence at each detection point in the air after a current emission source is arranged in a region to be detected; secondly, for any detection point, determining a plurality of step response segments corresponding to the emission state in the time sequence of the total magnetic field data at any detection point based on the emission state of the current emission source; then determining the magnetic field value at any detection point under the excitation condition of the current emission source based on the plurality of step response segments; and finally, determining the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the detection points. The method can reduce the noise interference of the alternating electromagnetic field in the environment, and the time sequence of the magnetic total field data is not influenced by the attitude error of the coil, thereby improving the accuracy and reliability of the detection result. Moreover, the receiver is not required to have a higher time sampling rate, and the complexity of circuit design is reduced.

Description

Geologic body resistivity detection method and system for ground transmission and air reception
Technical Field
The invention relates to the technical field of geophysical exploration, in particular to a surface-emitting aerial-receiving geologic body resistivity detection method and system.
Background
The magneto-resistivity method (MMR) is a geophysical electromagnetic prospecting method that artificially introduces non-induced current (usually direct current or low-frequency alternating current) between two points, measures the magnetic field excited by the induced current on the earth surface, and according to the change rule of the magnetic field, solves engineering geological problems such as underground water flow path, underground pipelines, prospecting and the like.
In the areas with strong fluctuation and covered plants, the ground traffic conditions are poor, and the conventional electric and electromagnetic exploration work is difficult to develop. Therefore, related researchers have proposed a semi-aerial electromagnetic detection method, that is, an electromagnetic emission source is placed on the ground, an unmanned aerial vehicle and other aircraft are enabled to carry an induction coil to fly in the air, and electromagnetic response signals generated by a geologic body in a large range are rapidly measured. The existing semi-aviation electromagnetic exploration method has the following defects:
(1) the induction coil measures the time change rate (dB/dt) of magnetic flux, and is easily interfered by alternating electromagnetic field noise existing in the environment;
(2) the aircraft typically suspends a horizontal induction coil below it, measuring the vertical component in dB/dt. However, the induction coil inevitably rotates and swings in flight, so that the observation data contains coil attitude errors, and the calculation and the elimination are difficult, thereby bringing difficulty to the correct analysis and the interpretation of the observation data;
(3) the electromagnetic detection method using an alternating electromagnetic signal requires a relatively precise time synchronization between the transmission source and the receiver, and also requires a relatively high time sampling rate of the receiver, so that the circuit design is relatively complicated.
Based on this, it is urgently needed to provide a geologic body resistivity detection method for ground transmission and air reception.
The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
The invention aims to provide a geologic body resistivity detection method and a geologic body resistivity detection system which are transmitted from the ground and received in the air, wherein alternating electromagnetic field noise interference existing in the environment can be reduced, the time sequence of magnetic total field data cannot be influenced by coil attitude errors, and the accuracy and the reliability of a detection result are improved. Moreover, the receiver is not required to have a higher time sampling rate, and the complexity of circuit design is reduced.
Another object of the present invention is to reduce the weight of the aircraft, further reduce the power consumption of the aircraft, make the aircraft empty and have a long endurance time, and ensure that the magnetic total field detection device mounted has less influence on the center of gravity of the aircraft, so that the aircraft can fly more stably and safely.
In order to achieve the above object, the present invention provides a geologic body resistivity detection method by ground transmission and aerial reception, comprising:
acquiring a magnetic total field data time sequence at each detection point in the air after a current emission source is distributed in a region to be detected on the ground;
for any detection point, determining a plurality of step response segments corresponding to the emission state in the time series of magnetic total field data at the any detection point based on the emission state of the current emission source;
determining a magnetic field value at any probe point under the excitation condition of the current emission source based on the plurality of step response segments;
and determining the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the detection points.
In an embodiment of the present invention, the current emission source includes a positive electrode, a negative electrode and a power supply, and the positive electrode and the negative electrode are uniformly distributed in the region to be detected or in a preset range around the region to be detected;
the positive electrode is connected with the positive electrode of the power supply, and the negative electrode is connected with the negative electrode of the power supply;
the power supply is used for generating current signals of positive and negative bipolar square wave waveforms.
In one embodiment of the present invention, the current emission source includes a first current emission source and a second current emission source;
the positive electrode of the first current emission source is arranged at one side of the region to be detected, and the negative electrode of the first current emission source is arranged at a first position at one side far away from the region to be detected;
the negative electrode of the second current emission source is arranged on the other side of the region to be detected, and the positive electrode of the second current emission source is arranged at a second position far away from the other side of the region to be detected.
In one embodiment of the invention, the time series of magnetic total field data is detected by a magnetic total field detection device mounted on the aircraft;
correspondingly, the determining, based on the emission state of the current emission source, a plurality of step response segments corresponding to the emission state in the time series of total magnetic field data at any detection point specifically includes:
determining a current emission frequency of the current emission source and a flight speed of the aircraft;
determining a target magnetic total field data time sequence in the magnetic total field data time sequence at any detection point based on the current emission frequency and the flying speed;
based on the transmit state, a plurality of step response segments in the target total magnetic field data time series corresponding to the transmit state are determined.
In one embodiment of the present invention, the transmitting state includes transmitting a first current of a positive polarity square waveform and transmitting a second current of a negative polarity square waveform;
correspondingly, the determining the magnetic field value at any detection point under the excitation condition of the current emission source based on the plurality of step response segments specifically includes:
negating the step response in the step response segment corresponding to the second current to obtain a target step response segment;
overlapping the step response segment corresponding to the first current with the step response in the target step response segment to obtain an average step response time sequence;
and calculating the average value of the step response in the time window of the target width in the average step response time sequence, and taking the average value as the magnetic field value.
In an embodiment of the present invention, the determining, based on the emission state of the current emission source, a plurality of step response segments in the time series of total magnetic field data at any detection point corresponding to the emission state further includes:
and filtering out fixed frequency noise and low frequency components of the magnetic total field data at each moment in the time sequence of the magnetic total field data at any detection point.
In an embodiment of the present invention, the determining the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the respective detection points further includes:
determining a type of the geological volume based on the resistivity of the geological volume;
the types of the geologic body at least include a vein, a fracture zone, a karst cave, and groundwater.
The invention provides a geologic body resistivity detection system with ground transmission and aerial reception, which comprises: the device comprises a geologic body resistivity detection device, a current emission source and a magnetic total field detection device, wherein the geologic body resistivity detection device is used for emitting and receiving in the air from the ground, and the magnetic total field detection device is mounted on an aircraft;
the geologic body resistivity detection device is used for executing the above geologic body resistivity detection method of ground transmission and aerial reception;
the current emission source is used for emitting a current signal with positive and negative bipolar square wave waveforms;
the magnetic total field detection device is used for detecting the magnetic total field data time sequence at each detection point in the air.
The invention provides an electronic device, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the program to realize the steps of the geologic body resistivity detection method received in the air by the ground transmission.
The invention provides a non-transitory computer-readable storage medium having stored thereon a computer program, wherein the computer program, when executed by a processor, implements the steps of the method for surface-transmitted aerial-reception resistivity detection of geological volume as described above.
Compared with the prior art, according to the geologic body resistivity detection method and system for ground emission and aerial reception, firstly, a magnetic total field data time sequence at each detection point in the air after a current emission source is distributed in a ground area to be detected is obtained; secondly, for any detection point, determining a plurality of step response segments corresponding to the emission state in the time sequence of the total magnetic field data at any detection point based on the emission state of the current emission source; then determining the magnetic field value at any detection point under the excitation condition of the current emission source based on the step response segments; and finally, determining the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the detection points. The method can reduce the noise interference of the alternating electromagnetic field in the environment, and the time sequence of the magnetic total field data is not influenced by the attitude error of the coil, thereby improving the accuracy and reliability of the detection result. Moreover, the receiver is not required to have a higher time sampling rate, and the complexity of circuit design is reduced.
Drawings
FIG. 1 is a schematic flow diagram of a method for surface-transmitted aerial-received geologic resistivity survey in accordance with one embodiment of the present invention;
FIG. 2 is a schematic view showing the arrangement position of a current emission source according to an embodiment of the present invention;
FIG. 3 is a graph illustrating the values of static magnetic fields generated by subsurface currents in a breaker zone according to one embodiment of the present invention;
FIG. 4 is a schematic structural diagram of a surface-emitting aerial-receiving geologic resistivity survey apparatus according to one embodiment of the present invention;
FIG. 5 is a schematic diagram of a ground-based transmit aerial receive geologic resistivity survey system in accordance with one embodiment of the present invention;
fig. 6 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
The following detailed description of the present invention is provided in conjunction with the accompanying drawings, but it should be understood that the scope of the present invention is not limited to the specific embodiments.
Throughout the specification and claims, unless explicitly stated otherwise, the word "comprise", or variations such as "comprises" or "comprising", will be understood to imply the inclusion of a stated element or component but not the exclusion of any other element or component.
Fig. 1 is a schematic flow chart of a method for detecting resistivity of a geologic body received in the air by ground transmission, which is provided in an embodiment of the present invention, and as shown in fig. 1, the method includes:
s1, acquiring a magnetic total field data time sequence at each detection point in the air after a current emission source is arranged in the region to be detected on the ground;
s2, for any detection point, determining a plurality of step response segments corresponding to the emission state in the time sequence of the total magnetic field data at any detection point based on the emission state of the current emission source;
s3, determining the magnetic field value at any detection point under the excitation condition of the current emission source based on the step response segments;
and S4, determining the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the detection points.
Specifically, in the method for detecting resistivity of a geologic body through ground transmission and aerial reception provided in the embodiment of the present invention, an execution main body is a server, the server may be a local server or a cloud server, and the local server may be a computer, a tablet computer, a smart phone, and the like.
Step S1 is executed first, and the time sequence of the total magnetic field data at each detection point in the air after the current emission source is arranged in the region to be detected on the ground is obtained first. The current emission source is used for providing a current signal with positive and negative bipolar square wave waveforms for the earth of a region to be detected on the ground, so that the region to be detected generates a low-frequency magnetic field, and a time sequence of magnetic total field data at each detection point in the air is conveniently acquired. The current emission sources may be disposed in the region to be detected or in a preset range around the region to be detected, and the number of the current emission sources may be one or more, and may be selected according to the requirement, which is not specifically limited in the embodiment of the present invention.
The region to be detected is a ground region in which the resistivity of the corresponding underground geologic body needs to be determined, a plurality of detection points can be arranged in the air, and the detection points can be uniformly distributed in the air region corresponding to the region to be detected.
The magnetic total field data time sequence refers to magnetic total field data arranged according to a time sequence, and the magnetic total field data time sequence can be obtained by detecting a magnetic total field detection device mounted on an aircraft. The aircraft can include unmanned planes, helicopters, fixed-wing airplanes and the like, and the magnetic total field detection device can include a three-axis fluxgate, an optical pump magnetometer, a three-component superconducting quantum interference magnetometer, a proton magnetometer and the like. In the embodiment of the invention, each detection point in the air is each point of the magnetic total field detection device in the moving track in the air.
When the current emission source emits a current signal, the aircraft flies in the air, and a magnetic total field data time sequence can be obtained through the magnetic total field detection device in the flying process. When the three-component superconducting quantum interference magnetometer is used as a magnetic total field detection device, three-component data, namely magnetic induction intensity data in three orthogonal directions, can be obtained through measurement of the three-component superconducting quantum interference magnetometer, and a magnetic total field data time sequence can be obtained through squaring and then squaring of the three-component data. In the embodiment of the invention, the current emission source and the magnetic total field detection device are synchronized through a GPS clock.
In this process, the position and the route time of each probe point may be recorded by a Global Positioning System (GPS) on the aircraft.
Then, step 2 is executed, and for any detection point in the region to be detected, that is, any detection point, that is, each detection point, according to the emission state of the current emission source, a plurality of step response segments corresponding to the emission state in the time series of the total magnetic field data at the any detection point can be determined. The emission state of the current emission source can comprise the emission of a first current signal with a positive square waveform and the emission of a second current signal with a negative square waveform, and accordingly, a step response segment corresponding to the first current signal and a step response segment corresponding to the second current signal in the time sequence of the total magnetic field data at any detection point can be determined. The step response segment corresponding to the first current signal is a time segment of the step response of the region to be detected to the first current signal, and the step response segment corresponding to the second current signal is a time segment of the step response of the region to be detected to the second current signal. The number of the step response segments corresponding to the first current signal and the number of the step response segments corresponding to the second current signal may be multiple, which is not specifically limited in the embodiment of the present invention.
It should be noted that, the step response of the region to be detected to the first current signal, that is, the magnetic total field data corresponding to the time period when the current emission source emits the first current signal, and the step response of the region to be detected to the second current signal, that is, the magnetic total field data corresponding to the time period when the current emission source emits the second current signal.
Then, step S3 is executed, and in combination with the plurality of step response segments, a magnetic field value at any detection point under the excitation condition of the current emission source can be determined, where the magnetic field value is the magnetic field response of the underground geologic body to the current emission source.
And finally, executing step S4, and determining the resistivity of the underground geologic body corresponding to the region to be detected according to the magnetic field values at the detection points in the air. In the embodiment of the invention, inversion can be carried out according to the magnetic field value at each detection point in the air, so that the three-dimensional spatial distribution of the resistivity of the underground geologic body corresponding to the region to be detected is obtained.
The geologic body resistivity detection method of ground emission aerial reception provided in the embodiment of the invention comprises the steps of firstly, acquiring a magnetic total field data time sequence at each detection point in the air after a current emission source is arranged in a region to be detected; secondly, for any detection point, determining a plurality of step response segments corresponding to the emission state in the time sequence of the total magnetic field data at any detection point based on the emission state of the current emission source; then determining the magnetic field value at any detection point under the excitation condition of the current emission source based on the step response segments; and finally, determining the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the detection points. The method can reduce the noise interference of the alternating electromagnetic field in the environment, and the time sequence of the magnetic total field data is not influenced by the attitude error of the coil, thereby improving the accuracy and reliability of the detection result. Moreover, the receiver is not required to have a higher time sampling rate, and the complexity of circuit design is reduced.
On the basis of the above embodiment, in the embodiment of the present invention, after the magnetic field value at any detection point is determined, inversion imaging may be performed based on the magnetic field values at the detection points to obtain a finer underground magnetic field distribution image, so that a user can more intuitively grasp the distribution condition of the hollow magnetic field.
On the basis of the above embodiment, in the geologic body resistivity detection method of ground emission and aerial reception provided in the embodiment of the present invention, the current emission source includes a positive electrode, a negative electrode and a power supply, and the positive electrode and the negative electrode are uniformly distributed in the area to be detected or in a preset range around the area to be detected;
the positive electrode is connected with the positive electrode of the power supply, and the negative electrode is connected with the negative electrode of the power supply;
the power supply is used for generating current signals of positive and negative bipolar square wave waveforms.
Specifically, the current emission source in the embodiment of the present invention may include a positive electrode, a negative electrode, and a power supply, where the positive electrode and the negative electrode are uniformly distributed in the region to be detected or in a predetermined range around the region to be detected. The positive electrode is connected with the positive electrode of the power supply through a cable, the negative electrode is connected with the negative electrode of the power supply through a cable, and the positive electrode, the ground area between the positive electrode and the negative electrode, the negative electrode and the power supply can directly form a closed loop. The power supply can generate current signals of positive and negative bipolar square wave waveforms so as to ensure that the current emission source can provide the current signals of the positive and negative bipolar square wave waveforms for the earth of a region to be detected.
In the embodiment of the invention, the specific structure of the current emission source is provided, so that the current emission source can provide current signals of positive and negative bipolar square wave waveforms for the ground of the region to be detected in real time, and the detection feasibility of the magnetic field value at each detection point in the region to be detected is further ensured.
On the basis of the above embodiments, the method for detecting resistivity of a geologic body received in the air by ground emission provided in the embodiments of the present invention includes a first current emission source and a second current emission source;
the positive electrode of the first current emission source is arranged at one side of the region to be detected, and the negative electrode of the first current emission source is arranged at a first position at one side far away from the region to be detected;
the negative electrode of the second current emission source is arranged on the other side of the region to be detected, and the positive electrode of the second current emission source is arranged at a second position far away from the other side of the region to be detected.
Specifically, in the embodiment of the present invention, the number of the current emission sources may be 2, that is, the current emission sources include a first current emission source and a second current emission source, and each of the first current emission source and the second current emission source includes a positive electrode, a negative electrode, and a power supply.
The positive electrode of the first current emission source is arranged at one side of the region to be detected, and the negative electrode of the first current emission source is arranged at a first position at one side far away from the region to be detected, wherein the first position can be an equivalent position at infinity at one side of the region to be detected.
The negative electrode of the second current emission source is arranged at the other side of the region to be detected, and the positive electrode of the second current emission source is arranged at a second position far away from the other side of the region to be detected, wherein the second position can be an equivalent position at infinity at the other side of the region to be detected.
As shown in FIG. 2, the area to be detected is a mountain 2, and the broken zone 21 in the mountain 2 is detected, wherein a positive electrode of a first current emission source is arranged at A, and a negative electrode of a second current emission source is arranged at B. The power supplies in the first current emission source and the second current emission source input current signals with certain amplitude to the cable and the mountain underground according to a certain fundamental frequency in a positive and negative bipolar square wave shape.
When a first current emission source provides a current signal, a magnetic total field detection device mounted on an aircraft receives magnetic total field data along a measuring line in the air above a region to be detected to obtain a first group of data; when the second current emission source provides current signals, the magnetic total field data are repeatedly received once at the same line measuring position, and a second group of data are obtained. And two groups of data are respectively subjected to subsequent processing to obtain two groups of magnetic field values, and the two groups of magnetic field values are superposed to obtain the magnetic field values equivalent to that the positive electrodes and the negative electrodes are directly arranged at A, B.
In the embodiment of the invention, two current emission sources can be introduced into the area to be detected, the power supply across the mountain is equivalent through repeated emission and observation, and the difficulty of wiring in the mountain area is greatly reduced.
In the prior art, in order to measure reliable dB/dt data, the equivalent area of an induction coil is increased by using a multi-turn metal loop, the weight is heavier, and the problem of short flight endurance time of an aircraft is caused.
Therefore, on the basis of the above embodiment, in the geologic body resistivity detection method by ground transmission and aerial reception provided in the embodiment of the present invention, the time sequence of the magnetic total field data is detected by a magnetic total field detection device mounted on an aircraft;
correspondingly, the determining, based on the emission state of the current emission source, a plurality of step response segments corresponding to the emission state in the time series of total magnetic field data at any detection point specifically includes:
determining a current emission frequency of the current emission source and a flight speed of the aircraft;
determining a target magnetic total field data time sequence in the magnetic total field data time sequence at any detection point based on the current emission frequency and the flying speed;
based on the transmit state, a plurality of step response segments in the target total magnetic field data time series corresponding to the transmit state are determined.
Specifically, the total magnetic field detection device adopted in the embodiment of the present invention may include a three-axis fluxgate, an optical pump magnetometer, a three-component superconducting quantum interference magnetometer, a proton magnetometer, and the like. The mass of the magnetic total field detection device is lighter than that of a coil type receiver, the air-staying and endurance time is longer under the condition that the power consumption of the aircraft is unchanged, and the influence of the magnetic total field detection device mounted on the aircraft on the gravity center of the aircraft is smaller, so that the aircraft flies more stably and safely.
When a plurality of step response segments corresponding to the emission state in the time sequence of the magnetic total field data at any detection point are determined through the emission state of the current emission source, the current emission frequency of the current emission source and the flight speed of the aircraft can be determined firstly. The transmission signal frequency represents the number of bipolar square waves transmitted by the current transmission source at any detection point in unit time (usually 1s), and the flight speed of the aircraft refers to the distance flown by the aircraft in unit time (usually 1 s).
Because the aircraft does not perform hover measurement after reaching the detection point but performs continuous flight measurement, in order to obtain accurate approximate hover measurement, the magnetic total field data obtained by measurement near the detection point needs to be collected, which needs 1) to improve the current emission frequency, so that the time for measuring the number of the same bipolar square waves is shortened, namely the distance for the aircraft to cross the detection point is shortened; or 2) reducing the flight speed of the aircraft. The current emission frequency and the flight speed of the aircraft may be set according to actual needs, which is not particularly limited in the embodiment of the present invention.
According to the current emission frequency and the flight speed, a target magnetic total field data time sequence in the magnetic total field data time sequence at any detection point can be determined, and the target magnetic total field data time sequence is a result of arranging the magnetic total field data measured nearby the detection point according to the time sequence.
And finally, determining a plurality of step response segments corresponding to the emission state in the time sequence of the target magnetic total field data according to the emission state.
In the embodiment of the invention, when a plurality of step response segments are determined, the time sequence of the target magnetic total field data can be determined firstly, so that the magnetic field value at any detection point can be more accurate.
On the basis of the above embodiment, the method for detecting resistivity of a geologic body received in the air by ground emission provided in the embodiment of the present invention includes the steps of transmitting a first current signal having a positive square waveform and transmitting a second current signal having a negative square waveform in an emission state;
correspondingly, the determining the magnetic field value at any detection point under the excitation condition of the current emission source based on the plurality of step response segments specifically includes:
negating the step response in the step response segment corresponding to the second current signal to obtain a target step response segment;
overlapping the step response segment corresponding to the first current signal with the step response in the target step response segment to obtain an average step response time sequence;
and calculating the average value of the step response in the time window of the target width in the average step response time sequence, and taking the average value as the magnetic field value.
Specifically, in the embodiment of the present invention, the emission state of the current emission source may include emission of a first current signal of a positive polarity square waveform and emission of a second current signal of a negative polarity square waveform. Accordingly, the plurality of step response segments includes a step response segment corresponding to the first current signal and a step response segment corresponding to the second current signal. Therefore, when the magnetic field value at any detection point is determined according to the plurality of step response segments, the step response in the step response segment corresponding to the second current signal may be inverted to obtain a target step response segment, and then the step response segment corresponding to the first current signal and the step response in the target step response segment are superimposed to obtain an average step response time sequence.
And finally, calculating the average value of the step response in the time window of the target width in the average step response time sequence, and taking the average value as the magnetic field value at any detection point. The target width may be set as required, and the time window of the target width may be generally selected in the late stage of the average step response time sequence. The later stage is a relative concept, and can be defined differently according to actual conditions, generally meaning that the electromagnetic induction signal generated by the step response to the electromagnetic field disturbance is dissipated, and only the electric field or the magnetic field generated by the steady field is left. The specific value of this late stage depends on many factors such as the fundamental frequency of the square wave, the earth conductivity, etc. For example, if the current emission frequency is less than 2.5Hz and the step response corresponds to a time period of 0.0001s to 0.1s, then the late phase may be selected to be 0.01s to 0.1 s.
In the embodiment of the invention, a time window with a target width is selected in the late stage of the average step response time sequence, and the obtained magnetic field value at any detection point can be equivalent to the response of the steady magnetic field under the excitation condition of the steady current source.
On the basis of the foregoing embodiments, the method for detecting resistivity of a geologic body received in the air by ground transmission provided in an embodiment of the present invention includes determining, based on a transmission state of the current transmission source, a plurality of step response segments corresponding to the transmission state in a time series of total magnetic field data at any detection point, where the step response segments include:
and filtering out fixed frequency noise and low frequency components of the magnetic total field data at each moment in the time sequence of the magnetic total field data at any detection point.
Specifically, in the embodiment of the present invention, after the time sequence of the total magnetic field data at any detection point is determined, the fixed frequency noise and the low frequency component in the data may be filtered, and then the application is performed, so that the accuracy of the magnetic field value at any detection point may be further ensured.
Fixed frequency noise can be filtered by filtering methods such as empirical mode decomposition and infinite impulse response digital filters, and the fixed frequency noise can comprise power frequency noise and the like. Low-frequency components can be filtered out through Kalman filtering, a least square method, envelope extraction and other methods, and the low-frequency components can comprise motion noise and the like.
On the basis of the foregoing embodiment, the method for detecting resistivity of a geologic body transmitted and received in the air on the ground according to the embodiments of the present invention determines the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the respective detection points, and then further includes:
determining a type of the geological volume based on the resistivity of the geological volume;
the types of the geologic body at least include a vein, a fracture zone, a karst cave, and groundwater.
In particular, in embodiments of the present invention, the magnetic field around the low resistivity object is stronger because the low resistivity object is more prone to sink larger currents, as shown in fig. 3. After the magnetic field values of the detection points in the air are determined, the types of the underground geologic bodies corresponding to the areas to be detected can be determined according to the magnetic field values of the detection points, so that geologic bodies such as ore veins, fault zones (fracture zones), fracture zones, karst caves, underground water and the like can be conveniently and quickly found. The distribution of the MMR anomaly of the steady magnetic field at the Fracture zone (Fracture zone) under excitation of the current-emitting source is shown only in FIG. 3, the magnetic field value being given in nT.
On the basis of the above embodiment, the geologic body resistivity detection method based on ground emission and aerial reception provided in the embodiments of the present invention can replace the position of the current emission source or the aircraft for the region to be detected, and perform multiple detections, so as to realize the omnibearing coverage of different positions and targets moving toward the target in the detection region.
In summary, the geologic body resistivity detection method based on ground transmission and aerial reception provided in the embodiments of the present invention has the following specific advantages:
(1) the observation efficiency is high. Compared with the magneto-resistivity method compared with the conventional ground transmitting and ground receiving, the aircraft can cover a large area in a short time in air receiving, and reaches mountainous areas and vegetation coverage areas where the ground is difficult to pass through. When the device meets the area where the cables are not easily laid on the ground, the device can also use one side of the mountain to emit and combine with data processing for many times, and equivalently power is supplied across two sides of the mountain, so that the cables do not need to cross or bypass the mountain, and the equipment laying time is effectively reduced.
(2) The observation precision is high. The measured data of the method is the magnetic total field under the steady condition. The magnetic total field is not changed due to the rotation of the sensor, so that the influence of attitude error is low; and if the magnetic total field is measured, a magnetic total field detection device with higher precision, such as a magnetic total field magnetometer, an optical pump and the like, can also be used. In addition, the observation object is a magnetic field under a stable condition, is not easy to be subjected to complex alternating electromagnetic fields and aircraft electromagnetic interference in the environment, and is relatively easy to obtain high-precision data through simple data processing.
(3) The takeoff weight is light. The general quality of magnetic total field detection device is lighter than coil formula receiver, under the unchangeable condition of aircraft consumption, has longer time of staying empty and endurance, and the sensor of hanging influences littleer to aircraft focus, and flight is more stable and safe.
In the embodiment of the invention, the magnetic total field detection device carried by the aircraft is adopted for the condition that the area to be detected is the land area, and when the exploration is carried out on water, the magnetic total field detection device carried by a ship can be used for realizing large-area rapid scanning.
As shown in fig. 4, on the basis of the above embodiments, the embodiment of the present invention provides a device for detecting resistivity of a geologic body through ground transmission and air reception, which is used to execute the method for detecting resistivity of a geologic body through ground transmission and air reception provided in the above method embodiments.
The geologic body resistivity detection device of ground transmission aerial reception can include:
an obtaining module 41, configured to obtain a magnetic total field data time sequence at each detection point in the air after a current emission source is arranged in a region to be detected;
a first determining module 42, configured to determine, for any detection point, a plurality of step response segments corresponding to an emission state of the current emission source in a time series of total magnetic field data at the any detection point based on the emission state;
a second determining module 43, configured to determine a magnetic field value at any of the detection points under the excitation condition of the current emission source based on the plurality of step response segments;
and the detection module 44 is configured to determine the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the detection points.
On the basis of the above embodiment, the geobody resistivity detection device for ground emission and aerial reception provided in the embodiment of the present invention includes a current emission source, a current detection unit, a power supply, and a control unit, wherein the current emission source includes a positive electrode, a negative electrode, and a power supply, and the positive electrode and the negative electrode are uniformly distributed in the region to be detected or in a preset range around the region to be detected;
the positive electrode is connected with the positive electrode of the power supply, and the negative electrode is connected with the negative electrode of the power supply;
the power supply is used for generating current signals of positive and negative bipolar square wave waveforms.
On the basis of the above embodiment, the geologic body resistivity detection device received in the air by ground emission provided in the embodiment of the present invention includes a first current emission source and a second current emission source;
the positive electrode of the first current emission source is arranged at one side of the region to be detected, and the negative electrode of the first current emission source is arranged at a first position at one side far away from the region to be detected;
the negative electrode of the second current emission source is arranged on the other side of the region to be detected, and the positive electrode of the second current emission source is arranged at a second position far away from the other side of the region to be detected.
On the basis of the above embodiment, in the geological resistivity detection device which is provided in the embodiment of the present invention and transmits and receives data in the air on the ground, the time sequence of the magnetic total field data is obtained by detecting the magnetic total field detection device mounted on the aircraft;
accordingly, the first determining module is specifically configured to:
determining a current emission frequency of the current emission source and a flight speed of the aircraft;
determining a target magnetic total field data time sequence in the magnetic total field data time sequence at any detection point based on the current emission frequency and the flying speed;
based on the transmit state, a plurality of step response segments in the target total magnetic field data time series corresponding to the transmit state are determined.
On the basis of the above embodiment, the ground emission aerial-reception geologic body resistivity detection device provided in the embodiment of the present invention, wherein the emission state includes emission of a first current signal having a positive square waveform and emission of a second current signal having a negative square waveform;
accordingly, the second determining module is specifically configured to:
negating the step response in the step response segment corresponding to the second current signal to obtain a target step response segment;
overlapping the step response segment corresponding to the first current signal with the step response in the target step response segment to obtain an average step response time sequence;
and calculating the average value of the step response in the time window of the target width in the average step response time sequence, and taking the average value as the magnetic field value.
On the basis of the above embodiment, the geologic body resistivity detection device received in the air by ground transmission provided in the embodiment of the present invention further includes a filtering module, configured to:
and filtering out fixed frequency noise and low frequency components of the magnetic total field data at each moment in the time sequence of the magnetic total field data at any detection point.
On the basis of the above embodiment, the geologic body resistivity detection apparatus received in the air by ground transmission provided in the embodiment of the present invention further includes a third determination module, configured to:
determining a type of the geological volume based on the resistivity of the geological volume;
the types of the geologic body at least include a vein, a fracture zone, a karst cave, and groundwater.
As shown in fig. 5, on the basis of the above embodiments, an embodiment of the present invention provides a system for detecting resistivity of a geologic body by ground transmission and air reception, including: the device comprises a geologic body resistivity detection device 51, a current emission source 52 and a magnetic total field detection device 53, wherein the geologic body resistivity detection device 51 is used for emitting and receiving ground surface and air, and the magnetic total field detection device 53 is mounted on an aircraft, and the current emission source 52 and the magnetic total field detection device 53 are connected with the geologic body resistivity detection device 51; in the embodiment of the invention, the current emission source and the magnetic total field detection device are synchronized through a GPS clock.
The geologic body resistivity detection device 51 is used for executing the geologic body resistivity detection method of ground transmission and air reception provided in the embodiments of the above methods;
the current emission source 52 is used for emitting a current signal with positive and negative bipolar square waveforms;
the total magnetic field detection device 53 is configured to detect a time series of total magnetic field data at each detection point in the region to be detected.
Fig. 6 illustrates a physical structure diagram of an electronic device, which may include, as shown in fig. 6: a processor (processor)610, a communication Interface (Communications Interface)620, a memory (memory)630 and a communication bus 640, wherein the processor 610, the communication Interface 620 and the memory 630 communicate with each other via the communication bus 640. The processor 610 may invoke logic instructions in the memory 630 to perform the method for surface-transmitted over-the-air received resistivity detection of a geologic volume provided in the various embodiments described above, the method comprising: acquiring a magnetic total field data time sequence at each detection point in the air after a current emission source is arranged in a region to be detected; for any detection point, determining a plurality of step response segments corresponding to the emission state in the time series of magnetic total field data at the any detection point based on the emission state of the current emission source; determining a magnetic field value at any probe point under the excitation condition of the current emission source based on the plurality of step response segments; and determining the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the detection points.
In addition, the logic instructions in the memory 630 may be implemented in software functional units and stored in a computer readable storage medium when the logic instructions are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, the present invention also provides a computer program product, the computer program product including a computer program, the computer program being stored on a non-transitory computer-readable storage medium, wherein when the computer program is executed by a processor, the computer is capable of executing the method for geological detection based on the magnetic resistivity method provided in the above embodiments, the method including: acquiring a magnetic total field data time sequence at each detection point in the air after a current emission source is arranged in a region to be detected; for any detection point, determining a plurality of step response segments corresponding to the emission state in the time series of magnetic total field data at the any detection point based on the emission state of the current emission source; determining a magnetic field value at any probe point under the excitation condition of the current emission source based on the plurality of step response segments; and determining the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the detection points.
In yet another aspect, the present invention also provides a non-transitory computer-readable storage medium having stored thereon a computer program which, when executed by a processor, is implemented to perform a method of magneto-resistivity based geological detection provided in the above embodiments, the method comprising: acquiring a magnetic total field data time sequence at each detection point in the air after a current emission source is arranged in a region to be detected; for any detection point, determining a plurality of step response segments corresponding to the emission state in the time series of magnetic total field data at the any detection point based on the emission state of the current emission source; determining a magnetic field value at any probe point under the excitation condition of the current emission source based on the plurality of step response segments; and determining the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the detection points.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing descriptions of specific exemplary embodiments of the present invention have been presented for purposes of illustration and description. It is not intended to limit the invention to the precise form disclosed, and obviously many modifications and variations are possible in light of the above teaching. The exemplary embodiments were chosen and described in order to explain certain principles of the invention and its practical application to enable one skilled in the art to make and use various exemplary embodiments of the invention and various alternatives and modifications as are suited to the particular use contemplated. It is intended that the scope of the invention be defined by the claims and their equivalents.

Claims (10)

1. A method for detecting the resistivity of a geologic body transmitted and received in the air by the ground is characterized by comprising the following steps:
acquiring a magnetic total field data time sequence at each detection point in the air after a current emission source is distributed in a region to be detected on the ground;
for any detection point, determining a plurality of step response segments corresponding to the emission state in the time series of magnetic total field data at the any detection point based on the emission state of the current emission source;
determining a magnetic field value at any probe point under the excitation condition of the current emission source based on the plurality of step response segments;
and determining the resistivity of the underground geologic body corresponding to the region to be detected based on the magnetic field values at the detection points.
2. The method according to claim 1, wherein the current emitter comprises a positive electrode, a negative electrode and a power source, and the positive electrode and the negative electrode are uniformly distributed in the area to be detected or in a preset range around the area to be detected;
the positive electrode is connected with the positive electrode of the power supply, and the negative electrode is connected with the negative electrode of the power supply;
the power supply is used for generating current signals of positive and negative bipolar square wave waveforms.
3. The method of claim 2, wherein the current emitting sources comprise a first current emitting source and a second current emitting source;
the positive electrode of the first current emission source is arranged at one side of the region to be detected, and the negative electrode of the first current emission source is arranged at a first position at one side far away from the region to be detected;
the negative electrode of the second current emission source is arranged on the other side of the region to be detected, and the positive electrode of the second current emission source is arranged at a second position far away from the other side of the region to be detected.
4. The method of claim 1, wherein the time series of magnetic total field data is detected by a magnetic total field detection device mounted on an aircraft;
correspondingly, the determining, based on the emission state of the current emission source, a plurality of step response segments corresponding to the emission state in the time series of total magnetic field data at any detection point specifically includes:
determining a current emission frequency of the current emission source and a flight speed of the aircraft;
determining a target magnetic total field data time sequence in the magnetic total field data time sequence at any detection point based on the current emission frequency and the flying speed;
based on the transmit state, a plurality of step response segments in the target total magnetic field data time series corresponding to the transmit state are determined.
5. The method of claim 1, wherein the transmitting comprises transmitting a first current signal having a positive square waveform and transmitting a second current signal having a negative square waveform;
correspondingly, the determining the magnetic field value at any detection point under the excitation condition of the current emission source based on the plurality of step response segments specifically includes:
negating the step response in the step response segment corresponding to the second current signal to obtain a target step response segment;
overlapping the step response segment corresponding to the first current signal with the step response in the target step response segment to obtain an average step response time sequence;
and calculating the average value of the step response in the time window of the target width in the average step response time sequence, and taking the average value as the magnetic field value.
6. The method for detecting resistivity of a geologic body received over the air by a ground based transmission according to any one of claims 1-5, wherein determining a plurality of step response segments in the time series of total magnetic field data at any one detection point corresponding to the transmission state based on the transmission state of the current transmitting source further comprises:
and filtering out fixed frequency noise and low frequency components of the magnetic total field data at each moment in the time sequence of the magnetic total field data at any detection point.
7. The method for detecting the resistivity of the geologic body received in the air by surface transmission according to any one of claims 1-5, wherein the determining the resistivity of the underground geologic body corresponding to the area to be detected based on the magnetic field values at each detection point further comprises:
determining a type of the geological volume based on the resistivity of the geological volume;
the types of the geologic body at least include a vein, a fracture zone, a karst cave, and groundwater.
8. A ground-based transmit-receive-over-the-air geologic resistivity survey system, comprising: the device comprises a geologic body resistivity detection device, a current emission source and a magnetic total field detection device, wherein the geologic body resistivity detection device is used for emitting and receiving in the air from the ground, and the magnetic total field detection device is mounted on an aircraft;
the geologic body resistivity detection device is used for executing the geologic body resistivity detection method of ground transmission and aerial reception according to any one of claims 1 to 7;
the current emission source is used for emitting a current signal with positive and negative bipolar square wave waveforms;
the magnetic total field detection device is used for detecting the magnetic total field data time sequence at each detection point in the air.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program performs the steps of the method for resistivity detection of a geological volume received over the air from a surface transmission according to any of claims 1 to 7.
10. A non-transitory computer readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements the steps of the method for surface-transmitted aerial-reception resistivity survey of geologic volume as claimed in any one of claims 1 to 7.
CN202111349047.5A 2021-11-15 2021-11-15 Ground emission air-received geological body resistivity detection method and system Active CN114002747B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111349047.5A CN114002747B (en) 2021-11-15 2021-11-15 Ground emission air-received geological body resistivity detection method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111349047.5A CN114002747B (en) 2021-11-15 2021-11-15 Ground emission air-received geological body resistivity detection method and system

Publications (2)

Publication Number Publication Date
CN114002747A true CN114002747A (en) 2022-02-01
CN114002747B CN114002747B (en) 2023-10-31

Family

ID=79929214

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111349047.5A Active CN114002747B (en) 2021-11-15 2021-11-15 Ground emission air-received geological body resistivity detection method and system

Country Status (1)

Country Link
CN (1) CN114002747B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100235100A1 (en) * 2009-03-16 2010-09-16 Bruce Alan Hobbs Method for determining resistivity anisotropy from earth electromagnetic responses
CN105425305A (en) * 2015-11-04 2016-03-23 中石化石油工程地球物理有限公司江汉分公司 Method for using time domain transient horizontal electric field down step response to improve resolution capability of electrical layer
CN111708099A (en) * 2020-06-17 2020-09-25 中国地质大学(武汉) Shallow earth surface geological comprehensive geophysical multi-parameter remote sensing composite system and acquisition method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100235100A1 (en) * 2009-03-16 2010-09-16 Bruce Alan Hobbs Method for determining resistivity anisotropy from earth electromagnetic responses
CN105425305A (en) * 2015-11-04 2016-03-23 中石化石油工程地球物理有限公司江汉分公司 Method for using time domain transient horizontal electric field down step response to improve resolution capability of electrical layer
CN111708099A (en) * 2020-06-17 2020-09-25 中国地质大学(武汉) Shallow earth surface geological comprehensive geophysical multi-parameter remote sensing composite system and acquisition method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘瑞德,梅岩辉,黄力军: "电性源瞬变电磁快速反演解释方法", 工程地球物理学报 *
袁桂琴;熊盛青;孟庆敏;周锡华;林品荣;王书民;高文利;徐明才;史大年;李秋生;: "地球物理勘查技术与应用研究", 地质学报 *
陈忠昌: "半航空磁电阻率法隧道透水灾害监测-以 7.15 珠海石景山隧道事故为例", 《中国地球科学联合学术年会2021》 *
陈忠昌等: "半航空磁电阻率法山区不良地质体探测方法研究", 《中国地球物理学会地球物理技术委员会第九届学术会议-全域地球物理探测与智能感知学术研讨会》 *

Also Published As

Publication number Publication date
CN114002747B (en) 2023-10-31

Similar Documents

Publication Publication Date Title
Palacky et al. Airborne electromagnetic methods
US10948620B2 (en) Electromagnetic receiver tracking and real-time calibration system and method
US20100026304A1 (en) Method and Apparatus for Analysing Geological Features
Becken et al. DESMEX: A novel system development for semi-airborne electromagnetic exploration
US11531134B2 (en) Multi-sensor system for airborne geophysical prospecting and method
AU2002328690A1 (en) Airborne geophysical measurements
EP1444536A1 (en) Airborne geophysical measurements
CN104062687A (en) Air ground integrated geomagnetic field combined observation method and system
Parvar Development and evaluation of unmanned aerial vehicle (UAV) magnetometry systems
WO2015085426A1 (en) Electromagnetic surveying at low frequencies using an airborne transmitter with receivers on the ground
US20180038983A1 (en) Earth surveying with aerial drones for improved drilling applications
CN111538106A (en) Unmanned archaeological exploration method suitable for ultra-shallow water area
AU2016203396B2 (en) Magnetometer signal sampling within time-domain EM transmitters and method
CN114002747B (en) Ground emission air-received geological body resistivity detection method and system
CN110146928A (en) Magnetic survey method and apparatus based on geological exploration unmanned plane
JP2018194454A (en) Underground electromagnetic survey device
Eröss et al. Three-component VLF using an unmanned aerial system as sensor platform
US2636924A (en) Method of geophysical exploration
FI127192B (en) SYSTEM FOR MEASURING THE ELECTRICAL PROPERTIES OF UNDERGROUND MATERIALS
CN219915951U (en) Receiver for electromagnetic acquisition system of ground-air unmanned aerial vehicle
Prospect Report on a Helicopter-Borne Versatile Time Domain Electromagnetic (VTEM Plus) and Horizontal Magnetic Gradiometer Geophysical Survey
CN114779352B (en) Geologic body survey system based on air-ground-well multi-dimensional electromagnetic method
CN116148933A (en) Receiver for electromagnetic acquisition system of ground-air unmanned aerial vehicle
Stoll et al. Airborne radioem using a remote piloted aircraft system
RU2649658C2 (en) Method and system of aeroelectric works

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
GR01 Patent grant
GR01 Patent grant