A method, device and electronic device for extracting signal parameters under noise interference

Abstract

The present application provides a signal parameter extraction method, device and electronic device under noise interference, which are used to perform correlation analysis based on signal superposition and cancellation to achieve accurate extraction of target signals. The method includes: obtaining an NMR signal under NMR processing, where the NMR signal is formed by superposition of a target signal and a noise signal; constructing a second formula when the target signal is represented by a first formula, where the first formula is: #imgabs0#E(t) s represents the target signal, E s is the initial amplitude of the signal, T s is the average decay time, f t is the Larmor frequency, #imgabs1# is the initial phase, and the second formula is: #imgabs2# Z1(t) and Z2(t) are E(t) s Two signals with correlation are constructed, S(t) represents the NMR signal, a and b are real numbers; a and b are taken, and the cross-correlation function of Z1(t) and Z2(t) is calculated by the third formula: #imgabs3#R(τ) represents the cross-correlation function; the cross-correlation function is E s , T s , solve #imgabs4#, and bring the solution into the first formula to restore the target signal.

Description

Technical Field

[0001] The present application relates to the field of signal detection technology, and in particular to a method, device and electronic equipment for extracting signal parameters under noise interference. Background Art

[0002] There are many kinds of useful signals in human society, such as communication signals, radar signals, non-destructive testing signals, etc. These signals often carry information that people want to know. For example, in communication signals, the signal corresponding to human voice carries sound information, in radar signals, the signal reflected from the target's position carries position information, and in non-destructive testing signals, there is abnormal position information of the target body.

[0003] However, the signal transmission process currently faces various types of noise interference. The type of interference is related to the corresponding signal receiving method. For example, the radar receives electromagnetic wave signals, so it is inevitably subject to electromagnetic noise interference. The seismic exploration method in non-destructive testing receives signals from the vibration of the target body, so it is inevitably subject to vibration signal interference caused by surrounding pedestrians and vehicles. The noise signal often covers up the target signal we want to obtain, making it impossible to accurately obtain the corresponding signal parameters, and it is also impossible to accurately restore the information carried in the signal.

[0004] In the face of such problems, the commonly used methods are often to deal with noise. For example, a signal can be repeatedly collected and then superimposed to reduce the interference of noise in the signal. For example, the frequency difference between the noise signal and the target signal can be used to design a filter for denoising.

[0005] In the process of studying the existing related technologies, the inventors found that the existing denoising technologies have certain limitations and cannot completely eliminate the noise signal, which means that it is still impossible to obtain an accurate target signal, but only a signal with a relatively low noise content, which is still not conducive to accurately obtaining the information carried by the target signal. Summary of the invention

[0006] The present application provides a method, device and electronic device for extracting signal parameters under noise interference, which are used to perform correlation analysis based on signal superposition and cancellation to achieve accurate extraction of target signals.

[0007] In a first aspect, the present application provides a method for extracting signal parameters under noise interference, the method comprising:

[0008] Acquiring a nuclear magnetic resonance (NMR) signal collected under NMR detection processing, wherein the NMR signal is formed by superposition of a target signal and a noise signal;

[0009] When the target signal is expressed by the first formula, the second formula is constructed. The first formula is: Among them, E(t) s represents the target signal, E s is the initial amplitude of the signal, T s is the average decay time, f t is the Larmor frequency, is the initial phase, and the second formula is: Where Z1(t) and Z2(t) are E(t) s Two signals with correlation are constructed, S(t) represents NMR signal, a and b are real numbers;

[0010] Take values ​​for a and b and bring them into the third formula to calculate the cross-correlation function of Z1(t) and Z2(t). The third formula is: Where R(τ) represents the cross-correlation function;

[0011] E s , T s , , and bring the solution into the first formula to restore the target signal.

[0012] In combination with the first aspect of the present application, in a first possible implementation of the first aspect of the present application, the cross-correlation function is E s , T s , The solution includes:

[0013] The cross-correlation function is E- s , T s , The preset mathematical operations include at least one of addition, subtraction, multiplication, division, absolute value, extreme value, norm, derivative and partial derivative.

[0014] In combination with the first aspect of the present application, in a second possible implementation of the first aspect of the present application, the values ​​of a and b are specifically two groups, and after the two groups of a and b are substituted into the third formula respectively, two cross-correlation functions R are obtained. 12 (τ) and R 34 (τ), the cross-correlation function is E s , T s , The solution includes:

[0015] R12 (τ) and R 34 (τ), and find the 2-norm respectively, and get ||R 12 (τ)||2 and||R 34 (τ)||2;

[0016] Ask||R 12 (τ)||2 and||R 34 (τ)||2, we get

[0017] The average decay time range of 0 to 1 second, the phase range of 0 to 2π, and the initial amplitude range extracted from the NMR signal are brought into G to solve the problem of E when G reaches its maximum value. s , T s , The value of .

[0018] In combination with the second possible implementation of the first aspect of the present application, in a third possible implementation of the first aspect of the present application, in the process of solving G, an improved artificial fish swarm algorithm, simulated annealing method or genetic algorithm is used for search processing.

[0019] In combination with the first aspect of the present application, in a fourth possible implementation of the first aspect of the present application, the NMR signal is specifically a signal collected by NMR detection during groundwater detection.

[0020] In combination with the fourth possible implementation of the first aspect of the present application, in a fifth possible implementation of the first aspect of the present application, before obtaining and processing the collected NMR signal under magnetic resonance NMR detection, the method further includes:

[0021] Initiate NMR detection processing in the target area.

[0022] In combination with the fourth possible implementation of the first aspect of the present application, in a sixth possible implementation of the first aspect of the present application, Es, Ts, After solving the problem and substituting the solution result into the first formula to restore the target signal, the method further includes:

[0023] Based on the changing patterns reflected by different target signals restored through different NMR detection processes, the groundwater location, groundwater burial depth and groundwater content are determined.

[0024] In a second aspect, the present application provides a signal parameter extraction device under noise interference, the device comprising:

[0025] An acquisition unit, used for acquiring an NMR signal collected under NMR detection processing, where the NMR signal is formed by superposition of a target signal and a noise signal;

[0026] The construction unit is used to construct a second formula when the target signal is expressed by the first formula, and the first formula is: Among them, E(t) s represents the target signal, E s is the initial amplitude of the signal, T s is the average decay time, f t is the preset Larmor frequency, is the initial phase, and the second formula is: Wherein, Z1(t) and Z2(t) are two signals with correlation constructed by E(t)s, S(t) represents NMR signal, and a and b are real numbers;

[0027] Substitute into the unit for taking values ​​of a and b, and substitute into the third formula to calculate the cross-correlation function of Z1(t) and Z2(t). The third formula is: Where R(τ) represents the cross-correlation function;

[0028] Solving unit, used to perform Es, Ts, , and bring the solution into the first formula to restore the target signal.

[0029] In conjunction with the second aspect of the present application, in a first possible implementation manner of the second aspect of the present application, the solving unit is specifically configured to:

[0030] The cross-correlation function is E- s , T s , The preset mathematical operations include at least one of addition, subtraction, multiplication, division, absolute value, extreme value, norm, derivative and partial derivative.

[0031] In combination with the second aspect of the present application, in a second possible implementation of the second aspect of the present application, the values ​​of a and b are specifically two groups, and after the two groups of a and b are substituted into the third formula respectively, two cross-correlation functions R are obtained. 12 (τ) and R 34 (τ), solving unit, specifically used for:

[0032] R 12 (τ) and R 34 (τ), and find the 2-norm respectively, and get ||R 12 (τ)||2 and||R 34 (τ)||2;

[0033] Ask||R 12 (τ)||2 and||R 34 (τ)||2, we get

[0034] The average decay time range of 0 to 1 second, the phase range of 0 to 2π, and the initial amplitude range extracted from the NMR signal are brought into G to solve the problem of E when G reaches its maximum value. s , T s , The value of .

[0035] In combination with the second possible implementation of the second aspect of the present application, in a third possible implementation of the second aspect of the present application, in the process of solving G, an improved artificial fish swarm algorithm, simulated annealing method or genetic algorithm is used for search processing.

[0036] In combination with the second aspect of the present application, in a fourth possible implementation of the second aspect of the present application, the NMR signal is specifically a signal collected by NMR detection during groundwater detection.

[0037] In combination with the fourth possible implementation manner of the second aspect of the present application, in a fifth possible implementation manner of the second aspect of the present application, the apparatus further includes an initiating unit, configured to:

[0038] Initiate NMR detection processing in the target area.

[0039] In combination with the fourth possible implementation manner of the second aspect of the present application, in a sixth possible implementation manner of the second aspect of the present application, the apparatus further includes a determining unit, which is configured to:

[0040] Based on the changing patterns reflected by different target signals restored through different NMR detection processes, the groundwater location, groundwater burial depth and groundwater content are determined.

[0041] In a third aspect, the present application provides an electronic device comprising a processor and a memory, wherein a computer program is stored in the memory, and when the processor calls the computer program in the memory, the method provided in the first aspect of the present application or any possible implementation method of the first aspect of the present application is executed.

[0042] In a fourth aspect, the present application provides a computer-readable storage medium, which stores multiple instructions, and the instructions are suitable for a processor to load to execute the method provided by the first aspect of the present application or any possible implementation of the first aspect of the present application.

