CN106198745B - Anchoring defect length recognition methods based on reflected energy than parameter - Google Patents
Anchoring defect length recognition methods based on reflected energy than parameter Download PDFInfo
- Publication number
- CN106198745B CN106198745B CN201610526205.2A CN201610526205A CN106198745B CN 106198745 B CN106198745 B CN 106198745B CN 201610526205 A CN201610526205 A CN 201610526205A CN 106198745 B CN106198745 B CN 106198745B
- Authority
- CN
- China
- Prior art keywords
- anchoring
- frequency coefficient
- layer
- low frequency
- signal
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention discloses a kind of anchoring defect length recognition methods based on reflected energy than parameter, the following steps are included: S1: reading in detection signal, nondestructive testing signal is read in from the Detection of Bolt Bonding Integrity nondestructive detecting instrument based on stess wave method, the signal that waveform variation is consistent, nothing outranges feature is screened out from it and is analyzed for data；S2: Multiscale Wavelet Decomposition selects wavelet basis function, then selects Decomposition order, carries out one-dimensional discrete Multiscale Wavelet Decomposition to nondestructive testing signal；S3: the low frequency coefficient obtained after wavelet decomposition is compared with high frequency coefficient of the detection signal under different scale, the place with common mutation feature is found out, and then obtain anchoring starting end position and the corresponding extreme value of anchoring section bottom position by waveform comparison method.The present invention can accurately identify anchoring defect length, and method is simple, easy to use, at low cost.
Description
Technical field
The present invention relates to mine bolt anchorage technology fields more particularly to a kind of anchoring based on reflected energy than parameter to lack
Fall into length recognition methods.
Background technique
In the evaluation of mine bolt anchoring quality, limit anchoring force is an important index；Scholars pass through a large amount of
Research, has found the relationship between anchorage length and limit anchoring force, therefore, as long as having measured anchoring by non-destructive testing technology
Length, so that it may calculate limit anchoring force using formula；The acquisition methods of anchorage length are mainly by lossless inspection at present
At the time of identifying that anchoring starting point and anchoring section bottom end correspond in survey signal, then it is calculated using formula；But it is this
Method can not determine that anchoring section whether there is defect, and the presence of defect will lead to practical anchorage length and reduce, and therefore, utilize this
The anchorage length that kind method is calculated is inaccurate, correspondingly, can not by the limit anchoring force that the length computation comes out
Accurately reflect the anchor situation of anchor pole.
Summary of the invention
Technical problems based on background technology, the invention proposes longer than the anchoring defect of parameter based on reflected energy
Spend recognition methods.
Anchoring defect length recognition methods proposed by the present invention based on reflected energy than parameter, comprising the following steps:
S1: detection signal is read in, reads in nothing from the Detection of Bolt Bonding Integrity nondestructive detecting instrument based on stess wave method
Damage detection signal is screened out from it the signal that waveform variation is consistent, nothing outranges feature and analyzes for data；
S2: Multiscale Wavelet Decomposition selects wavelet basis function, then selects Decomposition order, carries out to nondestructive testing signal
One-dimensional discrete Multiscale Wavelet Decomposition；
S3: waveform comparison method, by the high frequency of the low frequency coefficient obtained after wavelet decomposition and detection signal under different scale
Coefficient compares, and the place with common mutation feature is found out, and then obtains anchoring starting end position and anchoring section bottom end
The corresponding extreme value in position；
S4: theoretical calculation utilizes formula according to the design length of anchor pole and anchoring section(lfIndicate anchor pole certainly
By segment length, vfIndicate the speed that stress wave is propagated in free section) calculate the corresponding stress wave of anchoring starting end position
Propagation time utilizes formula(lmIndicate the length of anchorage part, vmIndicate consolidation velocity of wave) it calculates
The anchoring section bottom position corresponding stress wave propagation time, and then acquisition anchoring starting end position and anchoring section bottom position are corresponding
Signal waveform extreme value；
S5: calculating reflected energy ratio, and anchoring is originated end position, the corresponding signal waveform extreme value number of anchoring section bottom position
According to substitution formula(AsFor the corresponding signal waveform extreme value in anchoring starting point, AeFor anchoring section bottom end pair
The signal waveform extreme value answered), calculate reflected energy ratio；
S6: judgement anchoring defect length: by reflected energy ratio reSubstitute into the defect length obtained by laboratory test
With the corresponding relationship formula between reflected energy ratio(A (r in formulae), a, b, c are related with rock property
Coefficient, need to demarcate by testing), calculate anchoring defect specific length ld。
It preferably, by two groups of filter equalizers is first layer high frequency coefficient and the by nondestructive testing signal in the S2
One layer of low frequency coefficient；It then is that second layer high frequency coefficient and the second layer are low by two groups of filter equalizers by first layer low frequency coefficient
Second layer low frequency coefficient is then third layer high frequency coefficient and third layer low frequency system by two groups of filter equalizers by frequency coefficient
Nondestructive testing signal is disassembled according to frequency height after 5 times are decomposed, is carrying out Multiscale Wavelet Decomposition to signal by number
The wavelet basis function of Shi Suoxuan is respectively db5 and db6.
Preferably, in the S3, waveform comparison method includes paddy peak-paddy peak correspondent method, Feng Feng-paddy paddy correspondent method, peak
Paddy-peak peak correspondent method and Feng Feng-paddy peak correspondent method, the waveform being compared are respectively that layer 5 low frequency coefficient and layer 5 are high
Frequency coefficient, layer 5 low frequency coefficient and the 4th layer of high frequency coefficient.
In the present invention, the lossless inspection of stress wave is only applicable to than the anchoring defect length recognition methods of parameter based on reflected energy
Survey technology, stress wave nondestructive testing technology are to excite audio-frequency stress wave in termination point of anchor rod and receive reflection wave signal, by adopting
The reflection signal collected is analyzed to make the anchoring quality of anchor pole a kind of non-destructive testing technology of evaluation, more by small echo
Scale Decomposition can filter out the interference of nondestructive testing signal high frequency components, can extract anchoring starting end position and anchoring section
The corresponding extreme value of bottom position, the present invention can accurately identify anchoring defect length, and method is simple, easy to use, at low cost.
Detailed description of the invention
Fig. 1 is the schematic diagram proposed by the present invention based on reflected energy than the anchoring defect length recognition methods of parameter.
Specific embodiment
Combined with specific embodiments below the present invention is made further to explain.
Embodiment
With reference to Fig. 1, the present embodiment proposes the anchoring defect length recognition methods based on reflected energy than parameter, including with
Lower step:
S1: detection signal is read in, reads in nothing from the Detection of Bolt Bonding Integrity nondestructive detecting instrument based on stess wave method
Damage detection signal is screened out from it the signal that waveform variation is consistent, nothing outranges feature and analyzes for data；
S2: Multiscale Wavelet Decomposition selects wavelet basis function, then selects Decomposition order, carries out to nondestructive testing signal
One-dimensional discrete Multiscale Wavelet Decomposition；
S3: waveform comparison method, by the high frequency of the low frequency coefficient obtained after wavelet decomposition and detection signal under different scale
Coefficient compares, and the place with common mutation feature is found out, and then obtains anchoring starting end position and anchoring section bottom end
The corresponding extreme value in position；
S4: theoretical calculation utilizes formula according to the design length of anchor pole and anchoring section(lfIndicate anchor pole certainly
By segment length, vfIndicate the speed that stress wave is propagated in free section) calculate the corresponding stress wave of anchoring starting end position
Propagation time utilizes formula(lmIndicate the length of anchorage part, vmIndicate consolidation velocity of wave) it calculates
The anchoring section bottom position corresponding stress wave propagation time, and then acquisition anchoring starting end position and anchoring section bottom position are corresponding
Signal waveform extreme value；
S5: calculating reflected energy ratio, and anchoring is originated end position, the corresponding signal waveform extreme value number of anchoring section bottom position
According to substitution formula(AsFor the corresponding signal waveform extreme value in anchoring starting point, AeFor anchoring section bottom end pair
The signal waveform extreme value answered), calculate reflected energy ratio；
S6: judgement anchoring defect length: by reflected energy ratio reSubstitute into the defect length obtained by laboratory test
With the corresponding relationship formula between reflected energy ratio(A (r in formulae), a, b, c are related with rock property
Coefficient, need to demarcate by testing), calculate anchoring defect specific length ld。
It by two groups of filter equalizers is first layer high frequency coefficient and the by nondestructive testing signal in S2 in the present embodiment
One layer of low frequency coefficient；It then is that second layer high frequency coefficient and the second layer are low by two groups of filter equalizers by first layer low frequency coefficient
Second layer low frequency coefficient is then third layer high frequency coefficient and third layer low frequency system by two groups of filter equalizers by frequency coefficient
Nondestructive testing signal is disassembled according to frequency height after 5 times are decomposed, is carrying out Multiscale Wavelet Decomposition to signal by number
The wavelet basis function of Shi Suoxuan is respectively db5 and db6, and in S3, waveform comparison method includes paddy peak-paddy peak correspondent method, Feng Feng-paddy
Paddy correspondent method, peak valley-peak peak correspondent method and Feng Feng-paddy peak correspondent method, the waveform being compared are respectively layer 5 low frequency system
Several and layer 5 high frequency coefficient, layer 5 low frequency coefficient and the 4th layer of high frequency coefficient.
In specific application it is to be noted that either utilizing waveform comparison method or numerical method or two
Kind of method combines, all it cannot be guaranteed that accurately identifying anchoring starting point and the corresponding position in anchoring section bottom end, and if this
The inaccuracy of two positions identification will lead to occur using judging result of the reflected energy ratio method to defect length very big
Deviation.Therefore, in order to guarantee the accuracy of judging result, in practical applications, need to take that " signal carries out on multiple scales
Identify → calculate separately reValue → reValue is averaged " method.Specific processing step is as follows:
1. respectively using db5 and db6 as wavelet basis detection signal decomposition to the 5th layer.Take db5-d4 layers therein, db5-
D5 layers, db6-d4 layers and db6-d5 layers totally four groups of decomposed signal.
2. using anchoring starting point described in S3, S4 and the recognition methods of anchoring section bottom position on four groups of decomposed signals
Position identification is carried out respectively.In view of anchoring the identification knot of starting point and anchoring section bottom end corresponding position in different component solution signal
Fruit has certain deviation, and in order to guarantee the accuracy of recognition result, every group of signal identifies two groups of knots according to recognition rule
Fruit.
3. every group of recognition result is substituted into formula respectively(AsFor the corresponding signal wave in anchoring starting point
Shape extreme value, AeFor the corresponding signal waveform extreme value in anchoring section bottom end) reflected energy ratio is calculated, obtain 8 results.8 are tied
Fruit is averaged to obtain reAverage value, will the value substitute into S6 in calculate.
It is more by carrying out small echo to nondestructive testing signal than the anchoring defect length recognition methods of parameter based on reflected energy
Scale Decomposition can filter out the interference of upper frequency ingredient in signal, divide the relatively low frequency content in signal
Analysis, the present invention can accurately identify anchoring defect length, and method is simple, easy to use, at low cost.
The foregoing is only a preferred embodiment of the present invention, but scope of protection of the present invention is not limited thereto,
Anyone skilled in the art in the technical scope disclosed by the present invention, according to the technique and scheme of the present invention and its
Inventive concept is subject to equivalent substitution or change, should be covered by the protection scope of the present invention.
Claims (3)
1. the anchoring defect length recognition methods based on reflected energy than parameter, which comprises the following steps:
S1: detection signal is read in, reads in lossless inspection from the Detection of Bolt Bonding Integrity nondestructive detecting instrument based on stess wave method
Signal is surveyed, the signal that waveform variation is consistent, nothing outranges feature is screened out from it and is analyzed for data；
S2: Multiscale Wavelet Decomposition selects wavelet basis function, then selects Decomposition order, carries out to nondestructive testing signal one-dimensional
Discrete wavelet multi-resolution decomposition；
S3: waveform comparison method, by the high frequency coefficient of the low frequency coefficient obtained after wavelet decomposition and detection signal under different scale
It compares, the place with common mutation feature is found out, and then obtain anchoring starting end position and anchoring section bottom position
Corresponding extreme value；
S4: theoretical calculation utilizes formula according to the design length of anchor pole and anchoring sectionCalculate anchoring starting
The end position corresponding stress wave propagation time, utilize formulaIt is corresponding to calculate anchoring section bottom position
The stress wave propagation time, and then obtain the anchoring starting end position and corresponding signal waveform extreme value of anchoring section bottom position,
Wherein, lfIndicate free section length, vfIndicate the speed that stress wave is propagated in free section；
lmIndicate the length of anchorage part, vmIndicate consolidation velocity of wave；
S5: calculating reflected energy ratio, and anchoring is originated end position, anchoring section bottom position corresponding signal waveform extreme value data generation
Enter formulaReflected energy ratio is calculated,
Wherein, AsFor the corresponding signal waveform extreme value in anchoring starting point, AeFor the corresponding signal waveform extreme value in anchoring section bottom end；
S6: judgement anchoring defect length: by reflected energy ratio reSubstitute into the defect length and reflection obtained by laboratory test
Corresponding relationship formula between energy ratioA (r in above-mentioned relation formulae), a, b, c are to have with rock property
The coefficient of pass needs to demarcate by testing, and calculates the specific length l of anchoring defectd。
2. the anchoring defect length recognition methods according to claim 1 based on reflected energy than parameter, which is characterized in that
In the S2, it is first layer high frequency coefficient and first layer low frequency coefficient that nondestructive testing signal, which is passed through two groups of filter equalizers,；So
It is afterwards second layer high frequency coefficient and second layer low frequency coefficient by two groups of filter equalizers by first layer low frequency coefficient, then by the
Two layers of low frequency coefficient are third layer high frequency coefficient and third layer low frequency coefficient by two groups of filter equalizers, by third layer low frequency system
Number is the 4th layer of high frequency coefficient and the 4th layer of low frequency coefficient by two groups of filter equalizers, and the 4th layer of low frequency coefficient is passed through two groups
Filter equalizer be the without layer high frequency coefficient and the without layer low frequency coefficient, after 5 times are decomposed, by nondestructive testing signal according to frequency
Rate height disassembles, and when carrying out Multiscale Wavelet Decomposition to signal, selected wavelet basis function is respectively db5 and db6.
3. the anchoring defect length recognition methods according to claim 1 based on reflected energy than parameter, which is characterized in that
In the S3, waveform comparison method includes paddy peak-paddy peak correspondent method, Feng Feng-paddy paddy correspondent method, peak valley-peak peak correspondent method and peak
Peak-paddy peak correspondent method, the waveform being compared are respectively layer 5 low frequency coefficient and layer 5 high frequency coefficient, layer 5 low frequency
Coefficient and the 4th layer of high frequency coefficient.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610526205.2A CN106198745B (en) | 2016-07-05 | 2016-07-05 | Anchoring defect length recognition methods based on reflected energy than parameter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610526205.2A CN106198745B (en) | 2016-07-05 | 2016-07-05 | Anchoring defect length recognition methods based on reflected energy than parameter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106198745A CN106198745A (en) | 2016-12-07 |
| CN106198745B true CN106198745B (en) | 2019-03-08 |
Family
ID=57465471
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610526205.2A Expired - Fee Related CN106198745B (en) | 2016-07-05 | 2016-07-05 | Anchoring defect length recognition methods based on reflected energy than parameter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106198745B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106770651A (en) * | 2016-11-29 | 2017-05-31 | 山西大同大学 | A kind of layered roof in coal roadway anchoring body Failure Model survey meter and control method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007322401A (en) * | 2006-06-05 | 2007-12-13 | Aoki Asunaro Kensetsu Kk | Method for evaluating soundness of anchor |
| CN101571515A (en) * | 2009-05-31 | 2009-11-04 | 中国矿业大学(北京) | Nondestructive detection method and device of anchor rod based on variable-frequency pseudo-random signal |
| JP2010203810A (en) * | 2009-02-27 | 2010-09-16 | Shimizu Corp | Method and device for non-destructive inspection of concrete structure, and anchor bolt |
| KR20110066498A (en) * | 2009-12-11 | 2011-06-17 | 재단법인 포항산업과학연구원 | Apparatus for detecting defect in steel cable |
| CN102608223A (en) * | 2012-03-06 | 2012-07-25 | 张维平 | Digital signal analysis method based on spectral energy ratio |
| CN102890117A (en) * | 2012-02-15 | 2013-01-23 | 长江水利委员会长江科学院 | Method and device for detecting anchoring quality of long anchor cable |
-
2016
- 2016-07-05 CN CN201610526205.2A patent/CN106198745B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007322401A (en) * | 2006-06-05 | 2007-12-13 | Aoki Asunaro Kensetsu Kk | Method for evaluating soundness of anchor |
| JP2010203810A (en) * | 2009-02-27 | 2010-09-16 | Shimizu Corp | Method and device for non-destructive inspection of concrete structure, and anchor bolt |
| CN101571515A (en) * | 2009-05-31 | 2009-11-04 | 中国矿业大学(北京) | Nondestructive detection method and device of anchor rod based on variable-frequency pseudo-random signal |
| KR20110066498A (en) * | 2009-12-11 | 2011-06-17 | 재단법인 포항산업과학연구원 | Apparatus for detecting defect in steel cable |
| CN102890117A (en) * | 2012-02-15 | 2013-01-23 | 长江水利委员会长江科学院 | Method and device for detecting anchoring quality of long anchor cable |
| CN102608223A (en) * | 2012-03-06 | 2012-07-25 | 张维平 | Digital signal analysis method based on spectral energy ratio |
Non-Patent Citations (2)
| Title |
|---|
| The Application of HHT Method Used in the Anchor Bolt Testing Signal;Linsong Wang et al;《The 2nd International Conference on Computer Application and System Modeling》;20121231;第1099-1102页 * |
| 多点布测下锚固缺陷诊断的小波多尺度分析;孙冰 等;《煤炭学报》;20140731;第39卷(第7期);第1385-1390页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106198745A (en) | 2016-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104180754B (en) | Inversion method for biophysical property of forest based on LiDAR comprehensive waveform model | |
| US7383730B2 (en) | Methods for determining velocity of a stress wave within a material and homogeneity of properties within the material | |
| CN102928514B (en) | Frequency characteristic-based nondestructive detection method of stress waves of wood | |
| Socco et al. | Feasibility of ultrasonic tomography for nondestructive testing of decay on living trees | |
| JP6763394B2 (en) | Recording medium for storing the soil quality judgment device, soil quality judgment method and program | |
| CN109764986A (en) | A kind of steel member plane stress detection method based on ultrasonic shear waves phase spectrum | |
| CN106885846B (en) | Trees defect detecting device and detection method | |
| CN102183226B (en) | Based on the anchor rod nondestructive testing method of Multi-source Information Fusion | |
| CN106156503B (en) | A kind of multi-scale entropy characterizing method of anchor system internal flaw distribution | |
| CN114878582B (en) | Defect detection and analysis method and system for special steel | |
| CN103424420A (en) | Fitting-based nuclear magnetic resonance signal processing method | |
| Schultes et al. | Laser Optical Plankton Counter and Zooscan intercomparison in tropical and subtropical marine ecosystems | |
| CN105241428A (en) | Water depth retrieval method with hyper-spectrum | |
| CN103091342B (en) | Method for carrying out CT scanning analysis processing upon rock core sample | |
| EP3951381B1 (en) | Non-destructive testing method and device for testing and distinguishing internal and external defects of steel wire rope | |
| CN110954033A (en) | Concrete crack depth detection method and system | |
| CN103713051A (en) | Modeling method for spherical fruit texture prediction model | |
| JP3919015B2 (en) | Method and apparatus for nondestructive inspection of trees using acoustic tomography | |
| CN106644879A (en) | Method and device for determining permeability contribution values of different pore components of core | |
| CN112578435A (en) | Rock ultrasonic testing first arrival picking method and system | |
| CN106198745B (en) | Anchoring defect length recognition methods based on reflected energy than parameter | |
| CN104483389A (en) | Source array method based detection method of grouting quality of bridge prestressed pipeline | |
| CN110824007B (en) | Tubular pile crack detection method and system | |
| CN106680375B (en) | For determining the Air Coupling supersonic detection method of the elasticity modulus of material | |
| CN106645401B (en) | A kind of damage reason location and reconstructing method and system based on the estimation of frequency wave number |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190308 Termination date: 20200705 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |