Method and system for estimating volume stone content based on drill hole linear measurement and considering three-dimensional shape of block stone

Abstract

The invention provides a method and system for estimating the volume stone content based on drilling linear measurement by considering the three-dimensional shape of a block stone, and the method comprises the steps: analyzing the three-dimensional shape characteristics of the block stone through a CT scanning and digital image processing technology, and building a block stone three-dimensional model database containing the three-dimensional shape characteristic factors of the block stone; according to the maximum block stone particle size and the linear stone content obtained through field drilling, three-dimensional block stones in the database are randomly selected to generate a three-dimensional virtual model of the earth-rock aggregate; randomly generating a plurality of scanning lines in the three-dimensional virtual model, and calculating the relationship between the uncertainty factor of the soil-rock mixture and the sampling length ratio according to the linear stone content measured by the scanning lines; according to the on-site drilling sampling length ratio and the drilling linear stone content, the value of the uncertainty factor is determined, and the on-site volume stone content is estimated. And a scientific basis is provided for improving the safety and economic benefits of related engineering design of the earth-rock aggregate and reasonably selecting an excavation design and a supporting method.

Description

Technical Field

[0001] The present invention belongs to the field of geological engineering technology, and in particular relates to a method and system for estimating volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the block stone. Background Art

[0002] Soil-rock mixtures are widely distributed in natural geological bodies and human engineering activities, such as landslides, slope engineering, tunnel engineering, earth-rock dams, etc. A large number of studies in the past thirty years have shown that the structure and mechanical properties of soil-rock mixtures are mainly controlled by their volume rock content. VBP , accurately estimate the soil-rock mixture VBP It has the following important significance for related projects: ① It can accurately evaluate the shear strength of soil-rock mixture, effectively improve the safety and economic benefits of engineering design; ② It can provide a scientific basis for the reasonable selection of excavation equipment and support methods in related projects. VBP The prediction mainly relies on two methods: on-site screening and borehole catalog measurement. However, the above two methods have the following defects: ① On-site screening can only dig out the soil-rock mixture on the surface of the landslide for screening analysis, and the obtained stone content VBP It can only reflect the surface of the landslide VBP The accumulation layer of soil-rock mixture in landslides may be tens or even hundreds of meters thick. VBP As the burial depth increases, it is bound to change, and the relevant support engineering often passes through the entire accumulation layer, so it is necessary to master the entire soil-rock mixture accumulation layer. VBP It is difficult to achieve this only by screening. ② Since the diameter of the general drilling hole is only a dozen centimeters, for the entire landslide, the rock content obtained by drilling is usually considered to be a one-dimensional linear rock content. LBP ,A large number of studies have shown that if the number of holes or the total length of holes is insufficient, LBP and VBP There is a large error between the LBP Directly reflects the landslide accumulation layer VBP Therefore, the above traditional engineering methods are used to estimate the soil-rock mixture VBP There are defects in the process, and it is urgent to propose a method based on the limited number and length of drilling holes to accurately evaluate the LBP and VBP The error between the two and the measured value through drilling LBP Accurately estimate soil-rock mixture VBP The method and system provide a scientific basis for improving the design safety and economic benefits of soil-rock mixture related projects and for rationally selecting excavation design and support methods. Summary of the invention

[0003] The purpose of the present invention is to overcome the above-mentioned shortcomings and provide a method and system for estimating the volume stone content based on borehole linear measurement taking into account the three-dimensional shape of the block stone. The present invention constructs a three-dimensional virtual model of a soil-rock mixture that can take into account the real three-dimensional shape of the block stone, innovatively introduces uncertainty factors, and proposes a more accurate method and system for estimating the volume stone content of a soil-rock mixture based on on-site limited borehole linear measurement, which provides a scientific basis for improving the safety and economic benefits of engineering design related to soil-rock mixtures and rationally selecting excavation design and support methods.

[0004] In order to achieve the above technical features, the object of the present invention is achieved as follows: a method for estimating the volumetric stone content based on borehole linear measurement considering the three-dimensional shape of the block stone, characterized in that it includes the following steps: Step 1: Obtain several stone samples of different sizes in the soil-rock mixture study area, prepare a stone sample set and perform CT scanning, and use image digital processing technology to process the CT image of the stone set, separate and reconstruct the three-dimensional morphological model of the stone, analyze the three-dimensional morphological characteristics of the stone, and establish a three-dimensional stone model database containing the three-dimensional morphological characteristic factors of the stone; Step 2: Determine the area size of the soil-rock mixture 3D virtual model based on the maximum stone particle size in the 3D stone model database, and the linear stone content obtained by on-site drilling. LBP , assuming the volumetric stone content VBP 1 is equal to LBP By randomly selecting blocks from the block 3D model database, the volumetric stone content is generated as VBP 1. A three-dimensional virtual model of soil-rock mixture; Step 3: randomly generate parallel edges of the soil-rock mixture model cube and parallel to each other in the soil-rock mixture three-dimensional virtual model construction. n Scan lines, select different numbers of scan lines to calculate the linear stone content LBP , and determine the uncertainty factor of the virtual model UF and sampling length ratio N relationship; Step 4: Sampling length ratio according to on-site drilling N , Linear rock content obtained by drilling measurement LBP , and the uncertainty factors obtained through the virtual model UF and N relationship, determine UF The value of the estimated volumetric stone content on site VBP = LBP + / -( LBP × UF ).

[0005] Preferably, the specific steps in step 1 are: Step 11: Obtain sufficient rock samples of varying sizes in the soil-rock mixture study area, mix the rock with 1-2 mm foam particles in a transparent PVC cylindrical tube to prepare a rock sample, perform CT scanning, and obtain a three-dimensional CT image of the rock sample; Step 12: Convert the CT grayscale image into a three-dimensional CT binary image using a three-dimensional threshold segmentation method; Step 13: Using a three-dimensional median filtering method to reduce noise on the three-dimensional CT binary image; Step 14: Separate each 3D block stone model by using a 3D neighborhood connected component extraction method and a 3D watershed segmentation method; Step 15: mesh the surface of each separated stone block using a triangular mesh to analyze the three-dimensional morphological characteristics of the stone block; Step 16: Establish a three-dimensional model database of the block stone including the three-dimensional morphological characteristic factors of the block stone.

[0006] Preferably, the three-dimensional morphological characteristics of the stone block in step 15 include particle size, three-dimensional shape, sphericity, three-dimensional convexity and three-dimensional roundness.

[0007] Preferably, the specific steps in step 2 are: Step 21: According to the maximum stone particle size in the stone three-dimensional model database, determine that the soil-rock mixture three-dimensional virtual model area is a cubic area, and the area side length is the maximum stone particle size d max 1.5 times of Step 22: Linear rock content obtained from on-site drilling LBP , assuming that the volume of the three-dimensional virtual model of the soil-rock mixture contains VBP 1 is equal to LBP ; Step 23: randomly selecting a block rock model from the block rock three-dimensional model database and placing it within the soil-rock mixture model area; Step 24: Determine whether the block model being placed overlaps with the block model that has been placed, including block vertex overlap, block internal overlap and block edge overlap; the block vertex overlap and block internal overlap are determined by the volume calculated by the vector algorithm, and the edge overlap is determined by the Möller-Trumbore algorithm; Step 25: If it is determined in step 24 that the block being placed does not overlap with the block already placed, the block model is saved in the soil-rock mixture model, otherwise the block is placed again; Step 26: When the total volume of the blocks reaches VBP At 1, the stone placement is finished, and the volume of stone content generated is VBP 1 Three-dimensional virtual model of soil-rock mixture.

[0008] Preferably, the specific steps in step 3 are: Step 31: In the construction of the generated soil-rock mixture three-dimensional virtual model, randomly generate parallel edges of the soil-rock mixture model cube and parallel to each other. n Scan lines, each scan line length L n is equal to the side length of the soil-rock mixture model, where n It needs to be a multiple of 50; Step 32: Use the Möller-Trumbore algorithm to calculate the intersection of each scan line and the block stone model, and count the length of each scan line passing through the block stone model L B , divide this length by the length of the scan line L n , and obtain the linear stone content on the scan line LBP n ; Step 33: From n In the scan lines, randomly select n / 50 scan lines, calculate the average value of the linear stone content measured by each selected scan line LBP a , repeat this process λ times, you can get λ indivual LBP a , and calculate λ indivual LBP a Standard Deviation σ ( LBP a1 , LBP a2 ,…, LBP aλ ), based on the constructed three-dimensional virtual model of soil-rock mixture VBP 1, we can get the randomly selected n / 50 scan lines corresponding to UF = σ ( LBP a1 , LBP a2 ,…, LBP aλ ) / VBP 1; Step 34: Repeat the process of step 33 and randomly select n / 25,3 n / 50,2 n / 25, n / 10, n / 20, n / 5 and n / 2 scan lines, and calculate them separately n / 25,3 n / 50,2 n / 25, n / 10, n / 20, n / 5 and n / 2 scan lines UF value; Step 35: The randomly selected n / 50, n / 25,3 n / 50,2 n / 25, n / 10, n / 20, n / 5 and n / Uncertainty factor corresponding to 2 scan lines UF value, and determine UF Ratio to sampling length N The sampling length ratio N It is the total length of the selected scan line divided by the maximum stone particle size in the three-dimensional virtual model of the soil-rock mixture. d max get.

[0009] Preferably, the specific steps in step 4 are: Step 41: Based on the maximum visible stone particle size in the study area d maxL , and the total length of the on-site drilling L D , determine the sampling length ratio of the on-site drilling N = L D / d maxL ; Step 42: Sampling length ratio according to on-site drilling N , and the uncertainty factor obtained through the virtual model UF and N relationship, determine UF The value of Step 43: According to UF The linear rock content obtained by drilling measurement LBP , estimated volumetric stone content VBP = LBP + / -( LBP × UF ).

[0010] A system for estimating volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the block stone, the system is used to implement the method for estimating volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the block stone, comprising: The block stone three-dimensional model database module obtains a number of block stone samples of different sizes in the soil-rock mixture research area, prepares a block stone sample set and performs CT scanning, and uses image digital processing technology to process the block stone set CT image, separates and reconstructs the block stone three-dimensional morphological model, analyzes the block stone three-dimensional morphological characteristics, and establishes a block stone three-dimensional model database containing block stone three-dimensional morphological characteristic factors; The soil-rock mixture 3D virtual model module is based on the linear rock content obtained by on-site drilling. LBP , and assuming that the volumetric stone content VBP 1 is equal to LBP By randomly selecting blocks from the block 3D model database, the volumetric stone content is generated as VBP 1. A three-dimensional virtual model of soil-rock mixture; UF and N The relationship module randomly generates parallel edges of the cube of the soil-rock mixture model and parallel to each other in the three-dimensional virtual model construction of the soil-rock mixture. n Scan lines, select different numbers of scan lines to calculate the linear stone content LBP , and determine the uncertainty factor of the virtual model UF and sampling length ratio N relationship; The volume rock content estimation module is based on the sampling length ratio of the on-site drilling N , Linear rock content obtained by drilling measurement LBP , and the uncertainty factors obtained through the virtual model UF and N relationship, determine UF The value of the estimated volumetric stone content on site VBP = LBP + / -( LBP × UF ).

[0011] The soil-rock mixture three-dimensional virtual model module includes: The model area determination unit determines that the side length of the cube area of ​​the soil-rock mixture three-dimensional virtual model is 1.5 times the maximum block stone particle size according to the maximum block stone particle size in the block stone three-dimensional model database; The model volume rock content determination unit is based on the linear rock content obtained by on-site drilling LBP , assuming that the volume of the three-dimensional virtual model of the soil-rock mixture contains VBP 1 is equal to LBP ; The block stone random placement unit randomly selects block stone models from the block stone three-dimensional model database and places them within the soil-rock mixture model area; A block overlap determination unit determines whether the block model being placed has block vertex overlap, block internal overlap, and block edge overlap with the block model already placed; The block stone model saving unit saves the block stone model in the soil-rock mixture model if the block stone overlap determination unit determines that the block stone being placed does not overlap with the already placed blocks stone; The three-dimensional virtual model of the soil-rock mixture completes the unit. When the total volume of the blocks reaches VBP 1, the stone placement is finished, that is, the volume of stone content generated is VBP 1 Three-dimensional virtual model of soil-rock mixture.

[0012] The present invention has the following beneficial effects: 1. The present invention introduces uncertainty factors through innovation UF , a linear rock content method based on field drilling measurement with a limited number and length of boreholes is proposed. LBP Estimated volumetric rock content of soil-rock mixture VBP method and system.

[0013] 2. The present invention combines CT scanning technology with image digital processing technology to construct a three-dimensional virtual model of soil-rock mixture that can take into account the real three-dimensional shape of the block stone. Based on the innovative virtual drilling technology, the UF The relationship between the total length of the borehole and the maximum particle size of the rock mass can more accurately estimate the volumetric rock content of the soil-rock mixture considering the three-dimensional shape of the rock mass. VBP , providing a scientific basis for improving the design safety and economic benefits of soil-rock mixture related projects and for rationally selecting excavation design and support methods. BRIEF DESCRIPTION OF THE DRAWINGS

[0014] The present invention will be further described below in conjunction with the accompanying drawings and embodiments.

[0015] Figure 1 The present invention is a flow chart of a method for estimating volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the block of stone.

[0016] Figure 2 Schematic diagram of a grayscale histogram in an embodiment of the present invention.

[0017] Figure 3 It is a schematic diagram of a three-dimensional soil-rock mixture virtual model in an embodiment of the present invention.

[0018] Figure 4 It is a schematic diagram of a scanning line passing through a three-dimensional soil-rock mixture virtual model in an embodiment of the present invention.

[0019] Figure 5 Schematic diagram of the relationship between UF and N in an embodiment of the present invention.

[0020] Figure 6 The present invention is a schematic diagram of the structure of a system for estimating volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the rock block. DETAILED DESCRIPTION

[0021] The principles and features of the present invention are described below in conjunction with the accompanying drawings. The examples given are only used to explain the present invention and are not used to limit the scope of the present invention.

[0022] like Figure 1 As shown, a method for estimating volumetric stone content based on borehole linear measurement considering the three-dimensional shape of the block stone includes the following steps: Step 1: Obtain several stone samples of different sizes in the soil-rock mixture study area, prepare a stone sample set and perform CT scanning, and use image digital processing technology to process the CT image of the stone set, separate and reconstruct the three-dimensional morphological model of the stone, analyze the three-dimensional morphological characteristics of the stone, and establish a three-dimensional stone model database containing the three-dimensional morphological characteristic factors of the stone.

[0023] Step 2: Determine the area size of the soil-rock mixture 3D virtual model based on the maximum stone particle size in the 3D stone model database, and the linear stone content obtained by on-site drilling. LBP , assuming the volumetric stone content VBP 1 is equal to LBP By randomly selecting blocks from the block 3D model database, the volumetric stone content is generated as VBP 1 Three-dimensional virtual model of soil-rock mixture.

[0024] Step 3: randomly generate parallel edges of the soil-rock mixture model cube and parallel to each other in the soil-rock mixture three-dimensional virtual model construction. n Scan lines, select different numbers of scan lines to calculate the linear stone content LBP , and determine the uncertainty factor of the virtual model UF and sampling length ratio N relationship.

[0025] Step 4: Sampling length ratio according to on-site drilling N , Linear rock content obtained by drilling measurement LBP , and the uncertainty factors obtained through the virtual model UF and N relationship, determine UF The value of the estimated volumetric stone content on site VBP = LBP + / -( LBP× UF ).

[0026] In this embodiment, in step 1, the specific implementation of establishing a block stone model database through CT scanning is as follows: Step 11: Obtain sufficient rock samples of varying sizes in the soil-rock mixture study area, mix the rock with 1-2 mm foam particles in a transparent PVC cylindrical tube to prepare a rock sample, perform CT scanning, and obtain a three-dimensional CT image of the rock sample.

[0027] Several stone samples of different sizes in the soil-rock mixture study area were obtained and placed in a transparent PVC cylindrical tube with a diameter of 300 mm and a height of 300 mm. CT samples of the stone collection were prepared and scanned using the Philips Ingenuity CT system. The stone particle size in the prepared CT sample was between 10 and 100 mm. Lightweight foam particles (with a particle size of 1-2 mm) were used instead of soil to fill the gaps between the stones. From the perspective of sample preparation, the stones were separated so that they were kept as far away from each other as possible in the sample, reducing the difficulty of separating the stones in the later image processing. A transparent PVC cylindrical tube was used to hold the stones and lightweight foam particles, while allowing X-rays to pass through the tube wall, so that clear three-dimensional CT images of the stone collection could be obtained.

[0028] Step 12: Use the three-dimensional threshold segmentation method to convert the CT grayscale image into a three-dimensional CT binary image.

[0029] The CT grayscale image is converted into a three-dimensional CT binary image by a three-dimensional threshold segmentation method, so as to further extract the block stone information. First, by analyzing the grayscale histogram of the sample CT grayscale image sequence, the valley value between two grayscale value peaks is selected as the threshold, and the CT grayscale image is converted into a three-dimensional binary image (black and white image), that is, the grayscale level of pixels with a grayscale less than the threshold is changed to 0 (black), and the grayscale level of pixels with a grayscale greater than the threshold is changed to 255 (white); in this embodiment, the threshold value is selected. T is 61, such as Figure 2 shown.

[0030] Step 13: Use the 3D median filtering method to reduce noise on the 3D CT binary image.

[0031] In order to eliminate the noise in the binarized CT image, a three-dimensional median filter with a range of 3×3×3 is used to reduce the noise of the three-dimensional CT binary image.

[0032] Step 14: Use the three-dimensional neighborhood connected component extraction method and the three-dimensional watershed segmentation method to separate each three-dimensional stone model.

[0033] In the processing of CT images, the 26-neighborhood connected component extraction method is used to separate the individual rocks. However, some rocks that are closely in contact may not be successfully separated by the connected component method. In this case, the three-dimensional watershed segmentation method is further used to separate the unseparated rocks for the second time to ensure that each rock is completely separated.

[0034] Step 15: Use triangular mesh to mesh the surface of each separated stone block and analyze the three-dimensional morphological characteristics of the stone block (particle size, three-dimensional shape, sphericity, three-dimensional convexity, three-dimensional roundness, etc.).

[0035] For each block of stone that has been segmented, the Generate Surface function in AVIZO is used to generate a three-dimensional mesh on the surface of the block, and the geometric model of the block is saved as a ".stl" file type. The three-dimensional morphological characteristics of the block are analyzed through MATLAB programming, including but not limited to particle size, three-dimensional shape, sphericity, three-dimensional convexity, three-dimensional roundness and other morphological factors.

[0036] Step 16: Establish a three-dimensional model database of the block stone including the three-dimensional morphological characteristic factors of the block stone.

[0037] The geometric model of the reconstructed and separated stone blocks and the stone block morphological factors are all imported into MATLAB to establish a three-dimensional stone block model database containing the three-dimensional morphological characteristic factors of the stone blocks.

[0038] In this example, in step 2, constructing a three-dimensional virtual model of a soil-rock mixture is specifically implemented as follows: Step 21: According to the maximum stone particle size in the stone three-dimensional model database, determine that the soil-rock mixture three-dimensional virtual model area is a cubic area, and the side length of the area is 1.5 times the maximum stone particle size.

