Coal mine digital twinning dynamic modeling method and system based on multi-source heterogeneous data fusion

Abstract

The invention discloses a coal mine digital twinning dynamic modeling method and system based on multi-source heterogeneous data fusion, and relates to the technical field of model construction.The method comprises the following steps that an independent loop is constructed, a wind track topology is constructed according to a preset first constraint condition and the independent loop, and a coal mine three-dimensional model is obtained based on a preset method; according to preset ventilation parameters and wind track topology, virtual particle boundaries are calculated, and the coal mine three-dimensional model is updated. The design efficiency is improved through the preset construction method and constraint conditions, the flow path and distribution condition of the wind flow in the mine can be more accurately simulated by introducing the virtual particle boundary based on the preset method and the mine correlation graph, ventilation dead angles or wind speed abnormal areas can be identified, and the method is suitable for popularization and application. The coal mine three-dimensional model is dynamically updated according to the virtual particle boundary, the reliability of simulation analysis is improved, and a scientific basis can be provided for making an emergency plan through the three-dimensional model and wind flow simulation.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of model construction, in particular to a coal mine digital twin dynamic modeling method and system based on multi-source heterogeneous data fusion. BACKGROUND

[0002] In recent years, digital twin technology has developed rapidly, and digital twin has deepened its integration with various industries. Digital twin technology is applied in urban management (smart transportation, zero-carbon park, urban emergency), smart industry (process industry, discrete industry), autonomous driving testing, and medical scenarios (smart hospital, precision medicine, drug research and development). Through digital twin technology, various industry application scenarios are empowered.

[0003] At present, in the Chinese invention patent with the publication number CN119918148A, a coal mine tunnel parameterized drawing modeling and graph-model linkage method is disclosed. This method draws a two-dimensional tunnel based on two-dimensional attribute information, and constructs a three-dimensional tunnel model based on three-dimensional attribute information, making it possible to express the tunnel model in three dimensions. This method solves the problems of manual drawing of tunnel drawings, the need for manual intervention in converting two-dimensional drawings to three-dimensional models, high labor costs, low efficiency, and lack of support for collaboration. At the same time, the graph-model linkage method can realize bidirectional linkage update of two-dimensional tunnels and three-dimensional tunnel models. However, related technologies do not compensate for the main three-dimensional space from the ventilation perspective, which is not conducive to the rapidity of the model part of the commonly used and frequently changing space, and does not compensate for the dead angle area based on the non-dead angle area, which is not conducive to the accuracy of model construction. SUMMARY

[0004] The technical problem solved by the present application is that related technologies do not compensate for the main three-dimensional space from the ventilation perspective, which is not conducive to the rapidity of the model part of the commonly used and frequently changing space, and does not compensate for the dead angle area based on the non-dead angle area, which is not conducive to the accuracy of model construction.

[0005] To solve the above technical problems, the present application provides the following technical solutions. In a first aspect, a coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion includes the following steps: step S100, constructing a coordinate system, constructing an independent loop according to a predetermined construction method, and constructing a wind track topology according to a predetermined first constraint condition and the independent loop; Step S200, constructing a three-dimensional grid according to mine-related graphs, and replacing related parts of the three-dimensional grid through GTP algorithm and TIN algorithm to obtain a coal mine three-dimensional model; Step S300, setting fluid particle related parameters according to the wind track topology, setting virtual particle related parameters according to the fluid particle related parameters, and calculating the fluid particle related data and the virtual particle related data according to the Taylor expansion formula to obtain a virtual particle boundary. In step S400, a transformation relationship of coordinates is constructed according to the virtual particle boundary, coordinates of points of a corresponding part of the coal mine three-dimensional model are obtained according to the transformation relationship of coordinates, and the original coordinates in the coal mine three-dimensional model are replaced according to the coordinates of the points.

[0006] As a preferred scheme of the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion provided in the application, step S100 comprises the following sub-steps, including step S101, obtaining mine branch information; The mine branch information comprises a branch name, a branch number and a branch geometric center point position, and the branch geometric center point position is represented by a branch geometric center point longitude and a branch geometric center point latitude. In step S102, the branch geometric center point position is converted into a coordinate point in a three-dimensional Cartesian coordinate system, which is denoted as a first coordinate point. In step S103, an unordered topology is constructed according to the first coordinate point and the corresponding branch number. In step S104, simulation parameters and first constraint conditions are set, an independent loop is constructed according to a preset construction method, and a wind track topology is constructed according to the independent loop, an air inlet parameter, an air outlet parameter and the first constraint condition, and the wind track topology is an ordered topology.

[0007] As a preferred scheme of the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion provided in the application, an unordered sub-topology graph is constructed according to the air inlet, the air outlet and the return airway in the unordered topology graph, a distance minimum method is used to screen a path node, whether the path node is updated or jumps to a next unordered sub-topology graph is judged based on the return airway and the unordered sub-topology graph, and whether the path node is further updated or jumps to a next unordered sub-topology graph is judged according to the air outlet and the unordered sub-topology graph.

[0008] As a preferred scheme of the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion provided in the application, the construction method of the independent loop comprises the following steps, that is, obtaining a reserved ordered path graph, taking a reference starting node of any reserved ordered path graph as a starting point, taking a reference end node of the reserved ordered path graph as a midpoint, taking a path node of other reserved ordered path graphs except the reference starting node and the reference end node as a backup point, and constructing an independent loop according to the main point, the backup point and the terminal point, wherein the independent loop has the same meaning as the independent loop in circuit analysis, that is, the independent loop does not contain other loops, and the independent loop is configured with an acquisition constraint condition, the acquisition constraint condition is represented by that the number of nodes contained in the independent loop is the largest, the independent loop is numbered according to the branch number corresponding to the reference starting node, and the numbering of the independent loop is denoted as .

[0009] As the preferred scheme of the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion, after the independent loop is constructed, simulation parameters are set, the simulation parameters include air inlet parameters and air outlet parameters, the air inlet parameters include air inlet pressure and air inlet time point, the air pressure at each node in the independent loop is monitored, the air pressure at the air outlet is monitored, when the increase of the value of the air pressure at the air outlet is a first pressure value, the time point is obtained and recorded as an air outlet time point, the maximum value of the monitored air pressure at each node is obtained within the air inlet time point and the air outlet time point, and the maximum value of the monitored air pressure at each node is set as the air pressure at the corresponding node. The first constraint condition is that the air pressure reduction speed in the independent loop is minimum.

[0010] As the preferred scheme of the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion, the air pressure reduction speed of each independent loop is traversed, the independent loop with the minimum air pressure reduction speed is selected, and the independent loop with the minimum air pressure reduction speed is set as the air track topology, wherein the direction of the air track topology corresponds to the branch direction. Step S200 includes the following sub-steps, step S201, obtaining a mine related graph, and generating a three-dimensional grid; Step S202, according to the GTP algorithm, the stratum is constructed, according to the TIN algorithm, the fault is constructed, and according to the stratum and the fault, the part corresponding to the three-dimensional grid is replaced; Step S203, the three-dimensional geological model after replacement is recorded as a coal mine three-dimensional model.

[0011] As the preferred scheme of the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion, the mine related graph is obtained, each information in the mine related graph is converted into a coordinate point according to a three-dimensional Cartesian coordinate system, the coordinate points are supplemented according to the DIS interpolation algorithm, and the three-dimensional grid is obtained according to the supplemented coordinate points.

[0012] As the preferred scheme of the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion, according to the GTP algorithm, the stratum is constructed, and for the stratum, the construction method is configured, including taking the coordinate grid point as the reference point, generating the TIN model through the point-by-point interpolation method, obtaining the sequence of the stratum, obtaining the top elevation of each layer, setting the top elevation of each layer as the vertex coordinate of the corresponding triangle, setting the bottom elevation of each layer as the bottom coordinate of the corresponding triangle, setting the lithology as the internal attribute of the corresponding triangle, modifying the part corresponding to the three-dimensional grid according to each triangle with the internal attribute, and setting the part corresponding to the modified three-dimensional grid as the stratum. The part corresponding to the fault in the three-dimensional grid is acquired, a structural surface of the fault is constructed according to a TIN algorithm, and the three-dimensional grid is updated according to the structural surface of the fault.

[0013] As a preferred scheme of the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion, the virtual particle boundary is calculated according to a preset formula, the virtual particle boundary represents a virtual particle velocity of the wind track topology on the outermost boundary of the coal mine three-dimensional model, and the virtual particle velocity is obtained through a Taylor expansion formula. The speed point cloud of the virtual boundary is obtained according to the speed of the virtual particle b, denoted as the virtual particle boundary, and the transformation relationship between the speed of the virtual particle b and the coordinates of the corresponding point is obtained according to the virtual particle boundary and the virtual particle boundary corresponding to the first space at the entrance of the coal mine three-dimensional model, and the transformation relationship is represented as, ; Wherein, F1~F3 respectively represent the transformation rules of the virtual particle speed and the corresponding outermost boundary coordinate points x-axis coordinate, y-axis coordinate and z-axis coordinate. According to the new three-dimensional space region and the coal mine three-dimensional model, the coal mine three-dimensional model is updated, and the update is represented as replacing the original three-dimensional space region in the coal mine three-dimensional model with the new three-dimensional space region.

[0014] In the second aspect, the coal mine digital twin dynamic modeling system based on multi-source heterogeneous data fusion includes a construction module, a calculation module and an update module. The construction module is used to construct a coordinate system, construct an independent loop according to a preset construction method, construct a wind track topology according to a preset first constraint condition and the independent loop, construct a three-dimensional grid according to a mine related graph, and replace the related parts of the three-dimensional grid through a GTP algorithm and a TIN algorithm to obtain a coal mine three-dimensional model. The calculation module calculates the virtual particle boundary according to a preset ventilation parameter and the wind track topology. The update module constructs a coordinate transformation relationship according to the virtual particle boundary, obtains the coordinates of the points of the corresponding part of the coal mine three-dimensional model according to the coordinate transformation relationship, and replaces the original coordinates in the coal mine three-dimensional model according to the coordinates of the points.

[0015] The beneficial effects of the present application: through the preset construction method and constraint conditions, the wind track topology can be systematically constructed, the randomness of manual design is reduced, the design efficiency is improved, based on the preset method and mine related graph, the coal mine three-dimensional model can be quickly generated, the time of traditional manual modeling is saved, through the introduction of virtual particle boundary, the flow path and distribution of air flow in the mine can be more accurately simulated, which helps to identify the ventilation dead angle or abnormal wind speed area, the coal mine three-dimensional model is dynamically updated according to the virtual particle boundary, the model is closer to the actual ventilation condition, and the credibility of simulation analysis is improved, through the three-dimensional model and air flow simulation, the air flow change under the sudden events such as fire and gas explosion can be simulated, and scientific basis is provided for the formulation of emergency plan. BRIEF DESCRIPTION OF DRAWINGS

[0016] Figure 1 The basic flow diagram of the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion provided by an embodiment of the present application is shown. DETAILED DESCRIPTION

[0017] The technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

[0018] It should be understood that the step numbers used herein are only for the convenience of description, and are not limited to the execution sequence of the steps. It should be understood that the terms used in the specification of the present application are only for the purpose of describing specific embodiments and are not intended to limit the present application.

[0019] As used in the specification and the appended claims of the present application, the singular forms "a," "an," and "the" are intended to include plural forms unless the context clearly dictates otherwise.

[0020] The terms "comprise" and "include" indicate the presence of described features, whole, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, whole, steps, operations, elements, components and / or sets thereof

[0021] The term "and / or" means any combination of one or more of the associated listed items and all possible combinations, and includes these combinations.

[0022] With the development of digital twin technology, the coal mine industry can be empowered. Data-driven technology is the core, and three-dimensional geological automatic modeling, multi-dimensional data fusion analysis, and digital twin platform construction are the technical means. A comprehensive perception and virtual-real interaction coal mine digital twin system is constructed to provide a first-hand condition for intelligent decision-making and precise control. High-precision automatic modeling and dynamic updating of underground coal seams, faults, and roadways are achieved. Real-time visualization display and analysis of mine geological data are supported. Precise twin models can be used to further study disaster warning, ventilation simulation, production scheduling, etc. This significantly improves the efficiency and safety level of coal mine management and helps improve the intelligent level of auxiliary decision-making for coal mine production scheduling.

[0023] Based on this, the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion is provided.

[0024] The coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion provided in the embodiments of the present application is described in detail as follows. First, the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion in the embodiments of the present application is described.

[0025] The coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion provided in the embodiments of the present application is related to model construction. The coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion provided in the embodiments of the present application can be applied in a terminal, can be applied in a server end, and can also be software running in a terminal or a server end. In some embodiments, the terminal can be a smart phone, a tablet computer, a notebook computer, a desktop computer, etc. The server end can be configured as an independent physical server, can be configured as a server cluster or a distributed system composed of multiple physical servers, can also be configured as a cloud server providing cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, CDN, and basic cloud computing services such as big data and artificial intelligence platforms, etc. The software can be an application for implementing the coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion, etc., but is not limited to the above forms.

[0026] The application can also be implemented in a variety of general purpose or special purpose computer systems environments or configurations. Examples of well known computer systems, environments, or configurations that can be suitable for use with the application include, personal computers, server computers, handheld or laptop devices, tablet devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like. The application can be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, and the like, that perform particular tasks or implement particular abstract data types. The application can also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules can be located in local and remote computer storage media including memory storage devices.

[0027] Embodiments, with reference to Figure 1 For an embodiment of the application, a coal mine digital twin dynamic modeling method based on multi-source heterogeneous data fusion is provided, including the following steps: step S100, constructing a coordinate system, constructing an independent loop according to a preset construction method, and constructing a wind track topology according to a preset first constraint condition and the independent loop; Step S200, constructing a three-dimensional grid according to a mine-related graph, and replacing relevant parts of the three-dimensional grid through GTP algorithm and TIN algorithm to obtain a coal mine three-dimensional model; Step S300, setting fluid particle-related parameters according to the wind track topology, setting virtual particle-related parameters according to the fluid particle-related parameters, and calculating the fluid particle-related data and the virtual particle-related data according to the Taylor expansion formula to obtain a virtual particle boundary; Step S400, constructing a coordinate transformation relationship according to the virtual particle boundary, obtaining the coordinates of a point corresponding to a part of the coal mine three-dimensional model according to the coordinate transformation relationship, and replacing original coordinates in the coal mine three-dimensional model according to the coordinates of the point.

[0028] Further preferably, the wind track topology can be systematically constructed through the preset construction method and constraint condition, reducing the randomness of manual design and improving design efficiency. Based on the preset method and mine-related graph, the coal mine three-dimensional model can be quickly generated, saving the time of traditional manual modeling. By introducing the virtual particle boundary, the flow path and distribution of the air flow in the mine can be more accurately simulated, which helps to identify the ventilation dead angle or abnormal wind speed area. According to the virtual particle boundary, the coal mine three-dimensional model is dynamically updated, making the model closer to the actual ventilation condition and improving the credibility of simulation analysis. Through the three-dimensional model and air flow simulation, the air flow changes under sudden events such as fire and gas explosion can be simulated, providing a scientific basis for the development of emergency plans.

[0029] The step S100 comprises the following sub-steps, including, step S101, acquiring mine branch information; The mine branch information comprises branch name, branch number, branch geometric center point position, the branch geometric center point position is expressed as branch geometric center point longitude, branch geometric center point latitude; Step S102, convert the branch geometric center point position into a coordinate point in a three-dimensional Cartesian coordinate system, denoted as a first coordinate point; Step S103, according to the first coordinate point, the corresponding branch number, construct an unordered topology; Step S104, set simulation parameters and first constraint conditions, according to the preset construction method, construct independent loops, according to the independent loops, the air inlet parameter, the air outlet parameter, the first constraint condition, build a wind track topology, the wind track topology is an ordered topology.

