CN116910886A - An impact energy early warning method and system in karst areas - Google Patents
An impact energy early warning method and system in karst areas Download PDFInfo
- Publication number
- CN116910886A CN116910886A CN202311078824.6A CN202311078824A CN116910886A CN 116910886 A CN116910886 A CN 116910886A CN 202311078824 A CN202311078824 A CN 202311078824A CN 116910886 A CN116910886 A CN 116910886A
- Authority
- CN
- China
- Prior art keywords
- early warning
- impact energy
- karst
- prediction
- mine
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 235000019994 cava Nutrition 0.000 claims description 11
- 239000002689 soil Substances 0.000 claims description 10
- 238000004088 simulation Methods 0.000 claims description 9
- 238000012544 monitoring process Methods 0.000 claims description 8
- 238000007637 random forest analysis Methods 0.000 claims description 8
- 239000002245 particle Substances 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 abstract description 3
- 230000000704 physical effect Effects 0.000 abstract description 3
- 238000007405 data analysis Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 239000003245 coal Substances 0.000 description 24
- 238000005065 mining Methods 0.000 description 16
- 238000003066 decision tree Methods 0.000 description 14
- 239000011435 rock Substances 0.000 description 13
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000011161 development Methods 0.000 description 6
- 230000003628 erosive effect Effects 0.000 description 6
- 238000005755 formation reaction Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004090 dissolution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000009527 percussion Methods 0.000 description 2
- 230000035939 shock Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/08—Construction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Computer Hardware Design (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Architecture (AREA)
- Educational Administration (AREA)
- Civil Engineering (AREA)
- Game Theory and Decision Science (AREA)
- Mathematical Optimization (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Computational Mathematics (AREA)
- Structural Engineering (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
技术领域Technical field
本申请涉及冲击能量预警领域,具体涉及一种岩溶地区的冲击能量预警方法及系统。This application relates to the field of impact energy early warning, and specifically relates to an impact energy early warning method and system in karst areas.
背景技术Background technique
岩溶地区是指地层中富含溶洞和溶蚀形态的地区,由于岩石层的特殊结构和矿物组成,这类地区容易发生各种地质灾害,如地面塌陷、地下水突然涌出、地震等。当前的地质灾害监测方法大多采用传统的地质勘探手段,缺乏及时、准确的冲击能量预警方法,为有效防范和减轻岩溶地区灾害带来困难。Karst areas refer to areas rich in caves and dissolved forms in the strata. Due to the special structure and mineral composition of rock layers, such areas are prone to various geological disasters, such as ground collapse, sudden outflow of groundwater, earthquakes, etc. Most of the current geological disaster monitoring methods use traditional geological exploration methods and lack timely and accurate impact energy early warning methods, which brings difficulties to effectively prevent and mitigate disasters in karst areas.
发明内容Contents of the invention
本申请通过模拟岩溶地区的地质结构和物理性质的变化,识别冲击能量的释放和积聚过程,从而提前预警岩溶地区的地质灾害。This application simulates changes in the geological structure and physical properties of karst areas and identifies the release and accumulation process of impact energy, thereby providing early warning of geological disasters in karst areas.
为实现上述目的,本申请提供了一种岩溶地区的冲击能量预警方法,流程包括:In order to achieve the above purpose, this application provides an impact energy early warning method in karst areas. The process includes:
采集待开采矿井的历史工程地质资料;Collect historical engineering geological data of mines to be mined;
基于所述历史工程地质资料,布置勘探点,拟建建筑基础型式;Based on the historical engineering geological data, arrange exploration points and propose building foundation types;
基于所述建筑基础型式,构建模拟矿井;Based on the building foundation type, construct a simulated mine;
对所述模拟矿井进行初步预测,得到初步预测结果;Make preliminary predictions on the simulated mine and obtain preliminary prediction results;
基于所述初步预测结果,构建预测模型,得到最终预测结果,基于所述最终预测结果进行冲击能量的预警。Based on the preliminary prediction results, a prediction model is constructed to obtain the final prediction results, and an early warning of impact energy is performed based on the final prediction results.
优选的,拟建所述建筑基础型式的方法包括:Preferably, the proposed method for constructing the building foundation type includes:
通过所述历史工程地质资料,获取当地治理岩溶、土洞和塌陷的历史数据;Through the historical engineering geological data, obtain historical data on local control of karst, soil caves and collapse;
基于所述历史数据,确定所述建筑基础型式的荷载大小。Based on the historical data, the load magnitude of the building foundation type is determined.
优选的,当矿井条件复杂时,每个独立基础均布置勘探点,采用一柱一桩的方式逐柱布置勘探点。Preferably, when the mine conditions are complex, exploration points are arranged on each independent foundation, and the exploration points are arranged column by column in a column-by-pile manner.
优选的,构建所述模拟矿井的方法包括:基于勘探点的勘探结果,建立细观颗粒流模型,构建所述模拟矿井。Preferably, the method of constructing the simulated mine includes: establishing a mesoscopic particle flow model based on the exploration results of the exploration points, and constructing the simulated mine.
优选的,得到所述预测结果的方法包括:利用经验类比法和微震监测法得到所述预测结果。Preferably, the method for obtaining the prediction result includes: using an empirical analogy method and a microseismic monitoring method to obtain the prediction result.
优选的,构建所述预测模型的方法包括:基于所述预测结果,利用随机森林模型构建所述预测模型。Preferably, the method of constructing the prediction model includes: using a random forest model to construct the prediction model based on the prediction results.
本申请还提供了一种岩溶地区的冲击能量预警系统,包括:采集模块、拟建模块、定标模块、预测模块和构建模块;This application also provides an impact energy early warning system for karst areas, including: a collection module, a simulation module, a calibration module, a prediction module and a building module;
所述采集模块用于采集待开采矿井的历史工程地质资料;The collection module is used to collect historical engineering geological data of the mine to be mined;
所述拟建模块用于基于所述历史工程地质资料,布置勘探点,拟建建筑基础型式;The simulation module is used to arrange exploration points and the foundation type of the proposed building based on the historical engineering geological data;
所述定标模块用于基于所述建筑基础型式,构建模拟矿井;The calibration module is used to construct a simulated mine based on the building foundation type;
所述预测模块用于对所述模拟矿井进行初步预测,得到初步预测结果;The prediction module is used to make preliminary predictions on the simulated mine and obtain preliminary prediction results;
所述构建模块用于基于所述初步预测结果,构建预测模型,得到最终预测结果,基于所述最终预测结果进行冲击能量的预警。The building module is used to construct a prediction model based on the preliminary prediction results, obtain the final prediction results, and provide early warning of impact energy based on the final prediction results.
优选的,岩溶地区的冲击能量预警系统,其特征在于,所述拟建模块的工作流程包括:Preferably, the impact energy early warning system in karst areas is characterized in that the workflow of the simulation module includes:
通过所述历史工程地质资料,了解当地治理岩溶、土洞和塌陷的历史数据;Through the historical engineering geological data mentioned above, we can understand the historical data of local control of karst, soil caves and collapse;
基于所述历史数据,确定所述建筑基础型式的荷载大小。Based on the historical data, the load magnitude of the building foundation type is determined.
与现有技术相比,本申请的有益效果如下:Compared with the existing technology, the beneficial effects of this application are as follows:
本申请通过模拟地质结构和物理性质的变化,全面掌握冲击能量的释放和积聚过程,准确判断地质灾害风险。通过数据分析和处理,将原本需要监测信息转化为可行的预警标准,及时发布预警信息,提高地质灾害的预警准确性。This application simulates changes in geological structure and physical properties to comprehensively grasp the release and accumulation process of impact energy and accurately determine geological disaster risks. Through data analysis and processing, the originally required monitoring information is transformed into feasible early warning standards, early warning information is released in a timely manner, and the accuracy of early warning of geological disasters is improved.
附图说明Description of the drawings
为了更清楚地说明本申请的技术方案,下面对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to explain the technical solutions of the present application more clearly, the drawings required to be used in the embodiments are briefly introduced below. Obviously, the drawings in the following description are only some embodiments of the present application. For ordinary people in the art, Technical personnel can also obtain other drawings based on these drawings without exerting creative labor.
图1为本申请的方法流程示意图;Figure 1 is a schematic flow chart of the method of this application;
图2为本申请的系统结构示意图。Figure 2 is a schematic diagram of the system structure of this application.
具体实施方式Detailed ways
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only some of the embodiments of the present application, rather than all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by those of ordinary skill in the art without creative efforts fall within the scope of protection of this application.
为使本申请的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本申请作进一步详细的说明。In order to make the above objects, features and advantages of the present application more obvious and understandable, the present application will be described in further detail below in conjunction with the accompanying drawings and specific implementation modes.
实施例一Embodiment 1
在地质学研究上,冲击能量的预警本质上可以归结为对于冲击地压的预测,所以,本实施例利用对于预测冲击地压的方式来对于冲击能量进行预警。In geological research, the early warning of impact energy can essentially be attributed to the prediction of impact ground pressure. Therefore, this embodiment uses the method of predicting impact pressure to provide early warning of impact energy.
如图1所示,为本实施例方法流程示意图,流程包括:As shown in Figure 1, it is a schematic flow chart of the method in this embodiment. The process includes:
S1.采集待开采矿井的历史工程地质资料。S1. Collect historical engineering geological data of the mine to be mined.
岩溶地区即所谓的喀斯特地貌,喀斯特地貌以溶蚀作用为主,还包括流水的冲蚀、潜蚀,以及坍陷等机械侵蚀过程。所以因其特殊性,不能直接采用常规的方式来布置勘察点,容易发生危险,所以需要先采集拟建场地周边的工程地质资料(即历史工程地质资料)。Karst areas are so-called karst landforms. Karst landforms are dominated by dissolution and also include mechanical erosion processes such as water erosion, submerged erosion, and collapse. Therefore, due to its particularity, the conventional method cannot be used to directly arrange survey points, which is prone to danger. Therefore, it is necessary to collect engineering geological data around the proposed site (i.e., historical engineering geological data) first.
S2.基于历史工程地质资料,布置勘探点,拟建建筑基础型式。S2. Based on historical engineering geological data, arrange exploration points and propose building foundation types.
岩溶地区勘察之前,应充分收集拟建场地周边的工程地质资料,获取当地治理岩溶、土洞和塌陷的历史数据。在岩溶地区勘察工作量布置应注意以下几方面的问题:Before surveying karst areas, engineering geological data around the proposed site should be fully collected, and historical data on local control of karst, soil caves and collapse should be obtained. When arranging the survey workload in karst areas, attention should be paid to the following issues:
拟建建筑基础型式Proposed building foundation type
为了解拟建建筑荷载大小,应针对拟采用的基础型式布置勘探点。当条件复杂时每个独立基础均应布置勘探点,当采用一柱一桩基础时应逐柱布置勘探孔。In order to understand the load size of the proposed building, exploration points should be arranged according to the proposed foundation type. When conditions are complex, exploration points should be arranged on each independent foundation. When one column and one pile foundation is used, exploration holes should be arranged column by column.
岩溶稳定性的影响Effects of karst stability
勘察工作的布置应方法多样化和过程动态化山于地下岩溶的发育有隐蔽性、多样性和空间位置变化大的特点,在进行工作量布置时应注意采用多种手段相结合的方法。The layout of survey work should be diversified in methods and the process should be dynamic. The development of mountainous underground karst has the characteristics of concealment, diversity and large changes in spatial location. When arranging the workload, attention should be paid to using a combination of multiple methods.
对于埋藏较浅的覆盖型及裸露型岩溶地基,可以适当采用井探、槽探等方法查明浅部岩溶洞隙的形态。而对于埋藏型岩溶地基,在勘察过程中技术人员应该及时取得岩土层信息,如果出现岩溶发育较严重的场地,应根据实际情况,增加岩土测试、静探、物探等手段,进一步查明松散裂隙及隐蔽土洞等的埋藏位置。同时勘探点深度也应根据实际情况适当增减,以满足基础设计和稳定性验算的需要。For shallowly buried covered and exposed karst foundations, methods such as well exploration and trench exploration can be appropriately used to identify the shape of shallow karst caves. For buried karst foundations, technicians should obtain rock and soil layer information in a timely manner during the survey process. If there is a site with serious karst development, geotechnical testing, static exploration, geophysical exploration and other means should be added according to the actual situation to further identify Buried locations in loose fissures and hidden soil caves. At the same time, the depth of the exploration point should also be appropriately increased or decreased according to the actual situation to meet the needs of foundation design and stability verification.
S3.基于建筑基础型式,构建模拟矿井。S3. Construct a simulated mine based on the building foundation type.
之后,根据勘探点得到的结果,采用平行粘结模型建立矿井周边标准岩样的单轴压缩颗粒流数值模型(模拟矿井),通过试错法得到细观力学参数,得到标定结果。在本实施例中,标定结果可分为:岩层冲击倾向性和采场支承应力带特征,通过对于这两个细观参数进行预测,来得到后续预测结果。Afterwards, based on the results obtained from the exploration points, a parallel bonding model was used to establish a uniaxial compression particle flow numerical model (simulated mine) of standard rock samples around the mine. The mesoscopic mechanical parameters were obtained through a trial and error method, and the calibration results were obtained. In this embodiment, the calibration results can be divided into: rock formation impact tendency and stope support stress zone characteristics. By predicting these two mesoscopic parameters, subsequent prediction results are obtained.
S4.对模拟矿井进行初步预测,得到初步预测结果。S4. Make preliminary predictions on the simulated mine and obtain preliminary prediction results.
采用经验类比法对岩层冲击倾向进行估测Estimating rock formation impact tendency using empirical analogy method
通过对以往经验教训做出规律性的总结,并用于指导本矿或相似条件的其它矿井冲击地压煤层安全开采的方法称为经验类比法。The method of making regular summaries of past experiences and lessons and using them to guide the safe mining of percussion coal seams in this mine or other mines with similar conditions is called the experience analogy method.
在本实施例中,应着重考虑下列因素:本矿和邻矿的冲击地压现状和发展趋势;本煤层或邻层、邻区已发生过冲击地压;煤层老顶为厚5m以上、单轴抗压强度大于70MPa的坚硬岩层;岛形或半岛形煤柱;支承压力影响区;上部或下部遗留煤柱或回采边界,煤层厚度或倾角突然变化区;褶曲或断裂构造带等。In this embodiment, the following factors should be considered: the current status and development trend of rockbursts in this mine and adjacent mines; rockbursts have occurred in this coal seam or adjacent layers and adjacent areas; the old roof of the coal seam is more than 5m thick and single. Hard rock formations with axial compressive strength greater than 70MPa; island-shaped or peninsula-shaped coal pillars; support pressure affected areas; upper or lower remaining coal pillars or mining boundaries, areas with sudden changes in coal seam thickness or dip angle; fold or fault structural zones, etc.
就整个矿井而言,通过地质条件和开采技术条件的分析,可以圈定冲击地压可能发生的区域。随着开采深度的延深,当煤岩体应力满足强度条件时就可能发生冲击地压。始发冲击的深度通常称临界深度。从始发深度起,冲击地压就可能在煤柱、煤层凸出部位和邻近煤柱的上下煤层区段发生,并随着开采水平的延深,冲击地压发生地点和范围也随之扩大。所有靠近采掘工作面的区域,煤层厚度和倾角突然变化的区域,以及地质构造带都可能成为发生冲击地压的危险区域。As for the entire mine, through the analysis of geological conditions and mining technology conditions, the area where rockburst may occur can be delineated. As the mining depth increases, rockburst may occur when the stress of the coal and rock mass meets the strength condition. The depth of the initial impact is usually called the critical depth. From the initial depth, rockbursts may occur in coal pillars, coal seam protrusions, and upper and lower coal seam sections adjacent to coal pillars. As the mining level increases, the location and scope of rockbursts also expand. . All areas close to the mining face, areas with sudden changes in coal seam thickness and inclination, and geological structural zones may become dangerous areas for rockbursts.
采用微震监测法来对采场支承应力带特征进行估测Using microseismic monitoring method to estimate the characteristics of stope support stress zone
利用短周期地震仪监测记录0.5级以上冲击发生的次数及冲击地压释放的能量。利用此趋势预测预报近期冲击地压发生的趋势及应力释放情况。在定位系统建成之前,采用现在的地震仪现行监测。通过上述流程,得到初步预测结果。Use short-period seismometers to monitor and record the number of shocks above magnitude 0.5 and the energy released by the impact ground pressure. Use this trend prediction to predict the trend of recent geoburst occurrence and stress release. Until the positioning system is built, current seismometers will be used for current monitoring. Through the above process, preliminary prediction results are obtained.
经过初步预测,结合喀斯特地貌的特殊性,上述的临界深度不应超过668m;工作面超前支承压力集中范围为5~35m,应力集中系数为2.5,但上方砾岩层的超前压力影响范围达120m。因此,工作面采动集中应力对工作面影响较为明显。分层开采时上分层工作面周期来压强度最大达510kn/m2,来压较为强烈。After preliminary prediction, combined with the particularity of karst landforms, the above-mentioned critical depth should not exceed 668m; the advanced support pressure concentration range of the working face is 5 to 35m, and the stress concentration coefficient is 2.5, but the impact range of the advanced pressure of the upper conglomerate layer reaches 120m. Therefore, the concentrated stress of mining on the working face has a more obvious impact on the working face. During layered mining, the cyclic pressure intensity of the upper layered working face reaches a maximum of 510kn/m 2 , and the pressure is relatively strong.
S5.基于初步预测结果,构建预测模型,得到最终预测结果,基于最终预测结果进行冲击能量的预警。S5. Based on the preliminary prediction results, build a prediction model, obtain the final prediction results, and provide early warning of impact energy based on the final prediction results.
采用随机森林模型,基于初步预测结果,构建预测模型。A random forest model is used to build a prediction model based on the preliminary prediction results.
首先,收集由上述初步预测结果组成样本集。之后,采用均值代替法对收集到的样本集进行清洗;在本实施例中,主要为清洗缺失值、清洗与原数据类型不一致的内容、清洗掉不需要的数据,清洗逻辑错误数据。First, collect a sample set composed of the above preliminary prediction results. After that, the mean substitution method is used to clean the collected sample set; in this embodiment, the main steps are to clean missing values, clean content inconsistent with the original data type, clean unnecessary data, and clean logical error data.
之后,从样本集中筛选出特征数据构成特征集,利用随机森林模型基于建立的特征集采用bootstrap的方法从样本集的M个样本中随机抽取m个样本作为一个子训练集构建一颗决策树;采用相同的方法同步构建T(T>1)棵决策树;在特征集的n个特征中随机选择p个特征作为节点分裂的子集,根据平方误差选择p个特征误差最小的1个特征作为节点分裂特征,保持节点分裂直到该决策树不可再分裂;本实施例中,3个特征为根节点和内容节点,每个预测结果为输出也即叶子节点。对T棵决策树进行分裂组成随机森林;基于初步预测结果对分裂后的每颗决策树随机在M个样本中进行训练得到每颗决策树所对应的预测结果;对每颗决策树所对应的预测结果求取均值得到最终预测结果。After that, feature data are selected from the sample set to form a feature set, and a random forest model is used based on the established feature set and the bootstrap method is used to randomly select m samples from the M samples in the sample set as a sub-training set to build a decision tree; Use the same method to simultaneously construct T (T>1) decision trees; randomly select p features from the n features in the feature set as a subset of node splits, and select the feature with the smallest error among the p features as the subset based on the square error. The node split feature keeps the node split until the decision tree cannot be split anymore; in this embodiment, the three features are the root node and the content node, and each prediction result is the output, which is a leaf node. Split T decision trees to form a random forest; based on the preliminary prediction results, each split decision tree is randomly trained on M samples to obtain the prediction results corresponding to each decision tree; the prediction results corresponding to each decision tree are The prediction results are averaged to obtain the final prediction result.
最后基于最终预测结果进行冲击能量的预警。Finally, an early warning of impact energy is carried out based on the final prediction results.
实施例二Embodiment 2
如图2所示,为本实施例的系统结构示意图,包括:采集模块、拟建模块、定标模块、预测模块和构建模块。其中,采集模块用于采集待开采矿井的历史工程地质资料;拟建模块用于基于历史工程地质资料,布置勘探点,拟建建筑基础型式;定标模块用于基于建筑基础型式,构建模拟矿井;预测模块用于对模拟矿井进行初步预测,得到初步预测结果;构建模块用于基于初步预测结果,构建预测模型,得到最终预测结果,基于最终预测结果进行冲击能量的预警。As shown in Figure 2, it is a schematic diagram of the system structure of this embodiment, including: an acquisition module, a simulation module, a calibration module, a prediction module and a construction module. Among them, the acquisition module is used to collect historical engineering geological data of the mines to be mined; the simulation module is used to arrange exploration points and propose building foundation types based on historical engineering geological data; the calibration module is used to construct simulated mines based on the building foundation types. ; The prediction module is used to make preliminary predictions for the simulated mine and obtain preliminary prediction results; the construction module is used to build a prediction model based on the preliminary prediction results, obtain the final prediction results, and provide early warning of impact energy based on the final prediction results.
下面将结合本实施例,详细说明本申请如何解决实际生活中的技术问题。How this application solves technical problems in real life will be described in detail below in conjunction with this embodiment.
首先,利用采集模块采集待开采矿井的历史工程地质资料。First, use the acquisition module to collect historical engineering geological data of the mine to be mined.
岩溶地区即所谓的喀斯特地貌,喀斯特地貌以溶蚀作用为主,还包括流水的冲蚀、潜蚀,以及坍陷等机械侵蚀过程。所以因其特殊性,不能直接采用常规的方式来布置勘察点,容易发生危险,所以需要先采集拟建场地周边的工程地质资料(即历史工程地质资料)。Karst areas are so-called karst landforms. Karst landforms are dominated by dissolution and also include mechanical erosion processes such as water erosion, submerged erosion, and collapse. Therefore, due to its particularity, the conventional method cannot be used to directly arrange survey points, which is prone to danger. Therefore, it is necessary to collect engineering geological data around the proposed site (i.e., historical engineering geological data) first.
之后,拟建模块基于历史工程地质资料,布置勘探点,拟建建筑基础型式。After that, the simulation module arranges exploration points and proposes building foundation types based on historical engineering geological data.
岩溶地区勘察之前,应充分收集拟建场地周边的工程地质资料,获取当地治理岩溶、土洞和塌陷的历史数据。在岩溶地区勘察工作量布置应注意以下几方面的问题:Before surveying karst areas, engineering geological data around the proposed site should be fully collected, and historical data on local control of karst, soil caves and collapse should be obtained. When arranging the survey workload in karst areas, attention should be paid to the following issues:
拟建建筑基础型式Proposed building foundation type
为了解拟建建筑荷载大小,应针对拟采用的基础型式布置勘探点。当条件复杂时每个独立基础均应布置勘探点,当采用一柱一桩基础时应逐柱布置勘探孔。In order to understand the load size of the proposed building, exploration points should be arranged according to the proposed foundation type. When conditions are complex, exploration points should be arranged on each independent foundation. When one column and one pile foundation is used, exploration holes should be arranged column by column.
岩溶稳定性的影响Effects of karst stability
勘察工作的布置应方法多样化和过程动态化山于地下岩溶的发育有隐蔽性、多样性和空间位置变化大的特点,在进行工作量布置时应注意采用多种手段相结合的方法。The layout of survey work should be diversified in methods and the process should be dynamic. The development of mountainous underground karst has the characteristics of concealment, diversity and large changes in spatial location. When arranging the workload, attention should be paid to using a combination of multiple methods.
对于埋藏较浅的覆盖型及裸露型岩溶地基,可以适当采用井探、槽探等方法查明浅部岩溶洞隙的形态。而对于埋藏型岩溶地基,在勘察过程中技术人员应该及时取得岩土层信息,如果出现岩溶发育较严重的场地,应根据实际情况,增加岩土测试、静探、物探等手段,进一步查明松散裂隙及隐蔽土洞等的埋藏位置。同时勘探点深度也应根据实际情况适当增减,以满足基础设计和稳定性验算的需要。For shallowly buried covered and exposed karst foundations, methods such as well exploration and trench exploration can be appropriately used to identify the shape of shallow karst caves. For buried karst foundations, technicians should obtain rock and soil layer information in a timely manner during the survey process. If there is a site with serious karst development, geotechnical testing, static exploration, geophysical exploration and other means should be added according to the actual situation to further identify Buried locations in loose fissures and hidden soil caves. At the same time, the depth of the exploration point should also be appropriately increased or decreased according to the actual situation to meet the needs of foundation design and stability verification.
定标模块基于建筑基础型式,构建模拟矿井。The calibration module builds a simulated mine based on the basic building type.
之后,根据勘探点得到的结果,采用平行粘结模型建立矿井周边标准岩样的单轴压缩颗粒流数值模型(模拟矿井),通过试错法得到细观力学参数,得到标定结果。在本实施例中,标定结果可分为:岩层冲击倾向性和采场支承应力带特征,通过对于这两个细观参数进行预测,来得到后续预测结果。Afterwards, based on the results obtained from the exploration points, a parallel bonding model was used to establish a uniaxial compression particle flow numerical model (simulated mine) of standard rock samples around the mine. The mesoscopic mechanical parameters were obtained through a trial and error method, and the calibration results were obtained. In this embodiment, the calibration results can be divided into: rock formation impact tendency and stope support stress zone characteristics. By predicting these two mesoscopic parameters, subsequent prediction results are obtained.
预测模块对模拟矿井进行初步预测,得到初步预测结果。The prediction module makes preliminary predictions for the simulated mine and obtains preliminary prediction results.
采用经验类比法对岩层冲击倾向进行估测Estimating rock formation impact tendency using empirical analogy method
通过对以往经验教训做出规律性的总结,并用于指导本矿或相似条件的其它矿井冲击地压煤层安全开采的方法称为经验类比法。The method of making regular summaries of past experiences and lessons and using them to guide the safe mining of percussion coal seams in this mine or other mines with similar conditions is called the experience analogy method.
在本实施例中,应着重考虑下列因素:本矿和邻矿的冲击地压现状和发展趋势;本煤层或邻层、邻区已发生过冲击地压;煤层老顶为厚5m以上、单轴抗压强度大于70MPa的坚硬岩层;岛形或半岛形煤柱;支承压力影响区;上部或下部遗留煤柱或回采边界,煤层厚度或倾角突然变化区;褶曲或断裂构造带等。In this embodiment, the following factors should be considered: the current status and development trend of rockbursts in this mine and adjacent mines; rockbursts have occurred in this coal seam or adjacent layers and adjacent areas; the old roof of the coal seam is more than 5m thick and single. Hard rock formations with axial compressive strength greater than 70MPa; island-shaped or peninsula-shaped coal pillars; support pressure affected areas; upper or lower remaining coal pillars or mining boundaries, areas with sudden changes in coal seam thickness or dip angle; fold or fault structural zones, etc.
就整个矿井而言,通过地质条件和开采技术条件的分析,可以圈定冲击地压可能发生的区域。随着开采深度的延深,当煤岩体应力满足强度条件时就可能发生冲击地压。始发冲击的深度通常称临界深度。从始发深度起,冲击地压就可能在煤柱、煤层凸出部位和邻近煤柱的上下煤层区段发生,并随着开采水平的延深,冲击地压发生地点和范围也随之扩大。所有靠近采掘工作面的区域,煤层厚度和倾角突然变化的区域,以及地质构造带都可能成为发生冲击地压的危险区域。As for the entire mine, through the analysis of geological conditions and mining technology conditions, the area where rockburst may occur can be delineated. As the mining depth increases, rockburst may occur when the stress of the coal and rock mass meets the strength condition. The depth of the initial impact is usually called the critical depth. From the initial depth, rockbursts may occur in coal pillars, coal seam protrusions, and upper and lower coal seam sections adjacent to coal pillars. As the mining level increases, the location and scope of rockbursts also expand. . All areas close to the mining face, areas with sudden changes in coal seam thickness and inclination, and geological structural zones may become dangerous areas for rockbursts.
采用微震监测法来对采场支承应力带特征进行估测Using microseismic monitoring method to estimate the characteristics of stope support stress zone
利用短周期地震仪监测记录0.5级以上冲击发生的次数及冲击地压释放的能量。利用此趋势预测预报近期冲击地压发生的趋势及应力释放情况。在定位系统建成之前,采用现在的地震仪现行监测。通过上述流程,得到初步预测结果。Use short-period seismometers to monitor and record the number of shocks above magnitude 0.5 and the energy released by the impact ground pressure. Use this trend prediction to predict the trend of recent geoburst occurrence and stress release. Until the positioning system is built, current seismometers will be used for current monitoring. Through the above process, preliminary prediction results are obtained.
经过初步预测,结合喀斯特地貌的特殊性,上述的临界深度不应超过668m;工作面超前支承压力集中范围为5~35m,应力集中系数为2.5,但上方砾岩层的超前压力影响范围达120m。因此,工作面采动集中应力对工作面影响较为明显。分层开采时上分层工作面周期来压强度最大达510kn/m2,来压较为强烈。After preliminary prediction, combined with the particularity of karst landforms, the above-mentioned critical depth should not exceed 668m; the advanced support pressure concentration range of the working face is 5 to 35m, and the stress concentration coefficient is 2.5, but the impact range of the advanced pressure of the upper conglomerate layer reaches 120m. Therefore, the concentrated stress of mining on the working face has a more obvious impact on the working face. During layered mining, the cyclic pressure intensity of the upper layered working face reaches a maximum of 510kn/m 2 , and the pressure is relatively strong.
最后,构建模块基于初步预测结果,构建预测模型,得到最终预测结果,基于最终预测结果进行冲击能量的预警。Finally, the building module builds a prediction model based on the preliminary prediction results, obtains the final prediction results, and provides early warning of impact energy based on the final prediction results.
采用随机森林模型,基于初步预测结果,构建预测模型。A random forest model is used to build a prediction model based on the preliminary prediction results.
首先,收集由上述初步预测结果组成样本集。之后,采用均值代替法对收集到的样本集进行清洗;在本实施例中,主要为清洗缺失值、清洗与原数据类型不一致的内容、清洗掉不需要的数据,清洗逻辑错误数据。First, collect a sample set composed of the above preliminary prediction results. After that, the mean substitution method is used to clean the collected sample set; in this embodiment, the main steps are to clean missing values, clean content inconsistent with the original data type, clean unnecessary data, and clean logical error data.
之后,从样本集中筛选出特征数据构成特征集,利用随机森林模型基于建立的特征集采用bootstrap的方法从样本集的M个样本中随机抽取m个样本作为一个子训练集构建一颗决策树;采用相同的方法同步构建T(T>1)棵决策树;在特征集的n个特征中随机选择p个特征作为节点分裂的子集,根据平方误差选择p个特征误差最小的1个特征作为节点分裂特征,保持节点分裂直到该决策树不可再分裂;本实施例中,3个特征为根节点和内容节点,每个预测结果为输出也即叶子节点。对T棵决策树进行分裂组成随机森林;基于初步预测结果对分裂后的每颗决策树随机在M个样本中进行训练得到每颗决策树所对应的预测结果;对每颗决策树所对应的预测结果求取均值得到最终预测结果。After that, feature data are selected from the sample set to form a feature set, and a random forest model is used based on the established feature set and the bootstrap method is used to randomly select m samples from the M samples in the sample set as a sub-training set to build a decision tree; Use the same method to simultaneously construct T (T>1) decision trees; randomly select p features from the n features in the feature set as a subset of node splits, and select the feature with the smallest error among the p features as the subset based on the square error. The node split feature keeps the node split until the decision tree cannot be split anymore; in this embodiment, the three features are the root node and the content node, and each prediction result is the output, which is a leaf node. Split T decision trees to form a random forest; based on the preliminary prediction results, each split decision tree is randomly trained on M samples to obtain the prediction results corresponding to each decision tree; the prediction results corresponding to each decision tree are The prediction results are averaged to obtain the final prediction result.
最后基于最终预测结果进行冲击能量的预警。Finally, an early warning of impact energy is carried out based on the final prediction results.
以上所述的实施例仅是对本申请优选方式进行的描述,并非对本申请的范围进行限定,在不脱离本申请设计精神的前提下,本领域普通技术人员对本申请的技术方案做出的各种变形和改进,均应落入本申请权利要求书确定的保护范围内。The above-described embodiments are only descriptions of preferred modes of the present application and do not limit the scope of the present application. Without departing from the design spirit of the present application, those of ordinary skill in the art may make various modifications to the technical solutions of the present application. All modifications and improvements shall fall within the protection scope determined by the claims of this application.
Claims (8)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311078824.6A CN116910886B (en) | 2023-08-25 | 2023-08-25 | A method and system for early warning of impact energy in karst areas |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311078824.6A CN116910886B (en) | 2023-08-25 | 2023-08-25 | A method and system for early warning of impact energy in karst areas |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116910886A true CN116910886A (en) | 2023-10-20 |
CN116910886B CN116910886B (en) | 2024-05-07 |
Family
ID=88353306
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311078824.6A Active CN116910886B (en) | 2023-08-25 | 2023-08-25 | A method and system for early warning of impact energy in karst areas |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116910886B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200191994A1 (en) * | 2018-12-17 | 2020-06-18 | China University Of Petroleum (Beijing) | Modeling method and method for diagnosing lost circulation |
CN113486604A (en) * | 2021-06-16 | 2021-10-08 | 山东省鲁南地质工程勘察院(山东省地勘局第二地质大队) | Karst collapse early warning partitioning method based on numerical simulation |
CN116127752A (en) * | 2023-02-07 | 2023-05-16 | 中国矿业大学 | A method and system for predicting rock burst |
-
2023
- 2023-08-25 CN CN202311078824.6A patent/CN116910886B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200191994A1 (en) * | 2018-12-17 | 2020-06-18 | China University Of Petroleum (Beijing) | Modeling method and method for diagnosing lost circulation |
CN113486604A (en) * | 2021-06-16 | 2021-10-08 | 山东省鲁南地质工程勘察院(山东省地勘局第二地质大队) | Karst collapse early warning partitioning method based on numerical simulation |
CN116127752A (en) * | 2023-02-07 | 2023-05-16 | 中国矿业大学 | A method and system for predicting rock burst |
Non-Patent Citations (5)
Title |
---|
张以文;樊雷;付明勇;李桂森;马义华;: "基于GIS的矿山地质灾害管理系统构想", 岩土工程界, no. 06 * |
朱怀志;张峰;倪海明;: "鲍店煤矿保护层开采冲击诱因与防治体系", 煤矿安全, no. 04 * |
李宝富;刘永磊;: "冲击地压危险性等级识别的随机森林模型及应用", 科技导报, no. 01, pages 57 - 62 * |
王杰;罗周全;秦亚光;赵爽;: "基于随机森林理论的采场稳定性预测研究", 中国安全科学学报, no. 03 * |
罗利锐;刘志刚;闫怡冲;: "超前地质预报系统的提出及其发展方向", 岩土力学, no. 1 * |
Also Published As
Publication number | Publication date |
---|---|
CN116910886B (en) | 2024-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xie et al. | Study on the mechanical properties and mechanical response of coal mining at 1000 m or deeper | |
Xu et al. | Discrete element modeling of strata and surface movement induced by mining under open-pit final slope | |
Mazaira et al. | Intense rockburst impacts in deep underground construction and their prevention | |
Cao et al. | Risk assessment of fault water inrush during deep mining | |
Rahimi et al. | A comprehensive underground excavation design (CUED) methodology for geotechnical engineering design of deep underground mining and tunneling | |
CN111691872A (en) | Dynamic and accurate detection method for development range and degree of coal seam overlying strata fracture | |
Liu et al. | Analysis of support requirements for underground water-sealed oil storage cavern in China | |
Zhang et al. | Longwall mining–induced damage and fractures: Field measurements and simulation using FDM and DEM coupled method | |
Liu et al. | Numerical modelling of mining-induced inrushes from subjacent water conducting karst collapse columns in northern China | |
CN110159347A (en) | A kind of deep high stress tight roof stope dynamic disaster monitoring and pre-alarming method | |
CN106529762A (en) | New mine coal and gas outburst risk assessment method | |
Li et al. | Case study on the mining‐induced stress evolution of an extra‐thick coal seam under hard roof conditions | |
Zhu et al. | Simulation analysis of influencing factors of subsidence based on mining under huge loose strata: a case study of Heze mining area, China | |
Zhao et al. | Elastic wave prospecting of water-conducting fractured zones in coal mining | |
Jia et al. | Stability analysis of shallow goaf based on field monitoring and numerical simulation: A case study at an open-pit iron mine, China | |
Shi et al. | Study on the mechanisms and prevention of water inrush events in a deeply buried high-pressure coal seam—a case study of the Chensilou coal mine in China | |
Xue et al. | Numerical investigation on overburden migration behaviors in stope under thick magmatic rocks | |
Zhang et al. | Stability analysis and numerical simulation of foundation in old goaf under building load | |
Wei et al. | A case study of damage energy analysis and an early warning by microseismic monitoring for large area roof caving in shallow depth seams | |
Liu et al. | Study on stability of underlying room and pillar old goaf in close coal seam and mining of the upper coal seam | |
CN116910886B (en) | A method and system for early warning of impact energy in karst areas | |
Zhang et al. | Research on the underground water inrush mechanism based on the influence of fault | |
Wang et al. | A stress rotation–based method for improving roof stability of a deep longwall panel | |
Lu et al. | Floor Disturbance and Failure Characteristics of Super‐Large Mining Height Working Face | |
Zhou et al. | Impact of thick alluvial soil on a fractured water-conducting zone: an example from Huainan coal mine, China |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |