CN108226998B - Geological advance prediction method based on TSP system and 3D network of rock mass random discontinuities - Google Patents
Geological advance prediction method based on TSP system and 3D network of rock mass random discontinuities Download PDFInfo
- Publication number
- CN108226998B CN108226998B CN201711375524.9A CN201711375524A CN108226998B CN 108226998 B CN108226998 B CN 108226998B CN 201711375524 A CN201711375524 A CN 201711375524A CN 108226998 B CN108226998 B CN 108226998B
- Authority
- CN
- China
- Prior art keywords
- rock mass
- data
- discontinuous
- dimensional network
- discontinuous surface
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 48
- 230000008569 process Effects 0.000 claims abstract description 12
- 238000001514 detection method Methods 0.000 claims abstract description 4
- 238000004088 simulation Methods 0.000 claims description 15
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 8
- 238000009412 basement excavation Methods 0.000 claims description 7
- 238000013508 migration Methods 0.000 claims description 7
- 230000005012 migration Effects 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 238000001276 Kolmogorov–Smirnov test Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 3
- 230000000704 physical effect Effects 0.000 claims description 3
- 238000012552 review Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 2
- 238000012795 verification Methods 0.000 claims 1
- 238000010276 construction Methods 0.000 abstract description 11
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000005422 blasting Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000011160 research 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
- 235000019994 cava Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/612—Previously recorded data, e.g. time-lapse or 4D
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/641—Continuity of geobodies
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于TSP系统和岩体随机不连续面三维网络的地质超前预测方法,TSP系统现场采集数据,并对数据的有效性进行判别;对数据进行处理并复核,获取前方围岩三维网络模型数据,获取前期的勘探资料;划分岩体结构统计均质区和统计均质区内的不连续面优势组数;根据各组不连续面所对应的数据模拟生成不连续面空间三维网络数值模型;结合勘探资料、TSP分析结果和不连续面空间三维网络数值模型对前方岩体的稳定性进行综合超前预测。本发明以TSP的探测结果为原始数据,融合岩体随机不连续面三维网络对地下工程前方岩体稳定性进行准确预测预警,提高了预报的准确率,有力的指导地下工程的现场施工,保障工程一线的施工安全。
The invention discloses a geological advance prediction method based on a TSP system and a three-dimensional network of a random discontinuous surface of a rock mass. The TSP system collects data on site, and judges the validity of the data; processes and checks the data, and obtains the surrounding rock in front 3D network model data to obtain previous exploration data; divide the statistical homogeneous area of rock mass structure and the number of dominant groups of discontinuous surfaces in the statistical homogeneous area; generate 3D discontinuous surface space based on the data corresponding to each group of discontinuous surfaces Network numerical model; combined with exploration data, TSP analysis results and discontinuous surface space three-dimensional network numerical model to comprehensively predict the stability of the rock mass ahead. The invention uses the TSP detection results as the original data, integrates the three-dimensional network of the random discontinuous surface of the rock mass to accurately predict and warn the stability of the rock mass ahead of the underground project, improves the accuracy of the forecast, effectively guides the on-site construction of the underground project, and guarantees Construction safety on the front line of the project.
Description
技术领域technical field
本发明涉及一种基于TSP系统和岩体随机不连续面三维网络的地质超前预测方法。The invention relates to a geological advance prediction method based on a TSP system and a three-dimensional network of random discontinuous surfaces of rock masses.
背景技术Background technique
我国科技水平和综合国力的不断提高,公路、铁路和城市地铁等交通工程取得了举世瞩目的成就,与此同时,也带动起地下工程的建设进入跨越式的发展轨道。一批地质条件复杂、施工技术难度大的长大隧道相继建成。目前,中国已是世界上隧道和地下工程最多、发展速度最快、地质及结构形式最复杂的国家,在建和待建的山岭隧道、地铁隧道、水底隧道、水电压力隧道等地下工程不计其数。在地下工程开挖工程中,由于岩体被结构面切割成大小不同、形状不一的各种岩块,随着地下工程的持续开挖,围岩原有的静力平衡状态被打破,某些暴露在临空面上的不稳定块体会沿着软弱结构面滑移并导致局部掉块,严重的甚至产生连锁反应,在地下工程一定范围内造成围岩的局部坍塌、失稳,严重影响工程进度,人员、设备的安全。因此,为了保证地下工程施工安全,提高工程建设的效率,地下工程中的超前地质探测更显重要。With the continuous improvement of my country's scientific and technological level and comprehensive national strength, transportation projects such as highways, railways and urban subways have achieved world-renowned achievements. At the same time, it has also driven the construction of underground projects into a leapfrog development track. A number of long tunnels with complicated geological conditions and difficult construction techniques have been built one after another. At present, China is the country with the most tunnels and underground projects, the fastest development speed, and the most complicated geological and structural forms in the world. There are countless underground projects such as mountain tunnels, subway tunnels, underwater tunnels, and hydroelectric pressure tunnels under construction or to be built. number. In the excavation of underground engineering, because the rock mass is cut into various rock blocks of different sizes and shapes by the structural surface, with the continuous excavation of the underground engineering, the original static equilibrium state of the surrounding rock is broken, and a certain Some unstable blocks exposed on the air surface will slide along the weak structural surface and cause local block loss, and even cause a chain reaction in serious cases, causing local collapse and instability of surrounding rock within a certain range of underground engineering, seriously affecting Project progress, safety of personnel and equipment. Therefore, in order to ensure the safety of underground engineering construction and improve the efficiency of engineering construction, advanced geological detection in underground engineering is more important.
目前,TSP超前地质预报系统采用地震波反射原理,操作简单,能长距离地预报地下工程掌子面前方的地质情况。通过在掘进面后方一定距离内的钻孔中施以微型爆破来发射信号,爆破引发的地震波在岩体中以球面的形式向四周传播,其中一部分向地下工程前方传播,当地震波遇到岩石波阻抗差异界面(如断层、破碎带、岩性变化、溶洞和地下水等)时,部分地震信号反射回来,反射信号经接收传感器转换成电信号并放大,对反射信号进行相应的计算就可确定不良地质体相关参数。At present, the TSP advanced geological prediction system adopts the principle of seismic wave reflection, which is easy to operate and can predict the geological conditions in front of the underground engineering face for a long distance. The signal is transmitted by micro blasting in the borehole within a certain distance behind the excavation face. The seismic wave caused by the blasting propagates in the rock mass in the form of a spherical surface, and part of it propagates to the front of the underground project. When the seismic wave encounters the rock wave When the impedance difference interface (such as faults, broken zones, lithological changes, caves and groundwater, etc.), part of the seismic signal is reflected back, the reflected signal is converted into an electrical signal by the receiving sensor and amplified, and the corresponding calculation of the reflected signal can determine the fault. Geological related parameters.
岩体的稳定性是受岩体中不连续面的控制的。在进行岩体的稳定性分析前,首要的问题是搞清岩体的地质条件特别是不连续面在岩体中的展布特征。因为,不连续面反映了岩体不连续和不均匀的本质。然而岩体中不连续面的展布不具备有规可循的特征,它们在岩体中是随机分布的,而且其几何形态是测不准的,因此,采用统计推断和概率的方法进行研究。三维网络模拟正是建立在概率论和统计学的基础之上,它的研究对象是不连续面空间的几何特征。简化假定的不连续面几何形态,通过概率统计方法和空间解析几何等方法改正现场取样偏差,求解正确的迹长﹑大小﹑空间密度﹑产状分布等几何参数,通过蒙特-卡罗方法构筑不连续面在三维空间的组合形态,形成三维网络模型。The stability of a rock mass is controlled by discontinuities in the rock mass. Before the stability analysis of the rock mass, the primary problem is to find out the geological conditions of the rock mass, especially the distribution characteristics of the discontinuities in the rock mass. Because the discontinuous surface reflects the discontinuous and uneven nature of the rock mass. However, the distribution of discontinuities in the rock mass does not have the characteristics that can be followed. They are randomly distributed in the rock mass, and their geometric shapes are unpredictable. Therefore, statistical inference and probability methods are used for research. . The 3D network simulation is just based on probability theory and statistics, and its research object is the geometric characteristics of discontinuous surface space. Simplify the geometric form of the assumed discontinuity surface, correct the on-site sampling deviation through the method of probability statistics and spatial analytic geometry, solve the correct geometric parameters such as trace length, size, space density, and occurrence distribution, and construct the discontinuous surface through the Monte-Carlo method. The combination of continuous surfaces in three-dimensional space forms a three-dimensional network model.
传统地下工程超前地质预报对现场采集到的TSP数据,通过对应的后处理软件处理,获得P,SH和SV波的时间剖面、深度偏移剖面、提取的反射层、岩石物理力学参数、各反射层能量大小等成果,对获得的成果进行解译,可以预测前方的软弱带、破碎带、断层、含水情况等。TSP系统虽然在地下工程超前预报中广泛应用,但并没有结合岩体随机不连续面三维网络对地下工程前方岩体的稳定性进行预测。如何在地下工程中基于TSP系统和岩体随机不连续面三维网络对岩体稳定性超前预测,是目前亟待解决的一个技术难题。Advanced geological forecasting of traditional underground engineering processes the TSP data collected on site through corresponding post-processing software to obtain time profiles, depth migration profiles, extracted reflection layers, rock physical and mechanical parameters, and reflection profiles of P, SH, and SV waves. Interpretation of the obtained results can predict the weak zone, broken zone, fault, water content, etc. ahead. Although the TSP system is widely used in the advanced prediction of underground engineering, it does not combine the three-dimensional network of random discontinuities of the rock mass to predict the stability of the rock mass ahead of the underground engineering. How to predict the stability of rock mass in advance based on the TSP system and the 3D network of rock mass random discontinuities in underground engineering is a technical problem that needs to be solved urgently.
发明内容Contents of the invention
本发明为了解决上述问题,提出了一种基于TSP系统和岩体随机不连续面三维网络的地质超前预测方法,本发明在TSP系统基础上,以TSP的探测的前面围岩的产状信息(如倾向,倾角等)为原始数据,融合岩体随机不连续面三维网络对地下工程前方岩体稳定性进行的预测预警,准确、高效的判别前方围岩稳定性,降低由不稳定块体引发的隧道塌方失稳灾害的概率,有力的指导地下工程的现场施工,保障工程一线的施工安全。In order to solve the above-mentioned problems, the present invention proposes a geological advance prediction method based on the TSP system and the three-dimensional network of the random discontinuity surface of the rock mass. The present invention is based on the TSP system, with the occurrence information ( Such as inclination, inclination, etc.) as the original data, combined with the three-dimensional network of random discontinuity of the rock mass to predict and warn the stability of the rock mass ahead of the underground project, accurately and efficiently determine the stability of the surrounding rock ahead, and reduce the risk caused by unstable blocks. It can effectively guide the on-site construction of underground projects and ensure the construction safety of the front line of the project.
为了实现上述目的,本发明采用如下技术方案:In order to achieve the above object, the present invention adopts the following technical solutions:
一种基于TSP系统和岩体随机不连续面三维网络的地质超前预测方法,包括以下步骤:A geological advance prediction method based on a TSP system and a three-dimensional network of random discontinuities in rock mass, comprising the following steps:
(1)TSP系统现场采集数据,并对数据的有效性进行判别;(1) The TSP system collects data on site and judges the validity of the data;
(2)对数据进行处理并复核,获取前方围岩三维网络模型数据,获取前期的勘探资料;(2) Process and review the data, obtain the 3D network model data of the surrounding rock in front, and obtain the previous exploration data;
(3)划分岩体结构统计均质区和统计均质区内的不连续面优势组数;(3) divide the statistical homogeneous area of rock mass structure and the number of discontinuous surface dominant groups in the statistical homogeneous area;
(4)根据各组不连续面所对应的数据模拟生成不连续面空间三维网络数值模型;(4) According to the corresponding data simulation of each group of discontinuous surfaces, a three-dimensional network numerical model of discontinuous surfaces is generated;
(5)结合勘探资料、TSP分析结果和不连续面空间三维网络数值模型对前方岩体的稳定性进行综合超前预测。(5) Combining the exploration data, TSP analysis results and the numerical model of the three-dimensional network of the discontinuity surface, the stability of the rock mass in front is comprehensively predicted in advance.
进一步的,所述步骤(1)中,对数据有效性进行判别,是指在现场测试过程中选取对数据干扰小于阈值的信号。Further, in the step (1), judging the validity of the data refers to selecting a signal whose interference with the data is less than a threshold during the on-site test.
进一步的,所述步骤(2)中,TSP系统数据信息包括深度偏移图、三维显示图、产状信息和岩体2D物性图。Further, in the step (2), the TSP system data information includes depth migration map, three-dimensional display map, occurrence information and rock mass 2D physical property map.
进一步的,所述步骤(2)中,结合前期的勘探资料和TSP系统处理后生成的深度偏移图对开挖面前方围岩情况进行分析并做出初步预测。Further, in the step (2), analyze the surrounding rock conditions in front of the excavation face in combination with the previous exploration data and the depth migration map generated after processing by the TSP system, and make a preliminary prediction.
所述步骤(3)中,对岩体结构统计均质区划分,不同的岩体结构具有不同的岩体力学特征和岩体水力学特性,采用概率关联表施密特图对比法,结合实际的地质特点进行岩体结构统计均质区的划分。In the described step (3), the rock mass structure is statistically divided into homogeneous regions, and different rock mass structures have different rock mass mechanical characteristics and rock mass hydraulic properties, and the probability correlation table Schmidt diagram comparison method is adopted, combined with the actual According to the geological characteristics of the rock mass structure, the statistical homogeneity area is divided.
所述步骤(4)中,不连续面空间三维网络数值模型,过程包括:不连续面迹长的偏差校正及迹长大小模拟,不连续面间距、密度的模拟,不连续面产状测量偏差的校正或/和蒙特卡罗模拟及模型检验。In the step (4), the three-dimensional network numerical model of the discontinuous surface, the process includes: the deviation correction of the discontinuous surface trace length and the simulation of the trace length size, the simulation of the discontinuous surface spacing and density, and the discontinuous surface occurrence measurement deviation Calibration or/and Monte Carlo simulation and model checking.
所述步骤(4)中,采用建立在概率统计基础上的端点估值器方法进行不连续面迹长的纠正、圆盘直径大小的模拟、分布密度函数的拟合以及间距测量的偏差校正。In the step (4), an endpoint estimator method based on probability and statistics is used to correct the discontinuous track length, simulate the diameter of the disk, fit the distribution density function and correct the deviation of the distance measurement.
所述步骤(4)中,不连续面间距、密度的模拟过程为对不连续面的间距进行模拟,采用X2和K-S检验方法进行概率密度函数的最佳拟合,获得间距的概率分布密度函数及其对应的参数;采用Oda提出的张量法进行不连续面空间密度的求算。In described step (4), the simulation process of discontinuous surface distance, density is to simulate the distance of discontinuous surface, adopt X 2 and KS test method to carry out the best fitting of probability density function, obtain the probability distribution density of distance Functions and their corresponding parameters; the tensor method proposed by Oda is used to calculate the space density of discontinuous surfaces.
所述步骤(4)中,根据已经获得了各组不连续面在空间中的平均直径,结合权重分配公式,对每一个不连续面分配不同的权重,然后进行不连续面产状测量偏差校正。In the step (4), according to the average diameter of each group of discontinuous surfaces in space, combined with the weight distribution formula, each discontinuous surface is assigned a different weight, and then the discontinuous surface occurrence measurement deviation correction is performed .
所述步骤(4)中,蒙特卡罗模拟及模型检验的过程包括:In described step (4), the process of Monte Carlo simulation and model inspection comprises:
(4-1)确定生成网络模型的规模;(4-1) Determine the scale of the generated network model;
(4-2)确定生成网络模型的规模形态中不连续面的个数;(4-2) Determine the number of discontinuous surfaces in the scale form of the generated network model;
(4-3)确定每一条不连续面的空间位置,采用了蒙特卡罗模拟的方法,随机产生各不连续面的中心点坐标;(4-3) Determine the spatial position of each discontinuous surface, adopt the method of Monte Carlo simulation, randomly generate the center point coordinates of each discontinuous surface;
(4-4)确定各不连续面直径的大小:用蒙特卡罗法来模拟生成不连续面直径并使其服从已知的最佳概率分布;(4-4) Determine the size of each discontinuous surface diameter: use the Monte Carlo method to simulate and generate the discontinuous surface diameter and make it obey the known best probability distribution;
(4-5)确定不连续面的产状;(4-5) Determine the occurrence of the discontinuous surface;
(4-6)将所有模拟的结果进行组合,从而形成一个完整的模型,并对其进行检验。(4-6) Combine the results of all simulations to form a complete model and test it.
与现有技术相比,本发明的有益效果为:Compared with prior art, the beneficial effect of the present invention is:
本发明是基于TSP系统和岩体随机不连续面三维网络的岩体稳定性超前预测方法。这种方法以TSP的探测结果为原始数据,融合岩体随机不连续面三维网络对地下工程前方岩体稳定性进行准确的预测预警,提高了预报的准确率,有力的指导地下工程的现场施工,保障工程一线的施工安全。The invention is a rock mass stability advance prediction method based on the TSP system and the rock mass random discontinuous surface three-dimensional network. This method uses the detection results of TSP as the original data, and integrates the 3D network of random discontinuities of the rock mass to accurately predict and warn the stability of the rock mass ahead of the underground project, which improves the accuracy of the forecast and effectively guides the on-site construction of the underground project. , to ensure the construction safety of the front line of the project.
附图说明Description of drawings
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。The accompanying drawings constituting a part of the present application are used to provide further understanding of the present application, and the schematic embodiments and descriptions of the present application are used to explain the present application, and do not constitute improper limitations to the present application.
图1是本发明方法实现流程图。Fig. 1 is the flow chart of the realization of the method of the present invention.
具体实施方式:Detailed ways:
下面结合附图与实施例对本发明作进一步说明。The present invention will be further described below in conjunction with the accompanying drawings and embodiments.
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。It should be pointed out that the following detailed description is exemplary and intended to provide further explanation to the present application. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。It should be noted that the terminology used here is only for describing specific implementations, and is not intended to limit the exemplary implementations according to the present application. As used herein, unless the context clearly dictates otherwise, the singular is intended to include the plural, and it should also be understood that when the terms "comprising" and/or "comprising" are used in this specification, they mean There are features, steps, operations, means, components and/or combinations thereof.
在本发明中,术语如“上”、“下”、“左”、“右”、“前”、“后”、“竖直”、“水平”、“侧”、“底”等指示的方位或位置关系为基于附图所示的方位或位置关系,只是为了便于叙述本发明各部件或元件结构关系而确定的关系词,并非特指本发明中任一部件或元件,不能理解为对本发明的限制。In the present invention, terms such as "upper", "lower", "left", "right", "front", "rear", "vertical", "horizontal", "side", "bottom" etc. indicate The orientation or positional relationship is based on the orientation or positional relationship shown in the accompanying drawings, and is only a relative term determined for the convenience of describing the structural relationship of the various components or elements of the present invention, and does not specifically refer to any component or element in the present invention, and cannot be understood as a reference to the present invention. Invention Limitations.
本发明中,术语如“固接”、“相连”、“连接”等应做广义理解,表示可以是固定连接,也可以是一体地连接或可拆卸连接;可以是直接相连,也可以通过中间媒介间接相连。对于本领域的相关科研或技术人员,可以根据具体情况确定上述术语在本发明中的具体含义,不能理解为对本发明的限制。In the present invention, terms such as "fixed", "connected" and "connected" should be understood in a broad sense, which means that they can be fixedly connected, integrally connected or detachably connected; they can be directly connected or can be connected through the middle The medium is indirectly connected. For relevant researchers or technical personnel in the field, the specific meanings of the above terms in the present invention can be determined according to specific situations, and should not be construed as limitations on the present invention.
如图1所示,一种基于TSP系统和岩体随机不连续面三维网络的地质超前预测方法,包括以下步骤:As shown in Figure 1, a geological advance prediction method based on the TSP system and the three-dimensional network of random discontinuities in the rock mass includes the following steps:
步骤一:TSP系统现场采集数据,并对数据的有效性进行判别;Step 1: The TSP system collects data on site and judges the validity of the data;
步骤二:对数据进行处理并复核,获取前方围岩三维网络模型数据;Step 2: Process and review the data, and obtain the 3D network model data of the surrounding rock in front;
TSP系统数据处理后可得到开挖面前方围岩的深度偏移图、三维显示图、产状信息和岩体2D物性图。After the data processing of the TSP system, the depth migration map, three-dimensional display map, occurrence information and 2D physical property map of the rock mass in front of the excavation face can be obtained.
结合前期的勘探资料和TSP系统处理后生成的深度偏移图对开挖面前方围岩情况(断层、破碎带、含水情况等)进行分析并做出初步预测。Combining the previous exploration data and the depth migration map generated after processing by the TSP system, the surrounding rock conditions (faults, broken zones, water content, etc.) in front of the excavation face were analyzed and a preliminary prediction was made.
通过分析结论结合获取的产状信息提取相应的三维网络数值模型数据。The corresponding three-dimensional network numerical model data is extracted by combining the analysis conclusions with the obtained occurrence information.
步骤三:对岩体结构统计均质区和统计均质区内的不连续面优势组数进行划分;Step 3: dividing the statistically homogeneous area of the rock mass structure and the number of discontinuous surface dominant groups in the statistically homogeneous area;
不同的岩体结构具有不同的岩体力学特征和岩体水力学特性。采用概率关联表施密特图对比法,结合实际的地质特点进行岩体结构统计均质区的划分。Different rock mass structures have different rock mass mechanical characteristics and rock mass hydraulic properties. The statistical homogeneity area of rock mass structure is divided by using the Schmidt diagram comparison method of probability correlation table and combining with the actual geological characteristics.
对每一个统计均质区内的不连续面进行优势组数的划分。The number of dominant groups is divided for the discontinuous surfaces in each statistically homogeneous area.
步骤四:将各组不连续面所对应的数据模拟生成不连续面空间三维网络数值模型;Step 4: Simulate the data corresponding to each group of discontinuous surfaces to generate a three-dimensional network numerical model of the discontinuous surface space;
不连续面迹长的纠正、圆盘直径大小的模拟、分布密度函数的拟合以及间距测量的偏差校正。校正偏差的方法主要是采用建立在概率统计基础上的端点估值器方法。Correction of discontinuous track length, simulation of disc diameter, fitting of distribution density function, and deviation correction of spacing measurement. The method of correcting the deviation mainly adopts the endpoint estimator method based on probability and statistics.
不连续面间距、密度的模拟。首先对不连续面的间距进行模拟,采用X2和K-S检验方法进行概率密度函数的最佳拟合,可以获得间距的概率分布密度函数及其对应的参数;采用Oda提出的张量法进行不连续面空间密度的求算。Simulation of discontinuous face spacing, density. Firstly, simulate the spacing of the discontinuous surface, and use the X2 and K-S test methods to perform the best fitting of the probability density function, and obtain the probability distribution density function of the spacing and its corresponding parameters; use the tensor method proposed by Oda to carry out the discontinuous surface space Calculation of density.
不连续面产状测量偏差校正。根据已经获得了各组不连续面在空间中的平均直径,结合权重分配公式,对每一个不连续面分配不同的权重,然后进行校正。Discontinuity surface occurrence measurement deviation correction. According to the average diameter of each group of discontinuous surfaces in space, combined with the weight distribution formula, different weights are assigned to each discontinuous surface, and then correction is performed.
蒙特卡罗模拟及模型检验。生成3-D网络模型需要确定生成网络模型的规模,其形态是矩形盒子。过程如下:①确定方盒的规模;②确定方盒中不连续面的个数;③确定每一条不连续面的空间位置,采用了蒙特卡罗模拟的方法,随机产生各不连续面的中心点坐标;④确定各不连续面直径的大小:用蒙特卡罗法来模拟生成不连续面直径并使其服从已知的最佳概率分布;⑤确定不连续面的产状;⑥模型组合:把上述各步骤模拟的结果进行组合,从而形成一个完整的模型,并对其进行检验。Monte Carlo simulation and model checking. Generating a 3-D network model needs to determine the scale of the generated network model, and its shape is a rectangular box. The process is as follows: ① determine the scale of the box; ② determine the number of discontinuous surfaces in the box; ③ determine the spatial position of each discontinuous surface, and use the method of Monte Carlo simulation to randomly generate the center of each discontinuous surface Point coordinates; ④Determine the size of each discontinuous surface diameter: use Monte Carlo method to simulate and generate the discontinuous surface diameter and make it obey the known best probability distribution; ⑤Determine the occurrence of the discontinuous surface; ⑥Model combination: Combine the simulation results of the above steps to form a complete model and test it.
步骤五:综合解译。结合勘探资料、TSP分析结果和不连续面空间三维网络数值模型对前方岩体的稳定性进行综合超前预测。Step five: comprehensive interpretation. Combined with the exploration data, TSP analysis results and the numerical model of the three-dimensional network of the discontinuity surface, the stability of the front rock mass is comprehensively predicted in advance.
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。The above descriptions are only preferred embodiments of the present application, and are not intended to limit the present application. For those skilled in the art, there may be various modifications and changes in the present application. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of this application shall be included within the protection scope of this application.
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。Although the specific implementation of the present invention has been described above in conjunction with the accompanying drawings, it does not limit the protection scope of the present invention. Those skilled in the art should understand that on the basis of the technical solution of the present invention, those skilled in the art do not need to pay creative work Various modifications or variations that can be made are still within the protection scope of the present invention.
Claims (8)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711375524.9A CN108226998B (en) | 2017-12-19 | 2017-12-19 | Geological advance prediction method based on TSP system and 3D network of rock mass random discontinuities |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711375524.9A CN108226998B (en) | 2017-12-19 | 2017-12-19 | Geological advance prediction method based on TSP system and 3D network of rock mass random discontinuities |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108226998A CN108226998A (en) | 2018-06-29 |
CN108226998B true CN108226998B (en) | 2019-12-24 |
Family
ID=62652405
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711375524.9A Expired - Fee Related CN108226998B (en) | 2017-12-19 | 2017-12-19 | Geological advance prediction method based on TSP system and 3D network of rock mass random discontinuities |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108226998B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110096775A (en) * | 2019-04-20 | 2019-08-06 | 青岛理工大学 | Method for determining underground engineering excavation state |
CN111561934A (en) * | 2020-06-24 | 2020-08-21 | 平湖市中地测绘规划有限公司 | Geological exploration planning method based on unmanned aerial vehicle |
CN115310184B (en) * | 2022-08-23 | 2025-05-16 | 山东大学 | A method and system for analyzing discontinuous deformation of tunnel surrounding rock stability |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101743554A (en) * | 2007-03-30 | 2010-06-16 | 帕拉代姆地球物理有限公司 | Partitioning algorithm for building a stratigraphic grid |
CN103513293A (en) * | 2013-10-12 | 2014-01-15 | 广西大学 | Tunnel geology comprehensive advanced forecasting expert system and implementation method thereof |
CN103901503A (en) * | 2014-03-25 | 2014-07-02 | 中冶集团武汉勘察研究院有限公司 | Comprehensive detection method for front unfavorable geologic bodies in mine underground roadway tunneling |
CN105334548A (en) * | 2015-10-20 | 2016-02-17 | 中冶交通建设集团有限公司 | Geological forecasting method for tunnel construction in karst area |
CN106570287A (en) * | 2016-11-10 | 2017-04-19 | 中国人民解放军理工大学 | Method for predicting water inflow of tunnel based on three-dimensional discrete fracture network |
CN107589471A (en) * | 2017-07-13 | 2018-01-16 | 高军 | A kind of Railway Tunnel Synthetic Geological Prediction Ahead of Construction method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102221711B (en) * | 2011-05-23 | 2013-05-01 | 吉林大学 | Advanced prediction device for tunnel water inrush by using nuclear magnetic resonance differential detection and detection method |
-
2017
- 2017-12-19 CN CN201711375524.9A patent/CN108226998B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101743554A (en) * | 2007-03-30 | 2010-06-16 | 帕拉代姆地球物理有限公司 | Partitioning algorithm for building a stratigraphic grid |
CN103513293A (en) * | 2013-10-12 | 2014-01-15 | 广西大学 | Tunnel geology comprehensive advanced forecasting expert system and implementation method thereof |
CN103901503A (en) * | 2014-03-25 | 2014-07-02 | 中冶集团武汉勘察研究院有限公司 | Comprehensive detection method for front unfavorable geologic bodies in mine underground roadway tunneling |
CN105334548A (en) * | 2015-10-20 | 2016-02-17 | 中冶交通建设集团有限公司 | Geological forecasting method for tunnel construction in karst area |
CN106570287A (en) * | 2016-11-10 | 2017-04-19 | 中国人民解放军理工大学 | Method for predicting water inflow of tunnel based on three-dimensional discrete fracture network |
CN107589471A (en) * | 2017-07-13 | 2018-01-16 | 高军 | A kind of Railway Tunnel Synthetic Geological Prediction Ahead of Construction method |
Also Published As
Publication number | Publication date |
---|---|
CN108226998A (en) | 2018-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101739716B (en) | Three-dimensional space structure modeling and key block identifying method of engineering rock mass | |
CN101770038B (en) | Intelligent positioning method of mine microquake sources | |
WO2024169098A1 (en) | Geological determination method for delayed extremely-intense rockburst | |
CN104656124A (en) | Multi-parameter comprehensive rock burst predicting method based on geophysical exploration method | |
Zhang et al. | Determination of RVE with consideration of the spatial effect | |
CN109492262A (en) | A method of utilizing numerical simulation analysis non-uniform Distribution crack Drift stability | |
CN108226998B (en) | Geological advance prediction method based on TSP system and 3D network of rock mass random discontinuities | |
CN110174463B (en) | A Nondestructive Quantitative Testing Method for Three-Dimensional Mining Stress Field of Working Face | |
Fan et al. | Advanced stability analysis of the tunnels in jointed rock mass based on TSP and DEM | |
Zhang et al. | Determination of the optimal threshold and length measurements for RQD calculations | |
CN106980716B (en) | Underground cavern stability analysis method based on random block | |
Yang et al. | Numerical simulations of full-wave fields and analysis of channel wave characteristics in 3-D coal mine roadway models | |
Zhang et al. | Statistical analysis of stochastic blocks and its application to rock support | |
Dong et al. | Quantitative evaluation and classification method of the cataclastic texture rock mass based on the structural plane network simulation | |
CN116258025B (en) | Method for simulating landslide stability of structural mixed rock based on geophysical prospecting data | |
CN115510527A (en) | Tunnel surrounding rock stability discrimination and quantitative evaluation method based on safety index | |
CN110187379A (en) | A kind of test method based on TSP method tunnel geological forecast effect | |
CN117216935A (en) | Method for determining reserved deformation of horizontal broken stratum excavation driven by multiple processes | |
CN111239813B (en) | Seismic wave advanced prediction detection method for tunnel water-containing geological structure | |
CN106149770A (en) | Bored concrete pile pile foundation construction period synchronizes the large-section in-situ concrete pile hole wall Rock Mass Integrality detection method carried out | |
CN107832505B (en) | Advance prediction method of ill-conditioned blocks in underground engineering based on TSP system and block theory | |
Wang et al. | Research on multi-frequency ultrasonic scanning detecting technology of cavity in the test borehole | |
CN111538073A (en) | Advanced geological forecast detection simulation device and method for tunnel poor geologic body | |
CN117372639A (en) | Three-dimensional geological-geophysical modeling method combining sparse diaphysis section with auxiliary surface at any position | |
Wang et al. | Reliability analysis of rock slope excavation considering the stochasticity and finite persistence of wedges |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191224 Termination date: 20211219 |
|
CF01 | Termination of patent right due to non-payment of annual fee |