[0039] Since the maximum stone particle size in the stone 3D model database is 100 mm, the area size of the soil-rock mixture 3D virtual model is set to 150×150×150 mm.

[0040] Step 22: Linear rock content obtained from on-site drilling LBP , assuming that the volume of the three-dimensional virtual model of the soil-rock mixture contains VBP 1 is equal to LBP .

[0041] Linear rock content obtained by on-site drilling LBP is 30%, so the volume of the three-dimensional virtual model of the soil-rock mixture contains VBP 1 is 30%.

[0042] Step 23: Randomly select a block rock model from the block rock three-dimensional model database and place it within the soil-rock mixture model area.

[0043] The block stone models in the block stone three-dimensional model database are randomly selected and placed in the area of ​​150×150×150mm.

[0044] Step 24: Determine whether the block model being placed overlaps with the block model that has been placed (including block vertex overlap, block internal overlap and block edge overlap); the block vertex overlap and block internal overlap are determined by the volume size calculated by the vector algorithm, and the edge overlap is determined by the Möller-Trumbore algorithm.

[0045] The block stone model being generated in the area needs to be judged whether it overlaps with the block stone model that has been placed, and whether there is vertex overlap, internal overlap and edge overlap between the block stone model being placed and the block stone model that has been placed.

[0046] Vertex overlap and internal overlap are judged using a vector algorithm, by judging whether the vertex of the stone model being placed is inside the stone model that has been placed. If the volume formed by the vertex and the triangular faces of the stone model that has been placed is equal to the original volume of the stone model that has been placed, then it is judged that vertex overlap or internal overlap has occurred; otherwise, it is judged that no vertex overlap or internal overlap has occurred. Similarly, the same method is used to determine whether the stone that has been placed is inside the stone model that is being placed. If so, it is still judged that the two stones have vertex overlap or internal overlap.

[0047] Edge overlap is determined using the Möller-Trumbore algorithm. Through the two vertices of each edge of the stone model being placed, one of the vertices is used as the origin, and a ray is drawn through the direction of the line connecting the two points. If the ray intersects with the stone model that has been placed, it is determined whether the intersection is on the edge of the two vertices. If so, it is determined that the two stones have edge overlap, otherwise it is determined that there is no edge overlap.

[0048] Step 25: If it is determined in step 24 that the block being placed does not overlap with the block already placed, the block model is saved in the soil-rock mixture model, otherwise the block is placed again.

[0049] If it is determined in step 24 that the block being placed does not overlap with the blocks already placed, the block model is saved in the soil-rock mixture model, otherwise, the process returns to step 23 and random block placement is performed again.

[0050] Step 26: When the total volume of the blocks reaches VBP 1, the stone placement is finished, that is, the volume of stone content generated is VBP 1 Three-dimensional virtual model of soil-rock mixture.

[0051] When the total volume of the blocks reaches VBP When 1=30%, the placement of blocks of stone is terminated, and the virtual model of the soil-rock mixture with a volumetric stone content of 30% is saved. Figure 3 shown.

[0052] In this example, step 3 is specifically implemented as follows: Step 31: In the construction of the generated soil-rock mixture three-dimensional virtual model, randomly generate parallel edges of the soil-rock mixture model cube and parallel to each other. n Scan lines, each scan line length L n is equal to the side length of the soil-rock mixture model, where n Needs to be a multiple of 50.

[0053] In the construction of the generated three-dimensional virtual model of the soil-rock mixture, 100 scanning lines with a length of 150 mm and parallel to each other are randomly generated, parallel to the height side length of the soil-rock mixture model, such as Figure 4 shown.

[0054] Step 32: Use the Möller-Trumbore algorithm to calculate the intersection of each scan line and the block stone model, and count the length of each scan line passing through the block stone model L B , divide this length by the length of the scan line L n , and obtain the linear stone content on the scan line LBP n .

[0055] The Möller-Trumbore algorithm is used to calculate the intersection of each scan line and each block model, and the coordinates of the second intersection point are subtracted from the coordinates of the first intersection point. z The absolute value of the value, the coordinate value of the third intersection minus the coordinate value of the fourth intersection z The absolute value of the value, and so on, add up all these absolute values ​​and divide them by the length of the scanning line to get the linear stone content measured by the scanning line LBP 1. Similarly, 100 linear stone content values ​​are calculated ( LBP 1, LBP 2, …, LBP 100 ).

[0056] Step 33: From n In the scan lines, randomly select n / 50 scan lines, calculate the average value of the linear stone content measured by each selected scan line LBP a , repeat this process λtimes, you can get λ indivual LBP a , and calculate λ indivual LBP a Standard Deviation σ ( LBP a1 , LBP a2 ,…, LBP aλ ), based on the constructed three-dimensional virtual model of soil-rock mixture VBP 1, we can get the randomly selected n / 50 scan lines corresponding to UF = σ ( LBP a1 , LBP a2 ,…, LBP aλ ) / VBP 1.

[0057] Randomly select 2 scan lines from 100 scan lines and calculate the average linear stone content of each selected scan line LBP a ; Repeat this process 40 times to get 40 LBP a , calculate 40 average values LBP a Standard Deviation σ ( LBP a1 , LBP a2 ,…, LBP a40 ), and divide it by the three-dimensional virtual model of the soil-rock mixture VBP 1 value, we can get UF value.

[0058] Step 34: Repeat the process of step 33 and randomly select n / 25,3 n / 50,2 n / 25, n / 10, n / 20, n / 5 and n / 2 scan lines, and calculate them separately n / 25,3 n / 50,2 n / 25, n / 10, n / 20, n / 5 and n / 2 scan lines UF value.

[0059] Repeat the process of step 33, randomly select 4, 6, 8, 10, 15, 20 and 50 scan lines from 100 scan lines, and calculate the 4, 6, 8, 10, 15, 20 and 50 scan lines respectively. UF value.

[0060] Step 35: The randomly selected n / 50, n / 25,3 n / 50,2 n / 25, n / 10, n / 20, n / 5 and n / Uncertainty factor corresponding to 2 scan lines UF value, and determine UF Ratio to sampling length N The sampling length ratio N It is the total length of the selected scan line divided by the maximum stone particle size in the three-dimensional virtual model of the soil-rock mixture. d max Get, such as Figure 5 shown.

[0061] The uncertainty factors corresponding to the randomly selected 2, 4, 6, 8, 10, 15, 20 and 50 scan lines obtained in step 33 and step 34 are UF value and determined by data fitting UF Ratio to sampling length N In this example, the relationship is as follows: UF = -0.088 ln(N) + 0.381; Taking two randomly selected scanning lines as an example, the scanning line length is 150mm, the maximum stone particle size is 100mm, and the length ratio N =150*2 / 100=3.

[0062] In this example, the specific implementation of step 4 is as follows: Step 41: Based on the maximum visible stone particle size in the study area d maxL , and the total length of the on-site drilling L D , determine the sampling length ratio of the on-site drilling N = L D / d maxL .

[0063] Total length of drilling on siteL D The largest visible particle size of exposed rock blocks observed in the landslide study area is 40m. d maxL It is about 2m, and the sampling length ratio is calculated. N = L D / d maxL =40 / 2=20.

[0064] Step 42: Sampling length ratio according to on-site drilling N , and the uncertainty factor obtained through the virtual model UF and N relationship, determine UF The value of .

[0065] According to the sampling length ratio of the on-site drilling N= 20, Substitute into the relation UF = -0.088 ln(N) + 0.381, calculated UF =-0.088 ln (20) + 0.381=0.117.

[0066] Step 43: According to UF The linear rock content obtained by drilling measurement LBP , estimated volumetric stone content VBP = LBP + / -( LBP UF ).

[0067] The estimated volume of stone content on site is VBP =30%+ / -(30% 0.117)=26.5%~33.5%.

[0068] like Figure 6 As shown, a system for estimating volumetric stone content based on borehole linear measurement and considering the three-dimensional morphology of block stones comprises: a block stone three-dimensional model database module, which obtains a number of block stone samples of different sizes in the soil-rock mixture research area, prepares a block stone sample set and performs CT scanning, and uses image digital processing technology to process the block stone set CT image, separates and reconstructs the block stone three-dimensional morphological model, analyzes the three-dimensional morphological characteristics of the block stone, and establishes a block stone three-dimensional model database containing the block stone three-dimensional morphological characteristic factors; a soil-rock mixture three-dimensional virtual model module, which obtains the linear stone content obtained by on-site drilling according to the linear stone content obtained by on-site drilling. LBP , and assuming that the volumetric stone content VBP 1 is equal to LBPBy randomly selecting blocks from the block 3D model database, the volumetric stone content is generated VBP 1. A three-dimensional virtual model of soil-rock mixture; UF and N The relationship module randomly generates parallel edges of the cube of the soil-rock mixture model and parallel to each other in the three-dimensional virtual model construction of the soil-rock mixture. n Scan lines, select different numbers of scan lines to calculate the linear stone content LBP , and determine the uncertainty factor of the virtual model UF and sampling length ratio N The volume rock content estimation module uses the sampling length ratio of the on-site drilling holes to estimate the volume rock content. N , Linear rock content obtained by drilling measurement LBP , and the uncertainty factors obtained through the virtual model UF and N relationship, determine UF The value of the estimated volumetric stone content on site VBP = LBP + / -( LBP × UF ).

[0069] In this embodiment, the soil-rock mixture three-dimensional virtual model module includes: a model area determination unit, which determines that the side length of the cube area of ​​the soil-rock mixture three-dimensional virtual model is 1.5 times the maximum block stone size according to the maximum block stone size in the block stone three-dimensional model database; a model volume stone content determination unit, which determines the linear stone content obtained by on-site drilling. LBP , assuming that the volume of the three-dimensional virtual model of the soil-rock mixture contains VBP 1 is equal to LBP ; The block random placement unit randomly selects a block model from the block three-dimensional model database and places it within the soil-rock mixture model area; the block overlap judgment unit determines whether the block model being placed overlaps with the block model that has been placed in terms of block vertex overlap, block internal overlap, and block edge overlap; the block model preservation unit saves the block model in the soil-rock mixture model if the block overlap judgment unit determines that the block being placed does not overlap with the block that has been placed; the soil-rock mixture three-dimensional virtual model completion unit, when the total volume of the placed block reaches VBP 1, the stone placement is finished, that is, the volume of stone content generated is VBP 1 Three-dimensional virtual model of soil-rock mixture.

[0070] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for estimating volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the block, characterized in that: The following steps are involved: Step 1: Obtain several stone samples of different sizes in the soil-rock mixture study area, prepare a stone sample set and perform CT scanning, and use image digital processing technology to process the CT image of the stone set, separate and reconstruct the three-dimensional morphological model of the stone, analyze the three-dimensional morphological characteristics of the stone, and establish a three-dimensional stone model database containing the three-dimensional morphological characteristic factors of the stone; Step 2: Determine the area size of the soil-rock mixture 3D virtual model based on the maximum stone particle size in the 3D stone model database, and calculate the linear stone content obtained by on-site drilling. LBP , assuming the volumetric stone content VBP 1 is equal to LBP By randomly selecting blocks from the block 3D model database, the volumetric stone content is generated as VBP 1. A three-dimensional virtual model of soil-rock mixture; Step 3: randomly generate parallel edges of the soil-rock mixture model cube and parallel to each other in the soil-rock mixture three-dimensional virtual model construction. n Scan lines, select different numbers of scan lines to calculate the linear stone content LBP , and determine the uncertainty factor of the virtual model UF and sampling length ratio N relationship; Step 4: Sampling length ratio based on on-site drilling N , Linear rock content obtained by drilling measurement LBP , and the uncertainty factors obtained through the virtual model UF and N relationship, determine UF The value of the estimated volumetric stone content on site VBP = LBP + / -( LBP × UF ).

2. According to claim 1, a method for estimating volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the block, characterized in that: The specific steps in step 1 are: Step 11: Obtain sufficient rock samples of varying sizes in the soil-rock mixture study area, mix the rock with 1-2 mm foam particles in a transparent PVC cylindrical tube to prepare a rock sample, perform CT scanning, and obtain a three-dimensional CT image of the rock sample; Step 12: Convert the CT grayscale image into a three-dimensional CT binary image using a three-dimensional threshold segmentation method; Step 13: Using a 3D median filter method to reduce noise on the 3D CT binary image; Step 14: Separate each 3D block stone model by using a 3D neighborhood connected component extraction method and a 3D watershed segmentation method; Step 15: mesh the surface of each separated stone block using a triangular mesh to analyze the three-dimensional morphological characteristics of the stone block; Step 16: Establish a three-dimensional model database of the block stone including the three-dimensional morphological characteristic factors of the block stone.

3. According to claim 2, a method for estimating volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the block, characterized in that: The three-dimensional morphological characteristics of the stone blocks in step 15 include particle size, three-dimensional shape, sphericity, three-dimensional convexity and three-dimensional roundness.

4. The method for estimating volumetric stone content based on borehole linear measurement considering the three-dimensional shape of the block of stone according to claim 3, characterized in that: The specific steps in step 2 are: Step 21: According to the maximum stone particle size in the stone three-dimensional model database, determine that the soil-rock mixture three-dimensional virtual model area is a cubic area, and the area side length is the maximum stone particle size d max 1.5 times of Step 22: Linear rock content obtained from on-site drilling LBP , assuming that the volume of the three-dimensional virtual model of the soil-rock mixture contains VBP 1 is equal to LBP ; Step 23: randomly selecting a block rock model from the block rock three-dimensional model database and placing it within the soil-rock mixture model area; Step 24: Determine whether the block model being placed overlaps with the block model that has been placed, including block vertex overlap, block internal overlap and block edge overlap; the block vertex overlap and block internal overlap are determined by the volume calculated by the vector algorithm, and the edge overlap is determined by the Möller-Trumbore algorithm; Step 25: If it is determined in step 24 that the block being placed does not overlap with the blocks that have been placed, then the block model is saved in the soil-rock mixture model, otherwise the block is placed again; Step 26: When the total volume of the blocks reaches VBP 1, the stone placement is finished, that is, the volume of stone content generated is VBP 1 Three-dimensional virtual model of soil-rock mixture.

5. A method for estimating volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the block of stone according to claim 4, characterized in that: The specific steps in step 3 are: Step 31: In the construction of the generated soil-rock mixture three-dimensional virtual model, randomly generate parallel edges of the soil-rock mixture model cube and parallel to each other. n Scan lines, each scan line length L n is equal to the side length of the soil-rock mixture model, where n It needs to be a multiple of 50; Step 32: Use the Möller-Trumbore algorithm to calculate the intersection of each scan line and the block stone model, and count the length of each scan line passing through the block stone model L B , divide this length by the length of the scan line L n , and obtain the linear stone content on the scan line LBP n ; Step 33: From n In the scan lines, randomly select n / 50 scan lines, calculate the average value of the linear stone content measured by each selected scan line LBP a , repeat this process λ times, you can get λ indivual LBP a , and calculate λ indivual LBP a Standard Deviation σ ( LBP a1 , LBP a2 ,…, LBP aλ ), based on the constructed three-dimensional virtual model of soil-rock mixture VBP 1, we can get the randomly selected n / 50 scan lines corresponding to UF = σ ( LBP a1 , LBP a2 ,…, LBP aλ ) / VBP 1; Step 34: Repeat the process of step 33 and randomly select n / 25,3 n / 50,2 n / 25, n / 10, n / 20, n / 5 and n / 2 scan lines, and calculate them separately n / 25,3 n / 50,2 n / 25, n / 10, n / 20, n / 5 and n / 2 scan lines UF value; Step 35: The randomly selected n / 50, n / 25,3 n / 50,2 n / 25, n / 10, n / 20, n / 5 and n / Uncertainty factor corresponding to 2 scan lines UF value, and determine UF Ratio to sampling length N The sampling length ratio N It is the total length of the selected scan line divided by the maximum stone particle size in the three-dimensional virtual model of the soil-rock mixture. d max get.

6. A method for estimating volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the block of stone according to claim 5, characterized in that: The specific steps in step 4 are: Step 41: Based on the maximum visible stone particle size in the study area d maxL , and the total length of the on-site drilling L D , determine the sampling length ratio of the on-site drilling N = L D / d maxL ; step 42: According to the sampling length ratio of the on-site drilling N , and the uncertainty factor obtained through the virtual model UF and N relationship, determine UF The value of Step 43: According to UF The linear rock content obtained by drilling measurement LBP , estimated volumetric stone content VBP = LBP + / -( LBP × UF ).

7. A system for estimating the volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the block of stone, the system being used to implement the method for estimating the volumetric stone content based on borehole linear measurement taking into account the three-dimensional shape of the block of stone as claimed in any one of claims 1 to 6, characterized in that: include: The block stone three-dimensional model database module obtains a number of block stone samples of different sizes in the soil-rock mixture research area, prepares a block stone sample set and performs CT scanning, and uses image digital processing technology to process the block stone set CT image, separates and reconstructs the block stone three-dimensional morphological model, analyzes the block stone three-dimensional morphological characteristics, and establishes a block stone three-dimensional model database containing block stone three-dimensional morphological characteristic factors; The soil-rock mixture 3D virtual model module is based on the linear rock content obtained by on-site drilling. LBP , and assuming that the volumetric stone content VBP 1 is equal to LBP By randomly selecting blocks from the block 3D model database, the volumetric stone content is generated as VBP 1. A three-dimensional virtual model of soil-rock mixture; UF and N The relationship module randomly generates parallel edges of the cube of the soil-rock mixture model and parallel to each other in the three-dimensional virtual model construction of the soil-rock mixture. n Scan lines, select different numbers of scan lines to calculate the linear stone content LBP , and determine the uncertainty factor of the virtual model UF and sampling length ratio N relationship; The volume rock content estimation module is based on the sampling length ratio of the on-site drilling N , Linear rock content obtained by drilling measurement LBP , and the uncertainty factors obtained through the virtual model UF and N relationship, determine UF The value of the estimated volumetric stone content on site VBP = LBP + / -( LBP × UF ).

8. A system for estimating volumetric stone content based on borehole linear measurement considering the three-dimensional shape of a rock mass according to claim 7, characterized in that: The soil-rock mixture three-dimensional virtual model module includes: The model area determination unit determines that the side length of the cube area of ​​the soil-rock mixture three-dimensional virtual model is 1.5 times the maximum block stone particle size according to the maximum block stone particle size in the block stone three-dimensional model database; The model volume rock content determination unit is based on the linear rock content obtained by on-site drilling LBP , assuming that the volume of the three-dimensional virtual model of the soil-rock mixture contains VBP 1 is equal to LBP ; The block stone random placement unit randomly selects block stone models from the block stone three-dimensional model database and places them within the soil-rock mixture model area; A block overlap determination unit determines whether the block model being placed has block vertex overlap, block internal overlap, and block edge overlap with the block model already placed; The block stone model saving unit saves the block stone model in the soil-rock mixture model if the block stone overlap determination unit determines that the block stone being placed does not overlap with the already placed blocks stone; The three-dimensional virtual model of the soil-rock mixture completes the unit. When the total volume of the blocks reaches VBP 1, the stone placement is finished, that is, the volume of stone content generated is VBP 1 Three-dimensional virtual model of soil-rock mixture.

Images