Geological radar data enhancement processing method and device, electronic equipment and storage medium

Abstract

The invention provides a geological radar data enhancement processing method and device, electronic equipment and a storage medium, and belongs to the field of geological radar data processing, and the method comprises the steps: obtaining geological radar data of a target region; selecting a window size according to depth information of a target body in the geological radar data; performing full-coverage subdivision on the geological radar data by adopting the selected window; obtaining a window center filtering value of each window; wherein the window center filtering value is obtained by the following steps of: dividing a plurality of sub-domains in a neighborhood of a current window center point in a preset sub-domain form; respectively calculating a signal mean value and a signal mean square error in each sub-domain; and selecting the sub-domain with the minimum signal mean square error from all the sub-domains of the current window, and taking a signal mean value corresponding to the sub-domain as a window center filtering value of the current window. The method can sharpen the boundary and enhance the boundary features of the anomalous body, thereby improving the recognition effect of the anomalous body in the geological radar image.

Description

Technical Field

[0001] The invention relates to a geological radar data enhancement processing method, device, electronic equipment and storage medium, and belongs to the field of geological radar data processing. Background Art

[0002] Geological radar (GPR) is a high-frequency narrow pulse electromagnetic wave excited by a research instrument, with a frequency of usually 10 6 ~10 9 Hz, a non-destructive detection device that uses the received reflected waveform to determine the distribution of underground media. The important theoretical basis for using this method for field detection is that there are obvious dielectric differences inside the detection object, such as various unfavorable geological bodies (cavities, broken zones, water-filled caves, etc.) There are obvious physical differences between them and the surrounding rocks. When electromagnetic waves encounter these interfaces during propagation, they will produce strong reflections. Correlation analysis of the received reflected waveform can qualitatively infer important information such as the spatial position, burial depth and size of the abnormal body inside the detection object.

[0003] However, when the physical properties of the unfavorable geological body (cavity, fracture zone, water-filled cave, etc.) are not much different from those of the surrounding rocks, the received reflected waveform is not very obvious, and errors will occur when inferring important information such as the spatial position, burial depth and size of the abnormal body inside the detection object. Therefore, it is necessary to enhance the received geological radar data through technical means. Summary of the invention

[0004] The purpose of the present invention is to overcome the deficiencies in the prior art and provide a geological radar data enhancement processing method, device, electronic device and storage medium to solve the problem that the reflected waveform received by the geological radar is not obvious enough.

[0005] To achieve the above object, the present invention is implemented by adopting the following technical solutions:

[0006] In the first aspect, the present invention provides a method for enhancing and processing geological radar data, comprising: acquiring geological radar data of a target area; selecting a window size based on depth information of a target body in the geological radar data; using the selected window size to perform full coverage segmentation on the geological radar data to obtain a plurality of windows; and using small subdomain filtering to obtain window center filter values ​​of each window.

[0007] Furthermore, the process of obtaining the window center filter value includes: dividing a plurality of subdomains in the neighborhood of the current window center point in a preset subdomain form; respectively calculating the signal mean and signal mean square error in each of the subdomains; selecting the subdomain with the smallest signal mean square error among all the subdomains of the current window, and using the signal mean corresponding to the subdomain as the window center filter value of the current window.

[0008] Furthermore, the window size is A square window; is the number of signal points on each side of the square window; An odd number not less than 3.

[0009] Furthermore, the subdomain is in the shape of a parallelogram, and the center point of the window to which the subdomain belongs falls inside the parallelogram subdomain or on the edge of the parallelogram.

[0010] Furthermore, the calculation formula of the signal mean is:

[0011] ;

[0012] in, N = [ ( M + 1 ) / 2 ] 2

[0013] In the formula, For the The signal mean of each sub-domain; For the In the subdomain The signal value at each point; is the number of signal points in each subdomain.

[0014] Furthermore, the calculation formula of the signal mean square error is:

