CN105716545A - Structural surface roughness coefficient evaluation method based on Dice similarity measurement - Google Patents

Abstract

A structural surface roughness coefficient evaluation method based on Dice similarity measurement comprises the steps of: respectively selecting high-pixel images of m standard contour curves of Barton, and extracting coordinate data of each point of the image contour curves with the accuracy being 1[mu], according to the extracted x and y coordinate data, calculating fluctuating angles, and realizing characteristic vectors of the standard contour curves by means of normalization of adjacent fluctuating angle frequencies; selecting a structural surface to be measured in a field spot, counting the distribution of all adjacent fluctuating angles on a tested contour curve, recording the frequency of [theta]I in each interval, and forming tested contour curve characteristic vectors; carrying out normalization on the tested contour curve characteristic vectors; and according to the Dice similarity measurement method, determining the similarity of the tested curve with m standard curves respectively, and making changes to a similarity measurement result according to a formula (5), wherein r*k=1, and the corresponding roughness coefficient range is the roughness coefficient value of the tested curve. The structural surface roughness coefficient evaluation method is good in accuracy.

Description

Structural plane roughness coefficient evaluation methodology based on Dice similarity measure

Technical field

The present invention relates to a kind of structural plane roughness coefficient evaluation methodology, specifically the present invention is based on the relief feature analysis of Barton nominal contour curve, similarity analysis by nominal contour curve Yu sample contour curve characteristic vector, achieve the quantitative assessment of rock structural plane roughness coefficient, the invention belongs to field of engineering technology.

Background technology

Structural plane is the important composition of rock mass and widely distributed, and it exists the integrity and the seriality that destroy rock self, and therefore the feature of structural plane often determines mechanical property and the deformational behavior of rock mass.Structural plane roughness is the key factor affecting structural face shear strength parameter, and therefore, the research carrying out structural plane roughness coefficient (JRC) has highly important theory significance and engineer applied value.

At present, Chinese scholars has carried out a series of research for the structural plane coefficient of roughness from empirical method, statistical method, fractal theory three broad aspect.In engineering practice, contrast by Barton judge that the method for structural plane roughness coefficient is most widely used at the standard curve that 1981 propose, and be written with the recommendation method of ISRM (ISRM).The ultimate principle of the method is as follows: Barton rule of thumb proposes 10 structural plane nominal contour curves, and the structural plane coefficient of roughness is divided into 0～2, and 2～4 ..., 18～20, totally 10 ranks；The structural fece sample of collection and nominal contour curve are contrasted by people, find out immediate nominal contour curve according to its coarse fluctuating character, with the roughness value as this structural fece sample of the roughness value corresponding to this nominal contour curve.But, this method being determined structural plane roughness coefficient by the contrast of simple experience is very subjective, and its test result often changes along with the experience of tester is different, and therefore the method has obvious limitation.Experimental result can be produced much impacts to make this method clear actually, Beer et al. has carried out a network surveying research in 2002 for different testers, research shows that the roughness value result of the identical structural plane that different people judges has very big diversity, and the assembly average of different personnel's test result can tend towards stability along with the increase of tester's quantity, time more than more than 50 people, assembly average just tends towards stability.Alameda-Hern á ndezP et al. also carried out similar research in 2014, the difference is that the schooling considering tester that he is detailed, result shows the impact of this method test subject subjective judgment by artificial micro-judgment structural plane roughness, and test result is inaccurate.

Summary of the invention

For the deficiency that the accuracy overcoming existing structure surface roughness factor evaluation method is poor, the present invention provides the structural plane roughness coefficient evaluation methodology based on Dice similarity measure that a kind of accuracy is good.

The technical solution adopted for the present invention to solve the technical problems is:

A kind of structural plane roughness coefficient evaluation methodology based on Dice similarity measure, described evaluation methodology comprises the steps:

1) choose the high pixel photo of the m bar nominal contour curve of Barton respectively, with μ for precision, picture contour curve each point is carried out the extraction of coordinate data；

2) according to the x extracted, y coordinate data, adjacent two point (x on contour curve are calculated according to equation below_i,y_i)、(x_i+1,y_i+1) undulation angle θ_i,

${\mathrm{\θ}}_i=\arctan \left(\sqrt{{\left(\frac{y_{i+1}-y_i}{x_{i+1}-x_i}\right)}^2}\right)=\arctan \left(\sqrt{{\left(\frac{y_{i+1}-y_i}{\mathrm{\μ}}\right)}^2}\right)---\left(1\right)$

3) to set the number of degrees as interval, all adjacent contours angle, θ on statistical standard contour curve_iThe regularity of distribution, record θ_iOccur in the frequency in each interval, constitute nominal contour curvilinear characteristic vector S^*；

4) characteristic vector of nominal contour curve, by the normalization of adjacent contours angle frequency is realized, its formula is

$s_i=\frac{s_i^{*}-S_{\min }^{*}}{S_{max}^{*}-S_{min}^{*}},---\left(2\right)$

Wherein, s_iFor standard curve characteristic vector S^*The element of normalized vector；For the frequency on each angle interval of nominal contour curve；WithIt it is the minimum and maximum frequency of adjacent contours angle；

6) structural plane of required mensuration is selected at scene in the wild, adopts profilograph along the contour curve on measurement surface, rendering architecture face, direction, then passes through off-the-air picture and scans and data extraction, it is thus achieved that precision is the structural plane profile altitude data of μ；

7) to set the number of degrees as interval, all adjacent contours angle, θ on statistical test contour curve_iThe regularity of distribution, record θ_iOccur in the frequency in each interval, constitute test profile curvilinear characteristic vector T；

8) test profile curvilinear characteristic vector T is normalized by following formula (3)

$t_i=\frac{t_i^{*}-S_{\min }^{*}}{S_{max}^{*}-S_{min}^{*}}.---\left(3\right)$

t_iElement for the normalized vector of test curve characteristic vector；The frequency on number of degrees interval is often set for test contour curve；

9) according to Dice Similarity Measures, determining the similarity of test curve and m bar standard curve respectively, its formula is:

$D(T,S)=\frac{2T\·S}{\left|\right|T|{|}_2^2+|\left|S\right|{|}_2^2}=\frac{2\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{t}_i{s}_i}{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{t}_i^2+\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{s_i}^2},---\left(4\right)$

Similarity measure result is denoted as v_k=D (S_k, T), k=1,2 ..., m；M is the bar number of standard curve, S_kFor the normalized vector of kth bar standard curve, T is the normalized vector of test curve.

10) similarity measure result is done following change,

$r_k=\frac{2v_k-v_{min}-v_{max}}{v_{\max }-v_{\min }}.---\left(5\right)$

r_kNormalized value for similarity measure result；V_minIt it is the minima of similarity measure result；V_maxIt is the maximum of similarity measure result, r_k=1 illustrate the coarse relief feature of this standard curve and test curve closest to, its corresponding roughness value is exactly the roughness value value of this test curve.

Beneficial effects of the present invention is mainly manifested in: accuracy is good.

Accompanying drawing explanation

Fig. 1 is the schematic diagram of Barton nominal contour curve adjacent contours angle Statistical Distribution, wherein, a () represents standard curve 1, b () represents standard curve 2, c () represents standard curve 3, d () represents standard curve 4, (e) represents standard curve 5；F () represents standard curve 6, (g) represents standard curve 7, and (h) represents standard curve 8, and (i) represents standard curve 9, and (j) represents standard curve 10.

Fig. 2 is one embodiment of the invention test structural plane contour curve schematic diagram, and wherein, the unit of abscissa and vertical coordinate is: cm.

Detailed description of the invention

Below in conjunction with accompanying drawing, the invention will be further described.

Seeing figures.1.and.2, a kind of structural plane roughness coefficient evaluation methodology based on Dice similarity measure, described evaluation methodology comprises the steps:

1) choose the high pixel photo of 10 (taking m=10) bar nominal contour curve of Barton respectively, with μ=0.5mm for precision, picture contour curve each point is carried out the extraction of coordinate data；

2) according to the x extracted, y coordinate data, adjacent two point (x on contour curve are calculated according to equation below_i,y_i)、(x_i+1,y_i+1) undulation angle θ_i,

${\mathrm{\θ}}_i=\arctan \left(\sqrt{{\left(\frac{y_{i+1}-y_i}{x_{i+1}-x_i}\right)}^2}\right)=\arctan \left(\sqrt{{\left(\frac{y_{i+1}-y_i}{\mathrm{\μ}}\right)}^2}\right)---\left(1\right)$

3) with 2 ° (take and set the number of degrees as 2 °) for interval, all adjacent contours angle, θ on statistical standard contour curve_iThe regularity of distribution, record θ_iOccur in the frequency in each interval, constitute nominal contour curvilinear characteristic vector S^*.Statistical result shows adjacent contours angle, θ on 10 standard curves_iThe scope occurred is within 0～50 °；

4) characteristic vector of nominal contour curve, it is possible to by the normalization of adjacent contours angle frequency is realized, its formula is

$s_i=\frac{s_i^{*}-S_{\min }^{*}}{S_{max}^{*}-S_{min}^{*}},---\left(2\right)$

s_iFor standard curve characteristic vector S^*Normalized vector；For the frequency on the every 2 ° of intervals of nominal contour curve；WithIt it is the minimum and maximum frequency of adjacent contours angle；Article 10, the normalization characteristic vector of standard curve is as shown in table 1；

Table 1

6) structural plane of required mensuration is selected at scene in the wild, adopts profilograph along the contour curve on measurement surface, rendering architecture face, direction, then passes through off-the-air picture and scans and data extraction, it is thus achieved that precision is the structural plane profile altitude data of 0.5mm；

7) with 2 ° for interval, all adjacent contours angle, θ on statistical test contour curve_iThe regularity of distribution, record θ_iOccur in the frequency in each interval, constitute test profile curvilinear characteristic vector T；

8) test profile curvilinear characteristic vector T is normalized by following formula

$t_i=\frac{t_i^{*}-S_{\min }^{*}}{S_{max}^{*}-S_{min}^{*}}.---\left(3\right)$

t_iNormalized vector for test curve characteristic vector；For the frequency on the test every 2 ° of intervals of contour curve；

9) according to Dice Similarity Measures, determining test curve and the similarity of 10 standard curves respectively, its formula is:

$D(T,S)=\frac{2T\·S}{\left|\right|T|{|}_2^2+|\left|S\right|{|}_2^2}=\frac{2\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{t}_i{s}_i}{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{t}_i^2+\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{s_i}^2},---\left(4\right)$

Wherein, s_iFor standard curve characteristic vector S^*The element of normalized vector；For the frequency on the every 2 ° of intervals of nominal contour curve；WithIt it is the minimum and maximum frequency of adjacent contours angle；

10) similarity measure result is done following change by the similarity for the ease of determining test structural plane and normal structure face,

$r_k=\frac{2v_k-v_{min}-v_{max}}{v_{\max }-v_{\min }}.---\left(5\right)$

r_kNormalized value for similarity measure result；V_minIt it is the minima of similarity measure result；V_maxIt it is the maximum of similarity measure result.R_k=1 illustrate the coarse relief feature of this standard curve and test curve closest to, its corresponding roughness value scope is exactly the roughness value value of this test curve.

Example: the evaluation methodology of a kind of structural plane roughness coefficient based on Dice Similarity Measures, process is as follows:

First, on-site measurement granite structure face in the wild, its numeral topography is as shown in 2；

Secondly, with 2 ° for interval, all adjacent contours angle, θ on this test contour curve are added up_iThe regularity of distribution, record θ_iOccur in the frequency in each interval, constitute test profile curvilinear characteristic vector T, and after being normalized, obtain the normalization characteristic vector [0.00,0.26 of this test contour curve, 1.00,0.09,0.60,0.54,0.17,0.17,0.09,0.11,0.06,0.03,0.03,0.06,0.06,0.03,0.03,0.00,0.03,0.00,0.00,0.00,0.00,0.00,0.06].

Then, adopt Dice Similarity Measures to seek similarity with the normalization characteristic vector of table 1 standard curve, and it is as shown in table 2 to obtain the similarity after normalization.

Table 2

Finally, finding according to result of calculation, the Similarity value after the normalization of standard curve 6 (roughness value 10～12) and test curve is 1, and therefore the roughness value of this test curve is 10～12.

The core of the present invention is in that according to after Barton nominal contour curve and test profile curvilinear characteristic, find according to result of calculation, Similarity value after the normalization of standard curve 6 (roughness value 10～12) and test curve is 1, and therefore the roughness value of this test curve is 10～12.

The core of the present invention is in that the similarity quantification according to Barton nominal contour curve with test profile curvilinear characteristic vector determines the roughness value of test object, this method eliminates the impact of personal error, continuity simultaneously employs the short-cut method judging structural face shear strength based on Barton standard curve, when not losing simplicity, accuracy, it is achieved that the accurate evaluation of structural plane roughness coefficient.If changing in the scope that these characteristic vectors belong to the claims in the present invention and equivalent technologies thereof with the similar method for evaluating similarity of employing, then the present invention is also intended to comprise these changes and deformation.

Claims (1)

1. the structural plane roughness coefficient evaluation methodology based on Dice similarity measure, it is characterised in that: described evaluation methodology comprises the steps:

1) choose the high pixel photo of the m bar nominal contour curve of Barton respectively, with μ for precision, picture contour curve each point is carried out the extraction of coordinate data；

2) according to the x extracted, y coordinate data, adjacent two point (x on contour curve are calculated according to equation below_i,y_i)、(x_i+1,y_i+1) undulation angle θ_i,

${\mathrm{\θ}}_i=\arctan \left(\sqrt{{\left(\frac{y_{i+1}-y_i}{x_{i+1}-x_i}\right)}^2}\right)=\arctan \left(\sqrt{{\left(\frac{y_{i+1}-y_i}{\mathrm{\μ}}\right)}^2}\right)---\left(1\right)$

3) to set the number of degrees as interval, all adjacent contours angle, θ on statistical standard contour curve_iThe regularity of distribution, record θ_iOccur in the frequency in each interval, constitute nominal contour curvilinear characteristic vector S^*；

4) characteristic vector of nominal contour curve, by the normalization of adjacent contours angle frequency is realized, its formula is

$s_i=\frac{s_i^{*}-S_{\min }^{*}}{S_{max}^{*}-S_{min}^{*}},---\left(2\right)$

Wherein, s_iFor standard curve characteristic vector S^*The element of normalized vector；For the frequency on each angle interval of nominal contour curve；WithIt it is the minimum and maximum frequency of adjacent contours angle；

6) structural plane of required mensuration is selected at scene in the wild, adopts profilograph along the contour curve on measurement surface, rendering architecture face, direction, then passes through off-the-air picture and scans and data extraction, it is thus achieved that precision is the structural plane profile altitude data of μ；

7) to set the number of degrees as interval, all adjacent contours angle, θ on statistical test contour curve_iThe regularity of distribution, record θ_iOccur in the frequency in each interval, constitute test profile curvilinear characteristic vector T；

8) test profile curvilinear characteristic vector T is normalized by following formula (3)

$t_i=\frac{t_i^{*}-S_{\min }^{*}}{S_{max}^{*}-S_{min}^{*}}.---\left(3\right)$

t_iElement for the normalized vector of test curve characteristic vector；The frequency on number of degrees interval is often set for test contour curve；

9) according to Dice Similarity Measures, determining the similarity of test curve and m bar standard curve respectively, its formula is:

$D(T,S)=\frac{2T\·S}{\left|\right|T|{|}_2^2+|\left|S\right|{|}_2^2}=\frac{2\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{t}_i{s}_i}{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{t_i}^2+\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{s_i}^2},---\left(4\right)$

Similarity measure result is denoted as v_k=D (T, S_k), k=1,2 ..., m；M is the bar number of standard curve, S_kFor the normalized vector of kth bar standard curve, T is the normalized vector of test curve.

10) similarity measure result is done following change,

$r_k=\frac{2v_k-v_{min}-v_{max}}{v_{\max }-v_{\min }}.---\left(5\right)$

r_kNormalized value for similarity measure result；V_minIt it is the minima of similarity measure result；V_maxIt is the maximum of similarity measure result, r_k=1 illustrate the coarse relief feature of this standard curve and test curve closest to, its corresponding roughness value is exactly the roughness value value of this test curve.