A geological exploration method, device and equipment based on coded source magnetic charging method

Abstract

The present invention relates to a geological exploration method, device, and equipment based on a coded source magnetic charging method, belonging to the field of geological exploration technology. The present invention establishes a coded current field by inputting a coded signal current of a sixth-order inverse-repeating M sequence into a good conductor underground in the area to be explored; measuring the three-component magnetic field generated by the good conductor on the ground after charging to obtain a three-component magnetic field signal; sampling the three-component magnetic field signal and the coded signal current based on a preset sampling rate to obtain a corresponding three-component magnetic field response time series and a coded signal current time series; identifying the coded signal current time series and the three-component magnetic field response time series using a dual-phase FFT cyclic cross-correlation method to obtain a plurality of magnetic field frequency responses corresponding to different frequencies; and then calculating the total magnetic field intensity and magnetic amplitude frequency to obtain the distribution characteristics of the good conductor in the underground medium space, thereby solving the problem of low detection accuracy due to electromagnetic interference.

Description

Technical Field

[0001] The present invention belongs to the field of geological exploration technology, and in particular relates to a geological exploration method, device and equipment based on a coding source magnetic charging method. Background Art

[0002] During geological exploration, the charging method is often used to observe the electric field distribution generated by an artificial current field established in the earth to detect underground media. The charging method uses the fact that good conductors, after being charged, behave as equipotential bodies or quasi-equipotential bodies. The potential or potential gradient of the charged equipotential bodies at the surface is observed to understand the distribution of the charged equipotential bodies underground, thereby completing geological exploration.

[0003] However, ground currents always propagate along the most conductive paths. Areas of good conductivity have high current density, generating high potentials, while areas of high resistance have low current density, generating low potentials. Because measuring the electric field requires a grounding electrode, observations in areas where grounding is difficult can be ineffective. Because traditional charging methods use direct current or low-frequency alternating current to establish the current field, geological exploration using these methods is susceptible to electromagnetic interference and has low anti-interference capabilities, resulting in poor detection accuracy and results. Summary of the Invention

[0004] To this end, the present invention provides a geological exploration method, device and equipment based on the coded source magnetic charging method, which helps to solve the problem of low detection accuracy caused by electromagnetic interference during geological exploration using traditional charging methods.

[0005] To achieve the above objectives, the present invention adopts the following technical solutions:

[0006] In a first aspect, the present invention provides a geological exploration method based on a coded source magnetic charging method, comprising:

[0007] Inputting a coded signal current into a good conductor underground in the area to be explored to establish a coded current field; wherein the signal waveform of the coded signal current is a sixth-order inverse repeating M sequence;

[0008] Measuring the three-component magnetic field generated on the ground by the good conductor after the coded signal current is input, and obtaining a corresponding three-component magnetic field signal; wherein the three-component magnetic field includes an X-component magnetic field, a Y-component magnetic field, and a Z-component magnetic field;

[0009] The three-component magnetic field signal and the coded signal current are respectively collected based on a preset sampling rate to obtain a corresponding three-component magnetic field response time series and a coded signal current time series;

[0010] The coded signal current time series and the three-component magnetic field response time series are identified by a two-phase FFT cyclic cross-correlation method to obtain a plurality of magnetic field frequency responses corresponding to different frequencies; wherein the magnetic field frequency responses include a magnetic field frequency response of the X-th component magnetic field, a magnetic field frequency response of the Y-th component magnetic field, and a magnetic field frequency response of the Z-th component magnetic field corresponding to the plurality of different frequencies;

[0011] Based on the plurality of magnetic field frequency responses, calculating the corresponding total magnetic field intensity and magnetic amplitude frequency;

[0012] Based on the total magnetic field strength and the magnetic amplitude frequency, the distribution characteristics of the good electrical conductor in the underground medium space are obtained.

[0013] Furthermore, the calculating of the corresponding total magnetic field strength based on the plurality of magnetic field frequency responses includes:

[0014] The total magnetic field intensity is calculated based on the magnetic field frequency response of the X-th component magnetic field, the magnetic field frequency response of the Y-th component magnetic field, the magnetic field frequency response of the Z-th component magnetic field, and a preset total magnetic field intensity calculation formula, wherein the preset total magnetic field intensity calculation formula is as follows:

[0015]

[0016] Among them, H total (f i ) refers to the frequency f i The corresponding total magnetic field strength, H x (f i ) refers to the frequency f i The magnetic field frequency response corresponding to the X-component magnetic field, H y (f i ) refers to the frequency f i The corresponding magnetic field frequency response of the Y component magnetic field, H z (f i ) refers to the frequency f i The corresponding magnetic field frequency response of the Z-th component magnetic field.

[0017] Furthermore, the calculating and obtaining corresponding magnetic amplitude frequencies based on the plurality of magnetic field frequency responses includes:

[0018] The magnetic amplitude frequency is calculated according to the magnetic field frequency response and a preset magnetic amplitude frequency calculation formula, wherein the preset magnetic amplitude frequency calculation formula is as follows:

[0019]

[0020] Among them, H kfis the magnetic amplitude frequency corresponding to the kth component magnetic field, where k is one of X, Y, and Z; f D is the low frequency, f G is a high frequency, and 3f D ≤f G ≤11f D ;H k (f D ) refers to the frequency f D The corresponding magnetic field frequency response of the k-th component magnetic field, H k (f G ) refers to the frequency f G The corresponding magnetic field frequency response of the k-th component magnetic field.

[0021] Furthermore, obtaining the distribution characteristics of the good conductor in the underground medium space based on the total magnetic field strength and the magnetic amplitude frequency includes:

[0022] Based on the total magnetic field strength and the magnetic amplitude frequency, the distribution characteristics of the good conductor in the underground medium space are obtained by performing profile measurement and area measurement;

[0023] The distribution characteristics of the good conductor in the underground medium space include at least one of the following: a total magnetic field intensity profile, a magnetic amplitude frequency profile, and a plane contour map.

[0024] In a second aspect, the present invention provides a geological exploration device based on a coded source magnetic charging method, comprising: a coded signal current transmitter, a three-component fluxgate magnetometer, a data acquisition recorder, and a data processing terminal;

[0025] The coded signal current transmitter is used to input the coded signal current into the good conductor underground in the area to be explored to establish a coded current field; wherein the signal waveform of the coded signal current is a sixth-order inverse repeating M sequence;

[0026] The three-component fluxgate magnetometer is used to measure the three-component magnetic field generated on the ground by the good conductor after the coded signal current is input, and obtain a corresponding three-component magnetic field signal; wherein the three-component magnetic field includes an X-component magnetic field, a Y-component magnetic field, and a Z-component magnetic field;

[0027] The data acquisition recorder is used to respectively acquire the three-component magnetic field signal and the coded signal current based on a preset sampling rate to obtain a corresponding three-component magnetic field response time series and a coded signal current time series;

[0028] The data processing terminal is used to identify the coded signal current time series and the three-component magnetic field response time series through a two-phase FFT cyclic cross-correlation method to obtain a plurality of magnetic field frequency responses corresponding to different frequencies; wherein the magnetic field frequency responses include a magnetic field frequency response of the X-th component magnetic field, a magnetic field frequency response of the Y-th component magnetic field, and a magnetic field frequency response of the Z-th component magnetic field corresponding to a plurality of different frequencies; based on the plurality of magnetic field frequency responses, a corresponding total magnetic field intensity and magnetic amplitude frequency are calculated; based on the total magnetic field intensity and the magnetic amplitude frequency, a distribution characteristic of the good conductor in the underground medium space is obtained.

[0029] Furthermore, the coded signal current transmitter includes: a transmission controller, a high-power inverter and a generator rectifier power supply;

[0030] The transmitting controller drives the high-power inverter to invert the DC current output by the generator rectifier power supply into the coded signal current, and outputs the coded signal current to the good conductor through a preset electrode pair to establish the coded current field;

[0031] Wherein, the preset electrode pair includes a charging electrode A and an infinite electrode B.

[0032] Furthermore, the data acquisition recorder is used to:

[0033] storing the three-component magnetic field response time series as a corresponding magnetic field data file; and

[0034] The encoded signal current time series is stored as a corresponding current data file.

[0035] Furthermore, the data acquisition recorder is also used for:

[0036] The magnetic field data file and the current data file are transmitted to the data processing terminal for processing via a preset transmission method; wherein the preset transmission method includes at least one of the following: a wireless method and a wired method.

[0037] Furthermore, the data acquisition recorder has at least three measurement channels with the same consistency.

[0038] In a third aspect, the present invention provides a geological exploration device based on a coded source magnetic charging method, comprising:

[0039] One or more memories storing executable programs;

[0040] One or more processors are used to execute the executable program in the memory to implement the steps of any of the above methods.

[0041] The present invention adopts the above technical solution and has at least the following beneficial effects:

[0042] The present invention establishes a coded current field by inputting a coded signal current of a sixth-order inverse-repeated M sequence into a good conductor underground in an area to be explored; then, the three-component magnetic field generated on the ground by the good conductor after the coded signal current is input is measured to obtain a corresponding three-component magnetic field signal; wherein the three-component magnetic field includes an X-component magnetic field, a Y-component magnetic field, and a Z-component magnetic field; based on a preset sampling rate, the three-component magnetic field signal and the coded signal current are respectively collected to obtain a corresponding three-component magnetic field response time series and a coded signal current time series; the coded signal current time series and the three-component magnetic field response time series are identified by a two-phase FFT cyclic cross-correlation method to obtain magnetic field frequency responses corresponding to multiple different frequencies; based on the multiple magnetic field frequency responses, the corresponding total magnetic field strength and magnetic amplitude frequency are calculated; based on the total magnetic field strength and magnetic amplitude frequency, the distribution characteristics of the good conductor in the underground medium space are obtained, and anomaly analysis and geological interpretation of the distribution characteristics of the good conductor in the underground medium space are performed. By adopting the dual-phase FFT cyclic cross-correlation method, the coded signal current and the three-component magnetic field signal are identified, and the magnetic field frequency responses corresponding to multiple different frequencies are obtained. It can have a strong interference suppression capability against various types of random noise, thereby solving the problem of low detection accuracy caused by electromagnetic interference during geological exploration using traditional charging methods.

[0043] It is to be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0045] Figure 1 This is a flow chart of a geological exploration method based on a coded source magnetic charging method according to an embodiment of the present invention;

[0046] Figure 2 This is a schematic block diagram of a geological exploration device based on a coded source magnetic charging method according to an embodiment of the present invention;

[0047] Figure 3 is a block diagram of a coded signal current transmitter according to an embodiment of the present invention;

[0048] Figure 4This is a schematic block diagram of a geological exploration device based on a coded source magnetic charging method according to an embodiment of the present invention;

[0049] Figure 5 The present invention is a flowchart showing a data processing method of a geological exploration method based on a coded source magnetic charging method according to an embodiment of the present invention. DETAILED DESCRIPTION

[0050] To make the objectives, technical solutions, and advantages of the present invention more apparent, the technical solutions of the present invention will be described in detail below. Obviously, the embodiments described are only some of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other implementations obtained by those of ordinary skill in the art without inventive effort are within the scope of protection of the present invention.

[0051] See also Figure 1 , Figure 1 FIG. 1 is a flow chart of a geological exploration method based on a coding source magnetic charging method according to an embodiment of the present invention. Figure 1 As shown, the method includes the following steps:

[0052] Step S11: inputting a coded signal current into a good conductor underground in the area to be explored to establish a coded current field; wherein the signal waveform of the coded signal current is a sixth-order inverse repeating M sequence;

[0053] Step S12: measuring the three-component magnetic field generated on the ground by the good conductor after the coded signal current is input, and obtaining a corresponding three-component magnetic field signal; wherein the three-component magnetic field includes an X-component magnetic field, a Y-component magnetic field, and a Z-component magnetic field;

[0054] Step S13: collecting the three-component magnetic field signals and the coded signal current based on a preset sampling rate to obtain corresponding three-component magnetic field response time series and coded signal current time series;

[0055] Step S14: Identify the coded signal current time series and the three-component magnetic field response time series using a dual-phase FFT cyclic cross-correlation method to obtain a plurality of magnetic field frequency responses corresponding to different frequencies; wherein the magnetic field frequency responses include a magnetic field frequency response of the X-th component magnetic field, a magnetic field frequency response of the Y-th component magnetic field, and a magnetic field frequency response of the Z-th component magnetic field corresponding to the plurality of different frequencies;

[0056] Step S15: Based on the plurality of magnetic field frequency responses, calculating and obtaining the corresponding total magnetic field intensity and magnetic amplitude frequency;

[0057] Step S16: Based on the total magnetic field intensity and the magnetic amplitude frequency, obtain the distribution characteristics of the good conductor in the underground medium space.

[0058] It should be noted that the charging method is a geophysical electrical exploration method. When underground metal ore bodies have exposed outcrops, or underground rivers, underground pipelines, etc. have outcrops near the surface, we can establish an electric current field by charging the outcrops; since metal ore bodies, underground rivers, underground pipelines and other good conductors will behave as equipotential bodies after charging, they can generate a strong electric field on the surface. By observing the potential or potential gradient generated by this equipotential body on the surface, the distribution pattern of the potential and potential gradient can be used to infer the distribution characteristics of good conductors underground, so as to achieve the purpose of exploring underground geological bodies.

[0059] Specifically, in this embodiment, a coded signal current of a sixth-order inverse-repeating M sequence is used to excite the current field, charging the good conductor with a low-frequency coded current. At the same time, the magnetic field generated after the good conductor behaves as an equipotential body is observed on the ground, thereby inferring the spatial existence state of this good conductor underground.

[0060] The coding signal current is a sixth-order inverse-repeating M-sequence pseudo-random signal, which is output by a coding signal current transmitter. The period of the coding signal current is 6s, and the signal frequency range of the electromagnetic field is 1 / 6Hz to 10Hz.

[0061] After the coded signal current is input to the good conductor underground in the area to be explored, the good conductor is charged and becomes an equipotential body. By measuring the magnetic field generated by the equipotential body on the ground, the corresponding magnetic field signal is obtained. Since the magnetic field includes three components (i.e., H x , H y , H z ), therefore, the magnetic field is also called a three-component magnetic field.

[0062] The three-component magnetic field response time series and the coded signal current time series are obtained by collecting data of the three-component magnetic field signals and the coded signal current using a data acquisition and recorder at a preset sampling rate. The preset sampling rate is a pre-set rate value and is not specifically limited in this embodiment.

[0063] The FFT (Fast Fourier Transform) is an efficient algorithm for the DFT, also known as the Fast Fourier Transform. The Fourier transform is one of the most fundamental methods for time-to-frequency domain analysis. In digital processing applications, the DFT (Discrete Fourier Transform) forms the foundation of many digital signal processing methods.

[0064] Specifically, in this embodiment, coded signal currents are input to multiple locations of good conductors underground in the area to be explored, and three-component magnetic field signals are measured at multiple locations on the ground in the area to be explored. That is, there will be multiple coded signal currents and multiple three-component magnetic field signals at multiple test points. Figure 5 As shown, Figure 5 This is a flow chart illustrating data processing for a geological exploration method based on the coded source magnetic charging method, as illustrated in an embodiment of the present invention. After recording and collecting the three-component magnetic field response time series and the coded signal current time series at each measuring point, the frequency responses corresponding to the X, Y, and Z components of the magnetic field at each measuring point are identified using the dual-phase FFT cyclic cross-correlation method. In other words, multiple frequency responses of different frequencies are identified, namely, the magnetic field frequency responses of the X-th component magnetic field, the Y-th component magnetic field, and the Z-th component magnetic field corresponding to multiple frequencies. Based on the magnetic field frequency responses, the total magnetic field intensity and magnetic amplitude frequency at a certain frequency at each measuring point can be calculated. Parameters such as the total magnetic field intensity and magnetic amplitude frequency can then be used to indicate the spatial distribution of underground good conductors, thereby achieving the purpose of exploring underground geological bodies. Furthermore, the dual-phase FFT cyclic cross-correlation method is used to identify the frequency responses corresponding to each magnetic field component by performing cyclic cross-correlation on the recorded input coded signal current and the measured three-component magnetic field signal. Multi-cycle measurements can provide a strong interference suppression capability against various types of random noise.

[0065] It should be noted that when the coded signal current is used to continuously excite the low-frequency coded current field for multiple cycles, a wide-band three-component magnetic field frequency response can be obtained. At this time, the obtained frequencies are richer and more interpretation parameters can be obtained.

[0066] It can be understood that the present invention establishes a coded current field by inputting a coded signal current of a sixth-order inverse-repeated M sequence into a good conductor underground in the area to be explored; then, the three-component magnetic field generated by the good conductor on the ground after the coded signal current is input is measured to obtain a corresponding three-component magnetic field signal; wherein the three-component magnetic field includes an X-component magnetic field, a Y-component magnetic field, and a Z-component magnetic field; based on a preset sampling rate, the three-component magnetic field signal and the coded signal current are respectively collected to obtain a corresponding three-component magnetic field response time series and a coded signal current time series; the coded signal current time series and the three-component magnetic field response time series are identified by a two-phase FFT cyclic cross-correlation method to obtain a plurality of magnetic field frequency responses corresponding to different frequencies; based on the plurality of magnetic field frequency responses, the corresponding total magnetic field strength and magnetic amplitude frequency are calculated; according to the total magnetic field strength and magnetic amplitude frequency, the distribution characteristics of the good conductor in the underground medium space are obtained, and anomaly analysis and geological interpretation are performed on the distribution characteristics of the good conductor in the underground medium space. By adopting the dual-phase FFT cyclic cross-correlation method, the coded signal current and the three-component magnetic field signal are identified, and the magnetic field frequency responses corresponding to multiple different frequencies are obtained. It can have a strong interference suppression capability against various types of random noise, thereby solving the problem of low detection accuracy caused by electromagnetic interference during geological exploration using traditional charging methods.

[0067] Furthermore, the calculating of the corresponding total magnetic field strength based on the plurality of magnetic field frequency responses includes:

[0068] The total magnetic field intensity is calculated based on the magnetic field frequency response of the X-th component magnetic field, the magnetic field frequency response of the Y-th component magnetic field, the magnetic field frequency response of the Z-th component magnetic field, and a preset total magnetic field intensity calculation formula, wherein the preset total magnetic field intensity calculation formula is as follows:

[0069]

[0070] Among them, H total (f i ) refers to the frequency f i The corresponding total magnetic field strength, H x (f i ) refers to the frequency f i The magnetic field frequency response corresponding to the X-component magnetic field, H y (f i ) refers to the frequency f i The corresponding magnetic field frequency response of the Y component magnetic field, H z (f i ) refers to the frequency f i The corresponding magnetic field frequency response of the Z-th component magnetic field.

[0071] It should be noted that in this embodiment, coded signal currents will be input to multiple locations of good conductors underground in the area to be explored, and three-component magnetic field signals will also be measured at multiple locations on the ground in the area to be explored, that is, there will be multiple coded signal currents and multiple three-component magnetic field signals at multiple test points. After recording and collecting the three-component magnetic field response time series and coded signal current time series at each measuring point, the frequency response corresponding to the X, Y, and Z components of the magnetic field at each test point can be identified by the dual-phase FFT cyclic cross-correlation method. By using the frequency f i The corresponding magnetic field frequency response of the X-th component magnetic field, the magnetic field frequency response of the Y-th component magnetic field, the magnetic field frequency response of the Z-th component magnetic field and the preset total magnetic field intensity calculation formula can be used to calculate the frequency f i The corresponding total magnetic field strength.

[0072] Furthermore, the calculating and obtaining corresponding magnetic amplitude frequencies based on the plurality of magnetic field frequency responses includes:

[0073] The magnetic amplitude frequency is calculated according to the magnetic field frequency response and a preset magnetic amplitude frequency calculation formula, wherein the preset magnetic amplitude frequency calculation formula is as follows:

[0074]

[0075] Among them, H kf is the magnetic amplitude frequency corresponding to the kth component magnetic field, where k is one of X, Y, and Z; f D is the low frequency, f G is a high frequency, and 3f D ≤f G ≤11f D ;H k (f D ) refers to the frequency f D The corresponding magnetic field frequency response of the k-th component magnetic field, H k (f G ) refers to the frequency f G The corresponding magnetic field frequency response of the k-th component magnetic field.

[0076] It should be noted that the magnetic amplitude frequency can also be calculated based on the magnetic field frequency response and a preset magnetic amplitude frequency calculation formula. That is, at each test point, each component magnetic field can identify a magnetic field frequency response at several frequencies. Two frequencies are arbitrarily selected from these frequencies to calculate the magnetic amplitude frequency: a low frequency and a high frequency. The frequency difference between these two frequencies can be between 2 and 10 times, and the specific values ​​of the low frequency and high frequency are not limited.

