Method and device for inverting and reconstructing soil layer moisture content, electronic equipment and storage medium

Abstract

The invention provides a method and device for inverting the moisture content of a reconstructed soil layer, electronic equipment and a storage medium, and the method comprises the steps: detecting the structure of the reconstructed soil layer through a ground penetrating radar system, and obtaining the electromagnetic wave propagation data of each layer of the reconstructed soil layer; detecting the reconstructed soil layer provided with the reflector based on a ground penetrating radar system to obtain two-way travel time data, obtaining a relative dielectric constant of the dielectric layer based on the two-way travel time data and the burying depth of the reflector, and obtaining an empirical formula based on the relative dielectric constant; and correcting the empirical formula based on the actual moisture content of the reconstructed soil layer to obtain a moisture content model, and performing inversion based on the moisture content model to obtain the moisture content of the reconstructed soil layer. By combining the ground penetrating radar system and the corrected empirical formula, the method can significantly improve the accuracy of soil layer moisture content inversion, and is especially suitable for land reclamation and ecological restoration engineering in arid and semi-arid strip mine areas.

Description

Technical Field

[0001] The present invention belongs to the field of soil science and environmental monitoring technology, and particularly relates to a method, device, electronic equipment and storage medium for inverting and reconstructing soil layer moisture content. Background Art

[0002] The increase in mining activities, especially open-pit mining, has severely damaged the surface ecosystem, leading to increasingly severe water shortages and land degradation. Land reclamation and ecological restoration have become key components of ecological restoration in mining areas. Scientifically reconstructing the soil structure is a crucial prerequisite for ecological restoration. Soil moisture content, a key parameter for measuring soil moisture status, has a significant impact on vegetation growth and ecosystem recovery.

[0003] Existing soil moisture measurement methods, such as oven drying, time domain reflectometry (TDR), neutron analysis, and remote sensing, can provide soil moisture data to a certain extent. However, these methods suffer from shortcomings such as low measurement efficiency, limited measurement point distribution, insufficient monitoring depth, and potential environmental hazards. Ground penetrating radar (GPR), as a non-invasive, rapid, and continuous underground exploration technology, has broad application prospects, but its accuracy in inverting deep soil moisture still needs to be improved.

[0004] Therefore, the present application contemplates a method for detecting soil moisture content that combines GPR technology with a modified empirical formula. Summary of the Invention

[0005] In response to the problems existing in the prior art, the present invention provides a method, device, electronic equipment and storage medium for inverting and reconstructing the moisture content of a soil layer, which can efficiently and accurately detect the moisture content of a soil layer.

[0006] In a first aspect, an embodiment of the present invention provides a method for inverting and reconstructing soil moisture content, comprising:

[0007] Use a ground penetrating radar system to detect the structure of the reconstructed soil layer and obtain electromagnetic wave propagation data of each layer of the reconstructed soil layer;

[0008] Detecting the reconstructed soil layer provided with a reflector using a ground penetrating radar system to obtain two-way travel time data, obtaining a relative dielectric constant of the dielectric layer based on the two-way travel time data and the buried depth of the reflector, and obtaining an empirical formula based on the relative dielectric constant;

[0009] Based on the actual moisture content of the reconstructed soil layer, the empirical formula is modified to obtain a moisture content model, and the moisture content of the reconstructed soil layer is obtained by inversion based on the moisture content model.

[0010] Furthermore, the reconstructed soil layer provided with reflectors includes: a plurality of reflectors of different sizes sequentially spaced apart along the direction of emission of the radar electromagnetic waves.

[0011] Furthermore, the ground penetrating radar system preprocesses the detection data obtained by detection, and the preprocessing includes:

[0012] performing background denoising on the detection data to obtain denoised detection information;

[0013] Performing one-dimensional FIR filtering on the denoised detection information to obtain frequency-domain filtered detection information;

[0014] Gain control is performed on the frequency domain filtered detection information to obtain optimized detection data.

[0015] Furthermore, the empirical formula includes a Topp empirical formula, which includes:

[0016] θ v =-5.3×10 -2 +2.292×10 -2 ε r -5.5×10 -4 ε r 2 +4.3×10 -6 ε r 3 ;

[0017] In the formula, the θ v is the moisture content; ε r is the relative dielectric constant.

[0018] Furthermore, the relative dielectric constant is:

[0019]

[0020] Wherein, c is the speed of light; v is the propagation speed of radar electromagnetic waves, and the propagation speed of radar electromagnetic waves is calculated based on the two-way travel time data and the buried depth of the reflector.

[0021] Furthermore, the empirical formula is corrected using a polynomial regression analysis correction method, including:

[0022] Calculate the preliminary moisture content based on the empirical formula;

[0023] Based on the preliminary moisture content and the measured moisture content, the moisture content model was obtained by using the polynomial regression analysis correction method.

[0024] Furthermore, the moisture content model includes:

[0025] θ1=-24.80918+11.67223ε r -1.19808ε r 2 +0.04163ε r 3 ;

[0026] Where θ1 is the moisture content; ε r is the relative dielectric constant.

[0027] In a second aspect, an embodiment of the present disclosure provides a device for inverting and reconstructing soil moisture content, the device comprising:

[0028] The primary acquisition unit is configured to detect the structure of the reconstructed soil layer using a ground penetrating radar system to obtain electromagnetic wave propagation data of each layer of the reconstructed soil layer;

[0029] a secondary acquisition unit configured to detect the reconstructed soil layer provided with a reflector using a ground penetrating radar system to obtain two-way travel time data, obtain a relative dielectric constant of the dielectric layer based on the two-way travel time data and the burial depth of the reflector, and obtain an empirical formula based on the relative dielectric constant;

[0030] The output unit is configured to correct the empirical formula based on the actual moisture content of the reconstructed soil layer to obtain a moisture content model, and invert the moisture content of the reconstructed soil layer based on the moisture content model.

[0031] In a third aspect, an embodiment of the present disclosure provides an electronic device, the electronic device comprising:

[0032] at least one processor; and,

[0033] a memory communicatively connected to the at least one processor; wherein,

[0034] The memory stores instructions that can be executed by the at least one processor. The instructions are executed by the at least one processor to enable the at least one processor to perform the above-mentioned method for inverting and reconstructing the moisture content of the soil layer.

[0035] In a fourth aspect, an embodiment of the present disclosure provides a non-transitory computer-readable storage medium, which stores computer instructions for causing the computer to execute the above-mentioned method for inverting and reconstructing the moisture content of the soil layer.

[0036] Other optional features and technical effects of the embodiments of the present invention are partially described below, and partially can be understood by reading this document.

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

[0038] The present invention provides a method for inverting and reconstructing the moisture content of a soil layer, comprising: using a ground-penetrating radar system to detect the structure of a reconstructed soil layer to obtain electromagnetic wave propagation data for each layer of the reconstructed soil layer; detecting the reconstructed soil layer with a reflector disposed therein using the ground-penetrating radar system to obtain two-way travel time data; obtaining the relative dielectric constant of the dielectric layer based on the two-way travel time data and the burial depth of the reflector; and obtaining an empirical formula based on the relative dielectric constant; correcting the empirical formula based on the actual moisture content of the reconstructed soil layer to obtain a moisture content model; and inverting the moisture content of the reconstructed soil layer based on the moisture content model. By combining the ground-penetrating radar system with the corrected empirical formula, this method can significantly improve the accuracy of soil moisture inversion and is particularly suitable for land reclamation and ecological restoration projects in arid and semi-arid open-pit mining areas. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] Figure 1 A schematic flow chart of a method for inverting and reconstructing soil moisture content according to an embodiment of the present disclosure is shown;

[0040] Figure 2 A schematic diagram of the process of inverting and reconstructing the moisture content of a soil layer according to an embodiment of the present disclosure is shown;

[0041] Figure 3 A schematic diagram of the wave velocity calibration principle of the ground penetrating radar system GPR according to an embodiment of the present disclosure is shown;

[0042] Figure 4 A schematic diagram of measuring point layout and steel pipe burial according to an embodiment of the present disclosure is shown;

[0043] Figure 5 Shows an on-site detection diagram of a three-layer sponge structure area according to an embodiment of the present disclosure;

[0044] Figure 6 A radar detection diagram at a depth of 3 m according to an embodiment of the present disclosure is shown;

[0045] Figure 7 A radar detection image at a depth of 1.5 m according to an embodiment of the present disclosure is shown;

[0046] Figure 8 A radar detection diagram of a 3m mixed zone according to an embodiment of the present disclosure is shown;

[0047] Figure 9 A comparison curve of the TOPP formula calculated value and the measured value according to the embodiment of the present disclosure is shown;

[0048] Figure 10 A fitting curve diagram of relative dielectric constant and volume water content of an embodiment of the present disclosure is shown.

[0049] Figure 11A diagram of a device for inverting and reconstructing soil layer moisture content according to an embodiment of the present disclosure is shown. DETAILED DESCRIPTION

[0050] The embodiments of the present disclosure are described in detail below with reference to the accompanying drawings.

