文摘

提出了结构系统可靠性的定量方法研究随机岩石力学参数和负载的影响稳定的煤层底板承压含水层之上。获得地板的可靠性概率,提出了两种防水地板失败模式如下:(1)矿业完全移除地板的防水关键层。在这种情况下,失效模式的三种主要类型:采矿失败,承压水侵入,联合开采失败和入侵。(2)失效模式包括剪切和拉伸失败当关键层的厚度大于0。根据弹性薄板理论,计算的性能功能的可靠性概率模式可以获得。失效模式被视为系列系统。蒙特卡洛法计算每个失效模式的可靠性概率和系列系统。结果表明,随机关键层的岩石力学参数和负荷显著影响antiwater侵入的能力。此外,突水系数和可靠性概率可以同时作为突水风险的评价指标。这两个索引可以提高地板的可靠性的评估。

1。介绍

从挖掘上面层灰岩含水层突水,通常有较高的水压力,是一个有问题的中国煤田水文地质与工程地质领域,它严重威胁着煤矿的生产。因此,应该研究和水损害到我楼煤矿突水事故的发生,这是一个困难的问题要解决,应该有效地抑制(1- - - - - -3]。几项研究已经试图调查底板突水问题并提出了各种方法,比如突水系数(4- - - - - -6),板模式7- - - - - -9),三个地区(10- - - - - -12),和关键层(13- - - - - -17]。这些研究揭示了突水机理考虑各个方面,为矿山安全生产提供了积极的指导。然而,所有这些研究只考虑确定性方法不能反映真正的安全度。防水地板的地质条件复杂,导致岩石力学参数和负荷不确定,考虑到安全评价防水地板采用确定性方法是不确定的。相比之下,可靠性理论更适合在挖掘学习防水地板的稳定性。

使用结构可靠性方法研究的结果分析煤层的防水地板的稳定性高于缺乏承压含水层。地板的可靠性分析模型与突水系数作为性能指标建立了杜et al。18]。陆et al。19和卢和姚20.杜)扩展的模型等。18)利用剩余完整岩层为固支梁。在他们的模型中,限制水压力的梁的强度失效基于材料力学和弹性力学的理论被用作性能函数计算的可靠性防水地板的稳定性。然而,如上所述,这些研究只关注单个失效模式导致的突水和没有考虑多个失效模式的可能性也会引起突水。例如,拉伸或剪切破坏可能发生在有效防水地板或采矿失败类型和进步的承压水的入侵可能发生在一个无效的防水层。另一方面,突水系数只能反映两个因素,即水压力和厚度的防水层,在选择结构性能功能。然而,固定梁模型来源于材料力学很差精度分析的深梁相对较短的长度和高度。尽管弹性力学固支梁模型的准确性高,固支梁的应力分布时无法准确获得,导致相当大的错误的解决方案限制水压力。结构可靠性分析结果的准确性主要取决于性能函数的选择。因此,以前的计算结果并不完全可靠。

能够合理地反映防水地板的稳定的可靠性,本研究首先关键层理论适用于考虑地应力的影响。地板的关键层被视为一个矩形板和四个边夹,和弹性力学板理论是用来获得关键层破坏的力学判据,这被认为是一个函数的结构可靠性分析。在此基础上,分析了防水地板的失效模式,并计算每个失效模式的可靠性概率仿真采用蒙特卡罗方法。失效模式被视为系列系统的可靠性概率防水地板的稳定系统。随机变量的均值和方差,比如岩石强度参数和地面压力,通过案例讨论,影响可靠性的关键层的位置概率进行了分析。最后,各种失效模式和系统可靠性的防水地板某煤矿进行了分析和评价。

2。方法

2.1。机械故障形式的关键层

煤层的垂直距离的屋顶承压含水层叫做防水地板和厚度 (图1(一))。煤层开采之前,由于存在天然裂缝在地板上,沿着这些天然裂缝水上升到一定高度静水压力的作用下,原来被称为进步上升侵入承压水的高度。高度通常是小,如图1(一)

工作面开采后,采空区的形成和上覆地层的重量转移到采空区周围煤岩体形成采空区周围支承压力带。当支持压力达到或超过岩石强度极限限制的地板,一定范围内的岩体工作面的地板被摧毁和地板的失败区采矿形成(图1 (b))。同时,矿业和承压水的综合影响下,最初的入侵区扩展的裂纹进一步向上逐步入侵区承压水。最初的和进步的和入侵腰带被称为入侵区承压水,它有一个最大高度 目前,该参数的最大值是相对较大的结构断裂带和通常是通过实际测量测量。

对工作面,如图1,失败的最大深度 地板矿业在初始阶段可以根据断裂力学理论确定。采场被认为是一个无限岩石裂纹内部的一部分。采场的开采厚度比开采宽度小得多。因此,应力的计算可以简化为平面问题,和失败的最大深度 在采场的地板可以计算(1),如下所示21]: 在哪里 的平均体重是煤层上覆层(kN / m3), 是封面的煤层深度(米) 平均凝聚力(MPa)和内摩擦角(°)的岩石,分别和 的宽度是工作面(m)。

解决后的深度 地板的采矿失败的厚度 剩余的完整岩石地层可获得根据指令的高度 承压水(图2)。底板的突水主要取决于antiwater产能剩余的完整岩石地层压力。关键层理论指出,一个合理的岩体强度和耐水剩下完整岩石带,称为关键层的防水地板,有厚度 关键层的地板之间的距离和煤层的地板 (图1 (b))。承压水是否可以进入工作面主要取决于关键层的防水性能,这是畸形的类似于水压力下一盘。如图2正常的开采条件下,对长壁开采的脸,地板的关键层由一个矩形板和四个边缘夹轻轻接近水平或倾斜煤层时,和平板模型建立。

的设置x- - - - - -y协调与板的中面表面一致,协调的起源在于板的中心,盎司向中间平面轴是垂直向下的,板的厚度是一样的,防水的关键层 x方向是工作面的宽度和长度 y方向是工作面推进方向的长度

的平均体重地板的采矿失败带 ;剩下的完整岩石层的平均体重是多少 ;关键层的平均体重是多少 ; , , , 表示泊松比、弹性模量、内摩擦角,分别和粘结强度。板的上部受到自重 地板的矿业裂缝带和自重 ( )岩石的层,这是获得剩下的完整岩石带之间的区别,关键层。下部受到统一的水压力 自重的关键层被认为是表面力 在板上。因此,关键层的下部受到表面力 有趣的是,上部岩体的防水区是毁在这一点上;我们相信它不会承担任何额外的压力,将在和谐与关键层向上弯曲。因此,它是可行的和可控的考虑上部岩体的自重的关键层的外部负载模型。水平应力 徒在盘子里。假设一个均布荷载的计算。表达式显示为(2)。 在哪里 侧压力系数, 封面的煤层深度(米)认为,然后呢 的平均体重是煤层上覆层(kN / m3)。

考虑到厚板理论不成熟的力学、机械的分析关键层采用弹性薄板理论。鉴于关键层不被挖掘,他们仍然可以被视为一个连续介质假定为均匀各向同性。在中国,长壁工作面的宽度是150 - 240,首次屋顶的加权一步工作面大约是30 - 40米,和地板的突水通常发生在首次屋顶加权一步(8]。同时,含煤沉积岩体的层厚度通常低于10米,和防水板模型的关键层在初始开采可以满足薄板的一般条件(厚度,宽度之比不大于1/5)。

共同作用下的纵向载荷 和横向负载 在关键层的下部,薄板弯曲的微分方程可以表示为22,23]

的方程, 薄板的挠度, 薄板的弯曲刚度。

边界条件矩形薄板的四条边夹紧

上述微分方程准确是很难解决,由里兹方法可以解决。的挠度函数弯曲板的中面 在哪里 挠度函数的待定系数。

显然,挠度函数 可以满足固支板的边界条件。因为上面的方程收敛很快,方程可以取而代之的是第一项,可以满足采矿工程的精度要求;也就是说,

根据里茨法,挠度函数的系数 的关键层复合垂直和水平载荷下可以获得如下:

因此,挠度函数, ,关键层的开采层如下:

方程(8)表明,关键层的挠度增加水平地应力的作用下。根据弹性薄板理论,应力表达式关键层垂直和水平载荷下的水阻力可以获得如下(压应力是积极的和拉应力是负的):

在接下来的公式,

方程(9)表明,关键层的应力分量是线性分布沿板厚度在纵向和横向负载下,和最大值位于板的上下表面( )。讨论了关键层的应力分布规律。

的值 , , , , , , , , ,的分布 , , 在垂直和水平载荷下关键层可以通过使用(9)和(10),如图3(在此,只有上表面的状况 在分析了关键层下表面的压力是一个相反的分布规律)。

从图的分析3和(8点(中心), )在关键层的上表面拉应力最大值的经验σx= 5.08 MPa和σy= 23.34 MPa。中点(z =−0±15日hk/ 2)长边的关键层压应力最大值的经验σx= 13.34 MPa和σy= 32.81 MPa。剪切应力的分布 相反是对称中心的关键层中心点的值为零,最大剪切应力是中心的四个小矩形(图3 (c)),最大值为2.43 MPa。同样,压应力达到最大值的中心点防水关键层的下表面σx= 15.08 MPa和σy= 33.33 MPa,代表这些最大震级大于上表面的压应力。下表面的压应力达到最大值的中心点长边σx= 3.34 MPa和σy= 22.81 MPa,代表震级小于最大拉应力在表面。在相同的剪切应力xy坐标是一样的,在板的上下表面,但在相反的方向。

上述分析结果表明,底板关键层产生一个向上弯曲的共同作用下,横向载荷和纵向水平地应力越低,导致中间的拉上盘和下表面的边界地区,以及紧迫的边界区域的上表面和下盘的中部地区。根据弹性力学计算的结果,本文采用与数值相比,最大拉应力位于中心的关键层的上表面,最大压应力位于中心的下表面。这个结果表明关键层的中心点是突水的薄弱区域的地板上,应给予更多的关注。这个结果不同于使用固定梁分析获得的结论。

根据关键层的应力分布规律和岩石的失效特点,两种主要失效形式出现在地板上关键层在垂直和水平载荷:(1)剪切破坏发生在压缩关键层的一部分,(2)拉伸断裂发生在拉伸部分。不同的应用程序武力地板的关键层将改变地板的水控制策略。关键层的应力分量的表达式推导出在前面的部分中,和相应的屈服准则可以用来获得机械评价公式在关键层两种失效形式。

前者分析表明,最大压应力是中心点(x= 0,y= 0,z=hk/ 2)的下表面防水关键层。这个点的主应力可以通过(9)。

的公式,

针对摩尔-库仑采用的屈服准则。 在哪里

根据(11),(12)和(13),临界水压力(没有剪切破坏),发生在关键层

最大拉应力位于中心点(x= 0,y= 0,z=−hk/ 2)关键层的顶面。根据(9),此时的最大拉伸应力

为了防止关键层的拉伸断裂,最大拉应力(σ3)应小于允许的抗拉强度 岩石, 临界水压力(没有拉伸断裂),发生在关键层表示如下:

2.2。不透水的可靠性计算基本稳定

稳定系统的可靠性计算的防水地板主要可以归结为两个问题:(1)确定的主要失效模式防水地板和(2)计算的可靠性稳定性系统的防水地板上。

结构可靠性的概率是指工程结构将完成预定功能在给定的时间内,在指定的条件下24]。正如上面提到的,防水地板上的两种主要失效模式导致地上失去抵抗函数如下:(1)开采影响下,关键层影响地板的采矿失败深度和入侵承压水的高度,导致关键层的厚度为零,也就是说,关键层是防水的。失效模式可以分为三个subfailure模式在这种情况下。第一个是当采矿失败的深度 达到关键层的底部,也就是说, ;这就是所谓的采矿失败类型的失效模式。第二个是当入侵的高度 的承压水达到防水关键层的顶部,也就是说, ;这就是所谓的入侵类型的失效模式。第三个是当 ,和地板的采矿失败带接触承压水的入侵区,也就是说, ;这叫做采矿失败的组合类型的失效模式和入侵(2)当层关键层的厚度满足 ,地板上失败是主要关键层本身的强度不足造成的。关键层的失败可分为两种模式:剪切和拉伸断裂类型。图4总结了各种失效模式的防水地板上

代表的基本随机变量因素的数量 影响地板的防水功能,如岩体参数,故障层深度( ),入侵承压水的高度( ), ,和侧压力系数 当负载电阻(R- - - - - -年代)模式,可以列出每个失效模式的性能函数如下给出。

的性能函数失败模式挖掘和入侵类型如下:

,地板是一个可靠的状态。当 ,防水层失败。当 ,得到的表达式称为地板结构的极限状态方程。根据可靠性的计算方法,可获得不同失效模式的失效概率如下:

然后,失效模式的可靠性概率如下:

同样,合并后的失效模式的性能函数的采矿失败和入侵

采矿失败的组合类型的失效模式和入侵时可能发生的失效模式不发生。因此,可以计算失效概率的条件下 ,这一个条件概率。因此,失效模式的失效概率,这是采矿失败和入侵的组合类型,如下:

同样,地板的可靠性概率

当地板的厚度的关键层不是零,关键层的失败主要是由于缺乏关键层的力量。剪切破坏的表现功能

拉伸断裂的功能

同样,剪切或拉伸断裂失效模式的关键层可能会发生的情况 ,和失效模式的失效概率也是一个条件概率。失败的概率如下所示。

剪切破坏的失效模式

拉伸断裂的失效模式

同样,防止关键层的可靠性概率出现剪切或拉伸失效模式如下:

没有发生剪切破坏,

没有发生拉伸断裂,

计算可靠性的方法包括分析和仿真的方法。一阶二次矩方法是广泛使用的一种分析方法。最流行的仿真方法是蒙特卡洛方法,目前被认为是一个相对准确的方法计算的可靠性。考虑到各种失效模式的极限状态方程的复杂性在地板上突水推导出在这项研究中,计算导数当使用分析方法是困难的。因此,采用蒙特卡罗方法对可靠性计算。这个计算的基本原理如下。

在基本随机变量的联合概率密度函数 是设置为 ,这样地板结构的失效概率 在哪里 的指示功能吗 , , (25]。

方程(32)表明,失效概率的期望值是指示函数。不同的价值观 使用蒙特卡罗方法得到执行批量抽样的随机变量 ,平均值是地板失效概率的估计价值。根据地板的各种性能函数失败,每个失效模式的失效概率和可靠性概率可以用蒙特卡罗方法计算。

系统的可靠性是对多个性能函数结构可靠性问题。单一结构破坏主要是由于不同的失效模式。当发生任何故障模式时,整体结构损伤发生,和结构系统可以看作是一系列系统组成的各种失效模式;当所有失效模式发生结构性破坏,结构系统可以被视为一个并行系统。先前的分析表明,防水地板的任何失效模式的发生将导致突水的发生。因此,结构体系是一系列系统。

一系列系统的失效概率是每个失效模式的总和。性能的函数 th系列系统设置为失效模式 ,和一系列系统的失效概率的公式

当系列系统的失效概率计算的蒙特卡罗方法使用(33根据(),32),指示函数的级数系统成为[25]

系统的可靠性概率稳定防水地板是计算的方法用来计算单一失效模式的可靠性概率根据(34)。

3所示。结果与讨论

工作面的宽度 , 是最初的屋顶重量一步工作面,40米;封面的煤层深度H= 600,承压含水层的水压力p= 2.5 MPa,防水层的总厚度h= 45米,关键层的厚度hk= 8 m,其底部的距离煤层地板是35米。基本随机变量列在表中1。泊松比(μ)的岩石和土壤质量和重量 仅略有改变,因此它们可以被视为常数,用以下值用于计算: , , , 所有的随机变量服从正态分布的参数之间的相关性,和采样的数量计算的可靠性概率是100万次。在程序中,的值 从每个抽样计算应用于防水关键层来确定 需要改变。

3.1。在稳定可靠性概率随机变量的影响

每个失效模式的可靠性概率研究通过改变单个随机变量的平均值。数据56分别显示结果时平均凝聚力 地板和水平压力系数 是作为实例变量。

来自前分析,随机变量 影响价值的深度挖掘失败 的地板上,影响了关键层的实际厚度 地板上挖掘的失效模式的可靠性概率失败类型,结合入侵和采矿失败,关键层的剪切和拉伸断裂模式都是受变化的影响 (图5)。然而,当 高于6 MPa,每个失效模式的高水平和可靠性的概率系统延迟明显增加而增加的意思 这一现象表明,当防水地板的强度达到一定值,它不能有效改善防水地板的稳定性的可靠性概率通过不断改善地板强度。失效模式的入侵类型不受随机变量的影响 和可靠性概率是稳定在大约0.99998(图5(一个))。

在这项研究中,随机变量 影响故障模式发生概率的关键层的剪切和拉伸断裂系统。图6显示的可靠性概率关键层的剪切破坏模式随的均值的增加而减小 ,而拉伸失效模式的可靠性概率增加。当侧压力系数 较小,水平应力较小,关键层的拉伸断裂更可能发生剪切破坏。地应力的增加,产生的张应力横向弯曲偏移,部分由于水平压应力的存在;因此,在关键层拉应力发生的可能性减少,而剪切破坏的可能性增加。当水平应力的增加在一定程度上(如1.5),发生剪切破坏的可能性超过拉伸断裂的类型。在这个瞬间,系统可靠性概率达到最大值为0.96359。之后,增加的平均值 ,系统的可靠性概率不断减少。

此外,上述失效模式的可靠性概率基本稳定在大约0.99818,0.99998,和0.99994,由于采矿失败类型的地板上,和入侵类型和故障模式的两个是独立的

3.2。变异系数稳定可靠性概率的影响

改变单一的变异系数值随机变量使每个失效模式的可靠性概率的研究。图7显示结果当平均凝聚力 地板和入侵的承压水的高度 是作为例子。

7表明,变异系数增加 ,与相关的失效模式 和系统的可靠性概率在不同程度降低。当变异系数达到1.1,系统的可靠性概率大大降低。当变异系数 达到2,系统的失效模式从拉伸断裂类型更改为采矿失败。入侵破坏模式不受随机变量的影响 ,及其可靠性概率基本稳定在大约0.99998。

随机变量的值 影响关键层的实际厚度值,因此改变变异系数的值 和影响可靠性的概率入侵类型的失效模式,联合开采的失败和入侵类型和关键层的剪切和拉伸断裂类型(图8)。同样,变异系数 增加,每个失效模式的相关性和可靠性的概率系统不同程度的减少。可靠性概率减少明显的方式当变异系数达到2.6。拉伸断裂模式总是防水地板系统的主要故障模式如图所示的分析系统的可靠性概率。采矿失败不受随机变量的影响 ,因此,其可靠性概率基本稳定在大约0.99818。

上述分析表明,随机变量的变异系数是成反比的稳定概率防水地板上。各种失效模式的可靠性概率和系统可靠性概率和变异系数的增加减少。变异系数的影响在防水地板上稳定具有重要意义。

在第一个失效模式, , , 被定义为安全因素的失效模式挖掘故障类型,入侵类型,分别和两者的结合。相应地,在第二个主要失效模式,安全系数定义为临界水压力的比值(7)和(9)实际水压下关键层的剪切和拉伸失败。计算安全系数的平均值随机变量列在表中1,即平均安全系数。并给出了计算结果表2

2显示,当安全系数是用来评估的稳定性防水地板,只是价值的均值影响变量。岩体的物理力学参数和载荷的随机性不充分的考虑。数据78表明,变异系数增加,可靠性的概率减少下相同的安全系数。当安全系数很大,地板的可靠性概率仍相对较低。因此,地板的稳定性无法保证完全安全系数。可靠性理论是更客观和合理的使用来进行防水地板的稳定性研究。

3.3。关键层位置稳定可靠性概率的影响

分析关键层位置的影响地板的稳定性和可靠性概率,每个失效模式的概率在不同比率可以通过选择不同比率的计算 结果如图9,显示比例的增加,也就是说,当关键层的煤层开采的影响降低,采矿失败的失效模式的可靠性概率类型逐渐增加。失效模式的可靠性概率的入侵类型和组合类型逐渐从一个常数下降状态,这表明了关键层逐渐受承压水的侵入。关键层的剪切和拉伸断裂模式和系统的可靠性概率增加和减少的比率时,到达峰值 大约是0.7。上述发现表明关键层的位置显著影响防水地板的稳定性和可靠性概率;当关键层远离开采的影响范围和承压水的侵入,防水地板的稳定性和可靠性概率可以达到最大值。

4所示。应用程序

将研究成果应用于10-coal煤层的开采在第103 Yangliu煤矿矿区。Yangliu矿位于苏州城,中国安徽省。第103届矿区位于中间的我。封面10-coal缝深度是500 - 650,这是一个向东倾斜的单斜结构。地面高程+ 25米,顶部的距离太原组灰岩含水层的(即是45米。防水层的总厚度h= 45 m)和地层的岩性致密,主要是泥岩、粉砂岩、细砂岩(图10)。根据石灰石抽水试验的结果在太原形成矿区,静态的覆盖深度灰岩含水层的水位是−38米,单位涌水量= 0.0855 L / (s·m)和水丰富是虚弱的。然而,灰岩裂缝发展的异质性和当地富含水分差异很大,这表明一个潜在的水风险10-coal煤层的安全开采。采矿开始流失2014年石灰石矿区的水。当前的石灰岩覆盖深度340米后根据最新的观测结果脱水试验。

地板防水层的岩性分析表明,细砂层厚7.9米24.6米的煤层下出现一个稳定的分布。因此,岩石层可以被视为一个防水关键层,也就是说,hk= 7.9 m。矿山压力监测表明,屋顶的加权一步工作面在这个矿是40米。

的高度进步的入侵在结构层承压水很大断裂带和小裂缝发展的薄弱区域。据统计检测数据在矿区多年来,中值 是12米,标准差是2.1米,0.18,变异系数较大。侧压力系数的平均值 是1.2,标准偏差为0.08,变异系数在0.07并不大。

当计算地板的深度挖掘失败,均值和标准差的地板上岩体的物理力学参数计算使用(35)。

的公式, th分层厚度的地板, 力学参数的均值的吗 th层, 力学参数的标准偏差的 层。

3列出了岩石参数的随机变量的特征值(转换后的35根据测试数据)岩石的物理力学性质的实验室。泊松比和重量微妙的变化也视为常数的计算 , , , 各种随机变量在表3是正态分布,而不考虑计算过程中的参数之间的相关性。

根据封面的深度10-coal缝地板,可靠性概率计算从500 - 630米的深度,表中给出4。表4表明,水压力随深度,而安全系数的平均值和各失效模式的可靠性概率地板的减少和失效概率增加。

现场技术人员和我的员工经常使用突水系数法。这种方法被定义的煤矿防治水规定(2009)使用以下数学表达式: 在哪里 突水系数,MPa / m; 地板上的水压力防水层,MPa;和 地板防水层的总厚度,m。

根据中国的煤矿水控制的要求,结合淮北矿区的实际地质条件,整个地板的突水系数不超过0.06 MPa / m。利用这种方法有助于确定这个矿区的最高安全开采深度为610米。然而,计算结果(表4)表明,各种失效模式的平均安全系数较大时,开采深度是610米。此外,地板的可靠性概率仅为84.5%,而失效概率高达15.5%。报告的可靠性概率远远低于先前的研究在文献[26];具体来说,道路建设的可靠性概率不应低于90%。上面的突水系数的计算表明,如果采用传统的定值方法的最大挖掘深度10-coal seam是由安全系数,然后从地板下容易发生突水不利的外部因素。因此,它是不全面和不可靠的确定突水的风险从地板上只使用固定值的方法。结构系统可靠性的理论应该引入的风险预测突水的地板上。

结构系统可靠性理论的应用评价的煤矿的突水楼上承压含水层仍在试验阶段,在实践中需要进一步测试。合理分布模型的不确定因素,如岩体参数、负载,应该建立和模型大小,通过大量数据的积累更可靠的评价结果。

5。摘要和结论

(1)随机变量的均值的变化有重要影响地板的水压力抵抗能力。更大的平均值地板的内聚力对应于一个高可靠性的概率。然而,当地板内聚力达到一定值,地板的稳定性和可靠性概率不能有效地改善。一个更大的侧压力系数 显示更高水平压应力。的可靠性概率底板关键层的剪切破坏类型减少,拉伸断裂类型的可靠性概率增加,和系列的可靠性概率系统最初增加和减少。系统的可靠性概率达到最大值时的可靠性概率剪切破坏类型和强度失效模式是相同的(2)更大的随机变量的方差系数对应于一个较小的相关失效模式和系统的可靠性概率;然而,平均安全系数保持不变(3)关键层的位置有重要影响的稳定性和可靠性概率防水地板上。当关键层远离采矿和承压水的影响范围逐步入侵,防水地板的稳定性和可靠性概率可以达到最大值(4)的应用案例表明,监管和指定的安全系数计算的系统可靠性概率模型在这个研究可以作为评价指标来获得底板突水风险的评价结果更加准确和全面

数据可用性

使用的数据来支持本研究的发现可以从相应的作者。

的利益冲突

作者宣称没有利益冲突。

确认

本研究支持的部分研究经费来自中国的国家自然科学基金(41472235,41472235,41472235,41501570),煤炭资源安全开采的国家重点实验室,中国矿业大学&技术(没有。SKLCRSM16KFD01),安徽省自然科学基金(1508085 qe89),安徽大学优秀青年人才项目(gxyq2017004)、广东省重点实验室的土壤和地下水污染控制(2017号b030301012)和国家环境保护重点实验室综合表面Water-Groundwater污染控制。本文使用的数据可以通过联系相应的作者直接访问。