[0077] Furthermore, obtaining the distribution characteristics of the good conductor in the underground medium space based on the total magnetic field strength and the magnetic amplitude frequency includes:

[0078] Based on the total magnetic field strength and the magnetic amplitude frequency, the distribution characteristics of the good conductor in the underground medium space are obtained by performing profile measurement and area measurement;

[0079] The distribution characteristics of the good conductor in the underground medium space include at least one of the following: a total magnetic field intensity profile, a magnetic amplitude frequency profile, and a plane contour map.

[0080] It should be noted that at each test point, each component magnetic field can be identified as having a frequency response at several frequencies. Two frequencies are arbitrarily selected from these frequencies to calculate the magnetic amplitude frequency, and each frequency corresponds to a total magnetic field intensity. Through profile and area measurements, the distribution characteristics of the total magnetic field intensity, magnetic amplitude frequency profiles, and plane contour maps of good conductors in the underground medium space are drawn based on the obtained total magnetic field at multiple frequencies and multiple magnetic amplitude frequencies, providing more information for geological exploration interpretation.

[0081] See also Figure 2 , Figure 2 1 is a block diagram of a geological exploration device based on a coded source magnetic charging method according to an embodiment of the present invention. The geological exploration device 2 based on the coded source magnetic charging method may include: a coded signal current transmitter 21, a three-component fluxgate magnetometer 22, a data acquisition recorder 23, and a data processing terminal 24.

[0082] The coded signal current transmitter 21 is used to input the coded signal current into the good conductor underground in the exploration area to establish a coded current field; wherein the signal waveform of the coded signal current is a sixth-order inverse repeating M sequence;

[0083] The three-component fluxgate magnetometer 22 is used to measure the three-component magnetic field generated on the ground by the good conductor after the coded signal current is input, and obtain a corresponding three-component magnetic field signal; wherein the three-component magnetic field includes an X-component magnetic field, a Y-component magnetic field, and a Z-component magnetic field;

[0084] The data acquisition recorder 23 is used to respectively acquire the three-component magnetic field signal and the coded signal current based on a preset sampling rate to obtain a corresponding three-component magnetic field response time series and a coded signal current time series;

[0085] The data processing terminal 24 is used to identify the coded signal current time series and the three-component magnetic field response time series through a two-phase FFT cyclic cross-correlation method to obtain a plurality of magnetic field frequency responses corresponding to different frequencies; wherein the magnetic field frequency responses include a magnetic field frequency response of the X-th component magnetic field, a magnetic field frequency response of the Y-th component magnetic field, and a magnetic field frequency response of the Z-th component magnetic field corresponding to a plurality of different frequencies; based on the plurality of magnetic field frequency responses, the corresponding total magnetic field intensity and magnetic amplitude frequency are calculated; based on the total magnetic field intensity and the magnetic amplitude frequency, the distribution characteristics of the good conductor in the underground medium space are obtained.

[0086] It should be noted that the coded signal current transmitter 21 can invert direct current into a sixth-order inverse-repeating M-sequence pseudo-random signal, and transmit the sixth-order inverse-repeating M-sequence coded signal current to the good conductor underground in the area to be explored through an electrode pair to establish a coded current field.

[0087] The three-component fluxgate magnetometer 22 measures the three-component magnetic field generated by a good conductor (i.e., an equipotential body) on the ground after inputting the coded signal current. In actual geological exploration, one or more three-component fluxgate magnetometers can be used in an array configuration for observation. This array-based observation significantly improves accuracy and efficiency.

[0088] The data acquisition recorder 23 can be composed of a multi-channel signal conditioning module, a multi-channel high-speed A / D converter, an FPGA and an ARM processor, and can collect the three-component magnetic field signal output by the three-component fluxgate magnetometer 22 as a three-component magnetic field response time series. It can also record the coded signal current output by the coded signal current transmitter 21 as a coded signal current time series using the same preset sampling rate.

[0089] The data processing terminal 24 can use the dual-phase FFT cyclic cross-correlation method to identify the magnetic field frequency response corresponding to different frequencies at each measuring point using the synchronously collected coded signal current time series and three-component magnetic field response time series. Furthermore, the data processing terminal 24 can calculate the total magnetic field intensity and magnetic amplitude frequency for multiple frequencies at multiple measuring points based on preset total magnetic field intensity and magnetic amplitude frequency calculation formulas, and can also map the distribution characteristics of good conductors in the underground medium space.

[0090] It can be understood that the present invention uses a coded signal current transmitter 21 to input a coded signal current of a sixth-order inverse-repeating M sequence into a good conductor underground in the exploration area to establish a coded current field. Then, a three-component fluxgate magnetometer 22 measures the three-component magnetic field generated by the good conductor on the ground after the coded signal current is input, thereby obtaining a corresponding three-component magnetic field signal. The three-component magnetic field includes an X-component magnetic field, a Y-component magnetic field, and a Z-component magnetic field. A data acquisition recorder 23 collects the three-component magnetic field signal and the coded signal current at a preset sampling rate to obtain a corresponding three-component magnetic field response time series and a coded signal current time series. Finally, a data processing terminal 24 identifies the coded signal current time series and the three-component magnetic field response time series using a two-phase FFT cyclic cross-correlation method to obtain magnetic field frequency responses corresponding to multiple different frequencies. Based on the multiple magnetic field frequency responses, the corresponding total magnetic field strength and magnetic amplitude frequency are calculated. Based on the total magnetic field strength and magnetic amplitude frequency, the distribution characteristics of the good conductor in the underground medium space are obtained, and anomaly analysis and geological interpretation of the distribution characteristics of the good conductor in the underground medium space are performed. By adopting the dual-phase FFT cyclic cross-correlation method, the coded signal current and the three-component magnetic field signal are identified, and the magnetic field frequency responses corresponding to multiple different frequencies are obtained. It can have a strong interference suppression capability against various types of random noise, thereby solving the problem of low detection accuracy caused by electromagnetic interference during geological exploration using traditional charging methods.

[0091] Furthermore, the coded signal current transmitter 21 includes: a transmission controller 211, a high-power inverter 212 and a generator rectifier power supply 213;

[0092] The transmitting controller 211 drives the high-power inverter 212 to invert the DC current output by the generator rectifier power supply 213 into the coded signal current, and outputs the coded signal current to the good conductor through a preset electrode pair to establish the coded current field;

[0093] Wherein, the preset electrode pair includes a charging electrode A and an infinite electrode B.

[0094] It should be noted that if Figure 3 As shown, Figure 3 2 is a block diagram of a coded signal current transmitter according to an embodiment of the present invention. The coded signal current transmitter 21 includes a transmission controller 211 , a high-power inverter 212 and a generator rectifier power supply 213 .

[0095] The transmitting controller 211 divides the 7812000Hz clock crystal oscillator to obtain a 21Hz clock sequence, which stimulates the output of a 126-bit sixth-order inverse-repeating M-sequence. It also controls the high-power inverter 212 to modulate the high-voltage direct current to generate a sixth-order inverse-repeating M-sequence pseudo-random signal current. The generated coded signal current has a period of 6s, and the signal frequency range of the electromagnetic field is 1 / 6Hz to 10Hz.

[0096] Furthermore, the data acquisition recorder 23 is used to:

[0097] storing the three-component magnetic field response time series as a corresponding magnetic field data file; and

[0098] The encoded signal current time series is stored as a corresponding current data file.

[0099] Furthermore, the data acquisition recorder 23 is also used for:

[0100] The magnetic field data file and the current data file are transmitted to the data processing terminal for processing via a preset transmission method; wherein the preset transmission method includes at least one of the following: a wireless method and a wired method.

[0101] It should be noted that the data acquisition recorder 23 can collect the three-component magnetic field signals output by the three-component fluxgate magnetometer 22 as a three-component magnetic field response time series and store it as a magnetic field data file. It can also record the coded signal current output by the coded signal current transmitter 21 as a coded signal current time series using the same preset sampling rate and store it as a current data file. The magnetic field data file and the current data file can also be transmitted wirelessly or wired to a data processing terminal for processing.

[0102] Furthermore, the data acquisition recorder 23 has at least three measurement channels with the same consistency.

[0103] It should be noted that the data acquisition recorder 23 needs to collect the three-component magnetic field signal output by the three-component fluxgate magnetometer 22 as a three-component magnetic field response time series. Moreover, the three-component magnetic field includes the X-component magnetic field, the Y-component magnetic field, and the Z-component magnetic field, and the magnetic field has directionality. When collecting magnetic field signals, it is necessary to collect the three components of the magnetic field simultaneously. In order to ensure that the magnetic field strength of the three components under the same conditions is collected, it is necessary to ensure that the measurement channels are identical. In other words, at least three measurement channels with the same consistency are required to receive the magnetic field signals.

[0104] Furthermore, the data processing terminal 24 includes a total magnetic field strength calculation module, which is configured to calculate the corresponding total magnetic field strength based on the plurality of magnetic field frequency responses, including:

[0105] The total magnetic field intensity is calculated based on the magnetic field frequency response of the X-th component magnetic field, the magnetic field frequency response of the Y-th component magnetic field, the magnetic field frequency response of the Z-th component magnetic field, and a preset total magnetic field intensity calculation formula, wherein the preset total magnetic field intensity calculation formula is as follows:

[0106]

[0107] Among them, H totak (f i ) refers to the frequency f i The corresponding total magnetic field strength, H x (f i ) refers to the frequency f i The magnetic field frequency response corresponding to the X-component magnetic field, H y (f i ) refers to the frequency f i The corresponding magnetic field frequency response of the Y component magnetic field, H z (f i ) refers to the frequency f i The corresponding magnetic field frequency response of the Z-th component magnetic field.

[0108] Furthermore, the data processing terminal 24 includes a magnetic amplitude frequency calculation module, which is configured to calculate the corresponding magnetic amplitude frequency based on the plurality of magnetic field frequency responses, including:

[0109] The magnetic amplitude frequency is calculated according to the magnetic field frequency response and a preset magnetic amplitude frequency calculation formula, wherein the preset magnetic amplitude frequency calculation formula is as follows:

[0110]

[0111] Among them, H kf is the magnetic amplitude frequency corresponding to the kth component magnetic field, where k is one of X, Y, and Z; f D is the low frequency, f G is a high frequency, and 3f D ≤f G ≤11f D ;H k (f D ) refers to the frequency f D The corresponding magnetic field frequency response of the k-th component magnetic field, H k (fG ) refers to the frequency f G The corresponding magnetic field frequency response of the k-th component magnetic field.

[0112] Furthermore, the data processing terminal 24 includes a distribution feature generation module, which is used to obtain the distribution characteristics of the good conductor in the underground medium space based on the total magnetic field intensity and the magnetic amplitude frequency, including:

[0113] Based on the total magnetic field strength and the magnetic amplitude frequency, the distribution characteristics of the good conductor in the underground medium space are obtained by performing profile measurement and area measurement;

[0114] The distribution characteristics of the good conductor in the underground medium space include at least one of the following: a total magnetic field intensity profile, a magnetic amplitude frequency profile, and a plane contour map.

[0115] It should be noted that in the present invention, the geological exploration device 2 based on the coded source magnetic charging method does not require electrode grounding, can be set on the ground to measure the magnetic field, and can also be carried on a multi-rotor UAV platform for semi-aerial flight observation, thereby improving work efficiency.

[0116] It should be noted that the present invention can be used for metal mine exploration, as well as for groundwater flow monitoring and underground pipeline detection, which is not specifically limited in this embodiment.

[0117] Regarding the geological exploration device 2 based on the coded source magnetic charging method in the above embodiment, the specific manner in which each device performs operations has been described in detail in the embodiments of the above related methods and will not be elaborated on here.

[0118] See also Figure 4 , Figure 4 1 is a block diagram of a geological exploration device based on a coding source magnetic charging method according to an embodiment of the present invention. The geological exploration device 3 based on the coding source magnetic charging method includes:

[0119] One or more memories 31 on which executable programs are stored;

[0120] One or more processors 32 are configured to execute the executable program in the memory 31 to implement the steps of any of the above methods.

[0121] Regarding the geological exploration equipment 3 based on the coded source magnetic charging method in the above embodiment, the specific manner in which its processor 32 executes the program in the memory 31 has been described in detail in the embodiment of the method and will not be elaborated here.

[0122] It can be understood that the same or similar parts of the above embodiments can be referenced to each other, and the contents not described in detail in some embodiments can refer to the same or similar contents in other embodiments.

[0123] It should be noted that, in the description of the present invention, the terms "first," "second," etc. are used for descriptive purposes only and should not be understood as indicating or implying relative importance. In addition, in the description of the present invention, unless otherwise specified, the meaning of "plurality" or "multiple" is at least two.

[0124] It should be understood that when an element is referred to as being "fixed to" or "disposed on" another element, it can be directly on the other element or there may be an intermediate element at the same time; when an element is referred to as being "connected to" another element, it can be directly connected to the other element or there may be an intermediate element at the same time. In addition, the "connection" used here may include wireless connection; the wording "and / or" used includes any unit and all combinations of one or more associated listed items.

[0125] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, segment or portion of code comprising one or more executable instructions for implementing the steps of a specific logical function or process, and the scope of the preferred embodiments of the present invention includes alternative implementations in which functions may be performed out of the order shown or discussed, including performing functions in a substantially simultaneous manner or in the reverse order depending on the functions involved, which should be understood by those skilled in the art to which the embodiments of the present invention pertain.

[0126] It should be understood that various parts of the present invention can be implemented using hardware, software, firmware, or a combination thereof. In the above-described embodiments, multiple steps or methods can be implemented using software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented using hardware, as in another embodiment, any one of the following technologies known in the art or a combination thereof can be used: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application-specific integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0127] Those skilled in the art will understand that all or part of the steps in the method of the above embodiment can be completed by instructing related hardware through a program, and the program can be stored in a computer-readable storage medium. When the program is executed, it includes one or a combination of the steps of the method embodiment.

[0128] In addition, the functional units in the various embodiments of the present invention may be integrated into a single processing module, or each unit may exist physically separately, or two or more units may be integrated into a single module. The aforementioned integrated modules may be implemented in the form of hardware or in the form of software functional modules. If the integrated modules are implemented in the form of software functional modules and sold or used as independent products, they may also be stored in a computer-readable storage medium.

[0129] The storage medium mentioned above can be a read-only memory, a magnetic disk or an optical disk, etc.

[0130] Throughout this specification, reference to terms such as "one embodiment," "some embodiments," "examples," "specific examples," or "some examples" means that a specific feature, structure, material, or characteristic described in conjunction with that embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

[0131] Although the embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are illustrative and are not to be construed as limitations on the present invention. A person skilled in the art may change, modify, replace and modify the above embodiments within the scope of the present invention.

Claims

1. A geological exploration method based on a coded source magnetic charging method, characterized in that: include: Inputting a coded signal current into a good conductor underground in the area to be explored to establish a coded current field; wherein the signal waveform of the coded signal current is a sixth-order inverse repeating M sequence; Measuring the three-component magnetic field generated on the ground by the good conductor after the coded signal current is input, and obtaining a corresponding three-component magnetic field signal; wherein the three-component magnetic field includes an X-component magnetic field, a Y-component magnetic field, and a Z-component magnetic field; The three-component magnetic field signal and the coded signal current are respectively collected based on a preset sampling rate to obtain a corresponding three-component magnetic field response time series and a coded signal current time series; The coded signal current time series and the three-component magnetic field response time series are identified by a two-phase FFT cyclic cross-correlation method to obtain a plurality of magnetic field frequency responses corresponding to different frequencies; wherein the magnetic field frequency responses include a magnetic field frequency response of the X-th component magnetic field, a magnetic field frequency response of the Y-th component magnetic field, and a magnetic field frequency response of the Z-th component magnetic field corresponding to the plurality of different frequencies; Based on the plurality of magnetic field frequency responses, calculating the corresponding total magnetic field intensity and magnetic amplitude frequency; Based on the total magnetic field strength and the magnetic amplitude frequency, the distribution characteristics of the good electrical conductor in the underground medium space are obtained.

