发明名称 一种覆岩导水裂缝带高度的确定方法
摘要 本发明公开了一种覆岩导水裂缝带高度的确定方法,包括步骤:一、根据采煤工作面上覆各岩层的层厚、容重、弹性模量和抗拉强度,判断采煤工作面上覆各岩层中是否存在主关键层和/或亚关键层,并确定主关键层和/或亚关键层的位置;二、根据基载比与基采比初步判断主关键层在覆岩“三带”中的位置,初步判定覆岩导水裂缝带的高度范围;三、分三种情况计算采煤工作面上覆各岩层的总拉伸变形量,并比较采煤工作面上覆各岩层的总拉伸变形量与采煤工作面上覆岩层的极限拉伸变形量,判断采煤工作面上覆各岩层是否发生了破断;四、分三种情况确定覆岩导水裂缝带的高度。本发明方法步骤简单,简化了确定覆岩导水裂缝带高度的方法,能够获得较为精确的结果。
申请公布号 CN104933266B 申请公布日期 2016.06.15
申请号 CN201510400983.2 申请日期 2015.07.09
申请人 西安科技大学 发明人 赵兵朝;刘樟荣;同超;王守印;刘宾;刘飞;王春龙
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西安创知专利事务所 61213 代理人 谭文琰
主权项 一种覆岩导水裂缝带高度的确定方法,其特征在于该方法包括以下步骤:步骤一、根据采煤工作面上覆各岩层的层厚、容重、弹性模量和抗拉强度,判断采煤工作面上覆各岩层中是否存在主关键层和/或亚关键层,并确定主关键层和/或亚关键层的位置;步骤二、根据基载比J<sub>z</sub>与基采比J<sub>c</sub>初步判断主关键层在覆岩“三带”中的位置,初步判定覆岩导水裂缝带的高度范围;具体方法如下:当J<sub>z</sub><J<sub>z0</sub>且J<sub>c</sub><J<sub>c1</sub>时,主关键层位于“冒落带”,此时覆岩导水裂缝带的高度等于煤层的埋深,即确定出了覆岩导水裂缝带的高度;当J<sub>z</sub><J<sub>z0</sub>且J<sub>c1</sub><J<sub>c</sub><J<sub>c2</sub>时,主关键层位于“裂缝带”,此时覆岩导水裂缝带的高度大于主关键层与煤层的距离,执行步骤三;当J<sub>z</sub>&gt;J<sub>z0</sub>且J<sub>c</sub>&lt;J<sub>c2</sub>时,主关键层位于“裂缝带”,此时覆岩导水裂缝带的高度大于主关键层与煤层的距离,执行步骤三;当J<sub>z</sub><J<sub>z0</sub>且J<sub>c</sub>&gt;J<sub>c2</sub>时,主关键层位于“弯曲下沉带”,此时覆岩导水裂缝带的高度大于主关键层与煤层的距离,执行步骤三;当J<sub>z</sub>&gt;J<sub>z0</sub>且J<sub>c</sub>&gt;J<sub>c2</sub>时,主关键层位于“弯曲下沉带”,覆岩导水裂缝带的高度小于或等于主关键层与煤层的距离,执行步骤三;其中,<img file="FDA0000904704650000011.GIF" wi="246" he="134" />H<sub>基岩</sub>为基岩厚度,H<sub>载荷层</sub>为载荷层厚度;<img file="FDA0000904704650000012.GIF" wi="267" he="140" />m为煤层采高;J<sub>z0</sub>为基载比的临界值,J<sub>c1</sub>为基采比的下临界值,J<sub>c2</sub>为基采比的上临界值;步骤三、分三种情况计算采煤工作面上覆各岩层的总拉伸变形量,并比较采煤工作面上覆各岩层的总拉伸变形量与采煤工作面上覆岩层的极限拉伸变形量[ε],判断采煤工作面上覆各岩层是否发生了破断;具体方法为:情况一、当主关键层位于“裂缝带”时,此时采煤工作面开采尺寸大于等于主关键层的极限破断距,主关键层破断,覆岩导水裂缝带的高度大于主关键层的高度,根据公式<img file="FDA0000904704650000021.GIF" wi="363" he="135" />依次从下往上计算主关键层及其以上各岩层的总拉伸变形量ε,每计算一层,就比较一次ε与[ε],当ε&lt;[ε]时,判定为岩层未发生破断,当ε≥[ε]时,判定为岩层发生了破断,直到计算到相邻的两层岩层,下一层发生了破断,而上一层未发生破断后停止计算;情况二、当主关键层位于“弯曲下沉带”时,此时采煤工作面开采尺寸小于等于主关键层的极限破断距,主关键层未破断,覆岩导水裂缝带的高度小于主关键层的高度,当此时同时满足采煤工作面开采尺寸大于亚关键层的极限破断距时,亚关键层破断,覆岩导水裂缝带的高度介于亚关键层高度和主关键层高度之间,根据公式<img file="FDA0000904704650000022.GIF" wi="365" he="141" />依次从下往上计算亚关键层及其以上各岩层的总拉伸变形量ε,每计算一层,就比较一次ε与[ε],当ε&lt;[ε]时,判定为岩层未发生破断,当ε≥[ε]时,判定为岩层发生了破断,直到计算到相邻的两层岩层,下一层发生了破断,而上一层未发生破断后停止计算;情况三、当主关键层位于“弯曲下沉带”时,此时采煤工作面开采尺寸小于等于主关键层的极限破断距,主关键层未破断,覆岩导水裂缝带的高度小于主关键层的高度,当此时同时满足采煤工作面开采尺寸小于亚关键层的极限破断距时,亚关键层未破断,覆岩导水裂缝带的高度小于亚关键层的高度,根据公式<img file="FDA0000904704650000023.GIF" wi="364" he="136" />依次从上往下计算亚关键层及其以下各岩层的总拉伸变形量ε,每计算一层,就比较一次ε与[ε],当ε&lt;[ε]时,判定为岩层未发生破断,当ε≥[ε]时,判定为岩层发生了破断,直到计算到相邻的两层岩层,上一层发生了破断,而下一层未发生破断后停止计算;以上三种情况中,L<sub>0</sub>为岩层弯曲变形前的直线段长度且<img file="FDA0000904704650000024.GIF" wi="262" he="143" />H为岩层与煤层的垂直距离,β为岩层的移动角;L<sub>1</sub>为岩层弯曲变形后的曲线段弧长且<img file="FDA0000904704650000025.GIF" wi="437" he="119" />l<sub>a</sub>和l<sub>b</sub>的计算方法为:对采煤工作面上覆岩层下沉盆地的内边缘曲线和外边缘曲线分别采用两段椭圆曲线进行描述,两段椭圆曲线的短轴长度均为l<sub>a</sub>,且l<sub>a</sub>=mη;两段椭圆曲线的长轴长度均为l<sub>b</sub>,且<img file="FDA0000904704650000031.GIF" wi="349" he="142" />η为导水裂缝带发育高度达到最大时采煤工作面上覆岩层的下沉系数,且η=(m‑H(k<sub>s</sub>‑1)cosa)/m,m为煤层采高,k<sub>s</sub>为覆岩的残余碎胀系数,a为煤层倾角;步骤四、分三种情况确定覆岩导水裂缝带的高度,具体方法为:情况一、当主关键层位于“裂缝带”,且经步骤三判定主关键层上覆岩层第r层岩层破断而第r+1层岩层未破断时,覆岩导水裂缝带高度H<sub>f</sub>等于第r+1层岩层与煤层的距离H<sub>r+1</sub>;其中,r的取值为自然数;情况二、当主关键层位于“弯曲下沉带”,且采煤工作面开采尺寸大于亚关键层的极限破断距而小于主关键层的极限破断距,且经步骤三判定从亚关键层往上数第g层岩层破断而第g+1层岩层未破断时,覆岩导水裂缝带高度H<sub>f</sub>等于第g+1层岩层与煤层的距离H<sub>g+1</sub>;其中,g的取值为自然数;情况三、当主关键层位于“弯曲下沉带”,且采煤工作面开采尺寸小于亚关键层的极限破断距,且经步骤三判定从亚关键层往下数第j层岩层破断而第j+1层岩层未破断时,覆岩导水裂缝带高度H<sub>f</sub>等于第j+1层岩层与煤层的距离H<sub>j+1</sub>;其中,j的取值为自然数;步骤一中根据采煤工作面上覆各岩层的层厚、容重、弹性模量和抗拉强度,判断采煤工作面上覆各岩层中是否存在主关键层和/或亚关键层,并确定主关键层和/或亚关键层的位置的具体过程为:步骤101、根据公式(q<sub>n+1</sub>)<sub>1</sub>&lt;(q<sub>n</sub>)<sub>1</sub>判断采煤工作面上覆各岩层中是否存在坚硬岩层,当公式(q<sub>n+1</sub>)<sub>1</sub>&lt;(q<sub>n</sub>)<sub>1</sub>成立时,判定第n+1层岩层为坚硬岩层,并将第n+1层岩层作为第1层岩层,重新根据公式(q<sub>n+1</sub>)<sub>1</sub>&lt;(q<sub>n</sub>)<sub>1</sub>判断第n+1层岩层的上覆各岩层中是否存在坚硬岩层,重复以上过程,直至确定出采煤工作面上覆各岩层中的所有坚硬岩层位置,并将确定出的坚硬岩层自下而上依次编号为1、2、…、w;其中,(q<sub>n+1</sub>)<sub>1</sub>为从采煤工作面上覆第n+1层岩层往下看时第1层岩层所受的载荷,且<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mn>1</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>E</mi><mn>1</mn></msub><msubsup><mi>h</mi><mn>1</mn><mn>3</mn></msubsup><mrow><mo>(</mo><msub><mi>&gamma;</mi><mn>1</mn></msub><msub><mi>h</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&gamma;</mi><mn>2</mn></msub><msub><mi>h</mi><mn>2</mn></msub><mo>+</mo><mo>...</mo><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><msub><mi>h</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>E</mi><mn>1</mn></msub><msubsup><mi>h</mi><mn>1</mn><mn>3</mn></msubsup><mo>+</mo><msub><mi>E</mi><mn>2</mn></msub><msubsup><mi>h</mi><mn>2</mn><mn>3</mn></msubsup><mo>+</mo><mo>...</mo><mo>+</mo><msub><mi>E</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><msubsup><mi>h</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mn>3</mn></msubsup></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000904704650000032.GIF" wi="862" he="150" /></maths>(q<sub>n</sub>)<sub>1</sub>为从采煤工作面上覆第n层岩层往下看时第1层岩层所受的载荷,且<img file="FDA0000904704650000041.GIF" wi="812" he="166" />n、n+1为从煤层以上的第1层岩层从下往上数的岩层编号,且n的取值为自然数;E<sub>1</sub>、E<sub>2</sub>、…、E<sub>n+1</sub>分别为第1、2、…、n+1层岩层的弹性模量,h<sub>1</sub>、h<sub>2</sub>、…、h<sub>n+1</sub>分别为第1、2、…、n+1层岩层的厚度,γ<sub>1</sub>、γ<sub>2</sub>、…、γ<sub>n+1</sub>分别为第1、2、…、n+1层岩层的容重;w的取值为自然数;步骤102、根据公式<img file="FDA0000904704650000042.GIF" wi="301" he="158" />计算步骤101中确定出的各坚硬岩层中第k层坚硬岩层的破断距L<sub>k</sub>,并根据公式<img file="FDA0000904704650000043.GIF" wi="390" he="172" />计算步骤101中确定出的各坚硬岩层中第k+1层坚硬岩层的破断距L<sub>k+1</sub>;其中,h<sub>k</sub>、h<sub>k+1</sub>分别为第k、k+1层坚硬岩层的厚度,σ<sub>k</sub>、σ<sub>k+1</sub>分别为第k、k+1层坚硬岩层的抗拉强度,q<sub>k</sub>为第k层坚硬岩层与k+1层坚硬岩层之间所夹的岩层施加给第k层坚硬岩层的载荷且<img file="FDA0000904704650000044.GIF" wi="419" he="287" />E<sub>k</sub>为第k层坚硬岩层的弹性模量,h<sub>x</sub>为第k层坚硬岩层与k+1层坚硬岩层之间所夹的岩层中从第k层坚硬岩层往上数第x层岩层的厚度,γ<sub>x</sub>为第k层坚硬岩层与k+1层坚硬岩层之间所夹的岩层中从第k层坚硬岩层往上数第x层岩层的容重,E<sub>x</sub>为第k层坚硬岩层与k+1层坚硬岩层之间所夹的岩层中从第k层坚硬岩层往上数第x层岩层的弹性模量,s为第k层坚硬岩层与k+1层坚硬岩层之间所夹的岩层的层数,且s的取值为自然数,x为1~s的自然数;q<sub>k+1</sub>为第k+1层坚硬岩层与k+2层坚硬岩层之间所夹的岩层施加给第k+1层坚硬岩层的载荷且<img file="FDA0000904704650000045.GIF" wi="517" he="308" />E<sub>k+1</sub>为第k+1层坚硬岩层的弹性模量,h<sub>y</sub>为第k+1层坚硬岩层与k+2层坚硬岩层之间所夹的岩层中从第k+1层坚硬岩层往上数第y层岩层的厚度,γ<sub>y</sub>为第k+1层坚硬岩层与k+2层坚硬岩层之间所夹的岩层中从第k+1层坚硬岩层往上数第y层岩层的容重,E<sub>y</sub>为第k+1层坚硬岩层与k+2层坚硬岩层之间所夹的岩层中从第k+1层坚硬岩层往上数第y层岩层的弹性模量,s′为第k+1层坚硬岩层与k+2层坚硬岩层之间所夹的岩层的层数,且s′的取值为自然数,y为1~s′的自然数;步骤103、判断L<sub>k</sub>&lt;L<sub>k+1</sub>是否成立,当L<sub>k</sub>&lt;L<sub>k+1</sub>成立时,判定第k层坚硬岩层为主关键层,否则,当L<sub>k</sub>&lt;L<sub>k+1</sub>不成立时,根据公式<img file="FDA0000904704650000051.GIF" wi="302" he="169" />计算步骤101中确定出的各坚硬岩层中第k层坚硬岩层的破断距L′<sub>k</sub>,并根据公式<img file="FDA0000904704650000052.GIF" wi="407" he="167" />计算步骤101中确定出的各坚硬岩层中第k+2层坚硬岩层的破断距L<sub>k+2</sub>,然后判断L<sub>k</sub>&lt;L<sub>k+2</sub>是否成立,当L<sub>k</sub>&lt;L<sub>k+2</sub>成立时,判定第k层坚硬岩层为主关键层,否则,当L<sub>k</sub>&lt;L<sub>k+2</sub>不成立时,判定第k层坚硬岩层为亚关键层;其中,h<sub>k+2</sub>为第k+2层坚硬岩层的厚度,σ<sub>k+2</sub>为第k+2层坚硬岩层的抗拉强度,q′<sub>k</sub>为第k层坚硬岩层与k+2层坚硬岩层之间所夹的岩层施加给第k层坚硬岩层的载荷且<img file="FDA0000904704650000053.GIF" wi="461" he="295" />h<sub>λ</sub>为第k层坚硬岩层与k+2层坚硬岩层之间所夹的岩层中从第k层坚硬岩层往上数第λ层岩层的厚度,γ<sub>λ</sub>为第k层坚硬岩层与k+2层坚硬岩层之间所夹的岩层中从第k层坚硬岩层往上数第λ层岩层的容重,E<sub>λ</sub>为第k层坚硬岩层与k+2层坚硬岩层之间所夹的岩层中从第k层坚硬岩层往上数第λ层岩层的弹性模量,λ为1~(s+s′+1)的自然数;q<sub>k+2</sub>为第k+2层坚硬岩层与k+3层坚硬岩层之间所夹的岩层施加给第k+2层坚硬岩层的载荷且<img file="FDA0000904704650000054.GIF" wi="525" he="294" />E<sub>k+2</sub>为第k+2层坚硬岩层的弹性模量,h<sub>z</sub>为第k+2层坚硬岩层与k+3层坚硬岩层之间所夹的岩层中从第k+2层坚硬岩层往上数第z层岩层的厚度,γ<sub>z</sub>为第k+2层坚硬岩层与k+3层坚硬岩层之间所夹的岩层中从第k+2层坚硬岩层往上数第z层岩层的容重,E<sub>z</sub>为第k+2层坚硬岩层与k+3层坚硬岩层之间所夹的岩层中从第k+2层坚硬岩层往上数第z层岩层的弹性模量,s″为第k+2层坚硬岩层与k+3层坚硬岩层之间所夹的岩层的层数,且s″的取值为自然数,z为1~s″的自然数;重复步骤102和步骤103,直至判断完各坚硬岩层是否为主关键层或亚关键层,当将坚硬岩层判定为主关键层和/或亚关键层时,即确定了主关键层和/或亚关键层的位置。
地址 710054 陕西省西安市雁塔路中段58号