A time-domain airborne electromagnetic data filtering method, device and electronic equipment

Abstract

The application discloses a time domain airborne electromagnetic data filtering method and device and electronic equipment. The specific scheme is as follows: after a preset electromagnetic transmitting device sends an electromagnetic signal to the underground medium of a target exploration area, response electromagnetic data fed back by the underground medium based on the electromagnetic signal is acquired according to a preset sampling frequency and a preset sampling time length, wherein the response electromagnetic data is used for representing the electrical characteristic of the underground medium; filter parameters corresponding to the response electromagnetic data are determined according to the electromagnetic signal and the preset sampling frequency; and the response electromagnetic data is filtered based on the filter parameters to obtain target electromagnetic data corresponding to the response electromagnetic data. The application realizes effective filtering processing of the response electromagnetic data by the filter parameters determined according to the characteristics of the emitted electromagnetic signal, guarantees the accuracy of the target electromagnetic data, and improves the data signal-to-noise ratio.

Description

Technical Field

[0001] This invention relates to the field of data processing technology, and in particular to a time-domain airborne electromagnetic data filtering method, apparatus, and electronic device. Background Technology

[0002] Time-domain airborne electromagnetic exploration involves transmitting periodic electromagnetic signals into the subsurface and receiving the electromagnetic responses of the secondary fields from the subsurface medium. The electrical characteristics of the subsurface medium are then determined through inversion processing of the received electromagnetic waves. However, the received electromagnetic wave data may be affected by various noise sources, such as aircraft noise, instrument noise, or human noise, which can impact the reliability of the subsequent inversion results.

[0003] In existing technologies, filtering of time-domain airborne electromagnetic data is mostly based on the characteristics of noise. That is, different filtering methods are used for different types of noise. For example, notch comb filtering is used to remove power frequency interference, polynomial fitting or high-pass filtering is used to remove motion noise, and low-pass filtering is used to remove high-frequency noise. However, when there are many types of noise with different characteristics, the above filtering methods are not only complex, but also the simultaneous use of multiple filtering methods can easily have an adverse effect on the effective signal, which is not conducive to subsequent inversion processing. Summary of the Invention

[0004] This invention provides a time-domain airborne electromagnetic data filtering method, apparatus, and electronic device. By determining filter parameters based on the characteristics of the transmitted electromagnetic signal, it achieves effective filtering of the response electromagnetic data, ensuring the accuracy of the target electromagnetic data and improving the data signal-to-noise ratio.

[0005] According to one aspect of the present invention, a time-domain airborne electromagnetic data filtering method is provided, the method comprising:

[0006] After sending electromagnetic signals to the underground medium in the target exploration area through a preset electromagnetic transmission device, the response electromagnetic data of the underground medium based on the electromagnetic signals is obtained according to the preset sampling frequency and preset sampling duration. The response electromagnetic data is used to characterize the electrical characteristics of the underground medium.

[0007] Based on the electromagnetic signal and the preset sampling frequency, determine the filter parameters corresponding to the response electromagnetic data;

[0008] The response electromagnetic data is filtered based on the filter parameters to obtain the target electromagnetic data corresponding to the response electromagnetic data.

[0009] According to another aspect of the present invention, a time-domain airborne electromagnetic data filtering device is provided, the device comprising:

[0010] The electromagnetic data acquisition module is used to acquire the response electromagnetic data of the underground medium based on the electromagnetic signal after sending an electromagnetic signal to the underground medium in the target exploration area through a preset electromagnetic transmission device, according to a preset sampling frequency and a preset sampling duration. The response electromagnetic data is used to characterize the electrical characteristics of the underground medium.

[0011] The filter parameter determination module is used to determine the filter parameters corresponding to the response electromagnetic data based on the electromagnetic signal and the preset sampling frequency.

[0012] The electromagnetic data filtering module is used to filter the response electromagnetic data based on filter parameters to obtain the target electromagnetic data corresponding to the response electromagnetic data.

[0013] According to another aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0014] At least one processor; and

[0015] A memory that is communicatively connected to at least one processor; wherein,

[0016] The memory stores a computer program that can be executed by at least one processor, such that the at least one processor is able to perform the time-domain airborne electromagnetic data filtering method according to any embodiment of the present invention.

[0017] According to another aspect of the present invention, a computer-readable storage medium is provided that stores computer instructions for causing a processor to execute and implement the time-domain airborne electromagnetic data filtering method of any embodiment of the present invention.

[0018] According to another aspect of the present invention, a computer program product is provided, comprising a computer program, characterized in that the computer program, when executed by a processor, implements a time-domain airborne electromagnetic data filtering method as described in any embodiment of the present invention.

[0019] The technical solution of this invention involves transmitting electromagnetic signals to the underground medium of a target exploration area via a preset electromagnetic transmitting device. Then, based on a preset sampling frequency and a preset sampling duration, the response electromagnetic data fed back by the underground medium is acquired. Since the fed-back response electromagnetic data is susceptible to various noises, filter parameters corresponding to the response electromagnetic data can be determined based on the transmitted electromagnetic signal and the preset sampling frequency. These filter parameters are then used to filter the response electromagnetic data, yielding the target electromagnetic data corresponding to the response electromagnetic data. By determining the filter parameters based on the characteristics of the transmitted electromagnetic signal, the filtering process is simplified, achieving effective filtering of the response electromagnetic data. This amplifies the useful signal and suppresses noise signals, ensuring the accuracy of the target electromagnetic data, improving the data signal-to-noise ratio, and guaranteeing the reliability and accuracy of subsequent inversion processing based on the target electromagnetic data.

[0020] It should be understood that the description in this section is not intended to identify key or essential features of the embodiments of the present invention, nor is it intended to limit the scope of the invention. Other features of the invention will become readily apparent from the following description. Attached Figure Description

[0021] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0022] Figure 1 This is a flowchart of a time-domain airborne electromagnetic data filtering method provided in an embodiment of the present invention;

[0023] Figure 2 This is a flowchart of a time-domain airborne electromagnetic data filtering method provided in an embodiment of the present invention;

[0024] Figure 3 This is a schematic diagram of the structure of a time-domain airborne electromagnetic data filtering device provided in an embodiment of the present invention;

[0025] Figure 4 This is a schematic diagram of the structure of an electronic device that implements the time-domain airborne electromagnetic data filtering method according to an embodiment of the present invention. Detailed Implementation

[0026] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0027] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0028] Example 1

[0029] Figure 1 This is a flowchart of a time-domain airborne electromagnetic data filtering method provided in Embodiment 1 of the present invention. This embodiment is applicable to situations where filter parameters are determined based on the characteristics of the transmitted electromagnetic signal to filter the response electromagnetic data. This method can be executed by a time-domain airborne electromagnetic data filtering device, which can be implemented in hardware and / or software. This time-domain airborne electromagnetic data filtering device can be configured in electronic devices such as mobile phones, computers, or servers. Figure 1 As shown, the method includes:

[0030] S110. After sending an electromagnetic signal to the underground medium in the target exploration area through a preset electromagnetic transmission device, the response electromagnetic data of the underground medium based on the electromagnetic signal is obtained according to the preset sampling frequency and preset sampling duration.

[0031] The electromagnetic response data is used to characterize the electrical properties of the subsurface medium. The preset electromagnetic transmitting device can be a pre-set device used to transmit electromagnetic signals. Since this embodiment of the invention is applied in a time-domain airborne electromagnetic exploration scenario, the preset electromagnetic transmitting device can be an electromagnetic transmitting coil pre-deployed on an aircraft. The electromagnetic response data is the time-domain airborne electromagnetic data. The target exploration area can be the area currently requiring exploration processing. The subsurface medium can be understood as the subsurface medium body of the target exploration area. The electromagnetic signal can be understood as a signal used for geological exploration of the target exploration area, i.e., a signal used to detect the electrical properties of the subsurface medium.

[0032] Since the underground medium generates an induced electromagnetic field after an electromagnetic signal is sent to it, the response electromagnetic data corresponding to the induced electromagnetic field can be acquired according to a preset sampling frequency and a preset sampling duration. That is, the preset sampling frequency can be understood as the pre-set frequency for acquiring the response electromagnetic data. The preset sampling duration can be understood as the duration for acquiring the response electromagnetic data. The response electromagnetic data can be understood as the induced signal corresponding to the induced electromagnetic field generated by the electromagnetic signal in the underground medium. The response electromagnetic data is used to characterize the electrical characteristics of the underground medium. Electrical characteristics may include at least one of the following: electrical conductivity, magnetic permeability, and polarizability of the underground medium.

[0033] Specifically, after sending electromagnetic signals to the underground medium of the target exploration area through a preset electromagnetic transmission device, the response electromagnetic data corresponding to the electromagnetic signals fed back by the underground medium can be obtained according to the preset sampling frequency and preset sampling duration.

[0034] For example, in the application scenario of time-domain airborne electromagnetic exploration, the preset electromagnetic transmission device can be an electromagnetic transmitting coil deployed on an aircraft. That is, electromagnetic signals are transmitted to the target exploration area through the electromagnetic transmitting coil pre-deployed on the aircraft, so as to induce a corresponding electromagnetic field in the subsurface medium of the target exploration area. According to a preset sampling frequency and preset sampling duration, the induced electromagnetic signal corresponding to the induced electromagnetic field of the subsurface medium is acquired, i.e., the response electromagnetic data.

[0035] S120. Determine the filter parameters corresponding to the response electromagnetic data based on the electromagnetic signal and the preset sampling frequency.

[0036] The filter parameters can be those of a peak comb filter. Optionally, the filter parameters include at least one of the following: filter order, filter bandwidth, superimposed filter coefficients, and feedback filter coefficients. The filter order, i.e., the order of the peak comb filter, represents the degree of the highest-order term in the transfer function of the peak comb filter. The higher the filter order, the steeper the transition band (passband to stopband). That is, lower filter orders are suitable for denoising filtering requirements, while higher filter orders are suitable for scenarios such as frequency band separation. The filter bandwidth can be used to characterize the frequency range that the response electromagnetic data is allowed to pass through. In this embodiment of the invention, the filter bandwidth is the 3dB bandwidth of the peak comb filter. The superimposed filter coefficients are filter coefficients used for denoising, smoothing, and edge enhancement operations on the response electromagnetic data. The feedback filter coefficients are used to adjust the zero-point and pole-point positions of the response electromagnetic data and optimize baseband noise.

[0037] Specifically, the polarity and fundamental frequency information of the transmitted electromagnetic signal are determined. Based on the polarity and fundamental frequency information of the transmitted electromagnetic signal, and the constraint conditions or coefficient determination functions corresponding to the polarity information, at least one filter parameter among the filter order, filter bandwidth, superimposed filter coefficients, and feedback filter coefficients corresponding to the response electromagnetic data is determined, so as to perform filtering processing on the response electromagnetic data according to the filter parameters.

[0038] S130. Filter the response electromagnetic data based on the filter parameters to obtain the target electromagnetic data corresponding to the response electromagnetic data.

[0039] The target electromagnetic data can be understood as the response electromagnetic data after filtering.

[0040] Specifically, the response electromagnetic data is filtered according to the filter parameters to determine the target electromagnetic data corresponding to the response electromagnetic data. Based on the target electromagnetic data, the electrical characteristics of the underground medium in the target exploration area are determined, providing data support for subsequent exploration and development of the target exploration area.

[0041] In this embodiment of the invention, the method for filtering the response electromagnetic data to obtain the target electromagnetic data may be as follows: Based on the waveform periodic information of the response electromagnetic data, at least one time slice within a preset sampling duration and a portion of the electromagnetic data corresponding to each time slice are determined; for at least one time slice, the portion of the electromagnetic data corresponding to the target time slice is filtered according to filter parameters and a preset filtering function to obtain first output data; the next adjacent time slice of the target time slice is taken as the target time slice, and the portion of the electromagnetic data corresponding to the target time slice is filtered based on the first output data, filter parameters, and the preset filtering function to obtain second output data; the process of determining the target time slice is repeated, and the filtering process of determining the portion of the electromagnetic data corresponding to the target time slice is repeated based on the second output data, filter parameters, and the preset filtering function, until the filtering of the response electromagnetic data within the preset sampling duration is completed, and the target electromagnetic data is obtained.

[0042] Here, waveform periodicity information can be understood as the signal period of the electromagnetic wave data in response. A time slice can be a duration obtained by dividing a preset sampling time based on waveform periodicity information. Partial electromagnetic data can be understood as the response electromagnetic data corresponding to the time slice. Optionally, partial electromagnetic data within each time slice can be represented as:

[0043] X i ={x[iN],x[i-N+1],…,x[i+N-1]}

