CN111123369B - Geological exploration wave detection method, device, equipment and medium - Google Patents

Geological exploration wave detection method, device, equipment and medium Download PDF

Info

Publication number
CN111123369B
CN111123369B CN202010009649.5A CN202010009649A CN111123369B CN 111123369 B CN111123369 B CN 111123369B CN 202010009649 A CN202010009649 A CN 202010009649A CN 111123369 B CN111123369 B CN 111123369B
Authority
CN
China
Prior art keywords
electromagnetic
signal
signals
geological exploration
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010009649.5A
Other languages
Chinese (zh)
Other versions
CN111123369A (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.)
Hunan Non-Ferrous Metals Geology Investigation Bureau Rows 247
Original Assignee
Hunan Non-Ferrous Metals Geology Investigation Bureau Rows 247
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 Hunan Non-Ferrous Metals Geology Investigation Bureau Rows 247 filed Critical Hunan Non-Ferrous Metals Geology Investigation Bureau Rows 247
Priority to CN202010009649.5A priority Critical patent/CN111123369B/en
Publication of CN111123369A publication Critical patent/CN111123369A/en
Application granted granted Critical
Publication of CN111123369B publication Critical patent/CN111123369B/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/12Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
    • 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/38Processing data, e.g. for analysis, for interpretation, for correction

Abstract

The application relates to a geological exploration wave detection method, a device, equipment and a medium, wherein the method comprises the following steps: acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals; acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values; acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range; and acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal. According to the geological exploration electromagnetic detection signal processing method and device, the characteristic parameter value of the received geological exploration electromagnetic detection signal is compared with the predicted threshold range, the noise signal in the received geological exploration electromagnetic detection signal is filtered out, the geological exploration electromagnetic detection signal capable of reflecting the real geological structure is obtained, and the accuracy of geological exploration by using the electromagnetic signal is improved.

Description

Geological exploration wave detection method, device, equipment and medium
Technical Field
The application relates to the technical field of geological exploration, in particular to a geological exploration detection method, a device, equipment and a medium.
Background
The method has the advantages that the geology is surveyed and detected through various means and methods, mineral deposits with industrial significance are found in mineral product general survey, the quality and quantity of the mineral products are found, a proper bearing stratum and the technical conditions of mining and utilization are determined, mineral product reserves and geological data required by mine construction design are provided, and investigation and research work on geological conditions such as rocks, stratums, structures, mineral products, hydrology, landforms and the like in a certain area is facilitated.
In the frequency domain electromagnetic method for exploring geological conditions, besides a plurality of interferences exist in a non-detection space, an electric non-uniform medium may exist in a near-surface or a detection blind area, and the topographic relief of the surface can also generate interference on electromagnetic waves, namely a static effect. Electromagnetic fields with different frequencies in a detection space, such as an electromagnetic field generated by industrial discrete current or a terrestrial electromagnetic field and a secondary electromagnetic field generated by the electromagnetic field and induced to a penetrated medium, can generate a superposition effect with the electromagnetic field for measurement, and the electromagnetic fields with different frequencies in a non-detection space can also generate a superposition effect with the electromagnetic field for measurement, so that the electromagnetic field for measurement can be inhibited or eliminated, the accuracy of a measurement result is influenced, and even the effective measurement cannot be realized. Static effects can affect the observed data value of the electric field, and the superposition effect of electromagnetic fields with different frequencies and electromagnetic fields for measurement mainly affects the observed data value of the magnetic field, such as the Carniian resistivity (Camanird resistivity) value representing the distribution of the magnetic and electric structures of the detected space medium, generates false information or masks useful information, affects the objectivity and authenticity of the electromagnetic detection result, and affects the cognition and reliability of the distribution of the magnetic and electric structures of the detected space medium.
Disclosure of Invention
In view of the above, it is necessary to provide a geological exploration detection method, apparatus, device and medium capable of filtering a noise signal interfering with a measurement electromagnetic wave so that a geological exploration detection result can reflect the true condition of a geological structure.
One aspect of the present application provides a method of geophysical prospecting detection comprising:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
and acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal.
In the geological exploration wave detection method in the above embodiment, the characteristic parameter values of the induction electromagnetic signals of each frequency are obtained by receiving the induction electromagnetic signals of geological exploration electromagnetic detection signals of different frequencies, where the characteristic parameter values include a wavelength value, a wave amplitude value, and a first time value; when the characteristic parameter value of the induction electromagnetic signal belongs to a preset threshold range, taking the induction electromagnetic signal with the frequency as an intrinsic induction electromagnetic signal; and acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal. The transmitted geological exploration electromagnetic detection signal is interfered by terrain, different transmission medium layers or other electromagnetic signals on the ground in the transmission process, so that a plurality of noise signals are mixed in the induction electromagnetic signal of the received geological exploration electromagnetic detection signal, the induction electromagnetic signal of the received geological exploration electromagnetic detection signal can not reflect a real measurement result, and the authenticity or the accuracy of the detection result of the geological exploration electromagnetic detection signal is influenced. The noise signals in the received geological exploration electromagnetic detection signals are filtered by comparing the characteristic parameter values of the received geological exploration electromagnetic detection signals with the predicted threshold range, so that the geological exploration electromagnetic detection signals capable of reflecting real geological structures are obtained, and the accuracy of geological exploration by using the electromagnetic signals is improved.
In one embodiment, the obtaining the characteristic parameter values of the induced electromagnetic signals at the frequencies includes:
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000021
Figure GDA0003259773360000022
Obtaining the first time value based on the waviness value;
wherein v isnValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThickness value of the nth medium.
In the geological exploration wave detection method in the embodiment, the weighted average value of the propagation speeds of the electromagnetic signals of all frequencies in various medium layers propagating through is calculated to be used as the waviness value of the electromagnetic signals of all frequencies
Figure GDA0003259773360000023
Acquiring the first time value based on the waviness value to acquire characteristic parameter values of the electromagnetic signals of various frequencies, and screening out geological exploration electromagnetic detection signals capable of reflecting real geological conditions based on the acquired characteristic parameter values to improve utilization of electromagnetic informationThe accuracy of the geological exploration is carried out.
In one embodiment, the obtaining the characteristic parameter value of the induced electromagnetic signal at each frequency further includes:
acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000025
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000026
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000024
in the geological exploration wave detection method in the embodiment, the first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point is obtained1And a second distance S from the test point to a ground electromagnetic signal receiving point2(ii) a Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000027
And calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000028
To calculate a first time value T of said electromagnetic signal at each frequency, said frequencyThe first time value T of the electromagnetic signal is the sum of the propagation time of the electromagnetic signal from the emitting point to the test point and twice the propagation time of the electromagnetic signal from the induction point to the ground test point. And acquiring characteristic parameter values of the electromagnetic signals of all frequencies based on the acquired first time value, and screening out geological exploration electromagnetic detection signals capable of reflecting real geological conditions based on the acquired characteristic parameter values so as to improve the accuracy of geological exploration by utilizing the electromagnetic signals.
In one embodiment, the obtaining the characteristic parameter values of the induced electromagnetic signals at the frequencies includes:
calculating the wave amplitude value of the electromagnetic signal of each frequency according to the following formula, wherein the wave amplitude value of the electromagnetic signal comprises the wave amplitude value E of electric field signals which are vertical to each otherRXAmplitude H of the sum magnetic field signalRX
Figure GDA0003259773360000031
Figure GDA0003259773360000032
Figure GDA0003259773360000033
Figure GDA0003259773360000034
Figure GDA0003259773360000035
Wherein E isTxAmplitude, H, of electric field signals in the initial electromagnetic signal emitted by the source of the electromagnetic waveTxThe amplitude value of a magnetic field signal in an initial electromagnetic signal emitted by the electromagnetic wave emission source, alpha is an attenuation coefficient, f is frequency, gamma is conductivity, and gamma isnConductivity for propagation of electromagnetic waves in the nth medium, GRxFor amplification in the reception of signals, GTxFor amplification when transmitting signals, lambdaεIs the wavelength of an electromagnetic wave in a medium having a dielectric constant of epsilonnIs the dielectric constant, μ, of an electromagnetic wave in the nth mediumnTo be the permeability of the electromagnetic wave in the nth medium,
Figure GDA0003259773360000036
is the average value of the electrical conductivity of the electromagnetic wave in the n media that propagate through,
Figure GDA0003259773360000037
is the average value of the wavelength values of the electromagnetic wave in the n media propagating through,
Figure GDA0003259773360000038
is the average value of the magnetic permeability of the electromagnetic wave in the n media which are propagated through, and r is the distance of the electromagnetic wave propagation.
In one embodiment, the acquiring a geophysical survey pickup signal based on the intrinsic induced electromagnetic signal comprises:
fitting the intrinsic induction electromagnetic signal with a preset theoretical electromagnetic wave waveform curve;
and acquiring the intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal.
In the geological exploration wave detection method in the above embodiment, because the induced electromagnetic signal of the received geological exploration electromagnetic detection signal is mixed with the noise signal, the propagation time and the wavelength value of the induced electromagnetic signal mixed with the noise signal deviate from the geological exploration electromagnetic detection signal under the condition of the theoretical non-interference signal, and the intrinsic induced electromagnetic signal capable of reflecting the real geological condition is screened out, so that the intrinsic induced electromagnetic signal is fitted with the preset theoretical electromagnetic wave waveform curve; and acquiring the intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal, thereby acquiring a geological exploration detection result signal capable of reflecting real geological conditions to acquire a geological structure or a preset geological measurement result.
In one embodiment, the fitting the intrinsic induction electromagnetic signal to a preset theoretical electromagnetic wave waveform curve includes:
generating a preset theoretical electromagnetic wave waveform curve, wherein the theoretical electromagnetic wave waveform curve comprises a theoretical magnetic field curve and a theoretical electric field curve which change in a sine rule;
comparing the wavelength value and the amplitude value of the intrinsic induction electromagnetic signal with the wavelength value and the amplitude value of the theoretical electromagnetic wave waveform curve, respectively.
In one embodiment, the geological exploration detection method further comprises:
and acquiring a second time value between transmission of the geological exploration electromagnetic detection electromagnetic wave from the transmitting point to the test point through the inductive switch.
In the geological exploration wave detection method in the above embodiment, an inductive switch may be disposed at the detection point, and when the geological exploration electromagnetic detection electromagnetic wave propagates to the detection point, the inductive switch may detect a second time value between transmission of the geological exploration electromagnetic detection electromagnetic wave from the transmission point to the detection point, so as to reduce the calculation complexity of the first time value in the characteristic parameter values of the electromagnetic signal of each frequency.
In one embodiment, the geological exploration detection method further comprises:
acquiring the frequency f of the geological exploration result signal;
obtaining the electric field intensity E of the geological exploration result signalx
Obtaining the magnetic field intensity H of the geological exploration result signaly
Calculating the Carniya resistivity rho according to the following formulas
Figure GDA0003259773360000041
Geological exploration wave detection method in the above embodimentBased on the obtained frequency f and electric field intensity E of the geological exploration result signalxMagnetic field intensity HyAnd calculating the Carnia resistivity to obtain the geological exploration detection measurement result.
An aspect of the application provides a geological exploration detection device, includes: the device comprises an electromagnetic signal acquisition device, a signal processing device and a signal processing device, wherein the electromagnetic signal acquisition device is used for acquiring electromagnetic signals with different frequencies, and the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals; the tuning subsystem is used for acquiring characteristic parameter values of the electromagnetic signals of all frequencies and acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values; and the fitting subsystem is used for acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal.
In the geological exploration detection device in the above embodiment, the electromagnetic signal acquisition device is used to receive the induced electromagnetic signals of the geological exploration electromagnetic detection signals with different frequencies, and obtain characteristic parameter values of the induced electromagnetic signals with different frequencies, where the characteristic parameter values include a wavelength value, a wave amplitude value, and a first time value; when the characteristic parameter value of the induction electromagnetic signal belongs to a preset threshold range, the tuning subsystem takes the induction electromagnetic signal with the frequency as an intrinsic induction electromagnetic signal; and the fitting subsystem acquires a geological exploration detection result signal based on the intrinsic induction electromagnetic signal. The transmitted geological exploration electromagnetic detection signal is interfered by terrain, different transmission medium layers or other electromagnetic signals on the ground in the transmission process, so that a plurality of noise signals are mixed in the induction electromagnetic signal of the received geological exploration electromagnetic detection signal, the induction electromagnetic signal of the received geological exploration electromagnetic detection signal can not reflect a real measurement result, and the authenticity or the accuracy of the detection result of the geological exploration electromagnetic detection signal is influenced. The noise signals in the received geological exploration electromagnetic detection signals are filtered by comparing the characteristic parameter values of the received geological exploration electromagnetic detection signals with the predicted threshold range, so that the geological exploration electromagnetic detection signals capable of reflecting real geological structures are obtained, and the accuracy of geological exploration by using the electromagnetic signals is improved.
In one embodiment, the geological exploration detection device further comprises:
and the inductive switch is used for acquiring a second time value between the transmission of the geological exploration electromagnetic detection electromagnetic wave from the transmitting point to the test point.
In the geological exploration detecting device in the above embodiment, an inductive switch may be disposed at the detection point, and when the geological exploration electromagnetic detection wave propagates to the detection point, the inductive switch may detect a second time value between transmission of the geological exploration electromagnetic detection wave from the transmission point to the detection point, so as to reduce the calculation complexity of the first time value in the characteristic parameter values of the electromagnetic signal of each frequency.
In one embodiment, a Carniian resistivity acquisition module is arranged in the fitting subsystem and used for acquiring the frequency f and the electric field intensity E of the geological exploration result signalxAnd magnetic field strength Hy
Calculating the Carniya resistivity rho according to the following formulas
Figure GDA0003259773360000042
In the geological exploration detecting device in the above embodiment, the frequency f and the electric field strength E of the geological exploration result signal are obtainedxMagnetic field intensity HyAnd calculating the Carnia resistivity to obtain the geological exploration detection measurement result.
An aspect of the present application provides a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the method described in any of the embodiments of the present application when executing the computer program.
Another aspect of the present application provides a computer-readable storage medium, on which a computer program is stored, wherein the computer program, when executed by a processor, implements the steps of any of the methods described in the embodiments of the present application.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art to obtain drawings of other embodiments based on these drawings without any creative effort.
Fig. 1 is a diagram illustrating an application scenario of a method for detecting geological exploration according to a first embodiment of the present application.
FIG. 2 is a schematic flow chart of a method of seismic survey detection provided in a second embodiment of the present application.
FIG. 3 is a schematic flow chart of a method of seismic survey detection provided in a third embodiment of the present application.
FIG. 4 is a block diagram of a seismic survey detection device provided in a fourth embodiment of the present application.
FIG. 5 is a block diagram of a seismic survey detection device provided in a fifth embodiment of the present application.
FIG. 6 is a block diagram of a seismic survey detecting device according to a sixth embodiment of the present application.
Fig. 7 is an internal structural diagram of a computer device provided in an embodiment of the present application.
Detailed Description
To facilitate an understanding of the present application, the present application will now be described more fully with reference to the accompanying drawings. Preferred embodiments of the present application are illustrated in the accompanying drawings. This application may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used herein in the description of the present application is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
Where the terms "comprising," "having," and "including" are used herein, another element may be added unless an explicit limitation is used, such as "only," "consisting of … …," etc. Unless mentioned to the contrary, terms in the singular may include the plural and are not to be construed as being one in number.
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
The frequency domain electromagnetic method exploration is characterized in that except the existence of interference in a detection space, the interference exists in the electrical non-uniformity or the topographic relief of a shallow medium close to the earth surface or a detection blind area, namely, the static effect; electromagnetic fields with different frequencies in a detection space, such as an electromagnetic field generated by industrial discrete current and a terrestrial electromagnetic field and a secondary electromagnetic field generated by the electromagnetic fields and the terrestrial electromagnetic field induced to a penetrated medium, generate a superposition effect with a useful electromagnetic field and generate a superposition effect with the electromagnetic field and the useful electromagnetic field with different frequencies in a non-detection space, and the superposition effect is suppressed or eliminated by the existing device.
Generally speaking, the static effect only affects the electric field observation data value, and the superposition effect of the electromagnetic fields with different frequencies and the useful electromagnetic field mainly affects the magnetic field observation data value, such as the objective and real kania resistivity values characterizing the magnetic and dielectric medium structure distribution of the detection space are affected, false information is generated or the useful information is covered, wrong recognition is formed, and the cognition degree and the reliability of the magnetic and dielectric medium structure distribution of the detection space are affected.
In one embodiment of the present application, a method for detecting geological exploration is provided, which can be applied in the application environment as shown in fig. 1. Wherein the first terminal 102 communicates with the server 104 via a network. Specifically, the first terminal 102 may display an input interface of the geological exploration detection method, and the input interface may be a selective input mode by inputting a command, for example, a transmission command of geological exploration electromagnetic detection signals with different frequencies is selectively input. The first terminal 102 may be, but not limited to, various personal computers, notebook computers, smart phones, tablet computers, and portable wearable devices, and the server 104 may be implemented by an independent server or a server cluster formed by a plurality of servers. It should be noted that, in the present embodiment, the first terminal 102 communicates with the server 104 through a network, and data required for geological exploration measurement through electromagnetic signals or basic geographic data in a measurement area may be acquired through the server, including but not limited to speed values of electromagnetic waves propagating in different medium layers, known terrain parameters in a preset area, and the like. For example, the real-time information data running on each extension server and equipment connected to the switchboard server through the network is acquired. Then, the first terminal 102 acquires characteristic parameter values of the electromagnetic signals of each frequency by acquiring geological exploration electromagnetic detection signals or induction electromagnetic signals of the geological exploration electromagnetic detection signals in a preset frequency band, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values; acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range; and acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal, and further acquiring a geological exploration measurement result presented in a preset output form.
Further, in an embodiment of the present application, as shown in fig. 2, there is provided a method for detecting geological exploration, which is exemplified by the application of the method to the first terminal in fig. 1, and includes the following steps:
step 202, obtaining electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals.
The first terminal is a terminal with an operable interface, the first terminal interface displays an input window and can also display a webpage, and a user can input data through the input window and can also browse the webpage to acquire data stored locally or on a server connected with the first terminal through a network. The web application may be a browser or other application that can display page content (e.g., industrial control software applications, instant messaging applications, etc.). The frequency of the electromagnetic signal transmitted for geological exploration electromagnetic detection can be displayed in the webpage, and after a user selects a preset measuring area, the first terminal can acquire basic geographic data, common geological structure data and the like in the measuring area through the server. After the first terminal acquires the induction electromagnetic signals of the geological exploration electromagnetic detection signals with different frequencies, the induction electromagnetic signals can be processed or analyzed based on acquired basic geographic data and common geological structure data in the preset area.
And 204, acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values.
Due to the fact that characteristic parameter values such as a wavelength value, a wave amplitude value and a first time value of electromagnetic signals with different frequencies acquired by the first terminal are different, the geological exploration electromagnetic detection signals are influenced by landforms, geological stratum structures, dielectrics, noise signals and the like in the transmission process, the characteristic parameter values such as the wavelength value, the wave amplitude value and the first time value can change, and certain errors exist between the characteristic parameter values corresponding to the geological exploration electromagnetic detection signals or induction electromagnetic signals of the geological exploration electromagnetic detection signals acquired under the condition that no theoretical interference or noise exists.
More specifically, the acquiring characteristic parameter values of the induced electromagnetic signals at various frequencies includes:
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000061
Figure GDA0003259773360000062
Obtaining the first time value based on the waviness value;
wherein v isnValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThickness value of the nth medium.
In the geological exploration wave detection method in the embodiment, the weighted average value of the propagation speeds of the electromagnetic signals of all frequencies in various medium layers propagating through is calculated to be used as the waviness value of the electromagnetic signals of all frequencies
Figure GDA0003259773360000071
And acquiring the first time value based on the waviness value to acquire characteristic parameter values of the electromagnetic signals of various frequencies, and screening out geological exploration electromagnetic detection signals capable of reflecting real geological conditions based on the acquired characteristic parameter values to improve the accuracy of geological exploration by utilizing the electromagnetic signals.
More specifically, the acquiring characteristic parameter values of the induced electromagnetic signals at various frequencies includes:
acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000074
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000075
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000072
in the geological exploration wave detection method in the embodiment, the first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point is obtained1And a second distance S from the test point to a ground electromagnetic signal receiving point2(ii) a Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000076
And calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000077
And calculating a first time value T of the electromagnetic signal of each frequency, wherein the first time value T of the electromagnetic signal of each frequency is the sum of the propagation time of the electromagnetic signal from the emitting point to the test point and twice the propagation time of the electromagnetic signal from the induction point to the ground test point. And acquiring characteristic parameter values of the electromagnetic signals of all frequencies based on the acquired first time value, and screening out geological exploration electromagnetic detection signals capable of reflecting real geological conditions based on the acquired characteristic parameter values so as to improve the accuracy of geological exploration by utilizing the electromagnetic signals.
More specifically, in the geological exploration wave detection method in the above embodiment, the amplitude is divided into electric field amplitude and magnetic field amplitude. The amplitude of the electromagnetic wave of the same frequency is related to the energy, i.e., power, of the electromagnetic wave. The electromagnetic wave energy or power, the electric field amplitude and the magnetic field amplitude at the emission source can be calculated by corresponding formulas and can also be measured by instruments. The amplitudes of the electric field and the magnetic field at the receiving position or the information acquisition position are obtained through calculation, and the reasons are that the electromagnetic wave energy is attenuated due to factors such as distance, media and the like in the transmission process, and the static effect and the superposition effect are caused. In resource exploration, a medium in an electromagnetic wave propagation process is very complex, absorption and scattering of medium attenuation factors or energy attenuation calculation are simplified, and only distance and absorption factors are considered. Then the electromagnetic wave energy attenuation formula:
Figure GDA0003259773360000073
in the above formula, PRxElectromagnetic wave energy or power at a receiving place or an information acquisition place/point is W/mW; pTxIs the energy value or power value of the emitted electromagnetic wave, and the unit is W/mW; gTxAmplification or gain factor (dB) when transmitting electromagnetic signals; gTxAmplification or gain factor (dB) in receiving electromagnetic signals; alpha is an attenuation coefficient and has the unit of NP/m; λ is the electromagnetic wave wavelength, in m; r is the propagation distance of the electromagnetic wave in m. It should be noted that: pTxThe energy conservation law of the secondary field generated by sensing the depth medium by transmitting electromagnetic waves to the detection depth, and the energy of the electromagnetic waves transmitted to the detection depth is converted into the energy of the secondary field and is unchanged, so that P in the formulaTxThe remaining energy from the emission energy of the electromagnetic wave can be considered, and the propagation distance r is the sum of the distance from the emission position to the observation point, the distance from the observation point to the secondary field source and the distance from the secondary field source to the observation point.
More specifically, the acquiring characteristic parameter values of the induced electromagnetic signals at various frequencies includes:
calculating the wave amplitude value of the electromagnetic signal of each frequency according to the following formula, wherein the wave amplitude value of the electromagnetic signal comprises the wave amplitude value E of electric field signals which are vertical to each otherRXAmplitude H of the sum magnetic field signalRX
Figure GDA0003259773360000081
Figure GDA0003259773360000082
Figure GDA0003259773360000083
Figure GDA0003259773360000084
Figure GDA0003259773360000085
Figure GDA0003259773360000086
Wherein E isTxAmplitude, H, of electric field signals in the initial electromagnetic signal emitted by the source of the electromagnetic waveTxThe amplitude value of a magnetic field signal in an initial electromagnetic signal emitted by the electromagnetic wave emission source, alpha is an attenuation coefficient, f is frequency, gamma is conductivity, and gamma isnConductivity for propagation of electromagnetic waves in the nth medium, GRxFor amplification in the reception of signals, GTxFor amplification when transmitting signals, lambdaεIs the wavelength of an electromagnetic wave in a medium having a dielectric constant of epsilonnIs the dielectric constant, μ, of an electromagnetic wave in the nth mediumnTo be the permeability of the electromagnetic wave in the nth medium,
Figure GDA0003259773360000087
is the average value of the electrical conductivity of the electromagnetic wave in the n media that propagate through,
Figure GDA0003259773360000088
is the average value of the wavelength values of the electromagnetic wave in the n media propagating through,
Figure GDA0003259773360000089
the method is characterized in that the method is an average value of magnetic permeability of electromagnetic waves in n media which are propagated through, r is the distance of electromagnetic wave propagation, and the propagation distance r is the sum of the distance value from an electromagnetic signal emission position to an observation point, the distance value from the observation point to a secondary field source and the distance from the secondary field source to an electromagnetic wave signal receiving point.
And step 206, acquiring the intrinsic induction electromagnetic signal of which the characteristic parameter value belongs to a preset threshold range.
By obtaining characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values, geological exploration electromagnetic detection signals or induction electromagnetic signals of the geological exploration electromagnetic detection signals, of which the characteristic parameter values belong to a preset threshold range, are screened out.
And step 208, acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal.
The acquired geological exploration electromagnetic detection signals or induction electromagnetic signals thereof with various frequencies are possibly influenced by terrain, geological features, geological stratum structures, dielectrics, noise signals and the like in the propagation process, characteristic parameter values such as wavelength values, wave amplitude values and first time values of the acquired geological exploration electromagnetic detection signals or induction electromagnetic signals thereof change, certain errors exist between the characteristic parameter values corresponding to the geological exploration electromagnetic detection signals or the induction electromagnetic signals thereof acquired under the condition that the theory is free of interference or noise, and the fitting degree between the electromagnetic signals with various noise interferences filtered out and the characteristic parameter values corresponding to the geological exploration electromagnetic detection signals or the induction electromagnetic signals thereof acquired under the condition that the theory is free of interference or noise is higher. Therefore, by acquiring the intrinsic induction electromagnetic signal of which the characteristic parameter value belongs to the preset threshold range, and acquiring the geological exploration detection result signal capable of reflecting the real geological structure based on the intrinsic induction electromagnetic signal, the geological exploration measurement result in the preset form can be acquired. The geological exploration wave detection method provided by the embodiment of the application can inhibit or eliminate static effect and superposition effect, so that the acquired data information is close to objective. Even if the interference of a special case in the superposition effect is inhibited or eliminated possibly difficultly, so that the measurement effect is possibly poor, various medium samples in a measurement area can be scientifically collected and the physical properties of the samples can be measured, medium parameters are screened by utilizing big data, and the maximum and minimum electric field and magnetic field amplitude values approaching the objective of the measurement area are calculated by utilizing the corresponding formulas provided above, so that the amplitude of the electric field and magnetic field sine curves is tuned to fit and constrain the collected electric field and magnetic field amplitude values, and the interference phenomenon is effectively inhibited or eliminated.
Further, in an embodiment of the present application, a method for detecting geological exploration is provided, further comprising:
and acquiring a second time value between transmission of the geological exploration electromagnetic detection electromagnetic wave from the transmitting point to the test point through the inductive switch.
Specifically, an inductive switch may be provided at the probing point, and the inductive switch may detect a second time value between transmission of the geophysical prospecting electromagnetic detection waves from the transmission point to the testing point when the geophysical prospecting electromagnetic detection waves propagate to the probing point.
In the geological exploration wave detection method in the above embodiment, the inductive switch is arranged at the detection point, and when the geological exploration electromagnetic detection electromagnetic wave propagates to the detection point, the inductive switch can detect a second time value between transmission of the geological exploration electromagnetic detection electromagnetic wave from the transmission point to the detection point, so that the calculation complexity of the first time value in the characteristic parameter values of the electromagnetic signals of each frequency can be reduced.
Further, in an embodiment of the present application, as shown in fig. 3, there is provided a geological exploration detection method, which is described by taking the method as an example of being applied to the first terminal in fig. 1, and the geological exploration detection method further includes:
step 2081, fitting the intrinsic induction electromagnetic signal with a preset theoretical electromagnetic wave waveform curve.
Specifically, discrete points on the intrinsic induced electromagnetic signal with the characteristic parameter value falling within a preset threshold range may be obtained, and the discrete points may be fitted to the geological exploration detection result signal under the theoretical noise-free and interference condition, for example, the fitting may be performed by a minimum two-way method, so as to obtain a curve with the highest fitting degree to the geological exploration detection result signal curve under the theoretical noise-free and interference condition as the geological exploration detection result signal curve.
And 2082, obtaining the intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal.
By acquiring the curve with the highest fitting degree with the geological exploration detection result signal curve under the condition of no noise and interference with the theory as the geological exploration detection result signal curve, the acquired geological exploration detection result signal is closer to the geological exploration detection result signal acquired under the ideal condition, so that the geological exploration measurement result can be more real according to the acquired geological exploration detection result signal with the highest fitting degree, and the accuracy of the geological exploration electromagnetic detection result is improved.
Further, in a geological survey detection method provided in an embodiment of the present application, the method further includes:
acquiring the frequency f and the electric field intensity E of the geological exploration result signalxAnd magnetic field strength HyCalculating the Carniya resistivity rho according to the following formulas
Figure GDA0003259773360000091
In the geological exploration wave detection method in the above embodiment, the frequency f and the electric field strength E of the geological exploration result signal are obtainedxMagnetic field intensity HyAnd calculating the Carnia resistivity to obtain the geological exploration detection measurement result. Since the geophysical survey result signal has been compared to a geophysical survey pickup result signal that can be obtained under ideal conditions, the Carniya resistivity ρ calculated based on the geological survey result signalsAnd also more reflective of the true geological structure.
It should be understood that although the various steps in the flow charts of fig. 2-3 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least some of the steps in fig. 2-3 may include multiple sub-steps or multiple stages that are not necessarily performed at the same time, but may be performed at different times, and the order of performance of the sub-steps or stages is not necessarily sequential, but may be performed in turn or alternating with other steps or at least some of the sub-steps or stages of other steps.
In one embodiment of the present application, as shown in FIG. 4, there is provided a geophysical prospecting detection apparatus comprising: an electromagnetic signal acquisition device 20, a tuning subsystem 40, and a fitting subsystem 60, wherein:
the electromagnetic signal acquisition device 20 is used for acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
the tuning subsystem 40 is configured to obtain characteristic parameter values of the electromagnetic signals of each frequency, and obtain intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range, where the characteristic parameter values include a wavelength value, a wave amplitude value, and a first time value;
and a fitting subsystem 60 for obtaining a geophysical exploration demodulation result signal based on the intrinsic induced electromagnetic signal.
In particular, the tuning subsystem 40 is designed with a "wavelength, time" scale that cannot be calculated, i.e., tuned in wavelength, for natural sources. The fitting subsystem 60 resembles a "filter" that removes or suppresses the filtering superimposed on the "sinusoidally varying magnetic or electric field". A set of sub-system capable of generating a 'sine-law changing magnetic field or electric field theoretical curve' can be designed in the fitting sub-system 60, fitting degree is carried out on the 'sine curve by adjusting the wavelength and the amplitude' of the sine curve and an actually received magnetic field or electric field or discrete amplitude particle curve, the fitting degree is the best and is used as a magnetic field and electric field observation data value of the frequency point of the measuring point, and therefore the Carniya resistivity of the frequency point of the measuring point is calculated. The wavelengths required by the fitting subsystem 60 and the tuning subsystem are both taken as value sensing positions, an induced electromagnetic field induced by a transmitting electromagnetic field of a medium at a certain depth below the ground is a secondary field, and the secondary field is transmitted to a wavelength value calculated by a stratum medium which passes through a ground observation point and finally reaches the ground observation point. By setting the static effect correction originally performed in the data processing in the observation instrument, the workload of data processing and the subjectivity of the static effect correction method selection can be reduced. The superposition effect of industrial discrete current and earth electromagnetic field in a detection airspace and secondary electromagnetic field induced by the industrial discrete current and the earth electromagnetic field is inhibited or eliminated, and the signal-to-noise ratio is improved; data observation is possible to be carried out in the transition region, and the working efficiency is improved.
In the geological exploration detection device in the above embodiment, the electromagnetic signal acquisition device is used to receive the induced electromagnetic signals of the geological exploration electromagnetic detection signals with different frequencies, and obtain characteristic parameter values of the induced electromagnetic signals with different frequencies, where the characteristic parameter values include a wavelength value, a wave amplitude value, and a first time value; when the characteristic parameter value of the induction electromagnetic signal belongs to a preset threshold range, the tuning subsystem takes the induction electromagnetic signal with the frequency as an intrinsic induction electromagnetic signal; and the fitting subsystem acquires a geological exploration detection result signal based on the intrinsic induction electromagnetic signal. The transmitted geological exploration electromagnetic detection signal is interfered by terrain, different transmission medium layers or other electromagnetic signals on the ground in the transmission process, so that a plurality of noise signals are mixed in the induction electromagnetic signal of the received geological exploration electromagnetic detection signal, the induction electromagnetic signal of the received geological exploration electromagnetic detection signal can not reflect a real measurement result, and the authenticity or the accuracy of the detection result of the geological exploration electromagnetic detection signal is influenced. The noise signals in the received geological exploration electromagnetic detection signals are filtered by comparing the characteristic parameter values of the received geological exploration electromagnetic detection signals with the predicted threshold range, so that the geological exploration electromagnetic detection signals capable of reflecting real geological structures are obtained, and the accuracy of geological exploration by using the electromagnetic signals is improved.
Further, in a geological exploration detecting device provided in an embodiment of the present application, as shown in fig. 5, an inductive switch 50 is further included for obtaining a second time value between transmission of the geological exploration electromagnetic detection electromagnetic wave from the emission point to the test point.
In the geological exploration detecting device in the above embodiment, an inductive switch may be disposed at the detection point, and when the geological exploration electromagnetic detection wave propagates to the detection point, the inductive switch may detect a second time value between transmission of the geological exploration electromagnetic detection wave from the transmission point to the detection point, so as to reduce the calculation complexity of the first time value in the characteristic parameter values of the electromagnetic signal of each frequency.
Further, in a geological exploration detecting device provided in an embodiment of the present application, as shown in fig. 6, a carney resistivity obtaining module 61 is disposed in the fitting subsystem 60 for obtaining the frequency f and the electric field strength E of the geological exploration result signalxAnd magnetic field strength Hy
Calculating the Carniya resistivity rho according to the following formulas
Figure GDA0003259773360000111
In the geological exploration detecting device in the above embodiment, the frequency f and the electric field strength E of the geological exploration result signal are obtainedxMagnetic field intensity HyAnd calculating the Carnia resistivity to obtain the geological exploration detection measurement result.
For specific limitations of the geophysical prospecting detection means, reference is made to the above limitations of the geological prospecting detection method, which are not described in detail here.
In one embodiment of the present application, a computer device is provided, and the computer device may be a terminal, and the internal structure diagram thereof may be as shown in fig. 7. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to implement a method of prospecting detection. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on the shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.
Those skilled in the art will appreciate that the architecture shown in fig. 7 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment of the present application, there is provided a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the following steps when executing the computer program:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
and acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal.
In one embodiment of the application, the processor when executing the computer program further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000112
Figure GDA0003259773360000113
Obtaining the first time value based on the waviness value;
acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal;
wherein v isnValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThickness value of the nth medium.
In one embodiment of the application, the processor when executing the computer program further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000121
Figure GDA0003259773360000122
Acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000129
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA00032597733600001210
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000123
acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal;
wherein v isnValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThickness value of the nth medium.
In one embodiment of the application, the processor when executing the computer program further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000124
Figure GDA0003259773360000125
Acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining the test point to ground potentialSecond distance S of magnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000126
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000127
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000128
acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
fitting the intrinsic induction electromagnetic signal with a preset theoretical electromagnetic wave waveform curve;
obtaining an intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal;
wherein v isnValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThickness value of the nth medium.
In one embodiment of the application, the processor when executing the computer program further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000131
Figure GDA0003259773360000132
Acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000133
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000134
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000135
calculating the wave amplitude value of the electromagnetic signal of each frequency according to the following formula, wherein the wave amplitude value of the electromagnetic signal comprises the wave amplitude value E of electric field signals which are vertical to each otherRXAmplitude H of the sum magnetic field signalRX
Figure GDA0003259773360000136
Figure GDA0003259773360000137
Figure GDA0003259773360000138
Figure GDA0003259773360000139
Figure GDA00032597733600001310
Figure GDA00032597733600001311
Acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
fitting the intrinsic induction electromagnetic signal with a preset theoretical electromagnetic wave waveform curve;
obtaining an intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal;
wherein E isTxAmplitude, H, of electric field signals in the initial electromagnetic signal emitted by the source of the electromagnetic waveTxAmplitude, v, of a magnetic field signal in an initial electromagnetic signal emitted by said electromagnetic wave emission sourcenValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThe thickness value of the nth medium, alpha is attenuation coefficient, f is frequency, gamma is conductivity, gammanConductivity for propagation of electromagnetic waves in the nth medium, GRxFor amplification in the reception of signals, GTxFor amplification when transmitting signals, lambdaεIs the wavelength of an electromagnetic wave in a medium having a dielectric constant of epsilonnIs the dielectric constant, μ, of an electromagnetic wave in the nth mediumnTo be the permeability of the electromagnetic wave in the nth medium,
Figure GDA0003259773360000141
is the average value of the electrical conductivity of the electromagnetic wave in the n media that propagate through,
Figure GDA0003259773360000142
is the average value of the wavelength values of the electromagnetic wave in the n media propagating through,
Figure GDA0003259773360000143
is the average value of the magnetic permeability of the electromagnetic wave in the n media which are propagated through, and r is the distance of the electromagnetic wave propagation.
In one embodiment of the application, the processor when executing the computer program further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000144
Figure GDA0003259773360000145
Acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000146
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000147
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000148
calculating the wave amplitude value of the electromagnetic signal of each frequency according to the following formula, wherein the wave amplitude value of the electromagnetic signal comprises the wave amplitude value E of electric field signals which are vertical to each otherRXAmplitude H of the sum magnetic field signalRX
Figure GDA0003259773360000149
Figure GDA00032597733600001410
Figure GDA0003259773360000151
Figure GDA0003259773360000152
Figure GDA0003259773360000153
Figure GDA0003259773360000154
Acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
generating a preset theoretical electromagnetic wave waveform curve, wherein the theoretical electromagnetic wave waveform curve comprises a theoretical magnetic field curve and a theoretical electric field curve which change in a sine rule;
comparing the wavelength value and the amplitude value of the intrinsic induction electromagnetic signal with the wavelength value and the amplitude value of the theoretical electromagnetic wave waveform curve, respectively;
obtaining an intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal;
wherein E isTxAmplitude, H, of electric field signals in the initial electromagnetic signal emitted by the source of the electromagnetic waveTxAmplitude, v, of a magnetic field signal in an initial electromagnetic signal emitted by said electromagnetic wave emission sourcenValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThe thickness value of the nth medium, alpha is attenuation coefficient, f is frequency, gamma is conductivity, gammanConductivity for propagation of electromagnetic waves in the nth medium, GRxFor amplification in the reception of signals, GTxFor amplification when transmitting signals, lambdaεIs the wavelength of an electromagnetic wave in a medium having a dielectric constant of epsilonnIs the dielectric constant, μ, of an electromagnetic wave in the nth mediumnTo be the permeability of the electromagnetic wave in the nth medium,
Figure GDA0003259773360000155
is the average value of the electrical conductivity of the electromagnetic wave in the n media that propagate through,
Figure GDA0003259773360000156
is the average value of the wavelength values of the electromagnetic wave in the n media propagating through,
Figure GDA0003259773360000157
is the average value of the magnetic permeability of the electromagnetic wave in the n media which are propagated through, and r is the distance of the electromagnetic wave propagation.
In one embodiment of the application, the processor when executing the computer program further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
according to the followingCalculating the waviness value of the electromagnetic signal of each frequency by a formula
Figure GDA0003259773360000158
Figure GDA0003259773360000159
Acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA00032597733600001510
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA00032597733600001511
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000161
calculating the wave amplitude value of the electromagnetic signal of each frequency according to the following formula, wherein the wave amplitude value of the electromagnetic signal comprises the wave amplitude value E of electric field signals which are vertical to each otherRXAmplitude H of the sum magnetic field signalRX
Figure GDA0003259773360000162
Figure GDA0003259773360000163
Figure GDA0003259773360000164
Figure GDA0003259773360000165
Figure GDA0003259773360000166
Figure GDA0003259773360000167
Acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
generating a preset theoretical electromagnetic wave waveform curve, wherein the theoretical electromagnetic wave waveform curve comprises a theoretical magnetic field curve and a theoretical electric field curve which change in a sine rule;
comparing the wavelength value and the amplitude value of the intrinsic induction electromagnetic signal with the wavelength value and the amplitude value of the theoretical electromagnetic wave waveform curve, respectively;
obtaining an intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal;
acquiring the frequency f of the geological exploration result signal;
obtaining the electric field intensity E of the geological exploration result signalx
Obtaining the magnetic field intensity H of the geological exploration result signaly
Calculating the Carniya resistivity rho according to the following formulas
Figure GDA0003259773360000168
Wherein E isTxFor electromagnetic wave emission source emissionOf the initial electromagnetic signal, HTxAmplitude, v, of a magnetic field signal in an initial electromagnetic signal emitted by said electromagnetic wave emission sourcenValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThe thickness value of the nth medium, alpha is attenuation coefficient, f is frequency, gamma is conductivity, gammanConductivity for propagation of electromagnetic waves in the nth medium, GRxFor amplification in the reception of signals, GTxFor amplification when transmitting signals, lambdaεIs the wavelength of an electromagnetic wave in a medium having a dielectric constant of epsilonnIs the dielectric constant, μ, of an electromagnetic wave in the nth mediumnTo be the permeability of the electromagnetic wave in the nth medium,
Figure GDA0003259773360000171
is the average value of the electrical conductivity of the electromagnetic wave in the n media that propagate through,
Figure GDA0003259773360000172
is the average value of the wavelength values of the electromagnetic wave in the n media propagating through,
Figure GDA0003259773360000173
is the average value of the magnetic permeability of the electromagnetic wave in the n media which are propagated through, and r is the distance of the electromagnetic wave propagation.
In an embodiment of the application, a computer-readable storage medium is provided, on which a computer program is stored, which computer program, when being executed by a processor, carries out the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
and acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal.
In one embodiment of the application, the computer program when executed by the processor further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000174
Figure GDA0003259773360000175
Obtaining the first time value based on the waviness value;
acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal;
wherein v isnValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThickness value of the nth medium.
In one embodiment of the application, the computer program when executed by the processor further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000176
Figure GDA0003259773360000177
Acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000178
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000179
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000181
acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal;
wherein v isnValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThickness value of the nth medium.
In one embodiment of the application, the computer program when executed by the processor further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000182
Figure GDA0003259773360000183
Acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000184
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000185
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000186
acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
fitting the intrinsic induction electromagnetic signal with a preset theoretical electromagnetic wave waveform curve;
obtaining an intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal;
wherein v isnValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThickness value of the nth medium.
In one embodiment of the application, the computer program when executed by the processor further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000187
Figure GDA0003259773360000188
Acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000191
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000192
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000193
calculating the wave amplitude value of the electromagnetic signal of each frequency according to the following formula, wherein the wave amplitude value of the electromagnetic signal comprises electric field signals which are vertical to each otherAmplitude of wave ERXAmplitude H of the sum magnetic field signalRX
Figure GDA0003259773360000194
Figure GDA0003259773360000195
Figure GDA0003259773360000196
Figure GDA0003259773360000197
Figure GDA0003259773360000198
Figure GDA0003259773360000199
Acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
fitting the intrinsic induction electromagnetic signal with a preset theoretical electromagnetic wave waveform curve;
obtaining an intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal;
wherein E isTxAmplitude, H, of electric field signals in the initial electromagnetic signal emitted by the source of the electromagnetic waveTxAmplitude, v, of a magnetic field signal in an initial electromagnetic signal emitted by said electromagnetic wave emission sourcenValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThe thickness value of the nth medium, alpha is attenuation coefficient, f is frequency, gamma is conductivity, gammanConductivity for propagation of electromagnetic waves in the nth medium, GRxFor amplification in the reception of signals, GTxFor amplification when transmitting signals, lambdaεIs the wavelength of an electromagnetic wave in a medium having a dielectric constant of epsilonnIs the dielectric constant, μ, of an electromagnetic wave in the nth mediumnTo be the permeability of the electromagnetic wave in the nth medium,
Figure GDA00032597733600001910
is the average value of the electrical conductivity of the electromagnetic wave in the n media that propagate through,
Figure GDA00032597733600001911
is the average value of the wavelength values of the electromagnetic wave in the n media propagating through,
Figure GDA00032597733600001912
is the average value of the magnetic permeability of the electromagnetic wave in the n media which are propagated through, and r is the distance of the electromagnetic wave propagation.
In one embodiment of the application, the computer program when executed by the processor further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000201
Figure GDA0003259773360000202
Acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000207
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000208
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
calculating the wave amplitude value of the electromagnetic signal of each frequency according to the following formula, wherein the wave amplitude value of the electromagnetic signal comprises the wave amplitude value E of electric field signals which are vertical to each otherRXAmplitude H of the sum magnetic field signalRX
Figure GDA0003259773360000203
Figure GDA0003259773360000204
Figure GDA0003259773360000205
Figure GDA0003259773360000206
Acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
generating a preset theoretical electromagnetic wave waveform curve, wherein the theoretical electromagnetic wave waveform curve comprises a theoretical magnetic field curve and a theoretical electric field curve which change in a sine rule;
comparing the wavelength value and the amplitude value of the intrinsic induction electromagnetic signal with the wavelength value and the amplitude value of the theoretical electromagnetic wave waveform curve, respectively;
obtaining an intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal;
wherein E isTxAmplitude, H, of electric field signals in the initial electromagnetic signal emitted by the source of the electromagnetic waveTxAmplitude, v, of a magnetic field signal in an initial electromagnetic signal emitted by said electromagnetic wave emission sourcenValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThe thickness value of the nth medium, alpha is attenuation coefficient, f is frequency, gamma is conductivity, gammanConductivity for propagation of electromagnetic waves in the nth medium, GRxFor amplification in the reception of signals, GTxFor amplification when transmitting signals, lambdaεIs the wavelength of an electromagnetic wave in a medium having a dielectric constant of epsilonnIs the dielectric constant, μ, of an electromagnetic wave in the nth mediumnTo be the permeability of the electromagnetic wave in the nth medium,
Figure GDA0003259773360000211
is the average value of the electrical conductivity of the electromagnetic wave in the n media that propagate through,
Figure GDA0003259773360000212
is the average value of the wavelength values of the electromagnetic wave in the n media propagating through,
Figure GDA0003259773360000213
is the average value of the magnetic permeability of the electromagnetic wave in the n media which are propagated through, and r is the distance of the electromagnetic wave propagation.
In one embodiment of the application, the computer program when executed by the processor further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000214
Figure GDA0003259773360000215
Acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA00032597733600002112
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA00032597733600002113
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000216
calculating the wave amplitude value of the electromagnetic signal of each frequency according to the following formula, wherein the wave amplitude value of the electromagnetic signal comprises the wave amplitude value E of electric field signals which are vertical to each otherRXAmplitude H of the sum magnetic field signalRX
Figure GDA0003259773360000217
Figure GDA0003259773360000218
Figure GDA0003259773360000219
Figure GDA00032597733600002110
Figure GDA00032597733600002111
Acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
generating a preset theoretical electromagnetic wave waveform curve, wherein the theoretical electromagnetic wave waveform curve comprises a theoretical magnetic field curve and a theoretical electric field curve which change in a sine rule;
comparing the wavelength value and the amplitude value of the intrinsic induction electromagnetic signal with the wavelength value and the amplitude value of the theoretical electromagnetic wave waveform curve, respectively;
obtaining an intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal;
acquiring the frequency f of the geological exploration result signal;
obtaining the electric field intensity E of the geological exploration result signalx
Obtaining the magnetic field intensity H of the geological exploration result signaly
Calculating the Carniya resistivity rho according to the following formulas
Figure GDA0003259773360000221
Wherein E isTxAmplitude, H, of electric field signals in the initial electromagnetic signal emitted by the source of the electromagnetic waveTxAmplitude, v, of a magnetic field signal in an initial electromagnetic signal emitted by said electromagnetic wave emission sourcenValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThe thickness value of the nth medium, alpha is attenuation coefficient, f is frequency, gamma is conductivity, gammanConductivity for propagation of electromagnetic waves in the nth medium, GRxFor amplification in the reception of signals, GTxFor amplification when transmitting signals, lambdaεIs the wavelength of an electromagnetic wave in a medium having a dielectric constant of epsilonnIs the dielectric constant, μ, of an electromagnetic wave in the nth mediumnTo be the permeability of the electromagnetic wave in the nth medium,
Figure GDA0003259773360000222
is the average value of the electrical conductivity of the electromagnetic wave in the n media that propagate through,
Figure GDA0003259773360000223
is the average value of the wavelength values of the electromagnetic wave in the n media propagating through,
Figure GDA0003259773360000224
is the average value of the magnetic permeability of the electromagnetic wave in the n media which are propagated through, and r is the distance of the electromagnetic wave propagation.
In one embodiment of the application, the computer program when executed by the processor further performs the steps of:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values;
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure GDA0003259773360000225
Figure GDA0003259773360000226
Acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure GDA0003259773360000227
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure GDA0003259773360000228
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure GDA0003259773360000231
calculating the wave amplitude value of the electromagnetic signal of each frequency according to the following formula, wherein the wave amplitude value of the electromagnetic signal comprises the wave amplitude value E of electric field signals which are vertical to each otherRXAmplitude H of the sum magnetic field signalRX
Figure GDA0003259773360000232
Figure GDA0003259773360000233
Figure GDA0003259773360000234
Figure GDA0003259773360000235
Figure GDA0003259773360000236
Acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
generating a preset theoretical electromagnetic wave waveform curve, wherein the theoretical electromagnetic wave waveform curve comprises a theoretical magnetic field curve and a theoretical electric field curve which change in a sine rule;
comparing the wavelength value and the amplitude value of the intrinsic induction electromagnetic signal with the wavelength value and the amplitude value of the theoretical electromagnetic wave waveform curve, respectively;
obtaining an intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal;
acquiring the frequency f of the geological exploration result signal;
obtaining the electric field intensity E of the geological exploration result signalx
Obtaining the magnetic field intensity H of the geological exploration result signaly
Calculating the Carniya resistivity rho according to the following formulas
Figure GDA0003259773360000237
Wherein E isTxAmplitude, H, of electric field signals in the initial electromagnetic signal emitted by the source of the electromagnetic waveTxAmplitude, v, of a magnetic field signal in an initial electromagnetic signal emitted by said electromagnetic wave emission sourcenValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThe thickness value of the nth medium, alpha is attenuation coefficient, f is frequency, gamma is conductivity, gammanConductivity for propagation of electromagnetic waves in the nth medium, GRxFor amplification in the reception of signals, GTxFor amplification when transmitting signals, lambdaεIs the wavelength of an electromagnetic wave in a medium having a dielectric constant of epsilonnIs electricityDielectric constant, mu, of magnetic wave in the n-th mediumnTo be the permeability of the electromagnetic wave in the nth medium,
Figure GDA0003259773360000238
is the average value of the electrical conductivity of the electromagnetic wave in the n media that propagate through,
Figure GDA0003259773360000241
is the average value of the wavelength values of the electromagnetic wave in the n media propagating through,
Figure GDA0003259773360000242
is the average value of the magnetic permeability of the electromagnetic wave in the n media which are propagated through, and r is the distance of the electromagnetic wave propagation.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory, among others. Non-volatile memory can include read-only memory (ROM), Programmable ROM (PROM), Electrically Programmable ROM (EPROM), Electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Link DRAM (SLDRAM), Rambus Direct RAM (RDRAM), direct bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
The technical features of the embodiments described above may be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the embodiments described above are not described, but should be considered as being within the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the claims. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (10)

