CN104143115A - Technological method for achieving soil water content classified identification through geological radar technology - Google Patents
Technological method for achieving soil water content classified identification through geological radar technology Download PDFInfo
- Publication number
- CN104143115A CN104143115A CN201410286867.8A CN201410286867A CN104143115A CN 104143115 A CN104143115 A CN 104143115A CN 201410286867 A CN201410286867 A CN 201410286867A CN 104143115 A CN104143115 A CN 104143115A
- Authority
- CN
- China
- Prior art keywords
- data
- classification
- soil water
- identification
- geological radar
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 46
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005516 engineering process Methods 0.000 title claims abstract description 16
- 238000001228 spectrum Methods 0.000 claims abstract description 30
- 238000013528 artificial neural network Methods 0.000 claims abstract description 19
- 238000001914 filtration Methods 0.000 claims abstract description 13
- 230000009466 transformation Effects 0.000 claims abstract description 13
- 238000010606 normalization Methods 0.000 claims abstract description 9
- 239000013598 vector Substances 0.000 claims abstract description 9
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 238000012937 correction Methods 0.000 claims abstract description 6
- 238000000513 principal component analysis Methods 0.000 claims description 20
- 238000012549 training Methods 0.000 claims description 14
- 230000000704 physical effect Effects 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 4
- 238000010586 diagram Methods 0.000 claims description 3
- 238000000151 deposition Methods 0.000 claims description 2
- 230000001537 neural effect Effects 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 11
- 238000007781 pre-processing Methods 0.000 abstract 2
- 239000011159 matrix material Substances 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- FNMKZDDKPDBYJM-UHFFFAOYSA-N 3-(1,3-benzodioxol-5-yl)-7-(3-methylbut-2-enoxy)chromen-4-one Chemical compound C1=C2OCOC2=CC(C2=COC=3C(C2=O)=CC=C(C=3)OCC=C(C)C)=C1 FNMKZDDKPDBYJM-UHFFFAOYSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000005654 stationary process Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000007805 zymography Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a technological method for achieving soil water content classified identification through a geological radar technology. The technological method comprises the four parts of data preprocessing, soil water content parameter extraction, neural network classified identification and result display. Data preprocessing comprises the steps of null line correction, wavelet transformation and lowpass filtering. Soil water content parameter extraction comprises the steps that the power spectrum of reflection signals is solved with an autoregression moving average spectrum estimation method, data normalization is conducted, feature vectors and feature values of the power spectrum are extracted with a principle component analyzing method, and a sample feature vector database is constructed. Neural network classified identification comprises the steps that a neural network is trained through sample feature vectors, the trained mature network is used for conducting classified identification on data to be identified. Result display comprises the step that classified results are mapped and displayed. According to the technological method for achieving soil water content classified identification through the geological radar technology, automatic fast classified identification of the water content of soil is achieved, and a guiding function is performed on land detection and land rehabilitation.
Description
Technical field
The present invention relates to Georadar Data and process and interpretation technique, particularly a kind of geological radar technology is realized the technical method of soil water-containing Classification and Identification.
Background technology
Geological radar technology is a science of carrying out buried target detection by transmitting high-frequency impulse electromagnetic wave (frequency range is at 106 ~ 109Hz).The features such as that geological radar has is simple to operate, detection accuracy is high, not damaged, picking rate are fast, are current engineering detecting and the most active technical method of prospecting, and the application in Geotechnical Engineering is increasingly extensive.
Geological radar data is processed and is belonged in theory the category that digital signal processing is explained.Integrated interpretation method mainly contains direct interpretation procedure and indirect interpretation method.Direct method is by the raw data of radar detection being done to the processing of some routines, according to the power of radar appearance, phase characteristic, the characteristic informations such as the variation of lineups, again in conjunction with drilling data and other relevant geologic informations, directly reflected signal is made quantitatively and qualitative interpretation, but the geologic condition in the face of more complicated, be difficult to the correct external appearance characteristic of explaining after being complicated, and complicated and changeable due to underground medium, when often having, echoed signal becomes, the feature of non-stationary and randomness, so adopting the analysis and processing method of random signal analyzes data information, extract the proper vector of the self structure of surveying thing, choose principal ingredient wherein, and preserve complicated Nonlinear Mapping information with neural network model, realize accordingly automatic classification identification, this interpretation methods can be good at the principal character that thing is surveyed in reflection, and can realize automatic explanation, rapidly and efficiently.
Summary of the invention
The shortcoming that the object of the invention is to cannot accurately realize fast for geologic radar detection soil water-containing classification, provides a kind of geological radar technology to realize the technical method of soil water-containing Classification and Identification.
The technical scheme that the present invention solves its technical matters employing is:
Geological radar technology is realized the technical method of soil water-containing Classification and Identification, it is characterized in that comprising following concrete steps:
1) data pre-service:
The pretreated input data of data are the raw data that geological radar collects, and first Georadar Data are carried out to zero line correction, then carry out wavelet transformation, finally carry out low-pass filtering, and its output data are for completing pretreated data;
2) soil water-containing information extraction:
The step 1) of usining completes pre-service output data afterwards as the input data of this step, it is carried out to autoregressive moving average and ask power spectrum, carry out again data normalization, finally carry out principal component analysis (PCA), extract soil physical property infomation, deposit in sample database, it is soil water-containing characteristic that this step completes later output data;
3) neural network classification identification:
Using step 2) complete the soil water-containing characteristic of output after processing as the input data of this step, the data of depositing in sample database are inputted to neural metwork training neural network to be trained, then with the ripe network of training, complete the Classification and Identification of data to be sorted, the output data after this step completes are classification number under data to be identified;
4) result shows:
After the step 3) of usining completes and processes, under the data to be identified of output, classification number, as the input data of this step, carries out color range modulation to classification results, draws out search coverage classification of soils result figure.
Described zero line is proofreaied and correct: every track data that geological radar is collected carries out zero line correction, removes the signal drift noise of instrument self; Described wavelet transformation is: the data after zero line is proofreaied and correct are carried out wavelet transformation, filtering high frequency saltus step noise and extraneous high frequency interference noise again; Described low-pass filtering: the data after wavelet transformation are carried out to low-pass filtering again, near the energy direct current of filtered signal.
Described autoregressive moving average asks power spectrum to be: the every track data in Georadar Data after pre-service is carried out to the variation of autoregressive moving average power spectrum, try to achieve power spectrum information; Described data normalization is: power spectrum data is normalized; Described principal component analysis (PCA) is: the data after normalization are carried out to principal component analysis (PCA); Described extraction soil physical property infomation is: using the principal character vector after principal component analysis (PCA) as soil physical property infomation; Described sample database is: the soil physical property infomation obtaining is stored in property data base.
Described neural network training is: the characteristic input neural network in property data base is trained, until network training is ripe; Described neural network classification is identified as: the ripe neural network of radar data input training of the classification to be predicted after pretreated is carried out to Classification and Identification, provide affiliated classification number.
Described color range is modulated to: the classification results that step 3) is obtained carries out color range modulation, the corresponding color value of each class; Described drafting classification of soils result figure is: utilize the color range modulating that classification results one-tenth figure is shown, image can intuitively show search coverage soil water-containing classification schematic diagram.
The invention has the beneficial effects as follows: based on geological radar technology, a kind of technical method of soil water-containing Classification and Identification is being provided aspect soil water-containing Classification and Identification, compare with traditional boring with sampling method of testing, saved cost, realized the detection and identify on yardstick among a small circle, and by neural network, realized fast automatic identification, greatly improved efficiency.For actual detection provides a kind of technological means of precise and high efficiency.
Accompanying drawing explanation
Fig. 1 is that data is processed explanation FB(flow block).
Fig. 2 is geologic radar detection curve map, and wherein (a) is raw data curve; (b) for zero line is proofreaied and correct rear data and curves; (c) be data and curves after wavelet transformation; (d) be data and curves after low-pass filtering.
Fig. 3 is power spectrum conversion comparison diagram, the power spectrum chart that wherein (a) obtains for classical Fourier transform, the power spectrum chart (b) obtaining for autoregressive moving average power Spectral Estimation.
Fig. 4 is the corresponding power spectrum chart of different water cut, and wherein (a), for water percentage is the power spectrum chart of 20% correspondence, (b) for water percentage is the power spectrum chart of 30% correspondence, for water percentage, is (c) power spectrum chart of 40% correspondence.
Fig. 5 is geologic radar detection original section image.
Fig. 6 is the pretreated profile image of geological radar.
Fig. 7 is imaging color range table.
Fig. 8 is the classification sectional view of geologic radar detection soil water-containing.
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in further detail:
Fig. 1 is the data processing explanation FB(flow block) that geological radar technology of the present invention is realized soil water-containing Classification and Identification technical method.The concrete implementation detail of each step is as follows:
1. data pre-service
The pretreated object of data is before soil water-containing information extraction, to carry out the removal of various undesired signals, guarantees the reliability that soil water-containing response characteristic is extracted.Fig. 5 is the raw data profile image that geological radar gathers, and Fig. 6 is through pretreated Ground Penetrating Radar image, is specifically implemented as follows:
First the Georadar Data collecting is carried out to zero line correction, remove the signal drift noise of instrument self, Fig. 2 (a) original signal, Fig. 2 (b) is the signal after signal zero line is proofreaied and correct, the method of implementing is as follows: calculate signal average, by original signal value, deduct required average.
Secondly, the curve after zero line is proofreaied and correct carries out wavelet transformation, and compacting high frequency noise, uses Moret wavelet function, and scale parameter is 2, and transformation results is referring to Fig. 2 (c).
Finally, the signal after wavelet transformation is carried out to low-pass filtering, get rid of the interference of signal DC component, low-pass filtering parameter is 25MHz, and after filtering, result is referring to Fig. 2 (d).
2. soil physical property infomation extracts
First to pretreated Georadar Data signal
carrying out autoregressive moving average spectrum estimates.
for meeting the equation of difference
There is spectral density
, wherein
,
with
,
,
be called " rational expression spectral density ".Here adopt Cadzow Zymography in order to reduce number of parameters and to realize spectra calculation.Rational expression spectral density is rewritten into
, and be decomposed into
, wherein
, all the other
,
between pass be:
Stationary process power spectrum is written as:
Water cut has a great impact the dynamic characteristic of Electromagnetic Wave Propagation, on power spectrum, there is obvious performance, because water-bearing media electric conductivity strengthens, therefore medium electro-magnetic wave absorption coefficient strengthens, especially more obvious to the absorption of high-frequency signal, radio-frequency component occupation rate declines, and contrary low-frequency component occupancy volume increases, and can embody from the side media water-bearing rate feature.Utilize geological radar simulation tool (GPRMAX) to set up model herein, obtain the detection data of different water cut soil, reflection wave is analyzed, the (a) and (b) of Fig. 4, (c) are respectively that the water percentage that simulates is 20%, 30% and 40% power spectrum chart, and its low-frequency component accounting increases successively.
Secondly power spectrum information is normalized, first finds out maximin wherein, with treating that normalized numerical value deducts minimum value and obtains divided by the difference of maximin.Normalization can be eliminated the impact of data on neural network from dimension, reduces the fluctuation of network.
Then carry out principal component analysis (PCA), the data after normalization are extracted to major component.Specific implementation is: the sample matrix X to N dimension (experiment is 1024 dimensions), the covariance matrix S of first step compute matrix X; Second step calculates the latent vector e1 of covariance matrix S, e2 ... the eigenvalue of eN, eigenvalue, is chosen wherein latent vector corresponding to the larger eigenvalue of general impacts to little sequence by large, takes here overall accumulative total contribution degree is reached to 85% latent vector corresponding to eigenvalue; Among the space that the 3rd step data for projection is opened to eigenvector.Principal component analysis (PCA) has realized data compression effectively, kicks out of redundant data, at angle of statistics, holds the difference between signal, parses major influence factors from multiple factors.
Finally characteristic is deposited in to sample database storage administration.
3. neural network classification identification
First by the data input neural network training network in sample database, until network reaches frequency of training (setting 100000 times) or the systematic error of regulation, reach permissible value (setting 0.00001), think that network training is ripe.
Then by the ripe network of data input training of the search coverage to be identified of handling well, the result of record sort identification.
4. result shows
The classification results of output is 1,2,3,4,5,6 these 6 classes, corresponding soil moisture content is more than 0 ~ 10%, 10% ~ 20%, 20% ~ 30%, 30% ~ 40%, 40% ~ 50%, 50% this 6 kind respectively, stipulate the corresponding a kind of specific color of a classification number, modulation color range, as shown in Figure 7.
Then by the color range of modulation, become figure to show the Output rusults of search coverage, provide soil water-containing Classification and Identification sectional view, as shown in Figure 8.
Claims (5)
1. geological radar technology is realized the technical method of soil water-containing Classification and Identification, it is characterized in that comprising following concrete steps:
1) data pre-service
The pretreated input data of data are the raw data that geological radar collects, and first Georadar Data are carried out to zero line correction, then carry out wavelet transformation, finally carry out low-pass filtering, and its output data are for completing pretreated data;
2) soil water-containing information extraction
The step 1) of usining completes pre-service output data afterwards as the input data of this step, it is carried out to autoregressive moving average and ask power spectrum, carry out again data normalization, finally carry out principal component analysis (PCA), extract soil physical property infomation, deposit in sample database, it is soil water-containing characteristic that this step completes later output data;
3) neural network classification identification
Using step 2) complete the soil water-containing characteristic of output after processing as the input data of this step, the data of depositing in sample database are inputted to neural metwork training neural network to be trained, then with the ripe network of training, complete the Classification and Identification of data to be sorted, the output data after this step completes are classification number under data to be identified;
4) result shows
After the step 3) of usining completes and processes, under the data to be identified of output, classification number, as the input data of this step, carries out color range modulation to classification results, draws out search coverage classification of soils result figure.
2. geological radar technology according to claim 1 is realized the technical method of soil water-containing Classification and Identification, it is characterized in that, in step 1), described zero line is proofreaied and correct and is: every track data that geological radar is collected carries out zero line correction, removes the signal drift noise of instrument self; Described wavelet transformation is: the data after zero line is proofreaied and correct are carried out wavelet transformation, filtering high frequency saltus step noise and extraneous high frequency interference noise again; Described low-pass filtering: the data after wavelet transformation are carried out to low-pass filtering again, near the energy direct current of filtered signal.
3. geological radar technology according to claim 1 is realized the technical method of soil water-containing Classification and Identification, step 2) in, described autoregressive moving average asks power spectrum to be: the every track data in Georadar Data after pre-service is carried out to the variation of autoregressive moving average power spectrum, try to achieve power spectrum information; Described data normalization is: power spectrum data is normalized; Described principal component analysis (PCA) is: the data after normalization are carried out to principal component analysis (PCA); Described extraction soil physical property infomation is: using the principal character vector after principal component analysis (PCA) as soil physical property infomation; Described sample database is: the soil physical property infomation obtaining is stored in property data base.
4. geological radar technology according to claim 1 is realized the technical method of soil water-containing Classification and Identification, in step 3), described neural network training is: the characteristic input neural network in property data base is trained, until network training is ripe; Described neural network classification is identified as: the ripe neural network of radar data input training of the classification to be predicted after pretreated is carried out to Classification and Identification, provide affiliated classification number.
5. geological radar technology according to claim 1 is realized the technical method of soil water-containing Classification and Identification, and in step 4), described color range is modulated to: the classification results that step 3) is obtained carries out color range modulation, the corresponding color value of each class; Described drafting classification of soils result figure is: utilize the color range modulating that classification results one-tenth figure is shown, image can intuitively show search coverage soil water-containing classification schematic diagram.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410286867.8A CN104143115A (en) | 2014-06-25 | 2014-06-25 | Technological method for achieving soil water content classified identification through geological radar technology |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410286867.8A CN104143115A (en) | 2014-06-25 | 2014-06-25 | Technological method for achieving soil water content classified identification through geological radar technology |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104143115A true CN104143115A (en) | 2014-11-12 |
Family
ID=51852284
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410286867.8A Pending CN104143115A (en) | 2014-06-25 | 2014-06-25 | Technological method for achieving soil water content classified identification through geological radar technology |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104143115A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105628904A (en) * | 2015-12-22 | 2016-06-01 | 中国铁道科学研究院铁道建筑研究所 | Ground penetrating radar based water content detection method for railroad bed |
CN107632138A (en) * | 2017-09-14 | 2018-01-26 | 中国神华能源股份有限公司 | Calculate method, apparatus, medium and the system of sandy loam stratum porous media |
CN108038471A (en) * | 2017-12-27 | 2018-05-15 | 哈尔滨工程大学 | A kind of underwater sound communication signal type Identification method based on depth learning technology |
CN109164429A (en) * | 2018-10-17 | 2019-01-08 | 赣南师范大学 | Multi-fractal features Aircraft Targets classification method based on principal component analysis |
CN109670629A (en) * | 2018-11-16 | 2019-04-23 | 浙江蓝卓工业互联网信息技术有限公司 | Coal-burning boiler thermal efficiency forecast method based on shot and long term Memory Neural Networks |
CN111123263A (en) * | 2019-12-12 | 2020-05-08 | 创新奇智(合肥)科技有限公司 | Tunnel reinforcing steel bar identification and detection system and method based on geological radar detection data |
CN112085044A (en) * | 2019-06-14 | 2020-12-15 | 中南大学 | Slope dynamic classification method based on automatic monitoring data |
-
2014
- 2014-06-25 CN CN201410286867.8A patent/CN104143115A/en active Pending
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105628904A (en) * | 2015-12-22 | 2016-06-01 | 中国铁道科学研究院铁道建筑研究所 | Ground penetrating radar based water content detection method for railroad bed |
CN107632138A (en) * | 2017-09-14 | 2018-01-26 | 中国神华能源股份有限公司 | Calculate method, apparatus, medium and the system of sandy loam stratum porous media |
CN108038471A (en) * | 2017-12-27 | 2018-05-15 | 哈尔滨工程大学 | A kind of underwater sound communication signal type Identification method based on depth learning technology |
CN109164429A (en) * | 2018-10-17 | 2019-01-08 | 赣南师范大学 | Multi-fractal features Aircraft Targets classification method based on principal component analysis |
CN109164429B (en) * | 2018-10-17 | 2022-09-30 | 赣南师范大学 | Multi-fractal characteristic aircraft target classification method based on principal component analysis |
CN109670629A (en) * | 2018-11-16 | 2019-04-23 | 浙江蓝卓工业互联网信息技术有限公司 | Coal-burning boiler thermal efficiency forecast method based on shot and long term Memory Neural Networks |
CN109670629B (en) * | 2018-11-16 | 2021-09-07 | 浙江蓝卓工业互联网信息技术有限公司 | Coal-fired boiler thermal efficiency prediction method based on long-short term memory neural network |
CN112085044A (en) * | 2019-06-14 | 2020-12-15 | 中南大学 | Slope dynamic classification method based on automatic monitoring data |
CN112085044B (en) * | 2019-06-14 | 2023-11-24 | 中南大学 | Automatic monitoring data-based dynamic classification method for side slopes |
CN111123263A (en) * | 2019-12-12 | 2020-05-08 | 创新奇智(合肥)科技有限公司 | Tunnel reinforcing steel bar identification and detection system and method based on geological radar detection data |
CN111123263B (en) * | 2019-12-12 | 2022-04-19 | 创新奇智(浙江)科技有限公司 | Tunnel reinforcing steel bar identification and detection system and method based on geological radar detection data |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104143115A (en) | Technological method for achieving soil water content classified identification through geological radar technology | |
CN103941254A (en) | Soil physical property classification recognition method and device based on geological radar | |
CN105137498A (en) | Underground target detection and recognition system and method based on feature fusion | |
CN114564982B (en) | Automatic identification method for radar signal modulation type | |
Liu et al. | Deep learning and recognition of radar jamming based on CNN | |
CN111950503B (en) | Aviation transient electromagnetic data processing method and device and computing equipment | |
CN110110738A (en) | A kind of Recognition Method of Radar Emitters based on multi-feature fusion | |
CN104122588A (en) | Spectral decomposition based post-stack seismic data resolution ratio increasing method | |
CN103487832A (en) | Method for classifying supervised waveforms in three-dimensional seismic signal | |
CN104331583B (en) | A kind of Multifractal Modeling method based on Observed sea clutter | |
CN106249299A (en) | The lower thin layer weak reflected seismic energy restoration methods of strong reflection shielding and device | |
CN107179550B (en) | A kind of seismic signal zero phase deconvolution method of data-driven | |
CN107255831A (en) | A kind of extracting method of prestack frequency dispersion attribute | |
CN108710777A (en) | Abnormality recognition method is visited in the diversification that own coding neural network is accumulated based on multireel | |
CN103728662A (en) | Method for estimating stratum medium quality factors based on seismic signal envelope peak | |
CN102928875B (en) | Wavelet extraction method based on fractional number order Fourier | |
CN113887398A (en) | GPR signal denoising method based on variational modal decomposition and singular spectrum analysis | |
CN104268561A (en) | Hyperspectral image mixing eliminating method based on structure prior low rank representation | |
CN109655883A (en) | A kind of earthquake dividing method and system for target | |
CN107632326A (en) | Geophysical signal denoising method | |
CN109655904A (en) | A kind of multiple dimensioned solution cavity recognition methods and system | |
Bagheri et al. | Modeling the facies of reservoir using seismic data with missing attributes by dissimilarity based classification | |
Li et al. | Magnetotelluric noise suppression via convolutional neural network | |
CN107832798A (en) | Polarimetric SAR Image object detection method based on NSCT ladder pessimistic concurrency controls | |
Shang et al. | Automatic microseismic events detection using morphological multiscale top-hat transformation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20141112 |