CN115618700B - Mineral identification technology-based rock modeling method and system for containing physical surface - Google Patents

Abstract

The invention provides a method and a system for modeling a rock containing a physical surface based on a mineral identification technology, wherein the method comprises the following steps: manufacturing a rock sample slice; observing the rock sample slice, and acquiring a polarization microscope image of the rock sample slice with a clear observation of a schistosomie surface; processing the polarizing microscope image to obtain region area images of various minerals; vectorizing the area images of various minerals to obtain a complete rock area boundary vector diagram; importing the boundary vector diagram of the rock region into particle flow software, and enabling the Geometry corresponding to the imported boundary vector diagram of the rock region to coincide with a rock numerical model; by grouping mineral particles in the rock, setting different geometric, physical mechanical and contact parameters for the mineral particles with different components; and obtaining the numerical model of the rock containing the schistose surface by taking the change of the microscopic parameters of different minerals in the particle flow software. The method of the invention can accurately simulate the real schistose face shape of the rock in the nature.

Description

Mineral identification technology-based rock modeling method and system for containing physical surface

Technical Field

The invention relates to the technical field of rock mechanics, in particular to a method and a system for modeling a rock containing a physical surface based on a mineral identification technology.

Background

Schists are a common type of rock in rock mass engineering and are affected by regional metamorphism for a long time. The material contains a specific chip structure, so that the mechanical property of the material shows obvious anisotropy, and the influence on engineering safety is obvious. Therefore, the research on the sheet surface of the rock and the revealing of the physical mechanics law are important research contents of rock mechanics.

The particle flow software is numerical simulation software based on a discrete element method which is widely used in the field of rock engineering in recent years, and can carry out numerical simulation on rock materials. How to set a proper chip surface for simulating the chip structure of the rock is the basis of whether the data obtained by numerical simulation is accurate or not.

In the process of implementing the invention, the inventor finds that at least the following problems exist in the prior art:

at present, the modeling method for setting the rock schistosity surface in the particle flow software mainly comprises the following steps:

(1) The method comprises the steps of building rigid clusters in various shapes, randomly adding various rigid clusters into the built model in a random distribution mode, and simulating the schistose surface of the rock by utilizing the directional arrangement of the rigid clusters. In the method, because the shape and the size of the artificial component rigid cluster are specific and the type and the arrangement mode are limited, the real morphological characteristics of the rock physiological surface in a natural state cannot be accurately simulated.

(2) And cutting the rock model by using the strip joints with different inclination angles, and enabling the crack in the rock to expand towards the joint direction by utilizing the change of contact parameters of the joint, such as the change of normal bonding strength and tangential bonding strength. However, this method causes anisotropy of the rock through anisotropy of the joint, ignores mineral grains in the rock, and cannot depict the phenomenon of rock schistosity caused by directional development of the mineral grains.

(3) The particles in the whole model are divided according to the main mineral composition of the rock, the mode of directional arrangement of the particles in the mode of maximum mineral proportion is adopted to generate the chip surface, and then the model is filled in the mode of random distribution of other proportion minerals, so that the characteristic of the rock containing the chip surface is simulated. The method only focuses on the directional arrangement of particles with the most mineral proportion, and neglects the arrangement mode of other types of minerals.

Although the three modeling methods have characteristics, the real situation of the rock containing the physical surface in the indoor test cannot be accurately simulated. In a discrete element numerical simulation experiment about the rock containing the chip surface, how to establish a rock model capable of accurately reflecting the real directional arrangement condition of minerals in particle flow software is the key for determining the accuracy of the calculation result of the discrete element method.

Therefore, a need exists for a method and system for modeling rock containing physical surfaces based on mineral identification technology to at least partially solve the above technical problems.

Disclosure of Invention

In view of this, embodiments of the present invention provide a method and a system for modeling a rock containing a physical surface based on a mineral identification technology, so as to solve at least one of the problems in the prior art.

One aspect of the invention provides a mineral identification technology-based rock modeling method for a rock containing a geographical surface, which comprises the following steps:

collecting a rock sample with a sheet-like structure, and manufacturing a rock sample slice;

observing the rock sample slice through a polarization microscope, and acquiring a polarization microscope image of the rock sample slice with a clear observation of a schistose surface;

processing the polarization microscope image to obtain region area images of various minerals in the rock sample;

vectorizing the area images of various minerals to obtain a complete rock area boundary vector diagram;

establishing a rock numerical model through particle flow software, and determining geometric dimension parameters of the model;

the method comprises the steps that a geometer corresponding to a boundary vector diagram of a rock region is coincided with a rock numerical model by introducing the boundary vector diagram of the rock region into particle flow software;

the method comprises the steps of grouping mineral particles in the rock, setting different geometric, physical mechanics and contact parameters for the mineral particles with different components, and storing the different geometric, physical mechanics and contact parameters;

and finally obtaining the rock numerical model containing the chip surface by taking the change of different geometric, physical mechanics and contact parameters in the particle flow software.

In some embodiments of the invention, by collecting a rock sample in a sheet-like configuration and making a slice of the rock sample, comprises:

and respectively manufacturing the rock sample into rock sample slices in the north-south, east-west and horizontal orthogonal directions. And/or

The thickness of the rock sample slice is not more than 30um. This allows the rock sample slice to have good observation effect under a polarization microscope.

In some embodiments of the invention, the polarization microscope image is processed to obtain an area image of each type of mineral in the rock sample, including:

carrying out RGB color extraction treatment on the polarizing microscope image to obtain mineral component images of various minerals after the RGB color extraction treatment;

and carrying out gray processing on the mineral component images of the various minerals to obtain region area images of the various minerals.

In some embodiments of the invention, the method further comprises at least one of:

and carrying out tolerance rate processing on the region area images of the various minerals. Wherein, the regional area images of various minerals are processed by Tolerance (Tolerance) in order to make the junction of the boundary of the mineral region clearer and more accurate.

And coloring the area images of the various minerals. The coloring process is performed to facilitate the vectorization process of the subsequent image.