[0015] σ i = ∑ j = 1 N [ Δ Ez i − Ez i ( j ) ] 2 N ;

[0016] In the formula, For the The mean square error of the signal in each sub-domain.

[0017] In a second aspect, the present invention provides a geological radar data enhancement processing device, comprising:

[0018] Data acquisition module: used to obtain geological radar data of the target area;

[0019] Window selection module: used for selecting the window size according to the depth information of the target body in the geological radar data;

[0020] Window segmentation module: used to perform full coverage segmentation on the geological radar data using the selected window size to obtain a plurality of windows;

[0021] Window processing module: used to obtain the window center filter value of each window by using small sub-domain filtering.

[0022] Furthermore, the window processing module includes:

[0023] Subdomain division module: used to divide multiple subdomains in the neighborhood of the current window center point in a preset subdomain form;

[0024] Subdomain calculation module: used to calculate the signal mean and signal mean square error in each subdomain respectively;

[0025] Decision module: used to select the sub-domain with the smallest signal mean square error among all sub-domains of the current window, and use the signal mean corresponding to the sub-domain as the window center filtering value of the current window.

[0026] In a third aspect, the present invention provides an electronic device comprising: a processor and a memory, wherein the memory stores computer-readable instructions, and when the computer-readable instructions are executed by the processor, the geological radar data enhancement processing method as described in the first aspect is implemented.

[0027] In a fourth aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the geological radar data enhancement processing method as described in the first aspect is implemented.

[0028] Compared with the prior art, the present invention has the following beneficial effects:

[0029] The geological radar data enhancement processing method provided by the present invention divides the data into windows, and uses small subdomain filtering to divide the window into different subdomains. When the electric field strength data in a subdomain changes smoothly, it indicates that it is in a stable field strength area, that is, the medium properties in the area reflected by it are relatively uniform; on the contrary, when a subdomain is in an area with a large change in electric field strength, it indicates that the medium properties in the area reflected by it have a large change, and the central part of the window is most likely to belong to a subdomain with a relatively gentle abnormal change. The filtering value of the center point of the window is obtained by the subdomain with the gentlest abnormal change in a window, which highlights the boundary characteristics of the abnormal change, sharpens the boundary characteristics of the abnormal body, and can improve the recognition effect of the abnormal body in the geological radar image. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] Figure 1 This is a flow chart of the geological radar data enhancement processing method provided in Example 1 of the present invention;

[0031] Figure 2 1 is a diagram of the subdomain partitioning method in Example 1 of the present invention; wherein a is the partitioning method when the window is 5×5, and b is the partitioning method when the window is 9×9;

[0032] Figure 3The original profile of the soft soil layer with a burial depth of 2m in Example 2 of the present invention and the profiles after processing with different window sizes; wherein, a is a profile formed by the original geological radar data of the soft soil layer with a burial depth of 2m, b is a profile formed by enhancing the geological radar data of the soft soil layer with a burial depth of 2m using a 5×5 window, and c is a profile formed by enhancing the geological radar data of the soft soil layer with a burial depth of 2m using a 9×9 window;

[0033] Figure 4 These are the original profiles of the soft soil layer with a burial depth of 10 m and the profiles after processing with different window sizes in Example 2 of the present invention; wherein, a is a profile view formed by the original geological radar data of the soft soil layer with a burial depth of 10 m, b is a profile view formed by enhancing the geological radar data of the soft soil layer with a burial depth of 10 m using a 5×5 window, and c is a profile view formed by enhancing the geological radar data of the soft soil layer with a burial depth of 10 m using a 9×9 window. DETAILED DESCRIPTION

[0034] The terms "including" and "having" and any variations thereof in the specification and claims of this application and the above-mentioned drawings are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices. The technical solutions in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings 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.

[0035] Example 1

[0036] The present invention provides a method for enhancing and processing geological radar data. Figure 1 As shown, the method comprises the following steps:

