A charging method data denoising method and system in a strong interference environment

Abstract

The present application relates to the technical field of geophysical electromagnetic exploration, and discloses a charging method data denoising method and system in a strong interference environment, the method comprising: obtaining the transmission data transmitted by the transmission end and the response data received by the receiving end; determining the position range of the target potential difference in the main cycle according to the time difference between the starting time of the two; for the response data, dividing it into multiple main cycle segments according to the main cycle length of the transmission data waveform, superimposing the multiple main cycle segments and calculating the average value to obtain the superimposed main cycle segment; according to the step detection function, identifying the step of the superimposed main cycle segment, and determining the position of the target potential difference step in combination with the position range of the target potential difference; and calculating the target potential difference value according to the potential data before and after the target potential difference step. The present application can effectively suppress noise, improve the signal-to-noise ratio of the data, and extract the real and accurate potential difference value for the charging method data subjected to strong interference.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geophysical electromagnetic exploration, in particular to a charging method data denoising method and system in a strong interference environment. BACKGROUND

[0002] The statements in this section merely provide background information related to the present application and do not necessarily constitute prior art.

[0003] Charging method is a method of charging a good conductor such as ore body or aquifer, and understanding the size, occurrence state and even movement state (when the target good conductor is a flowing conductor such as underground water) of the target good conductor by observing the spatial distribution of its charging field. The observation methods of charging method mainly include potential method, potential gradient method and tracing equipotential line method. Usually, charging method uses square wave to excite signals, and two power supply electrodes are located in the measurement area and at infinity respectively. The data is acquired by potential observation (measuring the potential difference between the measurement points in the measurement area and the position at infinity) and potential difference observation (measuring the potential difference between two measurement points in the measurement area). The potential information of all measurement points in the measurement area is obtained by using a relatively small amount of potential data combined with a relatively large amount of potential difference data and the mathematical and physical relationship between potential difference and potential, and then the potential distribution map of the measurement area is drawn. In actual field work, potential difference data observation can be used instead of potential data observation. The transmitting end is fixed, and data is received at a plurality of measurement points in turn by arranging a plurality of measurement points. The potential difference between these measurement points and the position of the transmitting end is calculated. Not only the gradient curve with strong resolution can be obtained, but also the potential distribution of the measurement area can be calculated by potential difference. The device is light and easy to implement.

[0004] However, in the application process of charging method, it is always affected by various electromagnetic noise interference, especially in a strong interference environment, which leads to the inability to correctly extract potential difference information from the collected data. The error of potential difference information will lead to the error of subsequent calculation. The existing operation procedures and related documents mostly only mention the evaluation method of charging method data quality, and few mention the method or means for processing the charging method data affected by strong interference and extracting usable potential difference information. SUMMARY

[0005] The present application provides an automatic denoising processing method for charging method data in a strong interference environment. The charging method data affected by strong interference is processed to effectively suppress noise, improve data signal-to-noise ratio, and extract relatively true and accurate potential difference values.

[0006] To achieve the above purpose, one or more embodiments of the present application provide the following technical solutions:

[0007] A charging method data denoising method in a strong interference environment, characterized in that it comprises the following steps:

[0008] acquire the transmission data transmitted by the transmission end and the response data received by the reception end;

[0009] determine a position range of the target potential difference in a main period according to a time difference between the start time of the transmission data and the response data;

[0010] divide the response data into a plurality of main period segments according to the length of the main period of the transmission data waveform, superimpose the plurality of main period segments and calculate the average value to obtain the superimposed main period segment;

[0011] perform step detection on the superimposed main period segment according to a step detection function, and determine the position of the step of the target potential difference in combination with the position range of the target potential difference;

[0012] calculate the target potential difference value according to the potential data before and after the step of the target potential difference.

[0013] Further, the time difference between the start time of the transmission data and the response data is the length of the main period of the transmission data waveform is L, and the theoretical position of the target potential difference in the main period is:

[0014]

[0015] extend a certain time period forward and backward from the theoretical position to obtain the position range of the target potential difference in the main period.

[0016] Further, the step detection function is created according to the length L of the main period of the transmission data waveform, and specifically:

[0017]

[0018] Further, the step detection on the superimposed main period segment includes:

[0019] perform bidirectional continuation on the superimposed main period segment;

[0020] perform step detection on the superimposed main period segment after bidirectional continuation according to the step detection function to obtain a step detection result;

[0021] obtain a step indication function by taking the two ends of the step detection result and the parts corresponding to bidirectional continuation, and the wave peak of the step indication function corresponds to the step position in the superimposed main period segment.

[0022] Further, the step detection on the superimposed main period segment after bidirectional continuation includes:

[0023] perform convolution operation on the superimposed main period segment after bidirectional continuation and the step detection function.

[0024] The absolute value of the difference between the convolution operation result and the average value thereof is calculated to obtain a step detection result.

[0025] Further, the position of the target potential difference step is determined in combination with the position range in which the target potential difference is located.

[0026] According to the wave peak of the step indication function, one or more step positions in the superimposed main period section are determined.

[0027] The wave peak with the largest amplitude in the position range is recorded as the position of the target potential difference step.

[0028] Further, the calculation of the potential difference value includes:

[0029] The positions of the target sampling points are taken on the left and right sides of the position of the target potential difference step at a distance of a set time length from the center.

[0030] A plurality of data points are obtained around the two target sampling points, and the average value is calculated to obtain the average values of the potentials before and after the target potential difference step.

[0031] The average values of the potentials before and after the target potential difference step are subtracted to obtain the target potential difference value.

[0032] One or more embodiments provide a charging method data denoising system in a strong interference environment, comprising:

[0033] A data acquisition module is configured to acquire transmission data transmitted by a transmission end and response data received by a receiving end.

[0034] A position range determination module is configured to determine a position range of a target potential difference in a main period according to a time difference between the transmission data and the response data.

[0035] A period superposition module is configured to divide the response data into a plurality of main period sections according to the length of the main period of the transmission data waveform, superimpose the plurality of main period sections, and calculate the average value to obtain a superimposed main period section.

[0036] A step position determination module is configured to perform step identification on the superimposed main period section according to a step detection function, and determine the position of the target potential difference step in combination with the position range in which the target potential difference is located.

[0037] A potential difference value calculation module is configured to calculate a target potential difference value according to the potential data before and after the target potential difference step.

[0038] One or more embodiments provide an electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the method for data denoising of the charging method in a strong interference environment when executing the program.

[0039] One or more embodiments provide a computer readable storage medium, having stored thereon a computer program, which is executed by a processor to implement the method for data denoising of the charging method in a strong interference environment.

[0040] The above one or more technical solutions have the following beneficial effects:

[0041] By superimposing the received response data according to the period, various random noises in the received data can be suppressed, the proportion of the target potential difference in the data is relatively enhanced, which is helpful for subsequent accurate identification of the target potential difference, and the accuracy of the target potential difference value is ensured.

[0042] In the step detection stage, by bidirectional extension of the main period section, the boundary effect in step detection can be eliminated, and the interference peak in the subsequent step indicator function can be avoided, thereby improving the accuracy of target potential difference identification.

[0043] After obtaining the step position of the target potential difference, sampling points are taken before and after the step position, a plurality of potential data around each sampling point is obtained, the mean value of the potential data before and after the step position is obtained, and the potential difference value is calculated. Compared with using a point value as the potential data before and after the step position, the accuracy is higher, and the utilization rate of the charging method data in a strong interference environment is improved.

[0044] The application can effectively suppress the noise in the data, improve the signal-to-noise ratio of the data, and can be widely applied to the charging method exploration in a strong interference environment. Moreover, after the required parameters are input in advance, manual intervention is no longer needed, and all operations can be performed by computer software, which is conducive to rapid and automatic processing of a large amount of charging method data. BRIEF DESCRIPTION OF DRAWINGS

[0045] The drawings accompanying the specification of the present application form a part thereof, serve to provide further understanding of the present application, and together with the specification explain the exemplary embodiments of the present application and do not limit the present application.

[0046] Figure 1 The figure is a schematic diagram of an observation scheme for a charging method field work in a traditional case, wherein A is a power supply electrode in a measurement area, another power supply electrode connected to an infinite position, V with a circle represents a receiving instrument for recording measurement data, most of which is used to measure the potential difference between adjacent measurement points, and a small part is used to measure the potential of the measurement point (i.e. the potential difference between the measurement point and the infinite position);

[0047] Figure 2 FIG. 1 is a flowchart of a charging method data denoising method in a strong interference environment according to an embodiment of the present application;

[0048] Figure 3 FIG. 2 is a response data waveform example of a charging method in a strong interference environment according to an embodiment of the present application;

[0049] Figure 4 FIG. 3 is a waveform diagram of a main cycle segment obtained after period superposition of response data of a charging method according to an embodiment of the present application; Figure 3

[0050] Figure 5 FIG. 4 is a comparison diagram of a main cycle segment waveform and a step indicator function after superposition according to an embodiment of the present application; Figure 4

[0051] Figure 6 FIG. 5 is a diagram of extracting a potential difference value from a main cycle segment after superposition according to an embodiment of the present application. Figure 4 DETAILED DESCRIPTION

[0052] The present application will be further described below in conjunction with the accompanying drawings and embodiments.

[0053] It should be noted that the following detailed description is illustrative only and is intended to provide further description of the present application. Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which the present application pertains.

[0054] It should be noted that the terms used herein are only intended to describe specific embodiments and are not intended to limit exemplary embodiments according to the present application. As used herein, the singular form is intended to include the plural form unless the context clearly indicates otherwise, and it should also be understood that when the terms "comprise" and / or "include" are used in this specification, they indicate the presence of the features, steps, operations, devices, components, and / or combinations thereof.

[0055] The embodiments in the present application and the features in the embodiments can be combined with each other without conflict.

[0056] Embodiment One

[0057] The present embodiment discloses a charging method data denoising method in a strong interference environment, which uses period superposition, step identification, potential difference value extraction, and other means to process charging method data disturbed by noise. As shown in FIG. 1, the method comprises the following steps: Figure 2

[0058] Step 1: Obtain the transmission data transmitted by the transmission end and the response data received by the receiving end, and determine the position range of the target potential difference in the main cycle according to the difference between the start time of the transmission data and the response data.​​​​

[0059] Generally, the main period refers to the period corresponding to the minimum main frequency in the transmission waveform (there can be multiple main frequencies in the transmission waveform). Let the time when the transmission end starts to work be , the time when the reception end starts to work be , the difference between the two be , and the main period of the transmission waveform be , then the theoretical position of the step corresponding to the target potential difference on the superimposed main period waveform of the reception can be obtained by the following formula:

[0060] (1)

[0061] In actual cases, due to the inaccuracy of time recording, the range in which the target potential difference step exists can be appropriately enlarged. After the theoretical position is calculated by the above formula, according to the maximum error of the time recording, a certain time period is extended forward and backward from the theoretical position to obtain the position range of the target potential difference in the main period.

[0062] If there is no time recording, the approximate range in which the target step exists in all reception data can also be determined by observing and comparing the phase difference distribution of the transmission data and individual reception data according to the actual position of the target step in the individual reception data.

[0063] In the step 1, only the approximate range in which the target potential difference step exists in the main period needs to be determined, and it does not need to be very accurate. In the case where the data is weakly disturbed, the range can also be set as the entire main period.

[0064] Step 2: For the response data, divide it into multiple main period segments according to the main period length of the transmission data waveform, superimpose the amplitudes of the multiple main period segments and calculate the average value to obtain the superimposed main period segment.

[0065] By superimposing the received response data in cycles, various random noises in the reception data can be suppressed, and the proportion of the target potential difference in the data can be relatively enhanced.

[0066] Specifically, let the actually received charging method data be , and the main period length L of the transmission waveform be the time window length to truncate the reception waveform into a plurality of main period segments . Calculate the superimposed average value of all main period segments according to formula (2) to obtain the superimposed main period segment:

[0067] (2)

[0068] The original response data and the superimposed main period segment are shown in Figure 3 and Figure 4 , respectively.

[0069] Step 3: Based on the step detection function, perform step identification on the superimposed main period segment, and determine the position of the target potential difference step by combining the position range of the target potential difference.

[0070] Step 3 specifically includes:

[0071] Step 3.1: Based on the main period length of the transmitted waveform Construct the step detection function according to equation (3). .

[0072] (3)

[0073] Step 3.2: For the superimposed main period segment Perform bidirectional extension, in Each side extends outwards using the average of its first and last 10 data points as the extension value. The length of the superimposed principal period segment is obtained by calculating the length of the extended segment. .

[0074] Step 3.3: Extend the superimposed main periodic segment With step detection function Perform convolution and remove the absolute value after taking the mean, i.e., perform the processing according to equations (4) and (5) to obtain the detection result. :

[0075] (4)

[0076] (5)

[0077] Step 3.4: [The text appears to be incomplete and contains several grammatical errors. A more accurate translation would require Delete both ends The length of the step indicator function is determined to ensure the correspondence between the step indicator function and the superimposed main period segment; the remaining part is called the step indicator function. , The peak in the middle corresponds to In a step jump, the magnitude of the peak amplitude represents the height of the step jump.

[0078] Based on the predetermined range of the target potential difference, by selecting within that range The location of the target potential difference step can be automatically determined by the peak with the largest amplitude. Specifically, based on the peak of the step indicator function, one or more step positions in the superimposed main period segment are determined; the peak with the largest amplitude within the range of these positions is recorded as the location of the target potential difference step. Figure 5 This is a comparison diagram of the step indicator function and the superimposed main period segment.

[0079] Step 4: calculating the target potential difference value according to the potential data before and after the target potential difference step.

[0080] To ensure the accuracy of the acquired potential difference value, positions with a distance of the set time length from the center to the left and right of the center are taken as target sampling points, in this embodiment, positions with a distance of 5-20 time lengths from the center to the left and right of the center are taken as target sampling points. The average value of 5-20 data points around the two target sampling points is taken as the average value of the potential before and after the target potential difference step, and the target potential difference value is obtained by subtraction. Further, the purpose of suppressing noise, improving the data signal-to-noise ratio, and extracting a relatively true and accurate potential difference value is achieved.

[0081] Thus, the calculation of the potential difference between a measurement point and the position of the transmitting end is realized. Based on the above method, the potential difference between other measurement points and the position of the transmitting end is calculated.

[0082] Embodiment Two

[0083] Based on the method provided in Embodiment One, this embodiment provides a charging method data denoising system in a strong interference environment, the system comprising:

[0084] a data acquisition module, configured to acquire transmitting data transmitted by a transmitting end and response data received by a receiving end;

[0085] a position range determination module, configured to determine the position range of a target potential difference in a main period according to the time difference between the start time of the transmitting data and the response data;

[0086] a period superposition module, configured to divide the response data into a plurality of main period segments according to the main period length of the transmitting data waveform, superimpose the plurality of main period segments, and calculate the average value to obtain a superimposed main period segment;

[0087] a step position determination module, configured to perform step identification on the superimposed main period segment according to a step detection function, and determine the position of the target potential difference step in combination with the position range of the target potential difference;

[0088] a potential difference value calculation module, configured to calculate the target potential difference value according to the potential data before and after the target potential difference step.

[0089] Embodiment Three

[0090] The purpose of this embodiment is to provide an electronic device.

[0091] An electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor executes the program to implement the charging method data denoising method in a strong interference environment as described in Embodiment One.

[0092] Embodiment four

[0093] The purpose of this embodiment is to provide a computer readable storage medium.

[0094] A computer readable storage medium, having stored thereon a computer program, which, when executed by a processor, implements the method for data denoising of charging method in strong interference environment as described in embodiment one.

[0095] The steps and methods involved in the above embodiments two to four correspond to embodiment one, and the specific implementation can refer to the relevant description part of embodiment one.

[0096] Those skilled in the art should understand that the above-mentioned modules or steps of the present application can be realized by a general computer device.

[0097] Although the specific embodiments of the present application are described above in combination with the accompanying drawings, it is not a limitation on the protection scope of the present application, and those skilled in the art should understand that various modifications or changes made by those skilled in the art on the basis of the technical solutions of the present application without creative labor are still within the protection scope of the present application.

Claims

1. A method for data denoising of a charging method in a strong interference environment, characterized in that, The method comprises the following steps: acquiring transmission data transmitted by a transmission end and response data received by a receiving end; determining a position range of a target potential difference in a main period according to a time difference between the transmission data and the start time of the response data; dividing the response data into a plurality of main period segments according to the length of the main period of the transmission data waveform, superimposing the plurality of main period segments and calculating an average value to obtain superimposed main period segments; performing step detection on the superimposed main period segments according to a step detection function, and determining the position of the step of the target potential difference in combination with the position range of the target potential difference; specifically comprising: Step 3.1: Based on the main period length of the transmitted waveform According to the formula Construct a step detection function ; Step 3.2: to the superimposed main cycle section Bidirectional extension is performed, and the average values of the first and last 10 data points on both sides are taken as the extension values, respectively, and each is extended outward by the length to obtain the extended superimposed main cycle section . ; Step 3.3: The extended main period section of the superposition is convolved with the step detection function and the absolute value is taken, i.e. the processing of equation is performed, resulting in the detection result , :​ Step 3.4: The length of the step function is deleted at both ends to ensure the correspondence between the step function and the main period segment after superposition, and the remaining part is called the step function , The length of the step function is deleted at both ends to ensure the correspondence between the step function and the main period segment after superposition, and the remaining part is called the step function , The length of the step function is deleted at both ends to ensure the correspondence between the step function and the main period segment after superposition, and the remaining part is called the step function , According to the range where the predetermined target potential difference is located, the target potential difference is determined by selecting the range The position of the target potential difference step is automatically determined according to the wave peak with the largest amplitude in the range; one or more step positions in the main period section after superposition are determined according to the wave peak of the step indication function; the wave peak with the largest amplitude in the position range is recorded as the position of the target potential difference step. calculating the target potential difference value according to the potential data before and after the step of the target potential difference.

2. The method of claim 1, wherein the method is characterized by, Let the time difference between the start of the transmitted data and the response data be , the main cycle length of the transmitted data waveform is L, and the theoretical position of the target potential difference in the main cycle is: extending the theoretical position forward and backward by a certain time period to obtain the position range of the target potential difference in the main period.

3. The method of claim 1, wherein the method is characterized by, The step detection function is created according to the length L of the main period of the transmission data waveform, and specifically comprises: 。 4. The method of claim 1, wherein the method is characterized by, The calculation of the potential difference value comprises: taking the position of the step of the target potential difference as the center, and taking the positions of a certain length on the left and right sides of the center as target sampling points; acquiring a plurality of data points around the two target sampling points and calculating an average value to obtain the average value of the potential before and after the step of the target potential difference; subtracting the average value of the potential before and after the step of the target potential difference to obtain the target potential difference value.

5. A charging method data denoising system in a strong interference environment, characterized in that, comprising: a data acquisition module for acquiring transmission data transmitted by a transmission end and response data received by a receiving end; a position range determination module for determining a position range of a target potential difference in a main period according to a time difference between the transmission data and the start time of the response data; a period superposition module for dividing the response data into a plurality of main period segments according to the length of the main period of the transmission data waveform, superimposing the plurality of main period segments and calculating an average value to obtain superimposed main period segments; a step position determination module for performing step detection on the superimposed main period segments according to a step detection function, and determining the position of the step of the target potential difference in combination with the position range of the target potential difference; specifically comprising: Step 3.1: Based on the main period length of the transmitted waveform According to the formula Construct a step detection function ; Step 3.2: For the superimposed principal period segment Perform bidirectional extension, in Each side extends outwards using the average of its first and last 10 data points as the extension value. The length of the superimposed principal period segment is obtained by calculating the length of the extended segment. ; Step 3.3: The extended main period section of the superposition is convolved with the step detection function The absolute value of the result is taken, and the mean is removed, i.e. the processing of equation is performed, and the detection result is obtained :​​ Step 3.4: The length of the step function is deleted at both ends to ensure the correspondence between the step function and the main period segment after superposition, and the remaining part is called the step function , , , , , According to the range where the predetermined target potential difference is located, the target potential difference is determined by selecting the range The position of the target potential difference step is automatically determined according to the wave peak with the largest amplitude in the range; one or more step positions in the main period section after superposition are determined according to the wave peak of the step indication function; the wave peak with the largest amplitude in the position range is recorded as the position of the target potential difference step. a potential difference value calculation module for calculating the target potential difference value according to the potential data before and after the step of the target potential difference.

6. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor executes the program to realize the charging method data denoising method in a strong interference environment according to any one of claims 1-4.

7. A computer-readable storage medium having stored thereon a computer program, characterized in that The program is executed by the processor to realize the charging method data denoising method in a strong interference environment according to any one of claims 1-4.

Images