In some embodiments of the present invention, vectorizing the region area images of the various types of minerals to obtain a complete rock region boundary vector diagram, includes:

acquiring coordinate points on the boundary of the mineral area image based on the area image of each mineral;

forming coordinate points on the boundary into a data group, constructing a vector relation, and integrating in a vector space; according to the obtained vector data, deriving region boundary vector diagrams of various minerals;

and integrating and superposing the regional boundary vector diagrams of the various minerals to generate a complete rock regional boundary vector diagram.

In some embodiments of the invention, the method further comprises: and discarding repeated and similar vector data by adopting a lattice skipping method according to the obtained vector data in a close sequence. The effect of the processing is that the regional boundary vector diagram of the mineral can be smoothed, and the phenomenon that the boundary of the mineral is excessively hard due to excessive data is avoided.

In some embodiments of the invention, the rock numerical model is established by particle flow software, and the geometric dimension parameters of the model are determined, wherein the geometric dimension parameters comprise:

setting the range of the density and radius interval of a basic composition unit, namely particles, in the particle flow software according to the density and the particle size of the rock mineral;

and setting a boundary range of the rock sample, and randomly filling particles in the range, so that the particles are fully contacted.

In some embodiments of the present invention, the rock region boundary vector diagram is introduced into the particle flow software, and the Geometry corresponding to the introduced rock region boundary vector diagram coincides with the rock numerical model, where the specific content includes:

introducing the boundary vector diagram of the rock region into the particle flow software, converting the boundary vector diagram into Geometry, and determining the position relationship between the Geometry and a rock numerical model;

and positioning the Geometry corresponding to the boundary vector diagram of the rock region introduced into the particle flow software, and then scaling, translating and rotating the Geometry so that the Geometry corresponding to the boundary vector diagram of the introduced rock region coincides with the rock numerical model.

In some embodiments of the invention, by grouping mineral particles within a rock, different geometric, physico-mechanical and contact parameters are set for different compositions of mineral particles and stored, including:

dividing mineral particles in the rock into different components in a grouping mode according to a rock numerical model and a region divided by Geometry corresponding to the introduced rock region boundary vector diagram;

according to the geometrical morphological characteristics of the mineral represented by the geometrical morphological characteristics, different geometrical physical parameters are set for the mineral particles with different components;

setting different contact parameters for mineral particles of different components according to the physical and mechanical characteristics of the mineral represented by the contact parameters;

after the contact parameters are set, the Geometry for dividing different minerals is deleted, and all the geometric, physical mechanics and contact parameters of the rock numerical model are stored.

In another aspect of the present invention, a rock modeling system for a rock containing a geographical surface based on a mineral identification technology is provided, the system comprising:

the polarization microscope is used for observing the rock sample slice and acquiring a polarization microscope image of the rock sample slice with a clear observation of a schistose surface;

a memory for storing an executable program;

a processor for executing the executable program stored in the memory, causing the processor to perform actions comprising:

processing the polarization microscope image to obtain region area images of various minerals in the rock sample;

vectorizing the area images of various minerals to obtain a complete rock area boundary vector diagram;

establishing a rock numerical model through particle flow software, and determining geometric dimension parameters of the model;

the method comprises the steps that a Geometry corresponding to a boundary vector diagram of a rock region is coincided with a rock numerical model by leading the boundary vector diagram of the rock region into particle flow software;

the method comprises the steps of grouping mineral particles in the rock, setting different geometric, physical mechanics and contact parameters for the mineral particles with different components, and storing the different geometric, physical mechanics and contact parameters;

and finally obtaining the rock numerical model containing the chip surface by taking the change of different geometric, physical mechanics and contact parameters in the particle flow software.

According to the method and the system for modeling the rock containing the physiological surface based on the mineral identification technology, in the research process of the influence of the rock lamellar structure on the mechanical mechanism of the rock, compared with the traditional rock lamellar surface generation method adopted in discrete element software, the method can be used for setting the physiological surface of the rock of a specific type more accurately based on the mineral arrangement form of the rock body, the real physiological surface form of the rock in the nature can be accurately simulated, and the method and the system have a further promotion effect on the research of the micromechanical mechanism of the rock containing the physiological surface. In addition, compared with the traditional statistical-based method for randomly generating the chip surface, the chip surface of the rock is generated in the particle flow software by adopting the mineral identification-based technology, and a series of physical and mechanical properties reflected by a chip surface rock numerical model finally obtained by using the modeling method can more clearly represent a certain type of rock. Therefore, the method has higher universality and accuracy, and can improve the reliability and accuracy of the particle flow calculation result.

Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and drawings.

It will be appreciated by those skilled in the art that the objects and advantages that can be achieved with the present invention are not limited to the specific details set forth above, and that these and other objects that can be achieved with the present invention will be more clearly understood from the detailed description that follows.

Drawings

The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. The components in the figures are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the invention. For purposes of illustrating and describing some portions of the present invention, corresponding parts may be exaggerated in the drawings, i.e., may be larger relative to other components in an exemplary device actually made according to the present invention. In the drawings:

FIG. 1 is a flow chart of a method for modeling a rock across a face based on mineral identification technology in accordance with an embodiment of the present invention;

FIG. 2 is a partial flow chart of a method for modeling a rock across a face based on mineral identification technology in accordance with an embodiment of the present invention;

FIG. 3 is a partial flow chart of a method for modeling a rock across a face based on mineral identification in another embodiment of the present invention;

FIG. 4 is an image of a biotite quartz schist flake under a polarizing microscope in accordance with an embodiment of the present invention;

fig. 5 is a regional area image of the main minerals of the biotite quartz schist after RGB color-picking and gray-scale processing in an embodiment of the present invention, where the left side in the image is a regional area image of the biotite minerals and the right side is a regional area image of the quartz minerals;

FIG. 6 is an image of the area of a mineral region after coloring in accordance with an embodiment of the present invention;

fig. 7 is a mineral area boundary vector diagram in an embodiment of the present invention, in which the left side of the diagram is a region boundary vector diagram of a biotite mineral, and the right side is a region boundary vector diagram of a quartz mineral;

FIG. 8 is a regional boundary vector diagram of biotite quartz schist generated after mineral region boundary vector diagrams are overlaid in one embodiment of the present invention;

FIG. 9 is a rock numerical model built in the grain flow software in accordance with an embodiment of the present invention;

FIG. 10 is an image of a biotite quartz schist area boundary vector diagram integrated and superimposed with a rock numerical model in an embodiment of the present invention;

FIG. 11 is a graphical model of biotite quartz schist generated after grouping commands and setting the geometric, physical and mechanical and contact parameters of each mineral according to an embodiment of the present invention;

FIG. 12 is a schematic block diagram of a facies rock modeling system based on mineral identification techniques in an embodiment of the present invention.

Detailed Description

In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the following embodiments and accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.

It should be noted that, in order to avoid obscuring the present invention with unnecessary details, only the structures and/or processing steps closely related to the solution according to the present invention are shown in the drawings, and other details not so related to the present invention are omitted.

It should be emphasized that the term "comprises/comprising" when used herein, is taken to specify the presence of stated features, elements, steps or components, but does not preclude the presence or addition of one or more other features, elements, steps or components.

It is also noted herein that the term "coupled," if not specifically stated, may refer herein to not only a direct connection, but also an indirect connection in which an intermediate is present.

Hereinafter, embodiments of the present invention will be described with reference to the accompanying drawings. In the drawings, the same reference numerals denote the same or similar components, or the same or similar steps.

First, a method 100 for modeling a rock across a geographical surface based on a mineral identification technique according to an embodiment of the present application will be described with reference to fig. 1. As shown in fig. 1, a method 100 for modeling a rock containing a geographical surface based on a mineral identification technique may include the steps of:

in step S110, a rock sample having a sheet-like structure is collected, and a rock sample slice is produced.

In step S120, the rock sample sheet is observed by a polarization microscope, and a polarization microscope image of the rock sample sheet in which the schistose plane can be clearly observed is acquired.

In step S130, the polarization microscope image is processed to obtain area images of various types of minerals in the rock sample.

In step S140, vectorization processing is performed on the region area images of various types of minerals to obtain a complete rock region boundary vector diagram.

In step S150, a rock numerical model is established by the particle flow software, and geometric parameters of the model are determined.

In step S160, the Geometry corresponding to the boundary vector diagram of the rock region is overlapped with the rock numerical model by introducing the boundary vector diagram of the rock region into the particle flow software.

In step S170, by grouping the mineral particles in the rock, different geometric, physical mechanical and contact parameters are set for the mineral particles of different compositions, and the different geometric, physical mechanical and contact parameters are saved.

In step S180, a mathematical model of the rock containing the schistose face is finally obtained by taking the changes of the different geometric, physical mechanics and contact parameters in the particle flow software.

In the examples of the present application, based on mineral identification technology, first a rock sample containing sheet-like structures is collected and a slice of the rock sample is made. Then, the rock sample slice is observed by a polarization microscope, and a polarization microscope image of the rock sample slice with the schistose surface clearly observed is obtained. And processing the polarization microscope image to obtain the area image of each mineral in the rock sample. And then carrying out vectorization processing on the region area images of various minerals to obtain a complete rock region boundary vector diagram. And then establishing a rock numerical model through the particle flow software, determining geometric dimension parameters of the model, and leading the boundary vector diagram of the rock area into the particle flow software, so that the Geometry corresponding to the boundary vector diagram of the rock area is coincided with the rock numerical model. And then, by grouping the mineral particles in the rock, setting different geometric, physical and mechanical and contact parameters for the mineral particles with different components, and storing. And finally, obtaining the numerical model of the rock containing the chip surface by taking the change of different geometric, physical mechanics and contact parameters in the particle flow software.

As can be seen from the above description of the process, according to the mineral identification technology-based rock modeling method 100 for a rock slice plane in accordance with the embodiment of the present application, compared with the traditional rock slice plane generation method used in discrete element software, the method can be selected based on the mineral arrangement form of the rock itself, and set the slice plane of a specific type of rock more accurately, so as to accurately simulate the real slice plane form of the rock in nature. In addition, compared with the traditional statistical-based method for randomly generating the chip surface, the chip surface of the rock is generated in the particle flow software by adopting the mineral identification-based technology, and a series of physical and mechanical properties reflected by a chip surface rock numerical model finally obtained by using the modeling method can more clearly represent a certain type of rock. Therefore, the method has higher universality and accuracy.

The contents of the above steps of the method 100 for modeling a rock across a face based on mineral identification technology according to an embodiment of the present application will be described in detail below.

Containing schistose face rock, called schistose for short. Schist belongs to metamorphic rock-one of three types of rock: magma rock, sedimentary rock, metamorphic rock. Metamorphic rock refers to rock which is already in the earth crust and is subjected to the action of the internal force of the earth, the physicochemical properties of the rock are changed, and therefore the composition and the structural structure of the rock are changed. The common metamorphic effect is various, and the flaky structure in the schist is formed by the oriented arrangement of scaly, columnar and partially granular minerals under the metamorphic effect of a region due to complicated action factors such as abnormal heat flow in the region, pressure action and sometimes fluid phase addition. Schist is one of metamorphic rocks having a typical lamellar structure, and is a product of regional metamorphism. The mineral grain is characterized by having a flake structure, a plate structure and fibrous minerals which are arranged in parallel, have coarse grain sizes and can be distinguished by naked eyes. The common mineral is flaky mineral mica, and the granular mineral is mainly quartz and feldspar. Often contain characteristic metamorphic minerals such as andalusite, kyanite, garnet, cordierite, staurolite, echeverite and the like. Low strength, easy weathering and poor freezing resistance. For convenience of description in the embodiments of the present application, the method 100 for modeling the buccal rock based on the mineral identification technology in the embodiments of the present application is described below by taking biotite quartz schist as an example, but the method 100 for modeling the buccal rock based on the mineral identification technology in the embodiments of the present application is not intended to be applicable to only biotite quartz schist, and is also applicable to other schists.

