A coal seam gas content dynamic prediction method based on Bayesian optimization XGBoost

Abstract

The application discloses a coal seam gas content dynamic prediction method based on Bayesian optimization XGBoost, and belongs to the field of coal mine gas prevention and control, and comprises the following steps: collecting a sample set U composed of influence factors of coal seam gas occurrence, pre-processing the sample set U, obtaining a sample set, and proportionally dividing the sample set to obtain a training set, a verification set and a test set; the training set is divided into multiple sample sets; an XGBoost model is constructed, the nth sample set is taken as input, the actual coal seam gas content is taken as output, the XGBoost model is trained, and the nth coal seam gas content prediction model is obtained; the influence factors of coal seam gas occurrence are calculated through a three-dimensional geological model, the geological structure features and the coal seam roof and floor lithological features are comprehensively considered, and the hyperparameters of the XGBoost model are adaptively searched through a Bayesian optimization algorithm, so that the technical problem that the real-time prediction of the gas occurrence state under complex geological conditions cannot be realized by the traditional gas content prediction method which depends on static geological parameters is solved.

Description

Technical Field

[0001] This invention belongs to the field of coal mine gas control and relates to a dynamic construction method for coal seam gas content attribute models, specifically a dynamic prediction method for coal seam gas content based on Bayesian optimization XGBoost. Background Technology

[0002] In underground coal mining, accurate prediction of coal seam gas content is one of the core technical challenges in preventing coal and gas outbursts and ensuring safe mine production. With increasing mining depth and intensity, the occurrence of coal seam gas is significantly influenced by the combined effects of geological structure distribution and environmental parameters. Traditional gas content prediction models largely rely on static geological parameters, resulting in the following drawbacks: First, the rigidity of model hyperparameters leads to insufficient dynamic adaptability to gas emission characteristics under complex mining conditions; second, the nonlinear correlations between multi-source heterogeneous data are not fully integrated, resulting in insufficient generalization ability of the prediction model; and third, obtaining relevant geological parameters is difficult, leading to long update cycles for the prediction model.

[0003] Therefore, there is an urgent need to construct a coal seam gas content prediction method that integrates multi-source dynamic data driving and has adaptive parameter optimization capabilities, so as to realize real-time modeling of gas occurrence state under complex geological conditions. Summary of the Invention

[0004] To address the shortcomings of existing technologies, the present invention aims to provide a dynamic prediction method for coal seam gas content based on Bayesian optimization XGBoost, thereby solving the technical problem that traditional gas content prediction methods rely heavily on static geological parameters and cannot achieve real-time prediction of gas occurrence under complex geological conditions.

[0005] To solve the above-mentioned technical problems, the present invention adopts the following technical solution:

[0006] A dynamic prediction method for coal seam gas content based on Bayesian optimization XGBoost, using a three-dimensional geological model including a fault plane model and a stratigraphic model, includes the following steps:

[0007] Step 1: Collect a sample set U consisting of influencing factors of coal seam gas occurrence, and preprocess it to obtain the sample set. The dataset is then divided proportionally to obtain a training set, a validation set, and a test set; the training set is further divided into multiple sample sets. ;

[0008] The influencing factors of coal seam gas occurrence include geological structural features and the characteristics of the coal seam roof and floor.

[0009] The geological structural features include fault distribution and fold distribution;

[0010] The characteristics of the top and bottom plates of the coal seam include the thickness of the coal seam, the thickness of the overlying bedrock, the burial depth of the coal seam, the sand-mudstone ratio of the coal seam to the top interface of the Yan'an Formation, the mudstone thickness of the top and bottom plates of the coal seam, and the mudstone thickness within 20m of the top plate of the coal seam.

[0011] Step 2: Construct the XGBoost model using the nth sample set. Using the actual coal seam gas content as the input and the actual coal seam gas content as the output, the XGBoost model is trained to obtain the nth coal seam gas content prediction model; among them, the Bayesian method is used for optimization.

[0012] Step 3: Input the validation set into the nth coal seam gas content prediction model to obtain the predicted coal seam gas content. Use K-cross-validation to verify the error rate between the predicted coal seam gas content and the actual coal seam gas content corresponding to the validation set. Input the test set into the nth coal seam gas content prediction model to obtain the predicted coal seam gas content. Calculate the coefficient of determination using the predicted coal seam gas content and the actual coal seam gas content corresponding to the test set. ;

[0013] Step 4: Select the coefficient of determination with the smallest error rate. The coal seam gas content prediction model that is closest to 1 is used as the final coal seam gas content prediction model.

[0014] Step 5: Grid the coal seam space to obtain the spatial location set of the coal seam. ;

[0015] Step 6: Use a three-dimensional geological model to calculate the set of spatial locations. Spatial position Factors influencing coalbed methane occurrence;

[0016] The influencing factors of coal seam gas occurrence include geological structural features and the characteristics of the coal seam roof and floor.

[0017] The geological structural features include fault distribution and fold distribution;

[0018] The characteristics of the top and bottom plates of the coal seam include the thickness of the coal seam, the thickness of the overlying bedrock, the burial depth of the coal seam, the sand-mudstone ratio of the coal seam to the top interface of the Yan'an Formation, the mudstone thickness of the top and bottom plates of the coal seam, and the mudstone thickness within 20m of the top plate of the coal seam.

[0019] Step 7: Combine the spatial location set obtained in Step 6. Spatial position The influencing factors of coal seam gas occurrence are input into the final coal seam gas content prediction model to obtain spatial location. The content of coal seam gas;

[0020] Step 8: Measure the properties of the coal and rock strata in real time during the mining process, update the modeling parameters of the three-dimensional geological model based on the properties of the coal and rock strata, and return to Step 6.

[0021] This invention also includes the following technical features:

[0022] Between steps seven and eight, the method of harmonic averaging is also used to determine the spatial location. The coal seam gas content was processed to obtain a coal seam gas occurrence attribute model.

[0023] Step six specifically includes the following steps:

[0024] Step 6.1: Calculate the coal seam dip angle according to formula (6-1). And use it as a folding parameter;

[0025] (6-1)

[0026] in:

[0027] Indicates spatial location The tangent vector at the point;

[0028] Indicates spatial location The magnitude of the normal vector at that point;

[0029] Step 6.2: Use formula (6-2) to calculate the coal seam thickness, overlying bedrock thickness, coal seam burial depth, mudstone thickness of the top and bottom plates of the coal seam, sandstone thickness of the top and bottom plates of the coal seam, and mudstone thickness within 20m of the top plate of the coal seam.

[0030] (6-2)

[0031] in:

[0032] Indicates geological layers, The value range is 1 to 6, which respectively represent the coal seam, the overlying bedrock, the coal seam burial depth, the mudstone on the top and bottom of the coal seam, the sandstone on the top and bottom of the coal seam, and the mudstone within 20m of the top of the coal seam.

[0033] Indicates spatial location geological strata The thickness;

[0034] Indicates spatial location geological strata The corresponding equation for the top plate surface;

[0035] Indicates spatial location geological strata The corresponding equation for the base plate surface;

[0036] Step 6.3: Calculate the ratio of the thickness of mudstone at the top and bottom of the coal seam to the thickness of sandstone at the top and bottom of the coal seam obtained in Step 6.2, which gives the sandstone-mudstone ratio at the top interface of the coal seam to the Yan'an Formation.

[0037] Step 6.4: Calculate the fault parameters according to formula (6-3). This yields the fault distribution;

[0038] (6-3)

[0039] in:

[0040] Faults in coal seam space The difference in elevation;

[0041] Indicates spatial location Faults in the coal seam space The shortest distance between them.

[0042] In step five, the mesh generation method is the Delaunay tetrahedral mesh generation method.

[0043] In step one, the preprocessing specifically involves: using the maximum-minimum normalization method, mapping the influencing factors of gas occurrence in each coal seam to the [0, 1] interval through linear transformation to obtain the sample set. .

[0044] Compared with the prior art, the beneficial technical effects of this invention are:

[0045] (I) In this invention, the influencing factors of coal seam gas occurrence are calculated by a three-dimensional geological model. Taking into account the geological structure characteristics and the lithological characteristics of the roof and floor of the coal seam, the hyperparameters of the XGBoost model are adaptively searched using the Bayesian optimization algorithm to achieve dynamic prediction of coal seam gas content. This solves the technical problem that traditional gas content prediction methods rely heavily on static geological parameters and cannot achieve real-time prediction of gas occurrence under complex geological conditions.

[0046] (II) Based on prediction, this invention constructs a coal seam gas content attribute model. At the same time, by continuously updating the modeling parameters of the three-dimensional geological model and iteratively updating the three-dimensional geological model, the influencing factors affecting coal seam gas occurrence are calculated, thus realizing the continuous updating of the coal seam gas content attribute model. Attached Figure Description

[0047] Figure 1 This is a schematic diagram of the coal seam space after grid partitioning according to the present invention;

[0048] Figure 2 This is a comparison chart of the predicted and actual values ​​of the final coal seam gas content prediction model of this invention;

[0049] Figure 3 This is a schematic diagram of the coal seam gas content attribute model of the present invention.

[0050] The specific content of the present invention will be further explained in detail below with reference to the embodiments. Detailed Implementation

[0051] It should be noted that, unless otherwise specified, all components in this invention are those known in the art.

[0052] The following are specific embodiments of the present invention. It should be noted that the present invention is not limited to the following specific embodiments. All equivalent modifications made based on the technical solutions of this application fall within the protection scope of the present invention.

[0053] This invention presents a dynamic prediction method for coal seam gas content based on Bayesian optimization XGBoost. Based on a three-dimensional geological model, which includes a fault plane model and a stratigraphic model, the method comprises the following steps:

[0054] Step 1: Collect a sample set U consisting of influencing factors of coal seam gas occurrence, and preprocess it to obtain the sample set. The dataset is then divided proportionally to obtain a training set, a validation set, and a test set; the training set is further divided into multiple sample sets. ;

[0055] Factors influencing coal seam gas occurrence include geological structural features and the characteristics of the coal seam roof and floor;

[0056] Geological structural features include the distribution of faults and folds;

[0057] The characteristics of the top and bottom plates of the coal seam include the thickness of the coal seam, the thickness of the overlying bedrock, the burial depth of the coal seam, the sand-mudstone ratio of the coal seam to the top interface of the Yan'an Formation, the mudstone thickness of the top and bottom plates of the coal seam, and the mudstone thickness within 20m of the top plate of the coal seam.

[0058] Step 2: Construct the XGBoost model using the nth sample set. Using the actual coal seam gas content as the input and the actual coal seam gas content as the output, the XGBoost model is trained to obtain the nth coal seam gas content prediction model; among them, the Bayesian method is used for optimization.

[0059] Step 3: Input the validation set into the nth coal seam gas content prediction model to obtain the predicted coal seam gas content. Use K-cross-validation to verify the error rate between the predicted coal seam gas content and the actual coal seam gas content corresponding to the validation set. Input the test set into the nth coal seam gas content prediction model to obtain the predicted coal seam gas content. Calculate the coefficient of determination using the predicted coal seam gas content and the actual coal seam gas content corresponding to the test set. ;

[0060] Step 4: Select the coefficient of determination with the smallest error rate. The coal seam gas content prediction model that is closest to 1 is used as the final coal seam gas content prediction model.

[0061] Step 5: Grid the coal seam space to obtain the spatial location set of the coal seam. ;

[0062] Step 6: Use a three-dimensional geological model to calculate the set of spatial locations. Spatial position Factors influencing coalbed methane occurrence;

[0063] Factors influencing coal seam gas occurrence include geological structural features and the characteristics of the coal seam roof and floor;

[0064] Geological structural features include the distribution of faults and folds;

[0065] The characteristics of the top and bottom plates of the coal seam include the thickness of the coal seam, the thickness of the overlying bedrock, the burial depth of the coal seam, the sand-mudstone ratio of the coal seam to the top interface of the Yan'an Formation, the mudstone thickness of the top and bottom plates of the coal seam, and the mudstone thickness within 20m of the top plate of the coal seam.

[0066] Step 7: Combine the spatial location set obtained in Step 6. Spatial position The influencing factors of coal seam gas occurrence are input into the final coal seam gas content prediction model to obtain spatial location. The content of coal seam gas;

[0067] Step 8: Measure the properties of the coal and rock strata in real time during the mining process, update the modeling parameters of the three-dimensional geological model based on the properties of the coal and rock strata, and return to Step 6.

[0068] In the above scheme, the influencing factors of coal seam gas occurrence are calculated by a three-dimensional geological model. Taking into account the geological structural characteristics and the lithological characteristics of the roof and floor of the coal seam, the hyperparameters of the XGBoost model are adaptively searched using a Bayesian optimization algorithm to achieve dynamic prediction of coal seam gas content. This solves the technical problem that traditional gas content prediction methods rely heavily on static geological parameters and cannot achieve real-time prediction of gas occurrence under complex geological conditions.

[0069] Ideally, the ratio of training set, validation set, and test set should be 80%, 10%, and 10%, respectively.

[0070] The training set is used for training and learning feature parameters; the validation set is mainly used to tune the hyperparameters of the prediction model during training to confirm the optimal model; the test set is used for independent evaluation of the final model performance to verify the model's performance on unseen data; and K-cross-validation is used to improve the robustness and reliability of the evaluation results.

[0071] In this scheme, the XGBoost model based on the gradient boosting decision tree framework is adopted, which effectively constructs the complex mapping relationship between variables and avoids the possibility of highly nonlinear relationships between variables.

[0072] The spatial location calculated in this scheme Table 1 shows the influencing factors of coal seam gas occurrence and the corresponding actual coal seam gas content; the obtained spatial locations The influencing factors of coal seam gas occurrence are input into the final coal seam gas content prediction model to obtain spatial location. The predicted value of coal seam gas content, in this spatial location The comparison chart of predicted and actual coal seam gas content is shown below. Figure 2 As shown.

[0073] Between steps seven and eight, the method of harmonic averaging is also used to determine the spatial location. The coal seam gas content was processed to obtain a coal seam gas occurrence attribute model.

[0074] In the above scheme, a coal seam gas content attribute model is constructed based on prediction. At the same time, by continuously updating the modeling parameters of the three-dimensional geological model and iteratively updating the three-dimensional geological model, the influencing factors affecting coal seam gas occurrence are calculated, thus realizing the continuous updating of the coal seam gas content attribute model.

[0075] The constructed coal seam gas content attribute model is as follows Figure 3 As shown.

[0076] Step six specifically includes the following steps:

[0077] Step 6.1: Calculate the coal seam dip angle according to formula (6-1). And use it as a folding parameter;

[0078] (6-1)

[0079] in:

[0080] Indicates spatial location The tangent vector at the point;

[0081] Indicates spatial location The magnitude of the normal vector at that point;

[0082] Step 6.2: Use formula (6-2) to calculate the coal seam thickness, overlying bedrock thickness, coal seam burial depth, mudstone thickness of the top and bottom plates of the coal seam, sandstone thickness of the top and bottom plates of the coal seam, and mudstone thickness within 20m of the top plate of the coal seam.

[0083] (6-2)

[0084] in:

[0085] Indicates geological layers, The value range is 1 to 6, which respectively represent the coal seam, the overlying bedrock, the coal seam burial depth, the mudstone on the top and bottom of the coal seam, the sandstone on the top and bottom of the coal seam, and the mudstone within 20m of the top of the coal seam.

[0086] Indicates spatial location geological strata The thickness;

[0087] Indicates spatial location geological strata The corresponding equation for the top plate surface;

[0088] Indicates spatial location geological strata The corresponding equation for the base plate surface;

[0089] Step 6.3: Calculate the ratio of the thickness of mudstone at the top and bottom of the coal seam to the thickness of sandstone at the top and bottom of the coal seam obtained in Step 6.2, which gives the sandstone-mudstone ratio at the top interface of the coal seam to the Yan'an Formation.

[0090] Step 6.4: Calculate the fault parameters according to formula (6-3). This yields the fault distribution;

[0091] (6-3)

[0092] in:

[0093] Faults in coal seam space The difference in elevation;

[0094] Indicates spatial location Faults in the coal seam space The shortest distance between them;

[0095] Among them, the fault plane model is used to calculate spatial location. Faults in the coal seam space The shortest distance between them; the stratigraphic model is used to derive the virtual borehole, and the spatial location is obtained by fitting the interface of each geological layer in the virtual borehole. geological strata The corresponding equations for the top and bottom surfaces.

[0096] In step five, the mesh generation method is the Delaunay tetrahedral mesh generation method.

[0097] In the above technical solution, the coal seam space is divided into tetrahedral meshes after meshing, as shown in the figure. Figure 3 As shown.

[0098] In step one, the preprocessing specifically involves: using the maximum-minimum normalization method, mapping the influencing factors of gas occurrence in each coal seam to the [0, 1] interval through linear transformation to obtain the sample set. .

[0099] In the above technical solution, dimensional differences of all quantitative characteristics are eliminated through preprocessing.

[0100] Table 1

[0101]

Claims

1. A dynamic prediction method for coal seam gas content based on Bayesian optimization XGBoost, based on a three-dimensional geological model, wherein the three-dimensional geological model includes a fault plane model and a stratigraphic model; characterized in that, Includes the following steps: Step 1: Collect a sample set U consisting of influencing factors of coal seam gas occurrence, and preprocess it to obtain the sample set. The dataset is then divided proportionally to obtain a training set, a validation set, and a test set; the training set is further divided into multiple sample sets. ; The influencing factors of coal seam gas occurrence include geological structural features and the characteristics of the coal seam roof and floor. The geological structural features include fault distribution and fold distribution; The characteristics of the top and bottom plates of the coal seam include the thickness of the coal seam, the thickness of the overlying bedrock, the burial depth of the coal seam, the sand-mudstone ratio of the coal seam to the top interface of the Yan'an Formation, the mudstone thickness of the top and bottom plates of the coal seam, and the mudstone thickness within 20m of the top plate of the coal seam. Step 2: Construct the XGBoost model using the nth sample set. Using the actual coal seam gas content as the input and the actual coal seam gas content as the output, the XGBoost model is trained to obtain the nth coal seam gas content prediction model; among them, the Bayesian method is used for optimization. Step 3: Input the validation set into the nth coal seam gas content prediction model to obtain the predicted coal seam gas content. Use K-cross-validation to verify the error rate between the predicted coal seam gas content and the actual coal seam gas content corresponding to the validation set. Input the test set into the nth coal seam gas content prediction model to obtain the predicted coal seam gas content. Calculate the coefficient of determination using the predicted coal seam gas content and the actual coal seam gas content corresponding to the test set. ; Step 4: Select the coefficient of determination with the smallest error rate. The coal seam gas content prediction model that is closest to 1 is used as the final coal seam gas content prediction model. Step 5: Grid the coal seam space to obtain the spatial location set of the coal seam. ; Step 6: Use a three-dimensional geological model to calculate the set of spatial locations. Spatial position Factors influencing coalbed methane occurrence; The influencing factors of coal seam gas occurrence include geological structural features and the characteristics of the coal seam roof and floor. The geological structural features include fault distribution and fold distribution; The characteristics of the top and bottom plates of the coal seam include the thickness of the coal seam, the thickness of the overlying bedrock, the burial depth of the coal seam, the sand-mudstone ratio of the coal seam to the top interface of the Yan'an Formation, the mudstone thickness of the top and bottom plates of the coal seam, and the mudstone thickness within 20m of the top plate of the coal seam. Step 7: Combine the spatial location set obtained in Step 6. Spatial position The influencing factors of coal seam gas occurrence are input into the final coal seam gas content prediction model to obtain spatial location. The content of coal seam gas; Step 8: Measure the properties of the coal and rock strata in real time during the mining process, update the modeling parameters of the three-dimensional geological model based on the properties of the coal and rock strata, and return to Step 6.

2. The dynamic prediction method for coal seam gas content based on Bayesian optimized XGBoost as described in claim 1, characterized in that, Between steps seven and eight, the method of harmonic averaging is also used to determine the spatial location. The coal seam gas content was processed to obtain a coal seam gas occurrence attribute model.

3. The dynamic prediction method for coal seam gas content based on Bayesian optimized XGBoost as described in claim 1, characterized in that, Step six specifically includes the following steps: Step 6.1: Calculate the coal seam dip angle according to formula (6-1). And use it as a folding parameter; （6-1） in: Indicates spatial location The tangent vector at the point; Indicates spatial location The magnitude of the normal vector at that point; Step 6.2: Use formula (6-2) to calculate the coal seam thickness, overlying bedrock thickness, coal seam burial depth, mudstone thickness of the top and bottom plates of the coal seam, sandstone thickness of the top and bottom plates of the coal seam, and mudstone thickness within 20m of the top plate of the coal seam. （6-2） in: Indicates geological layers, The value range is 1 to 6, which respectively represent the coal seam, the overlying bedrock, the coal seam burial depth, the mudstone on the top and bottom of the coal seam, the sandstone on the top and bottom of the coal seam, and the mudstone within 20m of the top of the coal seam. Indicates spatial location geological strata The thickness; Indicates spatial location geological strata The corresponding equation for the top plate surface; Indicates spatial location geological strata The corresponding equation for the base plate surface; Step 6.3: Calculate the ratio of the thickness of mudstone at the top and bottom of the coal seam to the thickness of sandstone at the top and bottom of the coal seam obtained in Step 6.2, which gives the sandstone-mudstone ratio at the top interface of the coal seam to the Yan'an Formation. Step 6.4: Calculate the fault parameters according to formula (6-3). This yields the fault distribution; （6-3） in: Faults in coal seam space The difference in elevation; Indicates spatial location Faults in the coal seam space The shortest distance between them.

4. The dynamic prediction method for coal seam gas content based on Bayesian optimized XGBoost as described in claim 1, characterized in that, In step five, the mesh generation method is the Delaunay tetrahedral mesh generation method.

5. The dynamic prediction method for coal seam gas content based on Bayesian optimized XGBoost as described in claim 1, characterized in that, In step one, the preprocessing specifically involves: using the maximum-minimum normalization method, mapping the influencing factors of gas occurrence in each coal seam to the [0, 1] interval through linear transformation to obtain the sample set. .

Images