CN107169616B - Relative entropy prediction method for relative complexity of mine non-mining area structure - Google Patents
Relative entropy prediction method for relative complexity of mine non-mining area structure Download PDFInfo
- Publication number
- CN107169616B CN107169616B CN201710599393.6A CN201710599393A CN107169616B CN 107169616 B CN107169616 B CN 107169616B CN 201710599393 A CN201710599393 A CN 201710599393A CN 107169616 B CN107169616 B CN 107169616B
- Authority
- CN
- China
- Prior art keywords
- relative
- area
- index data
- complexity
- coal seam
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Abstract
The invention discloses a relatively complex construction process of a mine non-mining areaRelative entropy prediction method of degree. The prediction method comprises the steps of obtaining evaluation index data of a mined area, calculating a relative entropy actual value S of the relative complexity of the mined area structure according to three structure type index data, namely a flexure type, a fault number and a fault scale, and substituting the relative entropy actual value S of the relative complexity of the mined area structure and non-structure type index data into a formula S ═ beta0+β1x1+…+βpxpTo calculate a relation parameter beta p and a relation constant beta0(ii) a After the relation parameter beta p is verified, the relation parameter beta p and the relation constant beta are added0And non-structural type index data X of a non-mining area of a mine fieldpSubstituting S ═ β0+β1X1+…+βpXpCalculating the relative entropy S' of the relative construction complexity of the unexplored area of the mine field; and finally, judging the complexity level corresponding to the relative entropy S' of the relative complexity of the structure of the non-mining area of the mining area well field by utilizing a relation comparison table of the relative complexity and the relative entropy of the structure. The prediction method is simple to operate and has a good prediction effect.
Description
Technical Field
The invention belongs to the technical field of mine exploitation, and particularly relates to a relative entropy prediction method for the relative complexity of a mine unexplored area structure.
Background
The mechanical mining process of the coal mine needs to evaluate and predict mining geological conditions firstly, the mine geological structure is the most main factor in the mining geological conditions, the mine geological structure can influence the mining deployment and the safe mining of the coal mine, the complexity of the mine structure not only determines the selection of the design and the development mode of the coal mine well, but also has important influence on water inrush of the coal mine, surrounding rock deformation and coal and gas outburst. Therefore, the evaluation and prediction of the complexity of the mine construction are important components of modern mine construction.
The quantitative evaluation of the mine structure is to select proper quantitative indexes on the basis of the research of the structure rule, adopt a proper mathematical model and take the computer technology as a means to quantitatively judge the structure complexity of different sections of the mine and mark out corresponding structure grades. At present, a unique method for quantitative evaluation and research of geological structure complexity, such as an isocratic block index method, an artificial intelligence method, a fuzzy mathematics method, a geometric fractal method, a clustering method and the like, has certain limitations and is not perfect. The formation of geological formations in coal mines, typically faults and buckling, is well documented in the areas exposed by the mining of the coal mine, but not in the unexplored areas. How to predict the situation in the non-mining area according to the evaluation result of the geological structure complexity in the mining area still needs to carry out related research work.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a relative entropy prediction method for the relative complexity of the mine unexplored area structure. The prediction method can calculate the relative entropy of the relative complexity of the structure of the mined area of the same mine according to the relative entropy of the relative complexity of the structure of the mined area of the mine, so that the relative complexity of the structure of the mined area of the mine can be judged.
In order to achieve the purpose, the invention adopts the technical scheme that: the relative entropy prediction method for the relative complexity of the mine unexplored area structure is characterized by comprising the following steps of:
dividing a mined area of a mining area well field into a mined area latticed equal-performance block section consisting of a plurality of unit grids;
acquiring evaluation index data of the mined area of the mine field through the grid-shaped isocratic block of the mined area, wherein the evaluation index data comprises the fold type, the number of faults, the fault scale, the coal seam thickness abnormal index, the overlying sandstone thickness, the coal seam roof lithology, the coal seam floor elevation abnormal index and the coal seam floor contour number, the fold type, the number of faults and the fault scale are collectively called as construction type index data, and the coal seam thickness, the coal seam thickness abnormal index, the overlying sandstone thickness, the coal seam roof lithology, the coal seam floor elevation abnormal index and the coal seam floor contour number are collectively called as non-construction type index data;
thirdly, calculating a relative entropy actual value S of the relative complexity of the mined area of the mining area well field according to the construction type index data;
drawing contour lines on the equal-sex blocks according to the actual value S of the relative entropy of the relative complexity of the structure of the mined area of the mine field to form a relative complexity distinguishing diagram of the structure, and establishing a relation comparison table of the relative complexity of the structure and the relative entropy, wherein the relative complexity of the structure is divided into three levels of complexity, medium and simple;
step five, according to the relative entropy actual value S of the relative complexity of the structure of the mined area and the non-structure type index data x of the mined area and the mined areapObtaining the formula S ═ beta0+β1x1+…+βpxpThe relation parameter beta p and the relation constant beta in (1)0；
Step six, verifying a relation parameter beta p: because the relation parameter beta p is an estimated value, substituting the relation parameter beta p into a formula
Calculating the relative entropy of the relative complexity of the mined area structure of the mining area well field
Judging the calculation value of the relative entropy according to the relation comparison table of the relative complexity and the relative entropy
And whether the actual value S of the relative entropy belongs to the same construction relative complexity level, if so, adjusting the relation parameter beta p is not needed, and if not, adjusting the relation parameter beta p;
dividing the non-mining area of the mine field into non-mining area latticed equal-property blocks consisting of a plurality of cells;
step eight, acquiring non-structural type index data such as coal seam thickness, abnormal index of coal seam thickness, overlying sandstone thickness, lithology of a coal seam roof, elevation of a coal seam floor, abnormal index of elevation of the coal seam floor, number of contour lines of the coal seam floor and the like of the non-mined area of the mining area well field through the latticed isocratic block sections of the non-mined area;
ninthly, acquiring unstructured type index data X of the unexplored area of the mine fieldpSubstituting the formula S ═ β0+β1X1+…+βpXpCalculating the relative entropy S' of the relative complexity of the structure of the non-mining area of the mine field;
step ten, judging the complexity level corresponding to the relative entropy S' of the relative complexity of the structure of the non-mining area of the mining area well field by using a relation comparison table of the relative complexity of the structure and the relative entropy.
The relative entropy prediction method for the relative complexity of the mine non-mining area structure is characterized by specifically comprising the following steps of, when calculating the relative entropy actual value S of the relative complexity of the mine field and the mine field mined area structure in the third step:
301, making an entropy calculation assignment table of the construction type index data, wherein the entropy calculation assignment table assigns values to large faults, medium faults and small faults, and single-dip type flexure types, syncline flexure types, anticline flexure types and total fault numbers in fault scales according to weights;
step 302, making a statistical information table of evaluation index data of a mined area of a mining area well field, and assigning the structure type index data in the statistical information table according to the entropy value calculation assignment table;
step 303, passing a formula
Calculating a relative entropy actual value S of the relative complexity of the mined area structure of the mining area well field, wherein n is 5 to represent index numbers of five specific structure types, namely a fold type, the total number of faults, the number of large-scale faults, the number of medium-scale faults and the number of small-scale faultsAccording to the number of i, PiThe probability of occurrence of index data i of five specific structure types, namely fold type, total number of faults, number of large-scale faults, number of medium-scale faults and number of small-scale faults is specific
Wherein X is the assignment of the assignment table to the index data of the construction type according to the entropy value calculation, and i is the number corresponding to each construction type.
The relative entropy prediction method for the relative complexity of the mine unexplored area structure is characterized in that the unit grids are square, the size of each unit grid is 200m multiplied by 200m, and the direction of each unit grid is parallel to the advancing direction of a mine working face and perpendicular to the advancing direction of the mine working face.
The relative entropy prediction method for the relative complexity of the mine unexplored area structure is characterized by comprising the following steps of:
step S1, determining coordinates at the cell center point of the equal block, specifically: firstly, connecting the center points of each unit of the equal block segments by line segments by means of computer aided drafting (AutoCAD); then selecting the line segment, obtaining coordinates of the center point of the cell by adopting a query command, and storing the coordinates as a text document;
step S2, obtaining each evaluation index data contour, specifically: drawing each evaluation index data contour line according to the range of the horizontal coordinate and the vertical coordinate of the target area by using a kriging interpolation method according to the data disclosed by drilling in the target area by using geographic data gridding drawing software Surfer to obtain a corresponding grid file, and storing the grid file as a text document;
step S3, acquiring each evaluation index data at the center point of the cell of the equal sex block, specifically: and opening the text document obtained by the geographic data gridding drawing software Surfer in the step S2 by using a spatial analysis function in the geographic information system MapGIS by using the geographic information system MapGIS, saving the text document as a grid file, opening the grid file by using an image processing function in the geographic information system MapGIS, and selecting the text document data obtained by the AutoCAD software in the step S1 so as to obtain each evaluation index data at the central point of the cell.
Compared with the prior art, the invention has the following advantages:
1. the prediction method of the invention is based on the evaluation of the relative structural complexity of the area which is already mined and exposed in the mine field, and combined with the statistical evaluation index data in the equal block of the field division, the relative structural complexity of the area which is not mined and exposed in the same mine is predicted.
2. By adopting the prediction method, different mining modes and mining parameters can be adopted according to geological conditions disclosed by unexplored areas of the mine field in the mining area and aiming at different relative complexity, and a certain technical basis is provided for safe and efficient mining of the mine.
3. The prediction method of the invention calculates the relative entropy actual value S of the relative complexity of the structure of the mined area by acquiring the evaluation index data of the mined area and mainly according to the index data of three structure types, namely the fold type, the fault number and the fault scale, and then substitutes the relative entropy actual value S of the relative complexity of the structure of the mined area and the index data of the non-structure type into a formula S ═ beta0+β1x1+…+βpxpTo calculate a relation parameter beta p and a relation constant beta0(ii) a After the relation parameter beta p is verified, the relation parameter beta p and the relation constant beta are added0And non-structural type index data X of a non-mining area of a mine fieldpSubstituting S ═ β0+β1X1+…+βpXpCalculating the relative entropy S' of the relative complexity of the structure of the non-mining area of the mine field; and finally, judging the complexity level corresponding to the relative entropy S' of the relative complexity of the structure of the non-mining area of the mining area well field by utilizing a relation comparison table of the relative complexity and the relative entropy of the structure. The prediction method is simple to operate and has a good prediction effect.
The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is a schematic diagram of the grid-shaped equal block division of the mined area and the distribution of the structural type indexes in the mined area.
FIG. 3 is a schematic view of a contour line of coal seam thickness within a grid-like isotropic block of a mined area according to the present invention.
FIG. 4 is a schematic view of the thickness contour of overlying sandstone in a latticed equal block of a mined area according to the present invention.
FIG. 5 is a schematic diagram of the lithology distribution of the coal seam roof in a grid-like equal-performance block of a mined area according to the present invention.
FIG. 6 is a schematic view of the elevation contour of the floor of the coal seam in a grid-like equal-performance block of a mined area according to the present invention.
FIG. 7 is a graph of the relative complexity of the structure of the grid-like equal-quality block of the mined area.
Detailed Description
A method for predicting relative entropy of relative complexity of a mine unexplored area structure as shown in fig. 1 comprises the following steps:
dividing a mined area of a mining area well field into a mined area latticed equal-performance block section consisting of a plurality of unit grids;
specifically, a contour diagram of a coal seam floor of a certain well field is used as a base map to divide a grid-shaped equal block of a mined area, the directions of unit grids in the grid-shaped equal block of the mined area are determined to be a north-south direction and a east-west direction according to the advancing direction of a working face in the well field, the unit grids are square, and the size of the unit grids is 200m multiplied by 200 m;
acquiring evaluation index data of the mined area of the mine field through the grid-shaped isocratic block of the mined area, wherein the evaluation index data comprises the fold type, the number of faults, the fault scale, the coal seam thickness abnormal index, the overlying sandstone thickness, the coal seam roof lithology, the coal seam floor elevation abnormal index and the coal seam floor contour number, the fold type, the number of faults and the fault scale are collectively called as construction type index data, and the coal seam thickness, the coal seam thickness abnormal index, the overlying sandstone thickness, the coal seam roof lithology, the coal seam floor elevation abnormal index and the coal seam floor contour number are collectively called as non-construction type index data;
the fold type mainly reflects the complexity of a fold structure, and is distinguished by three types, namely anticline, syncline and monocline, the complexity of the fold structure reflected by the three types is reduced in sequence, the fold type is determined to be anticline or syncline according to the fold axis trace type in a cell during statistics, if no fold axis trace occurs in a certain cell, the cell is considered to be corresponding to the monocline, and the corresponding monocline is taken as a statistical value of the index;
the number of faults, the index mainly reflects the development density of the faults in a certain range, and the total number of the faults in the cells is used as a statistical value of the index during statistics;
the fault scale mainly reflects the development degree of fracture structures, the fault fall is used for distinguishing, the larger the fault fall is, the more the surface fracture structure develops, and the number of large faults (fall is more than 5m), the number of medium faults (fall is 1-5 m) and the number of small faults (fall is less than 1m) are respectively recorded as the statistical values of the index during statistics according to the condition of the fall of each fault in a cell;
the fault thickness mainly reflects the spatial distribution form of the geological structure to the coal seam thickness, the change of the coal seam thickness is controlled by various geological structures, and the coal seam thickness value of the central point in a unit cell is taken as a statistical value of the index during statistics;
the abnormal index of the coal seam thickness mainly reflects the change degree of the coal seam thickness relative to the average value, the development degree of the geological structure is evaluated by the index, and the ratio of the coal seam thickness of the central point in the unit cell to the average value of the coal seam thickness revealed by the drill hole is taken as the statistical value of the index during statistics;
the overlying sandstone thickness mainly reflects the difference of the condition change of the sedimentary medium, and the thickness of the rock stratum within 50cm above the coal seam roof is taken as a statistical range;
the lithology of the coal seam roof mainly reflects the difference of the condition change of the deposited medium, the lithology of the coal seam roof has a control effect on structural deformation in the coal seam, and according to the lithology of the coal seam roof disclosed by the drilling hole, the lithology corresponding to the central point in the unit cell is taken as a statistical value of the index during statistics;
the elevation of the coal seam floor mainly reflects the change of the spatial form of the coal seam floor, the fluctuation change of the coal seam floor is the result of the reconstruction of the geological structure, and the ratio of the elevation of the coal seam floor at the central point in a unit cell to the average value of the elevation of the coal seam floor revealed by a drill hole is taken as the statistical value of the index during statistics;
the abnormal index of the elevation of the coal seam floor mainly reflects the variability of the coal seam morphology, the development degree of a geological structure is evaluated by using the index, and the ratio of the elevation of the coal seam floor at the central point in a unit cell to the average value of the elevation of the coal seam bottom edge revealed by a drilling hole is taken as a statistical value of the index during statistics;
the number of contour lines of the coal seam floor mainly reflects the degree of urgency of the inclination of the coal seam floor, the denser the contour lines of the coal seam floor, the greater the inclination degree of the coal seam floor on the surface, namely the stronger the transformation effect of the geological structure, and the number of contour lines of the coal seam floor in a unit grid is taken as the statistical value of the index during statistics;
step three, calculating a relative entropy actual value S of the relative complexity of the mined area of the mine field according to the construction type index data, which specifically comprises the following steps:
301, making an entropy calculation assignment table of the construction type index data, wherein the entropy calculation assignment table assigns values to large faults, medium faults and small faults, and single-dip type flexure types, syncline flexure types, anticline flexure types and total fault numbers in fault scales according to weights;
in the implementation, the fault belongs to a large fault with the fall of more than 5m, the fault belongs to a medium fault with the fall of 1-5 m, and the fault belongs to a small fault with the fall of less than 5 m. When assigning, weights are assigned according to the sizes of faults and ruffles, and the assignment with large size is large, whereas the assignment with small size is small.
In this embodiment, the entropy value calculation assigned table is as follows:
table 1 entropy calculation assignment table
Step 302, making a statistical information table of evaluation index data of a mined area of a mining area well field, and assigning the structure type index data in the statistical information table according to the entropy value calculation assignment table;
in this embodiment, according to table 1, when the evaluation index data of the mined area and the mined field is quantitatively processed, the anticline is represented by 4, the syncline is represented by 3, the monocline is represented by 2, the information in the statistical data cells is shown in fig. 2, 3, 4, 5 and 6, the statistical result is shown in table 2,
TABLE 2 statistical information table of some cells of the grid-shaped equal block segments of the mined area
The coal seam thickness abnormality index and the coal seam floor elevation abnormality index are calculated according to the following formula:
coal seam thickness anomaly index being average coal seam thickness/coal seam thickness
Coal seam floor elevation abnormal index is equal to coal seam floor elevation/coal seam floor elevation average value
And according to the coal seam thickness and the coal seam floor elevation numerical values revealed by the drill holes in the well field, taking the arithmetic average value as the average value of the coal seam thickness and the average value of the coal seam floor elevation.
Step 303, passing a formula
Calculating a relative entropy actual value S of the relative complexity of the mined area structure of the mining area well field, wherein n is 5 to represent the number of index data i of five specific structure types including a fold type, the total number of faults, the number of large-scale faults, the number of medium-scale faults and the number of small-scale faults, and P isiIs of fold type, total number of fault and large sizeThe probability of occurrence of five specific structure type index data i, namely the number of fault, the number of medium-sized fault and the number of small-sized fault, is specific
Wherein X is the assignment of the assignment table to the index data of the construction type according to the entropy value calculation, and i is the number corresponding to each construction type.
In this embodiment, the specific calculation of the relative entropy actual value S of the relative complexity of the mined area and the shaft field is as follows:
drawing contour lines on the equal-sex blocks according to the relative entropy of the relative complexity of the structure of the mined area of the mine field to form a relative complexity distinguishing diagram of the structure, and establishing a relation comparison table of the relative complexity of the structure and the relative entropy, wherein the relative complexity of the structure is divided into three levels, namely complexity, medium and simple;
in this embodiment, the relationship between the relative complexity and the relative entropy is illustrated in the following table 3:
TABLE 3 construction of a relative complexity versus relative entropy
Step five, according to the relative entropy actual value S of the relative complexity of the structure of the mined area and the non-structure type index data x of the mined area and the mined areaiObtaining the formula S ═ beta0+β1x1+…+βpxpA relation parameter beta p and a relation constant beta0；
In this embodiment, the data (including the data in table 2) in the cells of the equal block segments divided by the mined area are counted according to the third step, the actual value S of the relative entropy is calculated, the actual value S of the relative entropy is used as a dependent variable, and the coal seam is processedThickness x1Coal seam thickness anomaly index x2Overlying sandstone thickness x3Lithology x of coal seam roof4Coal seam floor elevation x5And abnormal index x of elevation of coal seam floor6Number of contour lines x7As an independent variable, there is a correlation between a dependent variable and an independent variable, and the correlation is that the independent variable takes a fixed value and the dependent variable takes a value randomly within a certain error range. At the moment, the relation constant beta is estimated by adopting stepwise regression analysis based on the least square principle0Relation parameter betap。
The least squares principle, provided that a relationship constant β has been obtained0And a relation parameter betap(β1、β2、β3、β4、β5、β6、β7) According to the formula
Obtaining a calculated value of the relative entropy of a certain cell in the equal block
The numerical value is different from the actual relative entropy value S of the cell to obtain a calculated relative entropy value
The square of the difference with the actual value S of the relative entropy is Q, i.e.
The smaller Q is, the calculated value of the relative entropy of the cell is shown
The closer to the relative entropy actual value S, the minimum value of Q must exist, and the relationship constant beta is assumed at this time0And a relation parameter beta1、β2、β3、β4、β5、β6、β7I.e. the best estimate. According to the extreme value principle, Q is firstly respectively paired with beta0、β1、……、βpCalculating partial derivative, takingThe result is 0, an equation set is established, and then the equation set is solved to obtain a relation constant beta0And a relation parameter beta1、β2、β3、β4、β5、β6、β7The following are:
in the formula, i refers to the number of cells, and j refers to the number of independent variables.
For the cell statistics listed in table 2, specific calculations are given as follows:
Si＝0.8614
the stepwise regression analysis is to determine the formula S ═ beta0+β1x1+…+βpxpWhen considering each argument xpFor the importance degree of the formula, important variables which contribute a large amount to the formula are reserved, and secondary variables which contribute a small amount to the formula are removed.
Through a plurality of cells, the relation constant beta can be obtained0And a relation parameter betap(p＝1,2,…)；
Step six, verifying a relation parameter beta p: because the relation parameter beta p is an estimated value, substituting the relation parameter beta p into a formula
Calculating the relative entropy of the relative complexity of the mined area structure of the mining area well field
Judging the calculation value of the relative entropy according to the relation comparison table of the relative complexity and the relative entropy
And whether the actual value S of the relative entropy belongs to the same structure relative complexity degree, if so, adjusting the relation parameter beta p is not needed, and if not, adjusting the relation parameter beta p;
dividing the non-mining area of the mine field into non-mining area latticed equal-property blocks consisting of a plurality of cells;
specifically, a contour diagram of a coal bed floor of a certain well field is used as a base diagram to divide latticed equal blocks of an unexploited area, the directions of unit grids in the latticed equal blocks of the unexploited area are determined to be the north-south direction and the east-west direction according to the advancing direction of a working face in the well field, the unit grids are square, and the size of the unit grids is 200m multiplied by 200 m;
step eight, acquiring non-structural type index data such as coal seam thickness, abnormal index of coal seam thickness, overlying sandstone thickness, lithology of a coal seam roof, elevation of a coal seam floor, abnormal index of elevation of the coal seam floor, number of contour lines of the coal seam floor and the like of the non-mined area of the mining area well field through the latticed isocratic block sections of the non-mined area;
ninthly, acquiring unstructured type index data X of the unexplored area of the mine fieldpSubstituting the formula S ═ β0+β1X1+…+βpXpCalculating the relative entropy S' of the relative complexity of the structure of the non-mining area of the mine field;
step ten, judging the complexity level corresponding to the relative entropy S' of the relative complexity of the structure of the non-mining area of the mining area well field by using a relation comparison table of the relative complexity of the structure and the relative entropy.
In this embodiment, the following steps are adopted when obtaining the evaluation index data of the meshed isogenous block of the mined area and the evaluation index data of the meshed isogenous block of the non-mined area of the mining area well field:
step S1, determining coordinates at the cell center point of the equal block, specifically: firstly, connecting the center points of each unit of the equal block segments by line segments by means of computer aided drafting (AutoCAD); then selecting the line segment, obtaining coordinates of the center point of the cell by adopting a query command, and storing the coordinates as a text document;
step S2, obtaining each evaluation index data contour, specifically: drawing each evaluation index data contour line according to the range of the horizontal coordinate and the vertical coordinate of the target area by using a kriging interpolation method according to the data disclosed by drilling in the target area by using geographic data gridding drawing software Surfer to obtain a corresponding grid file, and storing the grid file as a text document;
step S3, acquiring each evaluation index data at the center point of the cell of the equal sex block, specifically: and opening the text document obtained by the geographic data gridding drawing software Surfer in the step S2 by using a spatial analysis function in the geographic information system MapGIS by using the geographic information system MapGIS, saving the text document as a grid file, opening the grid file by using an image processing function in the geographic information system MapGIS, and selecting the text document data obtained by the AutoCAD software in the step S1 so as to obtain each evaluation index data at the central point of the cell.
In this embodiment, the relative entropy prediction method for the relative complexity of the structure of the unexploded mine area predicts the relative complexity of the structure of the unexploded area in the same mine by combining statistical evaluation index data in the equal block segments divided by the mine field on the basis of evaluating the relative complexity of the structure of the area already mined and exposed by the mine field. By adopting the prediction method, different mining modes and mining parameters can be adopted according to geological conditions disclosed by unexplored areas of the mine field in the mining area and aiming at different relative complexity, and a certain technical basis is provided for safe and efficient mining of the mine.
In this embodiment, the prediction method calculates a relative entropy actual value S of the relative complexity of the structure of the mined area by acquiring evaluation index data of the mined area, mainly based on three structure type index data, namely a flexure type, a fault number and a fault scale, and substitutes the relative entropy actual value S of the relative complexity of the structure of the mined area and non-structure type index data into a formula S ═ β0+β1x1+…+βpxpTo calculate a relation parameter beta p and a relation constant beta0(ii) a After the relation parameter beta p is verified, the relation parameter beta p and the relation constant beta are added0And non-mined areas of the mine fieldStructure type index data XpSubstituting S ═ β0+β1X1+…+βpXpCalculating the relative entropy S' of the relative complexity of the structure of the non-mining area of the mine field; and finally, judging the complexity level corresponding to the relative entropy S' of the relative complexity of the structure of the non-mining area of the mining area well field by utilizing a relation comparison table of the relative complexity and the relative entropy of the structure. The prediction method is simple to operate and has a good prediction effect.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and all simple modifications, changes and equivalent structural changes made to the above embodiment according to the technical spirit of the present invention still fall within the protection scope of the technical solution of the present invention.
Claims (3)
1. The relative entropy prediction method for the relative complexity of the mine unexplored area structure is characterized by comprising the following steps of:
dividing a mined area of a mining area well field into a mined area latticed equal-performance block section consisting of a plurality of unit grids;
acquiring evaluation index data of the mined area of the mine field through the grid-shaped isocratic block of the mined area, wherein the evaluation index data comprises the fold type, the number of faults, the fault scale, the coal seam thickness abnormal index, the overlying sandstone thickness, the coal seam roof lithology, the coal seam floor elevation abnormal index and the coal seam floor contour number, the fold type, the number of faults and the fault scale are collectively called as construction type index data, and the coal seam thickness, the coal seam thickness abnormal index, the overlying sandstone thickness, the coal seam roof lithology, the coal seam floor elevation abnormal index and the coal seam floor contour number are collectively called as non-construction type index data;
thirdly, calculating a relative entropy actual value S of the relative complexity of the mined area of the mining area well field according to the construction type index data;
drawing isolines on the equal-performance blocks according to the relative entropy actual value S of the relative complexity of the mined area of the mining area well field to form a relative complexity partition diagram, and establishing a relation comparison table of the relative complexity and the relative entropy of the construction, wherein the relative complexity of the construction is divided into three levels, namely complexity, medium and simple;
step five, according to the relative entropy actual value S of the relative complexity of the structure of the mined area and the non-structure type index data x of the mined area and the mined areapObtaining the formula S ═ beta0+β1x1+...+βpxpThe relation parameter beta p and the relation constant beta in (1)0；
Step six, verifying a relation parameter beta p: because the relation parameter beta p is an estimated value, substituting the relation parameter beta p into a formula
Calculating the relative entropy of the relative complexity of the mined area structure of the mining area well field
Judging the calculation value of the relative entropy according to the relation comparison table of the relative complexity and the relative entropy
And whether the actual value S of the relative entropy belongs to the same construction relative complexity level, if so, adjusting the relation parameter beta p is not needed, and if not, adjusting the relation parameter beta p;
dividing the non-mining area of the mine field into non-mining area latticed equal-property blocks consisting of a plurality of cells;
step eight, acquiring the coal seam thickness, the coal seam thickness abnormal index, the overlying sandstone thickness, the coal seam roof lithology, the coal seam floor elevation abnormal index and the coal seam floor contour line number non-structural type index data of the mined area of the mined well field through the latticed equal-sex block section of the mined area;
ninthly, acquiring unstructured type index data X of the unexplored area of the mine fieldpSubstitute formula S'＝β0+β1X1+...+βpXpCalculating the relative entropy S' of the relative complexity of the structure of the non-mining area of the mine field;
step ten, judging the complexity level corresponding to the relative entropy S' of the relative complexity of the structure of the non-mining area of the mining area well field by utilizing a relation comparison table of the relative complexity of the structure and the relative entropy;
when calculating the relative entropy actual value S of the relative complexity of the mined area structure of the mining area well field in the third step, the method specifically comprises the following steps:
301, making an entropy calculation assignment table of the construction type index data, wherein the entropy calculation assignment table assigns values to large faults, medium faults and small faults, and single-dip type flexure types, syncline flexure types, anticline flexure types and total fault numbers in fault scales according to weights;
step 302, making a statistical information table of evaluation index data of a mined area of a mining area well field, and assigning the structure type index data in the statistical information table according to the entropy value calculation assignment table;
step 303, passing a formula
Calculating a relative entropy actual value S of the relative complexity of the mined area structure of the mining area well field, wherein n is 5 to represent the number of index data i of five specific structure types including a fold type, the total number of faults, the number of large-scale faults, the number of medium-scale faults and the number of small-scale faults, and P isiThe probability of occurrence of index data i of five specific structure types, namely fold type, total number of faults, number of large-scale faults, number of medium-scale faults and number of small-scale faults is specific
Wherein X is the assignment of the assignment table to the index data of the construction type according to the entropy value calculation, and i is the number corresponding to each construction type.
2. The method of claim 1, wherein the cells are square in shape and have a size of 200m x 200m, and are oriented in directions both parallel and perpendicular to the direction of advancement of the mine face.
3. The method for predicting the relative entropy of the relative complexity of the mine unexplored area structure according to claim 1, wherein the following steps are adopted when obtaining the evaluation index data of the meshed isosexual block of the mined area and the evaluation index data of the meshed isosexual block of the mine field unexplored area:
step S1, determining coordinates at the cell center point of the equal block, specifically: firstly, connecting the center points of each unit of the equal block segments by line segments by means of computer aided drafting (AutoCAD); then selecting the line segment, obtaining coordinates of the center point of the cell by adopting a query command, and storing the coordinates as a text document;
step S2, obtaining each evaluation index data contour, specifically: drawing each evaluation index data contour line according to the range of the horizontal coordinate and the vertical coordinate of the target area by using a kriging interpolation method according to the data disclosed by drilling in the target area by using geographic data gridding drawing software Surfer to obtain a corresponding grid file, and storing the grid file as a text document;
step S3, acquiring each evaluation index data at the center point of the cell of the equal sex block, specifically: and opening the text document obtained by the geographic data gridding drawing software Surfer in the step S2 by using a spatial analysis function in the geographic information system MapGIS by using the geographic information system MapGIS, saving the text document as a grid file, opening the grid file by using an image processing function in the geographic information system MapGIS, and selecting the text document data obtained by the AutoCAD software in the step S1 so as to obtain each evaluation index data at the central point of the cell.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710599393.6A CN107169616B (en) | 2017-07-21 | 2017-07-21 | Relative entropy prediction method for relative complexity of mine non-mining area structure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710599393.6A CN107169616B (en) | 2017-07-21 | 2017-07-21 | Relative entropy prediction method for relative complexity of mine non-mining area structure |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107169616A CN107169616A (en) | 2017-09-15 |
| CN107169616B true CN107169616B (en) | 2020-11-13 |
Family
ID=59818039
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710599393.6A Expired - Fee Related CN107169616B (en) | 2017-07-21 | 2017-07-21 | Relative entropy prediction method for relative complexity of mine non-mining area structure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107169616B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109409042B (en) * | 2018-08-23 | 2021-04-20 | 顺丰科技有限公司 | User authority distribution abnormity detection system, method, equipment and storage medium |
| CN110633557B (en) * | 2019-10-30 | 2023-04-14 | 太原理工大学 | Identification method for favorable area of coal bed gas structure |
| CN111582536B (en) * | 2020-03-23 | 2024-01-30 | 中国矿业大学(北京) | Feature learning-based hidden fault prediction method, device, equipment and medium |
| CN112380709A (en) * | 2020-11-18 | 2021-02-19 | 中国矿业大学 | Fold structure complexity identification method based on coal seam spreading characteristics |
| CN113821917A (en) * | 2021-08-31 | 2021-12-21 | 国能榆林能源有限责任公司 | Method for estimating shape of mined area top and bottom plates in fully mechanized mining equipment operation process |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101661114B (en) * | 2009-06-12 | 2010-10-13 | 中国矿业大学(北京) | Prediction method of small-scale structures in front of tunneled mine coal-shift based on ANN |
| CN101881182B (en) * | 2010-07-14 | 2011-10-26 | 中国矿业大学(北京) | Simulation and forecast evaluation method for delaying water inrush by mine fracture structure |
| CN102508941B (en) * | 2011-09-30 | 2014-10-29 | 中煤科工集团重庆研究院有限公司 | Method for accurate detection of geologic structures and division of coal and gas outburst areas |
| AU2012241140B2 (en) * | 2011-10-11 | 2016-01-21 | BELLAIRS, Jennifer Annette | A System and Method for Generating a Quantitative Factor of Safety and a Visual Representation Thereof |
| CN104200284B (en) * | 2014-09-01 | 2019-09-06 | 山东科技大学 | Mine fault tectonic prediction technique |
| CN104268647A (en) * | 2014-09-28 | 2015-01-07 | 山东科技大学 | Small fault strike extension length detecting method |
| CN104850897A (en) * | 2015-02-25 | 2015-08-19 | 中国矿业大学 | Prediction method for coal and gas outburst based on seismic information |
| CN104881571A (en) * | 2015-04-29 | 2015-09-02 | 西安科技大学 | Method for evaluating geological environment bearing capability of coal mine zone |
-
2017
- 2017-07-21 CN CN201710599393.6A patent/CN107169616B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN107169616A (en) | 2017-09-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107169616B (en) | Relative entropy prediction method for relative complexity of mine non-mining area structure | |
| RU2669948C2 (en) | Multistage oil field design optimisation under uncertainty | |
| CN110318744B (en) | Method for predicting shale gas resources | |
| WO2017058702A1 (en) | Ground support design tool | |
| EP2795527A1 (en) | Systems and methods for selecting facies model realizations | |
| EP2811108A1 (en) | Method for Assesing Production Strategy Plans | |
| CN101278301A (en) | Method of analyzing oil and gas production project description | |
| CN115293316A (en) | Prediction method for deep-buried thick coal seam water flowing fractured zone based on SSA-ELMAN | |
| AU2021103225A4 (en) | Embankment dam analysis to reduce flood impact using sensor modelling and gis techniques | |
| CN113033954B (en) | Intelligent decision-making method for coordinated development of coal and coalbed methane | |
| WO2012114146A2 (en) | Computerized method for the estimation of a value for at least a parameter of a hydrocarbon-producing region, for planning the operation and operating the region | |
| CN109640333B (en) | Underwater wireless sensor network coverage vulnerability remediation algorithm based on cluster division | |
| Wang et al. | Dynamic building energy demand modelling at urban scale for the case of Switzerland | |
| CN108595792A (en) | A kind of fuel tanker fuel measurement sensor placement optimization method | |
| CN113742962B (en) | Shale reservoir attribute three-dimensional modeling method based on horizontal well | |
| CN113221399B (en) | Decision method for dividing scheme of railway tunnel fabricated lining components | |
| CN116050877A (en) | Metal mine aquifer water enrichment evaluation method and evaluation system | |
| Yu et al. | Climatic formations: Evolutionary dynamics of urban morphologies | |
| CN114297865A (en) | Active utilization method and device for space stress field of three-dimensional well pattern | |
| CN109472867B (en) | Method for quantifying influence range of drilling position information | |
| Osterloh et al. | Probabilistic forecasting and model validation for the first-eocene large-scale pilot Steamflood, Partitioned Zone, Saudi Arabia and Kuwait | |
| CN113761699A (en) | Crack model construction method and device, electronic equipment and medium | |
| CN112734076A (en) | Large-scale positioning prediction method for hard rock type uranium ore resources | |
| CN114764546A (en) | Method for predicting residual oil of fault block oil reservoir based on generation of antagonistic neural network | |
| CN115576007B (en) | Semi-deterministic fracture modeling method and system based on disorder matrix |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Du Rongjun Inventor after: Xia Yucheng Inventor after: Wang Sherong Inventor after: Wei Zhaoxiang Inventor before: Du Rongjun Inventor before: Xia Yucheng |
|
| CB03 | Change of inventor or designer information | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180306 Address after: Yanta Road 710054 Shaanxi city of Xi'an province No. 58 Applicant after: XI'AN University OF SCIENCE AND TECHNOLOGY Applicant after: SHAANXI COAL MINING HANCHENG MINING Co.,Ltd. Address before: Yanta Road 710054 middle Shaanxi province Xi'an City No. 58 Applicant before: Xi'an University of Science and Technology |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201113 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |