CN111239726A - Method, device and readable storage medium for identifying underground pollutant area - Google Patents
Method, device and readable storage medium for identifying underground pollutant area Download PDFInfo
- Publication number
- CN111239726A CN111239726A CN202010120610.0A CN202010120610A CN111239726A CN 111239726 A CN111239726 A CN 111239726A CN 202010120610 A CN202010120610 A CN 202010120610A CN 111239726 A CN111239726 A CN 111239726A
- Authority
- CN
- China
- Prior art keywords
- dielectric constant
- identifying
- stratum
- abnormal
- zone
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 239000003344 environmental pollutant Substances 0.000 title claims abstract description 19
- 231100000719 pollutant Toxicity 0.000 title claims abstract description 19
- 230000002159 abnormal effect Effects 0.000 claims abstract description 20
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 15
- 238000001514 detection method Methods 0.000 claims abstract description 14
- 230000000149 penetrating effect Effects 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 7
- 230000006835 compression Effects 0.000 claims abstract description 5
- 238000007906 compression Methods 0.000 claims abstract description 5
- 230000010354 integration Effects 0.000 claims abstract description 4
- 238000011109 contamination Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 5
- 239000000356 contaminant Substances 0.000 claims 1
- 238000011084 recovery Methods 0.000 claims 1
- 241000282414 Homo sapiens Species 0.000 description 5
- 239000008346 aqueous phase Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003895 groundwater pollution Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000009439 industrial construction Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
- G01N27/06—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a liquid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/22—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance
- G01N27/221—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance by investigating the dielectric properties
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electromagnetism (AREA)
- Electrochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method, equipment and readable storage medium for identifying underground pollutant areas, which uses a ground penetrating radar to detect underground preset areas; restoring the amplitude of the section; performing compression wavelet processing on the detection signal through pulse deconvolution; obtaining the reflection coefficient of the stratum by a sparse deconvolution method based on an L1 norm; obtaining an estimated wave impedance value by performing bandwidth-limited integration on the reflection coefficient; calculating a global impedance value through global constraints; calculating the dielectric constant and the water content to generate a three-dimensional slice image, and observing the abnormal dielectric constant on the depth slice through time slicing; the abnormal region band is divided based on the abnormal dielectric constant. The method and the device realize finding out the position of the pollutant without damaging the stratum, avoid damaging the stratum structure and achieve the aim of treating the problem that the pollutant pollutes the environment.
Description
Technical Field
The invention relates to the field of radar data processing, in particular to a method, equipment and a readable storage medium for identifying an underground pollutant area.
Background
Underground water is a precious resource on which human beings live, and the guarantee of water safety is a prerequisite for human health, but due to the developed industrial construction of human beings, the underground water pollution prevention and control becomes a problem facing human beings at present while bringing convenience and progress to human beings and causing serious damage to the ecological environment. The underground water detection work is not suitable for adopting a well-digging sampling detection method, the cost is high, and the stratum structure is damaged.
Currently, Ground Penetrating Radar (GPR) is a mature geophysical exploration tool, and can be used for monitoring the activity of underground water in real time and detecting Light non-aqueous phase pollutants (LNAPL). People put higher requirements on the detection capability of the ground penetrating radar, hope to obtain physical parameters in hydrogeology and make quantitative explanation and analysis on underground water and pollution conditions.
The physical parameters of the underground medium have great influence on the detection of GPR, electromagnetic waves propagate in different underground media, the propagation speed and the amplitude of the electromagnetic waves are changed, and the physical parameters influencing the propagation of the electromagnetic waves are the dielectric constant and the conductivity of the medium, which respectively influence the speed and the amplitude of the waves. The conductivity range of the light non-aqueous phase fluid is generally 0-0.02S/m, the relative dielectric constant range is 1.8-3, and the light non-aqueous phase fluid belongs to a low-loss medium and is suitable for high-frequency electromagnetic wave detection. GPR as a new nondestructive detection technology has the advantages of high efficiency, high signal-to-noise ratio, high resolution and the like, so that GPR is very suitable for detecting and monitoring LNAPLs pollution.
How to realize effective measurement of groundwater pollution area based on GPR is a technical problem to be solved at present.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention provides a method for identifying an underground pollutant area, which is characterized by comprising the following steps:
s1: detecting an underground preset area by using a ground penetrating radar;
s2: restoring the amplitude of the section;
s3: performing compression wavelet processing on the detection signal through pulse deconvolution;
s4: obtaining the reflection coefficient of the stratum by a sparse deconvolution method based on an L1 norm;
s5: obtaining an estimated wave impedance value by performing bandwidth-limited integration on the reflection coefficient;
s6: calculating a global impedance value through global constraints;
s7: calculating the dielectric constant and the water content to generate a three-dimensional slice image, and observing the abnormal dielectric constant on the depth slice through time slicing; the abnormal region band is divided based on the abnormal dielectric constant.
The invention also provides a device for realizing the method for identifying the underground pollutant area, which comprises the following steps:
a memory for storing a computer program and a method of identifying a zone of subsurface contamination;
a processor for executing the computer program and the method for identifying a subterranean zone of contamination to carry out the steps of the method for identifying a subterranean zone of contamination.
The present invention also provides a readable storage medium having a method of identifying a zone of subsurface contamination, the readable storage medium having stored thereon a computer program for execution by a processor to perform the steps of the method of identifying a zone of subsurface contamination.
According to the technical scheme, the invention has the following advantages:
the method comprises the steps of detecting an underground preset area by using a ground penetrating radar; restoring the amplitude of the section; performing compression wavelet processing on the detection signal through pulse deconvolution; calculating the dielectric constant and the water content in the corresponding process to generate a three-dimensional slice image, and observing the abnormal dielectric constant on the depth slice through time slicing; the abnormal region band is divided based on the abnormal dielectric constant. The method and the device realize finding out the position of the pollutant without damaging the stratum, avoid damaging the stratum structure and achieve the aim of treating the problem that the pollutant pollutes the environment. The cost output is reduced to the maximum extent, the detection time is shortened, the detection range is expanded, and all-weather real-time monitoring can be carried out to meet the detection requirement of the environmental pollution problem caused by urban development.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings used in the description will be briefly introduced, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained based on these drawings without creative efforts.
FIG. 1 is a flow chart of a method of identifying a zone of subsurface contamination;
FIG. 2 is a flow chart of an embodiment of a method of identifying a subsurface pollutant zone.
Detailed Description
Those of ordinary skill in the art will appreciate that the elements and algorithm steps of the examples described in connection with the embodiments disclosed herein may be embodied in electronic hardware, computer software, or combinations of both, and that the components and steps of the examples have been described in a functional general in the foregoing description for the purpose of illustrating clearly the interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The block diagrams shown in the figures are functional entities only and do not necessarily correspond to physically separate entities. I.e. these functional entities may be implemented in the form of software, or in one or more hardware modules or integrated circuits, or in different networks and/or processor means and/or microcontroller means.
In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may also be an electric, mechanical or other form of connection.
The present invention provides a method for identifying a zone of subsurface contamination, as shown in fig. 1 and 2, comprising:
s1: detecting an underground preset area by using a ground penetrating radar;
in the first step, ground penetrating radar section data are prepared, and basic data preprocessing operations including exponential gain, prediction deconvolution, band-pass filtering and the like are performed.
The acquired data is processed, here data editing of the acquired GPR profiles. And (4) removing bad tracks, unifying track spacing and track number of each section, and unifying time sampling points and sampling intervals. And removing interference caused by the instrument and removing direct current drift. Renaming of files, etc. If the sections are many, the data can be conveniently arranged and edited for subsequent batch processing.
S2: restoring the amplitude of the section;
the attenuation of GPR electromagnetic waves in the stratum is fast, if the signal amplitude is not recovered, deep stratum reflection information cannot be obtained, and a plurality of gain methods such as linear gain, manual gain and the like exist.
S3: performing compression wavelet processing on the detection signal through pulse deconvolution;
in a received radar signal including multiples formed by multiple reflections between layers, occurring at fixed time intervals after the primary wave is formed, the effects of ringing caused by the multiple reflections can be eliminated by predictive deconvolution. In order to obtain the desired sharp pulse as the received signal, the wavelet is compressed by pulse deconvolution, thereby improving the longitudinal resolution of the signal.
S4: obtaining the reflection coefficient of the stratum by a sparse deconvolution method based on an L1 norm;
the GPR signal can be regarded as a convolution of a reflection coefficient and a radar wavelet, and the reflection coefficient of the stratum is obtained through a sparse deconvolution method based on an L1 norm.
S5: obtaining an estimated wave impedance value by performing bandwidth-limited integration on the reflection coefficient;
s6: calculating a global impedance value through global constraints;
the invention calculates a global impedance value based on global constraints. The wave impedance value inverted from the amplitude causes an error in the impedance amplitude range due to some uncertainty, and is generally combined with environmental hydrological data, such as well drilling data. Through well head data, the depth wave impedance value of a point can be obtained, namely the real impedance value, effective constraint can be carried out through the most calculated wave impedance of the well data under the environment that the stratum is layered and continuous, and the more the number of the wells, the better the constraint effect.
S7: calculating the dielectric constant and the water content to generate a three-dimensional slice image, and observing the abnormal dielectric constant on the depth slice through time slicing; the abnormal region band is divided based on the abnormal dielectric constant.
The method calculates the stratum dielectric constant and further calculates the stratum water content by mixing the hydrological medium model with the obtained wave impedance. And finally, generating a three-dimensional slice image, and observing the place with abnormal dielectric constant on the depth slice through the time slice so as to divide the distribution range of the place. The slice images were drawn using the slice function of MATLAB.
Specifically, a plurality of pieces of profile data are integrated into three-dimensional data. The time sampling points and the channel number of each section are required to be set to be uniform values so as to facilitate the three-dimensional mapping display.
And performing wave impedance inversion on the three-dimensional radar data. Firstly, the reflection coefficient of data is solved, the invention applies an L1 norm method, and experiments prove that the reflection coefficient is closer to a true value than the reflection coefficient solved by the damped least square method, and the error is smaller. And further solving the formation wave impedance value through the reflection coefficient, and adding a background value to make constraint. And finally, information such as the formation dielectric constant, the water content and the like is inverted through wave impedance.
And as the depth of the stratum increases, the dielectric constant value is gradually increased under the influence of water content, and when the dielectric constant value of the stratum in a certain area is found to be greatly different from the surrounding medium, the dielectric constant value is obviously smaller. The dielectric constant of the light non-aqueous phase fluid is between 1.8 and 3, but the light non-aqueous phase fluid is mixed with a stratum medium and is influenced by the water content, the mixed dielectric constant can be larger than 3, the specific situation refers to actual stratum background data, such as well drilling data, and an abnormal zone can be marked out by combining the analysis of the stratum background data.
The method is visually displayed in a three-dimensional mapping mode.
The slice image can be drawn by using the slice function of MATLAB in a specific mapping mode, and the slice image of the three-dimensional data at any point of any coordinate can be randomly checked by the slice function, so that an image with an obvious imaging effect can be conveniently screened out.
For example, calling the function: slice3D is plot as plot of plot3Dslice (data _3D, a, B, C, 1). Wherein data _3D is processed three-dimensional data, and A, B and C are coordinate vectors respectively corresponding to X, Y and Z coordinates. a, B and C are corresponding coordinate points, each point forms a slice image, and a, B and C can be a plurality of points as long as a belongs to A, B belongs to B and C belongs to C.
Based on the method, the invention also provides equipment for realizing the method for identifying the underground pollutant area, which comprises the following steps: a memory for storing a computer program and a method of identifying a zone of subsurface contamination; a processor for executing the computer program and the method for identifying a subterranean zone of contamination to carry out the steps of the method for identifying a subterranean zone of contamination.
The invention also provides a readable storage medium with a method for identifying a subsurface pollutant zone based on the above method, the readable storage medium having stored thereon a computer program, the computer program being executed by a processor to perform the steps of the method for identifying a subsurface pollutant zone.
The apparatus implementing the method for identifying a zone of subsurface contamination is the exemplary elements and algorithm steps described in connection with the embodiments disclosed herein, which can be implemented in electronic hardware, computer software, or combinations of both, the components and steps of the various examples having been described above generally in terms of their functionality for clarity of illustration of interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
Through the above description of the embodiments, those skilled in the art will readily understand that the apparatus for implementing the method for identifying a subsurface pollutant region described herein can be implemented by software, and can also be implemented by software in combination with necessary hardware. Therefore, the technical solution according to the disclosed embodiment of the apparatus for implementing the method for identifying an underground pollutant area can be embodied in the form of a software product, which can be stored in a non-volatile storage medium (which can be a CD-ROM, a usb disk, a mobile hard disk, etc.) or on a network, and includes several instructions to make a computing device (which can be a personal computer, a server, a mobile terminal, or a network device, etc.) execute the indexing method according to the disclosed embodiment.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (7)
1. A method of identifying a subterranean contaminant zone, comprising:
s1: detecting an underground preset area by using a ground penetrating radar;
s2: restoring the amplitude of the section;
s3: performing compression wavelet processing on the detection signal through pulse deconvolution;
s4: obtaining the reflection coefficient of the stratum by a sparse deconvolution method based on an L1 norm;
s5: obtaining an estimated wave impedance value by performing bandwidth-limited integration on the reflection coefficient;
s6: calculating a global impedance value through global constraints;
s7: calculating the dielectric constant and the water content to generate a three-dimensional slice image, and observing the abnormal dielectric constant on the depth slice through time slicing; the abnormal region band is divided based on the abnormal dielectric constant.
2. The method of claim 1,
step S2 further includes:
the exponential gain adopted for the amplitude recovery of the section is adopted, and the obtained stratum attenuation factor is used as a gain coefficient, so that the amplitude value of the original signal can be better restored.
3. The method of claim 1,
step S7 further includes:
the slice images were drawn using the slice function of MATLAB.
4. The method of claim 1,
step S7 further includes: and setting the time sampling point of each section and the number of the sections to be uniform values, and configuring the time sampling points and the number of the sections into a three-dimensional map for displaying.
5. The method of claim 1,
step S7 further includes:
the method comprises the steps that the dielectric constant value of a stratum is gradually increased along with the increase of the depth under the influence of water content, when the dielectric constant value of the stratum in a certain area is compared with the dielectric constant of the medium around the certain area and exceeds a threshold value, or when the dielectric constant value of the stratum in the certain area is larger than a preset threshold value, the area is divided into places with abnormal dielectric constants, and abnormal area bands are divided based on the abnormal dielectric constants.
6. An apparatus for implementing a method for identifying a zone of subsurface contamination, comprising:
a memory for storing a computer program and a method of identifying a zone of subsurface contamination;
a processor for executing the computer program and the method for identifying a subsurface pollutant zone to realize the steps of the method for identifying a subsurface pollutant zone according to any one of claims 1 to 5.
7. A readable storage medium having a method of identifying a subterranean zone of contamination, the readable storage medium having stored thereon a computer program for execution by a processor to perform the steps of the method of identifying a subterranean zone of contamination of any of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010120610.0A CN111239726A (en) | 2020-02-26 | 2020-02-26 | Method, device and readable storage medium for identifying underground pollutant area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010120610.0A CN111239726A (en) | 2020-02-26 | 2020-02-26 | Method, device and readable storage medium for identifying underground pollutant area |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111239726A true CN111239726A (en) | 2020-06-05 |
Family
ID=70875060
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010120610.0A Pending CN111239726A (en) | 2020-02-26 | 2020-02-26 | Method, device and readable storage medium for identifying underground pollutant area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111239726A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113092702A (en) * | 2021-04-06 | 2021-07-09 | 桂林理工大学 | Water body pollution detection method based on ground penetrating radar-seismic wave amplitude attribute |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107356985A (en) * | 2017-08-24 | 2017-11-17 | 浙江大学 | The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion |
| CN107544097A (en) * | 2017-06-27 | 2018-01-05 | 上海市环境科学研究院 | A kind of soil pollution based on geophysical probing technique is accurately positioned and accurate evaluation method |
| CN110618453A (en) * | 2019-08-07 | 2019-12-27 | 成都理工大学 | Wave impedance inversion method based on improved damping least square method |
-
2020
- 2020-02-26 CN CN202010120610.0A patent/CN111239726A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107544097A (en) * | 2017-06-27 | 2018-01-05 | 上海市环境科学研究院 | A kind of soil pollution based on geophysical probing technique is accurately positioned and accurate evaluation method |
| CN107356985A (en) * | 2017-08-24 | 2017-11-17 | 浙江大学 | The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion |
| CN110618453A (en) * | 2019-08-07 | 2019-12-27 | 成都理工大学 | Wave impedance inversion method based on improved damping least square method |
Non-Patent Citations (1)
| Title |
|---|
| 佘松盛: "探地雷达方法对近地表LNAPL污染物的探测研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113092702A (en) * | 2021-04-06 | 2021-07-09 | 桂林理工大学 | Water body pollution detection method based on ground penetrating radar-seismic wave amplitude attribute |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11409023B2 (en) | Methods to handle discontinuity in constructing design space using moving least squares | |
| Bergbauer et al. | Improving curvature analyses of deformed horizons using scale-dependent filtering techniques | |
| Lark et al. | Changes in variance and correlation of soil properties with scale and location: analysis using an adapted maximal overlap discrete wavelet transform | |
| Wang et al. | A wavelet-based baseline drift correction method for grounded electrical source airborne transient electromagnetic signals | |
| Bignardi et al. | OpenHVSR-Processing toolkit: Enhanced HVSR processing of distributed microtremor measurements and spatial variation of their informative content | |
| Foroozan et al. | Time-reversal ground-penetrating radar: Range estimation with Cramér–Rao lower bounds | |
| JP5719075B1 (en) | Cavity thickness exploration method | |
| Um et al. | Finite-element analysis of top-casing electric source method for imaging hydraulically active fracture zones | |
| Li et al. | Automatic fracture–vug identification and extraction from electric imaging logging data based on path morphology | |
| Lv et al. | Noise removal for semi-airborne data using wavelet threshold and singular value decomposition | |
| He et al. | GPR image denoising with NSST-UNET and an improved BM3D | |
| CN111239726A (en) | Method, device and readable storage medium for identifying underground pollutant area | |
| Giannakis | Realistic numerical modelling of ground penetrating radar for landmine detection | |
| Lei et al. | GPR detection localization of underground structures based on deep learning and reverse time migration | |
| Lu et al. | The tunnel seismic advance prediction method with wide illumination and a high signal‐to‐noise ratio | |
| Chen et al. | A data-driven multidimensional signal-noise decomposition approach for GPR data processing | |
| CN106610507A (en) | Method and device for removing ghost reflection interference | |
| Chen et al. | Extraction of gravity anomalies associated with gold mineralization: a comparison of singular value decomposition and bi-dimensional empirical mode decomposition | |
| Chen et al. | Ground penetrating radar signal processing improves mapping accuracy of underground voids and seawater table: an application in deteriorating coastal structure, Nanfangao Port, Taiwan | |
| Tang et al. | Geological disaster survey based on Curvelet transform with borehole Ground Penetrating Radar in Tonglushan old mine site | |
| Zhang et al. | A method for eliminating caprock thickness influence on anomaly intensities in geochemical surface survey for hydrocarbons | |
| Ferahtia et al. | A fuzzy logic-based filter for the removal of spike noise from 2D electrical resistivity data | |
| Reza et al. | Maximum entropy estimation of density function using order statistics | |
| CN113031072A (en) | Method, device and equipment for suppressing multiple between virtual homophase axial layers | |
| Sicking et al. | Pre-drill reservoir evaluation using passive seismic imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200605 |