一种耕地下采煤区域变形控制的地表沉叠加陷计算方法技术

本发明专利技术公开了一种耕地下采煤区域变形控制的地表沉陷叠加计算方法，属于地表沉陷预测领域。将地表沉陷问题简化为采空区下沉的高度，之后根据采空区间隔留下的垮落区空间以及充填体空间计算填满垮落区后的地表沉陷高度，由于充填体在受力过程中会出现压缩，因此根据充填体的充实率判断充填体在受力后的压缩空间，从而计算出充填体压实空间产生的地表沉陷高度，将前后两个地面沉陷高度相加即可获得地下开采空间产生的地表沉陷的总高度。其计算步骤简单，更符合实际情况，有效提高了区域岩层控制下地表移动和变形计算结果的精度，对实现煤粮复合区煤炭绿色开采和耕地协同修复治理有很强的指导作用。有很强的指导作用。有很强的指导作用。

【技术实现步骤摘要】
一种耕地下采煤区域变形控制的地表沉叠加陷计算方法


[0001]本专利技术涉及一种耕地下采煤区域变形控制的地表沉叠加陷计算方法，适用于煤粮复合区耕地下采煤设计，属于地表沉陷预测领域。
技术背景
[0002]耕地是国家粮食安全的重要基础，煤炭是我国能源战略安全的根本保证，煤粮复合区采煤沉陷耕地损毁与保护矛盾日趋凸显。长期以来，我国煤粮复合区煤炭开发主要采用长壁开采全部垮落法管理顶板，覆岩破坏和地表沉陷剧烈，对耕地破坏尤为严重。目前，虽然充填开采技术应用较为广泛，但全部充填成本高、效益低，而耕地的抗变形能力强，所以部分带状充填的区域性岩层控制必将成为未来耕地保护的主流方向。因此，如何计算煤粮复合区源头控制与耕地保护相结合的区域性地表移动与变形是亟待解决的瓶颈难题之一。
[0003]煤粮复合区采煤耕地协同保护的前提是需要实现区域性地表移动与变形的精准计算，保障其变形极值不超过耕地防控指标。目前，我国进行采煤地表沉陷和移动变形预计应用最广泛的方法是概率积分法，此外，在充填开采领域用于地表移动变形计算的等效采高理论经过不断发展研究，成为计算充填开采地表移动变形主流理论，但面向煤粮复合区耕地下采煤区域变形控制的地表沉陷计算方法尚缺乏相应的计算方法。

技术实现思路

[0004]针对上述技术中的问题，提供一种耕地下采煤区域变形控制的地表沉叠加陷计算方法，计算步骤简单，更符合实际情况，实现了如何计算煤粮复合区耕地下煤炭开采区域变形控制。
[0005]为实现上述技术目的，本专利技术的一种耕地下采煤区域变形控制的地表沉陷叠加计算方法，将地表沉陷问题简化为采空区下沉的高度，之后根据采空区间隔留下的垮落区空间以及充填体空间计算填满垮落区后的地表沉陷高度，由于充填体在受力过程中会出现压缩，因此根据充填体的充实率判断充填体在受力后的压缩空间，从而计算出充填体压实空间产生的地表沉陷高度，将前后两个地面沉陷高度相加即可获得地下开采空间产生的地表沉陷的总高度；其中实际采空区充填体越接近采空区中心区域压缩程度越高，越接近采空区两侧则压缩程度越低，因此在采用相同的充填材料和充填方法的情况下则认为所有采空区充填体的充实率相同，即压缩程度是相同的。
[0006]具体步骤如下：
[0007]针对煤粮复合区进行耕地下采煤区域，收集研究采煤区域地质资料、地层分布情况、钻孔柱状图等地质采矿资料，确定待开采煤层的开采深度H、开采厚度m、煤层倾角α信息；
[0008]确定该区域开采方法使用间隔垮落与回填交替设置的回填开采方式进行，然后根据采动程度和下沉系数的非线性关系，确定垮落工作面开采宽度为b，数量为n；充填开采工
作面宽度为a，数量为n
‑
1，填开采工作面的充实率为
[0009]基于等效采高理论计算充填工作面的等效采高等效采高m1为顶板压缩回填后的开采工作面的高度，为压缩后回填开采工作面充填体的高度；
[0010]将垮落工作面的原始开采范围为定义为开采范围R2，剖面范围R2等于采厚m乘以开采宽度b，将充填开采工作面顶部被下沉的顶板压缩空间的剖面定义为开采范围R1，开采范围R1等于等效采高m1乘以充填开采工作面宽度a；
[0011]将面向煤粮复合区耕地保护的地表变形区域控制的剖面开采范围R看作是n
‑
1长度为a、高度为m1的开采范围R1加上n个长度为b、高度为m的开采范围R2；
[0012]利用概率积分模型计算出开采范围R1引起的地表下沉W1、地表水平移动U1、地表倾斜I1、地表水平变形ε1和地表曲率K1，并且计算出开采范围R2引起的地表下沉W2、地表水平移动U2、地表倾斜I2、地表水平变形ε2和地表曲率K2；
[0013]最后计算出地表变形区域控制的地表下沉W＝W1+W2、地表水平移动U＝U1+U2、地表倾斜I＝I1+I2、地表水平变形ε＝ε1+ε2和地表曲率K＝K1+K2，即可准确预测出耕地下采煤区域变形控制的地表沉陷。
[0014]进一步，根据矿区已有的工作面开采沉陷实测参数确定地表沉陷概率积分模型预计参数，包括下沉系数q、水平移动系数B、主要影响角正切tanβ，开采影响传播角θ，拐点偏移距S，开采范围R1中单个垮落工作面沉陷预计参数与本矿区相同面宽工作面垮落法开采条件下的实测参数相同，根据这些参数和前文已经确定的开采深度H、开采厚度m、煤层倾角α即可进行地表沉陷计算。
[0015]进一步，开采范围R2中单个充填工作面沉陷预计参数参照《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采指南》，除下沉系数外其余预计参数采用本矿区薄煤层垮落法开采条件下实测参数，下沉系数应增大5％～10％。
[0016]有益效果：本专利技术以煤粮复合区耕地保护为前提，兼顾煤炭高效绿色开采，创造性的提出了一种耕地下采煤区域变形控制的地表沉陷叠加计算方法，解决了煤粮复合区耕地下采煤区域变形控制设计时地表移动变形如何计算的难题，其计算步骤简单，更符合实际情况，从科学合理的角度实现了地表沉陷的叠加计算，为煤粮复合区耕地下煤炭绿色开采设计提供了技术支撑，对实现煤粮复合区煤炭绿色开采和耕地保护协调可持续发展等具有重要的实际意义与应用价值。
附图说明
[0017]图1为本专利技术的耕地下采煤区域变形控制的地表沉陷叠加计算示意图。
具体实施方式
[0018]下面将结合图和具体实施过程对本专利技术做进一步详细说明：
[0019]如图1所示，本专利技术的一种耕地下采煤区域变形控制的地表沉陷叠加计算方法，步骤为：
[0020]首先通过研究区域地质采矿资料确定煤层的采深H、采厚m、倾角α等信息；然后根据采动程度和下沉系数的非线性关系，确定垮落工作面开采宽度b；之后设计充填开采工作面宽度a、充实率以及垮落工作面的个数n，接着可以确定充填工作面的个数为n
‑
1，并基于
等效采高理论计算充填工作面的等效采高这样面向煤粮复合区耕地保护的地表变形区域控制的剖面开采范围R看作是n
‑
1长度为a、高度为m1的开采范围R1加上n个长度为b、高度为m的开采范围R2。然后可以利用概率积分模型计算出开采范围R1引起的地表下沉W1、地表水平移动U1、地表倾斜I1、地表水平变形ε1和地表曲率K1，同样可以计算出开采范围R2引起的地表下沉W2、地表水平移动U2、地表倾斜I2、地表水平变形ε2和地表曲率K2，这样就可以计算出地表变形区域控制的地表下沉W＝W1+W2、地表水平移动U＝U1+U2、地表倾斜I＝I1+I2、地表水平变形ε＝ε1+ε2和地表曲率K＝K1+K2，进而实现了耕地下采煤区域变形控制的地表沉陷预测。
[0021]具体步骤为：
[0022]步骤1：针对煤粮复合区进行耕地下采煤区域，收集研究采煤区域地质资料、地层分布情况、钻孔柱状图等地质采矿资料，确定待开采煤层的开采深度H、开采厚度m、煤层倾角α等信息。
[0023]步骤2：根据矿区现有的工作面开采地表沉陷资料，建立采动程度和下沉系数的非线性关系，基于此可以确定本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种耕地下采煤区域变形控制的地表沉陷叠加计算方法，其特征在于：将地表沉陷问题简化为采空区下沉的高度，之后根据采空区间隔留下的垮落区空间以及充填体空间计算填满垮落区后的地表沉陷高度，由于充填体在受力过程中会出现压缩，因此根据充填体的充实率判断充填体在受力后的压缩空间，从而计算出充填体压实空间产生的地表沉陷高度，将前后两个地面沉陷高度相加即可获得地下开采空间产生的地表沉陷的总高度；其中实际采空区充填体越接近采空区中心区域压缩程度越高，越接近采空区两侧则压缩程度越低，因此在采用相同的充填材料和充填方法的情况下则认为所有采空区充填体的充实率相同，即压缩程度是相同的。2.根据权利要求1所述耕地下采煤区域变形控制的地表沉陷叠加计算方法，其特征在于具体步骤如下：针对煤粮复合区进行耕地下采煤区域，收集研究采煤区域地质资料、地层分布情况、钻孔柱状图等地质采矿资料，确定待开采煤层的开采深度H、开采厚度m、煤层倾角α信息；确定该区域开采方法使用间隔垮落与回填交替设置的回填开采方式进行，然后根据采动程度和下沉系数的非线性关系，确定垮落工作面开采宽度为b，数量为n；充填开采工作面宽度为a，数量为n
‑
1，填开采工作面的充实率为基于等效采高理论计算充填工作面的等效采高等效采高m1为顶板压缩回填后的开采工作面的高度，为压缩后回填开采工作面充填体的高度；将垮落工作面的原始开采范围为定义为开采范围R2，剖面范围R2等于采厚m乘以开采宽度b，将充填开采工作面顶部被下沉的顶板压缩空间的剖面定义为开采范围R1，开采范围R1等于等效...

【专利技术属性】
技术研发人员：郭广礼，李怀展，查剑锋，宫亚强，蒋乾，刘文锴，胡青峰，殷和建，汪铁宁，张学伟，王跃宗，陈延康，张鑫，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：