1. A method of seismic survey detection, comprising:
acquiring electromagnetic signals with different frequencies, wherein the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
acquiring characteristic parameter values of the electromagnetic signals of all frequencies, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values, and the first time values of the electromagnetic signals of all frequencies are the sum of the propagation time of the electromagnetic signals from the emitting points to the test points and the propagation time of twice the propagation time of the electromagnetic signals from the induction points to the ground test points;
acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range;
and acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal.
2. The method of claim 1, wherein said obtaining values of characteristic parameters of said induced electromagnetic signal at each frequency comprises:
calculating the waviness value of the electromagnetic signal of each frequency according to the following formula
Figure FDA0003259773350000011
Figure FDA0003259773350000012
Obtaining the first time value based on the waviness value;
wherein v isnValue of propagation velocity in the nth medium through which said electromagnetic signal passes, hnThickness value of the nth medium.
3. The method of claim 2, wherein said obtaining values of characteristic parameters of said induced electromagnetic signal at each frequency further comprises:
acquiring a first distance S from the emission point of the geological exploration electromagnetic detection signal to the test point1
Obtaining a second distance S from the test point to a ground electromagnetic signal receiving point2
Calculating a first waviness value of the electromagnetic signal between the emission point of the geological exploration electromagnetic detection signal and the test point
Figure FDA0003259773350000013
Calculating a second waviness value of the electromagnetic signal between the test point and the ground electromagnetic signal receiving point
Figure FDA0003259773350000014
Calculating a first time value T of the electromagnetic signal at each frequency according to the following formula:
Figure FDA0003259773350000015
4. the method of claim 3, wherein said obtaining values of parameters characteristic of said induced electromagnetic signal at each frequency comprises:
calculating the wave amplitude value of the electromagnetic signal of each frequency according to the following formula, wherein the wave amplitude value of the electromagnetic signal comprises the wave amplitude value E of electric field signals which are vertical to each otherRXAmplitude H of the sum magnetic field signalRX
Figure FDA0003259773350000016
Figure FDA0003259773350000017
Figure FDA0003259773350000018
Figure FDA0003259773350000021
Figure FDA0003259773350000022
Figure FDA0003259773350000023
Wherein E isTxAmplitude, H, of electric field signals in the initial electromagnetic signal emitted by the source of the electromagnetic waveTxThe amplitude value of a magnetic field signal in an initial electromagnetic signal emitted by the electromagnetic wave emission source, alpha is an attenuation coefficient, f is frequency, gamma is conductivity, and gamma isnConductivity for propagation of electromagnetic waves in the nth medium, GRxFor amplification in the reception of signals, GTxFor amplification when transmitting signals, lambdaεIs the wavelength of an electromagnetic wave in a medium having a dielectric constant of epsilonnIs the dielectric constant, μ, of an electromagnetic wave in the nth mediumnTo be the permeability of the electromagnetic wave in the nth medium,
Figure FDA0003259773350000024
is the average value of the electrical conductivity of the electromagnetic wave in the n media that propagate through,
Figure FDA0003259773350000025
is the average value of the wavelength values of the electromagnetic wave in the n media propagating through,
Figure FDA0003259773350000026
is the average value of the magnetic permeability of the electromagnetic wave in the n media which are propagated through, and r is the distance of the electromagnetic wave propagation.
5. The method of any of claims 1-4, wherein said obtaining a geophysical survey pickup signal based on the intrinsic induced electromagnetic signal comprises:
fitting the intrinsic induction electromagnetic signal with a preset theoretical electromagnetic wave waveform curve;
and acquiring the intrinsic induction electromagnetic signal with the best fitting degree as a geological exploration result signal.
6. The method of claim 5, further comprising:
acquiring the frequency f of the geological exploration result signal;
obtaining the electric field intensity E of the geological exploration result signalx
Obtaining the magnetic field intensity H of the geological exploration result signaly
Calculating the Carniya resistivity rho according to the following formulas
Figure FDA0003259773350000027
7. The method of claim 5, wherein fitting the intrinsic induced electromagnetic signal to a preset theoretical electromagnetic wave waveform curve comprises:
generating a preset theoretical electromagnetic wave waveform curve, wherein the theoretical electromagnetic wave waveform curve comprises a theoretical magnetic field curve and a theoretical electric field curve which change in a sine rule;
comparing the wavelength value and the amplitude value of the intrinsic induction electromagnetic signal with the wavelength value and the amplitude value of the theoretical electromagnetic wave waveform curve, respectively.
8. A geophysical exploration wave detection device, comprising:
the device comprises an electromagnetic signal acquisition device, a signal processing device and a signal processing device, wherein the electromagnetic signal acquisition device is used for acquiring electromagnetic signals with different frequencies, and the electromagnetic signals are induction electromagnetic signals of geological exploration electromagnetic detection signals;
the tuning subsystem is used for acquiring characteristic parameter values of the electromagnetic signals of all frequencies and acquiring intrinsic induction electromagnetic signals of which the characteristic parameter values belong to a preset threshold range, wherein the characteristic parameter values comprise wavelength values, wave amplitude values and first time values, and the first time values of the electromagnetic signals of all frequencies are the sum of the propagation time of the electromagnetic signals from the transmitting point to the test point and the propagation time of twice the electromagnetic signals from the sensing point to the ground test point;
and the fitting subsystem is used for acquiring a geological exploration detection result signal based on the intrinsic induction electromagnetic signal.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method of any of claims 1 to 7 are implemented when the computer program is executed by the processor.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 7.
CN202010009649.5A 2020-01-06 2020-01-06 Geological exploration wave detection method, device, equipment and medium Active CN111123369B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010009649.5A CN111123369B (en) 2020-01-06 2020-01-06 Geological exploration wave detection method, device, equipment and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010009649.5A CN111123369B (en) 2020-01-06 2020-01-06 Geological exploration wave detection method, device, equipment and medium