[0030] Further preferably, the branch name comprises main inclined shaft, auxiliary inclined shaft, pedestrian inclined shaft, track dark inclined shaft, track dark inclined shaft upper yard, track dark inclined shaft lower yard, track transfer lane, belt dark inclined shaft, two mining area belt down slope, two mining area track down slope, transportation crossheading, return air crossheading, ventilation opening, return air roadway, air outlet.

[0031] Further preferably, the branch number corresponds to the branch name, and the branch number is expressed as Q ij Wherein, i represents the i-th branch name, j represents the j-th branch site of the i-th branch name, i and j are natural numbers.

[0032] Further preferably, for the first coordinate point, based on the formula for converting latitude and longitude to a three-dimensional Cartesian coordinate system, the formula for converting latitude and longitude to a three-dimensional Cartesian coordinate system comprises, Wherein, x, y, z are the horizontal axis coordinate, vertical axis coordinate, vertical axis coordinate in the three-dimensional Cartesian coordinate system, is the latitude, is the longitude, is the curvature radius of the earth, is the altitude, is the eccentricity of the ellipsoid, that is, is a constant.

[0033] ​Further preferably, for the disordered topology graph, a construction strategy is configured, including, taking the first coordinate point of the air outlet, the first coordinate point of the air return lane, and the first coordinate point corresponding to the air outlet as reference points, connecting the first coordinate point of the air outlet with each of the first coordinate points other than the first coordinate point corresponding to the air outlet as the reference point, connecting the first coordinate point corresponding to the air outlet with each of the first coordinate points other than the first coordinate point of the air outlet as the reference point, and connecting the first coordinate point of the air return lane with each of the first coordinate points other than the first coordinate point of the air outlet and the first coordinate point corresponding to the air outlet as the reference point, to obtain a three-dimensional topology graph, which is denoted as the disordered topology graph.

[0034] According to the air outlet, the air outlet, and the air return lane in the disordered topology graph, a disordered sub-topology graph is constructed, a route node is screened according to the distance minimum method, and whether the route node is updated or jumps to the next disordered sub-topology graph is determined based on the air return lane and the disordered sub-topology graph. According to the air outlet and the disordered sub-topology graph, it is determined whether the route node is further updated or jumps to the next disordered sub-topology graph.

[0035] Further preferably, step S104 includes the following sub-steps. In step S1041, any air outlet is taken as a reference starting node, the air outlet with the minimum straight-line distance from the reference starting node is taken as a reference end node, the air return lane with the minimum straight-line distance from the reference starting node is taken as a reference relay node, and the reference starting node, the reference relay node, and the reference end node are constructed into a disordered sub-topology graph according to the construction logic of the disordered topology graph. In each disordered sub-topology graph, the air outlets are different, and the reference end nodes of different disordered sub-graphs can be the same. In step S1042, any disordered sub-topology graph is jumped to for analysis, the node with the closest straight-line distance from the reference starting node is obtained, the node is set as a route node, the reference starting node and the route node are connected in a straight line, and an arrow in the direction from the reference starting node to the route node is marked on the connected straight line. In step S1043, the route node is taken as a new reference starting node, and steps S1042-S1043 are repeated until the new reference starting node is the reference end node, the repetition of steps S1042-S1043 is terminated, an ordered path graph is obtained, and step S1044 is jumped to. In step S1044, it is determined whether the ordered path graph includes the reference relay node. When the ordered path graph includes the reference relay node, step S1045 is jumped to, and when the ordered path graph does not include the reference relay node, step S1046 is jumped to. In step S1045, it is determined whether there is an obtuse angle between each straight line from the relay node to the reference end node and the straight line corresponding to the reference starting node. When the angle between each straight line from the relay node to the reference end node and the corresponding straight line of the reference start node is an obtuse angle, jump to step S1047; When the angle between each straight line from the relay node to the reference end node and the corresponding straight line of the reference start node is not an obtuse angle, output the ordered path graph; In step S1046, the minimum distance in the screening method of any new reference start node in step S1043 is updated to the second minimum distance, or the minimum distance in the screening method of the passing node in step S1042 is updated to the second minimum distance, and steps S1042 to S1044 are cycled; When the minimum distance in the screening method of all new reference start nodes in step S1043 is updated to the second minimum distance, and the minimum distance in the screening method of the passing node in step S1042 is updated to the second minimum distance, and there is no reference relay node, delete each node in the unordered sub-topology graph, and jump to the next unordered sub-topology graph; When there is a reference relay node, stop the step of updating the minimum distance in the screening method of any new reference start node in step S1043 to the second minimum distance, or updating the minimum distance in the screening method of the passing node in step S1042 to the second minimum distance, and determine whether the angle between each straight line from the relay node to the reference end node and the corresponding straight line of the reference start node is an obtuse angle, When the angle between each straight line from the relay node to the reference end node and the corresponding straight line of the reference start node is an obtuse angle, continue the step of updating the minimum distance in the screening method of any new reference start node in step S1043 to the second minimum distance, or updating the minimum distance in the screening method of the passing node in step S1042 to the second minimum distance; When the angle between each straight line from the relay node to the reference end node and the corresponding straight line of the reference start node is not an obtuse angle, output the ordered path graph; When the minimum distance in the screening method of all new reference start nodes in step S1043 is updated to the second minimum distance, and the minimum distance in the screening method of the passing node in step S1042 is updated to the second minimum distance, and the angle between each straight line from the relay node to the reference end node and the corresponding straight line of the reference start node is an obtuse angle, delete each node in the unordered sub-topology graph, and jump to the next unordered sub-topology graph.

[0036] The construction method of the independent loop includes: obtaining the reserved ordered path diagram, taking the reference starting node of any reserved ordered path diagram as a starting point, taking the reference end node of the reserved ordered path diagram as a midpoint, and taking the path nodes of the other reserved ordered path diagrams except the reference starting node and the reference end node as backup points, and constructing an independent loop according to the main point, the backup point and the terminal point, wherein the independent loop has the same meaning as the independent loop in circuit analysis, and is represented as not containing other loops in the independent loop, and the independent loop is configured with an acquisition constraint condition, and the acquisition constraint condition is represented as that the number of nodes contained in the independent loop is the most, the independent loop is numbered according to the branch number corresponding to the reference starting node, and the number of the independent loop is represented as .

[0037] After the independent loop is constructed, simulation parameters are set, the simulation parameters include air inlet parameters and air outlet parameters, the air inlet parameters include air inlet pressure and air inlet time point, the air pressure at each node in the independent loop is monitored, the air pressure at the air outlet is monitored, when the increase of the value of the air pressure at the air outlet is a first pressure value, the time point is obtained and recorded as an air outlet time point, and the maximum value of the monitored air pressure at each node is obtained within the air inlet time point and the air outlet time point, and the maximum value of the monitored air pressure at each node is set as the air pressure at the corresponding node.

[0038] The first constraint condition is represented as that the air pressure reduction speed in the independent loop is the smallest.

[0039] Further preferably, a calculation method is configured for the air pressure reduction speed, and the calculation method of the air pressure reduction speed includes: splitting the independent loop into two branches according to the direction of the starting point and the terminal point, selecting any two adjacent nodes of the branch, the adjacent nodes being represented as the closest straight-line distance, recording as a node group, obtaining the air pressure of the node group, calculating the first difference value of the air pressure of the node group, traversing the first difference value of each node group of the branch, calculating the average value of the first difference value, obtaining the average value of the first difference value of the other branch, calculating the average value of the average value of the first difference value, and recording the average value of the average value of the first difference value as the average pressure drop of the independent loop; The difference value between the air outlet time point and the air inlet time point is calculated and recorded as the pressure drop time of the independent loop; The first ratio of the average pressure drop and the pressure drop time is calculated, and the first ratio is set as the air pressure reduction speed.

[0040] The air pressure reduction speeds corresponding to each independent loop are traversed, the independent loop with the smallest air pressure reduction speed is selected, and the independent loop with the smallest air pressure reduction speed is set as a wind track topology, wherein the direction of the wind track topology corresponds to the branch direction; Step S200 includes the following sub-steps: step S201, obtaining a mine related diagram to generate a three-dimensional grid; Step S202, according to the GTP algorithm to construct strata, according to the TIN algorithm to construct faults, according to the strata, faults, replace the corresponding part of the three-dimensional grid; Step S203, the replaced three-dimensional geological model is recorded as a coal mine three-dimensional model.

[0041] Further preferably, the mine-related map includes a mine structure outline map, a mine well-to-well profile contrast line design map, a mine stratum comprehensive columnar chart, a mine fault development feature, a sedimentary cycle profile map; In the mine structure outline map, the faults are represented by red solid lines, the anticlines are represented by green solid lines, the synclines are represented by green dashed lines, and the denudation area is represented by a blue closed curve; In the mine well-to-well profile contrast line design map, the borehole coordinates and the coal seam thickness are included; In the mine stratum comprehensive columnar chart, the stratum era system, the cross-sectional columnar chart, the rock thickness, and the lithology are included; The mine fault development feature includes a fault dip rose diagram, a fault dip distribution histogram, and a fault throw distribution histogram; The sedimentary cycle profile map includes a cycle profile map, a sedimentary facies, a sedimentary microfacies, and a sedimentary subfacies, wherein the sedimentary microfacies and the sedimentary subfacies do not appear at the same time, the cycle profile map includes a positive cycle profile map, a reverse cycle profile map, a symmetric cycle profile map, and a uniform cycle profile map, when the cycle profile map is a positive cycle profile map, or a reverse cycle profile map, or a uniform cycle profile map, the sedimentary microfacies exists and the sedimentary subfacies does not exist, and when the cycle profile map is a symmetric cycle profile map, the sedimentary subfacies exists and the sedimentary microfacies does not exist.

[0042] Further preferably, the stratum era system includes the Cambrian system, the Ordovician system, the Silurian system, the Devonian system, the Carboniferous system, the Permian system, the Triassic system, the Jurassic system, the Cretaceous system, the Paleogene system, the Neogene system, and the Quaternary system; The rock thickness and the coal seam thickness are represented as the average value of the vertical thickness of the rock in the rock area and the average value of the vertical thickness of the coal seam in the coal seam area, and are calculated by the calculation formula of the mean value theorem of integration; The borehole coordinates are represented as the central longitude and the central latitude of the circular face of the set borehole; The lithology includes magmatic rock (igneous rock), sedimentary rock, and metamorphic rock, wherein the magmatic rock (igneous rock) includes granite, granite porphyry, rhyolite, syenite, basalt, and trachyte, the sedimentary rock includes arenite, quartz gravel, quartz angular gravel, lime gravel, and marl, and the metamorphic rock includes gneiss, mica schist, chlorite schist, marble, and hornblende schist; The sedimentary facies includes a terrestrial group, a transitional group, and a marine group, wherein the terrestrial group includes residual facies, slope facies, and aeolian facies, the transitional group includes delta facies and estuary facies, and the marine group includes shore facies and shallow sea shelf facies.

