A ground penetrating radar data batch simulation method, system, device and medium

Abstract

The present application belongs to the field of ground penetrating radar forward simulation of underground empty area, and particularly relates to a ground penetrating radar data batch simulation method, system, device and medium. The present application obtains stratum parameters where an actual engineering underground empty area disease body is located, and constructs an underground empty area random generation algorithm. Different shapes, sizes and buried depths of underground empty area models are batch generated. A model storage module is established to store the batch generated underground empty area random shape models. The stored underground empty area random shape models are batch read and combined with model material parameters for corresponding ground penetrating radar forward simulation to obtain radar data of the corresponding underground empty area models. The present application can batch simulate ground penetrating radar data of underground empty area, can randomly batch generate underground empty area radar data of different sizes, shapes and buried depths, and can provide data basis for accurate identification of underground empty area and efficient data support for diagnosis and treatment of underground concealed diseases.

Description

Technical Field

[0001] This invention belongs to the field of forward modeling of ground-penetrating radar in underground airspace, and specifically relates to a method, system, equipment and medium for batch simulation of ground-penetrating radar data. Background Technology

[0002] In underground void identification, obtaining a large amount of underground void data is crucial for the identification, diagnosis, and maintenance of underground voids. However, the amount of actual detected void data is limited and insufficient to meet the needs of identification and diagnosis. To compensate for this data scarcity, numerical simulation has become an important supplementary method, capable of effectively simulating large amounts of data. However, current numerical simulation methods for ground-penetrating radar (GPR) data are quite cumbersome, requiring the design of computational models for different situations, and the computation time required to generate large amounts of void data is long. This approach is inefficient for generating batches of data required for underground void identification. Therefore, there is an urgent need for a method that can efficiently generate underground void data in batches to fill the gaps in actual detection data, improve the efficiency of underground void identification, provide strong support for batch simulation of underground void GPR data, and provide a reference for the prevention and treatment of hidden underground voids. Summary of the Invention

[0003] This invention proposes a method, system, equipment, and medium for batch simulation of ground-penetrating radar data to solve the problem of low efficiency in underground airspace identification in existing technologies.

[0004] To achieve the above objectives, the present invention proposes the following technical solution:

[0005] A method for batch simulation of ground-penetrating radar data includes the following steps:

[0006] Step 1: Obtain the stratigraphic parameters of the target underground void area;

[0007] Step 2: Based on the geological parameters, construct a random generation algorithm for underground voids to generate underground void models of different shapes, sizes and burial depths in batches;

[0008] Step 3: Store the underground void model;

[0009] Step 4: Obtain the model material parameters of the target underground void area;

[0010] Step 5: Read the stored underground void model and model material parameters, perform ground-penetrating radar forward modeling, and obtain radar data of the target underground void.

[0011] Preferably, the formation parameters in step 1 include formation type, formation thickness, formation dielectric constant, and formation conductivity.

[0012] Preferably, the algorithm for randomly generating underground voids in step 2 is as follows:

[0013] Step 21, construct the basic model function matrix;

[0014] Step 22: Establish an empty area model generation algorithm, which generates an empty area model based on the matrix function in step 21.

[0015] Preferably, the process of constructing the basic model function matrix in step 21 is as follows:

[0016] Determine the first corner point of the model, and set it as (x1, y1);

[0017] Define the length and width of the model as x. a and y b ;

[0018] Use the `plot` function to draw the base model of the required size. The matrix functions needed to build the base model in the Python environment are:

[0019] plt.plot([x1, x1, x1+x) a ,x1+x a ],[y1,y1+y b y1+y b ,y1]).

[0020] Preferably, the empty area model generation algorithm in step 22 includes a regular empty area model generation method and an irregular empty area model generation method.

[0021] Preferably, the method for generating the regular-shaped empty area model specifically includes:

[0022]

[0023] f represents the function that generates a regular-shaped empty region model, where p1, p2, ... p n These are the parameters required for model generation;

[0024] (x1, y1), (x2, y2),…, (x n y n ): Represents the coordinates of each point in the generated regular-shaped empty area model, i.e., the boundary points of the model;

[0025] x k and y k : Represents the rectangular coordinates of the k-th point, where k = 0, 1, 2, ..., n;

[0026] x c and y c : Represents the rectangular coordinates of the center point of the figure;

[0027] γ ak and γ bk Let represent the major axis diameter and minor axis diameter of the circle containing the k-th point, respectively;

[0028] θ k Let represent the angle between the k-th point and the center point of the figure, i.e., the polar angle, where k = 0, 1, 2, ..., n.

[0029] Preferably, the method for generating the irregular-shaped empty region model specifically includes:

[0030] A basic graphic function is established based on the regular shape empty area model generation method. Based on the basic graphic function, the shape size range and irregularity value are defined, and the parameter variable θ is generated based on a random generator. k and γ k This yielded empty areas of different sizes and burial depths;

[0031] The np.linspace function generates a regular shape of random size. The random parameters generated by the np.random.normal function are then superimposed on the shape generated by the np.linspace function, making the basic regular shape an irregular shape.

[0032] The edge angles are cyclically superimposed, and the variable parameters are modified by a Gaussian distribution function to continuously optimize the irregularity of the graphic edges.

[0033] Delete the points outside the shape, draw the shape and fill the points inside the shape to complete the generation of the irregular empty area model.

[0034] A ground-penetrating radar (GPR) data batch simulation system for implementing the aforementioned method for batch simulation of GPR data includes the following modules:

[0035] The system includes a stratigraphic parameter acquisition module, a subsurface void model building module, a model storage module, an integrated software module, a model forward modeling module, and a subsurface void ground-penetrating radar data storage module.

[0036] The formation parameter acquisition module is used to acquire the formation parameters of the target underground void area.

[0037] The underground void model building module is used to construct an underground void random generation algorithm based on the stratigraphic parameters of the target underground void, and to generate underground void models of different shapes, sizes and burial depths in batches.

[0038] The model storage module is used to store underground void models.

[0039] The integrated software module is used to establish the input program for model material parameters and to obtain model material parameters;

[0040] The model forward modeling module is used to read the underground void model and model material parameters, and perform ground-penetrating radar forward modeling to obtain radar data of the target underground void model;

[0041] The underground airspace ground-penetrating radar data storage module is used to store radar data of the target underground airspace model.

[0042] An electronic device, including a memory and a processor;

[0043] Memory, used to store computer programs;

[0044] A processor is used to execute the computer program, which, when executed by the processor, implements the steps of the method for batch simulation of ground-penetrating radar data.

[0045] A computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the method for batch simulation of ground-penetrating radar data.

[0046] The advantages of this invention are:

[0047] This invention proposes a batch simulation method for ground-penetrating radar (GPR) data. Based on stratigraphic parameters, a random generation algorithm for underground voids is constructed, which can randomly generate underground void models of different sizes, shapes, and burial depths in batches. After storage, forward modeling is performed to obtain radar data of the target underground void model. This invention realizes batch simulation of GPR data for underground voids, providing a data foundation for accurate identification of underground voids and providing efficient data support for the diagnosis and treatment of hidden underground diseases. Attached Figure Description

[0048] The accompanying drawings, which form part of this specification, are used to provide a further understanding of the invention. The illustrative embodiments of the invention and their descriptions are used to explain the invention and do not constitute an undue limitation of the invention. In the drawings:

[0049] Figure 1 This is a schematic diagram of a method for batch simulation of ground-penetrating radar data.

[0050] Figure 2 This is a schematic diagram of an electronic device for simulating batch data from ground-penetrating radar. Detailed Implementation

[0051] The present invention will now be described in detail with reference to the accompanying drawings and embodiments. It should be noted that, unless otherwise specified, the embodiments and features described herein can be combined with each other.

[0052] The following detailed description is exemplary and intended to provide further detailed explanation of the invention. Unless otherwise specified, all technical terms used in this invention have the same meaning as commonly understood by one of ordinary skill in the art. The terminology used in this invention is for describing particular embodiments only and is not intended to limit the scope of exemplary embodiments according to the invention.

[0053] Example 1:

[0054] Please see Figure 1 As shown, the present invention provides a method for batch simulation of ground penetrating radar data, comprising:

[0055] Step 1: Obtain the geological parameters of the underground voids in the actual project through field exploration, for subsequent numerical simulation; the geological parameters include geological type, geological thickness, geological dielectric constant and geological conductivity.

[0056] Step 2: Based on the obtained stratigraphic parameters, construct a random generation algorithm for underground voids, establish a multi-parameter underground void simulation module, and generate underground void models of different shapes, sizes, and burial depths in batches.

[0057] The multi-parameter underground void simulation module includes a regular void construction module and an irregular void construction module;

[0058] The specific construction steps include:

[0059] Step 21, Model Size and Stratigraphic Extent Construction: Use the plot function to generate the base model, thereby defining the size and extent of the model;

[0060] The specific implementation steps are as follows:

[0061] Determine the first corner point of the model, and set it as (x1, y1);

[0062] Define the length and width of the model as x. a and y b ;

[0063] Use the `plot` function to draw the base model of the required size, including the model extent and formation thickness. The matrix functions needed to build the base model in the Python environment are shown below:

[0064] plt.plot([x1, x1, x1+x) a ,x1+x a ,],[y1,y1+y b y1+y b ,y1])

[0065] Step 22: Random generation algorithms were established for both regular and irregular empty area models, and empty area models were generated based on the matrix function in Step 21.

[0066] The regular-shaped void model construction method can generate void models of different sizes and burial depths with the required shapes, including circles, ellipses, triangles, rectangles, trapezoids and polygons.

[0067] The method for constructing the regular-shaped empty region model is as follows:

[0068]

[0069] f represents the function that generates a regular-shaped empty region model, where p1, p2, ... p n These are the parameters required for model generation.

[0070] (x1, y1), (x2, y2),…, (x n y n ): Represents the coordinates of each point in the generated regular-shaped empty area model, i.e., the boundary points of the model. n has no range of values.

[0071] x k and y k : Represents the rectangular coordinates of the k-th point, where k = 0, 1, 2, ..., n.

[0072] x c and y c : Represents the rectangular coordinates of the center point of the figure.

[0073] γ ak and γ bk : Represents the major and minor axis diameters of the circle containing the k-th point.

[0074] θ k : represents the angle between the k-th point and the center point of the figure, i.e., the polar angle, where k = 0, 1, 2, ..., n.

[0075] Furthermore, the method for constructing the irregular-shaped empty region model includes the following specific steps:

[0076] A method for constructing regular-shaped empty regions is used to establish basic graphics functions;

[0077] Define the shape size range and the degree of irregularity;

[0078] The parameter variables are generated based on a random generator, and the parameter variables include θ. k γ ak and γ bk By adjusting the parameter variables, voids of different sizes and burial depths can be generated;

[0079] Use np.linspace to generate a graph of random size according to a set rule;

[0080] By overlaying the random parameters generated by np.random.normal onto the graph generated by np.linspace, the basic regular graph becomes an irregular shape;

[0081] The edge angles are cyclically superimposed, and the parameter variables are modified by a Gaussian distribution function to continuously optimize the irregularity of the graphic edges.

[0082] Calculate the distance from each point in the graphic to the center point of the graphic, keep the points inside the graphic, and delete the points outside the graphic;

[0083] Draw a shape and fill the points inside the shape.

[0084] Step 3: Establish a model storage module to store the generated batch of random shape models of underground voids. The random shape model storage module uses h5py to store the constructed random shape model data of underground voids as HDF5 files.

[0085] Step 4: Create an integrated software package for gprMax and Python, establish a Python environment for gprMax, and then run gprMax via a Python script. Also, create a program to input the model's material parameters.

[0086] Step 5: Using the established integrated software, batch read the stored random shape models and material parameters of underground voids to perform ground-penetrating radar forward modeling, obtaining the radar data of the corresponding underground void models. The integrated software reads the stored underground void models using the following program syntax:

[0087] #geometry_objects_read:f1f2f3file1file2

[0088] Where f1f2f3 are the coordinates (x, y, z) of the lower left corner of the field where the geometry array should be placed, file1 is the path and filename of the HDF5 file containing the integer array defining the geometry, and file2 is the path and filename of the text file containing the commands.

[0089] Step 6: Establish a ground-penetrating radar data storage module for underground airspace, perform a forward modeling process to obtain batches of underground airspace radar data, and store the batches of underground airspace radar data in a database.

[0090] Example 2:

[0091] This invention provides a system for batch simulation of ground-penetrating radar data, specifically comprising:

[0092] The system includes a geological parameter acquisition module, a subsurface void model building module, a model storage module, an integrated software module, a model forward modeling module, and a subsurface void ground-penetrating radar data storage module.

[0093] The geological parameter acquisition module is used to implement step 1 above to acquire the geological parameters of the underground void disease body in the actual project.

[0094] The underground void model building module is used to implement step 2 above. Based on the obtained stratum parameters, it constructs an underground void random generation algorithm, establishes a multi-parameter underground void simulation module, and generates underground void models of different shapes, sizes and burial depths in batches.

[0095] The model storage module is used to implement step 3 above and store the generated batch of random shape models of underground voids.

[0096] The integrated software module is used to implement step 4 above and establish the input program for model material parameters;

[0097] The model forward modeling module is used to implement step 5 above, batch read the stored random shape model and model material parameters of the underground void area, perform ground-penetrating radar forward modeling of the corresponding model, and obtain the radar data of the corresponding underground void area model.

[0098] The underground airspace ground-penetrating radar data storage module is used to implement step 6 above, storing the batch of underground airspace radar data obtained from the model forward modeling into the underground airspace ground-penetrating radar data storage module.

[0099] Example 3:

[0100] Please see Figure 2 As shown, the present invention also provides an electronic device 100; the electronic device 100 includes a memory 101, at least one processor 102, a computer program 103 stored in the memory 101 and executable on the at least one processor 102, and at least one communication bus 104.

[0101] The memory 101 can be used to store the computer program 103. The processor 102 implements the steps of the method for batch simulation of ground-penetrating radar data as described in Embodiment 1 by running or executing the computer program stored in the memory 101 and calling the data stored in the memory 101. The memory 101 may mainly include a program storage area and a data storage area. The program storage area may store the operating system, at least one application program required for a function (such as sound playback function, image playback function, etc.), etc.; the data storage area may store data created according to the use of the electronic device 100 (such as audio data), etc. In addition, the memory 101 may include non-volatile memory, such as hard disk, memory, plug-in hard disk, smart media card (SMC), secure digital (SD) card, flash card, at least one disk storage device, flash memory device, or other non-volatile solid-state storage device.

[0102] The at least one processor 102 may be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The processor 102 may be a microprocessor or any conventional processor. The processor 102 is the control center of the electronic device 100, connecting various parts of the electronic device 100 via various interfaces and lines.

[0103] The memory 101 in the electronic device 100 stores multiple instructions to implement a method for batch simulation of ground-penetrating radar data.

[0104] Example 4:

[0105] If the modules / units integrated in the electronic device 100 are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, all or part of the processes in the methods of the above embodiments can also be implemented by a computer program instructing related hardware. The computer program can be stored in a computer-readable storage medium, and when executed by a processor, it can implement the steps of the various method embodiments described above. The computer program includes computer program code, which can be in the form of source code, object code, executable files, or certain intermediate forms. The computer-readable medium can include: any entity or device capable of carrying the computer program code, a recording medium, a USB flash drive, a portable hard drive, a magnetic disk, an optical disk, a computer memory, and a read-only memory (ROM).

[0106] As is known from common technical knowledge, this invention can be implemented through other embodiments that do not depart from its spirit or essential characteristics. Therefore, the disclosed embodiments described above are merely illustrative in all respects and are not the only ones. All modifications within the scope of this invention or its equivalents are included in this invention.

[0107] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0108] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0109] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0110] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0111] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that modifications or equivalent substitutions can still be made to the specific implementation of the present invention. Any modifications or equivalent substitutions that do not depart from the spirit and scope of the present invention should be covered within the scope of protection of the claims of the present invention.

Claims

1. A method for batch simulation of ground-penetrating radar data, characterized in that, Includes the following steps: Step 1: Obtain the stratigraphic parameters of the target underground void area; Step 2: Based on the geological parameters, construct a random generation algorithm for underground voids to generate underground void models of different shapes, sizes and burial depths in batches; Step 3: Store the underground void model; Step 4: Obtain the model material parameters of the target underground void area; Step 5: Read the stored underground void model and model material parameters, perform ground-penetrating radar forward modeling, and obtain radar data of the target underground void. The specific steps of constructing the random generation algorithm for underground voids in step 2 are as follows: Step 21, construct the basic model matrix function; Step 22: Establish an empty area model generation algorithm, and generate an empty area model based on the matrix function in Step 21; The empty area model generation algorithm in step 22 includes a regular empty area model generation method and an irregular empty area model generation method; The method for generating the irregular-shaped void region model is as follows: A basic graphic function is established based on the regular shape void model generation method. Based on the basic graphic function, the shape size range and irregularity value are defined, and parameter variables are generated based on a random generator to obtain voids of different sizes and burial depths. The np.linspace function generates a regular shape of random size. The random parameters generated by the np.random.normal function are then superimposed on the shape generated by the np.linspace function to make the basic regular shape become an irregular shape. The edge angles are cyclically superimposed, and the variable parameters are modified through a Gaussian distribution function to continuously optimize the irregularity of the graphic edges; Delete the points outside the shape, draw the shape and fill the points inside the shape to complete the generation of the irregular empty area model.

2. The method for batch simulation of ground-penetrating radar data as described in claim 1, characterized in that, The formation parameters in step 1 include formation type, formation thickness, formation dielectric constant, and formation conductivity.

3. The method for batch simulation of ground-penetrating radar data as described in claim 1, characterized in that, The process of constructing the basic model matrix function in step 21 is as follows: Determine the first corner point of the model, and set it as (x1, y1); Define the length and width of the model as x. a and y b ; Use the `plot` function to draw the base model of the required size. The matrix functions needed to build the base model in the Python environment are: plt.plot([x1, x1 , x1+ x a , x1+ x a ], [y1, y1+ y b , y1+ y b , y1])。 4. The method for batch simulation of ground-penetrating radar data as described in claim 1, characterized in that, The method for generating the regular-shaped empty area model is as follows: f represents the function that generates a regular-shaped empty region model, where p1, p2, ... p n These are the parameters required for model generation; (x1, y1), (x2, y2), ⋯, (x n y n ): Represents the coordinates of each point in the generated regular-shaped empty area model, i.e., the boundary points of the model; x k and y k : Represents the rectangular coordinates of the k-th point, where k = 0, 1, 2, ..., n; x c and y c : Represents the rectangular coordinates of the center point of the figure; γ ak and γ bk Let these represent the major axis diameter and minor axis diameter of the ellipse containing the k-th point, respectively. θ k Let represent the angle between the k-th point and the center point of the figure, i.e., the polar angle, where k = 0, 1, 2, ..., n.

5. A ground-penetrating radar data batch simulation system for implementing the batch simulation method for ground-penetrating radar data according to any one of claims 1-4, characterized in that, Includes the following modules: The system includes a stratigraphic parameter acquisition module, a subsurface void model building module, a model storage module, an integrated software module, a model forward modeling module, and a subsurface void ground-penetrating radar data storage module. The formation parameter acquisition module is used to acquire the formation parameters of the target underground void area; The underground void model building module is used to construct an underground void random generation algorithm based on the stratigraphic parameters of the target underground void, and to generate underground void models of different shapes, sizes and burial depths in batches. The model storage module is used to store underground void models; The integrated software module is used to establish the input program for model material parameters and to obtain model material parameters; The model forward modeling module is used to read the underground void model and model material parameters, and perform ground-penetrating radar forward modeling to obtain radar data of the target underground void model; The underground airspace ground-penetrating radar data storage module is used to store radar data of the target underground airspace model.

6. An electronic device, characterized in that, Including memory and processor; Memory, used to store computer programs; A processor is configured to execute the computer program, which, when executed by the processor, implements the steps of a batch simulation method for ground-penetrating radar data as described in any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that, The storage medium stores a computer program, which, when executed by a processor, implements the steps of a batch simulation method for ground-penetrating radar data as described in any one of claims 1 to 4.

Images