[0044] Among them, X i This represents a portion of the electromagnetic data within a time slice. x[iN] represents the response electromagnetic data corresponding to the i-th sampling point within the time slice, x[i-N+1] represents the response electromagnetic data corresponding to the (i-N+1)-th sampling point within the time slice, and x[i+N-1] represents the response electromagnetic data corresponding to the (i+N-1)-th sampling point within the time slice. N represents the filter order. The preset filtering function can be a function used to filter the response electromagnetic data. The target time slice can be the time slice corresponding to the first portion of electromagnetic data to be filtered. The first output data can be the electromagnetic data after filtering the portion of electromagnetic data from the target time slice. The second output data can be the electromagnetic data after filtering the portion of electromagnetic data from the next adjacent time slice.

[0045] For a portion of the electromagnetic data within the current time slice, filtering can be performed using the following preset filtering function.

[0046] y[j]=b0x[j]+b1x[jN]-a N y[jN], where j=i,i+1,...,i+N-1;

[0047] Where, when j = i, y[j] represents the output data corresponding to the electromagnetic data of the i-th sampling point in the current time slice after filtering. x[j] represents the electromagnetic data of the i-th sampling point in the current time slice. x[jN] represents the electromagnetic data corresponding to the iN-th sampling point. y[jN] represents the output data of the electromagnetic data corresponding to the iN-th sampling point after filtering. b0 and b1 represent the superimposed filtering coefficients. N This represents the feedback filter coefficients.

[0048] It should be noted that since the aforementioned preset filtering function can filter the electromagnetic data of N sampling points each time, if x[jN] is not present in the current time slice, it can be obtained based on some electromagnetic data from the previous time slice adjacent to the current time slice. Similarly, if y[jN] is not present in the output data of the electromagnetic data corresponding to the current time slice, it can be obtained based on the output data corresponding to some electromagnetic data from the previous time slice adjacent to the current time slice. For example, if the electromagnetic data of the Nth sampling point before the current sampling point is not present in the second time slice, the corresponding electromagnetic data can be obtained from the first time slice. Correspondingly, if the output data corresponding to the electromagnetic data of the Nth sampling point before the current sampling point is not present in the second time slice, the corresponding first output data can be obtained from the first time slice.

[0049] Specifically, based on the waveform periodic information of the response electromagnetic data, at least one time slice corresponding to a preset sampling duration is determined, along with a portion of the electromagnetic data corresponding to each time slice. Based on the at least one time slice, the starting time slice is used as the target time slice, and the portion of the electromagnetic data within the target time slice is filtered using filter parameters and a preset filtering function to determine the first output data. The next adjacent time slice is used as the target time slice, and the portion of the electromagnetic data within this target time slice is filtered using the first output data, filter parameters, and the preset filtering function to obtain the second output data. The process of determining the target time slice is repeated, based on the second output data, filter parameters, and the preset filtering function. The filtering process for determining the portion of the electromagnetic data is repeated until the filtering of the portion of the electromagnetic data in all time slices within the preset sampling duration is completed. Based on the output data corresponding to all the portion of the electromagnetic data, the target electromagnetic data is determined.

[0050] Optionally, after obtaining the target electromagnetic data corresponding to the response electromagnetic data, the method further includes: performing inversion processing based on the target electromagnetic data within a preset sampling time to determine the electrical characteristics of the underground medium in the target exploration area.

[0051] Specifically, the target electromagnetic data within a preset sampling period is divided into several "channels" according to specific rules (such as time windows or frequency intervals). The target electromagnetic data in the same channel are superimposed and averaged to reduce random noise, thereby achieving the superimposed channel extraction process of the target electromagnetic data. The superimposed channel-extracted target electromagnetic data is then inverted to determine the electrical characteristics of the subsurface medium in the target exploration area.

[0052] The technical solution of this embodiment involves transmitting electromagnetic signals to the underground medium of the target exploration area via a preset electromagnetic transmitting device. Then, based on a preset sampling frequency and a preset sampling duration, the response electromagnetic data fed back by the underground medium is acquired. Since the fed-back response electromagnetic data is susceptible to various noises, filter parameters corresponding to the response electromagnetic data can be determined based on the transmitted electromagnetic signal and the preset sampling frequency. These filter parameters are then used to filter the response electromagnetic data, yielding the target electromagnetic data corresponding to the response electromagnetic data. By determining the filter parameters based on the characteristics of the transmitted electromagnetic signal, the filtering process is simplified, achieving effective filtering of the response electromagnetic data. This amplifies the useful signal and suppresses noise signals, ensuring the accuracy of the target electromagnetic data, improving the data signal-to-noise ratio, and guaranteeing the reliability and accuracy of subsequent inversion processing based on the target electromagnetic data.

[0053] Example 2

[0054] Figure 2 This is a flowchart of a time-domain airborne electromagnetic data filtering method provided in Embodiment 2 of the present invention. This embodiment further refines the step of "determining the filter parameters corresponding to the response electromagnetic data based on the electromagnetic signal and the preset sampling frequency" based on the above embodiments. For specific implementation details, please refer to the technical solution of this embodiment. Technical terms that are the same as or corresponding to those in the above embodiments will not be repeated here. Figure 2 As shown, the method includes:

[0055] S210. After sending an electromagnetic signal to the underground medium in the target exploration area through a preset electromagnetic transmission device, the response electromagnetic data of the underground medium based on the electromagnetic signal is obtained according to the preset sampling frequency and preset sampling duration.

[0056] Among them, the response electromagnetic data is used to characterize the electrical properties of the underground medium.

[0057] S220. Determine the target fundamental frequency information and polarity information corresponding to the electromagnetic signal.

[0058] The target fundamental frequency information can be the fundamental frequency information corresponding to the transmitted electromagnetic signal. Polarity information can be used to characterize the directional characteristics of the transmitted electromagnetic signal. Optionally, the transmitted electromagnetic signal can be a unipolar wave or a bipolar wave, that is, the polarity information can be unipolar or bipolar.

[0059] Specifically, the target fundamental frequency information and polarity information of the transmitted electromagnetic signal are determined, and the filter parameters adapted to the response electromagnetic data are determined based on the target fundamental frequency information and polarity information.

[0060] In this embodiment of the invention, the target fundamental frequency information can be determined by: performing frequency domain analysis on the response electromagnetic data to determine the offset data of the frequency information of the response electromagnetic data relative to the fundamental frequency information of the electromagnetic signal; adjusting the fundamental frequency information based on the offset data to determine the target fundamental frequency information.

[0061] The fundamental frequency information of the electromagnetic signal can be the frequency information of the electromagnetic signal transmitted by a preset electromagnetic transmitting device. Offset data can be understood as the deviation between the frequency information of the responding electromagnetic data and the fundamental frequency information of the electromagnetic signal. The target fundamental frequency information can be the adjusted fundamental frequency information corresponding to the electromagnetic signal.

[0062] Specifically, since the transmission frequency of the preset electromagnetic transmitting device may drift when sending electromagnetic signals, frequency domain analysis can be performed on the received response electromagnetic data to determine the offset of the frequency information of the response electromagnetic data relative to the fundamental frequency information of the electromagnetic signal. The fundamental frequency information of the electromagnetic signal can then be adjusted based on the offset data to determine the target fundamental frequency information. Optionally, the fundamental frequency information of the electromagnetic signal can be adjusted based on the offset information to obtain the target fundamental frequency range.

[0063] Optionally, adaptively adjusting the passband width of the peak comb filter (set to ±5% of the fundamental frequency information of the electromagnetic signal) can address frequency drift (<±2%) during electromagnetic signal transmission. Based on this, a better signal-to-noise ratio can be achieved in complex electromagnetic interference scenarios, and the scenario adaptability of this embodiment of the invention is improved.

[0064] S230. Based on the preset sampling frequency, target fundamental frequency information, and polarity information, determine the filter parameters corresponding to the polarity information.

[0065] The filter parameters include at least one of the following: filter order, filter bandwidth, superimposed filter coefficients, and feedback filter coefficients.

[0066] Specifically, the filter bandwidth can be determined based on the target fundamental frequency information. The filter order can be determined based on the preset sampling frequency, target fundamental frequency information, and polarity information. The superposition filter coefficients and feedback filter coefficients can be determined based on the filter bandwidth, filter order, and preset sampling frequency. Based on this, the filter order, filter bandwidth, superposition filter coefficients, and feedback filter coefficients can be obtained, and these parameters can be used to filter the response electromagnetic data.

[0067] In this embodiment of the invention, the filter parameters can be determined as follows: the filter bandwidth is determined based on the target fundamental frequency information and preset constraints; the filter order corresponding to the polarity information is determined based on the preset sampling frequency and the target fundamental frequency information; the superposition filter coefficient and feedback filter coefficient corresponding to the polarity information are determined based on the filter order, filter bandwidth, and preset sampling frequency; and the filter bandwidth, filter order, superposition filter coefficient, and feedback filter coefficient are determined as the filter parameters for the response electromagnetic data.

[0068] The preset constraints can be pre-set constraints corresponding to the filter bandwidth. Optionally, the preset constraints can be as follows:

[0069] Δf<f0 / 2

[0070] Where Δf represents the filter bandwidth and f0 represents the target fundamental frequency information.

[0071] Specifically, based on the target fundamental frequency information and preset constraints, the 3dB bandwidth of the peak comb filter, i.e., the filter bandwidth, is determined. Based on the polarity information of the electromagnetic signal, a filter order determination function corresponding to the polarity information is determined. Based on the filter order determination function, the preset sampling frequency, and the target fundamental frequency information, the filter order is determined. Based on the polarity information of the electromagnetic signal, preset superposition filter coefficient determination functions and preset feedback filter coefficient determination functions corresponding to the polarity information are determined. Based on the preset superposition filter coefficient determination function, the filter order, the filter bandwidth, and the preset sampling frequency, the superposition filter coefficients are determined. Based on the preset feedback filter coefficient determination function, the filter order, the filter bandwidth, and the preset sampling frequency, the feedback filter coefficients are determined. The filter bandwidth, filter order, superposition filter coefficients, and feedback filter coefficients are determined as the filter parameters used for filtering the response electromagnetic data.

[0072] Optionally, the filter order can be determined based on the preset sampling frequency and the target fundamental frequency information as follows: when the polarity information is bipolar, a first ratio of the preset sampling frequency to the target fundamental frequency information of a preset multiple is determined, and the first ratio is used as the filter order corresponding to bipolarity; when the polarity information is unipolar, a second ratio of the preset sampling frequency to the target fundamental frequency information is determined, and the second ratio is used as the filter order corresponding to unipolarity.

[0073] Bipolarity characterizes the alternating positive and negative amplitude changes of the electromagnetic signal. The preset multiplier can be a pre-set value multiplied by the target fundamental frequency information. Optionally, the preset multiplier can be 2. The first ratio can be the ratio of the preset sampling frequency to the target fundamental frequency information at the preset multiplier. Unipolarity characterizes the unidirectional change of the electromagnetic signal amplitude. The second ratio can be the ratio of the preset sampling frequency to the target fundamental frequency information.

[0074] Specifically, when the polarity information of the transmitted electromagnetic signal is determined to be bipolar (i.e., the waveform of the electromagnetic signal is a bipolar wave), a first ratio of the preset sampling frequency to the target fundamental frequency information (a preset multiple) is determined, and this first ratio is used as the filter order corresponding to bipolarity. Optionally, the filter order corresponding to bipolarity can be determined using the following filter order determination function.

[0075] N = f s / (2f0)

[0076] Where N represents the filter order, f s This indicates the preset sampling frequency, and f0 represents the target base frequency information, with a preset multiple of 2.

[0077] When the polarity information of the electromagnetic signal is determined to be unipolar, that is, when the waveform of the electromagnetic signal is a unipolar wave, a second ratio of the preset sampling frequency to the target fundamental frequency information is determined, and this second ratio is used as the filter order corresponding to unipolarity.

[0078] Alternatively, the filter order corresponding to unipolarity can be determined using the following filter order determination function.

[0079] N = f s / f0

[0080] Where N represents the filter order, f s f0 represents the preset sampling frequency and f0 represents the target base frequency information.

[0081] Optionally, the superposition filter coefficients and feedback filter coefficients can be determined based on the filter order, filter bandwidth, and preset sampling frequency as follows: When the polarity information is bipolar, a first parameter is determined based on the filter order, filter bandwidth, preset sampling frequency, and preset tangent function corresponding to bipolarity; the superposition filter coefficients and feedback filter coefficients corresponding to bipolarity are determined based on the first parameter and the preset superposition filter coefficient determination function and preset feedback filter coefficient determination function corresponding to bipolarity, respectively; when the polarity information is unipolar, a second parameter is determined based on the filter order, filter bandwidth, preset sampling frequency, and preset tangent function corresponding to unipolarity; the superposition filter coefficients and feedback filter coefficients corresponding to unipolarity are determined based on the second parameter and the preset superposition filter coefficient determination function and preset feedback filter coefficient determination function corresponding to unipolarity, respectively.

[0082] The preset tangent function can be a pre-defined tangent function used to process the filter order, filter bandwidth, and preset sampling frequency. Optionally, the preset tangent function can be expressed as:

[0083] A = tan(2NπΔf / fs)

[0084] Where N represents the filter order, Δf represents the filter bandwidth, and fs represents the preset sampling frequency. It should be noted that when N represents the filter order corresponding to bipolarity, A represents the first parameter. When N represents the filter order corresponding to unipolarity, A represents the second parameter. That is, the first parameter can be the result of substituting the filter order, filter bandwidth, and preset sampling frequency corresponding to bipolarity into the preset tangent function. The second parameter can be the result of substituting the filter order, filter bandwidth, and preset sampling frequency corresponding to unipolarity into the preset tangent function.

[0085] A preset superposition filter coefficient determination function corresponding to bipolarity can be used to determine the superposition filter coefficients corresponding to bipolarity. Based on this, when the polarity information of the transmitted electromagnetic signal is bipolar, the response electromagnetic data corresponding to the electromagnetic signal can be filtered using the superposition filter coefficients corresponding to bipolarity. Optionally, the preset superposition filter coefficient determination function corresponding to bipolarity can be expressed as:

[0086]

[0087] Where b0 and b1 represent the superimposed filter coefficients corresponding to bipolarity, N represents the filter order, Δf represents the filter bandwidth, and fs represents the preset sampling frequency.

[0088] A preset superposition filter coefficient determination function corresponding to unipolarity can be used to determine the superposition filter coefficients corresponding to unipolarity. Based on this, when the polarity information of the transmitted electromagnetic signal is unipolar, the response electromagnetic data corresponding to the electromagnetic signal can be filtered using the superposition filter coefficients corresponding to unipolarity. Optionally, the preset superposition filter coefficient determination function corresponding to unipolarity can be expressed as:

[0089]

[0090] Where b0' and b1' represent the superimposed filter coefficients corresponding to unipolarity, N represents the filter order, Δf represents the filter bandwidth, and fs represents the preset sampling frequency.

[0091] The feedback filter coefficients corresponding to bipolarity can be used to determine the feedback filter coefficients corresponding to bipolarity. Based on this, when the polarity information of the transmitted electromagnetic signal is bipolar, the response electromagnetic data corresponding to the electromagnetic signal can be filtered using the feedback filter coefficients corresponding to bipolarity. Optionally, the function for determining the preset feedback filter coefficients corresponding to bipolarity can be expressed as:

[0092]

[0093] Where N represents the filter order, Δf represents the filter bandwidth, fs represents the preset sampling frequency, and a N This represents the feedback filter coefficients corresponding to bipolarity.

[0094] The feedback filter coefficients corresponding to unipolarity can be used to determine the feedback filter coefficients corresponding to unipolarity. Based on this, when the polarity information of the transmitted electromagnetic signal is unipolar, the response electromagnetic data corresponding to the electromagnetic signal can be filtered using the feedback filter coefficients corresponding to unipolarity. Optionally, the function for determining the preset feedback filter coefficients corresponding to unipolarity can be expressed as:

[0095]

[0096] Among them, a N ' represents the feedback filter coefficient corresponding to unipolarity, N represents the filter order, Δf represents the filter bandwidth, and fs represents the preset sampling frequency.

[0097] It should be noted that when determining the superposition filter coefficients and feedback filter coefficients, the polarity information of the electromagnetic signal needs to be determined so that the filter order, filter bandwidth, and preset sampling frequency corresponding to the polarity information can be substituted into the preset superposition filter coefficient determination function and preset feedback filter coefficient determination function corresponding to the polarity information.

[0098] Specifically, when the polarity information is bipolar, the filter order, filter bandwidth, and preset sampling frequency corresponding to bipolarity are substituted into a preset tangent function to determine the first parameter corresponding to bipolarity. Based on the first parameter and a preset superposition filter coefficient determination function corresponding to bipolarity, the superposition filter coefficients corresponding to bipolarity are determined. Based on the first parameter and a preset feedback filter coefficient determination function corresponding to bipolarity, the feedback filter coefficients corresponding to bipolarity are determined.

[0099] When the polarity information is unipolar, the filter order, filter bandwidth, and preset sampling frequency corresponding to the unipolarity are substituted into a preset tangent function to determine the second parameter corresponding to the unipolarity. Based on the second parameter and a preset function for determining superimposed filter coefficients corresponding to the unipolarity, the superimposed filter coefficients corresponding to the unipolarity are determined. Based on the second parameter and a preset function for determining feedback filter coefficients corresponding to the unipolarity, the feedback filter coefficients corresponding to the unipolarity are determined.

[0100] S240. Filter the response electromagnetic data based on the filter parameters to obtain the target electromagnetic data corresponding to the response electromagnetic data.

[0101] The technical solution of this embodiment involves sending electromagnetic signals to the underground medium of the target exploration area via a preset electromagnetic transmitting device. Then, based on a preset sampling frequency and a preset sampling duration, the response electromagnetic data fed back by the underground medium is acquired. The target fundamental frequency information and polarity information corresponding to the electromagnetic signal are determined. Based on the preset sampling frequency, target fundamental frequency information, and polarity information, filter parameters corresponding to the polarity information are determined. These filter parameters are peak comb filter parameters. Filtering the response electromagnetic data using the peak comb filter parameters to obtain the target electromagnetic data avoids waveform distortion caused by time-domain smoothing. While improving the signal-to-noise ratio, it ensures the signal integrity of the target electromagnetic data, such as preserving key time-domain features like the signal rise edge and attenuation curve. Determining the filter parameters based on the characteristics of the transmitted electromagnetic signal allows for targeted enhancement of the effective signal energy in the response electromagnetic data and concentration of the signal energy within specific harmonic channels, thereby increasing the signal amplitude of the target electromagnetic data. Meanwhile, based on the multi-passband characteristics of the comb filter, a passband can be formed at integer multiples of the fundamental frequency, while a deep stopband (typically reaching -40dB or more) can be established in the non-harmonic frequency band. Based on this, various heterogeneous noise sources in the response electromagnetic data, such as power frequency interference (50 / 60Hz), high-frequency random noise (>1kHz), and motion noise (0.1-10Hz), can be filtered using the filter parameters determined in this embodiment of the invention. This simplifies the filtering process, achieves effective filtering of the response electromagnetic data, amplifies the useful signal, suppresses noise signals, ensures the accuracy of the target electromagnetic data, improves the data signal-to-noise ratio, and guarantees the reliability and accuracy of subsequent inversion processing based on the target electromagnetic data.

[0102] Example 3

[0103] Figure 3 This is a schematic diagram of a time-domain airborne electromagnetic data filtering device provided in Embodiment 3 of the present invention. Figure 3 As shown, the device includes: an electromagnetic data acquisition module 310, a filter parameter determination module 320, and an electromagnetic data filtering module 330.

[0104] The electromagnetic data acquisition module 310 is used to acquire the response electromagnetic data of the underground medium based on the electromagnetic signal after sending an electromagnetic signal to the underground medium in the target exploration area through a preset electromagnetic transmitting device, according to a preset sampling frequency and a preset sampling duration. The response electromagnetic data is used to characterize the electrical characteristics of the underground medium. The filter parameter determination module 320 is used to determine the filter parameters corresponding to the response electromagnetic data according to the electromagnetic signal and the preset sampling frequency. The electromagnetic data filtering module 330 is used to filter the response electromagnetic data based on the filter parameters to obtain the target electromagnetic data corresponding to the response electromagnetic data.

[0105] The technical solution of this embodiment involves transmitting electromagnetic signals to the underground medium of the target exploration area via a preset electromagnetic transmitting device. Then, based on a preset sampling frequency and a preset sampling duration, the response electromagnetic data fed back by the underground medium is acquired. Since the fed-back response electromagnetic data is susceptible to various noises, filter parameters corresponding to the response electromagnetic data can be determined based on the transmitted electromagnetic signal and the preset sampling frequency. These filter parameters are then used to filter the response electromagnetic data, yielding the target electromagnetic data corresponding to the response electromagnetic data. By determining the filter parameters based on the characteristics of the transmitted electromagnetic signal, the filtering process is simplified, achieving effective filtering of the response electromagnetic data. This amplifies the useful signal and suppresses noise signals, ensuring the accuracy of the target electromagnetic data, improving the data signal-to-noise ratio, and guaranteeing the reliability and accuracy of subsequent inversion processing based on the target electromagnetic data.

[0106] Based on the above embodiments, optionally, the filter parameter determination module includes: a fundamental frequency information and polarity information determination unit, used to determine the target fundamental frequency information and polarity information corresponding to the electromagnetic signal; and a filter parameter determination unit, used to determine the filter parameters corresponding to the polarity information according to the preset sampling frequency, the target fundamental frequency information and the polarity information, wherein the filter parameters include at least one of the following: filter order, filter bandwidth, superimposed filter coefficients and feedback filter coefficients.

[0107] Optionally, the fundamental frequency information and polarity information determination unit includes: a fundamental frequency information determination subunit, used to perform frequency domain analysis on the response electromagnetic data, determine the offset data of the frequency information of the response electromagnetic data relative to the fundamental frequency information of the electromagnetic signal; and adjust the fundamental frequency information based on the offset data to determine the target fundamental frequency information.

[0108] Optionally, the filter parameter determination unit includes: a filter bandwidth determination subunit, used to determine the filter bandwidth based on the target fundamental frequency information and preset constraints; a filter order determination subunit, used to determine the filter order corresponding to the polarity information based on the preset sampling frequency and the target fundamental frequency information; a filter coefficient determination subunit, used to determine the superimposed filter coefficient and feedback filter coefficient corresponding to the polarity information based on the filter order, filter bandwidth, and preset sampling frequency; and the filter parameter determination subunit is used to determine the filter bandwidth, filter order, superimposed filter coefficient, and feedback filter coefficient as filter parameters for the response electromagnetic data.

[0109] Optionally, the filter order determination subunit is used to determine a first ratio of the preset sampling frequency to the target fundamental frequency information by a preset multiple when the polarity information is bipolar, and to use the first ratio as the filter order corresponding to bipolarity; and to determine a second ratio of the preset sampling frequency to the target fundamental frequency information when the polarity information is unipolar, and to use the second ratio as the filter order corresponding to unipolarity.

[0110] Optionally, the filter coefficient determination subunit is used to determine the first parameter based on the filter order, filter bandwidth, preset sampling frequency, and preset tangent function corresponding to bipolarity when the polarity information is bipolar; and to determine the superimposed filter coefficient and feedback filter coefficient corresponding to bipolarity based on the first parameter and the preset superimposed filter coefficient determination function and preset feedback filter coefficient determination function corresponding to bipolarity, respectively; and to determine the second parameter based on the filter order, filter bandwidth, preset sampling frequency, and preset tangent function corresponding to unipolarity when the polarity information is unipolarity; and to determine the superimposed filter coefficient and feedback filter coefficient corresponding to unipolarity based on the second parameter and the preset superimposed filter coefficient determination function and preset feedback filter coefficient determination function corresponding to unipolarity, respectively.

[0111] Optionally, an electromagnetic data filtering module is used to determine at least one time slice and a portion of the electromagnetic data corresponding to each time slice based on the waveform period information of the response electromagnetic data; for at least one time slice, the portion of the electromagnetic data corresponding to the target time slice is filtered according to filter parameters and a preset filtering function to obtain first output data; the next adjacent time slice of the target time slice is taken as the target time slice, and the portion of the electromagnetic data corresponding to the target time slice is filtered based on the first output data, filter parameters, and the preset filtering function to obtain second output data; the process of determining the target time slice is repeated, and the filtering process of determining the portion of the electromagnetic data corresponding to the target time slice is repeated based on the second output data, filter parameters, and the preset filtering function, until the filtering process of the response electromagnetic data within the preset sampling time is completed to obtain the target electromagnetic data.

[0112] Optionally, the device further includes an electrical characteristic determination module, used to perform inversion processing based on the target electromagnetic data within a preset sampling period to determine the electrical characteristics of the underground medium in the target exploration area.

[0113] The time-domain airborne electromagnetic data filtering device provided in the embodiments of the present invention can execute the time-domain airborne electromagnetic data filtering method provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of executing the method.

[0114] Example 4

[0115] Figure 4 This is a schematic diagram of the structure of an electronic device provided in Embodiment 4 of the present invention. The electronic device 10 is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device may also represent various forms of mobile devices, such as personal digital processors, cellular phones, smartphones, wearable devices (such as helmets, glasses, watches, etc.), and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the invention described and / or claimed herein.

[0116] like Figure 4 As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12 or a random access memory (RAM) 13, communicatively connected to the at least one processor 11. The memory stores computer programs executable by the at least one processor. The processor 11 can perform various appropriate actions and processes based on the computer program stored in the ROM 12 or loaded from storage unit 18 into the RAM 13. The RAM 13 may also store various programs and data required for the operation of the electronic device 10. The processor 11, ROM 12, and RAM 13 are interconnected via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0117] Multiple components in electronic device 10 are connected to I / O interface 15, including: input unit 16, such as keyboard, mouse, etc.; output unit 17, such as various types of displays, speakers, etc.; storage unit 18, such as disk, optical disk, etc.; and communication unit 19, such as network card, modem, wireless transceiver, etc. Communication unit 19 allows electronic device 10 to exchange information / data with other devices through computer networks such as the Internet and / or various telecommunications networks.

[0118] Processor 11 can be a variety of general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special-purpose artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. Processor 11 performs the various methods and processes described above, such as time-domain airborne electromagnetic data filtering methods.

[0119] In some embodiments, the time-domain airborne electromagnetic data filtering method may be implemented as a computer program tangibly contained in a computer-readable storage medium, such as storage unit 18. In some embodiments, part or all of the computer program may be loaded and / or installed on electronic device 10 via ROM 12 and / or communication unit 19. When the computer program is loaded into RAM 13 and executed by processor 11, one or more steps of the time-domain airborne electromagnetic data filtering method described above may be performed. Alternatively, in other embodiments, processor 11 may be configured to perform the time-domain airborne electromagnetic data filtering method by any other suitable means (e.g., by means of firmware).

[0120] Various embodiments of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field-programmable gate arrays (FPGAs), application-specific integrated circuits (ASICs), application-specific standard products (ASSPs), system-on-a-chip (SoCs), complex programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments may include implementations in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which may be a dedicated or general-purpose programmable processor, capable of receiving data and instructions from a storage system, at least one input device, and at least one output device, and transmitting data and instructions to the storage system, the at least one input device, and the at least one output device.

[0121] Computer programs for implementing the time-domain airborne electromagnetic data filtering method of the present invention can be written in any combination of one or more programming languages. These computer programs can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when executed by the processor, the computer programs cause the functions / operations specified in the flowcharts and / or block diagrams to be implemented. The computer programs can be executed entirely on a machine, partially on a machine, as a standalone software package partially on a machine and partially on a remote machine, or entirely on a remote machine or server.

[0122] In particular, according to embodiments of the present invention, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, embodiments of the present invention include a computer program product comprising a computer program carried on a non-transitory computer-readable medium, the computer program containing program code for performing the methods shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via communication unit 19, or installed from storage unit 18, or installed from ROM 12. When the computer program is executed by processor 11, it performs the functions defined in the methods of the embodiments of the present invention.

[0123] Example 5

[0124] Embodiment 5 of the present invention also provides a computer-readable storage medium storing computer instructions for causing a processor to execute a time-domain airborne electromagnetic data filtering method, the method comprising:

[0125] After transmitting electromagnetic signals to the underground medium in the target exploration area through a preset electromagnetic transmitting device, the response electromagnetic data of the underground medium based on the electromagnetic signals is obtained according to the preset sampling frequency and preset sampling duration. The response electromagnetic data is used to characterize the electrical characteristics of the underground medium. Based on the electromagnetic signals and the preset sampling frequency, the filter parameters corresponding to the response electromagnetic data are determined. The response electromagnetic data is filtered based on the filter parameters to obtain the target electromagnetic data corresponding to the response electromagnetic data.

[0126] In the context of this invention, a computer-readable storage medium can be a tangible medium that may contain or store a computer program for use by or in conjunction with an instruction execution system, apparatus, or device. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination thereof. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. More specific examples of machine-readable storage media include electrical connections based on one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fibers, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0127] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user provides input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including sound input, voice input, or tactile input).

[0128] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as data servers), or computing systems that include middleware components (e.g., application servers), or computing systems that include frontend components (e.g., user computers with graphical user interfaces or web browsers through which users can interact with implementations of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., communication networks). Examples of communication networks include local area networks (LANs), wide area networks (WANs), blockchain networks, and the Internet.

[0129] A computing system can include clients and servers. Clients and servers are generally located far apart and typically interact through communication networks. The client-server relationship is created by computer programs running on the respective computers and having a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or cloud host, which is a hosting product within the cloud computing service system to address the shortcomings of traditional physical hosts and VPS services, such as high management difficulty and weak business scalability.

[0130] It should be understood that the various forms of processes shown above can be used, with steps reordered, added, or deleted. For example, the steps described in this invention can be executed in parallel, sequentially, or in different orders, as long as the desired result of the technical solution of this invention can be achieved, and this is not limited herein.

[0131] The specific embodiments described above do not constitute a limitation on the scope of protection of this invention. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of this invention should be included within the scope of protection of this invention.

Claims

1. A time-domain airborne electromagnetic data filtering method, characterized in that, include: After sending an electromagnetic signal to the underground medium in the target exploration area through a preset electromagnetic transmission device, the response electromagnetic data of the underground medium based on the electromagnetic signal is obtained according to a preset sampling frequency and a preset sampling duration. The response electromagnetic data is used to characterize the electrical characteristics of the underground medium. Based on the electromagnetic signal and the preset sampling frequency, determine the filter parameters corresponding to the response electromagnetic data; The response electromagnetic data is filtered based on the filter parameters to obtain the target electromagnetic data corresponding to the response electromagnetic data. The step of determining the filter parameters corresponding to the response electromagnetic data based on the electromagnetic signal and the preset sampling frequency includes: The target fundamental frequency information and polarity information corresponding to the electromagnetic signal are determined; based on the preset sampling frequency, the target fundamental frequency information and the polarity information, filter parameters corresponding to the polarity information are determined, wherein the filter parameters include at least one of the following: filter order, filter bandwidth, superposition filter coefficients and feedback filter coefficients; The step of determining the filter parameters corresponding to the polarity information based on the preset sampling frequency, the target fundamental frequency information, and the polarity information includes: Based on the target fundamental frequency information and preset constraints, the filter bandwidth is determined; based on the preset sampling frequency and the target fundamental frequency information, the filter order corresponding to the polarity information is determined; based on the filter order, the filter bandwidth, and the preset sampling frequency, the superposition filter coefficient and feedback filter coefficient corresponding to the polarity information are determined; the filter bandwidth, the filter order, the superposition filter coefficient, and the feedback filter coefficient are determined as the filter parameters of the response electromagnetic data. The step of determining the superposition filter coefficients and feedback filter coefficients corresponding to the polarity information based on the filter order, the filter bandwidth, and the preset sampling frequency includes: When the polarity information is bipolar, a first parameter is determined based on the filter order corresponding to the bipolarity, the filter bandwidth, the preset sampling frequency, and the preset tangent function; based on the first parameter and the preset superposition filter coefficient determination function and the preset feedback filter coefficient determination function corresponding to the bipolarity, the superposition filter coefficient and the feedback filter coefficient corresponding to the bipolarity are determined respectively. When the polarity information is unipolar, a second parameter is determined based on the filter order corresponding to the unipolarity, the filter bandwidth, the preset sampling frequency, and the preset tangent function; based on the second parameter and the preset superposition filter coefficient determination function and the preset feedback filter coefficient determination function corresponding to the unipolarity, the superposition filter coefficient and the feedback filter coefficient corresponding to the unipolarity are determined respectively.

2. The method according to claim 1, characterized in that, Determining the target fundamental frequency information corresponding to the electromagnetic signal includes: Frequency domain analysis is performed on the response electromagnetic data to determine the offset data of the frequency information of the response electromagnetic data relative to the fundamental frequency information of the electromagnetic signal; The baseband information is adjusted based on the offset data to determine the target baseband information.

3. The method according to claim 1, characterized in that, The step of determining the filter order corresponding to the polarity information based on the preset sampling frequency and the target fundamental frequency information includes: When the polarity information is bipolar, a first ratio of the preset sampling frequency to the target fundamental frequency information of a preset multiple is determined, and the first ratio is used as the filter order corresponding to the bipolarity. When the polarity information is unipolar, a second ratio of the preset sampling frequency to the target fundamental frequency information is determined, and the second ratio is used as the filter order corresponding to the unipolarity.

4. The method according to claim 1, characterized in that, The step of filtering the response electromagnetic data based on the filter parameters to obtain target electromagnetic data corresponding to the response electromagnetic data includes: Based on the waveform period information of the response electromagnetic data, at least one time slice within the preset sampling duration and a portion of the electromagnetic data corresponding to each time slice are determined. For the at least one time slice, the electromagnetic data corresponding to the target time slice is filtered according to the filter parameters and the preset filtering function to obtain the first output data; The next adjacent time slice of the target time slice is taken as the target time slice, and based on the first output data, the filter parameters and the preset filtering function, the electromagnetic data corresponding to the target time slice is filtered to obtain the second output data. The process of determining the target time slice is repeated, based on the second output data, the filter parameters, and the preset filter function, to repeatedly determine the filtering process of the partial electromagnetic data corresponding to the target time slice, until the filtering process of the response electromagnetic data within the preset sampling time is completed, and the target electromagnetic data is obtained.

5. The method according to claim 1, characterized in that, After obtaining the target electromagnetic data corresponding to the response electromagnetic data, the method further includes: Based on the target electromagnetic data within the preset sampling time, inversion processing is performed to determine the electrical characteristics of the underground medium in the target exploration area.

6. A time-domain airborne electromagnetic data filtering device, characterized in that, include: An electromagnetic data acquisition module is used to acquire response electromagnetic data fed back by the underground medium based on the electromagnetic signal after sending an electromagnetic signal to the underground medium in the target exploration area through a preset electromagnetic transmission device, according to a preset sampling frequency and a preset sampling duration. The response electromagnetic data is used to characterize the electrical characteristics of the underground medium. The filter parameter determination module is used to determine the filter parameters corresponding to the response electromagnetic data based on the electromagnetic signal and the preset sampling frequency. An electromagnetic data filtering module is used to filter the response electromagnetic data based on the filter parameters to obtain target electromagnetic data corresponding to the response electromagnetic data. The filter parameter determination module includes: a fundamental frequency information and polarity information determination unit, used to determine the target fundamental frequency information and polarity information corresponding to the electromagnetic signal; and a filter parameter determination unit, used to determine filter parameters corresponding to the polarity information based on the preset sampling frequency, the target fundamental frequency information, and the polarity information, wherein the filter parameters include at least one of: filter order, filter bandwidth, superimposed filter coefficients, and feedback filter coefficients. The filter parameter determination unit includes: a filter bandwidth determination subunit, used to determine the filter bandwidth based on the target fundamental frequency information and preset constraints; a filter order determination subunit, used to determine the filter order corresponding to the polarity information based on the preset sampling frequency and the target fundamental frequency information; a filter coefficient determination subunit, used to determine the superimposed filter coefficient and feedback filter coefficient corresponding to the polarity information based on the filter order, the filter bandwidth, and the preset sampling frequency; and a filter parameter determination subunit, used to determine the filter bandwidth, the filter order, the superimposed filter coefficient, and the feedback filter coefficient as the filter parameters of the response electromagnetic data. The filter coefficient determination subunit is configured to: determine a first parameter based on the filter order corresponding to the bipolarity, the filter bandwidth, the preset sampling frequency, and the preset tangent function when the polarity information is bipolar; determine the superimposed filter coefficient and the feedback filter coefficient corresponding to the bipolarity based on the first parameter and the preset superimposed filter coefficient determination function and the preset feedback filter coefficient determination function corresponding to the bipolarity, respectively; determine a second parameter based on the filter order corresponding to the unipolarity, the filter bandwidth, the preset sampling frequency, and the preset tangent function when the polarity information is unipolar; and determine the superimposed filter coefficient and the feedback filter coefficient corresponding to the unipolarity based on the second parameter and the preset superimposed filter coefficient determination function and the preset feedback filter coefficient determination function corresponding to the unipolarity, respectively.

7. An electronic device, characterized in that, The electronic device includes: At least one processor; and A memory communicatively connected to the at least one processor; wherein, The memory stores a computer program that can be executed by the at least one processor, the computer program being executed by the at least one processor to enable the at least one processor to perform the time-domain airborne electromagnetic data filtering method according to any one of claims 1-5.

Images