[0051] The following describes the embodiments of the present disclosure through specific examples, and those skilled in the art can easily understand other advantages and effects of the present disclosure from the contents disclosed in this specification. Obviously, the described embodiments are only a part of the embodiments of the present disclosure, rather than all of the embodiments. The present disclosure can also be implemented or applied through other different specific embodiments, and the details in this specification can also be modified or changed in various ways based on different viewpoints and applications without departing from the spirit of the present disclosure. It should be noted that, in the absence of conflict, the following embodiments and features in the embodiments can be combined with each other. Based on the embodiments in the present disclosure, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present disclosure.

[0052] It should be noted that various aspects of the embodiments within the scope of the appended claims are described below. It should be apparent that the aspects described herein can be embodied in a wide variety of forms, and any specific structure and / or function described herein is merely illustrative. Based on this disclosure, it should be understood by those skilled in the art that an aspect described herein can be implemented independently of any other aspect, and two or more of these aspects can be combined in various ways. For example, any number of aspects described herein can be used to implement an apparatus and / or practice a method. In addition, other structures and / or functionalities other than one or more of the aspects described herein can be used to implement this apparatus and / or practice this method.

[0053] It should also be noted that the illustrations provided in the following embodiments are only schematic illustrations of the basic concept of the present disclosure. The illustrations only show components related to the present disclosure and are not drawn according to the number, shape and size of components in actual implementation. In actual implementation, the type, quantity and proportion of each component can be changed at will, and the component layout type may also be more complicated.

[0054] Additionally, in the following description, specific details are provided to provide a thorough understanding of the examples. However, one skilled in the art will appreciate that the aspects described can be practiced without these specific details.

[0055] In the first aspect, Figure 1 A method 100 for inverting and reconstructing soil moisture content provided by an embodiment of the present invention is shown, comprising the following steps:

[0056] In step S101, a ground penetrating radar system is used to detect the structure of the reconstructed soil layer to obtain electromagnetic wave propagation data of each layer of the reconstructed soil layer. Specifically, Figure 3 This diagram shows the wave velocity calibration principle of a ground-penetrating radar (GPR) system. This system is a non-invasive method for reconstructing soil layer structure. The electromagnetic wave propagation data for each layer includes the propagation time and velocity of each layer. By processing the received reflected signals, GPR images are generated. These images intuitively demonstrate the layered structure of the reconstructed soil layer and the interfaces between the layers. By observing and analyzing GPR images, information such as soil layer thickness, stratigraphic sequence, and material composition can be obtained.

[0057] In some embodiments, the reconstructed soil layer is a three-layer sponge structure, specifically including a topsoil ecological layer, an aquifer, and an aquiclude, and is compacted. The ecological layer is made of loess, the aquifer is made of yellow sandstone, and the aquiclude is made of red clay. During the reconstruction process, it is necessary to set it up according to the actual conditions of the area to be tested. At the same time, it is necessary to pre-treat the reconstructed soil layer using an electronic scale, a drying box, etc. to approximate the actual soil layer parameters of the area to be tested.

[0058] The ground-penetrating radar system includes at least an antenna, a data acquisition module, and a signal processing unit. The antenna is used to transmit and receive electromagnetic waves. The transmitting antenna transmits high-frequency electromagnetic waves underground, while the receiving antenna captures reflected waves from the interfaces of different underground media. The data acquisition module is used to record and process the reflected signals captured by the receiving antenna. The signal processing unit is used to further process and analyze the digital signals output by the data acquisition module, including filtering, correction, and imaging. Specifically, the ground-penetrating radar system preprocesses the detection data obtained by detection, and the preprocessing includes:

[0059] Background denoising is performed on the detection data to obtain denoised detection information; the background denoising includes removing environmental noise, equipment noise, electromagnetic interference, etc., usually using a sliding window averaging method or an adaptive filtering method.

[0060] The de-noised detection information is subjected to one-dimensional FIR filtering to obtain frequency-domain filtered detection information. This one-dimensional FIR filtering, also known as finite impulse response filtering, removes unwanted frequency components. The parameters of the one-dimensional FIR filter, including the filter order and cutoff frequency, are designed based on the characteristics of the radar signal and filtering requirements. The radar signal is input into the designed one-dimensional FIR filter for filtering. The one-dimensional FIR filter performs a weighted summation on the input signal to obtain a filtered output signal, which is the frequency-domain filtered detection information.

[0061] Gain control is performed on the frequency-domain filtered detection information to obtain optimized detection data. The gain control is used to adjust the gain of the receiver to maintain the stability and reliability of the output signal. Common gain control methods include automatic gain control or sensitivity time gain control.

[0062] Next, go to step S102.

