Disclosure of Invention
The invention aims to solve the problems and provides a fault geometric modeling method for an intelligent mine, which can store and display fault data.
In a first aspect, the invention provides a fault geometric modeling method for an intelligent mine, which comprises the following steps:
1) Data preparation, namely surveying and acquiring fault position data and physical and mechanical attribute data;
2) The stratum modeling comprises the steps of extracting middle plane thickness data in GIS data, shifting middle plane nodes, wherein the offset is half of the thickness to form an upper boundary and a lower boundary;
3) Establishing a region on a horizontal plane, enabling projections of all stratum in GIS data on the horizontal plane to fall inside the region, and arranging uniform Delaunay triangulation grids in the region; because the grids and nodes of the stratum interface offset in the step are not in one-to-one correspondence in the vertical direction, interpolation is needed to be carried out on the stratum interface at the corresponding position according to the plane coordinates of the newly built grid nodes, and a new stratum interface is established in the process, wherein the interface has uniform control point grids and can be aligned in the vertical direction. The new stratum interface still does not have the stratum topological structure, and the next step is to determine the stratum interface on the basis of the new stratum interface so that the geometric boundaries of the upper stratum and the lower stratum form a joint relation.
4) Extracting stratum interfaces to establish stratum, namely interpolating the stratum interfaces according to the sequence from top to bottom for each grid node coordinate arranged on a horizontal plane, wherein the bottom plate of an upper stratum and the top plate of a lower stratum in all vertically adjacent two stratum adopt common nodes to form a final stratum interface, and finally connecting the interfaces belonging to the same stratum to express each stratum through a triangular prism set;
5) The method comprises the steps of adding fault units, generating triangular grids with uniform surface grids inwards according to geometrical information of faults, giving distances between layers and fault interaction information to the surface grid units according to survey data after the generation is finished, wherein the fault interaction information comprises slip directions and slip distance information on two sides of the faults, the information represents the fault to belong to the categories of forward faults, reverse faults or walk-slip faults, grid data comprise node information and triangular unit information, each triangular grid is formed by connecting three nodes, the data of the faults are stored in two storage formats, one of the data is stored on each node, namely, the given attribute is unique to a certain node, the geometrical parameter is given to the node, the other data is stored on the triangular grid, namely, the given attribute is unique to a certain triangular grid, and the physical property parameter is given to the grid.
Further, the geometrical modeling method of the fault oriented to the intelligent mine comprises the following types of data, namely a plane formed by three section point coordinates, known fault occurrence information and a plane passing through a certain point, and a curved surface formed by fault lines and fault inclination angle data.
Furthermore, according to the fault geometric modeling method for the intelligent mine, the surface grid unit and the body unit of the whole solving domain are mutually independent, and the independent fault grid is adopted, so that a complex geometric calculation process in geometric modeling is omitted, and a strict topological structure is not required. By adopting the method, a data storage format of faults in stratum geometry is established, and basic data is provided for stratum grid division and simulation calculation.
Further, according to the fault geometric modeling method for the intelligent mine, the data stored in the surface grid comprise but are not limited to geometric parameters including thickness and normal direction, and physical parameters including Young modulus, poisson's ratio, cohesive force, internal friction angle and permeability coefficient.
In a second aspect, the invention provides a fault geometric modeling system for an intelligent mine, which comprises an acquisition module, a modeling module, a stratum construction module and a fault data adding module;
The acquisition module is used for surveying and acquiring fault position data and physical and mechanical attribute data;
The modeling module is used for extracting middle plane thickness data in GIS data, offsetting the middle plane nodes, wherein the offset is half of the thickness to form an upper boundary and a lower boundary;
The stratum construction module is used for establishing an area on a horizontal plane, enabling projections of all stratum in GIS data on the horizontal plane to fall inside the area, arranging uniform Delaunay triangulation grids in the area, interpolating stratum interfaces according to the sequence from top to bottom for each grid node coordinate arranged on the horizontal plane, forming a final stratum interface by adopting common nodes on a bottom plate of an upper stratum and a top plate of a lower stratum in all vertically adjacent two strata, and finally connecting the interfaces belonging to the same stratum and expressing each stratum through a triangular prism set;
The fault data adding module is used for generating triangular grids with uniform surface grids inwards according to geometric information of faults, endowing distances among fault layers and fault interaction information on surface grid units according to survey data after the generation is finished, wherein the fault interaction information comprises slip directions and dislocation distance information on two sides of the faults, the information characterizes the faults to belong to categories of forward faults, reverse faults or walk-slip faults, grid data comprise node information and triangle unit information, each triangular grid is formed by connecting three nodes, the data of the faults are stored in two storage formats, one storage format is stored on each node, namely the endowed attribute is unique to a certain node, the geometric parameter is endowed on the node, the other storage format is stored on the triangular grid, namely the endowed attribute is unique to a certain triangular grid, and the physical parameter is endowed on the grid.
Further, the geometrical modeling system of the fault oriented to the intelligent mine comprises the following type data, namely a plane formed by three section point coordinates, known fault occurrence information and a plane passing through a certain point, and a curved surface formed by fault lines and fault inclination angle data.
Furthermore, according to the fault geometric modeling system for the intelligent mine, the surface grid unit and the body unit of the whole solving domain are mutually independent, and the independent fault grid is adopted, so that a complex geometric calculation process in geometric modeling is omitted, and a strict topological structure is not required. By adopting the method, a data storage format of faults in stratum geometry is established, and basic data is provided for stratum grid division and simulation calculation.
Further, the geometrical modeling system of fault oriented to intelligent mine comprises geometrical parameters including thickness, normal direction, physical parameters including Young modulus, poisson's ratio, cohesion, internal friction angle and permeability coefficient.
In a third aspect, the invention provides a fault geometry modeling device for an intelligent mine, comprising a memory and a processor, wherein the memory is used for storing a computer program, and the processor is used for realizing the fault geometry modeling method for the intelligent mine according to the first aspect when the computer program is executed.
In a fourth aspect, the present invention provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a method for geometrical modeling of a fault for a smart mine according to the first aspect.
According to the fault geometric modeling method for the intelligent mine, fault data are integrated into a geometric model based on the existing stratum geometric generation method, so that the problem that storage and display of fault data are lacking in the prior art is solved, compared with the prior art, the fault geometric modeling method not only comprises space position information of the fault, but also comprises attribute information of the fault at different space positions, due to the fact that the fault is not geometrically associated, a complex geometric calculation process is omitted, modeling efficiency is remarkably improved, the stored fault data are used as bases for grid division and grid attribute giving, more detailed geological data storage can be provided for a mine digital model, and accurate data support is provided for intelligent mine simulation.
Detailed Description
The geometrical modeling method of fault oriented to intelligent mine is described in detail below through the drawings and the embodiments.
Example 1
The embodiment discloses a fault geometric modeling method for intelligent mines, which comprises the following steps as shown in fig. 1:
1) The method comprises the steps of data preparation, namely surveying and obtaining position data and physical and mechanical attribute data of faults, wherein GIS data comprise position and thickness information of strata, each stratum consists of a triangulation plane with thickness attribute, and different stratum data are independent from each other and have no geometric topological relation. The data of the middle plane is stored by adopting a plane geographic coordinate and an elevation coordinate, the space position of the middle plane of a certain stratum is represented, and each point data stores the stratum thickness of the stratum under the geographic coordinate of the point;
2) The stratum modeling comprises the steps of extracting middle plane thickness data in GIS data, shifting middle plane nodes, wherein the offset is half of the thickness to form an upper boundary and a lower boundary;
3) The stratum interface reconstruction comprises the steps of establishing a region on a horizontal plane, enabling projections of all stratum in GIS data on the horizontal plane to fall inside the region, and arranging uniform Delaunay triangulation grids in the region, wherein the grids and nodes of the stratum interface offset in the step are not in one-to-one correspondence in the vertical direction, so that interpolation is needed to be carried out on the stratum interface at corresponding positions according to the plane coordinates of the newly established grid nodes, a new stratum interface is established in the process, and the new stratum interface is provided with uniform control point grids and can be aligned in the vertical direction, as shown in fig. 2. The new stratum interface still does not have the stratum topological structure, and the next step is to determine the stratum interface on the basis of the new stratum interface so that the geometric boundaries of the upper stratum and the lower stratum form a joint relation.
4) Extracting stratum interfaces to establish stratum, namely interpolating the stratum interfaces according to the sequence from top to bottom for each grid node coordinate arranged on a horizontal plane, wherein the bottom plate of an upper stratum and the top plate of a lower stratum in all vertically adjacent two stratum adopt common nodes to form a final stratum interface, and finally connecting the interfaces belonging to the same stratum to express each stratum through a triangular prism set;
5) The method comprises the steps of adding fault units, generating triangular grids with uniform surface grids inwards according to geometrical information of faults, giving distances between layers and fault interaction information to the surface grid units according to survey data after the generation is finished, wherein the fault interaction information comprises slip directions and slip distance information on two sides of the faults, the information represents the fault to belong to the categories of forward faults, reverse faults or walk-slip faults, grid data comprise node information and triangular unit information, each triangular grid is formed by connecting three nodes, the data of the faults are stored in two storage formats, one of the data is stored on each node, namely, the given attribute is unique to a certain node, the geometrical parameter is given to the node, the other data is stored on the triangular grid, namely, the given attribute is unique to a certain triangular grid, and the physical property parameter is given to the grid.
In the embodiment of the disclosure, as shown in fig. 2, the upper and lower interfaces of two strata are obtained by shifting the middle plane of the strata according to GIS data, and fig. 2 (a) shows a schematic diagram of the middle plane shift boundary of the strata. An area capable of containing the horizontal projection range of 4 formation interfaces is established on the horizontal plane and a uniform Delaunay triangulation grid is arranged, and fig. 2 (b) is a schematic diagram of Delaunay triangulation grid arrangement. The method comprises the steps of carrying out straight lines of vertical horizontal planes on each point in a grid, searching a stratum interface which is passed by the straight line from top to bottom on each straight line, calculating an intersection point of each stratum interface which is passed by the straight line, and taking the intersection point as a data point on the straight line, wherein the data point stores the top and bottom interface information of a stratum to which the straight line belongs, and (c) in fig. 2 is a stratum interface reconstruction illustration, and since the stratum interfaces are made in pairs of offset, a new interface point is calculated by taking an average value of the second point (the bottom plate of the first stratum) on the straight line and the elevation of the subsequent point (the top plate of the second stratum) on the straight line, wherein the point is a boundary point of the first stratum and the second stratum of the stratum, so that point distribution on the interfaces of all strata can be obtained by the method, and (d) in fig. 2 is an interface interpolation point illustration. And (e) in FIG. 2, connecting the interpolated points to form a stratum interface. According to the boundary range of the stratum interface, the boundaries around the stratum are connected up and down to form the boundary of the stratum side, the step establishes the complete stratum interface of each stratum and forms a topological structure among stratum geometries, and (f) in fig. 2 is a schematic of formation of the stratum structure. Finally, the fault geometry is added and fault physical parameters are assigned to the fault grid cells in accordance with the survey data, with (g) in fig. 2 being an illustration of the addition of fault data.
In the embodiment of the disclosure, the determination of the curved surface of the fault comprises the following types of data, namely a plane formed by three cross section point coordinates, known fault occurrence information and a plane passing through a certain point, and a curved surface formed by fault lines and fault inclination angle data.
In the embodiment of the disclosure, the surface grid unit and the volume unit of the whole solving domain are mutually independent, and the data stored in the surface grid comprises geometric parameters such as thickness, normal direction, physical parameters such as Young modulus, poisson ratio, cohesive force, internal friction angle and permeability coefficient.
According to the GIS modeling data improved by the fault geometric modeling method for the intelligent mine, not only is the space position information of the fault contained, but also the attribute information of the fault at different space positions is contained, and as the fault geometric modeling method is not related, a complex geometric calculation process is omitted, the fault information is mainly embodied on grid division and grid attribute configuration, and grids near the fault can be encrypted according to simulation requirements in specific application. When a volume grid is used for representing faults, smaller rigidity and strength parameters of fault units can be given to represent fault properties according to the width of the faults. Meanwhile, in the fluid seepage simulation, different permeability can be given according to water conductivity at different positions of faults, so that simulation analysis of a complex geological model of a mine is solved.
Example two
The embodiment discloses a fault geometric modeling system for an intelligent mine, which comprises an acquisition module, a modeling module, a stratum construction module and a fault data adding module;
The acquisition module is used for surveying and acquiring fault position data and physical and mechanical attribute data;
The modeling module is used for extracting middle plane thickness data in GIS data, offsetting the middle plane nodes, wherein the offset is half of the thickness to form an upper boundary and a lower boundary;
The stratum construction module is used for establishing an area on a horizontal plane, enabling projections of all stratum in GIS data on the horizontal plane to fall inside the area, arranging uniform Delaunay triangulation grids in the area, interpolating stratum interfaces according to the sequence from top to bottom for each grid node coordinate arranged on the horizontal plane, forming a final stratum interface by adopting common nodes on a bottom plate of an upper stratum and a top plate of a lower stratum in all vertically adjacent two strata, and finally connecting the interfaces belonging to the same stratum and expressing each stratum through a triangular prism set;
The fault data adding module is used for generating triangular grids with uniform surface grids inwards according to geometric information of faults, endowing distances among fault layers and fault interaction information on surface grid units according to survey data after the generation is finished, wherein the fault interaction information comprises slip directions and dislocation distance information on two sides of the faults, the information characterizes the faults to belong to categories of forward faults, reverse faults or walk-slip faults, grid data comprise node information and triangle unit information, each triangular grid is formed by connecting three nodes, the data of the faults are stored in two storage formats, one storage format is stored on each node, namely the endowed attribute is unique to a certain node, the geometric parameter is endowed on the node, the other storage format is stored on the triangular grid, namely the endowed attribute is unique to a certain triangular grid, and the physical parameter is endowed on the grid.
In the embodiment of the disclosure, the determination of the curved surface of the fault comprises the following types of data, namely a plane formed by three cross section point coordinates, known fault occurrence information and a plane passing through a certain point, and a curved surface formed by fault lines and fault inclination angle data.
In the embodiment of the disclosure, the surface grid unit and the volume unit of the whole solving domain are mutually independent, and the data stored in the surface grid comprises, but is not limited to, geometric parameters such as thickness, normal direction, physical parameters such as Young modulus, poisson ratio, cohesion, internal friction angle and permeability coefficient.
The process implemented by the fault geometry modeling system for intelligent mine in this embodiment is the same as the fault geometry modeling method for intelligent mine in the foregoing embodiment, and will not be described in detail here.
Example III
The embodiment discloses a fault geometry modeling device for an intelligent mine, which comprises a memory and a processor, wherein the memory is used for storing a computer program, and the processor is used for realizing the fault geometry modeling method for the intelligent mine according to the first embodiment when executing the computer program, and the specific modeling method steps are the same as those of the first embodiment and are not repeated.
Example IV
The embodiment discloses a computer readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the geometrical modeling method for fault of intelligent mine according to the foregoing embodiment is implemented, and the specific steps of the modeling method are the same as those of the foregoing embodiment, and are not repeated herein.
The computer of the embodiments of the present application may be a general purpose computer, a special purpose computer, a computer network, or other programmable apparatus. The computer instructions may be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. The computer readable storage medium may be any available medium that can be read by a computer or a data storage device such as a server, data center, etc. that contains an integration of one or more available media. The usable medium may be a magnetic medium (e.g., floppy disk, hard disk, magnetic tape), an optical medium (e.g., digital Versatile Disk (DVD)), or a semiconductor medium (e.g., solid State Disk (SSD)), etc. The software formed by the computer storage code can be located in random access memory, flash memory, read-only memory, programmable read-only memory or electrically erasable programmable memory, registers and other storage media which are mature in the field.
The functional modules in the embodiments of the present application may be integrated into one processing unit or module, or each module may exist alone physically, or two or more modules may be integrated into one unit or module. In the above embodiments, it may be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented in software, may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the processes or functions described in the embodiments of the present application are fully or partially implemented.
The foregoing is merely illustrative of the present application, and the present application is not limited thereto, and any person skilled in the art will readily recognize that variations or substitutions are within the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.