2. The method according to claim 1, characterized in that The calculating and obtaining the corresponding total magnetic field strength based on the plurality of magnetic field frequency responses includes: The total magnetic field intensity is calculated based on the magnetic field frequency response of the X-th component magnetic field, the magnetic field frequency response of the Y-th component magnetic field, the magnetic field frequency response of the Z-th component magnetic field, and a preset total magnetic field intensity calculation formula, wherein the preset total magnetic field intensity calculation formula is as follows: Among them, H total (f i ) refers to the frequency f i The corresponding total magnetic field strength, H x (f i ) refers to the frequency f i The magnetic field frequency response corresponding to the X-component magnetic field, H y (f i ) refers to the frequency f i The corresponding magnetic field frequency response of the Y component magnetic field, H z (f i ) refers to the frequency f i The corresponding magnetic field frequency response of the Z-th component magnetic field.

3. The method according to claim 2, characterized in that The calculating and obtaining corresponding magnetic amplitude frequencies based on the plurality of magnetic field frequency responses includes: The magnetic amplitude frequency is calculated according to the magnetic field frequency response and a preset magnetic amplitude frequency calculation formula, wherein the preset magnetic amplitude frequency calculation formula is as follows: Among them, H kf is the magnetic amplitude frequency corresponding to the kth component magnetic field, where k is one of X, Y, and Z; f D is the low frequency, f G is a high frequency, and 3f D ≤f G ≤11f D ;H k (f D ) refers to the frequency f D The corresponding magnetic field frequency response of the k-th component magnetic field, H k (f G ) refers to the frequency f G The corresponding magnetic field frequency response of the k-th component magnetic field.

4. The method according to claim 3, characterized in that The obtaining of the distribution characteristics of the good conductor in the underground medium space based on the total magnetic field strength and the magnetic amplitude frequency includes: Based on the total magnetic field strength and the magnetic amplitude frequency, the distribution characteristics of the good conductor in the underground medium space are obtained by performing profile measurement and area measurement; The distribution characteristics of the good conductor in the underground medium space include at least one of the following: a total magnetic field intensity profile, a magnetic amplitude frequency profile, and a plane contour map.

5. A geological exploration device based on the coded source magnetic charging method, characterized in that: include: Coded signal current transmitter, three-component fluxgate magnetometer, data acquisition recorder and data processing terminal; The coded signal current transmitter is used to input the coded signal current into the good conductor underground in the area to be explored to establish a coded current field; wherein the signal waveform of the coded signal current is a sixth-order inverse repeating M sequence; The three-component fluxgate magnetometer is used to measure the three-component magnetic field generated on the ground by the good conductor after the coded signal current is input, and obtain a corresponding three-component magnetic field signal; wherein the three-component magnetic field includes an X-component magnetic field, a Y-component magnetic field, and a Z-component magnetic field; The data acquisition recorder is used to respectively acquire the three-component magnetic field signal and the coded signal current based on a preset sampling rate to obtain a corresponding three-component magnetic field response time series and a coded signal current time series; The data processing terminal is used to identify the coded signal current time series and the three-component magnetic field response time series through a two-phase FFT cyclic cross-correlation method to obtain a plurality of magnetic field frequency responses corresponding to different frequencies; wherein the magnetic field frequency responses include a magnetic field frequency response of the X-th component magnetic field, a magnetic field frequency response of the Y-th component magnetic field, and a magnetic field frequency response of the Z-th component magnetic field corresponding to a plurality of different frequencies; based on the plurality of magnetic field frequency responses, a corresponding total magnetic field intensity and magnetic amplitude frequency are calculated; based on the total magnetic field intensity and the magnetic amplitude frequency, a distribution characteristic of the good conductor in the underground medium space is obtained.

6. The device according to claim 5, characterized in that The coded signal current transmitter includes: a transmission controller, a high-power inverter and a generator rectifier power supply; The transmitting controller drives the high-power inverter to invert the DC current output by the generator rectifier power supply into the coded signal current, and outputs the coded signal current to the good conductor through a preset electrode pair to establish the coded current field; Wherein, the preset electrode pair includes a charging electrode A and an infinite electrode B.

7. The device according to claim 5, characterized in that The data acquisition recorder is used to: storing the three-component magnetic field response time series as a corresponding magnetic field data file; and The encoded signal current time series is stored as a corresponding current data file.

8. The device according to claim 7, characterized in that The data acquisition recorder is also used for: The magnetic field data file and the current data file are transmitted to the data processing terminal for processing via a preset transmission method; wherein the preset transmission method includes at least one of the following: a wireless method and a wired method.

9. The device according to claim 8, characterized in that The data acquisition recorder has at least three measurement channels with the same consistency.

10. A geological exploration device based on the coded source magnetic charging method, characterized in that: include: One or more memories storing executable programs; One or more processors, configured to execute the executable program in the memory to implement the steps of the method according to any one of claims 1 to 4.

Images