Publications (2)

Publication Number Publication Date
CN111123369A CN111123369A (en) 2020-05-08
CN111123369B true CN111123369B (en) 2021-11-30

Family

ID=70486885

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010009649.5A Active CN111123369B (en) 2020-01-06 2020-01-06 Geological exploration wave detection method, device, equipment and medium

Country Status (1)

Country Link
CN (1) CN111123369B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104090302A (en) * 2014-07-02 2014-10-08 中国海洋石油总公司 Method for anomaly analysis of frequency domain of underground medium of work area
WO2015152955A1 (en) * 2014-04-03 2015-10-08 Halliburton Energy Services, Inc. Multi-component induction logging systems and methods using selected frequency inversion
CN105204073A (en) * 2015-09-18 2015-12-30 中南大学 Tensor apparent conductivity measurement method
CN108828681A (en) * 2018-08-01 2018-11-16 中国石油天然气集团有限公司 The determination method and apparatus of formation resistivity and polarizability
CN108983298A (en) * 2018-05-31 2018-12-11 湖南继善高科技有限公司 A kind of four-dimensional electromagnetic monitoring method and system in real time of oil gas pressure break

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9335434B2 (en) * 2012-11-02 2016-05-10 Pgs Geophysical As Method and system for processing data acquired in an electromagnetic survey
EP2985634B1 (en) * 2014-08-11 2019-02-27 Services Pétroliers Schlumberger Method and apparatus for determining resistivity of a formation

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015152955A1 (en) * 2014-04-03 2015-10-08 Halliburton Energy Services, Inc. Multi-component induction logging systems and methods using selected frequency inversion
CN104090302A (en) * 2014-07-02 2014-10-08 中国海洋石油总公司 Method for anomaly analysis of frequency domain of underground medium of work area
CN105204073A (en) * 2015-09-18 2015-12-30 中南大学 Tensor apparent conductivity measurement method
CN108983298A (en) * 2018-05-31 2018-12-11 湖南继善高科技有限公司 A kind of four-dimensional electromagnetic monitoring method and system in real time of oil gas pressure break
CN108828681A (en) * 2018-08-01 2018-11-16 中国石油天然气集团有限公司 The determination method and apparatus of formation resistivity and polarizability

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"强电磁干扰区大地电磁远参考技术试验效果分析";徐志敏 等;《物探与化探》;20180630;第42卷(第3期);第560-568页 *

Also Published As

Publication number Publication date
CN111123369A (en) 2020-05-08

Similar Documents

Publication Publication Date Title
Wang Frequencies of the Ricker wavelet
Arai et al. S-wave velocity profiling by inversion of microtremor H/V spectrum
Sinha et al. Spectral decomposition of seismic data with continuous-wavelet transform
Li et al. Prediction of fracture density using genetic algorithm support vector machine based on acoustic logging data
Perron et al. Robustness of kappa (κ) measurement in low‐to‐moderate seismicity areas: Insight from a site‐specific study in Provence, France
Luo et al. Analysis of group-velocity dispersion of high-frequency Rayleigh waves for near-surface applications
Dasios et al. Compressional‐wave Q estimation from full‐waveform sonic data
Zhou et al. Large-depth hydrogeological detection in the North China-type coalfield through short-offset grounded-wire TEM
Bergamo et al. On the relation between empirical amplification and proxies measured at Swiss and Japanese stations: Systematic regression analysis and neural network prediction of amplification
Wang The W transform
Xu et al. Method for calculating the fracture porosity of tight-fracture reservoirs
Hallal et al. An H/V geostatistical approach for building pseudo-3D Vs models to account for spatial variability in ground response analyses Part I: Model development
Akin et al. Site characterization using surface wave methods in the Arsin-Trabzon province, NE Turkey
Farrugia et al. Validation of 1D earthquake site characterization methods with observed earthquake site amplification in Alberta, Canada
Xue et al. Estimation of seismic quality factor in the time-frequency domain using variational mode decomposition
Xu et al. Fusion of multiple time‐domain GPR datasets of different center frequencies
Wang et al. Relational database for horizontal‐to‐vertical spectral ratios
Kimura et al. Systematic monitoring of instrumentation health in high-density broadband seismic networks
CN113703044B (en) Correction method and device for ancient river channel width, electronic equipment and storage medium
CN111123369B (en) Geological exploration wave detection method, device, equipment and medium
Liu et al. Quantitative analysis of near‐surface seismological complexity based on the generalized Lippmann–Schwinger matrix equation
Jiang et al. Hydrocarbon detection based on empirical mode decomposition, teager-kaiser energy, and the cepstrum
Tian et al. Seismic depositional sequence characterization by using enhanced multichannel variational-mode decomposition
Grégoire et al. GPR abilities for the detection and characterisation of open fractures in a salt mine
Fan et al. Inversion based on deep learning of logging-while-drilling directional resistivity measurements

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