一种底板隔水岩层稳定性计算方法与流程

1.本发明属于岩土灾害技术领域，具体涉及一种底板隔水岩层稳定性计算方法。


背景技术：

2.煤炭的发展在经济建设中有着举足轻重的地位，其在国民经济和人民生活中的起着中坚力量的作用。现阶段，很多煤炭均分布于地下水层附近，很多煤炭资源受到水害的严重威胁。近些年，煤岩层底板突水事故时有发生，以灰岩承压水带来的煤岩层底板突水隐患尤为突出，合理科学的提前排查底板突水隐患、降低甚至消除底板突水威胁就显得极为重要。
3.在现有技术中，对底板隔水层稳定性方面在以往并没有系统性的理论求解方法，同时，对底板隔水岩层安全的评价方法均未兼顾底板采动破坏深度与承压水导升高度双重影响，且现有的方法未求得最小隔水层厚度和最大极限水压力值。另外，现阶段，现场实际判断底板安全性依然参考突水系数法进行，突水系数法是采用统计归纳法所得到经验结果，针对具体的某一个煤矿或采面底板稳定性安全评价没有现成的理论求解方法，缺乏理论支撑，且没有考虑底板岩石力学性质，无系统性的理论求解方法，在实际操作过程中，无法给工程实践人员带来科学的指导。


技术实现要素：