In the embodiment of the present application, a rock sample with a sheet-like structure is collected in step S110, and a rock sample slice is produced. Taking the biotite quartz schist as an example, a biotite quartz schist sample can be collected in a strong structure activity area (namely, a schist forming area), such as China's Hetian area, and the rock sample has an obvious schist structure under visual observation. And preparing the collected rock sample into a cylindrical sample, and manufacturing a microscope slice according to the development trend of the schistose surface.

In order to ensure that the rock sample slice has good subsequent observation effect under a polarizing microscope, the rock sample can be respectively manufactured into the rock sample slice according to three orthogonal directions of north and south, east and west and horizontal, wherein the three orthogonal directions of north and south, east and west and horizontal are only exemplary directions set for obtaining good observation effect, and do not represent limitation to the rock sample slice, and the rock sample slice can be adaptively adjusted according to actual conditions. Further, it is preferable that the thickness of the rock-like thin section is not more than 30um, and for example, the thickness of the rock-like thin section may be limited to 20um.

In the embodiment of the present application, the rock sample sheet is observed by the polarization microscope in step S120, and a polarization microscope image of the rock sample sheet in which the schistose plane can be clearly observed is acquired. A complete image of the field coverage area covering the rock specimen sheet can be obtained, for example, by a laika DM4P smart polarization microscope, for better observation of the microstructure in the rock specimen sheet. A polarizing microscope image of a rock sample sheet having a physically oriented arrangement of distinct scale-like, columnar and partially granular minerals was observed and selected. As shown in fig. 4, the biotite quartz schist is taken as an example, which is an image of a biotite quartz schist flake under a polarization microscope.

In the embodiment of the present application, the polarization microscope image is processed in step S130 to obtain the area image of each type of mineral in the rock sample. Referring to fig. 1 and 2, step S130 may include the following steps:

in step S131, RGB color extraction is performed on the polarization microscope image to obtain mineral component images of various minerals after RGB color extraction.

In step S132, the mineral component images of the various types of minerals are subjected to gray scale processing to obtain area images of the various types of minerals.

In step S133, the tolerance ratio processing is performed on the region area images of the respective types of minerals.

In step S134, the region area images of the respective types of minerals are subjected to coloring processing.

Specifically, biotite quartz schist is still taken as an example. In step S131, RGB color sampling processing is performed on the polarization microscope image to obtain mineral component images of various minerals after RGB color sampling processing. Because the color gamut of biotite and quartz mineral is greatly different under the polarizing microscope (the color gamut of various minerals in the actual schist under the polarizing microscope is generally greatly different), the reflected RGB value range also has certain difference, so that various minerals in the rock can be respectively subjected to color taking treatment according to the RGB interval range of various minerals, thereby dividing different mineral components in the area of the rock, and obtaining mineral component images of various minerals subjected to RGB color taking treatment. According to the RGB interval range of the biotite mineral and the quartz mineral, the two main minerals of the biotite quartz schist are respectively subjected to color extraction treatment, different mineral components in the area of the area are divided, and mineral component images of the biotite and the quartz mineral subjected to RGB color extraction treatment are obtained. In step S132, the mineral component images of the various minerals are subjected to gray scale processing to obtain area images of the various minerals. The range of gray colors from pure white to pure black is specified to be 0-225 levels, 225 levels for white and 0 levels for black. The area image of the mineral of the type can be obtained by setting the gray value of the selected mineral to be 255 (white) and the gray values of the other minerals to be 0 (black). According to the method, the area image of the biotite mineral and the area image of the quartz mineral are sequentially obtained. As shown in fig. 5, the image is a regional area image of the main mineral of biotite quartz schist after RGB color matching and gray scale processing, and the left side of the image is a regional area image of biotite mineral and the right side is a regional area image of quartz mineral. In step S133, the area images of the various minerals are subjected to tolerance ratio processing, and the boundary between the biotite and the quartz mineral area is more accurate after the processing. In the step S134, the area images of various minerals are colored, and when the images of biotite minerals and quartz minerals can be colored respectively, the areas of the same minerals are filled with colors in a random color filling manner, so that the vectorization processing in the step S140 is performed on the area images of the areas conveniently. Referring to fig. 6, it is a mineral area image of the main mineral of biotite quartz schist after coloring treatment.

In the illustrated embodiment, although step S130 includes step S133 and step S134, it does not mean that step S130 necessarily includes step S133 and step S134, which is only one example of a preferred embodiment. Step S130 may not include any of step S133 and step S134, or may include only any of step S133 and step S134. In addition, in the illustrated embodiment, the order of step S133 is prior to step S134, which is also only an example and does not mean that the order of step S133 is necessarily prior to step S134, and the order of step S133 may be arranged after step S134.

In the embodiment of the present application, vectorization processing is performed on the region area images of various types of minerals in step S140 to obtain a complete rock region boundary vector diagram. Referring to fig. 1 and 3, step S140 may include the following steps:

in step S141, based on the region area images of each type of minerals, coordinate points on the boundary of the region area images of the minerals are acquired.

In step S142, the coordinate points on the boundary are formed into a data set, a vector relationship is constructed, and the data set is integrated in a vector space; and according to the obtained vector data, deriving the region boundary vector diagram of various minerals.

In step S143, repeated, similar vector data is discarded by using the trellis method by putting the resulting vector data in close order.

In step S144, the region boundary vector maps of the various types of minerals are integrated and superimposed to generate a complete rock region boundary vector map.

Specifically, biotite quartz schist is still taken as an example. In step S141, based on the area image of each type of mineral, a coordinate point on the boundary of the area image of the mineral area is obtained. Coordinate points on the boundaries of the images of the areas of biotite and quartz mineral areas can be quickly obtained by using Opencv software. In step S142, a large number of coordinate points on the boundary form each data set, a vector relationship is constructed, and integration is performed in a vector space; and (4) deriving a regional boundary vector diagram of the biotite and the quartz mineral according to the obtained vector data. Referring to fig. 7, it is a region boundary vector diagram of biotite and quartz minerals, the left side of the diagram is a region boundary vector diagram of biotite minerals, and the right side is a region boundary vector diagram of quartz minerals. In step S143, repeated and similar vector data are discarded by using the trellis-skipping method according to the obtained vector data in a close order (for example, according to the arrangement order, or in an inverse arrangement order). The lattice skipping method can simplify repeated and similar vector data (namely discard processing) by adopting a mode of 'saving one group and deleting one group', or a mode of 'saving two groups and deleting one group', so as to smooth the regional boundary vector diagram of the mineral and avoid transition hardening of the boundary of the mineral caused by excessive data. In step S144, the regional boundary vector diagrams of various minerals are integrated and superimposed to generate a complete rock regional boundary vector diagram. Taking the biotite quartz schist as an example, the area boundary vector diagram of the biotite mineral and the area boundary vector diagram of the quartz mineral are superposed and integrated to obtain the area boundary vector diagram of the biotite quartz schist. As shown in fig. 8, it is a regional boundary vector diagram of biotite quartz schist generated after superimposing a regional boundary vector diagram of biotite and quartz mineral.

In addition, in the illustrated embodiment, although step S140 includes step S143, it does not mean that step S140 necessarily includes step S143, which is only one example of a more preferred embodiment. Step S143 may not be included in step S140.

In the embodiment of the present application, a rock numerical model is established by the particle flow software in step S150, and geometric parameters of the model are determined. Step S150 specifically includes: setting the range of the density and radius interval of a basic composition unit, namely particles, in the particle flow software according to the density and particle size of rock minerals; and setting a boundary range of the rock sample, randomly filling particles in the range, and fully contacting the particles.

The biotite quartz schist is still used as an example for illustration. Firstly, establishing a rock numerical model of the biotite quartz schist in particle flow software, wherein the radius of particles is 0.02cm-0.04cm, the overall size of the rock numerical model is 10cm multiplied by 15cm, then randomly filling the particles in the range, and fully contacting the particles, and making a basis for grouping and setting contact parameters in the step S170. As shown in fig. 9, which is a rock numerical model built in the grain flow software. The particle flow software may be PFC, UDEC, or 3DEC, etc. invented by ITASCA. The particle flow software according to the embodiment of the present application is described by taking PFC software as an example, but the present invention is also applicable to other particle flow software.

In the embodiment of the present application, in step S160, the rock region boundary vector diagram is introduced into the particle flow software, and the Geometry corresponding to the introduced rock region boundary vector diagram coincides with the rock numerical model, and the content of the Geometry may include: importing the boundary vector diagram of the rock region into the particle flow software, converting the boundary vector diagram into Geometry, and determining the position relationship between the Geometry and the rock numerical model; and positioning the Geometry corresponding to the boundary vector diagram of the rock region introduced into the particle flow software, and then scaling, translating and rotating the Geometry so as to ensure that the Geometry corresponding to the boundary vector diagram of the introduced rock region is matched with the rock numerical model in a fitting manner.

The biotite quartz schist is taken as an example. And (3) introducing the regional boundary vector diagram of the biotite quartz schist into particle flow software, converting the regional boundary vector diagram into Geometry, and determining the position relationship between the Geometry and a rock numerical model. And (3) zooming, translating and rotating the Geometry corresponding to the regional boundary vector diagram of the imported biotite quartz schist so as to coincide with a rock numerical model generated in particle flow software. As shown in fig. 10, it is an image obtained by integrating and overlaying the biotite quartz schist area boundary vector diagram and the rock numerical model. The method comprises the steps of carrying out scaling, translation and rotation on Geometry corresponding to a region boundary vector diagram of imported biotite quartz schist, wherein the scaling, translation and rotation are all operations carried out for the coincidence of the Geometry and a rock numerical model generated in particle flow software, so that the scaling, translation and rotation are not all necessary operations, and only one step or two steps of the scaling, translation and rotation may be required according to actual conditions.

In the embodiment of the present application, in step S170, by grouping the mineral particles in the rock, different geometric, physical mechanical and contact parameters are set for the mineral particles of different compositions and stored. The specific content in step S170 may include: dividing mineral particles in the rock into different components in a grouping mode according to a rock numerical model and a region divided by Geometry corresponding to the introduced rock region boundary vector diagram; setting different geometrical physical parameters for mineral particles with different components according to the geometrical morphological characteristics of the mineral represented by the geometrical morphological characteristics; according to the physical and mechanical characteristics of the mineral represented by the contact parameters, different contact parameters are set for the mineral particles with different components; after the contact parameters are set, the Geometry for dividing different minerals is deleted, and various geometric, physical mechanics and contact parameters of the rock numerical model are saved.

Specifically, according to the region divided by Geometry generated by a rock numerical model and a region boundary vector diagram of the imported biotite quartz schist, the internal particles of the rock are divided into the components of the biotite minerals and the components of the quartz minerals in a grouping mode. Different geometrical and physical parameters are set for the mineral particles with different components according to the geometrical and morphological characteristics of the biotite mineral and the quartz mineral. Different contact parameters are set for mineral particles with different components according to the physical and mechanical characteristics of the biotite mineral and the quartz mineral. After the microscopic parameters are set, the Geometry for dividing different minerals is deleted, and the geometric, physical mechanics and contact parameters of the rock numerical model are stored. Referring to fig. 11, a numerical model of biotite quartz schist is generated after grouping and setting the geometric, physico-mechanical and contact parameters of the minerals.

In the embodiment of the application, the effect of generating the rock schistose surface is achieved by changing the values of the microscopic parameters of different minerals in the particle flow software in step S180, and finally the rock numerical model containing the schistose surface is obtained.

Based on the above description, according to the mineral identification technology-based facet-containing rock modeling method 100 in the embodiment of the present application, compared with the traditional rock slice facet generation method adopted in discrete element software, the method can set the facet of a specific type of rock more accurately based on the mineral arrangement form of the rock itself, and can accurately simulate the real facet form of the rock in nature. In addition, compared with the traditional statistical-based method for randomly generating the schistose surface, the method for generating the schistose surface of the rock in the particle flow software by adopting the mineral identification-based technology can more clearly represent a certain type of rock by using a series of physical and mechanical properties reflected by a schistose surface rock numerical model finally obtained by the modeling method. Therefore, the method has higher universality and accuracy.

The above exemplarily illustrates the method 100 for modeling the rock across the facies based on the mineral identification technology according to the embodiment of the present application. A rock modeling system 200 for a rock across a geographical surface based on mineral identification technology provided by another aspect of the present application is described below in conjunction with fig. 12.

An example modeling system 200 for implementing a method of modeling a faceted rock based on mineral identification techniques according to an embodiment of the present invention is described with reference to FIG. 12. The system 200 includes at least one polarization microscope 210 and at least one example electronic device 220. The connection relationship between the polarization microscope 210 and the electronic device 220 may be only that the polarization microscope 210 transmits the acquired polarization microscope image of the rock sample slice with the clearly observable schistose surface to the electronic device 220, or the electronic device 220 may control the polarization microscope 210 to be turned on and off, control the polarization microscope 210 to acquire the time and frequency of the polarization microscope image of the rock sample slice with the clearly observable schistose surface, and control the polarization microscope 210 to transmit the polarization microscope image of the rock sample slice with the clearly observable schistose surface to the electronic device 220 at a set time.

As shown in fig. 12, the polarization microscope 210 is used for observing a rock sample slice and acquiring a polarization microscope image of the rock sample slice in which a physiological plane can be clearly observed.

The electronic device 220 may include one or more processors 221, one or more memories 222, an input device 223, and an output device 224, which may be interconnected by a bus system 225 and/or other form of connection mechanism (not shown). It should be noted that the components and configuration of the electronic device 220 shown in fig. 12 are exemplary only, and not limiting, and the electronic device may have other components and configurations as desired.

The processor 221 may be a Central Processing Unit (CPU) or other form of processing unit having data processing capabilities and/or instruction execution capabilities, and may control other components in the electronic device 220 to perform desired functions.

The memory 222 may include one or more computer program products that may include various forms of computer-readable storage media, such as volatile memory and/or non-volatile memory. The volatile memory may include, for example, random Access Memory (RAM), cache memory (cache), and/or the like. The non-volatile memory may include, for example, read Only Memory (ROM), hard disk, flash memory, etc. On which one or more computer program instructions may be stored that may be executed by processor 221 to implement client-side functionality (implemented by the processor) and/or other desired functionality in embodiments of the invention described herein. Various applications and various data, such as various data used and/or generated by the applications, may also be stored in the computer-readable storage medium.

The input device 223 may be a device used by a user to input instructions, and may include one or more of a keyboard, a mouse, a microphone, a touch screen, and the like. In addition, the input device 223 may be any interface for receiving information.

The output device 224 may output various information (e.g., images or sounds) to an outside (e.g., a user), and may include one or more of a display, a speaker, and the like. The output device 224 may be any other device having an output function.

For example, the example electronic device 220 for implementing the method 100 for modeling the rocks containing the facies according to the mineral identification technology of the embodiment of the present invention may be applied to electronic devices such as a terminal device (e.g., a mobile phone), a tablet computer, a notebook computer, an ultra-mobile personal computer (UMPC), a handheld computer, a netbook, a Personal Digital Assistant (PDA), a wearable device (e.g., a smart watch, a smart glasses, or a smart helmet), an Augmented Reality (AR) \\ Virtual Reality (VR) device, a smart home device, and a vehicle-mounted computer, and the embodiment of the present invention does not limit this.

Referring to fig. 12, the electronic device 220 in the mineral identification technology-based multi-faceted rock modeling system 200 according to an embodiment of the present application includes a processor 221 and a memory 222, the memory 222 storing an executable program executed by the processor 221, the executable program, when executed by the processor 221, causing the processor 221 to perform the aforementioned multi-faceted rock modeling method 100 according to an embodiment of the present application based on mineral identification technology. The specific operation of the mineral identification technology-based faceted rock modeling system according to the embodiments of the present application can be understood by those skilled in the art with reference to the foregoing description, and for the sake of brevity, specific details are not repeated here, and only some main operations of the processor 221 are described.

In one embodiment of the application, the executable program, when executed by the processor 221, causes the processor 221 to perform the steps of: processing the polarizing microscope image to obtain region area images of various minerals in the rock sample; vectorizing the region area images of various minerals to obtain a complete rock region boundary vector diagram; establishing a rock numerical model through particle flow software, and determining geometric dimension parameters of the model; introducing a rock region boundary vector diagram into particle flow software, and enabling a Geometry corresponding to the introduced rock region boundary vector diagram to coincide with a rock numerical model; mineral particles in the rock are grouped, and different geometric, physical mechanics and contact parameters are set for the mineral particles with different components and are stored; and finally obtaining the rock numerical model containing the chip surface by taking the change of different geometric, physical mechanics and contact parameters in the particle flow software.

In one embodiment of the application, the executable program, when executed by the processor 221, causes the processor 221 to further perform the steps of: carrying out RGB color extraction treatment on the polarizing microscope image to obtain mineral component images of various minerals after the RGB color extraction treatment; and carrying out gray level processing on the mineral component images of various minerals to obtain region area images of various minerals.

In one embodiment of the application, the executable program, when executed by the processor 221, causes the processor 221 to further perform the steps of: carrying out tolerance rate processing on the region area images of various minerals; and coloring the area images of various minerals.

In one embodiment of the application, the executable program, when executed by the processor 221, causes the processor 221 to further perform the steps of: acquiring coordinate points on the boundary of the mineral area image based on the area image of each mineral; forming coordinate points on the boundary into a data group, constructing a vector relation, and integrating in a vector space; according to the obtained vector data, deriving region boundary vector diagrams of various minerals; and integrating and superposing the regional boundary vector diagrams of various minerals to generate a complete rock regional boundary vector diagram.

In one embodiment of the application, the executable program, when executed by the processor 221, causes the processor 221 to further perform the steps of: the repeated and similar vector data are discarded by adopting a lattice skipping method according to the similar sequence of the obtained vector data.

In one embodiment of the application, the executable program, when executed by the processor 221, causes the processor 221 to further perform the steps of: setting the range of the density and radius interval of a basic composition unit, namely particles, in the particle flow software according to the density and the particle size of the rock mineral; and setting a boundary range of the rock sample, and randomly filling particles in the range, so that the particles are fully contacted.

In one embodiment of the application, the executable program, when executed by the processor 221, causes the processor 221 to further perform the steps of: importing the boundary vector diagram of the rock region into particle flow software, converting the boundary vector diagram into Geometry, and determining the position relationship between the Geometry and a rock numerical model; and positioning the Geometry corresponding to the boundary vector diagram of the rock region introduced into the particle flow software, and then scaling, translating and rotating the Geometry so as to ensure that the Geometry corresponding to the boundary vector diagram of the introduced rock region is matched with the rock numerical model in a matching way.

In one embodiment of the application, the executable program, when executed by the processor 221, causes the processor 221 to further perform the steps of: dividing mineral particles in the rock into different components in a grouping mode according to a rock numerical model and a region divided by Geometry corresponding to a guided rock region boundary vector diagram; according to the geometrical morphological characteristics of the mineral represented by the geometrical morphological characteristics, different geometrical physical parameters are set for the mineral particles with different components; according to the physical and mechanical characteristics of the mineral represented by the contact parameters, different contact parameters are set for the mineral particles with different components; after the contact parameters are set, the Geometry for dividing different minerals is deleted, and various geometric, physical mechanics and contact parameters of the rock numerical model are saved.

Furthermore, according to an embodiment of the present application, there is also provided a storage medium on which a computer program is stored, which when executed by a processor is configured to execute the corresponding steps of the method for modeling the facies rocks based on the mineral identification technology according to an embodiment of the present application. The storage medium may include, for example, a memory card of a smart phone, a storage component of a tablet computer, a hard disk of a personal computer, a Read Only Memory (ROM), an Erasable Programmable Read Only Memory (EPROM), a portable compact disc read only memory (CD-ROM), a USB memory, or any combination of the above storage media. The computer-readable storage medium may be any combination of one or more computer-readable storage media.

Based on the above description, according to the method and the system for modeling the rock containing the physical face based on the mineral identification technology in the embodiment of the application, compared with the traditional method for generating the physical face of the rock slice adopted in the discrete element software, the method can set the physical face of the rock of a specific type more accurately based on the mineral arrangement form of the rock body, and can accurately simulate the real physical face form of the rock in the nature. In addition, compared with the traditional statistical-based method for randomly generating the chip surface, the chip surface of the rock is generated in the particle flow software by adopting the technology based on mineral identification, and a series of physical and mechanical properties reflected by a chip surface rock numerical model finally obtained by using the modeling method can more clearly represent a certain type of rock, so that the method has higher universality and accuracy.

Although the example embodiments have been described herein with reference to the accompanying drawings, it is to be understood that the above-described example embodiments are merely illustrative and are not intended to limit the scope of the present application thereto. Various changes and modifications may be effected therein by one of ordinary skill in the pertinent art without departing from the scope or spirit of the present application. All such changes and modifications are intended to be included within the scope of the present application as claimed in the appended claims.

Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.

In the several embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other manners. For example, the above-described device embodiments are merely illustrative, and for example, the division of the units is only one logical functional division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another device, or some features may be omitted, or not executed.

In the description provided herein, numerous specific details are set forth. However, it is understood that embodiments of the application may be practiced without these specific details. In some instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure an understanding of this description.

Similarly, it should be appreciated that in the description of exemplary embodiments of the present application, various features of the present application are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the application and aiding in the understanding of one or more of the various inventive aspects. However, the method of the present application should not be construed to reflect the intent: this application is intended to cover such departures from the present disclosure as come within known or customary practice in the art to which this invention pertains. Rather, as the following claims reflect, inventive aspects lie in less than all features of a single disclosed embodiment. Thus, the claims following the detailed description are hereby expressly incorporated into this detailed description, with each claim standing on its own as a separate embodiment of this application.

It will be understood by those skilled in the art that all of the features disclosed in this specification (including any accompanying claims, abstract and drawings), and all of the processes or elements of any method or apparatus so disclosed, may be combined in any combination, except combinations where such features are mutually exclusive. Each feature disclosed in this specification (including any accompanying claims, abstract and drawings) may be replaced by alternative features serving the same, equivalent or similar purpose, unless expressly stated otherwise.

Furthermore, those skilled in the art will appreciate that while some embodiments described herein include some features included in other embodiments, rather than other features, combinations of features of different embodiments are meant to be within the scope of the application and form different embodiments. For example, in the claims, any of the claimed embodiments may be used in any combination.

The various component embodiments of the present application may be implemented in hardware, or in software modules running on one or more processors, or in a combination thereof. Those skilled in the art will appreciate that a microprocessor or Digital Signal Processor (DSP) may be used in practice to implement some or all of the functionality of some of the modules according to embodiments of the present application. The present application may also be embodied as apparatus programs (e.g., computer programs and computer program products) for performing a portion or all of the methods described herein. Such programs implementing the present application may be stored on a computer readable medium or may be in the form of one or more signals. Such a signal may be downloaded from an internet website or provided on a carrier signal or in any other form.

It should be noted that the above-mentioned embodiments illustrate rather than limit the application, and that those skilled in the art will be able to design alternative embodiments without departing from the scope of the appended claims. In the claims, any reference signs placed between parentheses shall not be construed as limiting the claim. The word "comprising" does not exclude the presence of elements or steps not listed in a claim. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements. The application may be implemented by means of hardware comprising several distinct elements, and by means of a suitably programmed computer. In the unit claims enumerating several means, several of these means can be embodied by one and the same item of hardware. The usage of the words first, second and third, etcetera do not indicate any ordering. These words may be interpreted as names.

The above description is only for the specific embodiments of the present application or descriptions thereof, and the protection scope of the present application is not limited thereto, and any person skilled in the art can easily think of the changes or substitutions within the technical scope of the present application, and should be covered by the protection scope of the present application. The protection scope of the present application shall be subject to the protection scope of the claims.

Claims (10)

1. A mineral identification technology-based method for modeling a rock containing a physical surface is characterized by comprising the following steps:

collecting a rock sample with a sheet-like structure, and manufacturing a rock sample slice;

observing the rock sample slice through a polarization microscope, and acquiring a polarization microscope image of the rock sample slice with a clear observation of a schistose surface;

processing the polarization microscope image to obtain region area images of various minerals in the rock sample;

vectorizing the area images of various minerals to obtain a complete rock area boundary vector diagram;

establishing a rock numerical model through particle flow software, and determining geometric dimension parameters of the model;

the method comprises the steps that a Geometry corresponding to a boundary vector diagram of a rock region is coincided with a rock numerical model by leading the boundary vector diagram of the rock region into particle flow software;

the method comprises the steps of grouping mineral particles in the rock, setting different geometric, physical mechanics and contact parameters for the mineral particles with different components, and storing the different geometric, physical mechanics and contact parameters;

and finally obtaining the rock numerical model containing the chip surface by taking the change of different geometric, physical mechanics and contact parameters in the particle flow software.

2. The mineral identification technology-based rock modeling method for the geographical surface of the rock as claimed in claim 1, wherein the step of collecting rock samples with a sheet-like structure and making rock sample slices comprises:

respectively manufacturing the rock sample into rock sample slices in the north-south, east-west and horizontal orthogonal directions; and/or

The thickness of the rock sample slice is not more than 30um.

3. The method for modeling a faceted rock based on mineral identification technology of claim 1, wherein processing the polarization microscope image to obtain a regional area image of each type of mineral in the rock sample comprises:

carrying out RGB color extraction treatment on the polarizing microscope image to obtain mineral component images of various minerals after the RGB color extraction treatment;

and carrying out gray level processing on the mineral component images of the various minerals to obtain region area images of the various minerals.

4. The mineral identification technology-based method for modeling a faceted rock according to claim 1 or 3, wherein said method further comprises at least one of:

carrying out tolerance rate processing on the region area images of the various minerals;

and coloring the area images of the various minerals.

5. The mineral identification technology-based geographical face-specific rock modeling method according to claim 1, wherein vectorizing the region area images of each type of mineral to obtain a complete rock region boundary vector diagram comprises:

acquiring coordinate points on the boundary of the mineral area image based on the area image of each mineral;

forming coordinate points on the boundary into a data group, constructing a vector relation, and integrating in a vector space; according to the obtained vector data, deriving region boundary vector diagrams of various minerals;

and integrating and superposing the regional boundary vector diagrams of the various minerals to generate a complete rock regional boundary vector diagram.

6. The mineral identification technology-based method for modeling a geographical face-bearing rock according to claim 5, further comprising:

and discarding repeated and similar vector data by adopting a lattice skipping method according to the obtained vector data in a close sequence.

7. The mineral identification technology-based method for modeling the rock with the physical face according to claim 1, wherein the rock numerical model is established through particle flow software, and the determination of the geometric size parameters of the model comprises the following steps:

setting the range of the density and radius interval of a basic composition unit, namely particles, in the particle flow software according to the density and particle size of rock minerals;

and setting a boundary range of the rock sample, and randomly filling particles in the range, so that the particles are fully contacted.

8. The method for modeling a facies-containing rock based on a mineral identification technology as claimed in claim 1 wherein the step of registering the Geometry corresponding to the introduced rock region boundary vector diagram with a rock numerical model by introducing the rock region boundary vector diagram into the particle flow software comprises:

introducing the boundary vector diagram of the rock region into the particle flow software, converting the boundary vector diagram into Geometry, and determining the position relationship between the Geometry and a rock numerical model;

and positioning the Geometry corresponding to the boundary vector diagram of the rock region introduced into the particle flow software, and then scaling, translating and rotating the Geometry so that the Geometry corresponding to the boundary vector diagram of the introduced rock region coincides with the rock numerical model.

9. The mineral identification technology-based rock modeling method for the geographical surface of the rock as claimed in claim 8, wherein the method comprises the steps of grouping mineral particles in the rock, setting different geometric, physical mechanical and contact parameters for the mineral particles of different compositions, and storing the different geometric, physical mechanical and contact parameters, and comprises:

dividing mineral particles in the rock into different components in a grouping mode according to a rock numerical model and a region divided by Geometry corresponding to the introduced rock region boundary vector diagram;

according to the geometrical morphological characteristics of the mineral represented by the geometrical morphological characteristics, different geometrical physical parameters are set for the mineral particles with different components;

setting different contact parameters for mineral particles of different components according to the physical and mechanical characteristics of the mineral represented by the contact parameters;

after the contact parameters are set, the Geometry for dividing different minerals is deleted, and various geometric, physical mechanics and contact parameters of the rock numerical model are saved.

10. A rock modeling system for containing physical surfaces based on mineral identification technology, the system comprising:

the polarization microscope is used for observing the rock sample slice and acquiring a polarization microscope image of the rock sample slice with a clear observation of a schistose surface;

a memory for storing an executable program;

a processor for executing the executable program stored in the memory, causing the processor to perform actions comprising:

processing the polarization microscope image to obtain region area images of various minerals in the rock sample;

vectorizing the area images of various minerals to obtain a complete rock area boundary vector diagram;

establishing a rock numerical model through particle flow software, and determining geometric dimension parameters of the model;

the method comprises the steps that a Geometry corresponding to a boundary vector diagram of a rock region is coincided with a rock numerical model by leading the boundary vector diagram of the rock region into particle flow software;

the method comprises the steps of grouping mineral particles in the rock, setting different geometric, physical mechanics and contact parameters for the mineral particles with different components, and storing the different geometric, physical mechanics and contact parameters;

and finally obtaining the rock numerical model containing the chip surface by taking the change of different geometric, physical mechanics and contact parameters in the particle flow software.

Images