[0063] In step S102, the reconstructed soil layer provided with a reflector is detected by a ground penetrating radar system to obtain two-way travel time data, the relative dielectric constant of the dielectric layer is obtained based on the two-way travel time data and the buried depth of the reflector, and an empirical formula is obtained based on the relative dielectric constant.

[0064] In some embodiments, the reconstructed soil layer provided with reflectors includes: a plurality of reflectors of different sizes spaced sequentially along the direction of radar electromagnetic wave emission, to ensure spatial coverage and representativeness of the data, thereby enhancing the signal reflection effect of the ground penetrating radar system and facilitating subsequent wave velocity calibration. Specifically, the reflectors are stainless steel reflectors, which can be stainless steel pipes, by selecting stainless steel pipes of different diameters. Furthermore, in some embodiments, the stainless steel pipes are provided at four locations, with diameters of 2 cm, 4 cm, 6 cm, and 8 cm, respectively, and the burial depths can be 20 cm, 40 cm, 60 cm, and 80 cm. Figure 4 A schematic diagram of the measurement point layout and steel pipe burial according to an embodiment of the present disclosure is shown. The steel pipes in the figure are respectively arranged at 25 cm and 45 cm of the reconstructed soil layer.

[0065] For example, two-way travel time data of each measuring point is obtained by ground penetrating radar, and the electromagnetic wave propagation speed in different dielectric layers is calculated in combination with the depth information of the buried reflector. Based on the electromagnetic wave propagation speed, the relative dielectric constant can be calculated. The relative dielectric constant is:

[0066]

[0067] Wherein, c is the speed of light; v is the propagation speed of radar electromagnetic waves, and the propagation speed of radar electromagnetic waves is calculated based on the two-way travel time data and the buried depth of the reflector.

[0068] Exemplarily, the empirical formula includes a Topp empirical formula, which includes:

[0069] θ v =-5.3×10 -2 +2.292×10 -2 ε r -5.5×10 -4 ε r 2+4.3×10 -6 ε r 3 ;

[0070] In the formula, the θ v is the moisture content; ε r is the relative dielectric constant.

[0071] Next, go to step S103.

[0072] In step S103, based on the actual moisture content of the reconstructed soil layer, the empirical formula is modified to obtain a moisture content model, and the moisture content of the reconstructed soil layer is obtained by inversion based on the moisture content model. In some embodiments, the empirical formula is modified using a polynomial regression analysis correction method, including:

[0073] Calculate the preliminary moisture content based on the empirical formula;

[0074] Based on the preliminary moisture content and the measured moisture content, a polynomial regression analysis correction method is used to obtain a moisture content model. Generally speaking, there is a certain difference between the preliminary moisture content and the measured moisture content. In this case, the empirical formula is corrected based on the actual measurement data through polynomial regression analysis.

[0075] Furthermore, the moisture content model includes:

[0076] θ1=-24.80918+11.67223ε r -1.19808ε r 2 +0.04163ε r 3 ;

[0077] Where θ1 is the moisture content; ε r is the relative dielectric constant, and the goodness of fit R 2 =0.93618, specifically, Figure 10 A fitting curve diagram of relative dielectric constant and volume water content of an embodiment of the present disclosure is shown.

[0078] In some embodiments, data processing software is also included, and the data processing software includes GR software and SPSS software. The data processing software is used to process the electromagnetic wave propagation data, data analysis, wave velocity calculation, empirical formula correction and water content inversion.

[0079] For example, Figure 2 As shown, the present application discloses a process 200 for inverting and reconstructing soil moisture content, comprising the following steps:

[0080] Step S201: Reconstructing the soil layer:

[0081] At the 100-acre (approximately 166 acres) demonstration site of the Heidaigou open-pit mine dump in Inner Mongolia, a three-layer sponge structure was constructed in the 3-meter and 1.5-meter zones, respectively. The 3-meter zone employed a 0.9-meter topsoil ecological layer, a 1.0-meter aquifer, and a 1.1-meter aquiclude. The 1.5-meter zone employed a 0.4-meter topsoil ecological layer, a 0.5-meter aquifer, and a 0.6-meter aquiclude. During construction, rollers were used to compact the layers layer by layer, ensuring clear stratification and distinct positions.

[0082] Step S202: Ground penetrating radar detection and wave velocity calibration:

[0083] Figure 5 The on-site detection diagram of the three-layer sponge structure area of ​​the embodiment of the present disclosure is shown as follows: Figure 5 As shown, a 400MHz GPR system was used to detect along the survey lines B1-B2 and C1-C2. Stainless steel pipes with diameters of 2cm, 4cm, 6cm, and 8cm were buried at different depths below points B1 and B2, respectively, at depths of 20cm, 40cm, 60cm, and 80cm, as reflectors. The reconstructed soil layer was detected using a ground penetrating radar system. Figure 6 FIG4 shows a radar detection diagram at a depth of 3 m according to an embodiment of the present disclosure. Figure 7 FIG4 shows a radar detection diagram at a depth of 1.5 m according to an embodiment of the present disclosure. Figure 8 A radar detection image of a 3m mixed zone in an embodiment of the present disclosure is shown. Using a ground-penetrating radar (GPR) system to measure two-way travel time and combine it with the burial depth, the electromagnetic wave propagation velocity and relative dielectric constant of each layer were calculated, resulting in 22 sets of dielectric constant data, as shown in Table 1.

[0084] Table 1. Error between the moisture content calculated using the Topp formula and the measured value in area B

[0085]

[0086]

[0087] Step S203: Water content inversion and model correction:

[0088] The moisture content of the soil layer in area B was preliminarily calculated using the Topp formula and compared with the measured moisture content of the soil layer. Figure 9A comparison curve of the TOPP formula calculated values ​​and the measured values ​​of the embodiment of the present disclosure is shown. Table 1 shows the data of the TOPP formula calculated values ​​and the measured values. It can be seen from the figure and the table that the difference between the TOPP formula calculated values ​​and the measured values ​​is large, that is, the error obtained by the TOPP formula is large. Therefore, polynomial regression analysis is used to correct the Topp formula and establish a model for the relationship between volumetric water content and dielectric constant applicable to this area: θ1 = -24.80918 + 11.67223εr - 1.19808εr 2 +0.04163εr 3 .

[0089] Step S204: Soil layer moisture content inversion output:

[0090] The revised formula, that is, the moisture content model, was applied to the data of area C to calculate the soil moisture content. The results were compared with the measured values. As shown in Table 2, the error range was significantly reduced, with an average error of only 0.523%, verifying the high accuracy and reliability of the inversion method.

[0091] Table 2. Error between the moisture content measured by the moisture content model and the measured value

[0092]

[0093] The second embodiment of the present invention further provides a device for inverting and reconstructing soil moisture content, the device comprising:

[0094] A collection unit is configured to collect and organize cathodoluminescence (CL) images of zircon, wherein characteristic information of the zircon in the CL images includes at least growth zoning of the zircon;

[0095] The zircon age prediction unit is configured to pre-process the zircon CL image using a preset target zircon age prediction model and obtain the zircon CL image age; wherein the preset target zircon age prediction model includes a ring band extraction module, a multi-scale feature fusion module and an age prediction module, the ring band extraction module is used to extract the ring band features in the CL image information, the multi-scale feature fusion module is used to perform multi-scale feature extraction and fusion on the ring band features to obtain multi-scale spatial information, and the age prediction module is used to obtain the final zircon age prediction result based on the multi-scale spatial information.

[0096] A third embodiment of the present invention further provides an electronic device, comprising:

[0097] at least one processor; and,

[0098] a memory communicatively connected to the at least one processor; wherein,

[0099] The memory stores instructions that can be executed by the at least one processor. The instructions are executed by the at least one processor to enable the at least one processor to execute the method for inverting and reconstructing soil layer moisture content in any of the aforementioned embodiments.

[0100] The fourth embodiment of the present invention also provides a non-transitory computer-readable storage medium, which stores computer instructions, and the computer instructions are used to enable the computer to execute the method for inverting and reconstructing the moisture content of the soil layer described in any of the aforementioned embodiments.

[0101] The fifth embodiment of the present invention also provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium. The computer program includes program instructions, which, when executed by a computer, enable the computer to execute the method for inverting and reconstructing the moisture content of the soil layer of any of the aforementioned embodiments.

[0102] Figure 11 A schematic diagram of a method or device 1000 that can implement an embodiment of the present invention is shown. In some embodiments, the method or device 1000 may include more or fewer devices than shown. In some embodiments, the method or device 1000 may be implemented using a single device or multiple devices. In some embodiments, the method or device 1000 may be implemented using cloud-based or distributed devices.

[0103] like Figure 11 As shown, device 1000 includes a processor 1001, which can perform various appropriate operations and processes according to the programs and / or data stored in a read-only memory (ROM) 1002 or the programs and / or data loaded from a storage portion 1008 into a random access memory (RAM) 1003. Processor 1001 can be a multi-core processor or can include multiple processors. In some embodiments, processor 1001 can include a general-purpose main processor and one or more special coprocessors, such as a central processing unit (CPU), a graphics processing unit (GPU), a neural network processor (NPU), a digital signal processor (DSP), etc. Various programs and data required for the operation of device 1000 are also stored in RAM 1003. Processor 1001, ROM 1002, and RAM 1003 are connected to each other via a bus 1004. An input / output (I / O) interface 1005 is also connected to bus 1004.

[0104] The processor and memory are used together to execute the program stored in the memory. When the program is executed by the computer, the methods, steps or functions described in the above embodiments can be implemented.

[0105] The following components are connected to the I / O interface 1005: an input section 1006 including a keyboard, a mouse, a touch screen, and the like; an output section 1007 including devices such as a cathode ray tube (CRT), a liquid crystal display (LCD), and a speaker; a storage section 1008 including a hard disk and the like; and a communication section 1009 including a network interface card such as a LAN card or a modem. The communication section 1009 performs communication processing via a network such as the Internet. A drive 1010 is also connected to the I / O interface 1005 as needed. A removable medium 1011, such as a magnetic disk, an optical disk, a magneto-optical disk, or a semiconductor memory, is installed in the drive 1010 as needed, so that computer programs read therefrom can be installed into the storage section 1008 as needed. Figure 11 Only some components are shown schematically, which does not mean that the device 1000 only includes Figure 11 Components shown.

[0106] The systems, devices, modules, or units described in the above embodiments may be implemented by a computer or its associated components. The computer may be, for example, a mobile terminal, a smartphone, a personal computer, a laptop computer, an in-vehicle human-computer interaction device, a personal digital assistant, a media player, a navigation device, a game console, a tablet computer, a wearable device, a smart TV, an Internet of Things system, a smart home, an industrial computer, a server, or a combination thereof.

[0107] Although not shown, in an embodiment of the present invention, a computer-readable storage medium is provided, on which a computer program / instruction is stored. When the computer program / instruction is executed by a processor, the method for inverting and reconstructing the moisture content of a soil layer described in the embodiment is implemented.

[0108] Storage media in embodiments of the present invention include permanent and non-permanent, removable and non-removable items that can be used to store information using any method or technology. Examples of storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device.

[0109] Although not shown, an embodiment of the present invention further provides a computer program product, including: a computer program / instruction, which, when executed by a processor, implements the method for inverting and reconstructing soil layer moisture content described in the embodiment.

[0110] The methods, programs, systems, and apparatuses of the embodiments of the present invention may be executed or implemented in a single or multiple networked computers, or may be practiced in a distributed computing environment. In the embodiments of this specification, in these distributed computing environments, tasks may be performed by remote processing devices connected via a communication network.

[0111] Those skilled in the art will appreciate that the embodiments of this specification may be provided as methods, systems, or computer program products. Therefore, those skilled in the art will appreciate that the functional modules / units or controllers and related method steps described in the above embodiments may be implemented using software, hardware, or a combination of software / hardware.

[0112] Unless explicitly stated, the actions or steps of the methods, procedures, and methods described in accordance with the embodiments of the present invention do not have to be performed in a specific order and can still achieve the desired results. In some implementations, multitasking and parallel processing are also possible or may be advantageous.

[0113] In this document, multiple embodiments of the present invention are described, but for the sake of brevity, the description of each embodiment is not exhaustive, and the same or similar features or parts between the embodiments may be omitted. In this document, "one embodiment", "some embodiments", "example", "specific example", or "some examples" are intended to apply to at least one embodiment or example according to the present invention, but not all embodiments. The above terms do not necessarily mean to refer to the same embodiment or example. Those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, unless they are mutually contradictory.

[0114] While the exemplary systems and methods of the present invention have been specifically shown and described with reference to the foregoing embodiments, these are merely examples of the best modes for implementing the present systems and methods. Those skilled in the art will appreciate that various changes may be made to the embodiments of the systems and methods described herein when implementing the present systems and / or methods without departing from the spirit and scope of the present invention as defined in the appended claims.

[0115] In addition, the method, system, device, and medium for inverting and reconstructing soil moisture content according to the present invention can also be implemented in the following ways:

[0116] (1) A method for inverting and reconstructing soil moisture content, characterized by comprising:

[0117] Use a ground penetrating radar system to detect the structure of the reconstructed soil layer and obtain electromagnetic wave propagation data of each layer of the reconstructed soil layer;

[0118] Detecting the reconstructed soil layer provided with a reflector using a ground penetrating radar system to obtain two-way travel time data, obtaining a relative dielectric constant of the dielectric layer based on the two-way travel time data and the buried depth of the reflector, and obtaining an empirical formula based on the relative dielectric constant;

[0119] Based on the actual moisture content of the reconstructed soil layer, the empirical formula is modified to obtain a moisture content model, and the moisture content of the reconstructed soil layer is obtained by inversion based on the moisture content model.

[0120] (2) The method according to (1) is characterized in that the reconstructed soil layer provided with reflectors includes: a plurality of reflectors of different sizes arranged in sequence and spaced apart along the direction of emission of the radar electromagnetic wave.

[0121] (3) The method according to (1), characterized in that the ground penetrating radar system preprocesses the detection data obtained by detection, and the preprocessing includes:

[0122] performing background denoising on the detection data to obtain denoised detection information;

[0123] Performing one-dimensional FIR filtering on the denoised detection information to obtain frequency-domain filtered detection information;

[0124] Gain control is performed on the frequency domain filtered detection information to obtain optimized detection data.

[0125] (4) The method according to (1), characterized in that the empirical formula includes a Topp empirical formula, and the Topp empirical formula includes:

[0126] θ v =-5.3×10 -2 +2.292×10 -2 ε r -5.5×10 -4 ε r 2 +4.3×10 -6 ε r 3 ;

[0127] In the formula, the θ v is the moisture content; ε r is the relative dielectric constant.

[0128] (5) The method according to (1), characterized in that the relative dielectric constant is:

[0129]

[0130] Wherein, c is the speed of light; v is the propagation speed of radar electromagnetic waves, and the propagation speed of radar electromagnetic waves is calculated based on the two-way travel time data and the buried depth of the reflector.

[0131] (6) The method according to (1) is characterized in that the correction of the empirical formula adopts a polynomial regression analysis correction method, including:

[0132] Calculate the preliminary moisture content based on the empirical formula;

[0133] Based on the preliminary moisture content and the measured moisture content, the moisture content model was obtained by using the polynomial regression analysis correction method.

[0134] (7) The method according to (1), characterized in that the moisture content model includes:

[0135] θ1=-24.80918+11.67223ε r -1.19808ε r 2 +0.04163ε r 3 ;

[0136] Where θ1 is the moisture content; ε r is the relative dielectric constant.

[0137] (8) A device for inverting and reconstructing soil moisture content, characterized in that the device comprises:

[0138] The primary acquisition unit is configured to detect the structure of the reconstructed soil layer using a ground penetrating radar system to obtain electromagnetic wave propagation data of each layer of the reconstructed soil layer;

[0139] a secondary acquisition unit configured to detect the reconstructed soil layer provided with a reflector using a ground penetrating radar system to obtain two-way travel time data, obtain a relative dielectric constant of the dielectric layer based on the two-way travel time data and the burial depth of the reflector, and obtain an empirical formula based on the relative dielectric constant;

[0140] The output unit is configured to correct the empirical formula based on the actual moisture content of the reconstructed soil layer to obtain a moisture content model, and invert the moisture content of the reconstructed soil layer based on the moisture content model.

[0141] (9) The device according to (8) is characterized in that the primary acquisition unit is further configured so that the reconstructed soil layer provided with reflectors includes: a plurality of reflectors of different sizes arranged in sequence and spaced apart along the direction of emission of the radar electromagnetic wave.

[0142] (10) The device according to (8), characterized in that the primary acquisition unit is further configured to pre-process the detection data obtained by the ground penetrating radar system, and the pre-processing includes:

[0143] performing background denoising on the detection data to obtain denoised detection information;

[0144] Performing one-dimensional FIR filtering on the denoised detection information to obtain frequency-domain filtered detection information;

[0145] Gain control is performed on the frequency domain filtered detection information to obtain optimized detection data.

[0146] (11) The device according to (8), characterized in that the secondary acquisition unit is further configured such that the empirical formula includes a Topp empirical formula, and the Topp empirical formula includes:

[0147] θ v =-5.3×10 -2 +2.292×10 -2 ε r -5.5×10 -4 ε r 2 +4.3×10 -6 ε r 3 ;

[0148] In the formula, the θ v is the moisture content; ε r is the relative dielectric constant.

[0149] (12) The device according to (8), characterized in that the secondary acquisition unit is further configured such that the relative dielectric constant is:

[0150]

[0151] Wherein, c is the speed of light; v is the propagation speed of radar electromagnetic waves, and the propagation speed of radar electromagnetic waves is calculated based on the two-way travel time data and the buried depth of the reflector.

[0152] (13) The device according to (8), characterized in that the output unit is further configured to correct the empirical formula using a polynomial regression analysis correction method, including:

[0153] Calculate the preliminary moisture content based on the empirical formula;

[0154] Based on the preliminary moisture content and the measured moisture content, the moisture content model was obtained by using the polynomial regression analysis correction method.

[0155] (14) The device according to (8), characterized in that the output unit is further configured such that the moisture content model includes:

[0156] θ1=-24.80918+11.67223ε r -1.19808ε r 2 +0.04163ε r 3 ;

[0157] Where θ1 is the moisture content; ε r is the relative dielectric constant.

[0158] (15) An electronic device, characterized in that the electronic device comprises:

[0159] at least one processor; and,

[0160] a memory communicatively connected to the at least one processor; wherein,

[0161] The memory stores instructions that can be executed by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute the method for inverting and reconstructing the moisture content of the soil layer as described in any one of (1) to (14).

[0162] (16) A non-transitory computer-readable storage medium, characterized in that the non-transitory computer-readable storage medium stores computer instructions, which are used to enable the computer to execute the method for inverting and reconstructing soil layer moisture content described in any one of (1) to (14).

[0163] (17) A computer program product comprising a computer program stored on a non-transitory computer-readable storage medium, the computer program comprising program instructions which, when executed by a computer, cause the computer to execute the method for inverting and reconstructing soil layer moisture content as described in any one of (1) to (14).

Claims

1. A method for inverting and reconstructing soil moisture content, characterized in that: include: Use a ground penetrating radar system to detect the structure of the reconstructed soil layer and obtain electromagnetic wave propagation data of each layer of the reconstructed soil layer; Detecting the reconstructed soil layer provided with a reflector using a ground penetrating radar system to obtain two-way travel time data, obtaining a relative dielectric constant of the dielectric layer based on the two-way travel time data and the buried depth of the reflector, and obtaining an empirical formula based on the relative dielectric constant; Based on the actual moisture content of the reconstructed soil layer, the empirical formula is modified to obtain a moisture content model, and the moisture content of the reconstructed soil layer is obtained by inversion based on the moisture content model.

2. The method according to claim 1, characterized in that The reconstructed soil layer provided with reflectors includes: a plurality of reflectors of different sizes sequentially spaced apart along the direction of emission of radar electromagnetic waves.

3. The method according to claim 1, characterized in that The ground penetrating radar system pre-processes the detection data obtained by detection, and the pre-processing includes: performing background denoising on the detection data to obtain denoised detection information; Performing one-dimensional FIR filtering on the denoised detection information to obtain frequency-domain filtered detection information; Gain control is performed on the frequency domain filtered detection information to obtain optimized detection data.

4. The method according to claim 1, wherein The empirical formula includes the Topp empirical formula, which includes: i v =-5.3×10 -2 +2.292×10 -2 e r -5.5×10 -4 e r 2 +4.3×10 -6 e r 3 ; In the formula, the θ v is the moisture content; ε r is the relative dielectric constant.

5. The method according to claim 1, wherein The relative dielectric constant is: Wherein, c is the speed of light; v is the propagation speed of radar electromagnetic waves, and the propagation speed of radar electromagnetic waves is calculated based on the two-way travel time data and the buried depth of the reflector.

6. The method according to claim 1, wherein The empirical formula is corrected using a polynomial regression analysis correction method, including: Calculate the preliminary moisture content based on the empirical formula; Based on the preliminary moisture content and the measured moisture content, the moisture content model was obtained by using the polynomial regression analysis correction method.

7. The method according to claim 1, characterized in that The moisture content model includes: θ1=-24.80918+11.67223ε r -1.19808e r 2 +0.04163e r 3 ; Where θ1 is the moisture content; ε r is the relative dielectric constant.

8. A device for inverting and reconstructing soil moisture content, characterized in that: Based on the method for inverting and reconstructing soil moisture content according to any one of claims 1 to 7, the device comprises: The primary acquisition unit is configured to detect the structure of the reconstructed soil layer using a ground penetrating radar system to obtain electromagnetic wave propagation data of each layer of the reconstructed soil layer; a secondary acquisition unit configured to detect the reconstructed soil layer provided with a reflector using a ground penetrating radar system to obtain two-way travel time data, obtain a relative dielectric constant of the dielectric layer based on the two-way travel time data and the burial depth of the reflector, and obtain an empirical formula based on the relative dielectric constant; The output unit is configured to correct the empirical formula based on the actual moisture content of the reconstructed soil layer to obtain a moisture content model, and invert the moisture content of the reconstructed soil layer based on the moisture content model.

9. An electronic device, characterized in that: The electronic device includes: at least one processor; and, a memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute the method for inverting and reconstructing soil layer moisture content as described in any one of claims 1 to 7.

10. A non-transitory computer-readable storage medium, characterized in that The non-transitory computer-readable storage medium stores computer instructions, which are used to enable the computer to execute the method for inverting and reconstructing soil layer moisture content according to any one of claims 1 to 7.

Images