[0037] Step 1: Obtain geological radar data of the target area.

[0038] In geological exploration, ground penetrating radar data helps analyze underground structures and detect potential geological problems. Ground penetrating radar data is usually collected by scanners, probes, or arrays. These devices send high-frequency electromagnetic waves and receive the reflected signals. The collected data includes the amplitude and arrival time of the electromagnetic waves. The arrival time can be used to calculate the depth information of the target geological layer.

[0039] In field exploration, geological radar equipment is usually scanned by vehicle or manual means.

[0040] Step 2: Select the window size according to the depth information of the target body in the geological radar data.

[0041] Specifically, the size of the window in filtering has an obvious impact on the filtering result. The larger the window, the more beneficial it is to compress the stepped band. However, when it increases to a certain size, the enhancement effect will not be significantly improved. The ability of a small window to compress the stepped band is limited. In some cases, a small window may not achieve the compression effect, but its advantage is that the aliasing and distortion generated are small. Therefore, when selecting a window, both good filtering effect and aspects such as aliasing and distortion should be considered.

[0042] In a specific embodiment, the window size is a square window; where is the number of signal points on each side length of the square window; is an odd number not less than 3.

[0043] In a specific embodiment, the selection criterion for the window size is: for geological anomalies with an estimated target body burial depth less than 10 m, a 5×5 window is used; for geological anomalies with an estimated target body burial depth greater than 10 m and above, a 9×9 window is used.

[0044] Step 3: Perform full-coverage dissection on the ground penetrating radar data with the selected window size to obtain multiple windows.

[0045] In a specific embodiment, the ground penetrating radar data is subjected to full-coverage dissection in a sequentially seamless arrangement manner.

[0046] Step 4: Obtain the window center filtering value of each window by using small sub-domain filtering.

[0047] In a specific embodiment, the process of obtaining the window center filtering value includes:

[0048] Step 41: Dissect multiple sub-domains in the neighborhood of the current window center point in a preset sub-domain form.

[0049] In the prior art, common sub-domain dissection methods in the filtering process include the "rice" - shaped sub-domain dissection method, Tomita - Tsuji - shaped sub-domain dissection method, Nagao - Matsuyama - shaped sub-domain dissection method, symmetric - shaped sub-domain dissection method, "field" - shaped sub-domain dissection method, "cross" - shaped sub-domain dissection method, "×" - shaped sub-domain dissection method, improved "field" - shaped sub-domain dissection method, sliding - shaped sub-domain dissection method, and pentagon - hexagon sub-domain dissection method.

[0050] In a specific embodiment of the present invention, a parallelogram sub-domain form is used to dissect each window. As Figure 2 shown, the center point of the window to which the sub-domain belongs falls inside the parallelogram sub-domain or on the parallelogram edge, where Figure 2In figure a, the parallelogram subdomain partitioning method is shown when the window size is 5×5, and in figure b, the parallelogram subdomain partitioning method is shown when the window size is 9×9.

[0051] Step 42: Calculate the signal mean and signal mean square error in each of the sub-domains respectively.

[0052] In a specific embodiment, the calculation formula of the signal mean is:

[0053] ;

[0054] in, N = [ ( M + 1 ) / 2 ] 2 ;

[0055] In the formula, For the The signal mean of each sub-domain; For the In the subdomain The signal value at each point; is the number of signal points in each subdomain.

[0056] The calculation formula of signal mean square error is:

[0057] σ i = ∑ j = 1 N [ Δ Ez i − Ez i ( j ) ] 2 N ;

[0058] In the formula, For the The mean square error of the signal in each sub-domain.

[0059] Step 43: Select a sub-domain with the smallest signal mean square error among all sub-domains of the current window, and use the signal mean value corresponding to the sub-domain as the window center filtering value of the current window.

[0060] When the electric field strength data in a subdomain changes slowly, it indicates that it is in a stable field strength zone, that is, the medium properties reflected in the subdomain are relatively uniform. On the contrary, when a subdomain is in an area with large changes in electric field strength, the medium properties reflected in the area have large changes, and the center of the window is most likely to belong to the subdomain with relatively slow abnormal changes. The subdomain with the smallest mean square error is selected from all subdomains, and the mean of the subdomain is used as the filtering value of the center point of the window, which highlights the boundary characteristics of the abnormal changes and sharpens the boundary characteristics of the abnormal body, which is conducive to improving the recognition effect of the abnormal body in the geological radar image.

[0061] Example 2

[0062] This embodiment provides specific performances of the method of the present invention in the scenarios of surveying a soft soil layer at a burial depth of 2 m and a soft soil layer at a burial depth of 10 m.

[0063] In actual survey practice, geological radar equipment is first used for scanning to obtain geological radar data of soft soil layers with a burial depth of 2m and soft soil layers with a burial depth of 10m respectively.

[0064] Specifically, Figure 3 The original profile of the soft soil layer with a burial depth of 2m and the profile after processing with different window sizes, where a is the profile formed by the original geological radar data of the soft soil layer with a burial depth of 2m.

[0065] Figure 4 The original profile of the soft soil layer with a burial depth of 10m and the profile after processing with different window sizes, where a is the profile formed by the original geological radar data of the soft soil layer with a burial depth of 10m.

[0066] Then, the original geological radar data of the soft soil layer with a burial depth of 2m and the soft soil layer with a burial depth of 10m were segmented using a 5×5 window and a 9×9 window, respectively, and the window center filtering value of each segmentation window was obtained.

[0067] Specifically, Figure 3 b shows the profile of the enhanced geological radar data of the soft soil layer with a burial depth of 2m using a 5×5 window, while Figure 3 c shows the profile of the enhanced geological radar data of the soft soil layer at a depth of 2 m after being processed using a 9×9 window.

[0068] Figure 4 b shows the profile of the enhanced processing of the soft soil layer geological radar data with a burial depth of 10m using a 5×5 window, while Figure 4 c shows the profile of the enhanced geological radar data of the soft soil layer at a depth of 10 m after being processed using a 9×9 window.

[0069] By comparing a, b and c in the two figures respectively, it can be seen that in the geological radar image after small subdomain filtering, the reflected wave is more obvious, the shape is more continuous and complete, the boundary of the anomaly is sharpened, the boundary characteristics of the anomaly are enhanced, and the recognition effect of the anomaly is also improved.

[0070] By comparing b and c in the two figures, it can be seen that the larger the window is, the more obvious the enhancement effect on the reflected wave is, and the better the recognition effect of the abnormal body is.

[0071] from Figure 3 and Figure 4 It can also be seen that when the soft soil layer is shallow, the 5×5 small window can also enhance the reflected wave well. Although it is not as obvious as the 9×9 large window enhancement effect, it can also better reflect the boundary characteristics of the anomaly. When the soft soil layer is deep, the 9×9 large window can better reflect the boundary characteristics of the anomaly compared to the 5×5 small window.

[0072] In summary, the present invention proposes to use small subdomain filtering technology to enhance the geological radar data. This technology can sharpen the boundaries and enhance the boundary features of the anomalies, thereby improving the recognition effect of the anomalies in the geological radar images.

[0073] Example 3

[0074] This embodiment provides a geological radar data enhancement processing device, including:

[0075] Data acquisition module: used to obtain geological radar data of the target area;

[0076] Window selection module: used for selecting the window size according to the depth information of the target body in the geological radar data;

[0077] Window segmentation module: used to perform full coverage segmentation on the geological radar data using the selected window size to obtain a plurality of windows;

[0078] Window processing module: used to obtain the window center filter value of each window by using small sub-domain filtering;

[0079] Example 4

[0080] Based on Example 3, this example provides a specific composition of a window processing module.

[0081] The window processing module includes:

[0082] Subdomain division module: used to divide multiple subdomains in the neighborhood of the current window center point in a preset subdomain form;

[0083] Subdomain calculation module: used to calculate the signal mean and signal mean square error in each subdomain respectively;

[0084] Decision module: used to select the sub-domain with the smallest signal mean square error among all sub-domains of the current window, and use the signal mean corresponding to the sub-domain as the window center filtering value of the current window.

[0085] Example 5

[0086] This embodiment provides an electronic device, including: a processor and a memory, wherein the memory stores computer-readable instructions, and when the computer-readable instructions are executed by the processor, the geological radar data enhancement processing method as described in Example 1 is implemented.

[0087] Example 6

[0088] This embodiment provides a computer storage medium, on which a computer program is stored. When the computer program is executed by a processor, the geological radar data enhancement processing method as described in Example 1 is implemented.

[0089] The storage medium may include, for example, a storage component of a tablet computer, a hard disk of a computer, a read-only memory (ROM), an erasable programmable read-only memory (EPROM), a portable compact disk read-only memory (CD-ROM), a USB memory, or any combination of the above storage media. The computer storage medium may be any combination of one or more computer storage media.

[0090] The above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A geological radar data enhancement processing method, characterized in that: include: Acquire geological radar data of the target area; Selecting a window size according to the depth information of the target body in the geological radar data; Using the selected window size to perform full coverage segmentation on the geological radar data to obtain a plurality of windows; Small subdomain filtering is used to obtain the window center filter value of each window.

2. The geological radar data enhancement processing method according to claim 1, characterized in that: The process of obtaining the window center filter value includes: Divide the neighborhood of the current window center point into multiple subdomains in a preset subdomain form; Calculating the signal mean and signal mean square error in each of the sub-domains respectively; A subdomain with the smallest signal mean square error is selected from all subdomains of the current window, and the signal mean corresponding to the subdomain is used as the window center filtering value of the current window.

3. The geological radar data enhancement processing method according to claim 1, characterized in that: The window size is A square window; is the number of signal points on each side of the square window; An odd number not less than 3.

4. The geological radar data enhancement processing method according to claim 2, characterized in that: The subdomain is in the shape of a parallelogram, and the center point of the window to which the subdomain belongs falls inside the parallelogram subdomain or on the edge of the parallelogram.

5. The geological radar data enhancement processing method according to claim 2, characterized in that: The calculation formula of the signal mean is: ； in, ; In the formula, For the The signal mean of each sub-domain; For the In the subdomain The signal value at each point; is the number of signal points in each subdomain.

6. The geological radar data enhancement processing method according to claim 2, characterized in that: The calculation formula of the signal mean square error is: ； In the formula, For the The mean square error of the signal in each sub-domain.

7. A geological radar data enhancement processing device, characterized in that: include: Data acquisition module: used to obtain geological radar data of the target area; Window selection module: used for selecting the window size according to the depth information of the target body in the geological radar data; Window segmentation module: used to perform full coverage segmentation on the geological radar data using the selected window size to obtain a plurality of windows; Window processing module: Use small subdomain filtering to obtain the window center filter value of each window.

8. The geological radar data enhancement processing device according to claim 7, characterized in that: The window processing module includes: Subdomain division module: used to divide multiple subdomains in the neighborhood of the current window center point in a preset subdomain form; Subdomain calculation module: used to calculate the signal mean and signal mean square error in each subdomain respectively; Decision module: used to select the sub-domain with the smallest signal mean square error among all sub-domains of the current window, and use the signal mean corresponding to the sub-domain as the window center filtering value of the current window.

9. An electronic device, characterized in that: include: A processor and a memory, wherein the memory stores computer-readable instructions, and when the computer-readable instructions are executed by the processor, the geological radar data enhancement processing method according to any one of claims 1 to 6 is implemented.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by the processor, the geological radar data enhancement processing method according to any one of claims 1 to 6 is implemented.

Images