4.针对上述现有技术存在的问题，本发明提供一种底板隔水岩层稳定性计算方法，该方法能推导得出底板隔水岩层厚度安全阈值和最大极限水压力值，可科学合理的判断底板整体安全性，为高承压水体上开采煤岩层提供安全保障和理论指导。
5.为了实现上述目的，本发明提供一种底板隔水岩层稳定性计算方法，若煤岩层底板岩性相近，则执行步骤一；若煤岩层底板岩性存在明显的上下两层分层结构，则执行步骤二；
6.步骤一：以单一岩性结构底板隔水岩层进行稳定性分析，具体方法如下：
7.s1：建立力学模型，进行受力分析；具体方法如下：
8.s10：将岩煤开采层至承压开采层之间的部分简化为一种底板岩性结构的底板岩层，再将一种底板岩性结构的由上到下依次分为底板采动裂隙层、有效隔水层和底板承压水导升裂隙层；
9.s11：根据公式(1)计算底板的有效隔水层高度h1；
10.h1＝h1‑
h0‑
c
ꢀꢀꢀ
(1)；
11.式中，h1为底板岩层的厚度，单位m；h0为底板采动裂隙层高度；c为底板承压水导升裂隙层高度；
12.s12：根据公式(2)计算底板承压水导升裂隙层顶部的等效水压力q2；根据公式(3)计算采空区冒落矸石的等效荷载q1；
13.q2＝q0‑
(h1 h0)
·
γ1ꢀꢀꢀ
(2)
[0014][0015]
式中，q0为高承压水实际作用的水压力；γ1为底板岩层的等效重度，单位kn/m3；l为采空区沿走向长度；γ为采空区冒落矸石的等效重度，单位kn/m3；
[0016]
s2：引入弹塑性破坏理论进行分析；具体方法如下：
[0017]
s20：将有效隔水层和底板采动裂隙层简化成理想的线性弹塑性体，并使有效隔水层的变形满足公式(4)，使底板采动裂隙层的变形满足公式(5)；
[0018]
σ
x
＝k
·
σ
y
σ
n
ꢀꢀꢀ
(4)；
[0019][0020]
式中，σ
x
为底板岩层中微单元体在x方向的水平应力，单位mpa；σ
y
为底板岩层中微单元体在y方向的垂直应力，单位mpa；σ
n
为底板岩层单轴抗压强度，单位mpa，或者通过实测获得；c
n
为底板岩层的粘聚力，单位mpa；为底板岩层的内摩擦角，单位
°
；
[0021]
为底板岩层残余强度，单位mpa，或者通过实测获得；为底板岩层的残余粘聚力，单位mpa；为底板岩层的残余内摩擦角，单位
°
；当底板岩层在弹性和残余阶段的内摩擦角不变时，即则有k
*
＝k；
[0022]
s3：选取微分单元，建立微分平衡方程；具体方法如下：
[0023]
从底板隔离层中选取一段厚度为dy的微单元体，以水平方向为x轴、垂直向下为y轴建立底板岩层微单元体进行受力分析，其中dσ
y
、dσ
x
分别为微单元体在垂直和水平方向的主应力；使微单元体在y方向上的应力平衡方程满足公式(6)；并根据公式(7)计算微单元体与两侧岩层接触面上的摩擦阻力τ；根据公式(6)建立y方向应力平衡方程得到公式(8)；
[0024]
∑f＝0
ꢀꢀꢀꢀꢀꢀꢀꢀ
(6)；
[0025]
τ＝f
·
σ
x
c
ꢀꢀꢀꢀꢀꢀꢀ
(7)；
[0026]
σ
y
·
l 2τ
·
dy
‑
(σ
y
dσσ
y
)
·
l＝0
ꢀꢀꢀꢀ
(8)；
[0027]
式中，c为微单元体与两侧岩层接触面上的粘聚力，单位mpa，当微单元体处于底板采动裂隙层h0时，当微单元体处于有效隔水层h1时，c＝c
n
；f为微单元体与两侧岩层接触面上的摩擦系数，当微单元体处于底板采动裂隙层h0时该参数为f
*
，当微单元体处于有效隔水层h1时为f，
[0028]
s4：解微分方程，推导底板岩层所能承受的最大极限水压力值或最小隔水岩层厚度安全阈值；
[0029]
当微单元体处于底板采动裂隙层时，按如下方法进行求解：
[0030]
s41：联立公式(6)、公式(7)和公式(8)获得公式(9)，并求解公式(9)得到公式(10)；
[0031][0032][0033]
式中，为未知量；
[0034]
s42：先使y＝0时，得到公式(11)；并联立公式(3)和公式(11)得到公式(12)；再推导出公式(13)；
[0035][0036][0037][0038]
s42：使y＝h0，得到公式(14)；
[0039][0040]
当微单元体处于有效隔水层时，按如下方法进行求解：
[0041]
s43：联立公式(5)、公式(7)和公式(8)获得公式(15)，并求解公式(15)得到公式(16)；
[0042][0043][0044]
式中，c1为未知量；
[0045]
s44：先使y＝h0 h1，得到公式(17)；并联立公式(2)和公式(17)得到公式(18)；再联立公式(16)和公式(18)得到(19)；
[0046][0047][0048]
[0049]
s45：使y＝h0，得到公式(20)；在y＝h0交界面处，公式(20)等于公式(14)，进而得到公式(21)，再求解公式(21)得到公式(22)；
[0050][0051][0052][0053]
s5：根据实际问题，判断隔水层厚度是否安全，具体过程如下：
[0054]
当底板隔水岩层厚度h1已知时，通过公式(22)给出底板岩层所能承受的极限水压力值，实际水压力q与q0存在公式(23)中的关系：当底板隔水岩层受水压力值已知时，对公式(22)求解得到底板岩层最小厚度值h
min
，底板岩层实际厚度值h1与最小厚度值h
min
存在公式(24)中的关系：
[0055][0056][0057]
步骤二：以双岩性结构底板隔水岩层进行稳定性分析，具体方法如下：
[0058]
a1：建立力学模型，进行受力分析；具体方法如下：
[0059]
a10：将岩煤开采层至承压开采层之间的部分简化为两种底板岩性结构的底板岩层，再将两种底板岩性结构的底板岩层由上到下分别记为底板隔水岩层a和底板隔水岩层b，并将底板隔水岩层a由上到下依次分为底板采动裂隙层和上部有效隔水层，将底板隔水岩层b由上到下依次分为下部有效隔水层和承压水导升裂隙层，其中上部有效隔水层和下部有效隔水层形成总有效隔水层；
[0060]
a11：根据公式(25)计算底板的总有效隔水层高度；
[0061]
h2 h3＝h2‑
h0‑
c
ꢀꢀꢀ
(25)；
[0062]
式中，h2为底板岩层的总厚度，单位m，当底板隔水层实际厚度确定来求解承压水压力最大值时，h2＝h
实
；h0为底板采动裂隙层高度；c为底板承压水导升裂隙层高度；
[0063]
a12：根据公式(26)计算承压水导升裂隙层顶部的等效水压力q3；根据公式(27)计算采空区冒落矸石的等效荷载q1；
[0064]
q3＝q0‑
(h2 h0)
·
γ2‑
h3·
y3ꢀꢀꢀ
(26)；
[0065][0066]
式中，q0为高承压水实际作用的水压力；h2为上部有效隔水层的高度；h0为底板采动裂隙层高度；h3为下部有效隔水层的高度；γ2为底板隔水岩层a的等效重度，单位kn/m3；γ3为底板隔水岩层b的等效重度，单位kn/m3；l为采空区沿走向长度；γ为采空区冒落矸石的等效重度，单位kn/m3；
[0067]
a2：引入弹塑性破坏理论进行分析；具体方法如下：
[0068]
a20：将总有效隔水层和底板采动裂隙层简化成理想的线性弹塑性体，并使底板隔水岩层a中的上部有效隔水层的变形满足公式(27)，使底板隔水岩层b中的下部有效隔水层的变形满足公式(28)；使底板采动裂隙层的变形满足公式(29)；
[0069]
σ
2x
＝k2·
σ
2y
σ
2n
ꢀꢀꢀꢀ
(27)；
[0070]
σ
3x
＝k3·
σ
3y
σ
3n
ꢀꢀꢀꢀꢀ
（28)；
[0071][0072]
式中，σ
2x
为底板隔水岩层a中微单元体在x方向的水平应力，单位mpa；σ
2y
为底板隔水岩层a中微单元体在y方向的垂直应力，单位mpa；σ
2n
为底板隔水岩层a的单轴抗压强度，单位mpa；或者通过实测获得；c
2n
为底板隔水岩层a的粘聚力，单位mpa；为底板隔水岩层a的内摩擦角，单位
°
；
[0073]
σ
3x
为底板隔水岩层b中微单元体在x方向的水平应力，单位mpa；σ
3y
为底板隔水岩层b中微单元体在y方向的垂直应力，单位mpa；σ
3n
为底板隔水岩层b的单轴抗压强度，单位mpa；或者通过实测获得；σ
3n
为底板隔水岩层b的粘聚力，单位mpa；为底板隔水岩层b的内摩擦角，单位
°
；
[0074]
为底板隔水层顶部岩层a的残余强度，单位mpa，或者通过实测获得；为隔水层顶部岩层a的残余粘聚力，单位mpa；为隔水层顶部岩层a的残余内摩擦角，单位
°
；隔水层顶部岩层a在弹性和残余阶段的内摩擦角不变，即则有k
*
＝k；
[0075]
a3：选取微分单元，建立微分平衡方程；具体方法如下：
[0076]
从底板隔离层中选取一段厚度为dy的微单元体，以水平方向为x轴、垂直向下为y轴建立底板隔水层微单元体进行受力分析，其中dσ
y
、dσ
x
分别为微单元体在垂直和水平方向
的主应力；使微单元体在y方向上的微单元体应力的平衡方程满足公式(30)；并根据公式(31)计算微单元体与两侧岩层接触面上的摩擦阻力τ；根据公式(31)建立y方向应力平衡方程得到公式(32)；
[0077]
∑f＝0
ꢀꢀꢀꢀꢀꢀꢀꢀ
(30)；
[0078]
τ＝f
·
σ
x
c
ꢀꢀꢀꢀꢀꢀꢀ
(31)；
[0079]
σ
y
·
l 2τ
·
dy
‑
(σ
y
dσ
y
)
·
l＝0
ꢀꢀꢀꢀ
(32)；
[0080]
式中，c为微单元体与两侧岩层接触面上的粘聚力，单位mpa，当微单元体处于底板采动裂隙层h0时，当微单元体处于上部有效隔水层h2时，c＝c
2n
；当微单元体处于下部有效隔水层下部时，c＝c
3n
；f为微单元体与两侧岩层接触面上的摩擦系数，当微单元体处于底板采动裂隙层h0时该参数为f
*
，当微单元体处于上部有效隔水层h2时时为f2，当微单元体处于上部有效隔水层h3时为f3，
[0081]
a4：解微分方程，推导底板岩层所能承受的最大极限水压力值或最小隔水岩层厚度安全阈值；
[0082]
当微单元体处于底板采动裂隙层时，按如下方法进行求解：
[0083]
a41：联立公式(29)、公式(31)和公式(32)获得公式(33)，并求解公式(33)得到公式(34)；
[0084][0085][0086]
式中，为未知量；
[0087]
a42：先使y＝0时，得到公式(35)；并联立公式(27)和公式(35)得到公式(36)；再推导出公式(37)；
[0088][0089][0090][0091]
a42：使y＝h0，得到公式(38)；
[0092][0093]
当微单元体处于上部有效隔水层h2时，按如下方法进行求解：
[0094]
a43：联立公式(27)、公式(31)和公式(32)获得公式(39)，并求解公式(39)得到公式(40)；
[0095][0096][0097]
式中，c2为未知量；
[0098]
a44：使y＝h0，得到公式(41)；在y＝h0交界面处，公式(38)等于公式(41)，进而得到公式(42)，再联立公式(42)与公式(40)得到公式(43)；
[0099][0100][0101][0102]
a45：把y＝h0 h2代入公式(43)得到公式(44)；
[0103][0104]
当微单元体处于下部有效隔水层h3时，按如下方法进行求解：
[0105]
a46：把公式(28)和公式(31)代入公式(32)得到公式(45)；并求解公式(45)得到公式(46)；
[0106][0107][0108]
式中，c3为未知量；
[0109]
a47：使y＝h0 h2 h3得到公式(47)；并联立公式(24)和公式(47)得到公式(48)，联立公式(28)和公式(46)得到公式(49)；
[0110][0111][0112][0113]
a48：在y＝h0 h2交界面处，公式(44)等于公式(49)，进而得到公式(50)，再求解公式(50)得到公式(51)；
[0114][0115][0116]
a5：根据实际问题，判断隔水层厚度是否安全，具体过程如下：
[0117]
当底板隔水岩层的总厚度h2已知时，通过公式(51)给出底板隔水层所能承受的极限水压力值，实际水压力q与q0存在公式(52)中的关系：当底板隔水岩层受水压力值已知时，对公式(51)求解得到底板隔水层最小厚度值h
min
，底板隔水层实际厚度值h2与最小厚度值h
min
存在公式(53)中的关系：
[0118][0119][0120]
在步骤一中，底板采动裂隙层和承压水导升裂隙层的高度由现场实测或相关理论计算得出。
[0121]
在步骤二中，底板采动裂隙层和承压水导升裂隙层的高度由现场实测或相关理论计算得出。
[0122]
本发明通过求解给出了一种底板隔水层稳定性判别方法，并推导得出底板隔水岩层厚度安全阈值和最大极限水压力值，可科学合理的判断底板整体安全性，为高承压水体上采煤提供安全保障和理论指导。通过构建煤岩层开采层至承压含水层之间的底板隔水岩层整体稳定性力学模型，建立了一种底板隔水岩层承受极限承载水压力和采动应力耦合作用的理论计算方法，并给出了底板隔水岩层极限水压力和最小隔水层厚度阈值理论解，从而为研判底板安全稳定及突水风险性提供技术保障。
[0123]
相比较于其他传统的底板破坏深度分析理论，本发明更加切合实际，分析角度新颖，具有明显优势。本发明不仅可以推导得出单层、双层煤岩层底板所能承受的临界水压力，更能推广至多层煤岩层底板所能承受的最大极限水压力和最小底板安全厚度阈值；进而本发明可用计算单层、双层结构底板隔水岩层稳定性，还可以推广到多层底板结构稳定性的验算。其适用于综采工作面、综放工作面、炮采工作面、连采工作面等一切开采条件，适用于煤岩层开采和岩石保护层开采中的奥灰水防治、寒灰水防治和太灰水防治，也适用岩石层开采、金属矿开采、半煤岩层开采、巷道掘进、隧道开挖、地铁开挖等一切地下工程。同时，还适用于求解巷道底板安全厚度等相关煤矿防治水等一切矿井灾害防治领域。
附图说明
[0124]
图1是本发明中底板隔水层在承压水压力作用下的截面应力分析图；
[0125]
图2本发明中底板隔水层单元体受力分析图；
[0126]
图3是本发明中双岩性结构底板隔水层在承压水压力作用下的截面应力分析图；
[0127]
图4是本发明中双岩性结构底板隔水层单元体受力分析图；
[0128]
图5是本发明中双岩性结构底板极限水压力随工作面推进度变化曲线。
具体实施方式
[0129]
下面结合附图对本发明作进一步说明。
[0130]
本发明提供了一种底板隔水岩层稳定性计算方法，若煤岩层底板岩性相近，则执行步骤一；若煤岩层底板岩性存在明显的上下两层分层结构，则执行步骤二；
[0131]
步骤一：以单一岩性结构底板隔水岩层进行稳定性分析，具体方法如下：
[0132]
s1：考虑到回采工作面采用综采或综放工艺，开采宽度较大，可将其沿工作面回采走向方向看作平面应变问题处理，沿开采层走向做剖面，建立力学模型，则底板隔水层力学模型图如图1所示，进行受力分析；具体方法如下：
[0133]
s10：将岩煤开采层至承压开采层之间的部分简化为一种底板岩性结构的底板岩层，再将一种底板岩性结构的由上到下依次分为底板采动裂隙层、有效隔水层和底板承压水导升裂隙层；
[0134]
s11：根据公式(1)计算底板的有效隔水层高度h1；
[0135]
h1＝h1‑
h0‑
c
ꢀꢀꢀ
(1)；
[0136]
式中，h1为开采层底板至寒武系灰岩顶部之间底板岩层的厚度，单位m；当底板隔水层实际厚度确定来求解承压水压力最大值时，h1＝h
实
；h0为底板采动裂隙层高度；c为底板
承压水导升裂隙层高度；
[0137]
s12：考虑到承压水导升裂隙的存在，承压水将沿着导升裂隙发育至有效隔水层底部，根据公式(2)计算底板承压水导升裂隙层顶部的等效水压力q2；根据公式(3)计算采空区冒落矸石的等效荷载q1；
[0138]
q2＝q0‑
(h1 h0)
·
γ1ꢀꢀꢀ
(2)；
[0139][0140]
式中，q0为高承压水实际作用的水压力；γ1为底板岩层的等效重度，单位kn/m3；l为采空区沿走向长度；γ为采空区冒落矸石的等效重度，单位kn/m3；
[0141]
s2：引入弹塑性破坏理论进行分析；具体方法如下：
[0142]
s20：隔水层中实际受到承压水压力作用是有效隔水层和底板采动裂隙层，可将有效隔水层和底板采动裂隙层简化成理想的线性弹塑性体；隔水岩层中实际受到承压水压力作用是有效隔水岩层和底板采动裂隙层，可将有效隔水岩层和底板采动裂隙层简化成理想的线性弹塑性体。有效隔水层变形符合虎克定律，屈服时满足mohr
‑
coulomb准则，有效隔水层的变形满足公式(4)；底板采动裂隙层处于塑性破坏阶段，其强度值降低至残余强度，屈服条件仍满足mohr
‑
coulomb准则，使底板采动裂隙层的变形满足公式(5)；
[0143]
σ
x
＝k
·
σ
y
σ
n
ꢀꢀꢀ
(4)；
[0144][0145]
式中，σ
x
为底板岩层中微单元体在x方向的水平应力，单位mpa；σ
y
为底板岩层中微单元体在y方向的垂直应力，单位mpa；σ
n
为底板岩层单轴抗压强度，单位mpa，或者通过实测获得；c
n
为底板岩层的粘聚力，单位mpa；为底板岩层的内摩擦角，单位
°
；
[0146]
为底板岩层残余强度，单位mpa，或者通过实测获得；为底板岩层的残余粘聚力，单位mpa；为底板岩层的残余内摩擦角，单位
°
；当底板岩层在弹性和残余阶段的内摩擦角不变时，即则有k
*
＝k；
[0147]
s3：选取微分单元，建立微分平衡方程；具体方法如下：
[0148]
从底板隔离层中选取一段厚度为dy的微单元体，建立底板隔水层微单元体受力分析，如图2所示，以图2的o点为原点，水平方向为x轴、垂直向下为y轴建立底板岩层微单元体进行受力分析，其中dσ
y
、dσ
x
分别为微单元体在垂直和水平方向的主应力；要使得底板隔水层保持稳定，则微单元在y方向上的微单元体应力的平衡方程满足根据mohr
‑
coulomb准则，即满足公式(6)；并根据公式(7)计算微单元体与两侧岩层接触面上的摩擦阻力τ；根据公式(6)建立y方向应力平衡方程得到公式(8)；
[0149]
∑f＝0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)；
[0150]
τ＝f
·
σ
x
c
ꢀꢀꢀꢀꢀꢀꢀ
(7)；
[0151]
σ
y
·
l 2τ
·
dy
‑
(σ
y
dσ
y
)
·
l＝0
ꢀꢀꢀꢀ
(8)；
[0152]
式中，c为微单元体与两侧岩层接触面上的粘聚力，单位mpa，当微单元体处于底板采动裂隙层h0时，当微单元体处于有效隔水层h1时，c＝c
n
；f为微单元体与两侧岩层接触面上的摩擦系数，当微单元体处于底板采动裂隙层h0时该参数为f
*
，当微单元体处于有效隔水层h1时为f，
[0153]
s4：解微分方程，推导底板岩层所能承受的最大极限水压力值或最小隔水岩层厚度安全阈值；
[0154]
当微单元体处于底板采动裂隙层时，按如下方法进行求解：
[0155]
s41：联立公式(6)、公式(7)和公式(8)获得公式(9)，并求解公式(9)得到公式(10)；
[0156][0157][0158]
式中，为未知量；
[0159]
s42：先使y＝0时，得到公式(11)；并联立公式(3)和公式(11)得到公式(12)；再推导出公式(13)；
[0160][0161][0162][0163]
s42：使y＝h0，得到公式(14)；
[0164][0165]
当微单元体处于有效隔水层时，按如下方法进行求解：
[0166]
s43：联立公式(5)、公式(7)和公式(8)获得公式(15)，并求解公式(15)得到公式(16)；
[0167][0168][0169]
式中，c1为未知量；
[0170]
s44：先使y＝h0 h1，得到公式(17)；并联立公式(2)和公式(17)得到公式(18)；再联立公式(16)和公式(18)得到(19)；
[0171][0172][0173][0174]
s45：使y＝h0，得到公式(20)；在y＝h0交界面处，公式(20)等于公式(14)，进而得到公式(21)，再求解公式(21)得到公式(22)；
[0175][0176][0177][0178]
s5：根据实际问题，判断隔水层厚度是否安全，具体过程如下：
[0179]
当底板隔水岩层厚度h1已知时，通过公式(22)给出底板岩层所能承受的极限水压力值，实际水压力q与q0存在公式(23)中的关系：当底板隔水岩层受水压力值已知时，对公式(22)求解得到底板岩层最小厚度值h
min
，底板岩层实际厚度值h1与最小厚度值h
min
存在公式(24)中的关系：
[0180]
[0181][0182]
步骤二：以双岩性结构底板隔水岩层进行稳定性分析，具体方法如下：
[0183]
a1：建立力学模型，进行受力分析；具体方法如下：
[0184]
a10：如图3所示，将岩煤开采层至承压开采层之间的部分简化为两种底板岩性结构的底板岩层，再将两种底板岩性结构的底板岩层由上到下分别记为底板隔水岩层a和底板隔水岩层b，并将底板隔水岩层a由上到下依次分为底板采动裂隙层和上部有效隔水层，将底板隔水岩层b由上到下依次分为下部有效隔水层和承压水导升裂隙层，其中上部有效隔水层和下部有效隔水层形成总有效隔水层；
[0185]
a11：根据公式(25)计算底板的总有效隔水层高度；
[0186]
h2 h3＝h2‑
h0‑
c
ꢀꢀꢀꢀꢀꢀ
(25)；
[0187]
式中，h2为底板岩层的总厚度，单位m，当底板隔水层实际厚度确定来求解承压水压力最大值时，h2＝h
实
；h0为底板采动裂隙层高度；c为底板承压水导升裂隙层高度；
[0188]
a12：根据公式(26)计算承压水导升裂隙层顶部的等效水压力q3；根据公式(27)计算采空区冒落矸石的等效荷载q1；
[0189]
q3＝q0‑
(h2 h0)
·
γ2‑
h3·
γ3ꢀꢀꢀ
(26)；
[0190][0191]
式中，q0为高承压水实际作用的水压力；h2为上部有效隔水层的高度；h0为底板采动裂隙层高度；h3为下部有效隔水层的高度；γ2为底板隔水岩层a的等效重度，单位kn/m3；γ3为底板隔水岩层b的等效重度，单位kn/m3；l为采空区沿走向长度；γ为采空区冒落矸石的等效重度，单位kn/m3；
[0192]
a2：引入弹塑性破坏理论进行分析；具体方法如下：
[0193]
a20：底板隔水层中实际受到承压水压力作用于有效隔水岩层和底板采动裂隙层，可将总有效隔水层和底板采动裂隙层简化成理想的线性弹塑性体，总有效隔水岩层变形符合虎克定律，屈服时满足mohr
‑
coulomb准则，底板隔水岩层a中的上部有效隔水层的变形满足公式(27)，使底板隔水岩层b中的下部有效隔水层的变形满足公式(28)；底板采动裂隙层处于塑性破坏阶段，其强度值降低至残余强度，屈服条件仍满足mohr
‑
coulomb准则，使底板采动裂隙层的变形满足公式(29)；
[0194]
σ
2x
＝k2·
σ
2y
σ
2n
ꢀꢀꢀ
(27)；
[0195]
σ
3x
＝k3·
σ
3y
σ
3n
ꢀꢀꢀꢀ
(28)；
[0196][0197]
式中，σ
2x
为底板隔水岩层a中微单元体在x方向的水平应力，单位mpa；σ
2y
为底板隔水岩层a中微单元体在y方向的垂直应力，单位mpa；σ
2n
为底板隔水岩层a的单轴抗压强度，单位mpa；或者通过实测获得；c
2n
为底板隔水岩层a的粘聚
力，单位mpa；为底板隔水岩层a的内摩擦角，单位
°
；
[0198]
σ
3x
为底板隔水岩层b中微单元体在x方向的水平应力，单位mpa；σ
3y
为底板隔水岩层b中微单元体在y方向的垂直应力，单位mpa；σ
3n
为底板隔水岩层b的单轴抗压强度，单位mpa；或者通过实测获得；σ
3n
为底板隔水岩层b的粘聚力，单位mpa；为底板隔水岩层b的内摩擦角，单位
°
；
[0199]
为底板隔水层顶部岩层a的残余强度，单位mpa，或者通过实测获得；为隔水层顶部岩层a的残余粘聚力，单位mpa；为隔水层顶部岩层a的残余内摩擦角，单位
°
；隔水层顶部岩层a在弹性和残余阶段的内摩擦角不变，即则有k
*
＝k；
[0200]
a3：选取微分单元，建立微分平衡方程；具体方法如下：
[0201]
从底板隔离层中选取一段厚度为dy的微单元体，建立底板隔水层微单元体受力分析，如图4所示，以水平方向为x轴、垂直向下为y轴建立底板隔水层微单元体进行受力分析，其中dσ
y
、dσ
x
分别为微单元体在垂直和水平方向的主应力；使微单元体在y方向上的微单元体应力的平衡方程满足公式(30)；要使得底板隔水层保持稳定，则微单元在y方向上的微单元体应力的平衡方程满足根据mohr
‑
coulomb准则，根据公式(31)计算微单元体与两侧岩层接触面上的摩擦阻力τ，根据公式(31)建立y方向应力平衡方程得到公式(32)；
[0202]
∑f＝0
ꢀꢀꢀ
(30)；
[0203]
τ＝f
·
σ
x
c
ꢀꢀꢀꢀꢀꢀ
(31)；
[0204]
σ
y
·
l 2τ
·
dy
‑
(σ
y
dσ
y
)
·
l＝0
ꢀꢀꢀ
(32)；
[0205]
式中，c为微单元体与两侧岩层接触面上的粘聚力，单位mpa，当微单元体处于底板采动裂隙层h0时，当微单元体处于上部有效隔水层h2时，c＝c
2n
；当微单元体处于下部有效隔水层下部时，c＝c
3n
；为微单元体与两侧岩层接触面上的摩擦系数，当微单元体处于底板采动裂隙层h0时该参数为f
*
，当微单元体处于上部有效隔水层h2时时为f2，当微单元体处于上部有效隔水层h3时为f3，
[0206]
a4：解微分方程，推导底板岩层所能承受的最大极限水压力值或最小隔水岩层厚度安全阈值；
[0207]
当微单元体处于底板采动裂隙层时，按如下方法进行求解：
[0208]
a41：联立公式(29)、公式(31)和公式(32)获得公式(33)，并求解公式(33)得到公式(34)；
[0209]
[0210][0211]
式中，为未知量；
[0212]
a42：先使y＝0时，得到公式(35)；并联立公式(27)和公式(35)得到公式(36)；再推导出公式(37)；
[0213][0214][0215][0216]
a42：使y＝h0，得到公式(38)；
[0217][0218]
当微单元体处于上部有效隔水层h2时，按如下方法进行求解：
[0219]
a43：联立公式(27)、公式(31)和公式(32)获得公式(39)，并求解公式(39)得到公式(40)；
[0220][0221][0222]
式中，c2为未知量；
[0223]
a44：使y＝h0，得到公式(41)；在y＝h0交界面处，公式(38)等于公式(41)，进而得到公式(42)，再联立公式(42)与公式(40)得到公式(43)；
[0224][0225]
[0226][0227]
a45：把y＝h0 h2代入公式(43)得到公式(44)；
[0228][0229]
当微单元体处于下部有效隔水层h3时，按如下方法进行求解：
[0230]
a46：把公式(28)和公式(31)代入公式(32)得到公式(45)；并求解公式(45)得到公式(46)；
[0231][0232][0233]
式中，c3为未知量；
[0234]
a47：使y＝h0 h2 h3得到公式(47)；并联立公式(24)和公式(47)得到公式(48)，联立公式(28)和公式(46)得到公式(49)；
[0235][0236][0237][0238]
a48：在y＝h0 h2交界面处，公式(44)等于公式(49)，进而得到公式(50)，再求解公式(50)得到公式(51)；
[0239][0240][0241]
a5：根据实际问题，判断隔水层厚度是否安全，具体过程如下：
[0242]
当底板隔水岩层的总厚度h2(h2＝h0 h2 h3 c)已知时，通过公式(51)给出底板隔水层所能承受的极限水压力值，实际水压力q与q0存在公式(52)中的关系：当底板隔水岩层受水压力值已知时，对公式(51)求解得到底板隔水层最小厚度值h
min
，底板隔水层实际厚度值h2与最小厚度值h
min
存在公式(53)中的关系：
[0243][0244][0245]
在步骤一中，底板采动裂隙层和承压水导升裂隙层的高度由现场实测或相关理论计算得出。
[0246]
在步骤二中，底板采动裂隙层和承压水导升裂隙层的高度由现场实测或相关理论计算得出。
[0247]
本发明中，首先在高等岩石力学理论基础上，建立了单一结构底板隔水岩层与双岩性结构底板隔水层力学模型，求解得出关键层位(有效隔水层)的底板隔水岩层解析式；再选取微分单元，并结合mohr
‑
coulomb准则和塑性破坏强度理论，建立了采动应力和承压水压力耦合作用的底板微元微分平衡方程，推导得到底板隔水岩层最大极限水压力和最小隔水层厚度理论解，能推导得出底板隔水岩层所能承受的最大极限水压力与最小隔水层厚度安全阈值，可科学合理的判断底板整体安全性，为高承压水体上采煤提供安全保障和理论指导，为煤矿安全开采提供了切实可行的理论指导。同时，结合煤岩层底板实际厚度(h
实
)、承压水压力(p
实
)、底板破坏深度、承压水导升高度、岩石容重γ等已知参数，求得底板隔水层所能承受的极限水压力值p0，或利用科学计算软件求解底板隔水层最小厚度值hmin，若保证现场不发生突水灾害，需满足p实≤pmax或h实≥hmin基本条件。
[0248]
实施例：
[0249]
结合华北某矿承压开采工作面底板各岩层厚度及物理力学参数(表1)，将底板隔水层整体划分为较符合实际的两层进行分析，上层厚度56m，主要由灰岩和砂质泥岩相互交错而成；下层厚度8m，主要岩性为铝质泥岩。根据现场地质情况计算每层岩体物理力学参数，求解底板由两种不同岩性组成的情况下所能承受的最大极限水压力，判断当前底板隔水层稳定性及突水危险性，参数取值如表1所示。
[0250]
表1复合底板基本参数
[0251][0252]
对于承压水导升裂隙带高度取c＝1m。将参数代入式(52)，得到双岩性结构底板极限水压力随工作面推进度变化曲线(图5)。从图5中可以看出：在工作面开采初期，底板所能承受的极限水压力q0与工作面推进长度l呈反比关系，底板所能承受的极限水压力随工作面的推进逐渐降低，这是因为煤岩层被采出后底板岩层应力平衡状态发生改变，底板岩层上部荷载的消失使得其所能承受的极限水压力减小。但在l＝0～50m范围内，底板所能承受的极限水压力q0均大于8.5mpa，大于《规程》中突水系数所对应的水压力值，说明在开采初期，采空区底板裸露面积较小，从结构力学角度来看，有效隔水层岩梁弯矩较小，所能承受的水压力较大，所以在水压力一定的情况下，开采初期底板隔水层结构稳定性较强，最不易发生突水事故。
[0253]
取该矿承压开采工作面的开采长度l为155.9m，确定了有效隔水层所能承受的最大极限水压力q0＝3.46mpa。结合该矿现场工程地质条件，该工作面下部承压水压力为1.8mpa，小于底板所能承受的极限水压力3.46mpa(图5)，根据判别公式(53)可知底板突水“底板隔水岩层稳定性好、不突水”。因此在工作面回采过程中，若不遇较大构造情况下，底板能够保持安全稳定。