[0043] From the above content, it can be concluded that the present application has the following beneficial effects:

[0044] Regarding the extraction of the real target signal in the NMR signal, compared with the prior art that estimates and eliminates the noise signal of the NMR signal, the present application bypasses the noise signal, overcomes the problem that the noise signal is difficult to be approximately completely eliminated, and directly targets the characteristics of the target signal itself, performs correlation analysis on the basis of signal superposition and cancellation, to achieve accurate extraction of the target signal, which can promote a higher-precision NMR detection effect. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For those skilled in the art, other drawings can be obtained based on these drawings without creative work.

[0046] Figure 1 A schematic diagram of a flow chart of a signal parameter extraction method under noise interference of the present application;

[0047] Figure 2 A signal schematic diagram of the NMR signal of the present application;

[0048] Figure 3 This is another signal schematic diagram of the NMR signal of the present application;

[0049] Figure 4 A schematic diagram of a scenario of a signal parameter extraction method under noise interference in this application;

[0050] Figure 5 A schematic diagram of the structure of a signal parameter extraction device under noise interference in this application;

[0051] Figure 6 This is a structural schematic diagram of the electronic device of this application. DETAILED DESCRIPTION

[0052] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of this application.

[0053] The terms "first", "second", etc. in the specification and claims of this application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments described here can be implemented in an order other than that illustrated or described here. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or modules is not necessarily limited to those steps or modules clearly listed, but may include other steps or modules that are not clearly listed or inherent to these processes, methods, products or devices. The naming or numbering of steps in this application does not mean that the steps in the method flow must be executed in the time / logical sequence indicated by the naming or numbering. The process steps that have been named or numbered can change the execution order according to the technical purpose to be achieved, as long as the same or similar technical effects can be achieved.

[0054] The division of modules in this application is a logical division. There may be other division methods when it is implemented in actual applications. For example, multiple modules can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, and the indirect coupling or communication connection between modules can be electrical or other similar forms, which are not limited in this application. In addition, the modules or submodules described as separate components may or may not be physically separated, may or may not be physical modules, or may be distributed in multiple circuit modules, and some or all of the modules may be selected according to actual needs to achieve the purpose of the present application.

[0055] Before introducing the signal parameter extraction method under noise interference provided by the present application, the background content involved in the present application is first introduced.

[0056] The signal parameter extraction method, device and computer-readable storage medium under noise interference provided in the present application can be applied to electronic devices to perform correlation analysis based on signal superposition and cancellation to achieve accurate extraction of target signals.

[0057] The signal parameter extraction method under noise interference mentioned in the present application may be executed by a signal parameter extraction device under noise interference, or a server, a physical host, or a user equipment (UE) and other different types of electronic devices integrated with the signal parameter extraction device under noise interference. Among them, the signal parameter extraction device under noise interference may be implemented in hardware or software, and the UE may be a terminal device such as a smart phone, a tablet computer, a laptop computer, a desktop computer, or a personal digital assistant (PDA), and the electronic device may be set in a device cluster.

[0058] Next, the signal parameter extraction method under noise interference provided by the present application is introduced.

[0059] First, see Figure 1 , Figure 1 A flow chart of a method for extracting signal parameters under noise interference in the present application is shown. The method for extracting signal parameters under noise interference provided in the present application may specifically include the following steps S101 to S104:

[0060] Step S101, acquiring an NMR signal collected under NMR detection processing, where the NMR signal is formed by superposition of a target signal and a noise signal;

[0061] It can be understood that the extraction of real signals and target signals to be achieved in this application is aimed at the specific application scenarios of NMR detection.

[0062] NMR detection, that is, nuclear magnetic resonance processing, can also be understood as the processing involved in NMR signals in the magnetic resonance sounding (MRS) method. It has many application scenarios in detection. Generally speaking, the principle of nuclear magnetic resonance is applied. Since the atomic nuclei in different substances are different and have different magnetic gyrometry ratios, the nuclear magnetic resonance effect caused by the differences in the relaxation properties of the atomic nuclei of different substances can be used to achieve the purpose of detection.

[0063] In this regard, the corresponding NMR signal can be collected through the receiving coil configured by the NMR detection processing. In the NMR signal, the required target signal and the noise signal are superimposed together, which can be expressed as:

[0064] S(t)=E(t)+w(t),

[0065] Where S(t) represents the NMR signal, E(t) represents the target signal, and w(t) represents the noise signal.

[0066] Therefore, there is a need for denoising, which requires filtering out noise signals to obtain the target signals, provide data support for subsequent detection and processing of related objects, and for use in data analysis.

[0067] As for step S101, in a specific application, the NMR signal collected by the NMR detection and processing can be retrieved from a local storage space or from other devices, or the NMR signal can be collected locally in real time through NMR detection and processing, or the NMR signal can be collected in real time through NMR detection and processing by other devices. The specific details can be adjusted according to the actual situation and are not specifically limited here.

[0068] Step S102, constructing a second formula when the target signal is expressed by the first formula, where the first formula is: Among them, E(t) s represents the target signal, E s is the initial amplitude of the signal, T s is the average decay time, f t is the preset Larmor frequency, is the initial phase, and the second formula is: Where Z1(t) and Z2(t) are E(t) s Two signals with correlation are constructed, S(t) represents NMR signal, a and b are real numbers;

[0069] It should be understood that, in order to meet the denoising needs of NMR signals, compared with the prior art which estimates and eliminates noise signals, the present application bypasses the noise signals, overcomes the problem that it is difficult to completely eliminate the noise signals, and directly determines the target signals.

[0070] Specifically, in the following content, the present application uses the groundwater detection scenario based on NMR detection processing as an example.

[0071] In the NMR detection process, the hydrogen nuclei in the groundwater are excited by the alternating current pulses of a certain frequency to form a macroscopic magnetic moment. During the relaxation of the magnetic moment, a nuclear magnetic resonance signal is generated. The receiving coil located on the ground receives a free induction decay signal, namely the NMR signal.

[0072] refer to Figure 2 A signal schematic diagram of the NMR signal of the present application is shown, in which 0-t1 is the pure noise acquisition time before the instrument is excited; t1-t2 is the duration of the excitation pulse of the instrument transmitting the Larmor frequency, and the pulse envelope is a rectangle; t2-t4 is the transmission to reception time, called the dead time; t4-t5 is the NMR signal receiving time, and the NMR signal envelope decays according to the exponential law.

[0073] It can be understood that the present application believes that when the expression of the target signal can be written, a signal with multiple unknown parameters and the same form as the target function can be constructed. The value range of each unknown parameter can be given through known prior information to ensure that the real parameters of the target signal are within the corresponding range. Specifically, the target signal can be configured as follows:

[0074]

[0075] in,

[0076] E(t) s Indicates the target signal,

[0077] E s is the initial amplitude of the signal, which is related to the underground water content.

[0078] T s is the average decay time, i.e., the spin-spin relaxation time, which is closely related to the pore size of the aquifer. The common T s The approximate relationship with the corresponding aquifer can also be referred to the following Table 1.

[0079] f t is the preset Larmor frequency, which is related to the strength of the geomagnetic field and can be obtained through field measurements. It can be understood as a relatively unchanged factor in this NMR detection process.

[0080] is the initial phase, which is the phase shift of the secondary field relative to the excitation current and is related to the conductivity of the underground rock formation.

[0081] Table 1 - Measured T s Approximate relationship between values ​​and aquifer types

[0082]

[0083] NMR signals are alternating electromagnetic field signals, and are very weak, usually between tens and thousands of nanovolts, and are easily interfered by external electromagnetic noise. Therefore, when performing measurements, repeated experiments are often performed for later signal superposition to improve the signal-to-noise ratio, such as Figure 3 Another signal schematic diagram of the NMR signal of the present application is shown in Figure 3 In the figure, from top to bottom are the results of 5, 10, 15, and 25 repeated measurement data superpositions. The left side is the time domain image and the right side is the frequency domain image. The result of 5 superpositions does not show the exponential decay law. Only after multiple superpositions can the exponential decay law be seen, but the noise still submerges the real target signal E(t). s , the target signal E cannot be known s , Ts , No matter what the parameters are, it is impossible to accurately reflect the real target signal E(t) s , it is impossible to obtain the corresponding groundwater detection results.

[0084] Therefore, the signal extraction idea of ​​the present application can be applied for processing, that is, bypassing the noise signal, directly targeting the characteristics of the target signal itself, and performing correlation analysis on the basis of signal superposition and cancellation to achieve accurate extraction of the target signal.

[0085] First, it involves constructing E(t) s The two signals Z1(t) and Z2(t) with a correlation relationship have the following relationship, that is, the cross-correlation function expression obtained by performing cross-correlation calculation on the two variables:

[0086]

[0087] Among them, a and b are real numbers to be taken.

[0088] After constructing the cross-correlation function equation group, the specific parameters of the target signal can be solved through subsequent mathematical operations, so as to restore the target signal.

[0089] Step S103, taking values ​​of a and b, and bringing them into the third formula to calculate the cross-correlation function of Z1(t) and Z2(t), the third formula is: Where R(τ) represents the cross-correlation function;

[0090] After constructing the cross-correlation function equation group, we can take the values ​​of a and b real numbers and bring them into the cross-correlation function equation group, and then solve the cross-correlation function of Z1(t) and Z2(t), that is:

