CN110675495B

CN110675495B - Coal face coal seam three-dimensional display method based on geologic body modeling and application thereof

Info

Publication number: CN110675495B
Application number: CN201910898457.1A
Authority: CN
Inventors: 张群磊; 李东印; 袁瑞甫; 李化敏; 张旭和; 刘军磊; 王文; 王伸
Original assignee: Henan University of Technology
Current assignee: Henan University of Technology
Priority date: 2019-09-23
Filing date: 2019-09-23
Publication date: 2023-03-24
Anticipated expiration: 2039-09-23
Also published as: CN110675495A

Abstract

The invention provides a geological modeling-based coal face coal seam three-dimensional display method and application thereof, and solves the problem that three-dimensional coal seam space occurrence parameter changes cannot be intuitively, effectively and in real time in the automatic mining process of an extra-thick coal seam working face; the method comprises the steps of processing drilling column diagrams at different positions of an upper gate and a lower gate of a coal face through AutoCAD software, establishing a coal seam roof elevation line and a coal seam floor elevation line, extracting relative coordinates of elevation control points on the coal seam roof elevation line and the coal seam floor elevation line, carrying out interpolation homogenization processing on the elevation control points through geological modeling SUFER software, generating a three-dimensional geological geometry body of the coal face through ANSYS software, and carrying out three-dimensional display through CDEM software. The method can quickly and accurately generate the three-dimensional numerical model of the extra-thick coal seam of the fully mechanized caving face, conveniently display the three-dimensional structure of the extra-thick coal seam, acquire occurrence data of the coal seam at any pushing position of the working face in real time, and provide real-time and accurate geological reference data for dynamically adjusting automatic coal caving control parameters of the extra-thick coal seam.

Description

Coal face coal seam three-dimensional display method based on geologic body modeling and application thereof

Technical Field

The invention relates to the technical field of coal seam numerical modeling, in particular to a coal face coal seam three-dimensional display method based on geological body modeling and application thereof, which are used for realizing the modeling and display of a working face coal seam three-dimensional structure and solving the thickness and the inclination angle of a coal seam.

Background

The safe and efficient production of the fully mechanized top coal caving working face of the extra-thick coal seam depends on the geological occurrence conditions of the working face of the extra-thick coal seam to a great extent, and the occurrence conditions of the top coal of the extra-thick coal seam have great influence on the mining method and the production benefit of the whole mine. The design method of the top coal caving is influenced by the thickness of the coal seam, for example, the design and the model selection of the hydraulic support in the aspect of the support height are influenced by the setting of mining-discharging ratio; the thickness of the top coal affects the selection of the mining mode, for example, the top coal with different thicknesses determines whether the mining mode is a one-time mining full-height mining method, a top coal caving mining method or a layered top coal caving mining method. The coal seam gangue inclusion occurrence characteristics influence the caving characteristic of the top coal and the gangue content of the top coal; the top coals with different thicknesses influence the lumpiness of top coal crushing, the coal discharging time of the bracket, the working resistance of the bracket and the like; the fluctuation of the elevation of the floor of the coal bed affects the position of the roadway excavation. And structural characteristics of each coal rock layer above the coal bed influence the breaking process of the coal rock layer, the settlement of the earth surface, the breaking degree of top coal and the working resistance display characteristics of the bracket. Especially, automatic coal mining and automatic top coal discharging of the current coal mining working face are the main trends of development of the coal industry in future, and real-time dynamic acquisition of the thickness, the inclination angle and the gangue inclusion distribution of an automatic mining coal layer has important guiding significance for determining and dynamically adjusting control parameters of automatic mining equipment.

Therefore, the working face three-dimensional coal rock stratum modeling is carried out through engineering field geological analysis, the three-dimensional space information of the three-dimensional coal rock stratum is obtained in real time to study the occurrence condition of the top coal of the fully mechanized caving working face of the extra-thick coal seam and the occurrence characteristics of the overlying strata, and visual and reasonable geological analysis and modeling basis can be provided for studying the top coal caving rule of the working face of the extra-thick coal seam and an automatic coal caving process. The method has the advantages that the geological occurrence condition of the coal seam of the working face can be intuitively and thoroughly known, the selection and the determination of the coal mining and releasing method of the mine can be better guided, the coal mining process is timely adjusted, and the important influence is generated on the safe and efficient production of the mine.

The application number is 201610643584.3, the Chinese patent application is a working face coal seam three-dimensional modeling method based on geological data, geological data of a coal layer top plate of a coal face is processed through an interpolation method to generate the upper surface of a coal layer structure of the coal face, and coal seam geological data are provided for adjusting the height of a roller of a coal mining machine. However, only the upper surface of the coal seam structure is generated, the elevation fluctuation of the bottom plate of the extra-thick coal seam and the occurrence condition of gangue inclusion in the extra-thick coal seam cannot be considered, the coal seam thickness of the extra-thick coal seam cannot be effectively obtained, and geological occurrence parameters such as the coal seam strike inclination angle, the inclination angle and the like of a certain hydraulic support position of the coal face cannot be reasonably obtained, so that effective and reliable coal seam geological data cannot be well provided for the automatic emission of extra-thick top coal on the coal face of the extra-thick coal seam. The automatic coal mining and the automatic top coal discharging of the coal face are the main trend of the development of the coal industry, and the determination of the automatic control parameters of the automatic mining of the ultra-thick coal seam is very important.

Disclosure of Invention

The invention provides a three-dimensional modeling and three-dimensional model dynamic display method for a working face of an extra-thick coal seam, and a calculation method for calculating the thickness of the coal seam at any advance position of the working face of the extra-thick coal seam and the dip angle of the coal seam, which can correctly guide the coal mining and discharging of a mine, timely adjust the coal mining process, and dynamically acquire the thickness, the dip angle and the dip distribution of an automatic exploitation coal seam in real time, have important guiding significance for determining and dynamically adjusting the control parameters of automatic exploitation equipment, and improve the safe and efficient production of the mine.

In order to achieve the purpose, the technical scheme of the invention is realized as follows: a coal face coal seam three-dimensional display method based on geologic body modeling comprises the following steps:

the method comprises the following steps: drilling and coring are carried out on the top and the bottom of the coal bed at different positions of a transportation gateway and a return air gateway of the coal face, and drilling bar charts at different positions of the upper and lower gateways of the coal face are obtained according to the thickness of a core sample of the coal bed taken at different positions and the height of the gateway;

step two: according to the obtained drilling column diagrams at different positions of the upper and lower crossheading of the coal face, connecting elevation control points on the level of a coal seam roof in AutoCAD software to generate a coal seam roof elevation line, connecting elevation control points on the level of a coal seam floor to generate a coal seam floor elevation line, and extracting relative coordinates of the elevation control points on the coal seam roof elevation line and the coal seam floor elevation line;

step three: utilizing geological modeling SUFER software to respectively carry out interpolation homogenization treatment on elevation control points on a coal seam roof elevation line and a coal seam floor elevation line in a working face propulsion direction and a working face arrangement direction to generate a plurality of control points;

step four: uniformly processing the interpolation by the SUFER software to obtain elevation control points on a coal seam roof and forming the upper surface of a three-dimensional coal seam by a plurality of control points generated in the advancing direction and the arrangement direction of a working surface; uniformly processing the control points on the coal seam floor by using the SUFER software interpolation, and forming the lower surface of the three-dimensional coal seam by using a plurality of control points generated in the crossheading direction and the working face arrangement direction;

step five: respectively generating an upper geometric surface and a lower geometric surface of the coal seam on the basis of elevation control point data of the upper surface and the lower surface of the three-dimensional coal seam obtained after uniform interpolation processing of the SUFER software, establishing a side surface by taking the boundary of the upper geometric surface and the lower geometric surface as a reference to generate a closed space surface, and generating a three-dimensional geological geometry body of a working surface by the closed space surface;

step six: dividing the three-dimensional geological geometry of the working face into small finite element unit grids in ANSYS software by taking a calculation unit body containing spatial information as a basic data unit, and performing three-dimensional display on the model section and the three-dimensional geological geometry of the working face by using CDEM software.

The third step is that the realization method for carrying out interpolation homogenization treatment by the geological modeling SUFER software is a kriging interpolation method; the kriging interpolation method is based on a regionalized variable theory, takes a variation function as a main tool, and obtains an estimated value on the premise of ensuring that the estimated value meets the unbiased condition and the minimum variance condition.

The method for realizing the kriging interpolation method comprises the following steps: if the regionalized variable f (x) satisfies the second-order stationary assumption or the eigen assumption, the estimated value of the point P to be interpolated is

Wherein f is _i Is a function of n known points, w _i Is a full-index and->

And then according to the condition that the estimated variance is minimum: />

Finding an estimated value of the point P to be interpolated, where mu is Lagrangian and gamma (x) _j -x _i ) Is a known point x _i And x _i Function of variation between, gamma (x) _p -x _i ) Is a known point x _i And point x to be interpolated _p I =1,2, … …, n.

The method for forming the upper surface or the lower surface of the three-dimensional coal seam by the multiple control points in the fourth step comprises the following steps: connecting adjacent control points by using a curve to sequentially generate a plurality of line segments, sequentially generating a plurality of small surface areas by the plurality of line segments, and finally, forming the upper surface or the lower surface of each small surface area; the method for generating the closed space surface comprises the following steps: and sequentially generating side surfaces of the three-dimensional coal seam structure based on the space boundary conditions of the upper surface and the lower surface, and generating a closed space surface of the three-dimensional coal seam by the upper surface, the lower surface and the peripheral side surfaces of the three-dimensional coal seam structure.

The application of the geological modeling-based coal face coal seam three-dimensional display method comprises the following steps of: the method comprises the steps of dynamically obtaining space coordinates of any position in a three-dimensional geological geometric body of a working face in CDEM software, and calculating the thickness of the coal seam and coal seam address occurrence parameters of a trend average inclination angle and a trend average inclination angle at any advancing position of the working face by combining control point space coordinates of different positions of the coal seam, which are at a certain distance from a leading working face, at a certain distance from a lagging working face, at a certain distance from a gateway above the working face and at a certain distance from a gateway below the working face, at any advancing position of the working face.

On the basis of the established working face three-dimensional coal seam model, obtaining the plane coordinate of any point in the working face arrangement direction of the working face pushed to a specific position, obtaining the elevation point coordinates of the top surface and the bottom surface of the coal seam at the position through the plane coordinate of the position, and further obtaining the thickness of the coal seam at any point in the working face arrangement direction when the working face is pushed to the specific plane position; the specific implementation method comprises the following steps:

(1) When the coal face is pushed to a specific position X1 along the coal seam trend, the coal face is arranged in the direction, and the plane coordinates of the position of the face at the distance of Y1 from the working face gateway are (X1, Y1);

(2) In a space coordinate point database of the three-dimensional coal seam model, the maximum coordinate value of a space vertical coordinate corresponding to a (X1, Y1) plane coordinate is a point coordinate on the coal seam top surface corresponding to any position (X1, Y1), and the three-dimensional space coordinate of any position (X1, Y1) on the coal seam top surface of the working face is obtained: z1;

(3) In a space coordinate point database of the three-dimensional coal seam model, the minimum coordinate value of a space vertical coordinate corresponding to (X1, Y1) plane coordinates is a point coordinate on the coal seam bottom surface corresponding to any position (X1, Y1), and the three-dimensional space coordinate at any position (X1, Y1) of the coal seam bottom surface of the working face is obtained: z2;

(4) When the thickness of the coal seam at the position where the working face is advanced to the specific plane position (X1, Y1) is obtained by subtracting the minimum coordinate value of the spatial vertical coordinate corresponding to the (X1, Y1) plane coordinate from the maximum coordinate value of the spatial vertical coordinate corresponding to the (X1, Y1) plane coordinate, namely: (Z2-Z1) m.

The working face is pushed to a specific position, three space point coordinates of a coal seam floor in the pushing direction at a lagging pushing position, a pushing current position and a leading pushing position at the top of the coal seam are obtained, the relative coal seam dip angle of the coal seam floor at the pushing position is further obtained, the relative coal seam dip angles of the middle of the coal seam and the bottom of the coal seam at the pushing position are sequentially obtained, and the average coal seam strike dip angle at any pushing position in the coal seam strike direction is obtained by averaging the relative coal seam dip angles at the top of the coal seam, the middle of the coal seam and the bottom of the coal seam.

The method for calculating the average coal seam strike dip angle at any propulsion position in the coal seam strike direction comprises the following steps:

(1) Acquiring a three-dimensional space coordinate of a coal seam top surface corresponding to a position where a working surface is pushed to any position X1 and a working surface crossheading certain distance Y1 in the arrangement direction of the working surface: (X1, Y1, Z1);

(2) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire the three-dimensional spatial coordinates of the top surface of the coal seam corresponding to the position where the working surface is pushed to any position X1+ m and the working surface is arranged at a certain distance Y1 away from the working surface on one side along the trough in the direction of arrangement of the working surface: (Xm, Y1, zm 1); wherein m represents the distance of the leading coal face;

(3) The strike inclination angle of the top surface of the coal seam at the position where the working face is pushed to any position X1 relative to the top of the coal seam at the position m meters ahead of the advanced coal face is as follows: theta _X1 ＝tan ^-1 (|Zm1-Z1|/|Xm-X1|；

(4) The three-dimensional coal seam model acquires that the working face is pushed to any position X1, and the three-dimensional space coordinate of the middle part of the coal seam corresponding to a certain distance Y1 from the working face crossheading on one side in the arrangement direction of the working face is as follows: (X1, Y1, Z2);

(5) According to the method for acquiring the spatial coordinates of one point of the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire the three-dimensional spatial coordinates of the middle part of the coal seam corresponding to the position where the working face is pushed to any position X1+ m and the working face is crossheading at a certain distance Y1 from the arrangement direction of the working face to one side: (Xm, Y1, zm 2);

(6) The strike dip angle of the top surface of the coal seam at any position X1 pushed by the working face relative to the middle part of the coal seam at a distance of m meters in front of the advanced coal face is as follows: theta _X2 ＝tan ^-1 (|Zm2-Z2|/|Xm-X1|；

(7) The three-dimensional coal seam model acquires that the working face is pushed to any position X1, and the three-dimensional space coordinate of the coal seam bottom face corresponding to the position, which is away from the working face crossheading on one side by a certain distance Y1, in the arrangement direction of the working face is as follows: (X1, Y1, Z3);

(8) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire the three-dimensional spatial coordinates of the bottom surface of the coal seam corresponding to the position where the working surface is pushed to any position X1+ m and the working surface is arranged at a certain distance Y1 away from the working surface on one side along the trough in the direction of arrangement of the working surface: (Xm, Y1, zm 3);

(9) The dip angle of the coal seam bottom at the position where the working face is pushed to any position X1 relative to the distance m meters in front of the advanced coal face is as follows: theta _X2 ＝tan ^-1 (|Zm3-Z3|/|Xm-X1|；

(10) And then the average dip angle of the coal seam strike at any position (X1, Y1) is obtained by averaging the dip angle of the coal seam bottom strike, the dip angle of the coal seam middle strike and the dip angle of the coal seam top at any position (X1, Y1) according to the average dip angle of the coal seam strike at any position (X1, Y1) which is pushed by the working surface:

θ _X ＝(θ _X1 +θ _X2 +θ _X3 )/3。

the method comprises the steps that a working face is pushed to a specific position, three space point coordinates of a coal seam roof at any position of a coal seam inclination, at a certain distance from the position close to an upper gate way of the working face and at a certain distance from the position close to a lower gate way of the working face and in the arrangement direction of the working face are obtained, the relative coal seam inclination angle of the coal seam roof at the position is obtained, the relative coal seam inclination angles of the middle portion of the coal seam and the bottom portion of the coal seam at the position are obtained in sequence, and the average coal seam inclination angle of the coal seam at any position in the inclination direction at any pushing position is obtained by averaging the relative coal seam inclination angles of the top portion of the coal seam, the middle portion of the coal seam and the bottom portion of the coal seam.

The calculation method of the average inclined angle of the coal seam at any position of the coal seam at any propulsion position in the inclined direction comprises the following steps:

(1) Acquiring a three-dimensional space coordinate of a coal seam top corresponding elevation point corresponding to a position where a working surface is pushed to any position X1 and a working surface crossheading certain distance Y1 is arranged in the working surface arrangement direction from one side by a three-dimensional coal seam model: (X1, Y1, Z1);

(2) Acquiring a three-dimensional space coordinate of the coal seam top surface at a position (Y1 + n 1) of a hydraulic support position with a certain distance n1 from the Y1 position in the arrangement direction of the working surface and propelling the working surface to any position X1 by using a three-dimensional coal seam model, wherein the three-dimensional space coordinate is as follows: (X1, yn, zn 1); n1 represents the distance from the hydraulic support at the Y1 position in the arrangement direction of the working face;

(3) Then the working face is pushed to any position X1, and at a certain distance n1 position relative to the position of the working face Y1, the inclined dip angle of the top of the coal seam at the position of the working face Y1 is as follows:

θ _Y1 ＝tan ^-1 (|Zn1-Z1|/|Yn-Y1|)；

(4) Acquiring a three-dimensional space coordinate of a coal seam middle corresponding elevation point corresponding to a position where a working surface is pushed to any position X1 and a working surface crossheading certain distance Y1 in the arrangement direction of the working surface from a three-dimensional coal seam model: (X1, Y1, Z2);

(5) The three-dimensional space coordinate of the middle part of the coal seam at the position (Y1 + n 1) of the hydraulic support at a certain distance n1 from the Y1 position in the arrangement direction of the working face is obtained by the three-dimensional coal seam model, wherein the working face is pushed to any position X1: (X1, yn, zn 2);

(6) Then the working face is pushed to any position X1, and at a certain distance n from the position of the working face Y1, the inclined dip angle of the middle part of the coal seam at the position of the working face Y1 is as follows:

θ _Y2 ＝tan ^-1 (|Zn2-Z2|/|Yn-Y1|)；

(7) Acquiring a three-dimensional space coordinate of a coal bed bottom corresponding elevation point corresponding to a position where a working surface is pushed to any position X1 and a working surface crossheading certain distance Y1 is arranged in the working surface arrangement direction from a three-dimensional coal bed model: (X1, Y1, Z3);

(8) The three-dimensional space coordinates of the bottom of the coal seam at the position (Y1 + n 1) of the hydraulic support with a certain distance n1 from the Y1 position in the arrangement direction of the working face are obtained by the three-dimensional coal seam model, wherein the three-dimensional space coordinates are as follows: (X1, yn, zn 3);

(9) Then the working surface is pushed to any position X1, and at a certain distance n1 relative to the position Y1 of the working surface, the inclined inclination angle of the bottom of the coal seam at the position Y1 of the working surface is as follows:

θ _Y3 ＝tan ^-1 (|Zn3-Z3|/|Yn-Y1|)；

(10) The working surface is pushed to any position X1 and a position with a certain distance n1 relative to the position of the working surface Y1, the average dip angle of the coal seam inclination at the position of the working surface Y1 is obtained by averaging the dip angle of the coal seam bottom, the dip angle of the coal seam middle and the dip angle of the coal seam top:

θ _Y ＝(θ _Y1 +θ _Y2 +θ _Y3 )/3。

compared with the prior art, the invention has the beneficial effects that: processing drilling column diagrams at different positions of an upper gate and a lower gate of a coal face through AutoCAD software, establishing a coal seam roof elevation line and a coal seam floor elevation line, extracting relative coordinates of elevation control points on the coal seam roof elevation line and the coal seam floor elevation line, performing interpolation homogenization processing on the elevation control points through geological modeling SUFER software, generating a three-dimensional geological geometry body of the coal face through ANSYS software, performing three-dimensional display through CDEM software, and conveniently obtaining space coordinates of any position in the three-dimensional geological geometry body in the follow-up process; the thickness of the coal seam at any advancing position of the working face, the coal seam address occurrence parameters such as the average trend inclination angle and the average inclination angle of the inclined coal seam can be conveniently calculated, and the thickness of the coal seam at any supporting position of the hydraulic support at a specific position where the working face is advanced, the inclination angle of the coal seam in the advancing direction relative to the working face and the inclination angle of the coal seam in the arrangement direction relative to the working face can be obtained. The invention can quickly and accurately generate the three-dimensional numerical model of the ultra-thick coal bed of the fully mechanized caving face, conveniently display the three-dimensional structure of the ultra-thick coal bed layer, acquire occurrence data of the coal bed at any advancing position of the working face in real time, and provide real-time and reliable coal bed geological occurrence parameters for realizing the automatic coal mining and coal discharging operations of dynamically adjusting the coal discharging opening time of a specific bracket, the advancing layout of a hydraulic bracket, the swing amplitude and the swing frequency of a bracket tail beam, adjusting the initial supporting force of the hydraulic bracket, the opening pressure of a safety valve and the like.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the prior art descriptions will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.

FIG. 1 is a flow chart of the present invention.

FIG. 2 is a schematic diagram of an elevation line generated by using a coal seam drilling histogram of the present invention, wherein (a) is a 5202 roadway drilling histogram, and (b) is a 2202 roadway drilling histogram.

FIG. 3 is a schematic view of the surface structure of a coal seam of a working face of the present invention, wherein (a) is the upper surface and (b) is the lower surface.

Fig. 4 is a schematic diagram of a three-dimensional geometric model of a working surface according to the present invention, wherein (a) is an upper surface and a lower surface of a three-dimensional geometric body, (b) is a lateral surface of the three-dimensional geometric body, (c) is a top view of the three-dimensional geometric body, and (d) is a bottom view of the three-dimensional geometric body.

FIG. 5 is a schematic diagram of a three-dimensional numerical model of a working surface according to the present invention, wherein (a) is a three-dimensional coal seam grid model and (b) is a three-dimensional coal seam grid node.

Fig. 6 is a schematic diagram of a three-dimensional display model of a working surface according to the present invention, in which (a) one-point information acquisition of a three-dimensional coal seam space is performed, and (b) overall information acquisition of a three-dimensional coal seam is performed.

Detailed Description

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without inventive effort based on the embodiments of the present invention, are within the scope of the present invention.

As shown in fig. 1, a method for three-dimensionally displaying a coal seam structure of a coal face based on geologic body modeling and calculating the thickness and inclination angle of the coal seam,

the method comprises the following steps: drilling and coring are carried out on the top and the bottom of the coal bed at different positions of the transportation gateway and the return air gateway of the coal face, and drilling bar charts at different positions of the upper and lower gateways of the coal face are obtained according to the thickness of the core sample of the coal bed taken at different positions and the height of the gateway.

Drilling and coring are carried out on the top and the bottom of the coal seam at different positions of the crossheading of the coal face, and drilling bar charts at different positions of the crossheading of the coal face are obtained according to the thickness of a coal seam core sample taken at different positions of the crossheading and the height of the crossheading.

Drilling and coring are carried out on the top and the bottom of a coal bed at different positions of the same Xin coal mine second panel 8202 coal face 5202 crossheading 2-4 and the same Xin coal mine second panel 2202 crossheading 2-9, and drilling column diagrams at different positions of the crossheading of the working face are obtained according to the thickness of a coal bed core sample taken at different positions of the crossheading and the height of the crossheading at different positions, so that a CAD format file is established.

The drilling histogram is a CAD file, after core drilling and sampling on site, the heights of different rock core samples are respectively measured and then drawn in CAD software, and the coal bed top plate contour line and the coal bed bottom plate contour line are formed by connecting the centers of the coal rock interfaces in each histogram in sequence by using a plurality of lines in the CAD drilling histogram.

Step two: according to the obtained drilling column diagrams at different positions of the upper and lower crossheading of the coal face, connecting elevation control points on the horizon of the coal seam roof in the AutoCAD software to generate a coal seam roof elevation line, connecting elevation control points on the horizon of the coal seam floor to generate a coal seam floor elevation line, and extracting relative coordinates of the elevation control points on the coal seam roof elevation line and the coal seam floor elevation line.

The drilled histogram and top and bottom plate elevation lines for the 5202 lane and the 2202 lane are shown in fig. 2. Corresponding upper and lower elevation control points 2-10 and 2-11 of the coal seam are found out according to the position of the coal seam shown in the drilling histogram 2-1 of the upper crossheading, similarly, corresponding elevation control points are found out through analysis of all drilling histograms, elevation control points on the surface of the coal seam roof are connected in the AutoCAD software to generate elevation lines 2-2 and 2-5 of the coal seam roof, as shown in figure 2 (a), elevation control points on the surface of the coal seam floor are connected to generate elevation lines 2-3 and 2-6 of the coal seam floor, as shown in figure 2 (b). And connecting elevation control points on the surfaces of the top plate and the bottom plate of the coal seam in the AutoCAD software to generate elevation control lines of the bottom plate and the top plate of the coal seam, and extracting to obtain relative coordinates of the elevation control points of the top plate and the bottom plate of the coal seam of the upper and lower crossheading seams.

In order to research the thickness change of a three-dimensional coal seam on a working face, the left lower corner of a drilling histogram of the 5202 lane in the first step is moved to the origin of coordinates, other layers are hidden, only the layer of the elevation control line of the top floor and the bottom floor of the coal seam is reserved, the elevation control line of the top floor and the bottom floor is independently stored as a dxf format file of CAD software, and the spatial coordinates of the intersection points of the top floor and the bottom floor of the coal seam in the elevation line of the top floor of the coal seam and all drilling histograms of the 5202 lane are obtained through a dxf format multi-line control point coordinate extraction program. And moving the lower left corner of the 2202 th lane drilling histogram in the first step to the origin of coordinates, hiding other layers, only reserving the layer of the coal seam top and bottom plate elevation control line, independently storing the top and bottom plate elevation control line as a dxf format file of CAD software, and obtaining the space coordinates of the intersection points of the coal seam top and bottom plate levels in all the 2202 th lane drilling histogram through a dxf format multi-line control point coordinate extraction program. The relative elevation coordinates of the elevation line control points of the 5202 lane and the 2202 lane are shown in tables 1 and 2.

Table 1 5202 elevation control points of roadway coal seam roof and floor

Table 2 elevation control points of 2202 lane coal seam roof and floor

From tables 1 and 2, it can be seen that the coal seam roof elevation control points generated by 5202 lane drilling are 14 data, the coal seam roof elevation control points generated by 2202 lane drilling are 21 data, the coal seam floor elevation control points generated by 5202 lane drilling are 14 data, and the coal seam floor elevation control points generated by 2202 lane drilling are 21 data.

Step three: and utilizing geological modeling SUFER software to respectively carry out interpolation homogenization treatment on the elevation control points on the coal seam roof elevation line and the coal seam floor elevation line in the working face propulsion direction and the working face arrangement direction to generate a plurality of control points.

Interpolation processing is carried out on the control points of the top and bottom of the coal bed obtained by utilizing drilling data in the transportation crossheading and the return air crossheading of the geological modeling SUFER software working face, interpolation homogenization processing of the height control points of the top and bottom of the coal bed between two crossheadings is realized, namely, interpolation homogenization processing is respectively carried out on the height control points on the height lines of the top and bottom of the coal bed in the advancing direction and the arrangement direction of the working face, and a plurality of height control points in the crossheading direction and the arrangement direction of the coal face are further generated by supplementing and encrypting, so that the height control points of the coal bed bottom and the coal bed trend of the top of the upper and lower crossheading and the height control points of the coal bed trend between the upper and lower crossheading in the arrangement direction of the working face are enabled to present a change trend of smooth transition on the whole. Carrying out interpolation homogenization treatment on elevation control points on a coal seam roof elevation line and a coal seam floor elevation line in the axial direction of each crossheading by geological modeling SUFER software, and increasing and generating a plurality of control points in the crossheading direction, so that the coal seam trend elevation control points of the coal seam floor and the coal seam roof of the upper crossheading and the lower crossheading are in a smoothly-transitional change trend on the whole; according to the elevation control points of the upper and lower crossheading at the same working face propulsion position, interpolation homogenization processing is carried out in the coal face arrangement direction to generate a plurality of control points in the working face arrangement direction, so that the inclined elevation control points of the coal seam roof and floor layers between the upper and lower crossheading in the working face arrangement direction are in a smooth transition change trend on the whole.

The realization method for carrying out interpolation homogenization treatment by geological modeling SUFER software is a kriging interpolation method; the kriging interpolation method is based on a regionalized variable theory, takes a variation function as a main tool, and obtains an estimated value on the premise of ensuring that the estimated value meets the unbiased condition and the minimum variance condition.

Wherein f is _i Is a function of n known points, w _i Is a full-index and->

And according to the condition that the estimated variance is minimum: />

For a coal face with the push length of about 2000m, in order to accurately know the three-dimensional occurrence condition of a coal bed on the whole face, the data volume of only 35 control points on the surface of the top floor of the coal bed is obviously insufficient, so that the interpolation homogenization treatment is performed by using geological modeling SUFER software, and the number of the three-dimensional coal bed control points is reasonably increased.

Reading in coal seam roof elevation line control point Data interpolation of 5202 gate and 2202 gate in Excel table 1 and Excel table 2 by a Grid-Data module of the SUFER software to generate the three-dimensional coal seam upper surface of the 8202 coal face. The concrete implementation process of the geological modeling SUFER software is that interpolation processing is firstly carried out in the coal seam trend direction according to a Krigin interpolation method, and the concrete difference parameters are as follows: the working face advancing direction is 2000m, the interpolation interval is 20m, and the interpolation intervals are 100; secondly, interpolation processing is carried out in the coal seam inclination direction, and specific difference parameters are as follows: because the length of a single working face is only 200m, the change of the coal bed is relatively smooth, and two control point linear interpolation values are adopted in the coal bed trend direction; similarly, a Grid-Data module of the SUFER software reads in coal seam floor elevation line control Data of 5202 gateway and 2202 gateway in table 1 and Excel table 2 in Excel format to generate the three-dimensional coal seam lower surface of the 8202 coal face by interpolation, the specific implementation process comprises the following steps of firstly carrying out interpolation processing in the coal seam trend direction according to a Kriging interpolation method, and the specific difference parameters are as follows: the working face advancing direction is 2000m, the interpolation interval is 20m, and the interpolation intervals are 100; secondly, interpolation processing is carried out in the coal seam inclination direction, and specific difference parameters are as follows: because the length of a single working face is only 200m, the change of the coal bed is relatively gentle, and therefore two control points are linearly interpolated in the coal bed inclination direction. 35 coal seam control point data extracted by drilling of AutoCAD software are subjected to interpolation homogenization processing by utilizing SUFER software, the upper surface and the lower surface of the coal seam can be subdivided into 100 control intervals after being processed by the SUFER software, and 70 control points on the surface of the top floor plate of the coal seam are increased to 400 control points, so that the generated surface model of the top floor plate of the coal seam is closer to the real coal seam condition.

Step four: based on an elevation control point on a coal seam roof and a plurality of control points generated in the working face advancing direction and the working face arrangement direction after uniform processing by a SURER software Krigin interpolation method, automatically connecting adjacent control points to sequentially generate a plurality of line segments in a Map-New-3D surface module of the surfer software, sequentially generating a plurality of small surface areas by the plurality of line segments, and finally forming the upper surface of the three-dimensional coal seam by all the small surface areas; the method comprises the steps that control points on a coal seam floor and a plurality of control points generated in the crossheading direction and the working face arrangement direction are obtained after interpolation uniform processing based on the SUFFER software, a plurality of line segments are generated in sequence by automatically connecting adjacent control points in a Map-New-3D surface module of the surfer software, a plurality of small surface areas are generated in sequence by the line segments, and all the small surface areas finally form a three-dimensional coal seam to form the lower surface of the three-dimensional coal seam. The spatial forms of the upper surface and the lower surface of the three-dimensional coal seam structure subjected to interpolation processing by the SUFER software are shown in fig. 3 (a) and 3 (b).

Step five: based on elevation control point data of the upper surface and the lower surface of the three-dimensional coal seam obtained after the SUFER software interpolation is uniformly processed, the upper geometric surface and the lower geometric surface of the coal seam are respectively generated in ANSYS software, the side face is established by taking the boundary of the upper geometric surface and the lower geometric surface as the reference to generate a closed space surface, and the three-dimensional geological geometry body of the working face is generated by the closed space surface.

The occurrence state of the working face coal bed three-dimensional space is researched, so that a working face coal bed three-dimensional space geometrical body needs to be established. Converting three-dimensional coordinate data of 400 elevation control points of the upper surface and the lower surface of a three-dimensional coal seam structure obtained after interpolation processing of the SUFFER software into K-point geometric key point data of ANSYS software through the Sufer-to-ANSYS software, sequentially generating a plurality of line segments through APDL language programming processing, sequentially generating a plurality of small surface areas by the plurality of line segments, finally respectively forming the upper surface and the lower surface of the coal seam structure area by all the small surface areas, establishing a side surface by taking the boundary of the upper surface and the lower surface as a reference to generate a closed space surface, and generating a three-dimensional geological geometry of a working surface by the closed space surface.

The method for generating the closed space surface comprises the following steps: based on the space boundary conditions of the upper surface and the lower surface, side surfaces of a three-dimensional coal seam structure are sequentially generated in a Preferences-Create-Area module of ANSYS software through a height control line of a top floor and a bottom floor of a coal seam of the same gate, and a space geometry of the three-dimensional coal seam is generated in the Preferences-Create-Volume module of the ANSYS software through a closed space surface of the three-dimensional coal seam, wherein the closed space surface of the three-dimensional coal seam is composed of the upper surface, the lower surface and the peripheral side surfaces of the three-dimensional coal seam structure.

Based on the relative coordinate data of the elevation control points of the upper surface and the lower surface of the three-dimensional coal seam obtained after interpolation processing of the SUFFER software, the upper geometric surface and the lower geometric surface of the coal seam are respectively generated through APDL language programming processing of ANSYS software, the boundary of the upper geometric surface and the lower geometric surface is used as a standard for establishing a side face to generate a closed space surface, and the three-dimensional geological geometry of the working face is obtained.

In ANSYS software, the upper surface and the lower surface of a 8202-surface coal layer 8202 coal face three-dimensional coal layer generated by 400 elevation control points of the upper surface and the lower surface of the coal layer subjected to interpolation processing by the SUFFER software are shown in fig. 4 (a), the side surface of a three-dimensional coal layer structure established by taking the upper geometric surface and the lower geometric surface of the 8202-surface coal layer as the reference is shown in fig. 4 (b), and the upper view and the lower view of a working face three-dimensional geological body geometric body generated by closing the upper surface and the lower surface and the peripheral side surface of the coal layer are respectively shown in fig. 4 (c) and (d).

Step six: dividing the three-dimensional geological geometry body of the working face into small finite element unit grids in ANSYS software by taking a computing unit body containing spatial information as a basic data unit, and performing three-dimensional display on the cross section of the model and the three-dimensional geological geometry body of the working face by CDEM software.

In order to accurately obtain coordinate information of any point space position in the three-dimensional geological body of the working face, the three-dimensional geological geometric body of the working face needs to be subjected to grid division by using ANSYS software and taking a calculation unit body containing space information as a basic data unit, the unit type selected in the process of dividing the units by the ANSYS software is SOLID185, the unit types are tetrahedral units, the number of the units is 30461, the three-dimensional structural grid of the coal seam of the working face containing the coordinate information after the units are divided is shown in a graph (a) of figure 5, and the node distribution of the unit grid is shown in a graph (b) of figure 5.

In order to accurately and conveniently acquire the geometric information of the three-dimensional geological body of the working face, the CDEM software is used as three-dimensional geological body display software, and the CDEM software can conveniently acquire the coordinate information of the model unit and carry out three-dimensional display on the model section and the space body. The method specifically comprises the steps that a coal seam three-dimensional structure grid model after the unit division is converted into an ANSYS-to-CDEM.txt File through a File-Read Input from module of ANSYS software, the ANSYS-to-CDEM.txt File is converted into an ansys.dat format text File which can be identified by CDEM, the ANSYS-to-CDEM software unit generation-ANSYS module reads in the text File, three-dimensional display of a model section and a spatial body can be conducted, and spatial coordinates of any position in the model can be picked up through a mouse at will. The three-dimensional coal seam structure morphology and any point coordinate information displayed in the CDEM software are shown in fig. 6 (a). From figure 6 (b), it can be seen that said coal seam on 8202 face is about 13-27m in thickness, and the coal seam thickness in local area is about 30 m. The three-dimensional model display and unit coordinate acquisition functions based on the CDEM software can conveniently and visually acquire the coal and gangue stratum distribution and the space coordinates of any position in the three-dimensional coal seam model.

For the coal face of the ultra-thick coal seam, at different pushing positions, the thickness and the inclination angle of the coal seam can affect the working states of the hydraulic support and the coal mining machine and can also significantly affect the discharging effect of top coal, so that after the working face coal seam three-dimensional model is built, the thickness, the gangue inclusion occurrence and the coal seam inclination angle of the coal seam at any pushing position of the working face can be efficiently and automatically obtained, and the method has important significance for safely and efficiently carrying out coal mining. And calculating occurrence parameters such as coal seam thickness, dip angle and the like of the working face of the extra-thick coal seam at different propulsion positions based on the space coordinates of any position in the three-dimensional model of the space body obtained in the step six, and providing more accurate geological data information for automatic coal mining and coal caving of the extra-thick coal seam.

The method for acquiring the thickness of the coal seam pushed to a specific position by the working face comprises the following steps:

(1) On the basis of the established working face three-dimensional coal seam model, when the coal face is pushed to a specific position X1 along the coal seam trend, the arrangement direction of the coal face is that the plane coordinate of the working face position at the distance Y1 from the working face crossheading is (X1, Y1).

(2) In a space coordinate point database of the three-dimensional coal seam model, the maximum coordinate value of a space vertical coordinate corresponding to a (X1, Y1) plane coordinate is a point coordinate on the coal seam top surface corresponding to any position (X1, Y1), so that the three-dimensional space coordinate at any position (X1, Y1) on the coal seam top surface of the working surface can be efficiently obtained by indexing: and Z1.

(3) In a space coordinate point database of the three-dimensional coal bed model, the minimum coordinate value of a space vertical coordinate corresponding to a (X1, Y1) plane coordinate is a coordinate of one point on the coal bed bottom surface corresponding to any position (X1, Y1), so that the three-dimensional space coordinate of any position (X1, Y1) on the coal bed bottom surface of the working surface can be efficiently obtained by indexing: and Z2.

The method for acquiring the coal seam inclination angle in the advancing direction of the working face when the working face is advanced to a specific position comprises the following steps:

(1) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional spatial coordinates of the top surface of the coal seam corresponding to the position where the working surface is pushed to any position X1 and the working surface is at a certain distance Y1 from the working surface on one side in the arrangement direction of the working surface can be efficiently acquired by the three-dimensional coal seam model: (X1, Y1, Z1).

(2) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire the three-dimensional spatial coordinates of the top surface of the coal seam corresponding to the position where the working surface is pushed to any position X1+ m and the working surface is arranged at a certain distance Y1 away from the working surface on one side along the trough in the direction of arrangement of the working surface: (Xm, Y1, zm 1).

The coal seam strike inclination refers to the coal seam inclination of the advancing direction of the coal face when the working face is advanced to the position X1, and is obtained by the ratio of the elevation difference value of the advancing direction position at the advancing position of the current working face and a certain distance ahead of the coal face to the horizontal distance of the advancing direction. m represents a certain distance m meters from the leading coal face, and the coal seam inclination angle at a specific position Y1 (namely the position of the hydraulic support with a specific number of the working face) in the arrangement direction of the working face is obtained, so that the coordinate of Y1 is unchanged.

(3) The strike dip angle of the top surface of the coal seam at the position where the working face is pushed to any position X1 relative to the top of the coal seam at the position m meters ahead of the advanced coal face is as follows: theta.theta. _X1 ＝tan ^-1 (|Zm1-Z1|/|Xm-X1|)。

(4) According to the method for acquiring the spatial coordinates of one point of the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire that the working face is pushed to any position X1, and the three-dimensional spatial coordinates of the middle part of the coal seam corresponding to a position with a certain distance Y1 from the working face crossheading on one side in the arrangement direction of the working face are as follows: (X1, Y1, Z2).

(5) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire the three-dimensional spatial coordinates of the middle part of the coal seam corresponding to the position where the working face is pushed to any position X1+ m and the working face is crossheading at a certain distance Y1 from the arrangement direction of the working face: (Xm, Y1, zm 2).

(6) The strike inclination angle of the top surface of the coal seam at any position X1 pushed by the working face relative to the middle part of the coal seam at a distance of m meters in front of the advanced coal face is as follows: theta.theta. _X2 ＝tan ^-1 (|Zm2-Z2|/|Xm-X1|)。

(7) According to the method for acquiring the spatial coordinates of one point of the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire that the working face is pushed to any position X1, and the three-dimensional spatial coordinates of the bottom surface of the coal seam corresponding to a position with a certain distance Y1 from the working face crossheading on one side in the arrangement direction of the working face are as follows: (X1, Y1, Z3).

(8) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire the three-dimensional spatial coordinates of the bottom surface of the coal seam corresponding to the position where the working surface is pushed to any position X1+ m and the working surface is arranged at a certain distance Y1 away from the working surface on one side along the trough in the direction of arrangement of the working surface: (Xm, Y1, zm 3).

(9) The dip angle of the coal bed bottom, at the position where the working face is pushed to any position X1, of the coal bed top face relative to the distance m meters in front of the advanced coal face is as follows: theta.theta. _X2 ＝tan ^-1 (|Zm3-Z3|/|Xm-X1|)。

(10) And then the coal seam strike average dip angle of the working face at any position (X1, Y1) is pushed to the coal seam bottom strike dip angle at any position (X1, Y1) according to the working face, and the coal seam middle strike dip angle and the coal seam top strike dip angle are averaged to obtain:

θ _X ＝(θ _X1 +θ _X2 +θ _X3 )/3。

the method for acquiring the inclination angle of the coal seam in the working face arrangement direction when the working face is pushed to a specific position comprises the following steps:

(1) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional spatial coordinates of the corresponding elevation point at the top of the coal seam corresponding to the position where the working surface is pushed to any position X1 and the working surface is arranged at a certain distance Y1 away from the working surface on one side along the trough in the arrangement direction of the working surface can be efficiently acquired by the three-dimensional coal seam model: (X1, Y1, Z1).

(2) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire the three-dimensional spatial coordinates of the top surface of the coal seam at a position (Y1 + n 1) of the hydraulic support, wherein the position (Y1 + n 1) is a certain distance n1 from the position Y1 in the arrangement direction of the working surface, and the three-dimensional spatial coordinates of the top surface of the coal seam are as follows: (X1, yn, zn 1).

The dip angle of the coal seam inclination refers to the dip angle of the coal seam in the arrangement direction of the coal face when the working face is pushed to the position X1, and is obtained by the ratio of the elevation difference value to the horizontal distance at different positions (the positions of the hydraulic supports with different numbers on the working face) on the working face. n1 represents the position of the hydraulic support at a certain distance n1 from the position Y1 in the arrangement direction of the working face, and the inclined inclination angle of the coal bed at the position Y1 in the arrangement direction of the working face (namely the position of the hydraulic support with a specific number of the working face) is obtained when the working face is pushed to the position X1, so that the coordinate of X1 is not changed.

(3) Then the working surface is pushed to any position X1, and at a certain distance n1 relative to the position of the working surface Y1, the inclined inclination angle of the top of the coal seam at the position of the working surface Y1 is as follows:

θ _Y1 ＝tan ^-1 (|Zn1-Z1|/|Yn-Y1|)；

(4) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional spatial coordinates of the corresponding elevation point in the middle of the coal seam corresponding to the position where the working face is pushed to any position X1 and the working face is at a certain distance Y1 from the working face on one side in the arrangement direction of the working face can be efficiently acquired by the three-dimensional coal seam model: (X1, Y1, Z2).

(5) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire the three-dimensional spatial coordinates of the middle part of the coal seam at the position (Y1 + n 1) of the hydraulic support with a certain distance n1 from the position Y1 in the arrangement direction of the working face to the position X1 of the working face, wherein the three-dimensional spatial coordinates are as follows: (X1, yn, zn 2).

(6) Then the working surface is pushed to any position X1, and at a certain distance n1 relative to the position Y1 of the working surface, the inclined inclination angle of the middle part of the coal seam at the position Y1 of the working surface is as follows:

θ _Y2 ＝tan ^-1 (|Zn2-Z2|/|Yn-Y1|)。

(7) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional spatial coordinates of the corresponding elevation point at the bottom of the coal seam corresponding to the position where the working surface is pushed to any position X1 and the working surface is at a certain distance Y1 from the working surface on one side along the trough in the arrangement direction of the working surface can be efficiently acquired by the three-dimensional coal seam model: (X1, Y1, Z3).

(8) According to the method for acquiring the spatial coordinates of one point of the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire the three-dimensional spatial coordinates of the bottom of the coal seam at the position (Y1 + n 1) where the working face is pushed to any position X1 and the position of the hydraulic support is a certain distance n away from the position Y1 in the arrangement direction of the working face, wherein the three-dimensional spatial coordinates are as follows: (X1, yn, zn 3).

θ _Y3 ＝tan ^-1 (|Zn3-Z3|/|Yn-Y1|)；

(10) Then the working surface is pushed to any position X1, the position with a certain distance n1 relative to the position of the working surface Y1 is provided, the coal seam inclined average inclination angle at the position of the working surface Y1 is the coal seam bottom inclined inclination angle, and the coal seam middle inclined inclination angle and the coal seam top inclined inclination angle are obtained by averaging:

θ _Y ＝(θ _Y1 +θ _Y2 +θ _Y3 )/3。

the thickness of the coal bed at the supporting position of any hydraulic support of the working face, the coal bed inclination angle of the coal bed relative to the advancing direction of the working face and the coal bed inclination angle of the coal bed relative to the arrangement direction of the working face, which can be obtained by the method, of the working face can be obtained, and the geological occurrence information parameters of the coal bed at the specific position are input into a sensing controller arranged on the hydraulic support of the working face in real time, so that the opening time of the coal discharge port of the specific support, the advancing arrangement of the hydraulic support, the swing amplitude and the swing frequency of the tail beam of the support, the initial supporting force of the hydraulic support, the opening pressure of a safety valve and the like can be dynamically adjusted.

The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims

1. A coal face coal seam three-dimensional display method based on geological body modeling is characterized by comprising the following steps:

2. The geological modeling based coal face coal seam three-dimensional display method as claimed in claim 1, wherein the geological modeling refer software in step three is used for interpolation homogenization treatment by a kriging interpolation method; the kriging interpolation method is based on a regionalized variable theory, takes a variation function as a main tool, and obtains an estimated value on the premise of ensuring that the estimated value meets the unbiased condition and the minimum variance condition.

3. The geological modeling-based coal face coal seam three-dimensional display method according to claim 2, wherein the kriging interpolation method is implemented by the following steps: if the regionalized variable f (x) satisfies the second-order stationary assumption or the eigen assumption, the estimated value of the point P to be interpolated is

Wherein f is _i Is a function of n known points, w _i Is the full coefficient of n known points and

and then according to the condition that the estimated variance is minimum:

4. The geological modeling based coal face coal seam three-dimensional display method according to claim 3, wherein the method for combining the plurality of control points in the fourth step into the upper surface or the lower surface of the three-dimensional coal seam is as follows: connecting adjacent control points by using a curve to sequentially generate a plurality of line segments, sequentially generating a plurality of small surface areas by the plurality of line segments, and finally, forming the upper surface or the lower surface of each small surface area; the method for generating the closed space surface comprises the following steps: and sequentially generating side surfaces of the three-dimensional coal bed structure based on the space boundary conditions of the upper surface and the lower surface, and generating a closed space surface of the three-dimensional coal bed by the upper surface, the lower surface and the peripheral side surfaces of the three-dimensional coal bed structure.

5. The application of the geological modeling based coal face coal seam three-dimensional display method according to claim 1 or 4, characterized in that the thickness and inclination angle of the coal seam are calculated as follows: the method comprises the steps of dynamically obtaining space coordinates of any position in a three-dimensional geological geometric body of a working face in CDEM software, and calculating the thickness of the coal seam and coal seam address occurrence parameters of a trend average inclination angle and a trend average inclination angle at any advancing position of the working face by combining control point space coordinates of different positions of the coal seam, which are at a certain distance from a leading working face, at a certain distance from a lagging working face, at a certain distance from a gateway above the working face and at a certain distance from a gateway below the working face, at any advancing position of the working face.

6. The application of the geological modeling-based coal face coal seam three-dimensional display method is characterized in that on the basis of the established working face three-dimensional coal seam model, the plane coordinates of any point in the working face arrangement direction of the working face pushed to a specific position are obtained, the elevation point coordinates of the top face and the bottom face of the coal seam at the position are obtained through the plane coordinates of the position, and the thickness of the coal seam at any point in the working face arrangement direction when the working face is pushed to the specific plane position is further obtained; the specific implementation method comprises the following steps:

7. The application of the geological modeling-based coal face coal seam three-dimensional display method is characterized in that the working face is pushed to a specific position, three space point coordinates of a coal seam roof in the pushing direction at the lagging pushing position, the pushing current position and the advancing pushing position of the top of the coal seam are obtained, the relative coal seam inclination angle of the bottom of the coal seam at the pushing position is further obtained, the relative coal seam inclination angles of the middle of the coal seam and the bottom of the coal seam at the pushing position are sequentially obtained, and the average coal seam strike inclination angle at any pushing position in the coal seam strike direction is obtained by averaging the relative coal seam inclination angles of the top of the coal seam, the middle of the coal seam and the bottom of the coal seam.

8. The application of the geological modeling based coal face coal seam three-dimensional display method according to claim 7, wherein the average coal seam strike dip at any advance position of the coal seam strike direction is calculated by:

(5) According to the method for acquiring the spatial coordinates of one point on the top surface of the coal seam, the three-dimensional coal seam model can efficiently acquire the three-dimensional spatial coordinates of the middle part of the coal seam corresponding to the position where the working face is pushed to any position X1+ m and the working face is crossheading at a certain distance Y1 from the arrangement direction of the working face: (Xm, Y1, zm 2);

(6) The working face is pushed to the top surface of the coal seam at any position X1 relative to the middle trend of the coal seam at the position m meters ahead of the advanced coal mining working faceThe inclination angle is: theta _X2 ＝tan ^-1 (|Zm2-Z2|/|Xm-X1|；

(10) And then the average dip angle of the coal seam strike at any position (X1, Y1) pushed by the working face is obtained by averaging the dip angle of the coal seam bottom strike, the dip angle of the coal seam middle strike and the strike dip angle of the coal seam top at any position (X1, Y1) pushed by the working face:

θ _X ＝(θ _X1 +θ _X2 +θ _X3 )/3。

9. the application of the geological modeling-based coal face coal seam three-dimensional display method is characterized in that the working face is pushed to a specific position, three spatial point coordinates of a coal seam roof in the arrangement direction of the working face at any position of the coal seam inclination, a position close to an upper crossheading of the working face with respect to the position and a position close to a lower crossheading of the working face with respect to the position are obtained, the relative coal seam inclination angle of a coal seam floor at the position is obtained, the relative coal seam inclination angles of the middle part of the coal seam and the bottom part of the coal seam at the position are obtained in sequence, and the average coal seam inclination angle of the coal seam at any position in the inclination direction at any pushing position is obtained by averaging the relative coal seam inclination angles of the top part of the coal seam, the middle part of the coal seam and the bottom part of the coal seam.

10. The application of the geological modeling-based coal face coal seam three-dimensional display method is characterized in that the calculation method of the average coal seam inclination angle of the coal seam at any position in the inclination direction at any advancing position is as follows:

θ _Y1 ＝tan ^-1 (|Zn1-Z1|/|Yn-Y1|；

(6) Then the working surface is pushed to any position X1, and at a certain distance n relative to the position Y1 of the working surface, the inclined dip angle of the middle part of the coal seam at the position Y1 of the working surface is as follows:

θ _Y2 ＝tan ^-1 (|Zn2-Z2|/|Yn-Y1|；

(7) Acquiring a three-dimensional space coordinate of a coal bed bottom corresponding elevation point corresponding to a position where a working surface is pushed to any position X1 and a working surface crossheading certain distance Y1 in the arrangement direction of the working surface by using a three-dimensional coal bed model: (X1, Y1, Z3);

(8) Acquiring a three-dimensional space coordinate of the bottom of the coal seam at a position (Y1 + n 1) of a hydraulic support at a certain distance n1 from the Y1 position in the arrangement direction of the working surface and propelling the working surface to any position X1 by using a three-dimensional coal seam model, wherein the three-dimensional space coordinate is as follows: (X1, yn, zn 3);

θ _Y3 ＝tan ^-1 (|Zn3-Z3|/|Yn-Y1|；

(10) Then the working surface is pushed to any position X1, and at a certain distance n1 relative to the position Y1 of the working surface, the average inclination angle of the coal seam trend at the position Y1 of the working surface is obtained by averaging the inclination angle of the bottom of the coal seam, the inclination angle of the middle of the coal seam and the inclination angle of the top of the coal seam:

θ _Y ＝(θ _Y1 +θ _Y2 +θ _Y3 )/3。