CN110346791A - A kind of search of multichannel Ground Penetrating Radar layer position and optimization algorithm - Google Patents
A kind of search of multichannel Ground Penetrating Radar layer position and optimization algorithm Download PDFInfo
- Publication number
- CN110346791A CN110346791A CN201910785175.0A CN201910785175A CN110346791A CN 110346791 A CN110346791 A CN 110346791A CN 201910785175 A CN201910785175 A CN 201910785175A CN 110346791 A CN110346791 A CN 110346791A
- Authority
- CN
- China
- Prior art keywords
- consecutive points
- layer position
- position information
- layer
- site
- 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
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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a kind of search of multichannel Ground Penetrating Radar layer position and optimization algorithms, comprising the following steps: establishes the three-dimensional system of coordinate in region to be searched, and carries out collection point definition；The layer position information in a reference layer site is inputted in region to be searched, and marking reference layer site is A；On the basis of the A of reference layer site, the layer position information of consecutive points is acquired using amplitude tolerance threshold value and peak valley along X-axis and Y-axis external diffusion respectively；On the basis of consecutive points in the hope of layer position information, the layer position index that said three-dimensional body information C-Scan is acquired using amplitude tolerance threshold value and peak valley is repeated.Through the above scheme, the present invention has many advantages, such as that logic is easy, accurate and reliable, has very high practical value and promotional value in multichannel wavelet transform technical field.
Description
Technical field
The present invention relates to multichannel wavelet transform technical field, especially a kind of multichannel Ground Penetrating Radar layer position
Search and optimization algorithm.
Background technique
Ground Penetrating Radar, is to emit frequency electromagnetic waves to underground by transmitting antenna, is reflected back ground by receiving antenna reception
The electromagnetic wave in face, electromagnetic wave reflect when encountering the interface there are electrical property difference when propagating in underground medium, according to connecing
The features such as waveform, oscillator intensity and the variation of time of the electromagnetic wave received infer the spatial position of underground medium, structure, form
And buried depth.Traditional detection processing data are complicated, and carry out independent place to the data of detection acquisition after needing to detect one by one
There is a possibility that error in reason, any point.Therefore, it is necessary to the data to independent process to carry out later period optimization processing.
Summary of the invention
In view of the above-mentioned problems, the purpose of the present invention is to provide a kind of search of multichannel Ground Penetrating Radar layer position and optimization to calculate
Method, The technical solution adopted by the invention is as follows:
A kind of search of multichannel Ground Penetrating Radar layer position and optimization algorithm, comprising the following steps:
Step S1 establishes the three-dimensional system of coordinate in region to be searched, and carries out collection point definition.
Step S2, inputs the layer position information in a reference layer site in region to be searched, and marking reference layer site is A；
Layer position information is layer bit index value.
Step S3 uses amplitude tolerance threshold value and peak along X-axis and Y-axis external diffusion respectively on the basis of the A of reference layer site
Paddy acquires the layer position information of consecutive points.
Step S4 on the basis of the consecutive points in the hope of layer position information, repeats to acquire three using amplitude tolerance threshold value and peak valley
Tie up the layer position index of body information C-Scan.
Further, in the step S1, it is acquired a definition, comprising the following steps:
Step S11, label Z-direction are A-Scan；The index value range of the A-Scan be more than or equal to 0 and be less than etc.
In 255 positive integer；
Step S12, several A-Scan form a B-scan along the y axis；The Bscan composition one of n channel radar
Batch；Several Batch form a said three-dimensional body information C-Scan along the x axis；The n is the positive integer more than or equal to 1；
Preferably, in the step S3, the layer position information of consecutive points is acquired using amplitude tolerance threshold value and peak valley, it is specific to walk
It is rapid as follows:
Step S31, using the consecutive points A2 [i] of traversal retrieval reference layer site A；The i is greater than 0 and less than 255
Positive integer；
Step S32, it is preset whether the tolerance of the amplitude and preset reference value iRef that judge consecutive points A2 [i] is less than
Amplitude tolerance threshold value threshold；If more than, then return step S31, the consecutive points next point of traversal search consecutive points A2 [i]
A2[i+1]；Otherwise S33 is entered step；
Step S33 judges the kurtosis of consecutive points A2 [i], if the amplitude slope dA2 [i]≤0 of consecutive points A2 [i]
And dA2 [i+1] >=0, then consecutive points A2 [i] is wave crest, and enters step S34；If the amplitude slope dA2 of consecutive points A2 [i]
[i] >=0 and dA2 [i+1]≤0, then consecutive points A2 [i] is trough, and enters step S34；Otherwise return step S31, traversal
Search for the consecutive points next point A2 [i+1] of consecutive points A2 [i]；
Step S34, judge consecutive points A2 [i] for the range difference of the point of wave crest or trough layer position information corresponding to this point,
If the reference layer position information of a upper point of consecutive points A2 [i] be equal to -1, stop the retrieval of consecutive points A2 [i], otherwise into
Enter step S35；
Step S35 judges the layer potential difference of reference layer site A and consecutive points A2 [i] whether less than 5, then by consecutive points A2 [i]
Layer position iLast be labeled as i, and reference value iRef is modified as to the amplitude of consecutive points A2 [i].
It further, further include the transmitting of reference layer site A in the step S3, if upper one of consecutive points A2 [i]
The reference layer position information of point is not equal to -1, then consecutive points A2 [i] is used as reference layer site A, continues to test consecutive points A2's [i]
The layer position information of consecutive points.
Compared with prior art, the invention has the following advantages:
The present invention dexterously passes through the layer position information an of reference point, obtains the radar data of multichannel Ground Penetrating Radar acquisition
Entire three-dimension layer position data can simplify logic step without carrying out later period optimization to the data of detection and improve and examine
The accuracy of survey.The present invention acquires consecutive points by amplitude tolerance threshold value and the successional constraint information of peak valley and adjacent layer position
Layer position information, judge more accurate.In conclusion the present invention has many advantages, such as that high-efficient, logic is easy, accurate and reliable,
There is very high practical value and promotional value in multichannel wavelet transform technical field.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to the attached drawing used required in embodiment
It is briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as to protection
The restriction of range to those skilled in the art without creative efforts, can also be attached according to these
Figure obtains other relevant attached drawings.
Fig. 1 is logical flow chart of the invention.
Fig. 2 is some region of test result figure of the invention.
Specific embodiment
To keep the purposes, technical schemes and advantages of the application apparent, with reference to the accompanying drawings and examples to the present invention
It is described further, embodiments of the present invention include but is not limited to the following example.Based on the embodiment in the application, ability
Domain those of ordinary skill every other embodiment obtained without making creative work, belongs to the application
The range of protection.
Embodiment
As shown in Figure 1 to Figure 2, a kind of search of multichannel Ground Penetrating Radar layer position and optimization algorithm, packet are present embodiments provided
Include following steps:
The first step establishes the three-dimensional system of coordinate in region to be searched, and carries out collection point definition.Specifically:
(1) label Z-direction is A-Scan；Each A-Scan forms a waveform, i.e. index value range by 256 points
[0,255]。
(2) a series of A-Scan in the same direction y form a B-Scan, are labeled as A [i]；14 channel radar B-
Scan forms a Batch；C-Scan represents whole three-dimensional informations in the region of acquisition.
Second step, inputs the layer position information in a reference layer site in region to be searched, and marking reference layer site is A；
Layer position information is layer bit index value.Wherein, iRefDepth indicates the reference layer position of A point, and iRefValue indicates the reference of A point
Peak value (wave crest or trough)；ILast represents the reference layer position information of a point on A point；If last point does not retrieve layer position
Information remembers iLast=-1；If last point detection channel layer position information, iLast=iLastDepth；ILastDepth is
The layer place value of last point.StP, endP: the entire position C-Scan layers of initial search point index value, confining layers position fluctuation range are indicated.
Such as in stP=30, endP=50, represents and carry out the search of layer position within the scope of 30-50 index value.Amplitude tolerance threshold value
Threshold is preset as 6000.
Third step uses amplitude tolerance threshold value and peak along X-axis and Y-axis external diffusion respectively on the basis of the A of reference layer site
Paddy acquires the layer position information of consecutive points.
(31) using the consecutive points A2 [i] of traversal retrieval reference layer site A；The i is just whole greater than 0 and less than 255
Number；
(32) whether the tolerance of the amplitude and preset reference value iRef that judge consecutive points A2 [i] is less than preset amplitude
Tolerance threshold value threshold；If more than, then return step (31), the consecutive points next point A2 [i of traversal search consecutive points A2 [i]
+1]；Otherwise (33) are entered step；
(33) judge the kurtosis of consecutive points A2 [i], if the amplitude slope dA2 [i]≤0 of consecutive points A2 [i] and
DA2 [i+1] >=0, then consecutive points A2 [i] is wave crest, and enters step (34)；If the amplitude slope dA2 of consecutive points A2 [i]
[i] >=0 and dA2 [i+1]≤0, then consecutive points A2 [i] is trough, and enters step (34)；Otherwise return step (31), time
Go through the consecutive points next point A2 [i+1] of search consecutive points A2 [i]；
(31) consecutive points A2 [i] is judged for the range difference of the point of wave crest or trough layer position information corresponding to this point, if phase
The reference layer position information of a upper point of adjoint point A2 [i] is equal to -1, then stops the retrieval of consecutive points A2 [i], otherwise enters step
Suddenly (35)；
(35), the layer potential difference of reference layer site A and consecutive points A2 [i] are judged whether less than 5, then by consecutive points A2's [i]
Layer position iLast is labeled as i, and reference value iRef is modified as to the amplitude of consecutive points A2 [i].
(36), need to detect the layer position information of next point after the layer position information for acquiring consecutive points, therefore, it is necessary to carry out
The transmitting of reference layer site A, if the reference layer position information of a upper point of consecutive points A2 [i] is not equal to -1, by consecutive points A2
[i] is used as reference layer site A, continues to test the layer position information of the consecutive points of consecutive points A2 [i].
4th step on the basis of the consecutive points in the hope of layer position information, repeats step (31)~(36) and uses amplitude tolerance threshold
Value and peak valley acquire the layer position index of said three-dimensional body information C-Scan.
As soon as the radar of multichannel Ground Penetrating Radar acquisition can be obtained in conclusion the present invention only passes through a layer potential reference point
Entire C-Scan layers of position data of data.It is detected compared to single-pass channel layer position, greatly improves layer position search efficiency and detection essence
Degree.Compared with prior art, there is substantive distinguishing features outstanding and significant progress, in multichannel wavelet transform skill
Art field has very high practical value and promotional value.
Above-described embodiment is merely a preferred embodiment of the present invention, and it is not intended to limit the protection scope of the present invention, as long as using
Design principle of the invention, and the non-creative variation worked and made is carried out on this basis, it should belong to of the invention
Within protection scope.
Claims (4)
1. a kind of multichannel Ground Penetrating Radar layer position search and optimization algorithm, which comprises the following steps:
Step S1 establishes the three-dimensional system of coordinate in region to be searched, and carries out collection point definition；
Step S2, inputs the layer position information in a reference layer site in region to be searched, and marking reference layer site is A；It is described
Layer position information is layer bit index value；
Step S3 is asked along X-axis and Y-axis external diffusion using amplitude tolerance threshold value and peak valley respectively on the basis of the A of reference layer site
Obtain the layer position information of consecutive points；
Step S4 on the basis of the consecutive points in the hope of layer position information, repeats to acquire said three-dimensional body using amplitude tolerance threshold value and peak valley
The layer position index of information C-Scan.
2. a kind of multichannel Ground Penetrating Radar layer position search according to claim 1 and optimization algorithm, which is characterized in that described
In step S1, it is acquired a definition, comprising the following steps:
Step S11, label Z-direction are A-Scan；The index value range of the A-Scan is more than or equal to 0 and is less than or equal to 255
Positive integer；
Step S12, several A-Scan form a B-scan along the y axis；The Bscan of n channel radar forms a Batch；Number
A Batch forms a said three-dimensional body information C-Scan along the x axis；The n is the positive integer more than or equal to 1.
3. a kind of multichannel Ground Penetrating Radar layer position search according to claim 2 and optimization algorithm, which is characterized in that described
In step S3, the layer position information of consecutive points is acquired using amplitude tolerance threshold value and peak valley, the specific steps are as follows:
Step S31, using the consecutive points A2 [i] of traversal retrieval reference layer site A；The i is just whole greater than 0 and less than 255
Number；
Whether the tolerance of step S32, the amplitude for judging consecutive points A2 [i] and preset reference value iRef are less than preset amplitude
Tolerance threshold value threshold；If more than then return step S31, the consecutive points next point A2 [i+ of traversal search consecutive points A2 [i]
1]；Otherwise S33 is entered step；
Step S33 judges the kurtosis of consecutive points A2 [i], if the amplitude slope dA2 [i]≤0 of consecutive points A2 [i] and
DA2 [i+1] >=0, then consecutive points A2 [i] is wave crest, and enters step S34；If the amplitude slope dA2 of consecutive points A2 [i]
[i] >=0 and dA2 [i+1]≤0, then consecutive points A2 [i] is trough, and enters step S34；Otherwise return step S31, traversal
Search for the consecutive points next point A2 [i+1] of consecutive points A2 [i]；
Step S34 judges consecutive points A2 [i] for the range difference of the point of wave crest or trough layer position information corresponding to this point, if phase
The reference layer position information of a upper point of adjoint point A2 [i] is equal to -1, then stops the retrieval of consecutive points A2 [i], otherwise enters step
Rapid S35；
Step S35 judges the layer potential difference of reference layer site A and consecutive points A2 [i] whether less than 5, then by the layer of consecutive points A2 [i]
Position iLast is labeled as i, and reference value iRef is modified as to the amplitude of consecutive points A2 [i].
4. a kind of multichannel Ground Penetrating Radar layer position search according to claim 3 and optimization algorithm, which is characterized in that described
Further include the transmitting of reference layer site A in step S3, if the reference layer position information of a upper point of consecutive points A2 [i] not equal to-
1, then consecutive points A2 [i] is used as reference layer site A, continues to test the layer position information of the consecutive points of consecutive points A2 [i].
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910785175.0A CN110346791A (en) | 2019-08-23 | 2019-08-23 | A kind of search of multichannel Ground Penetrating Radar layer position and optimization algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910785175.0A CN110346791A (en) | 2019-08-23 | 2019-08-23 | A kind of search of multichannel Ground Penetrating Radar layer position and optimization algorithm |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110346791A true CN110346791A (en) | 2019-10-18 |
Family
ID=68181181
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910785175.0A Pending CN110346791A (en) | 2019-08-23 | 2019-08-23 | A kind of search of multichannel Ground Penetrating Radar layer position and optimization algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110346791A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111551927A (en) * | 2020-05-19 | 2020-08-18 | 上海圭目机器人有限公司 | Underground pipeline diameter measuring method based on three-dimensional ground penetrating radar |
| CN113253264A (en) * | 2021-05-10 | 2021-08-13 | 西安西测测试技术股份有限公司 | Ground penetrating radar data acquisition and recombination method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1272870A1 (en) * | 1999-12-08 | 2003-01-08 | STN ATLAS Marine Electronics GmbH | Method for determining the depth values of a body of water |
| US20080185903A1 (en) * | 2006-06-29 | 2008-08-07 | Igor Bausov | Look-ahead radar and horizon sensing for coal cutting drums and horizontal directional drills |
| CN102298160A (en) * | 2011-05-18 | 2011-12-28 | 西安交通大学 | Method for tracking data layer position of ground penetrating radar |
| CN104181596A (en) * | 2014-08-27 | 2014-12-03 | 中国石油集团东方地球物理勘探有限责任公司 | Geologic horizon automatic tracking method and device |
| CN105974407A (en) * | 2016-07-26 | 2016-09-28 | 中国科学院电子学研究所 | Ground penetrating radar underground horizon detection method |
| US20180335517A1 (en) * | 2016-11-10 | 2018-11-22 | China University Of Mining & Technology, Beijing | Method and Device for Detecting Discontinuous Body with Ground Penetrating Radar |
-
2019
- 2019-08-23 CN CN201910785175.0A patent/CN110346791A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1272870A1 (en) * | 1999-12-08 | 2003-01-08 | STN ATLAS Marine Electronics GmbH | Method for determining the depth values of a body of water |
| US20080185903A1 (en) * | 2006-06-29 | 2008-08-07 | Igor Bausov | Look-ahead radar and horizon sensing for coal cutting drums and horizontal directional drills |
| CN102298160A (en) * | 2011-05-18 | 2011-12-28 | 西安交通大学 | Method for tracking data layer position of ground penetrating radar |
| CN104181596A (en) * | 2014-08-27 | 2014-12-03 | 中国石油集团东方地球物理勘探有限责任公司 | Geologic horizon automatic tracking method and device |
| CN105974407A (en) * | 2016-07-26 | 2016-09-28 | 中国科学院电子学研究所 | Ground penetrating radar underground horizon detection method |
| US20180335517A1 (en) * | 2016-11-10 | 2018-11-22 | China University Of Mining & Technology, Beijing | Method and Device for Detecting Discontinuous Body with Ground Penetrating Radar |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111551927A (en) * | 2020-05-19 | 2020-08-18 | 上海圭目机器人有限公司 | Underground pipeline diameter measuring method based on three-dimensional ground penetrating radar |
| CN111551927B (en) * | 2020-05-19 | 2023-01-31 | 上海圭目机器人有限公司 | Underground pipeline diameter measuring method based on three-dimensional ground penetrating radar |
| CN113253264A (en) * | 2021-05-10 | 2021-08-13 | 西安西测测试技术股份有限公司 | Ground penetrating radar data acquisition and recombination method |
| CN113253264B (en) * | 2021-05-10 | 2023-09-22 | 西安西测测试技术股份有限公司 | Ground penetrating radar data acquisition and recombination method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6823266B2 (en) | Method for performing object-based connectivity analysis in 3-D seismic data volumes | |
| Liu et al. | Numerical simulation of subsurface radar for detecting buried pipes | |
| CN104656157B (en) | A kind of method and device in identification shale gas dessert area | |
| CN110346791A (en) | A kind of search of multichannel Ground Penetrating Radar layer position and optimization algorithm | |
| CN104422959A (en) | Method for detecting curvature property of reservoir boundary | |
| CA2464429A1 (en) | Global classification of sonic logs | |
| CN104570078B (en) | A kind of solution cavity detection method of the similitude cross directional variations rate based on frequency domain inclination angle | |
| CN103901467A (en) | Method for tracking positions of three-dimensional seismic data | |
| CN103163555B (en) | Middle-shallow Buried Gases gas pool identification method | |
| CN105242311A (en) | Method and apparatus for predicting time window of first arrival pickup | |
| CN111158050B (en) | Data acquisition system and method and tunnel seismic wave advanced prediction method | |
| CN105259578A (en) | Formation velocity determining method based on vertical seismic profile data | |
| CN112114358B (en) | Underground volcanic channel identification method based on three-dimensional seismic data representation | |
| CN111474592A (en) | Interwell electromagnetic detection system and method | |
| CN101452082A (en) | First arrival picking -up method for fractal seismic waves | |
| Yao et al. | Using seismic methods to detect connectivity of fracture networks controlled by strike-slip faults in ultra-deep carbonate reservoirs: A case study in northern tarim basin, China | |
| AU2002312550B2 (en) | Method for performing object-based connectivity analysis in 3-D seismic data volumes | |
| Wang et al. | Characteristic parameters extraction method of hidden Karst Cave from Borehole radar signal | |
| AU2002312550A1 (en) | Method for performing object-based connectivity analysis in 3-D seismic data volumes | |
| Chen et al. | Intelligent data fitting technique for 3D velocity reconstruction | |
| CN107015271B (en) | Surface structure survey method and apparatus | |
| CN110579798A (en) | Seismic acquisition observation method and system with equal reflection angle intervals | |
| CN108802826A (en) | Method, device and system for evaluating development condition of abnormal body of slotted hole in formation beside well | |
| CN114859414B (en) | Method and device for automatically extracting stratum dip angle information from seismic data | |
| CN109143354B (en) | Method and device for decomposing seismic waveform characteristics |
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: 20191018 |