CN112965120B - Transient electromagnetic signal processing method and device for geological exploration and storage medium - Google Patents
Transient electromagnetic signal processing method and device for geological exploration and storage medium Download PDFInfo
- Publication number
- CN112965120B CN112965120B CN202110136488.0A CN202110136488A CN112965120B CN 112965120 B CN112965120 B CN 112965120B CN 202110136488 A CN202110136488 A CN 202110136488A CN 112965120 B CN112965120 B CN 112965120B
- Authority
- CN
- China
- Prior art keywords
- signal
- secondary field
- point
- transient electromagnetic
- determining
- 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
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 89
- 238000003672 processing method Methods 0.000 title claims description 14
- 238000000034 method Methods 0.000 claims abstract description 22
- 230000011218 segmentation Effects 0.000 claims abstract description 17
- 238000005070 sampling Methods 0.000 claims description 43
- 238000004590 computer program Methods 0.000 claims description 9
- 230000005672 electromagnetic field Effects 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 abstract description 9
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/10—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices using induction coils
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a method, a device and a storage medium for processing transient electromagnetic signals for geological exploration, wherein the method comprises the following steps: acquiring a transient electromagnetic signal; determining a corresponding secondary field part signal according to the segmentation processing of the transient electromagnetic signal; selecting a corresponding secondary field ideal data segment according to the linearity of the secondary field part signal; performing linear fitting on ideal data segments of the secondary field under logarithmic coordinates; and determining an extended late signal of the secondary field part signal according to the linear fitting result of the ideal secondary field data segment, wherein the time starting point of the extended late signal is consistent with the time ending point of the secondary field part signal. In the inversion process of the transient electromagnetic method, the invention mainly adopts the data of the secondary field to analyze, and the ideal section of the secondary field signal is linearly fitted under the logarithmic coordinate, so that the signal-to-noise ratio of the signal can be greatly improved when the late signal is interfered or the signal is weaker, the algorithm model is simple, the calculated amount is less, and the calculation can be rapidly completed.
Description
Technical Field
The present invention relates to the field of signal processing technologies, and in particular, to a method and an apparatus for processing transient electromagnetic signals for geological exploration, and a storage medium.
Background
The transient electromagnetic method is a time domain electromagnetic method for solving the geological problems by using an ungrounded loop to emit a pulsed electromagnetic field into the ground and using a coil to observe the spatial and temporal distribution of the electromagnetic field generated by the underground eddy current induced by the pulsed electromagnetic field. The vortex alternating electromagnetic field with continuously changing structure and frequency in the conductive medium is a secondary field, the strength, direction, distribution rule and phase relation with the primary field are dependent on the electromagnetic property, size and shape of the conductive medium, and the conditions of the strength, frequency, inductive coupling relation between the primary field and the conductor and the like of the primary field.
The secondary field can be divided into an early signal stage and a late signal stage according to the transient process, and the signal of the secondary field is weak and is easy to interfere, so that the late signal is more difficult to capture. Common signal processing methods focus on denoising the secondary field, however, interference in the secondary field cannot be completely eliminated, when a late signal arrives in the secondary field, the signal attenuation is still obvious, and the acquisition of a complete and effective secondary field signal cannot be ensured, so that effective geological problem analysis cannot be performed. In summary, how to efficiently expand the secondary field signal is a problem to be solved.
Disclosure of Invention
In view of the foregoing, it is desirable to provide a method, apparatus and storage medium for processing transient electromagnetic signals for geological exploration, which solve the problem of how to efficiently expand the secondary field signals.
The invention provides a transient electromagnetic signal processing method for geological exploration, which comprises the following steps:
acquiring a transient electromagnetic signal;
determining a corresponding secondary field part signal according to the segmentation processing of the transient electromagnetic signal;
selecting a corresponding secondary field ideal data segment according to the linearity of the secondary field part signal;
performing linear fitting on the secondary field ideal data segment under logarithmic coordinates;
and determining an extended late signal of the secondary field part signal according to the linear fitting result of the secondary field ideal data segment, wherein the time starting point of the extended late signal is consistent with the time ending point of the secondary field part signal.
Further, the result of the segmentation processing of the transient electromagnetic signal sequentially includes a primary field part signal, a primary field and secondary field mixed signal, and the secondary field part signal, and determining the corresponding secondary field part signal according to the segmentation processing of the transient electromagnetic signal includes:
according to the linear slope difference of the transient electromagnetic signal under a logarithmic coordinate system, determining a first change point and a second change point of linear slope change;
determining the primary field part signal by taking the time starting point of the transient electromagnetic signal as the time starting point of the primary field part signal and the first change point as the time ending point of the primary field part signal;
the first change point is used as a time starting point of the primary field and secondary field mixed signal, and the second change point is used as a time end point of the primary field and secondary field mixed signal, so that the primary field and secondary field mixed signal is determined;
and determining the secondary field part signal by taking the second change point as a time starting point of the secondary field part signal and taking a time ending point of the transient electromagnetic signal as a time ending point of the primary field part signal.
Further, determining the first change point and the second change point of the slope change of the straight line according to the slope difference of the straight line of the transient electromagnetic signal in the logarithmic coordinate system comprises:
periodically sampling the transient electromagnetic signal to determine a plurality of sampling points;
determining the slope difference of the straight line corresponding to each sampling point according to the slope of the straight line between each sampling point and the adjacent sampling point;
comparing the slope differences of the straight lines corresponding to the sampling points, and determining the first change point and the second change point according to comparison results.
Further, determining the first change point and the second change point of the slope change of the straight line according to the slope difference of the straight line of the transient electromagnetic signal in the logarithmic coordinate system further comprises: and determining the first change point and the second change point according to a selection result of the manual work in the transient electromagnetic signal.
Further, the selecting the corresponding secondary field ideal data segment according to the linearity of the secondary field part signal includes:
periodically sampling the secondary field part signal to determine a plurality of sampling points;
taking each sampling point as a interception starting point, and determining each sampling point as an interception signal formed by the interception starting point according to a preset interception length;
calculating the linearity of each intercepted signal;
and determining the corresponding ideal data segment of the secondary field according to the intercepted signal with the minimum linearity.
Further, the linear fitting of the second order field ideal data segment in logarithmic coordinates includes:
converting the secondary field ideal data segment into a logarithmic coordinate, and determining a normalized data value corresponding to each data point in the secondary field ideal data segment;
and performing linear fitting according to the time starting point of the ideal data segment of the secondary field, the time ending point of the ideal data segment of the secondary field and the normalized data value, and determining a linear fitting equation under a logarithmic coordinate.
Further, the determining the extended late signal of the secondary field portion signal based on the linear fit of the secondary field ideal data segment comprises:
taking the time end point of the secondary field part signal as the time start point of the extended late signal;
determining a time end point of the extended late signal according to the time start point of the extended late signal and a preset extended duration;
and determining the extended late signal according to the time starting point of the extended late signal, the time ending point of the extended late signal and the linear fitting equation.
The invention also provides a transient electromagnetic signal processing device, comprising:
an acquisition unit for acquiring a transient electromagnetic signal;
the processing unit is used for determining a corresponding secondary field part signal according to the segmentation processing of the transient electromagnetic signal; the method is also used for selecting a corresponding secondary field ideal data segment according to the linearity of the secondary field part signal; the linear fitting method is also used for carrying out linear fitting on the ideal data segment of the secondary field under the logarithmic coordinate;
and the expansion unit is used for determining an expansion end signal of the secondary field part signal according to the linear fitting result of the secondary field ideal data segment, wherein the time starting point of the expansion end signal is consistent with the time ending point of the secondary field part signal.
The invention also provides a transient electromagnetic signal processing device, which comprises a processor and a memory, wherein the memory is stored with a computer program, and the computer program realizes the transient electromagnetic signal processing method for geological exploration when being executed by the processor.
The invention also provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a method of transient electromagnetic signal processing for geological exploration as described above.
Compared with the prior art, the invention has the beneficial effects that: firstly, effectively acquiring transient electromagnetic signals; then, segmenting by transient electromagnetic signals, and performing effective signal segmentation processing according to the distribution relation of the primary field and the secondary field so as to determine secondary field part signals in the signals; further, the linearity of the secondary field part signal is utilized to judge the corresponding linear relation, in general, the smaller the linearity is, the smaller the deviation between the description and the fitting straight line is, the corresponding secondary field ideal data segment is selected to determine and fit the signal segment with small deviation between the straight line and the signal segment with small interference; furthermore, under the logarithmic coordinates, linear fitting is carried out on an ideal section of the secondary field part signal, and when the late signal is interfered or the signal is weaker, the signal-to-noise ratio of the signal can be greatly improved; and finally, effectively expanding the secondary field part signal by utilizing the linear fitting result of the secondary field ideal data segment, avoiding that the secondary field part signal cannot acquire the related late signal due to too strong interference, effectively and accurately recovering the late signal by fitting the secondary field ideal data segment, and correspondingly expanding to obtain the complete signal information. In summary, the invention carries out linear fitting on the ideal section of the secondary field part signal under the logarithmic coordinate, and can greatly improve the signal-to-noise ratio of the signal when the late signal is interfered or the signal is weaker; in addition, the linear fitting is carried out on the ideal section of the secondary field part signal, the algorithm model is simple, the calculated amount is small, the calculation can be rapidly completed, the algorithm complexity is reduced, the extension of the late signal of the secondary field part signal is ensured to be rapidly and accurately carried out, and the integrity and the high efficiency of the transient electromagnetic signal are greatly improved. In addition, the transient electromagnetic signal processing method for geological exploration provided by the invention has simple algorithm, has no strict requirement on the computing capacity of hardware equipment, can be used for a PC (personal computer) and embedded system equipment with weaker computing power, and greatly improves the practicability of the algorithm.
Drawings
FIG. 1 is a schematic flow chart of a transient electromagnetic signal processing method for geological exploration;
FIG. 2 is a schematic diagram of a process for determining a secondary field portion signal according to the present invention;
FIG. 3 is a schematic flow chart of determining a first change point and a second change point according to the present invention;
FIG. 4 is a schematic flow chart of selecting ideal data segments of a secondary field according to the present invention;
FIG. 5 is a schematic flow chart of linear fitting of ideal data segments of a secondary field according to the present invention;
FIG. 6 is a flow chart of determining an extended late signal according to the present invention;
FIG. 7 is a signal diagram of a signal segment according to the present invention;
FIG. 8 is a signal diagram of an extended late signal according to the present invention;
fig. 9 is a schematic structural diagram of a transient electromagnetic signal processing device provided by the invention.
Detailed Description
Preferred embodiments of the present invention will now be described in detail with reference to the accompanying drawings, which form a part hereof, and together with the description serve to explain the principles of the invention, and are not intended to limit the scope of the invention.
Example 1
The embodiment of the invention provides a transient electromagnetic signal processing method for geological exploration, and referring to fig. 1, fig. 1 is a schematic flow chart of the transient electromagnetic signal processing method for geological exploration, which is provided by the invention, and comprises steps S1 to S5, wherein:
in step S1, a transient electromagnetic signal is acquired;
in step S2, determining a corresponding secondary field part signal according to the segmentation processing of the transient electromagnetic signal;
in step S3, selecting a corresponding secondary field ideal data segment according to the linearity of the secondary field part signal;
in step S4, linear fitting is carried out on ideal data segments of the secondary field under logarithmic coordinates;
in step S5, an extended late signal of the secondary field portion signal is determined according to the linear fitting result of the secondary field ideal data segment, wherein the time start point of the extended late signal coincides with the time end point of the secondary field portion signal.
In the embodiment of the invention, firstly, effective acquisition of transient electromagnetic signals is carried out; then, segmenting by transient electromagnetic signals, and performing effective signal segmentation processing according to the distribution relation of the primary field and the secondary field so as to determine secondary field part signals in the signals; further, the linearity of the secondary field part signal is utilized to judge the corresponding linear relation, in general, the smaller the linearity is, the smaller the deviation between the description and the fitting straight line is, the corresponding secondary field ideal data segment is selected to determine and fit the signal segment with small deviation between the straight line and the signal segment with small interference; furthermore, under the logarithmic coordinates, linear fitting is carried out on an ideal section of the secondary field signal, and when the late signal is interfered or the signal is weak, the signal-to-noise ratio of the signal can be greatly improved; and finally, effectively expanding the secondary field part signal by utilizing the linear fitting result of the secondary field ideal data segment, avoiding that the secondary field part signal cannot acquire the related late signal due to too strong interference, effectively and accurately recovering the late signal by fitting the secondary field ideal data segment, and correspondingly expanding to obtain the complete signal information.
Preferably, as seen in fig. 2, fig. 2 is a schematic flow chart of determining a secondary field part signal according to the present invention, and the step S2 includes steps S21 to S24, wherein:
in step S21, according to the difference of the slope of the line of the transient electromagnetic signal in the logarithmic coordinate system, determining a first change point and a second change point of the slope change of the line;
in step S22, the time start point of the transient electromagnetic signal is taken as the time start point of the primary field part signal, the first change point is taken as the time end point of the primary field part signal, and the primary field part signal is determined;
in step S23, the first change point is used as a time start point of the mixed signal of the primary field and the secondary field, and the second change point is used as a time end point of the mixed signal of the primary field and the secondary field, so as to determine the mixed signal of the primary field and the secondary field;
in step S24, the second change point is used as a time start point of the secondary field portion signal, and the time end point of the transient electromagnetic signal is used as a time end point of the primary field portion signal, so as to determine the secondary field portion signal.
As a specific embodiment, the first change point and the second change point with the largest change of the slope of the straight line are determined according to the slope difference of the straight line, the transient electromagnetic signal is divided into three sections by taking the first change point and the second change point as demarcation points, and therefore the primary field part signal, the primary field and secondary field mixed signal and the secondary field part signal are effectively distinguished.
Preferably, as seen in connection with fig. 3, fig. 3 is a schematic flow chart of determining a first change point and a second change point according to the present invention, and step S21 further includes steps S211 to S213, where:
in step S211, periodically sampling the transient electromagnetic signal to determine a plurality of sampling points;
in step S212, determining a difference of the slope of the line corresponding to each sampling point according to the slope of the line between each sampling point and the adjacent sampling points;
in step S213, the slope differences of the straight lines corresponding to the sampling points are compared, and the first change point and the second change point are determined according to the comparison result.
As a specific embodiment, the present embodiment compares the slope of a straight line by using a plurality of sampling points, and when the slope change between two adjacent sampling points is larger, that is, the slope difference of the straight line is larger, the signal change points corresponding to the sampling points, that is, the first change point and the second change point, are described, and the signal change is judged according to the slope of the straight line, so as to distinguish 3 sections, and form corresponding primary field part signals, primary field and secondary field mixed signals and secondary field part signals.
Preferably, the determining the first change point and the second change point further includes: and determining a first change point and a second change point according to a selection result of the human in the transient electromagnetic signal. As a specific embodiment, the first change point and the second change point can be manually determined according to the difference of the linear slope of the transient electromagnetic signal under the logarithmic coordinates, so that 3 sections are distinguished, and the time starting points and the end points of the 3 sections are determined, namely, a primary field part signal, a primary field and secondary field mixed signal and a secondary field part signal are formed.
Preferably, as seen in fig. 4, fig. 4 is a schematic flow chart of selecting ideal data segments of a secondary field according to the present invention, and step S3 further includes steps S31 to S34, wherein:
in step S31, the secondary field part signal is periodically sampled to determine a plurality of sampling points;
in step S32, each sampling point is taken as a interception start point, and according to a preset interception length, an interception signal formed by each sampling point corresponding to the interception start point is determined;
in step S33, the linearity of each intercepted signal is calculated;
in step S34, a corresponding secondary field ideal data segment is determined according to the truncated signal with the minimum linearity.
As a specific embodiment, the example circularly calculates the linearity of the intercepted signals formed correspondingly by taking a plurality of different sampling points as starting points, and the intercepted signals with the minimum linearity, namely the signal section with the minimum error with the fitting straight line, have the minimum linearity, so that the condition that the intercepted signals of the section are least and most ideal in interference degree is illustrated, and the intercepted signals are correspondingly expanded as ideal data sections of a secondary field, so that the effectiveness of the later-stage signals of the subsequent expansion is ensured.
Preferably, as seen in conjunction with fig. 5, fig. 5 is a schematic flow chart of linear fitting of the second-field ideal data segment according to the present invention, and step S4 includes steps S41 to S42, where:
in step S41, converting the secondary field ideal data segment to a logarithmic coordinate, and determining a normalized data value corresponding to each data point in the secondary field ideal data segment;
in step S42, a linear fitting equation in logarithmic coordinates is determined by performing linear fitting based on the time start point of the ideal data segment of the secondary field, the time end point of the ideal data segment of the secondary field, and the normalized data value.
As a specific embodiment, the present example expands the ideal data segment of the secondary field under the logarithmic coordinate, greatly improves the signal-to-noise ratio of the signal by using the logarithmic coordinate, and simultaneously reduces the data redundancy by using the normalized data value of each data point constituting the ideal data segment of the secondary field in the logarithmic coordinate system, thereby performing effective linear fitting, and determining the linear fitting equation correspondingly formed by the ideal data segment of the secondary field. It should be noted that, in the embodiment of the present invention, the fitting manner of determining the linear fitting equation by the normalized data values between the time start point and the time end point of the ideal data segment of the secondary field includes, but is not limited to, polynomial curve fitting, gaussian fitting, and power exponent fitting.
Preferably, as seen in fig. 6, fig. 6 is a schematic flow chart of determining an extended late signal according to the present invention, and step S5 includes steps S51 to S53, wherein:
in step S51, the time end of the secondary field part signal is used as the time start of the extended late signal;
in step S52, determining a time end point of the extended late signal according to the time start point of the extended late signal and a preset extended duration;
in step S53, the extended late signal is determined according to the time start of the extended late signal, the time end of the extended late signal, and the linear fit equation.
As a specific embodiment, the time starting point of the extended late signal and the time ending point of the extended late signal are brought into a linear fitting equation, so that the extended late signal corresponding to the secondary field part signal can be effectively determined, and the signal integrity is ensured.
In a specific embodiment of the present invention, as shown in fig. 7 and 8, fig. 7 is a signal diagram of signal segmentation provided by the present invention, and fig. 8 is a signal diagram of extended late signal provided by the present invention, where in fig. 7, a signal received after transmitting a pulsed electromagnetic field into the ground is segmented, a first change point is a signal point at 1000s, a second change point is a signal point at 2200s, and the first change point and the second change point are taken as boundaries, and are divided into 3 intervals, namely, the 1 st segment: primary field part signal, section 2: primary and secondary field mixed signal, section 3: and selecting the 4 th section with the minimum linearity from the secondary field part signals of the 3 rd section: the secondary field is an ideal data segment. In fig. 8, after selecting the ideal data segment of the secondary field (the 4 th segment), performing linear fitting on the ideal data segment of the secondary field (corresponding to the ideal secondary signal in fig. 8) to obtain a corresponding linear fitting equation, using the time end point of the ideal data segment of the secondary field as the time start point of the extended late signal, determining the time end point of the extended late signal according to the preset duration of the extended late signal, and finally bringing the time start point and the time end point of the extended late signal into the linear fitting equation to effectively determine any signal point of the extended late signal. After the extended late signal in fig. 8 is obtained, a more complete analysis and judgment of the geological condition is made by the signal.
Example 2
An embodiment of the present invention provides a transient electromagnetic signal processing device, and as shown in fig. 9, fig. 9 is a schematic structural diagram of the transient electromagnetic signal processing device provided by the present invention, where the transient electromagnetic signal processing device 900 includes:
an acquisition unit 901 for acquiring a transient electromagnetic signal;
a processing unit 902, configured to determine a corresponding secondary field portion signal according to the segmentation processing of the transient electromagnetic signal; the method is also used for selecting a corresponding secondary field ideal data segment according to the linearity of the secondary field part signal; the method is also used for carrying out linear fitting on the ideal data segment of the secondary field under the logarithmic coordinates;
and the expansion unit 903 is configured to determine an expanded late signal of the secondary field part signal according to a linear fitting result of the secondary field ideal data segment, where a time start point of the expanded late signal is consistent with a time end point of the secondary field part signal.
Example 3
The embodiment of the invention provides a transient electromagnetic signal processing device, which comprises a processor and a memory, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the transient electromagnetic signal processing method for geological exploration is realized.
Example 4
An embodiment of the present invention provides a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements a transient electromagnetic signal processing method for geological exploration as described above.
The invention discloses a transient electromagnetic signal processing method, a device and a storage medium for geological exploration, wherein firstly, effective acquisition of a transient electromagnetic signal is carried out; then, segmenting by transient electromagnetic signals, and performing effective signal segmentation processing according to the distribution relation of the primary field and the secondary field so as to determine secondary field part signals in the signals; further, the linearity of the secondary field part signal is utilized to judge the corresponding linear relation, in general, the smaller the linearity is, the smaller the deviation between the description and the fitting straight line is, the corresponding secondary field ideal data segment is selected to determine and fit the signal segment with small deviation between the straight line and the signal segment with small interference; furthermore, under the logarithmic coordinates, linear fitting is carried out on an ideal section of the secondary field part signal, and when the late signal is interfered or the signal is weaker, the signal-to-noise ratio of the signal can be greatly improved; and finally, effectively expanding the secondary field part signal by utilizing the linear fitting result of the secondary field ideal data segment, avoiding that the secondary field part signal cannot acquire the related late signal due to too strong interference, effectively and accurately recovering the late signal by fitting the secondary field ideal data segment, and correspondingly expanding to obtain the complete signal information.
According to the technical scheme, linear fitting is carried out on an ideal section of the secondary field part signal under a logarithmic coordinate, and when the late signal is interfered or the signal is weak, the signal-to-noise ratio of the signal can be greatly improved; in addition, the linear fitting is carried out on the ideal section of the secondary field part signal, the algorithm model is simple, the calculated amount is small, the calculation can be rapidly completed, the algorithm complexity is reduced, the extension of the late signal of the secondary field part signal is ensured to be rapidly and accurately carried out, and the integrity and the high efficiency of the transient electromagnetic signal are greatly improved. In addition, the transient electromagnetic signal processing method for geological exploration provided by the invention has simple algorithm, has no strict requirement on the computing capacity of hardware equipment, can be used for a PC (personal computer) and embedded system equipment with weaker computing power, and greatly improves the practicability of the algorithm.
The present invention is not limited to the above-mentioned embodiments, and any changes or substitutions that can be easily understood by those skilled in the art within the technical scope of the present invention are intended to be included in the scope of the present invention.
Claims (7)
1. A method of processing transient electromagnetic signals for geological exploration, comprising:
acquiring a transient electromagnetic signal, wherein the transient electromagnetic signal is a signal received after a pulse electromagnetic field is emitted to the ground;
determining a corresponding secondary field part signal according to the segmentation processing of the transient electromagnetic signal;
selecting a corresponding secondary field ideal data segment according to the linearity of the secondary field part signal;
performing linear fitting on the secondary field ideal data segment under logarithmic coordinates;
determining an extended late signal of the secondary field partial signal according to a linear fitting result of the secondary field ideal data segment, wherein the time starting point of the extended late signal is consistent with the time ending point of the secondary field ideal data segment;
the result of the segmentation processing of the transient electromagnetic signal sequentially comprises a primary field part signal, a primary field and secondary field mixed signal and a secondary field part signal, and the determining of the corresponding secondary field part signal according to the segmentation processing of the transient electromagnetic signal comprises the following steps:
according to the linear slope difference of the transient electromagnetic signal under a logarithmic coordinate system, determining a first change point and a second change point of linear slope change;
determining the primary field part signal by taking the time starting point of the transient electromagnetic signal as the time starting point of the primary field part signal and the first change point as the time ending point of the primary field part signal;
the first change point is used as a time starting point of the primary field and secondary field mixed signal, and the second change point is used as a time end point of the primary field and secondary field mixed signal, so that the primary field and secondary field mixed signal is determined;
determining the secondary field part signal by taking the second change point as a time starting point of the secondary field part signal and taking a time end point of the transient electromagnetic signal as a time end point of the primary field part signal;
the determining a first change point and a second change point of the slope change of the straight line according to the slope difference of the straight line of the transient electromagnetic signal under a logarithmic coordinate system comprises:
periodically sampling the transient electromagnetic signal to determine a plurality of sampling points;
determining the slope difference of the straight line corresponding to each sampling point according to the slope of the straight line between each sampling point and the adjacent sampling point;
comparing the slope differences of the straight lines corresponding to each sampling point, and determining the first change point and the second change point according to comparison results;
the selecting the corresponding ideal data segment of the secondary field according to the linearity of the secondary field part signal comprises:
periodically sampling the secondary field part signal to determine a plurality of sampling points;
taking each sampling point as a interception starting point, and determining each sampling point as an interception signal formed by the interception starting point according to a preset interception length;
calculating the linearity of each intercepted signal;
and determining the corresponding ideal data segment of the secondary field according to the intercepted signal with the minimum linearity.
2. The method for processing transient electromagnetic signals for geological exploration according to claim 1, wherein determining a first change point and a second change point of a change in slope of a straight line according to a difference in slope of the straight line of the transient electromagnetic signals in a logarithmic coordinate system further comprises: and determining the first change point and the second change point according to a selection result of the manual work in the transient electromagnetic signal.
3. The method of transient electromagnetic signal processing for geological exploration of claim 1, wherein said linearly fitting said secondary field ideal data segment in logarithmic coordinates comprises:
converting the secondary field ideal data segment into a logarithmic coordinate, and determining a normalized data value corresponding to each data point in the secondary field ideal data segment;
and performing linear fitting according to the time starting point of the ideal data segment of the secondary field, the time ending point of the ideal data segment of the secondary field and the normalized data value, and determining a linear fitting equation under a logarithmic coordinate.
4. A method of transient electromagnetic signal processing for geological exploration according to claim 3, wherein said determining an extended late signal of said secondary field portion signal from a linear fit of said secondary field ideal data segment comprises:
taking the time end point of the secondary field ideal data segment as the time start point of the extended late signal;
determining a time end point of the extended late signal according to the time start point of the extended late signal and a preset extended duration;
and determining the extended late signal according to the time starting point of the extended late signal, the time ending point of the extended late signal and the linear fitting equation.
5. A transient electromagnetic signal processing apparatus, comprising:
an acquisition unit for acquiring a transient electromagnetic signal;
the processing unit is used for determining a corresponding secondary field part signal according to the segmentation processing of the transient electromagnetic signal; the method is also used for selecting a corresponding secondary field ideal data segment according to the linearity of the secondary field part signal; the linear fitting method is also used for carrying out linear fitting on the ideal data segment of the secondary field under the logarithmic coordinate;
the expansion unit is used for determining an expansion late signal of the secondary field partial signal according to the linear fitting result of the secondary field ideal data segment, wherein the time starting point of the expansion late signal is consistent with the time ending point of the secondary field ideal data segment;
the result of the segmentation processing of the transient electromagnetic signal sequentially comprises a primary field part signal, a primary field and secondary field mixed signal and a secondary field part signal, and the determining of the corresponding secondary field part signal according to the segmentation processing of the transient electromagnetic signal comprises the following steps:
according to the linear slope difference of the transient electromagnetic signal under a logarithmic coordinate system, determining a first change point and a second change point of linear slope change;
determining the primary field part signal by taking the time starting point of the transient electromagnetic signal as the time starting point of the primary field part signal and the first change point as the time ending point of the primary field part signal;
the first change point is used as a time starting point of the primary field and secondary field mixed signal, and the second change point is used as a time end point of the primary field and secondary field mixed signal, so that the primary field and secondary field mixed signal is determined;
determining the secondary field part signal by taking the second change point as a time starting point of the secondary field part signal and taking a time end point of the transient electromagnetic signal as a time end point of the primary field part signal;
the determining a first change point and a second change point of the slope change of the straight line according to the slope difference of the straight line of the transient electromagnetic signal under a logarithmic coordinate system comprises:
periodically sampling the transient electromagnetic signal to determine a plurality of sampling points;
determining the slope difference of the straight line corresponding to each sampling point according to the slope of the straight line between each sampling point and the adjacent sampling point;
comparing the slope differences of the straight lines corresponding to each sampling point, and determining the first change point and the second change point according to comparison results;
the selecting the corresponding ideal data segment of the secondary field according to the linearity of the secondary field part signal comprises:
periodically sampling the secondary field part signal to determine a plurality of sampling points;
taking each sampling point as a interception starting point, and determining each sampling point as an interception signal formed by the interception starting point according to a preset interception length;
calculating the linearity of each intercepted signal;
and determining the corresponding ideal data segment of the secondary field according to the intercepted signal with the minimum linearity.
6. A transient electromagnetic signal processing device comprising a processor and a memory, said memory having stored thereon a computer program which, when executed by said processor, implements a transient electromagnetic signal processing method for geological exploration as claimed in any one of claims 1-4.
7. A computer-readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements a method for transient electromagnetic signal processing for geological exploration according to any of claims 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110136488.0A CN112965120B (en) | 2021-02-01 | 2021-02-01 | Transient electromagnetic signal processing method and device for geological exploration and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110136488.0A CN112965120B (en) | 2021-02-01 | 2021-02-01 | Transient electromagnetic signal processing method and device for geological exploration and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112965120A CN112965120A (en) | 2021-06-15 |
| CN112965120B true CN112965120B (en) | 2024-04-12 |
Family
ID=76273568
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110136488.0A Active CN112965120B (en) | 2021-02-01 | 2021-02-01 | Transient electromagnetic signal processing method and device for geological exploration and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112965120B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102176064A (en) * | 2011-02-14 | 2011-09-07 | 福州华虹智能科技开发有限公司 | Method for eliminating electromagnetic interference signals of underground coalmine transient electromagnetic instrument |
| CN105589108A (en) * | 2015-12-14 | 2016-05-18 | 中国科学院电子学研究所 | Rapid three-dimensional inversion method for transient electromagnetism based on different constraint conditions |
| CN107621656A (en) * | 2017-08-31 | 2018-01-23 | 湖南五维地质科技有限公司 | Window method when taking out of transient electromagnetic method |
| CN109350098A (en) * | 2018-08-27 | 2019-02-19 | 苏州瑞派宁科技有限公司 | Determination method, method for reconstructing and the device of the fit approach of signal |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11122992B2 (en) * | 2018-08-14 | 2021-09-21 | Biosense Webster (Israel) Ltd. | Determining a location and an orientation of a magnetic field sensor |
-
2021
- 2021-02-01 CN CN202110136488.0A patent/CN112965120B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102176064A (en) * | 2011-02-14 | 2011-09-07 | 福州华虹智能科技开发有限公司 | Method for eliminating electromagnetic interference signals of underground coalmine transient electromagnetic instrument |
| CN105589108A (en) * | 2015-12-14 | 2016-05-18 | 中国科学院电子学研究所 | Rapid three-dimensional inversion method for transient electromagnetism based on different constraint conditions |
| CN107621656A (en) * | 2017-08-31 | 2018-01-23 | 湖南五维地质科技有限公司 | Window method when taking out of transient electromagnetic method |
| CN109350098A (en) * | 2018-08-27 | 2019-02-19 | 苏州瑞派宁科技有限公司 | Determination method, method for reconstructing and the device of the fit approach of signal |
Non-Patent Citations (1)
| Title |
|---|
| 基于趋势面分析的瞬变电磁弱异常提取方法;程建远;李博凡;范涛;刘磊;;煤炭学报;20151215;第40卷(第12期);第2856-2864页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112965120A (en) | 2021-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105044793B (en) | A kind of inversion method and device of multiple tracks transient electromagnetic detecting data | |
| US20110037473A1 (en) | Data acquisition method with a three dimensional small bin electromagnetic consecutive array | |
| CN105353373B (en) | One kind is based on Hough transform Ground Penetrating Radar target extraction method and device | |
| CN106772632B (en) | A kind of tunnel advance geology exploration method based on delay transient electromagnetic | |
| CN103941254A (en) | Soil physical property classification recognition method and device based on geological radar | |
| CN113275341B (en) | Pipe cleaner tracking and positioning method based on distributed optical fiber vibration sensing | |
| CN105223481A (en) | Based on the partial-discharge ultrahigh-frequency signal initial time defining method of difference energy function | |
| US20050270037A1 (en) | Method and system for non-destructive evaluation of conducting structures | |
| CN104143115A (en) | Technological method for achieving soil water content classified identification through geological radar technology | |
| CN106137184B (en) | Electrocardiosignal QRS complex detection method based on wavelet transformation | |
| CN112965120B (en) | Transient electromagnetic signal processing method and device for geological exploration and storage medium | |
| CN108459353B (en) | Faint magnetic resonance signal extracting method and device under a kind of electromagnetic noise background | |
| CN111505465A (en) | Method and device for building database, looking up table and positioning local discharge characteristics of winding equipment | |
| CN115293219B (en) | Pulse signal denoising method integrating wavelet and kurtosis | |
| US11300603B2 (en) | Method and system for condition monitoring electrical equipment | |
| CN113608190B (en) | Sea surface target detection method and system based on three characteristics of singular space | |
| CN110927781A (en) | Seismic data optimization method and system for high-precision seismic coherence attribute | |
| CN104820213B (en) | A kind of search and track radar directional diagram method of estimation rejected based on dynamic | |
| CN109061379B (en) | Grounding grid topological structure based on induction voltage differential method and breakpoint identification method | |
| CN112363019A (en) | Grounding grid topological structure fault detection method and device | |
| JPH0843539A (en) | Processing method for received signal used for detection of buried object | |
| CN111177886B (en) | Marine ranging planning and soil thickness prediction method based on geophysical prospecting data analysis | |
| CN109520493A (en) | A kind of earth magnetism filtering method based on multiple dimensioned empirical mode decomposition | |
| CN118131229B (en) | Method for detecting geological radar in hole capable of distinguishing abnormal azimuth | |
| CN117872039B (en) | Line fault location method and system based on improved RBF network |
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 |