[0091]

[0092] The cross-correlation function is expanded as:

[0093]

[0094] At this time, the cross-correlation function R(τ) leaves the unknown E s , T s , Three.

[0095] Step S104, performing E on the cross-correlation function s , T s , , and bring the solution into the first formula to restore the target signal.

[0096] Based on the cross-correlation function R(τ) above, it can be solved based on the mathematical operation of the configuration and combined with the specific conditions of the configuration, so that E can be solved. s , T s , Third, substitute the above target signal expression into the first formula: The target signal can be restored and separated from the NMR signal mixed with noise signals.

[0097] It is worth mentioning that the formula expression of the cross-correlation function R(τ) is a theoretical formula specially configured for this application. In actual work, there is no signal with infinite duration, but when T is not infinite, the formula is still valid.

[0098] from Figure 1 It can be seen from the illustrated embodiments that, with respect to the extraction of the real target signal in the NMR signal, compared with the prior art in which the noise signal of the NMR signal is estimated and eliminated, the present application bypasses the noise signal, overcomes the problem that the noise signal is difficult to be approximately completely eliminated, and directly targets the characteristics of the target signal itself, performs correlation analysis on the basis of signal superposition and cancellation, to achieve accurate extraction of the target signal, thereby promoting a higher-precision NMR detection effect.

[0099] Next, we will continue to Figure 1 Each step of the illustrated embodiment and possible implementation methods thereof in practical applications are described in detail.

[0100] Regarding the solution process in step S104, as a practical implementation method of the present application, after performing E s , T s , The solution can include the following:

[0101] The cross-correlation function is E- s , T s , The preset mathematical operation may specifically include at least one of the mathematical operation strategies such as addition, subtraction, multiplication, division, taking absolute value, taking extreme value, taking norm, taking derivative and partial derivative.

[0102] Taking the norm as an example, as another practical implementation method, the values ​​of a and b in the previous step S103 can be specifically two groups (that is, after taking the values, one group of a1 and b1 is obtained, and another group of a2 and b2 is obtained). After the two groups of a and b are respectively substituted into the third formula, two cross-correlation functions R can be obtained. 12 (τ) and R 34 (τ), and in step S104 the cross-correlation function is Es , T s , The solution process can include the following contents:

[0103] R 12 (τ) and R 34 (τ), and find the 2-norm respectively, and get ||R 12 (τ)||2 and||R 34 (τ)||2;

[0104] Ask||R 12 (τ)||2 and||R 34 (τ)||2, we get

[0105] The average decay time range from 0 to 1 second (prior information, corresponding to T s ), the phase range from 0 to 2π (prior information, corresponding to ) and the initial amplitude range extracted from the NMR signal (it can also be extracted from the content shown in the signal waveform. The initial amplitude range is an inherent characteristic of the signal, corresponding to E s ) are brought into G to solve, and the solution is E when G reaches the maximum value. s , T s , The value of .

[0106] It can be seen that under this setting, E s , T s , The solution process provides a set of specific implementation solutions, which is not only easy to apply, but also found in the actual research and development process that it can effectively obtain E s , T s , In other words, the target signal can be restored with high precision.

[0107] In addition, it can be understood that the values ​​of a and b involved can be selected randomly and arbitrarily in principle, and can be selected according to the specific situation of the noisy signal, which can be used for mathematical operations of subsequent correlation analysis to solve E s , T s , That's it.

[0108] Furthermore, it can be seen from the exemplary embodiment here that when a and b are taken as values, not only one group but multiple groups can be taken. How to combine multiple groups of a and b to provide mathematical operations for subsequent related analysis can be determined based on the computational requirements of the mathematical operations for subsequent related analysis.

[0109] In the case of Es, Ts, When solving , the present application takes two sets of values, one set can be: a = 1, b = 2, then:

[0110]

[0111] The cross-correlation function is:

[0112]

[0113] After expansion, there are:

[0114]

[0115] Another group can be: a=-1, b=1, then we have:

[0116]

[0117] The cross-correlation function is:

[0118]

[0119] Later you can bring in Conduct E s , T s , The solution.

[0120] As an example, the Larmor frequency f t The preset value is 2136.8Hz, and the initial amplitude E s like Figure 3 The content shown must be less than 400nV. At this time, the ranges of each parameter are determined as follows: Es∈(0,400), Ts∈(0,1), After solving the three problems, the result is E s =263.2774nV, T s =0.1766s, is 1.71π.

[0121] In addition, in the process of solving G, that is, in E s , T s , In the process of solving the three, related solving algorithms and search algorithms may also be involved. As another practical implementation method, in the specific application of this application, an improved artificial fish swarm algorithm, a simulated annealing method or a genetic algorithm may be used for search processing. Of course, other different types of search algorithms may also be used to solve G, that is, E s , T s , The solution of the three.

[0122] Back to the specific application scenario, the signal parameter extraction method under noise interference provided in the present application may also include data processing procedures related to the application scenario.

[0123] Taking the application scenario of non-destructive detection of groundwater as an example, NMR detection can utilize the changing law of nuclear magnetic resonance signals generated by hydrogen nuclei (protons) in underground aquifers to detect the existence of underground aquifers and the burial depth and water content of groundwater, thereby achieving the purpose of depth measurement.

[0124] Correspondingly, the NMR signal involved in the above step S101 may specifically be a signal collected by NMR detection during groundwater detection.

[0125] In this case, as another practical implementation, before step S101, the following contents may also be involved:

[0126] Initiate NMR detection processing in the target area.

[0127] It can be understood that the target area is an area where groundwater exists. Thus, under this setting, when the electronic device executes the signal parameter extraction method under noise interference of the present application, Figure 4 A scenario schematic diagram of the signal parameter extraction method under noise interference of the present application is shown, which may also involve initiating NMR detection processing through the NMR equipment configured by itself, so that the NMR signal can be obtained by itself, and then a high-precision target signal can be determined through subsequent data processing for use in groundwater detection and analysis in the target area.

[0128] Furthermore, the electronic device can also perform the final groundwater detection and analysis task, that is, after performing E on the cross-correlation function in step S104, s , T s , After solving the problem and substituting the solution result into the first formula to restore the target signal, the signal parameter extraction method under noise interference of the present application may also include:

[0129] Based on the changing patterns reflected by different target signals restored through different NMR detection processes, the groundwater location, groundwater burial depth and groundwater content are determined.

[0130] It can be understood that the groundwater detection analysis (groundwater location, groundwater burial depth and groundwater water content) is determined in a dynamic manner, that is, it is carried out on the basis of different rounds of NMR detection processing results (corresponding to the extraction of multiple target signals, involving multiple executions of steps S101 to S104). In this way, based on the changing laws reflected by different target signals, the groundwater location, groundwater burial depth and groundwater content are ultimately accurately determined.

[0131] As for the changing pattern of the target signal and the corresponding relationship between the groundwater position, groundwater burial depth and groundwater water content, it is similar to the prior art and will not be elaborated here. Of course, the corresponding relationship adopted in this application can also be optimized and improved in specific applications, which will not be elaborated here.

[0132] The above is an introduction to the signal parameter extraction method under noise interference provided by the present application. In order to facilitate better implementation of the signal parameter extraction method under noise interference provided by the present application, the present application also provides a signal parameter extraction device under noise interference from the perspective of functional modules.

[0133] See also Figure 5 , Figure 5 Schematic diagram of a signal parameter extraction device under noise interference in this application. In this application, the signal parameter extraction device 500 under noise interference may specifically include the following structure:

[0134] An acquisition unit 501 is used to acquire an NMR signal collected under NMR detection processing, where the NMR signal is formed by superposition of a target signal and a noise signal;

[0135] The constructing unit 502 is used to construct a second formula when the target signal is expressed by the first formula, where the first formula is: Among them, E(t) s represents the target signal, E s is the initial amplitude of the signal, T s is the average decay time, f t is the preset Larmor frequency, is the initial phase, and the second formula is: Wherein, Z1(t) and Z2(t) are two signals with correlation constructed by E(t)s, S(t) represents NMR signal, and a and b are real numbers;

[0136] Substitute into unit 503, used to take values ​​of a and b, and substitute into the third formula to calculate the cross-correlation function of Z1(t) and Z2(t), the third formula is: Where R(τ) represents the cross-correlation function;

[0137] Solving unit 504, used to perform Es, Ts, , and bring the solution into the first formula to restore the target signal.

[0138] In an exemplary implementation, the solving unit 504 is specifically configured to:

[0139] The cross-correlation function is E- s , T s , The preset mathematical operations include at least one of addition, subtraction, multiplication, division, absolute value, extreme value, norm, derivative and partial derivative.

[0140] In another exemplary implementation, the values ​​of a and b are specifically two groups. After the two groups of a and b are substituted into the third formula, two cross-correlation functions R are obtained. 12 (τ) and R 34 (τ), the solving unit 504 is specifically used for:

[0141] R 12 (τ) and R 34 (τ), and find the 2-norm respectively, and get ||R 12 (τ)||2 and||R 34 (τ)||2;

[0142] Ask||R 12 (τ)||2 and||R 34 (τ)||2, we get

[0143] The average decay time range of 0 to 1 second, the phase range of 0 to 2π, and the initial amplitude range extracted from the NMR signal are brought into G to solve the problem of E when G reaches its maximum value. s , T s , The value of .

[0144] In yet another exemplary implementation, in the process of solving G, an improved artificial fish swarm algorithm, a simulated annealing method or a genetic algorithm is used for search processing.

[0145] In yet another exemplary implementation, the NMR signal is specifically a signal collected by NMR detection during groundwater detection.

[0146] In yet another exemplary implementation, the apparatus further includes an initiating unit 505, configured to:

[0147] Initiate NMR detection processing in the target area.

[0148] In yet another exemplary implementation, the apparatus further includes a determining unit 506, configured to:

[0149] Based on the changing patterns reflected by different target signals restored through different NMR detection processes, the groundwater location, groundwater burial depth and groundwater content are determined.

[0150] The present application also provides an electronic device from the perspective of hardware structure, see Figure 6 , Figure 6 601, a memory 602, and an input / output device 603. The processor 601 is used to execute the computer program stored in the memory 602 to implement the following Figure 1 The steps of the signal parameter extraction method under noise interference in the corresponding embodiment; or, the processor 601 is used to execute the computer program stored in the memory 602 to implement the following Figure 5 Corresponding to the functions of each unit in the embodiment, the memory 602 is used to store the processor 601 executing the above Figure 1 The computer program required by the signal parameter extraction method under noise interference in the corresponding embodiment.

[0151] Exemplarily, the computer program may be divided into one or more modules / units, one or more modules / units are stored in the memory 602, and executed by the processor 601 to complete the present application. One or more modules / units may be a series of computer program instruction segments capable of completing specific functions, and the instruction segments are used to describe the execution process of the computer program in the computer device.

[0152] The electronic device may include, but is not limited to, a processor 601, a memory 602, and an input / output device 603. Those skilled in the art will appreciate that the illustration is merely an example of an electronic device and does not constitute a limitation on the electronic device, and may include more or fewer components than shown in the illustration, or a combination of certain components, or different components. For example, the electronic device may also include a network access device, a bus, etc., and the processor 601, the memory 602, the input / output device 603, etc. are connected via a bus.

[0153] The processor 601 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc. The processor is the control center of the electronic device, and uses various interfaces and lines to connect various parts of the entire device.

[0154] The memory 602 can be used to store computer programs and / or modules. The processor 601 implements various functions of the computer device by running or executing the computer programs and / or modules stored in the memory 602 and calling the data stored in the memory 602. The memory 602 can mainly include a program storage area and a data storage area, wherein the program storage area can store an operating system, at least one application required for a function, etc.; the data storage area can store data created according to the use of the electronic device, etc. In addition, the memory can include a high-speed random access memory, and can also include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (SecureDigital, SD) card, a flash card (Flash Card), at least one disk storage device, a flash memory device, or other volatile solid-state storage devices.

[0155] When the processor 601 is used to execute the computer program stored in the memory 602, the following functions can be implemented:

[0156] Acquiring an NMR signal collected under NMR detection processing, where the NMR signal is formed by superposition of a target signal and a noise signal;

[0157] When the target signal is expressed by the first formula, the second formula is constructed. The first formula is: Among them, E(t) s represents the target signal, E s is the initial amplitude of the signal, T s is the average decay time, f t is the preset Larmor frequency, is the initial phase, and the second formula is: Wherein, Z1(t) and Z2(t) are two signals with correlation constructed by E(t)s, S(t) represents NMR signal, and a and b are real numbers;

[0158] Take values ​​for a and b and bring them into the third formula to calculate the cross-correlation function of Z1(t) and Z2(t). The third formula is: Where R(τ) represents the cross-correlation function;

[0159] Es, Ts, , and bring the solution into the first formula to restore the target signal.

[0160] Those skilled in the art can clearly understand that for the convenience and simplicity of description, the specific working process of the signal parameter extraction device under noise interference, the electronic device and its corresponding units described above can refer to the following. Figure 1 The description of the signal parameter extraction method under noise interference in the corresponding embodiment will not be repeated here.

[0161] A person of ordinary skill in the art will appreciate that all or part of the steps in the various methods of the above embodiments may be completed by instructions, or by controlling related hardware through instructions. The instructions may be stored in a computer-readable storage medium and loaded and executed by a processor.

[0162] To this end, the present application provides a computer-readable storage medium, in which a plurality of instructions are stored, and the instructions can be loaded by a processor to execute the present application as follows: Figure 1 The steps of the signal parameter extraction method under noise interference in the corresponding embodiment, the specific operation can refer to the following Figure 1 The description of the signal parameter extraction method under noise interference in the corresponding embodiment will not be repeated here.

[0163] The computer-readable storage medium may include: a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, etc.

[0164] Due to the instructions stored in the computer-readable storage medium, the present application can be executed. Figure 1 The steps of the signal parameter extraction method under noise interference in the corresponding embodiment, therefore, the present application can be implemented as follows Figure 1 The beneficial effects that can be achieved by the signal parameter extraction method under noise interference in the corresponding embodiment are detailed in the previous description and will not be repeated here.

[0165] The signal parameter extraction method, device, electronic device and computer-readable storage medium under noise interference provided by the present application are introduced in detail above. Specific examples are used in this article to illustrate the principles and implementation methods of the present application. The description of the above embodiments is only used to help understand the method of the present application and its core idea; at the same time, for technical personnel in this field, according to the idea of ​​the present application, there will be changes in the specific implementation method and application scope. In summary, the content of this specification should not be understood as a limitation on the present application.

Claims

1. A method for extracting signal parameters under noise interference, characterized in that: The method comprises: Acquiring an NMR signal collected under a nuclear magnetic resonance (NMR) detection process, wherein the NMR signal is formed by superposition of a target signal and a noise signal; When the target signal is expressed by the first formula, a second formula is constructed, where the first formula is: Among them, E(t) s represents the target signal, E s is the initial amplitude of the signal, T s is the average decay time, f t is the preset Larmor frequency, is the initial phase, the second formula is: Wherein, Z1(t) and Z2(t) are two signals with correlation constructed by E(t)s, S(t) represents the NMR signal, and a and b are real numbers; Take values ​​for a and b, and bring them into the third formula to calculate the cross-correlation function of Z1(t) and Z2(t). The third formula is: Wherein, R(τ) represents the cross-correlation function; The cross-correlation function is subjected to Es, Ts, and bring the solution into the first formula to restore the target signal.

2. The method according to claim 1, characterized in that: The cross-correlation function is subjected to E s , T s , The solution includes: The cross-correlation function is subjected to E based on a preset mathematical operation. s 、T s , The preset mathematical operation includes at least one of addition, subtraction, multiplication, division, absolute value, extreme value, norm, derivative and partial derivative.

3. The method according to claim 1, characterized in that The values ​​of a and b are specifically two groups. After the two groups of a and b are substituted into the third formula, the two cross-correlation functions R are obtained. 12 (τ) and R 34 (τ), the cross-correlation function is subjected to E s , T s , The solution includes: R 12 (τ) and R 34 (τ), and find the 2-norm respectively, and get ||R 12 (τ)||2 and||R 34 (τ)||2; Ask||R 12 (τ)||2 and||R 34 (τ)||2, we get The average decay time range of 0 to 1 second, the phase range of 0 to 2π, and the initial amplitude range extracted from the NMR signal are brought into G to solve, and the solution is E when G reaches the maximum value. s 、T s , The value of .

4. The method according to claim 3, characterized in that In the process of solving G, an improved artificial fish swarm algorithm, simulated annealing method or genetic algorithm is used for search processing.

5. The method according to claim 1, characterized in that The NMR signal is specifically a signal collected by NMR detection during groundwater detection.

6. The method according to claim 5, characterized in that Before obtaining and processing the NMR signal acquired under NMR detection, the method further comprises: The NMR detection process is initiated at a target region.

7. The method according to claim 5, characterized in that The cross-correlation function is subjected to Es, Ts, After solving the problem and bringing the solution into the first formula to restore the target signal, the method further includes: Based on the changing patterns reflected by the different target signals restored through different NMR detection processes, the groundwater location, groundwater burial depth and groundwater content are determined.

8. A signal parameter extraction device under noise interference, characterized in that: The device comprises: An acquisition unit, used for acquiring an NMR signal collected in a nuclear magnetic resonance (NMR) detection process, wherein the NMR signal is formed by superposition of a target signal and a noise signal; A constructing unit is used to construct a second formula when the target signal is expressed by a first formula, wherein the first formula is: Among them, E(t) s represents the target signal, E s is the initial amplitude of the signal, T s is the average decay time, f t is the preset Larmor frequency, is the initial phase, the second formula is: Wherein, Z1(t) and Z2(t) are two signals with correlation constructed by E(t)s, S(t) represents the NMR signal, and a and b are real numbers; The unit is used to take values ​​of a and b, and is used to calculate the cross-correlation function of Z1(t) and Z2(t) in the third formula. The third formula is: Wherein, R(τ) represents the cross-correlation function; A solving unit is used to perform Es, Ts, and bring the solution into the first formula to restore the target signal.

9. An electronic device, characterized in that: The method comprises a processor and a memory, wherein the memory stores a computer program, and the processor executes the method according to any one of claims 1 to 7 when calling the computer program in the memory.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a plurality of instructions, and the instructions are suitable for being loaded by a processor to execute the method according to any one of claims 1 to 7.

Images