[0043] Obtain a mine-related graph, convert each information in the mine-related graph into a coordinate point according to a three-dimensional Cartesian coordinate system, supplement the coordinate points according to a DIS interpolation algorithm, and obtain a three-dimensional grid according to the supplemented coordinate points.

[0044] Further preferably, obtain a boundary line according to a mine structure outline graph, the boundary line includes each solid line and dashed line, obtain an x-axis coordinate distribution range and a y-axis coordinate distribution range according to the boundary line, and the coordinate points are configured with a supplementary constraint condition, the supplementary constraint condition includes that a straight-line distance between any two coordinate points is less than or equal to a first distance, the coordinate points are interpolated along a direction of a unit vector of a z-axis direction, and the coordinate points are interpolated within a coordinate distribution range corresponding to the boundary line.

[0045] Further preferably, the x-axis coordinate distribution range and the y-axis coordinate distribution range are both double-closed intervals, and the coordinate points interpolated within the coordinate distribution range corresponding to the boundary line are expressed as that an x-coordinate of the interpolated coordinate point is distributed in the x-axis coordinate distribution range, a y-coordinate of the interpolated coordinate point is distributed in the y-axis coordinate distribution range, and a z-coordinate of the interpolated coordinate point is unrestricted.

[0046] Construct a stratum according to a GTP algorithm, and the stratum is configured with a construction method including taking a coordinate grid point as a reference point, generating a TIN model through a point-by-point interpolation method, obtaining a sequence of the stratum, obtaining a top elevation of each layer, setting the top elevation of each layer as a vertex coordinate of a corresponding triangle, setting a bottom elevation of each layer as a bottom coordinate of the corresponding triangle, setting a lithology as an internal attribute of the corresponding triangle, modifying a corresponding part of the three-dimensional grid according to each triangle with the internal attribute set, and setting the modified corresponding part of the three-dimensional grid as the stratum.

[0047] Obtain a part corresponding to a fault in the three-dimensional grid, construct a structural surface of the fault according to a TIN algorithm, and update the three-dimensional grid according to the structural surface of the fault.

[0048] Further preferably, the fault includes a high-angle normal fault, a low-angle normal fault, a high-angle reverse fault, and a low-angle reverse fault, obtain fault parameters including a fault throw and a fault dip angle, and the fault includes a footwall top surface and a hanging wall top surface.

[0049] Further preferably, obtain any two coordinate points on a fault line of the footwall top surface, denoted as A(x1, y1, z1) and B(x2, y2, z2), calculate a coordinate of a projection point of A on the hanging wall top surface , and the calculation expression of the coordinate of the projection point is ; wherein, is a fault dip angle, is a formation thickness, is a fault throw.

[0050] The intersection of the fault surface and the formation surface is obtained by a Boolean operation, and is recorded as a quasi-edge line. According to the calculation expression of the coordinates of the projection point , the upper edge line of the fault is updated so that any point in the upper edge line coincides with the projection of the corresponding point in the lower edge line. The updated quasi-edge line is recorded as the actual edge line. According to the actual edge line, the corresponding part in the three-dimensional grid is modified to generate a fault that conforms to the actual geological phenomenon. The corresponding part in the modified three-dimensional grid is recorded as the structural surface of the fault.

[0051] According to a preset formula, a virtual particle boundary is calculated, which represents the virtual particle velocity of the wind track topology on the outermost boundary of the coal mine three-dimensional model. The virtual particle velocity is obtained through a Taylor expansion formula.

[0052] Further preferably, for the outermost boundary of the wind track topology in the coal mine three-dimensional model, an acquisition method is configured to acquire a three-dimensional space region of the wind track topology on the coal mine three-dimensional model. The three-dimensional space region does not include coal mine mechanical equipment, personnel, rock, and coal seams. The curved surface near the inside of the three-dimensional space region is set as the outermost boundary.

[0053] Further preferably, one virtual particle b of the outermost boundary is acquired, and a fluid particle with a straight-line distance of a second distance from the virtual particle b is acquired. The wind pressure of the virtual particle b and the wind pressure of the fluid particle are acquired. The wind pressure of the fluid particle is monitored by a wind pressure sensor. The wind pressure of the virtual particle is Taylor expanded at the fluid particle m. The calculation expression of the Taylor expansion is ; ; ; ; wherein, represents the wind pressure estimate value of the fluid particle in the calculation domain, i.e., the wind pressure value detected at the node closest to the virtual particle b, represents the unit vector in the normal direction of the face element corresponding to the virtual particle b, represents the unit vector in the tangent direction of the face element corresponding to the virtual particle b, represents the projection vector of the line connecting the fluid particle coordinate point and the origin in the plane formed by the normal and tangent of the face element, represents the projection vector of the line connecting the virtual particle coordinate point and the origin in the plane formed by the normal and tangent of the face element, represents the modulus of the orientation vector, is expressed as a high-order infinitesimal, is expressed as the acceleration of gravity, is expressed as the wind energy density, is the wind speed of the fluid particle, that is, the wind pressure value detected at the node closest to the fluid particle b, is a constant, M is the Mth ventilation, U m is expressed as the air volume detected at the node closest to the fluid particle in the first time period, is the straight-line distance between the virtual particle b and the fluid particle, is the velocity compensation value of the fluid particle to the virtual particle.

[0054] Further preferably, according to the calculation expression of Taylor expansion, the velocity of the virtual particle b is obtained, and the calculation expression of the velocity of the virtual particle b is, ; wherein, is the velocity of the virtual particle b.

[0055] According to the velocity of the virtual particle b, the velocity point cloud of the virtual boundary is obtained, denoted as virtual particle boundary, according to the virtual particle boundary, according to the virtual particle boundary corresponding to the first space at the entrance of the coal mine three-dimensional model, the conversion relationship between the velocity of the virtual particle b and the coordinates of the corresponding point is obtained, and the conversion relationship is expressed as, ; wherein, F1~F3 respectively represent the conversion rule of the velocity of the virtual particle and the coordinate points x-axis coordinate, y-axis coordinate and z-axis coordinate of the corresponding outermost boundary.

[0056] Further preferably, according to the conversion relationship and the virtual particle boundary, the outermost boundary corresponding to the wind track topology is reconstructed to obtain a new three-dimensional space region.

[0057] According to the new three-dimensional space region and the coal mine three-dimensional model, the coal mine three-dimensional model is updated, and the update is expressed as replacing the original three-dimensional space region in the coal mine three-dimensional model with the new three-dimensional space region.

[0058] Further preferably, by the preset construction method and constraint condition, the wind track topology can be systematically constructed, the randomness of human design is reduced, the design efficiency is improved, based on the preset method and the mine related graph, the coal mine three-dimensional model can be quickly generated, the time of traditional manual modeling is saved, by introducing the virtual particle boundary, the flow path and distribution of the air flow in the mine can be more accurately simulated, which helps to identify the ventilation dead angle or the wind speed abnormal area, according to the virtual particle boundary, the coal mine three-dimensional model is dynamically updated, so that the model is closer to the actual ventilation condition, the credibility of the simulation analysis is improved, through the three-dimensional model and the air flow simulation, the air flow change under the sudden events such as fire and gas explosion can be simulated, which provides a scientific basis for the formulation of emergency plan.

[0059] Those skilled in the art will appreciate that embodiments of the present application can be readily used as a method, a system or a computer program product. Accordingly, the present application can take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application can take the form of a computer program product on one or more computer-usable storage media (or computer- readable storage media) having computer-usable program code embodied in the medium. The medium can be any available storage media that can be accessed by a computer. By way of example, and not limitation, such computer-usable storage media can include RAM, ROM, EEPROM, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other storage medium(s) that can be used to carry or store desired computer program code in the form of instructions or data structures and that can be accessed by a computer. Also, the present application can be embodied in a computer program product that can be traded as goods or merchandise, through the based on any such medium, transactional medium, or physical medium. Figure 1 one or more flows and / or blocks Figure 1 one or more flows and / or blocks

[0060] It should be noted that the above-mentioned embodiments are only used to illustrate but not to limit the technical solutions of the present application, and although the present application is described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present application can be modified or replaced equivalently without departing from the spirit and scope of the technical solutions of the present application, and they should be covered in the protection scope of the present application.

Claims

1. A dynamic modeling method for coal mine digital twins based on multi-source heterogeneous data fusion, characterized in that, The steps include: Step S100, constructing a coordinate system, constructing independent loops according to a preset construction method, and constructing a wind track topology according to a preset first constraint condition and independent loops. Step S200: Construct a three-dimensional mesh based on the relevant mine map, and replace the relevant parts of the three-dimensional mesh using the GTP algorithm and TIN algorithm to obtain a three-dimensional coal mine model; Step S300: Set fluid particle related parameters according to the wind track topology, set virtual particle related parameters according to the fluid particle related parameters, and calculate the fluid particle related data and virtual particle related data according to the Taylor expansion formula to obtain the virtual particle boundary. Step S400: Construct the coordinate transformation relationship based on the virtual particle boundary, obtain the coordinates of the corresponding points in the three-dimensional coal mine model based on the coordinate transformation relationship, and replace the original coordinates in the three-dimensional coal mine model based on the coordinates of the points.

2. The dynamic modeling method for coal mine digital twins based on multi-source heterogeneous data fusion as described in claim 1, characterized in that, Step S100 includes the following sub-steps, including step S101, obtaining mine branch information; Mine branch information includes branch name, branch number, and the location of the branch geometric center point. The location of the branch geometric center point is expressed as the longitude and latitude of the branch geometric center point. Step S102: Convert the position of the branch geometry center point into a coordinate point in a three-dimensional Cartesian coordinate system, and denot it as the first coordinate point; Step S103: Construct an unordered topology based on the first coordinate point and the corresponding branch number; Step S104: Set simulation parameters and first constraint conditions. Construct independent loops according to the preset construction method. Based on the independent loops, inlet parameters, outlet parameters, and first constraint conditions, construct the air duct topology, which is an ordered topology.

3. The dynamic modeling method for coal mine digital twins based on multi-source heterogeneous data fusion as described in claim 2, characterized in that, Based on the ventilation openings, air outlets, and return airways in the unordered topology graph, construct an unordered sub-topology graph. Filter path nodes using the minimum distance method. Based on the return airway and the unordered sub-topology graph, determine whether the path nodes should be updated or jump to the next unordered sub-topology graph. Based on the air outlet and the unordered sub-topology graph, determine whether the path nodes should be further updated or jump to the next unordered sub-topology graph.

4. The dynamic modeling method for coal mine digital twins based on multi-source heterogeneous data fusion as described in claim 1, characterized in that, The method for constructing independent loops includes: obtaining a preserved ordered path graph; taking the reference start node of any preserved ordered path graph as the starting point; taking the reference end node of any preserved ordered path graph as the midpoint; and taking the path nodes of other preserved ordered path graphs (excluding the reference start node and reference end node) as backup points; constructing an independent loop based on the main point, backup points, and end point. Here, "independent loop" has the same meaning as "independent loop" in circuit analysis, indicating that the independent loop does not contain other loops. For each independent loop, constraints are configured, which are expressed as: the independent loop contains the maximum number of nodes. The independent loops are then numbered according to the branch number corresponding to the reference start node. The numbering of the independent loop is represented as follows: .

5. The dynamic modeling method for coal mine digital twins based on multi-source heterogeneous data fusion as described in claim 4, characterized in that, After the independent loop is constructed, simulation parameters are set, including inlet parameters and outlet parameters. Inlet parameters include inlet pressure and inlet time point. The air pressure at each node in the independent loop is monitored, and the air pressure at the outlet is monitored. When the increase in the air pressure at the outlet is equal to the first pressure value, the time point is obtained and recorded as the outlet time point. Within the inlet and outlet time points, the maximum value of the air pressure at each node is obtained, and the maximum value of the air pressure at each node is set as the air pressure at the corresponding node. The first constraint is expressed as follows: the rate of decrease in wind pressure in an independent loop is minimized.

6. The dynamic modeling method for coal mine digital twins based on multi-source heterogeneous data fusion as described in claim 4, characterized in that, Iterate through the wind pressure reduction rate corresponding to each independent loop, select the independent loop with the smallest wind pressure reduction rate, and set the independent loop with the smallest wind pressure reduction rate as the wind rail topology, where the direction of the wind rail topology corresponds to the direction of the branch. Step S200 includes the following sub-steps: Step S201, obtain relevant mine maps and generate a three-dimensional mesh; Step S202: Construct the strata according to the GTP algorithm, construct the faults according to the TIN algorithm, and replace the corresponding parts of the three-dimensional mesh according to the strata and faults. Step S203: Record the replaced three-dimensional geological model as the coal mine three-dimensional model.

7. The dynamic modeling method for coal mine digital twins based on multi-source heterogeneous data fusion as described in claim 1, characterized in that, Obtain relevant maps of the mine, convert each piece of information in the relevant maps into coordinate points according to the three-dimensional Cartesian coordinate system, supplement the coordinate points according to the DIS interpolation algorithm, and obtain a three-dimensional mesh based on the supplemented coordinate points.

8. The dynamic modeling method for coal mine digital twins based on multi-source heterogeneous data fusion as described in claim 7, characterized in that, The stratigraphy is constructed based on the GTP algorithm. For the stratigraphy, a construction method is configured, including: using coordinate grid points as reference points, generating a TIN model through point-by-point interpolation, obtaining the stratigraphic sequence, obtaining the top elevation of each layer, setting the top elevation of each layer as the vertex coordinates of the corresponding triangle, setting the bottom elevation of each layer as the bottom coordinates of the corresponding triangle, setting lithology as the internal attribute of the corresponding triangle, modifying the corresponding part of the 3D mesh based on each triangle with set internal attributes, and setting the corresponding part of the modified 3D mesh as the stratigraphy. Obtain the portion corresponding to the fault layer in the 3D mesh, construct the fault surface according to the TIN algorithm, and update the 3D mesh based on the fault surface.

9. The dynamic modeling method for coal mine digital twins based on multi-source heterogeneous data fusion as described in claim 8, characterized in that, According to the preset formula, the virtual particle boundary is calculated. The virtual particle boundary is represented by the virtual particle velocity on the outermost boundary of the wind track topology in the three-dimensional coal mine model. The virtual particle velocity is obtained by Taylor expansion formula. The velocity point cloud of the virtual boundary is obtained based on the velocity of virtual particle b, denoted as the virtual particle boundary. Based on the virtual particle boundary, and the virtual particle boundary corresponding to the first space at the entrance of the 3D coal mine model, the transformation relationship between the velocity of virtual particle b and the coordinates of the corresponding point is obtained. This transformation relationship is expressed as follows: ; Where F1~F3 represent the conversion rules between the velocity of the virtual particle and the x-axis, y-axis, and z-axis coordinates of the corresponding outermost boundary point, respectively; Based on the new three-dimensional spatial region and the three-dimensional coal mine model, the three-dimensional coal mine model is updated. The update is represented by replacing the original three-dimensional spatial region in the three-dimensional coal mine model with the new three-dimensional spatial region.

10. A dynamic modeling system for coal mine digital twins based on multi-source heterogeneous data fusion, the system being used to execute the dynamic modeling method for coal mine digital twins based on multi-source heterogeneous data fusion as described in claim 1, characterized in that, It includes a building module, a computing module, and an update module; The construction module is used to construct a coordinate system, construct independent loops according to a preset construction method, construct wind track topology according to a preset first constraint condition and independent loops, construct a three-dimensional mesh according to the mine-related diagram, and replace the relevant parts of the three-dimensional mesh through the GTP algorithm and TIN algorithm to obtain a three-dimensional coal mine model. The calculation module calculates the virtual particle boundary based on preset ventilation parameters and air duct topology; The update module constructs a coordinate transformation relationship based on the virtual particle boundary, obtains the coordinates of the corresponding points in the three-dimensional coal mine model based on the coordinate transformation relationship, and replaces the original coordinates in the three-dimensional coal mine model based on the coordinates of the points.

Images