CN114187423A - Surrounding rock fracture reconstruction method, electronic equipment and storage medium for three-dimensional simulation test - Google Patents
Surrounding rock fracture reconstruction method, electronic equipment and storage medium for three-dimensional simulation test Download PDFInfo
- Publication number
- CN114187423A CN114187423A CN202111509800.2A CN202111509800A CN114187423A CN 114187423 A CN114187423 A CN 114187423A CN 202111509800 A CN202111509800 A CN 202111509800A CN 114187423 A CN114187423 A CN 114187423A
- Authority
- CN
- China
- Prior art keywords
- coal mining
- fracture
- grid
- model
- solid phase
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 73
- 238000004088 simulation Methods 0.000 title claims abstract description 67
- 239000011435 rock Substances 0.000 title claims abstract description 44
- 238000012360 testing method Methods 0.000 title claims abstract description 39
- 238000003860 storage Methods 0.000 title claims abstract description 13
- 239000003245 coal Substances 0.000 claims abstract description 234
- 238000005065 mining Methods 0.000 claims abstract description 231
- 239000007790 solid phase Substances 0.000 claims abstract description 104
- 238000009826 distribution Methods 0.000 claims abstract description 44
- 238000010146 3D printing Methods 0.000 claims abstract description 13
- 239000002245 particle Substances 0.000 claims description 28
- 238000005259 measurement Methods 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 16
- 230000002093 peripheral effect Effects 0.000 claims description 10
- 238000012544 monitoring process Methods 0.000 abstract description 13
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000002591 computed tomography Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 3
- 241001502050 Acis Species 0.000 description 2
- 238000000627 alternating current impedance spectroscopy Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000003708 edge detection Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 239000012071 phase Substances 0.000 description 1
- 238000007639 printing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a surrounding rock fracture reconstruction method, electronic equipment and a storage medium for a three-dimensional simulation test, wherein the method comprises the following steps: acquiring porosity of a simulated coal mining model in a three-dimensional simulation experiment as the porosity after coal mining; acquiring the distribution probability of solid phase growing nuclei statistically obtained by each grid on the coal mining model and the growing probability of the solid phase growing nuclei of each grid in multiple directions; generating a fracture digital model after coal mining according to the porosity after coal mining, the distribution probability of solid phase growth nuclei of each grid and the growth probability of the solid phase growth nuclei of each grid in multiple directions; and outputting the coal mining fracture digital model to a three-dimensional printing device, and reconstructing the three-dimensional printing device into a fracture three-dimensional model according to the coal mining fracture digital model. The method is combined with the model parameters of the three-dimensional simulation model to establish the fracture three-dimensional model, so that the reconstruction of the model is realized, and the problems of low monitoring precision of the surrounding rock fracture field and inaccuracy of the three-dimensional reconstruction model in the three-dimensional simulation experiment are solved.
Description
Technical Field
The invention relates to the technical field of coal mine correlation, in particular to a surrounding rock fracture reconstruction method, electronic equipment and a storage medium for a three-dimensional simulation test.
Background
Currently, there are two types of devices for monitoring the fracture field of the surrounding rock, drilling television and CT imaging.
Borehole television imagers, such as KDVJ-400, use high resolution color television cameras for monitoring and display on a liquid crystal display screen the borehole wall structure, monitoring the overall conditions within the borehole, and recording. The digital panoramic drilling shooting system adopts a new concept of panoramic images, utilizes a digital technology for processing, has good real-time performance, large observation range, high measurement precision, automatic detection process and strong digital processing capability, can accurately obtain parameters such as structural surface occurrence, fracture width and the like in a drill hole, and can form a plane development picture and a three-dimensional drill hole core picture which are processed by seamless splicing.
The CT scanning imaging system is used as a nondestructive testing means, and can rapidly and nondestructively acquire the distribution characteristics of the internal structure of the rock by using the scanning imaging principle and the visualization function of three-dimensional reconstruction. In an industrial CT scanning imaging experiment, when X rays penetrate through a rock test piece, the internal structure of the rock test piece is composed of mineral components with different densities, the absorption coefficients of all points to the X rays are different, the intensity of the X rays is attenuated in the scanning process, the rays contain the internal structure density information of the scanned rock test piece, and an imaging system receives and converts the internal structure density information into a digital image form, so that the internal structure distribution condition of the rock test piece is conveniently observed. And interfacing with an industrial CT scanning device through VG Studio max software, processing the original data and deriving a CT scanning slice image. Vectorizing a CT slice image stack by using three-dimensional visual image processing software Avizo, performing image threshold segmentation by using an interactive threshold module to extract a crack part in an image, and displaying the appearance and distribution characteristics of the crack by using a Volume Rendering module.
However, although the existing monitoring device can monitor the fracture field inside the model, the accuracy is low because the interference of contact monitoring is not considered, and the result deviation is large. The existing model reconstruction method can complete reconstruction work, but the algorithm is complicated, and the model result error is caused by insufficient perfection on a connection contact point when a multi-phase problem is encountered.
Disclosure of Invention
Based on this, it is necessary to provide a surrounding rock fracture reconstruction method, an electronic device and a storage medium for use in a three-dimensional simulation test, aiming at the technical problems of low monitoring precision of a surrounding rock fracture field and inaccuracy of a three-dimensional reconstruction model in a three-dimensional simulation test in the prior art.
The invention provides a surrounding rock fracture reconstruction method used in a three-dimensional simulation test, which comprises the following steps:
acquiring porosity of a simulated coal mining model in a three-dimensional simulation experiment as the porosity after coal mining;
acquiring the distribution probability of solid phase growing nuclei statistically obtained by each grid on the coal mining model and the growing probability of the solid phase growing nuclei of each grid in multiple directions;
generating a fracture digital model after coal mining according to the porosity after coal mining, the distribution probability of solid phase growth nuclei of each grid and the growth probability of the solid phase growth nuclei of each grid in multiple directions;
and outputting the coal mining crack digital model to a three-dimensional printing device, and reconstructing the coal mining crack digital model into a crack three-dimensional model by the three-dimensional printing device.
Further:
in the three-dimensional simulation experiment, the porosity of the simulated coal mining model after coal mining is used as the porosity after coal mining, and the method specifically comprises the following steps: acquiring a coal mining model image of a simulated coal mining model in a three-dimensional simulation experiment, and calculating the porosity of the simulated coal mining model according to the coal mining model image to serve as the porosity after coal mining;
the method for acquiring the distribution probability of the solid phase growth nuclei obtained by each grid statistic on the coal mining model and the growth probability of the solid phase growth nuclei of each grid in multiple directions specifically comprises the following steps: dividing the coal mining model image into a plurality of grids, and according to the image information of the coal mining model image, counting the distribution probability of solid phase growth nuclei of each grid and the growth probability of the solid phase growth nuclei of each grid in a plurality of directions.
Furthermore, in the three-dimensional simulation experiment, the acquiring of the post-coal mining model image of the simulated post-coal mining model specifically includes:
in the three-dimensional simulation experiment, images of a coal mining model after simulated coal mining, which are shot by a plurality of camera devices in the coal mining model, are obtained, and the images shot by the camera devices are fused to obtain a model image after coal mining.
Still further, still include:
identifying a fracture from the model image after coal mining as a real-time measurement fracture after coal mining, and calculating a fracture change value of the real-time measurement fracture after coal mining as a real-time measurement fracture change value after coal mining;
identifying a fracture from the coal mining fracture digital model as a fracture to be verified, and calculating a fracture change value to be verified of the fracture to be verified;
and comparing the actual measured fracture change value after coal mining with the fracture change value to be verified, if the difference value between the actual measured fracture change value after coal mining and the fracture change value to be verified is smaller than a preset difference threshold value, judging that the fracture digital model after coal mining is effective, otherwise, judging that the fracture digital model after coal mining is ineffective, and generating a new fracture digital model after coal mining according to the porosity after coal mining, the distribution probability of solid phase growth nuclei of each grid and the growth probability of the solid phase growth nuclei of each grid in multiple directions.
Still further, still include:
acquiring the porosity of a coal mining model before coal mining as the porosity before coal mining;
dividing the coal mining model into a plurality of grids, and acquiring the initial distribution probability of solid phase growth nuclei of each grid and the initial growth probability of the solid phase growth nuclei of each grid in a plurality of directions;
generating an initial fracture digital model according to the porosity before coal mining, the initial distribution probability of the solid phase growth nuclei of each grid and the initial growth probability of the solid phase growth nuclei of each grid in multiple directions;
identifying a fracture from the initial fracture digital model as an initial fracture;
the calculating of the change value of the crack to be verified specifically comprises the following steps: and comparing the initial crack with the crack to be verified, and determining the change value of the crack to be verified.
Still further, still include:
acquiring a pre-coal mining model image of a simulated pre-coal mining model in a three-dimensional simulation experiment;
identifying a fracture from the pre-coal mining model image as a pre-coal mining actual measurement fracture;
the method for calculating the fracture change value of the actual-measured fracture after coal mining as the actual-measured fracture change value after coal mining specifically comprises the following steps:
and comparing the coal mining actual measurement crack with the coal mining actual measurement crack to obtain a crack change value of the coal mining actual measurement crack as a coal mining actual measurement crack change value.
Further, still include:
acquiring a coal mining process model image of a coal mining model in a simulated coal mining process in a three-dimensional simulation experiment;
and identifying a crack from the coal mining process model image, and generating and displaying change data of the crack about a time axis.
Further, the generating a coal mining fracture digital model according to the porosity after coal mining, the distribution probability of the solid phase growth nuclei of each grid, and the growth probability of the solid phase growth nuclei of each grid in multiple directions specifically includes:
constructing a digital model based on the coordinates of the grid;
traversing all grids of the digital model, acquiring a grid node random number and the distribution probability of a solid phase growth core of each grid, if the grid node random number is smaller than the distribution probability of the solid phase growth core of each grid, generating a solid phase particle core by each grid at a corresponding coordinate position in the digital model, and generating a grid of the solid phase particle cores as an effective grid;
performing directional growth calculation on each effective grid until the porosity of the digital model reaches the porosity after coal mining;
forming contour equipotential lines based on an effective grid of a digital model, and generating a fracture digital model after coal mining based on the contour equipotential lines;
the calculation of the directional growth specifically comprises:
for an effective grid, acquiring a random number of peripheral nodes and the growth probability of a solid phase growth core of the effective grid in multiple directions, taking the direction in which the growth probability of the solid phase growth core of the effective grid is greater than the random number of the peripheral nodes as a growth direction, taking an adjacent grid positioned in the growth direction of the effective grid as a growth grid, generating a solid phase particle kernel at a corresponding coordinate position of the growth grid in the digital model, updating the growth grid generating the solid phase particle kernel into the effective grid, and calculating the porosity of the digital model, if the porosity of the digital model reaches the porosity after coal mining, ending the directional growth calculation, otherwise, selecting the next effective grid to execute the directional growth calculation.
The present invention provides an electronic device, including:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to at least one of the processors; wherein,
the memory stores instructions executable by at least one of the processors to enable the at least one of the processors to perform the method for reconstructing a wall rock fracture in a three-dimensional simulation test as described above.
The present invention provides a storage medium storing computer instructions for performing all the steps of the method for reconstructing a fracture of a surrounding rock in a three-dimensional simulation test as described above, when the computer instructions are executed by a computer.
According to the method, a fracture three-dimensional model is established by combining model parameters of a three-dimensional simulation model, the reconstruction of the model is realized, and finally an entity is manufactured by using a 3D printing technology to be compared with an experimental model, so that the problems of low monitoring precision of a surrounding rock fracture field and inaccuracy of the three-dimensional reconstruction model in a three-dimensional simulation experiment are solved. Meanwhile, a real-time image is transmitted by designing a mark point camera to obtain the development and expansion conditions of the fracture, so that the aim of monitoring the surrounding rock fracture field in a three-dimensional simulation test is fulfilled.
Drawings
FIG. 1 is a working flow chart of a method for reconstructing a surrounding rock fracture in a three-dimensional simulation test according to the invention;
FIG. 2 is a workflow diagram of a workflow diagram for generating a post-coal mining fracture digital model in a preferred embodiment of the present invention;
FIG. 3 is a schematic view of the growth direction in the grid according to the preferred embodiment of the present invention;
fig. 4 is a schematic diagram of a hardware structure of an electronic device according to the present invention.
Detailed Description
The following further describes embodiments of the present invention with reference to the accompanying drawings. In which like parts are designated by like reference numerals. It should be noted that the terms "front," "back," "left," "right," "upper" and "lower" used in the following description refer to directions in the drawings, and the terms "inner" and "outer" refer to directions toward and away from, respectively, the geometric center of a particular component.
Fig. 1 shows a work flow chart of a method for reconstructing a surrounding rock fracture in a three-dimensional simulation test according to the present invention, which includes:
s101, acquiring porosity of a coal mining model after simulated coal mining in a three-dimensional simulation experiment as the porosity after coal mining;
in S102, acquiring the distribution probability of solid phase growth nuclei statistically obtained by each grid on the coal mining model and the growth probability of the solid phase growth nuclei of each grid in multiple directions;
step S103, generating a fracture digital model after coal mining according to the porosity after coal mining, the distribution probability of solid phase growth nuclei of each grid and the growth probability of the solid phase growth nuclei of each grid in multiple directions;
and S104, outputting the coal mining fracture digital model to a three-dimensional printing device, and reconstructing the coal mining fracture digital model into a fracture three-dimensional model by the three-dimensional printing device.
Specifically, simulated coal mining is carried out on a coal mining model of the three-dimensional simulation experiment, and after the simulated coal mining, the porosity of the coal mining model is measured and calculated. The porosity after the simulated coal mining is then acquired as the porosity after coal mining by performing step S101. Then, the coal mining model is divided into a plurality of grids, and in step S102, the distribution probability of the solid phase growth sum and the growth probabilities in different directions of each grid are obtained. The distribution probability and the growth probability are obtained by adopting the existing statistical method on the measured data of the coal mining model.
Then, step S103 is executed to generate a coal mining fracture digital model, i.e., a three-dimensional porous medium model, and in step S104, the generated coal mining fracture digital model is output to a three-dimensional (3D) printing device, and the 3D printing device produces an entity to be compared with the coal mining model used in the experiment.
According to the method, a fracture three-dimensional model is established by combining model parameters of a three-dimensional simulation model, the reconstruction of the model is realized, and finally an entity is manufactured by using a 3D printing technology to be compared with an experimental model, so that the problems of low monitoring precision of a surrounding rock fracture field and inaccuracy of the three-dimensional reconstruction model in a three-dimensional simulation experiment are solved. Meanwhile, a real-time image is transmitted by designing a mark point camera to obtain the development and expansion conditions of the fracture, so that the aim of monitoring the surrounding rock fracture field in a three-dimensional simulation test is fulfilled.
In one embodiment:
in the three-dimensional simulation experiment, the porosity of the simulated coal mining model after coal mining is used as the porosity after coal mining, and the method specifically comprises the following steps: acquiring a coal mining model image of a simulated coal mining model in a three-dimensional simulation experiment, and calculating the porosity of the simulated coal mining model according to the coal mining model image to serve as the porosity after coal mining;
the method for acquiring the distribution probability of the solid phase growth nuclei obtained by each grid statistic on the coal mining model and the growth probability of the solid phase growth nuclei of each grid in multiple directions specifically comprises the following steps: dividing the coal mining model image into a plurality of grids, and according to the image information of the coal mining model image, counting the distribution probability of solid phase growth nuclei of each grid and the growth probability of the solid phase growth nuclei of each grid in a plurality of directions.
Specifically, images may be captured of the coal mining model, and the porosity after coal mining, the distribution probability of solid phase growth nuclei of each mesh, and the growth probability of solid phase growth nuclei of each mesh in a plurality of directions may be counted from the captured images.
According to the method and the device, the relevant parameters are determined by shooting the images of the coal mining model, so that the parameters are more accurate.
In one embodiment, the acquiring of the post-coal mining model image of the post-coal mining coal model in the three-dimensional simulation experiment specifically includes:
in the three-dimensional simulation experiment, images of a coal mining model after simulated coal mining, which are shot by a plurality of camera devices in the coal mining model, are obtained, and the images shot by the camera devices are fused to obtain a model image after coal mining.
Specifically, a plurality of high-definition camera devices are placed in a model of the simulation experiment platform. And after the model device is stable, performing simulated coal mining and collecting images. And then uploading the acquired images to a server, and fusing the images acquired by the plurality of camera devices by adopting the conventional image processing method to obtain a coal mining model image. Specifically, a binarization and edge detection method can be adopted to obtain a characteristic region outline of an original image, the pixel coordinates of the central point of the region are calculated through positioning, segmentation and decoding, and finally the difference values of different pixel points are converted into actual distances according to a scale, so that a model image after coal mining is obtained.
In one embodiment, the method further comprises the following steps:
identifying a fracture from the model image after coal mining as a real-time measurement fracture after coal mining, and calculating a fracture change value of the real-time measurement fracture after coal mining as a real-time measurement fracture change value after coal mining;
identifying a fracture from the coal mining fracture digital model as a fracture to be verified, and calculating a fracture change value to be verified of the fracture to be verified;
and comparing the actual measured fracture change value after coal mining with the fracture change value to be verified, if the difference value between the actual measured fracture change value after coal mining and the fracture change value to be verified is smaller than a preset difference threshold value, judging that the fracture digital model after coal mining is effective, otherwise, judging that the fracture digital model after coal mining is ineffective, and generating a new fracture digital model after coal mining according to the porosity after coal mining, the distribution probability of solid phase growth nuclei of each grid and the growth probability of the solid phase growth nuclei of each grid in multiple directions.
In the embodiment, the actual measured fracture change value after coal mining is compared with the fracture change value to be verified, so that whether the fracture digital model after coal mining is consistent with the actual coal mining model or not is judged, and the accuracy of the fracture digital model after coal mining is improved. The method for identifying the crack from the model image can be realized by adopting the existing crack identification method.
In one embodiment, the method further comprises the following steps:
acquiring the porosity of a coal mining model before coal mining as the porosity before coal mining;
dividing the coal mining model into a plurality of grids, and acquiring the initial distribution probability of solid phase growth nuclei of each grid and the initial growth probability of the solid phase growth nuclei of each grid in a plurality of directions;
generating an initial fracture digital model according to the porosity before coal mining, the initial distribution probability of the solid phase growth nuclei of each grid and the initial growth probability of the solid phase growth nuclei of each grid in multiple directions;
identifying a fracture from the initial fracture digital model as an initial fracture;
the calculating of the change value of the crack to be verified specifically comprises the following steps: and comparing the initial crack with the crack to be verified, and determining the change value of the crack to be verified.
In this embodiment, a fracture digital model is generated before coal mining, an initial fracture is identified from the fracture digital model generated before coal mining, and a fracture variation value is obtained by comparing the initial fracture with a fracture to be verified, where the fracture variation value includes, but is not limited to, fracture length, fracture width, or fracture depth.
In one embodiment, the method further comprises the following steps:
acquiring a pre-coal mining model image of a simulated pre-coal mining model in a three-dimensional simulation experiment;
identifying a fracture from the pre-coal mining model image as a pre-coal mining actual measurement fracture;
the method for calculating the fracture change value of the actual-measured fracture after coal mining as the actual-measured fracture change value after coal mining specifically comprises the following steps:
and comparing the coal mining actual measurement crack with the coal mining actual measurement crack to obtain a crack change value of the coal mining actual measurement crack as a coal mining actual measurement crack change value.
Specifically, if the coal mining model does not design a fracture before simulated coal mining, the actually measured fracture change value after coal mining is the relevant parameter of the actually measured fracture after coal mining. If the coal mining model is provided with a fracture before simulating coal mining, the actually measured fracture change value after coal mining is a comparison value of the related parameters of the actually measured fracture after coal mining and the actually measured fracture before coal mining.
In one embodiment, the method further comprises the following steps:
acquiring a coal mining process model image of a coal mining model in a simulated coal mining process in a three-dimensional simulation experiment;
and identifying a crack from the coal mining process model image, and generating and displaying change data of the crack about a time axis.
Specifically, the fracture variation data, including fracture length variation graphs and catastrophe trend graphs, are shown on the coordinate system in the order of the time axis. The length variation graph mainly shows the change process of the crack, and the sudden change trend graph mainly shows the change amplitude in the process.
The embodiment adds the change data of the fracture about the time axis, and is convenient for dynamically monitoring the three-dimensional simulation test surrounding rock fracture field.
In one embodiment, the generating a coal mining fracture digital model according to the porosity after coal mining, the distribution probability of solid phase growth nuclei of each grid, and the growth probability of solid phase growth nuclei of each grid in multiple directions specifically includes:
constructing a digital model based on the coordinates of the grid;
traversing all grids of the digital model, acquiring a grid node random number and the distribution probability of a solid phase growth core of each grid, if the grid node random number is smaller than the distribution probability of the solid phase growth core of each grid, generating a solid phase particle core by each grid at a corresponding coordinate position in the digital model, and generating a grid of the solid phase particle cores as an effective grid;
performing directional growth calculation on each effective grid until the porosity of the digital model reaches the porosity after coal mining;
forming contour equipotential lines based on an effective grid of a digital model, and generating a fracture digital model after coal mining based on the contour equipotential lines;
the calculation of the directional growth specifically comprises:
for an effective grid, acquiring a random number of peripheral nodes and the growth probability of a solid phase growth core of the effective grid in multiple directions, taking the direction in which the growth probability of the solid phase growth core of the effective grid is greater than the random number of the peripheral nodes as a growth direction, taking an adjacent grid positioned in the growth direction of the effective grid as a growth grid, generating a solid phase particle kernel at a corresponding coordinate position of the growth grid in the digital model, updating the growth grid generating the solid phase particle kernel into the effective grid, and calculating the porosity of the digital model, if the porosity of the digital model reaches the porosity after coal mining, ending the directional growth calculation, otherwise, selecting the next effective grid to execute the directional growth calculation.
Specifically, as shown in fig. 2, a workflow diagram for generating a coal mining fracture digital model is shown, which includes:
step S201, inputting parameters including porosity after coal mining, distribution probability of solid phase growth nuclei of each grid and growth probability of the solid phase growth nuclei of each grid in multiple directions;
step S202, traversing all grids;
step S203, if the random number of the grid nodes is smaller than the distribution probability of the solid phase growing kernel of the grid, the grid generates a solid phase particle kernel, otherwise, step S204 is executed;
for example, all the grids in the digital model are initialized to zero, and when the grids generate solid-phase particle kernels, positions corresponding to the grids in the digital model are set to be ones;
step S204, whether the traversal is finished or not is judged, if yes, the step S205 is executed, and if not, the step S202 is executed;
step S205, traversing all solid-phase particles and growing the solid-phase particles to the periphery;
step S206, if the random number of the peripheral nodes is smaller than the growth probability of the solid phase growth nucleus of the grid in a certain direction, taking the direction in which the growth probability of the solid phase growth nucleus of the grid is larger than the random number of the peripheral nodes as the growth direction, taking the adjacent grid positioned in the growth direction of the effective grid as the growth grid, and generating a solid phase particle kernel at the corresponding coordinate position of the growth grid in the digital model, otherwise, executing the step S205;
for example, when a grid generates a solid-phase particle kernel, a position corresponding to the grid in the digital model is set as one;
step S207, if the porosity of the digital model reaches the porosity after coal mining, executing step S208, otherwise executing step S205;
and step S208, generating a particle contour equipotential line according to the digital model, saving contour line data, generating a particle contour line graph, and saving the particle contour line graph as a DXF file.
Specifically, the coal mining model is divided into multiple layers, and each layer is divided into multiple grids. For example, the model image after coal mining is divided into multiple layers, and each layer is divided into multiple grids. The digital model also comprises a plurality of layers, and the grid of the coal mining model corresponds to the coordinate position of the corresponding layer in the digital model. As shown in fig. 3, each grid includes a growth direction 1, a growth direction 2, a growth direction 3, a growth direction 4, a growth direction 5, a growth direction 6, a growth direction 7, and a growth direction 8, wherein the growth directions 1 to 4 vary in length and the growth directions 5 to 8 vary in angle. The mesh in the digital model is initialized to zero. And when the grid is judged to generate the solid-phase particle kernel, placing one corresponding position of the grid in the digital model. When the porosity of the digital model reaches the post-coal porosity, a contour line of the solid phase particles is generated for each layer, for example, based on the coordinate position of one. And then programming the contour line data of all the layers through Auto CAD VBA, and sequentially overlapping and stacking the contour lines of all the layers to form a three-dimensional image. And finally, selecting ACIS (x. dxf) for file type output, and importing finite element analysis software to complete fracture reconstruction.
In addition, according to the porosity before coal mining, the initial distribution probability of the solid phase growth nuclei of each grid and the initial growth probability of the solid phase growth nuclei of each grid in multiple directions, the generated initial fracture digital model is consistent with the generated fracture digital model after coal mining, namely:
constructing a digital model based on the coordinates of the grid;
traversing all grids of the digital model, acquiring a grid node random number and the distribution probability of a solid phase growth core of each grid, if the grid node random number is smaller than the distribution probability of the solid phase growth core of each grid, generating a solid phase particle core by each grid at a corresponding coordinate position in the digital model, and generating a grid of the solid phase particle cores as an effective grid;
performing directional growth calculation on each effective grid until the porosity of the digital model reaches the porosity before coal mining;
forming contour equipotential lines based on an effective grid of a digital model, and generating a fracture digital model before coal mining based on the contour equipotential lines;
the calculation of the directional growth specifically comprises:
for an effective grid, acquiring a random number of peripheral nodes and the growth probability of a solid phase growth core of the effective grid in multiple directions, taking the direction in which the growth probability of the solid phase growth core of the effective grid is greater than the random number of the peripheral nodes as a growth direction, taking an adjacent grid positioned in the growth direction of the effective grid as a growth grid, generating a solid phase particle kernel at a corresponding coordinate position of the growth grid in the digital model, updating the growth grid generating the solid phase particle kernel into the effective grid, and calculating the porosity of the digital model, if the porosity of the digital model reaches the porosity before coal mining, ending the directional growth calculation, otherwise, selecting the next effective grid to execute the directional growth calculation.
As the best embodiment of the invention, the method for reconstructing the surrounding rock fracture in the three-dimensional simulation test comprises the following steps:
1. a plurality of high-definition camera devices are arranged in a model of the simulation experiment platform. And after the model device is stable, performing simulated coal mining and collecting images.
2. Uploading the collected image to a server, obtaining the characteristic area outline of the original image by adopting a binarization and edge detection method, calculating the pixel coordinate of the central point of the area through positioning, segmentation and decoding, and finally converting the difference value of different pixel points into an actual distance according to a scale so as to obtain the change value of the fracture.
3. And displaying the fracture change data including a fracture length change graph and an abrupt change trend graph in a coordinate system in a time axis sequence. The length variation graph mainly shows the change process of the crack, and the sudden change trend graph mainly shows the change amplitude in the process.
4. After the coal mining process is finished, measuring the porosity of the model, and obtaining the distribution probability P of the solid phase growing nucleus by each grid through actual measurement data statisticscdGiven probability of growth P in different directionsdi(i represents directions, and the values of 1, 2, 3, 4, 5, 6, 7 and 8 are shown in figure 3).
5. The entire computational domain is divided into multiple layers of grids, and a four parameter generation method (QSGS) generation program is programmed using Matlab according to the flow chart shown in fig. 3.
4. The porous medium image can be generated through a scatter diagram command scatter of Mat lab, the porous medium image is processed to generate a geometric model format which can be read by finite elements, contour lines of solid-phase particles are generated through a contour line command of Mat lab, and contour line data are stored in a data file. Because the contour line points are many, the data size is large, direct modeling is difficult by adopting general finite element analysis software, and Auto CAD VBA programming is selected to realize generation of the fracture model.
7. And finally, selecting ACIS (x. dxf) for file type output, and importing finite element analysis software to complete fracture reconstruction.
8. And (3) comparing the obtained image manufacturing entity with the original model by using a 3D printing technology, and verifying the accuracy of the three-dimensional reconstruction model.
Fig. 4 is a schematic diagram of a hardware structure of an electronic device according to the present invention, which includes:
at least one processor 401; and the number of the first and second groups,
a memory 402 communicatively coupled to at least one of the processors 401; wherein,
the memory 402 stores instructions executable by at least one of the processors to enable the at least one of the processors to perform the method for reconstructing a wall rock fracture in a three-dimensional simulation test as described above.
In fig. 4, one processor 401 is taken as an example.
The electronic device may further include: an input device 403 and a display device 404.
The processor 401, the memory 402, the input device 403, and the display device 404 may be connected by a bus or other means, and are illustrated as being connected by a bus.
The memory 402, which is a non-volatile computer-readable storage medium, may be used to store non-volatile software programs, non-volatile computer-executable programs, and modules, such as program instructions/modules corresponding to the method for reconstructing a surrounding rock fracture in a three-dimensional simulation test in the embodiment of the present application, for example, the method flow shown in fig. 1. The processor 401 executes various functional applications and data processing by running the nonvolatile software programs, instructions and modules stored in the memory 402, namely, the method for reconstructing the surrounding rock fracture in the three-dimensional simulation test in the above embodiments is implemented.
The memory 402 may include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required for at least one function; the storage data area may store data created according to the use of the surrounding rock fracture reconstruction method for the three-dimensional simulation test, and the like. Further, the memory 402 may include high speed random access memory, and may also include non-volatile memory, such as at least one magnetic disk storage device, flash memory device, or other non-volatile solid state storage device. In some embodiments, the memory 402 optionally includes a memory remotely located from the processor 401, and the remote memory may be connected via a network to a device for performing the method for reconstructing the wall rock fractures in the three-dimensional simulation test. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The input device 403 may receive input of user clicks and generate signal inputs related to user settings and functional controls for the method of reconstruction of the wall rock fractures in the three-dimensional simulation test. The display device 404 may include a display screen or the like.
The one or more modules stored in the memory 402, when executed by the one or more processors 401, perform the method for reconstructing a wall rock fracture in a three-dimensional simulation test in any of the method embodiments described above.
According to the method, a fracture three-dimensional model is established by combining model parameters of a three-dimensional simulation model, the reconstruction of the model is realized, and finally an entity is manufactured by using a 3D printing technology to be compared with an experimental model, so that the problems of low monitoring precision of a surrounding rock fracture field and inaccuracy of the three-dimensional reconstruction model in a three-dimensional simulation experiment are solved. Meanwhile, a real-time image is transmitted by designing a mark point camera to obtain the development and expansion conditions of the fracture, so that the aim of monitoring the surrounding rock fracture field in a three-dimensional simulation test is fulfilled.
An embodiment of the present invention provides a storage medium storing computer instructions for performing all the steps of the method for reconstructing a fracture of a surrounding rock in a three-dimensional simulation test as described above, when the computer executes the computer instructions.
The above-mentioned embodiments only express several embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the present invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (10)
1. A method for reconstructing a surrounding rock fracture in a three-dimensional simulation test is characterized by comprising the following steps:
acquiring porosity of a simulated coal mining model in a three-dimensional simulation experiment as the porosity after coal mining;
acquiring the distribution probability of solid phase growing nuclei statistically obtained by each grid on the coal mining model and the growing probability of the solid phase growing nuclei of each grid in multiple directions;
generating a fracture digital model after coal mining according to the porosity after coal mining, the distribution probability of solid phase growth nuclei of each grid and the growth probability of the solid phase growth nuclei of each grid in multiple directions;
and outputting the coal mining crack digital model to a three-dimensional printing device, and reconstructing the coal mining crack digital model into a crack three-dimensional model by the three-dimensional printing device.
2. The method for reconstructing the surrounding rock fractures in the three-dimensional simulation test according to claim 1, wherein the method comprises the following steps:
in the three-dimensional simulation experiment, the porosity of the simulated coal mining model after coal mining is used as the porosity after coal mining, and the method specifically comprises the following steps: acquiring a coal mining model image of a simulated coal mining model in a three-dimensional simulation experiment, and calculating the porosity of the simulated coal mining model according to the coal mining model image to serve as the porosity after coal mining;
the method for acquiring the distribution probability of the solid phase growth nuclei obtained by each grid statistic on the coal mining model and the growth probability of the solid phase growth nuclei of each grid in multiple directions specifically comprises the following steps: dividing the coal mining model image into a plurality of grids, and according to the image information of the coal mining model image, counting the distribution probability of solid phase growth nuclei of each grid and the growth probability of the solid phase growth nuclei of each grid in a plurality of directions.
3. The method for reconstructing the surrounding rock fracture in the three-dimensional simulation test according to claim 2, wherein the acquiring of the post-coal mining model image of the post-coal mining model in the three-dimensional simulation test specifically comprises:
in the three-dimensional simulation experiment, images of a coal mining model after simulated coal mining, which are shot by a plurality of camera devices in the coal mining model, are obtained, and the images shot by the camera devices are fused to obtain a model image after coal mining.
4. The method for reconstructing the surrounding rock fractures in the three-dimensional simulation test according to claim 3, further comprising the following steps of:
identifying a fracture from the model image after coal mining as a real-time measurement fracture after coal mining, and calculating a fracture change value of the real-time measurement fracture after coal mining as a real-time measurement fracture change value after coal mining;
identifying a fracture from the coal mining fracture digital model as a fracture to be verified, and calculating a fracture change value to be verified of the fracture to be verified;
and comparing the actual measured fracture change value after coal mining with the fracture change value to be verified, if the difference value between the actual measured fracture change value after coal mining and the fracture change value to be verified is smaller than a preset difference threshold value, judging that the fracture digital model after coal mining is effective, otherwise, judging that the fracture digital model after coal mining is ineffective, and generating a new fracture digital model after coal mining according to the porosity after coal mining, the distribution probability of solid phase growth nuclei of each grid and the growth probability of the solid phase growth nuclei of each grid in multiple directions.
5. The method for reconstructing the surrounding rock fractures in the three-dimensional simulation test according to claim 4, further comprising the following steps of:
acquiring the porosity of a coal mining model before coal mining as the porosity before coal mining;
dividing the coal mining model into a plurality of grids, and acquiring the initial distribution probability of solid phase growth nuclei of each grid and the initial growth probability of the solid phase growth nuclei of each grid in a plurality of directions;
generating an initial fracture digital model according to the porosity before coal mining, the initial distribution probability of the solid phase growth nuclei of each grid and the initial growth probability of the solid phase growth nuclei of each grid in multiple directions;
identifying a fracture from the initial fracture digital model as an initial fracture;
the calculating of the change value of the crack to be verified specifically comprises the following steps: and comparing the initial crack with the crack to be verified, and determining the change value of the crack to be verified.
6. The method for reconstructing the surrounding rock fractures in the three-dimensional simulation test according to claim 5, further comprising the following steps of:
acquiring a pre-coal mining model image of a simulated pre-coal mining model in a three-dimensional simulation experiment;
identifying a fracture from the pre-coal mining model image as a pre-coal mining actual measurement fracture;
the method for calculating the fracture change value of the actual-measured fracture after coal mining as the actual-measured fracture change value after coal mining specifically comprises the following steps:
and comparing the coal mining actual measurement crack with the coal mining actual measurement crack to obtain a crack change value of the coal mining actual measurement crack as a coal mining actual measurement crack change value.
7. The method for reconstructing the surrounding rock fractures in the three-dimensional simulation test according to claim 1, further comprising the following steps of:
acquiring a coal mining process model image of a coal mining model in a simulated coal mining process in a three-dimensional simulation experiment;
and identifying a crack from the coal mining process model image, and generating and displaying change data of the crack about a time axis.
8. The method for reconstructing the surrounding rock fractures in the three-dimensional simulation test according to claim 1, wherein the generating of the digital model of the fractures after coal mining according to the porosity after coal mining, the distribution probability of the solid phase growth nuclei of each grid, and the growth probability of the solid phase growth nuclei of each grid in multiple directions specifically comprises:
constructing a digital model based on the coordinates of the grid;
traversing all grids of the digital model, acquiring a grid node random number and the distribution probability of a solid phase growth core of each grid, if the grid node random number is smaller than the distribution probability of the solid phase growth core of each grid, generating a solid phase particle core by each grid at a corresponding coordinate position in the digital model, and generating a grid of the solid phase particle cores as an effective grid;
performing directional growth calculation on each effective grid until the porosity of the digital model reaches the porosity after coal mining;
forming contour equipotential lines based on an effective grid of a digital model, and generating a fracture digital model after coal mining based on the contour equipotential lines;
the calculation of the directional growth specifically comprises:
for an effective grid, acquiring a random number of peripheral nodes and the growth probability of a solid phase growth core of the effective grid in multiple directions, taking the direction in which the growth probability of the solid phase growth core of the effective grid is greater than the random number of the peripheral nodes as a growth direction, taking an adjacent grid positioned in the growth direction of the effective grid as a growth grid, generating a solid phase particle kernel at a corresponding coordinate position of the growth grid in the digital model, updating the growth grid generating the solid phase particle kernel into the effective grid, and calculating the porosity of the digital model, if the porosity of the digital model reaches the porosity after coal mining, ending the directional growth calculation, otherwise, selecting the next effective grid to execute the directional growth calculation.
9. An electronic device, comprising:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to at least one of the processors; wherein,
the memory stores instructions executable by at least one of the processors to enable the at least one processor to perform the method for reconstruction of a wall rock fracture in a three-dimensional simulation test as claimed in any one of claims 1 to 8.
10. A storage medium storing computer instructions for performing all the steps of the method for reconstructing a wall rock fracture in a three-dimensional simulation test according to any one of claims 1 to 8 when executed by a computer.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111509800.2A CN114187423A (en) | 2021-12-10 | 2021-12-10 | Surrounding rock fracture reconstruction method, electronic equipment and storage medium for three-dimensional simulation test |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111509800.2A CN114187423A (en) | 2021-12-10 | 2021-12-10 | Surrounding rock fracture reconstruction method, electronic equipment and storage medium for three-dimensional simulation test |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114187423A true CN114187423A (en) | 2022-03-15 |
Family
ID=80604459
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111509800.2A Pending CN114187423A (en) | 2021-12-10 | 2021-12-10 | Surrounding rock fracture reconstruction method, electronic equipment and storage medium for three-dimensional simulation test |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114187423A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116645602A (en) * | 2023-04-28 | 2023-08-25 | 长安大学 | Method, system and storage medium for intelligent identification and three-dimensional reconstruction of grotto weathered cracks |
-
2021
- 2021-12-10 CN CN202111509800.2A patent/CN114187423A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116645602A (en) * | 2023-04-28 | 2023-08-25 | 长安大学 | Method, system and storage medium for intelligent identification and three-dimensional reconstruction of grotto weathered cracks |
| CN116645602B (en) * | 2023-04-28 | 2024-05-07 | 长安大学 | Method, system and storage medium for intelligent identification and three-dimensional reconstruction of grotto weathered cracks |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Omar et al. | Towards an automated photogrammetry-based approach for monitoring and controlling construction site activities | |
| US7272264B2 (en) | System and method for hole filling in 3D models | |
| US9189862B2 (en) | Outline approximation for point cloud of building | |
| CN112633657B (en) | Construction quality management method, device, equipment and storage medium | |
| Truong-Hong et al. | Flying voxel method with Delaunay triangulation criterion for façade/feature detection for computation | |
| CN112991424B (en) | Fractal dimension calculation method and device based on octree algorithm | |
| CN113838005B (en) | Intelligent identification and three-dimensional reconstruction method and system for rock mass fracture based on dimension conversion | |
| Lai et al. | Surface roughness of rock faces through the curvature of triangulated meshes | |
| CN111870952B (en) | Altitude map generation method, device, equipment and storage medium | |
| TW201344633A (en) | System and method for measuring curved surfaces of objects | |
| CN111744183B (en) | Illumination sampling method and device in game and computer equipment | |
| CN112231874A (en) | Method and device for establishing underground pipeline model, computer equipment and storage medium | |
| CN113538668A (en) | Geological model construction method, geological model construction system, storage medium and electronic equipment | |
| CN114332134A (en) | Building facade extraction method and device based on dense point cloud | |
| CN114565722A (en) | Three-dimensional model monomer realization method | |
| CN114187423A (en) | Surrounding rock fracture reconstruction method, electronic equipment and storage medium for three-dimensional simulation test | |
| CN114821055A (en) | House model construction method and device, readable storage medium and electronic equipment | |
| Marino et al. | Measuring rock slope damage on rubble mound breakwater through digital photogrammetry | |
| Kim et al. | Automated two-dimensional geometric model reconstruction from point cloud data for construction quality inspection and maintenance | |
| CN111812670B (en) | Single photon laser radar space transformation noise judgment and filtering method and device | |
| CN117808852A (en) | Construction period digital monitoring method and system based on mapping point cloud | |
| Budianti et al. | Background blurring and removal for 3d modelling of cultural heritage objects | |
| Bueno et al. | Quantitative evaluation of CHT and GHT for column detection under different conditions of data quality | |
| Bajaj et al. | 3D Dense Reconstruction for Structural Defect Quantification | |
| CN116645602B (en) | Method, system and storage medium for intelligent identification and three-dimensional reconstruction of grotto weathered cracks |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |