CN115510776A - 采空区覆岩渗流能力计算与修复方法、装置及电子设备 - Google Patents

采空区覆岩渗流能力计算与修复方法、装置及电子设备 Download PDF

Info

Publication number
CN115510776A
CN115510776A CN202211113009.4A CN202211113009A CN115510776A CN 115510776 A CN115510776 A CN 115510776A CN 202211113009 A CN202211113009 A CN 202211113009A CN 115510776 A CN115510776 A CN 115510776A
Authority
CN
China
Prior art keywords
fracture
zone
mining area
channel
overlying strata
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.)
Granted
Application number
CN202211113009.4A
Other languages
English (en)
Other versions
CN115510776B (zh
Inventor
李全生
郭俊廷
贺安民
滕腾
朱艳照
李鹏
李军
张村
阎跃观
李召龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology Beijing CUMTB
China Energy Investment Corp Ltd
National Institute of Clean and Low Carbon Energy
Original Assignee
China University of Mining and Technology Beijing CUMTB
China Energy Investment Corp Ltd
National Institute of Clean and Low Carbon Energy
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology Beijing CUMTB, China Energy Investment Corp Ltd, National Institute of Clean and Low Carbon Energy filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN202211113009.4A priority Critical patent/CN115510776B/zh
Publication of CN115510776A publication Critical patent/CN115510776A/zh
Application granted granted Critical
Publication of CN115510776B publication Critical patent/CN115510776B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21FSAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00Methods or devices for use in mines or tunnels, not covered elsewhere
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Computing Systems (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geology (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了采空区覆岩渗流能力计算与修复方法、装置及电子设备,其方法包括:首先得到研究矿区覆岩的三带结构,其次利用超声波测试法对研究矿区覆岩的三带结构进行测量并得到修复前的渗水能力,最后对研究矿区覆岩进行修复并计算渗水能力,评估修复后研究矿区覆岩渗透改善情况。本发明结合覆岩岩性硬度与超声波测试法进行裂隙测量并能计算得到覆岩渗流能力,利用高压水射流对主要竖向裂隙通道进行水力化碎胀处理,实现覆岩结构裂隙带人工修复与治理,为煤炭安全生产、含水层保护以及环境保护提供技术保障。

Description

采空区覆岩渗流能力计算与修复方法、装置及电子设备
技术领域
本发明涉及采空区覆岩结构渗流能力修复与治理领域,尤其涉及一种采空区覆岩渗流能力计算与修复方法、装置及电子设备。
背景技术
85%以上的煤矿采用井工开采,煤炭地下开采后工作面后方采空区上覆岩层由下向上依次发育形成“三带”结构,包含冒落带、裂隙带和弯曲下沉带。其中弯曲下沉带由于没有形成连通的裂隙或仅仅存在一些因离层产生的横向分布裂隙,对水体的流动影响不大,但冒落带和裂隙带中裂隙发育纵横交错,尤其是裂隙带中贯穿岩层的竖向裂隙,为地下水的跃层流动提供了主要通道,不仅对煤矿安全生产带来严重的安全隐患,同时造成地下水系统的破坏,对矿区生态环境和水土保持具有重要影响。特别是,一些矿区(比如中国西部矿区)处于干旱气候地区,矿区煤矿上覆岩的弱胶结特征显著,岩层多为砂岩、砂质泥岩等水理效应显著的岩性,例如,经研究表明,布尔台矿上覆岩层中砂质泥岩吸水饱和后抗压强度和抗拉强度下降60%以上,岩石塑性增大30%以上,裂隙发育所需消耗的能量增大50%以上,生态环境极其脆弱,煤炭资源开采与生态、水环境保护的协调发展问题十分突出。采空区覆岩结构中冒落带、裂隙带的竖向裂隙渗流能力对采空区覆岩稳定性具有重大影响,对采空区覆岩渗流能力计算及修复治理是亟需解决的技术难点,也对煤炭安全生产、含水层保护以及环境保护具有重要意义。
发明内容
针对背景技术所指出的技术问题,本发明的目的在于提供一种采空区覆岩渗流能力计算与修复方法、装置及电子设备,结合覆岩岩性硬度分类计算三带结构高度,通过超声波测试法进行裂隙测量并能计算得到覆岩渗流能力,利用高压水射流对主要竖向裂隙通道进行水力化碎胀处理,实现覆岩结构裂隙带人工修复与治理,为煤炭安全生产、含水层保护以及环境保护提供技术保障。
本发明的目的通过下述技术方案实现:
采空区覆岩渗流能力计算与修复方法,其方法包括:
A、根据覆岩岩性硬度依次划分坚硬、中硬、软弱、极软四个等级,将研究矿区覆岩从下至上依次划分为冒落带、裂隙带和弯曲下沉带的三带结构,根据研究矿区的地质参数得出工作面覆岩的冒落带高度、裂隙带高度、弯曲下沉带高度;
B、利用超声波测试法对研究矿区覆岩的三带结构进行测量,得到裂隙带竖向裂隙通道、弯曲下沉带的分布位置数据,裂隙带竖向裂隙通道的分布位置数据包括竖向裂隙通道位置数据与尺寸数据,弯曲下沉带的分布位置数据包括弯曲下沉带的位置、高度数据;从裂隙带竖向裂隙通道的分布位置数据中筛选出主要的竖向裂隙通道;
C、按照如下公式计算研究矿区覆岩的渗水能力K:
K=βD3/12,其中β为裂隙长度,D为主要竖向裂隙通道开度;
D、对筛选出的主要竖向裂隙通道按照如下方法进行修复:
D1、以筛选出的主要竖向裂隙通道作为目标裂隙通道,由地面向位于目标裂隙通道上方的弯曲下沉带、目标裂隙通道竖向钻孔,在孔中布设高压水管,高压水管位于目标裂隙通道的部分开有出水喷头;
D2、利用高压水管进行水压控制在目标裂隙通道中进行高压水射流处理,将目标裂隙通道的大裂隙通道打碎成为n个开度更小的裂缝;
E、计算修复后研究矿区覆岩的渗水能力K
Figure BDA0003844311000000021
1≤i≤n,其中Di表示裂缝i的开度,βi为裂隙i的长度;
F、评估修复后研究矿区覆岩渗透改善情况,渗透改善情况采用包括渗流能力降低量K、渗流能力降低倍率B两个指标,K=K-K
Figure BDA0003844311000000031
本发明步骤E可以按照如下方法进行替换:方法E采用如下方法计算修复后研究矿区覆岩的渗水能力K
E1、目标裂隙通道两壁岩土的吸水膨胀变形公式表示为Δυ=αΔx/(xl+x),Δυ表示岩体体积变化量,α表示吸水膨胀系数,xl表示吸水膨胀的含水率修正参数,x表示总的含水率,Δx表示吸水量.
E2、岩土体材料可近似看成为各向同性材料,其某一方向的线应变约为体应变的三分之一,具体的层状裂隙通道的开度减小系数定义为
Figure BDA0003844311000000032
E3、计算修复后研究矿区覆岩的渗水能力计算公式为:
Figure BDA0003844311000000033
考虑到岩体碎裂吸水膨胀,修复后研究矿区覆岩的渗水能力K
Figure BDA0003844311000000034
其中Di表示裂缝i的开度,βi为裂隙i的长度,λ为具体的层状裂隙通道的开度减小系数。
本发明提供全面针对覆岩岩性硬度为坚硬、中硬、软弱、极软四个等级的分类计算,以本发明主要针对的覆岩岩性硬度为软弱为优选方案并提供如下计算方法:
覆岩岩性为软弱的冒落带高度公式如下:
Hm=100∑M/(6.2∑M+32)±1.5,其中M表示采高,单位m;
覆岩岩性为软弱的裂隙带高度公式如下:
Hl=100∑M/(3.1∑M+5)±4.0,其中M表示采高,单位m;
覆岩岩性为软弱的弯曲下沉带高度等于采深减去裂隙带高度。
本发明优选的步骤B方法包括:
向研究矿区覆岩发出超声波,根据超声波在研究矿区覆岩的三带结构中传播速度、振幅的区别识别三带结构的具体带层及带层的高度;若超声波在研究矿区覆岩三带结构某个带层具有明显波速变化,则判断该带层为裂隙带,同时根据超声波的波速变化记录情况得出裂隙带竖向裂隙通道的分布位置数据。
本发明优选的步骤D1中钻孔包括竖直孔和L型孔,竖直孔中布设的高压水管配合为竖直形状的高压水管,L型孔中布设的高压水管配合为L型高压水管。
本发明优选的步骤D2中单个目标裂隙通道的高压水射流处理后打碎成为n 个开度为Di、长度为βi的裂缝,则该目标裂隙通道的渗流能力降低倍率
Figure BDA0003844311000000041
其中β为目标裂隙通道裂隙长度,D为目标裂隙通道开度,Di表示裂缝i的开度,βi为裂隙i的长度。
本发明优选的目标裂隙通道的渗流能力降低量
Figure BDA0003844311000000042
一种采空区覆岩渗流能力计算与修复装置,所述装置包括:
获取模块,用于获取研究矿区覆岩的裂隙带竖向裂隙通道、弯曲下沉带的分布位置数据,获取模块内部按照覆岩岩性硬度等级构建有冒落带高度、裂隙带高度、弯曲下沉带高度的公式数据库,获取模块还用于获取研究矿区覆岩所对应的覆岩岩性硬度;
修复系统,筛选出的主要竖向裂隙通道进行竖向钻孔,在孔中布设高压水管,利用高压水管进行水压控制在目标裂隙通道中进行高压水射流处理,将目标裂隙通道的大裂隙通道打碎成为n个开度更小的裂缝;
处理模块,用于计算研究矿区覆岩修复前后的渗水能力,并计算得到渗流能力降低量、渗流能力降低倍率评估修复后研究矿区覆岩渗透改善情况。
一种电子设备,包括:至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器执行采空区覆岩渗流能力计算与修复方法的步骤。
本发明较现有技术相比,具有以下优点及有益效果:
(1)本发明结合覆岩岩性硬度分类计算三带结构高度,通过超声波测试法进行裂隙测量并能计算得到覆岩渗流能力,利用高压水射流对主要竖向裂隙通道进行水力化碎胀处理,实现覆岩结构裂隙带人工修复与治理,为煤炭安全生产、含水层保护以及环境保护提供技术保障。
(2)本发明能够全面针对覆岩岩性硬度四个等级进行渗流能力计算与修复,尤其是对覆岩为弱胶结特征显著,岩层多为砂岩、砂质泥岩等水理效应显著的岩性,利用高压水射流技术对裂隙带和弯曲下沉带主要竖向通道裂隙进行水力化碎胀处理,有效阻抑导水裂隙带纵向传导过程。
附图说明
图1为实施例中研究矿区地质参数示意及修复高压水管布设示意图;
图2为实施例中岩体体积变化率与吸水率的关系图;
图3为本发明高压水射流处理的结构示意图;
图4为实施例中举例大裂隙通道打碎为小裂隙的对比示意图。
具体实施方式
下面结合实施例对本发明作进一步地详细说明:
实施例
如图1~图4所示,采空区覆岩渗流能力计算与修复方法,其方法包括:
A、根据覆岩岩性硬度依次划分坚硬、中硬、软弱、极软四个等级,将研究矿区覆岩从下至上依次划分为冒落带、裂隙带和弯曲下沉带的三带结构,根据研究矿区的地质参数(包括覆岩岩性硬度、工作面布置、开采参数)得出工作面覆岩的冒落带高度、裂隙带高度、弯曲下沉带高度,地质参数为研究矿区的覆岩岩性硬度并按照覆岩岩性硬度划分坚硬、中硬、软弱、极软的四个等级确定研究矿区的覆岩岩性硬度等级,研究矿区的工作面布置和开采参数如图1所示,研究矿区覆岩从下至上依次划分为冒落带3、裂隙带4和弯曲下沉带5的三带结构,采空区1大部分位于冒落带3下方位置,采空区1部分处于煤层2中,弯曲下沉带5上方可能具有离层7。经研究总结,覆岩层岩性按照覆岩岩性硬度呈现不同分布特性,为了应对所有覆岩岩性硬度的计算,本发明提供如下高度公式:覆岩岩性为坚硬的冒落带高度公式如下: Hm=100×∑M/(2.1×∑M+16)±2.5,其中M表示采高,覆岩岩性为坚硬的裂隙带高度公式如下:Hl=100∑M/(1.2∑M+2)±8.9,其中M表示采高,单位m;覆岩岩性为坚硬的弯曲下沉带高度等于采深减去裂隙带高度。覆岩岩性为中硬的冒落带高度公式如下:Hm=100∑M/(4.7∑M+19)±2.2,其中M表示采高,单位m;覆岩岩性为中硬的裂隙带高度公式如下: Hl=100∑M/(1.6∑M+3.6)±5.6,其中M表示采高,单位m;覆岩岩性为中硬的弯曲下沉带高度等于采深减去裂隙带高度。覆岩岩性为软弱的冒落带高度公式如下:Hm=100∑M/(6.2∑M+32)±1.5,其中M表示采高,单位m;覆岩岩性为软弱的裂隙带高度公式如下: Hl=100∑M/(3.1∑M+5)±4.0,其中M表示采高,单位m;覆岩岩性为软弱的弯曲下沉带高度等于采深减去裂隙带高度。覆岩岩性为极软弱的冒落带高度公式如下:Hm=100∑M/(7.0M+63)±1.2,其中M表示采高;覆岩岩性为极软弱的裂隙带高度公式如下: Hl=100∑M/(5M+8)±3.0,其中M表示采高;覆岩岩性为极软弱的弯曲下沉带高度等于采深减去裂隙带高度。本实施例选取中国西部矿区岩层为例,矿区岩层弱胶结特征显著,覆岩岩性为软弱,采高M选取8.6m,采深H为 200m。土层厚度范围为0~27m。根据表中公式计算得冒落带高度Hm范围为8.6~ 11.6m,裂隙带高度Hl范围为23.2~31.2m,弯曲下沉带高度范围为142~178m。
B、利用超声波测试法对研究矿区覆岩的三带结构进行测量(重点是针对三带结构的竖向裂隙通道进行分布测量并筛选出主要的竖向裂隙通道,在筛选出主要的竖向裂隙通道后,可以对冒落带、裂隙带和弯曲下沉带进行精准测量),得到裂隙带竖向裂隙通道、弯曲下沉带的分布位置数据,裂隙带竖向裂隙通道的分布位置数据包括竖向裂隙通道位置数据与尺寸数据,弯曲下沉带的分布位置数据包括弯曲下沉带的位置、高度数据;从裂隙带竖向裂隙通道的分布位置数据中筛选出主要的竖向裂隙通道。在一些实施例中,步骤B方法包括:
向研究矿区覆岩发出超声波,根据超声波在研究矿区覆岩的三带结构中传播速度、振幅的区别识别三带结构的具体带层及带层的高度;若超声波在研究矿区覆岩三带结构某个带层具有明显波速变化,则判断该带层为裂隙带,同时根据超声波的波速变化记录情况得出裂隙带竖向裂隙通道的分布位置数据。以中国西部矿区岩层为例,获得覆岩裂隙带主要竖向裂隙通道高度为27~33.2m,开度为10cm,弯曲下沉带离层高度为157.2~157.8m离层高度为60cm。本发明采用按照覆岩层岩性进行理论计算和重点区域实际探测相结合,对于三带位置及结构有较为全面的了解,为打钻孔和导水管的铺设提供保障,省时省力。
C、按照如下公式计算研究矿区覆岩的渗水能力K:
K=βD3/12,其中β为裂隙长度,D为主要竖向裂隙通道开度;
D、对筛选出的主要竖向裂隙通道按照如下方法进行修复:
D1、以筛选出的主要竖向裂隙通道作为目标裂隙通道,由地面向位于目标裂隙通道上方的弯曲下沉带、目标裂隙通道竖向钻孔,在孔中布设高压水管,高压水管位于目标裂隙通道的部分开有出水喷头。在一些实施例中,步骤D1中钻孔包括竖直孔和L型孔,竖直孔中布设的高压水管配合为竖直形状的高压水管,L型孔中布设的高压水管配合为L型高压水管,竖直形状的高压水管与L型高压水管布设情况如图1所示,图1中标号6表示竖直孔或L型孔,标号8表示高压水管。
D2、利用高压水管进行水压控制在目标裂隙通道中进行高压水射流处理,将目标裂隙通道的大裂隙通道打碎成为n个开度更小的裂缝。如图3所示,在目标裂隙通道中布设高压水管8,高压水管8上的高压水喷头10置于目标裂隙通道中,高压水通过高压水管8传至目标裂隙通道内部,由高压水喷头10喷射而出,射入目标裂隙通道的侧壁岩层11,目标裂隙通道内部岩石碎胀12,将大裂隙通道打碎成为开度更小的裂缝13。
E、计算修复后研究矿区覆岩的渗水能力K
Figure BDA0003844311000000081
1≤i≤n,其中Di表示裂缝i的开度,βi为裂隙i的长度。简单举例如下:若步骤D2中将目标裂隙通道的大裂隙通道打碎成为n个开度为D/n(如图4打碎对比所示)、裂隙长度βi的裂缝,目标裂隙通道原始的渗水能力K=βD3/12,
Figure BDA0003844311000000082
一般来说βi会小于目标裂隙通道的长度,故此可以看出修复后目标裂隙通道的渗水能力至少降低了n2倍。明显减少主裂缝通道,增加复杂的次生小裂缝,能够实现有效抑制地下水渗流目的。以中国西部矿区岩层为例,通过高压水射流处理的一条开度为10cm的主裂隙变为10个开度为0.7cm 的次生裂隙,渗透能力由K减小为343K/100000,大大缩小了,因此,对裂隙进行水力化打散之后,碎裂岩石吸水膨胀,可有效阻抑导水裂隙带纵向传导过程。
在一些实施例中,方法E采用如下方法计算修复后研究矿区覆岩的渗水能力K
E1、目标裂隙通道两壁岩土的吸水膨胀变形公式表示为Δυ=αΔx/(xl+x),Δυ表示岩体体积变化量,α表示吸水膨胀系数,xl表示吸水膨胀的含水率修正参数,x表示总的含水率,Δx表示吸水量.
E2、岩土体材料可近似看成为各向同性材料,其某一方向的线应变约为体应变的三分之一,如图2所示,具体的层状裂隙通道的开度减小系数定义为
Figure BDA0003844311000000091
E3、计算修复后研究矿区覆岩的渗水能力计算公式为:
Figure BDA0003844311000000092
考虑到岩体碎裂吸水膨胀,修复后研究矿区覆岩的渗水能力K
Figure BDA0003844311000000093
其中Di表示裂缝i的开度,βi为裂隙i的长度,λ为具体的层状裂隙通道的开度减小系数。
F、评估修复后研究矿区覆岩渗透改善情况,渗透改善情况采用包括渗流能力降低量K、渗流能力降低倍率B两个指标,K=K-K
Figure BDA0003844311000000094
在一些实施例中,步骤D2中单个目标裂隙通道的高压水射流处理后打碎成为n个开度为Di、长度为βi的裂缝,则该目标裂隙通道的渗流能力降低倍率
Figure BDA0003844311000000095
其中β为目标裂隙通道裂隙长度,D为目标裂隙通道开度,Di表示裂缝i的开度,βi为裂隙i的长度。简单举例如下:若将目标裂隙通道的大裂隙通道打碎成为n个开度为(1-Δυ)D/n、裂隙长度β(假如打碎的裂缝长度取目标裂隙通道开始的裂隙长度,实际上是远小于的)的裂缝,目标裂隙通道原始的渗水能力K=βD3/12,
Figure BDA0003844311000000101
根据
Figure BDA0003844311000000102
公式计算得到,修复后目标裂隙通道的渗水能力至少降低了n2倍。目标裂隙通道的渗流能力降低量
Figure BDA0003844311000000103
渗流能力降低量也大大降低了。本发明修复方法对裂隙进行水力化打散之后,碎裂岩石吸水膨胀,可有效阻抑导水裂隙带纵向传导过程。
一种采空区覆岩渗流能力计算与修复装置,所述装置包括:
获取模块,用于获取研究矿区覆岩的裂隙带竖向裂隙通道、弯曲下沉带的分布位置数据,获取模块内部按照覆岩岩性硬度等级构建有冒落带高度、裂隙带高度、弯曲下沉带高度的公式数据库,获取模块还用于获取研究矿区覆岩所对应的覆岩岩性硬度;
修复系统,筛选出的主要竖向裂隙通道(从研究矿区覆岩的所有裂隙进行筛选)进行竖向钻孔,在孔中布设高压水管,利用高压水管进行水压控制在目标裂隙通道中进行高压水射流处理,将目标裂隙通道的大裂隙通道打碎成为n 个开度更小的裂缝;
处理模块,用于计算研究矿区覆岩修复前后的渗水能力,并计算得到渗流能力降低量、渗流能力降低倍率评估修复后研究矿区覆岩渗透改善情况。
一种电子设备,包括:至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器执行本发明采空区覆岩渗流能力计算与修复方法的步骤。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (9)

1.一种采空区覆岩渗流能力计算与修复方法,其特征在于:其方法包括:
A、根据覆岩岩性硬度依次划分坚硬、中硬、软弱、极软四个等级,将研究矿区覆岩从下至上依次划分为冒落带、裂隙带和弯曲下沉带的三带结构,根据研究矿区的地质参数得出工作面覆岩的冒落带高度、裂隙带高度、弯曲下沉带高度;
B、利用超声波测试法对研究矿区覆岩的三带结构进行测量,得到裂隙带竖向裂隙通道、弯曲下沉带的分布位置数据,裂隙带竖向裂隙通道的分布位置数据包括竖向裂隙通道位置数据与尺寸数据,弯曲下沉带的分布位置数据包括弯曲下沉带的位置、高度数据;从裂隙带竖向裂隙通道的分布位置数据中筛选出主要的竖向裂隙通道;
C、按照如下公式计算研究矿区覆岩的渗水能力K:
K=βD3/12,其中β为裂隙长度,D为主要竖向裂隙通道开度;
D、对筛选出的主要竖向裂隙通道按照如下方法进行修复:
D1、以筛选出的主要竖向裂隙通道作为目标裂隙通道,由地面向位于目标裂隙通道上方的弯曲下沉带、目标裂隙通道竖向钻孔,在孔中布设高压水管,高压水管位于目标裂隙通道的部分开有出水喷头;
D2、利用高压水管进行水压控制在目标裂隙通道中进行高压水射流处理,将目标裂隙通道的大裂隙通道打碎成为n个开度更小的裂缝;
E、计算修复后研究矿区覆岩的渗水能力K
Figure FDA0003844310990000011
1≤i≤n,其中Di表示裂缝i的开度,βi为裂隙i的长度;
F、评估修复后研究矿区覆岩渗透改善情况,渗透改善情况采用包括渗流能力降低量K、渗流能力降低倍率B两个指标,K=K-K
Figure FDA0003844310990000012
2.按照权利要求1所述的采空区覆岩渗流能力计算与修复方法,其特征在于:方法E采用如下方法计算修复后研究矿区覆岩的渗水能力K
E1、目标裂隙通道两壁岩土的吸水膨胀变形公式表示为Δυ=αΔx/(xl+x),Δυ表示岩体体积变化量,α表示吸水膨胀系数,xl表示吸水膨胀的含水率修正参数,x表示总的含水率,Δx表示吸水量.
E2、岩土体材料可近似看成为各向同性材料,其某一方向的线应变约为体应变的三分之一,具体的层状裂隙通道的开度减小系数定义为
Figure FDA0003844310990000021
E3、计算修复后研究矿区覆岩的渗水能力计算公式为:
Figure FDA0003844310990000022
考虑到岩体碎裂吸水膨胀,修复后研究矿区覆岩的渗水能力K
Figure FDA0003844310990000023
其中Di表示裂缝i的开度,βi为裂隙i的长度。
3.按照权利要求1所述的采空区覆岩渗流能力计算与修复方法,其特征在于:
覆岩岩性为软弱的冒落带高度公式如下:
Hm=100∑M/(6.2∑M+32)±1.5,其中M表示采高,单位m;
覆岩岩性为软弱的裂隙带高度公式如下:
Hl=100∑M/(3.1∑M+5)±4.0,其中M表示采高,单位m;
覆岩岩性为软弱的弯曲下沉带高度等于采深减去裂隙带高度。
4.按照权利要求1所述的采空区覆岩渗流能力计算与修复方法,其特征在于:步骤B方法包括:
向研究矿区覆岩发出超声波,根据超声波在研究矿区覆岩的三带结构中传播速度、振幅的区别识别三带结构的具体带层及带层的高度;若超声波在研究矿区覆岩三带结构某个带层具有明显波速变化,则判断该带层为裂隙带,同时根据超声波的波速变化记录情况得出裂隙带竖向裂隙通道的分布位置数据。
5.按照权利要求1所述的采空区覆岩渗流能力计算与修复方法,其特征在于:步骤D1中钻孔包括竖直孔和L型孔,竖直孔中布设的高压水管配合为竖直形状的高压水管,L型孔中布设的高压水管配合为L型高压水管。
6.按照权利要求2所述的采空区覆岩渗流能力计算与修复方法,其特征在于:步骤D2中单个目标裂隙通道的高压水射流处理后打碎成为n个开度为Di、长度为βi的裂缝,则该目标裂隙通道的渗流能力降低倍率
Figure FDA0003844310990000031
其中β为目标裂隙通道裂隙长度,D为目标裂隙通道开度,Di表示裂缝i的开度,βi为裂隙i的长度。
7.按照权利要求6所述的采空区覆岩渗流能力计算与修复方法,其特征在于:目标裂隙通道的渗流能力降低量
Figure FDA0003844310990000032
8.一种采空区覆岩渗流能力计算与修复装置,其特征在于:所述装置包括:
获取模块,用于获取研究矿区覆岩的裂隙带竖向裂隙通道、弯曲下沉带的分布位置数据,获取模块内部按照覆岩岩性硬度等级构建有冒落带高度、裂隙带高度、弯曲下沉带高度的公式数据库,获取模块还用于获取研究矿区覆岩所对应的覆岩岩性硬度;
修复系统,筛选出的主要竖向裂隙通道进行竖向钻孔,在孔中布设高压水管,利用高压水管进行水压控制在目标裂隙通道中进行高压水射流处理,将目标裂隙通道的大裂隙通道打碎成为n个开度更小的裂缝;
处理模块,用于计算研究矿区覆岩修复前后的渗水能力,并计算得到渗流能力降低量、渗流能力降低倍率评估修复后研究矿区覆岩渗透改善情况。
9.一种电子设备,其特征在于:包括:至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器执行如权利要求1~7任一所述的方法的步骤。
CN202211113009.4A 2022-09-14 2022-09-14 采空区覆岩渗流能力计算与修复方法、装置及电子设备 Active CN115510776B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211113009.4A CN115510776B (zh) 2022-09-14 2022-09-14 采空区覆岩渗流能力计算与修复方法、装置及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211113009.4A CN115510776B (zh) 2022-09-14 2022-09-14 采空区覆岩渗流能力计算与修复方法、装置及电子设备

Publications (2)

Publication Number Publication Date
CN115510776A true CN115510776A (zh) 2022-12-23
CN115510776B CN115510776B (zh) 2023-04-28

Family

ID=84504919

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211113009.4A Active CN115510776B (zh) 2022-09-14 2022-09-14 采空区覆岩渗流能力计算与修复方法、装置及电子设备

Country Status (1)

Country Link
CN (1) CN115510776B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN206557207U (zh) * 2017-02-28 2017-10-13 贵州理工学院 一种基于数值模拟研究重复采动覆岩裂隙变化实验装置
RU2646956C1 (ru) * 2017-05-31 2018-03-12 Общество с ограниченной ответственностью "Научно-исследовательский институт природных газов и газовых технологий - Газпром ВНИИГАЗ" Способ определения трещинной пористости горных пород
CN108119142A (zh) * 2017-11-09 2018-06-05 中国矿业大学 “五图-三带-两分区”保水采煤方法
US20200370433A1 (en) * 2018-06-08 2020-11-26 China University Of Mining And Technology Risk evaluation method of overburden bed-separation water disaster in mining area
CN112377259A (zh) * 2020-11-10 2021-02-19 煤炭科学技术研究院有限公司 一种瓦斯圈闭形成的判定方法及圈闭瓦斯资源量的计算方法
CN113863987A (zh) * 2021-10-21 2021-12-31 无锡亚冠微电子有限公司 一种隧道渗水动态化智能监测系统及方法
CN114109502A (zh) * 2020-08-25 2022-03-01 神华神东煤炭集团有限责任公司 采动区域覆岩裂隙自修复程度的探测方法及其应用
CN114818236A (zh) * 2021-01-27 2022-07-29 辽宁工程技术大学 一种重建煤体孔裂隙网络结构的渗流模拟方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN206557207U (zh) * 2017-02-28 2017-10-13 贵州理工学院 一种基于数值模拟研究重复采动覆岩裂隙变化实验装置
RU2646956C1 (ru) * 2017-05-31 2018-03-12 Общество с ограниченной ответственностью "Научно-исследовательский институт природных газов и газовых технологий - Газпром ВНИИГАЗ" Способ определения трещинной пористости горных пород
CN108119142A (zh) * 2017-11-09 2018-06-05 中国矿业大学 “五图-三带-两分区”保水采煤方法
US20200370433A1 (en) * 2018-06-08 2020-11-26 China University Of Mining And Technology Risk evaluation method of overburden bed-separation water disaster in mining area
CN114109502A (zh) * 2020-08-25 2022-03-01 神华神东煤炭集团有限责任公司 采动区域覆岩裂隙自修复程度的探测方法及其应用
CN112377259A (zh) * 2020-11-10 2021-02-19 煤炭科学技术研究院有限公司 一种瓦斯圈闭形成的判定方法及圈闭瓦斯资源量的计算方法
CN114818236A (zh) * 2021-01-27 2022-07-29 辽宁工程技术大学 一种重建煤体孔裂隙网络结构的渗流模拟方法
CN113863987A (zh) * 2021-10-21 2021-12-31 无锡亚冠微电子有限公司 一种隧道渗水动态化智能监测系统及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JU J ET AL: "Self-healing effect of water-conducting fractures due to water-rock interactions in undermined rock strata and its mechanisms" *
陈善乐: "综放工作面采空区渗流场域及渗透系数研究" *

Also Published As

Publication number Publication date
CN115510776B (zh) 2023-04-28

Similar Documents

Publication Publication Date Title
Zhang et al. Enhancement of gas drainage efficiency in a special thick coal seam through hydraulic flushing
Kang et al. Mechanical behaviors of coal measures and ground control technologies for China's deep coal mines–A review
Zhang et al. Prediction of the height of the water-conducting zone above the mined panel in solid backfill mining
CN108665105B (zh) 一种导水裂缝带发育高度预计方法
Sun et al. Multi-field coupling of water inrush channel formation in a deep mine with a buried fault
Wei et al. Formation and height of the interconnected fractures zone after extraction of thick coal seams with weak overburden in Western China
CN110735663A (zh) 一种垮落带内地面注浆充填浆体形态相似模拟方法
CN113622913B (zh) 一种全垮落法开采下井上下一体化隧道围岩变形控制方法
Xue et al. The influence of magmatic rock thickness on fracture and instability law of mining surrounding rock
Parise Rock failures in karst
CN104850695A (zh) 一种确定有效采空区放顶爆破方案的方法
Wang et al. Study on blasting technology for open-pit layering of complex mine adjacent to high and steep slope
Zhang et al. Research on roof cutting pressure relief of the gob-side entry retaining with roadside backfilling
Bezruchko et al. Prognosis for free methane traps of structural and tectonic type in Donbas
Cheng et al. ASSESSMENT MODEL OF STRATA PERMEABILITY CHANGE DUE TO UNDERGROUND LONGWALL MINING.
Nan et al. Migration law of different top coal thicknesses in top coal caving
CN114352355B (zh) 深部矿井冲击地压控制方法
CN115510776A (zh) 采空区覆岩渗流能力计算与修复方法、装置及电子设备
Tian-xuan et al. Optimization of goaf gas drainage parameters based on numerical simulation studying fracture in overlying strata
Li et al. Water inrush mechanism and prevention for thick coal mining under an extremely thick glutenite layer: a case study in the southwest of the Ordos Basin
CN117027802B (zh) 地面水平井分段压裂区域超前防治煤矿冲击地压的方法
Sun et al. Research and application of hydraulic fracturing axial roof cutting technology for gently inclined hard roof based on abra-sive jet
Wan-rong Experimental and Numerical Study of Rock Stratum Movement Characteristics in Longwall Mining
Lei et al. Failure evolution of collapse column area in deep mine based on numerical calculation and microseismic monitoring
Duan et al. Pressure Relief Gas Drainage in a Fully Mechanised Mining Face Based on the Comprehensive Determination of ‘Three Zones’ Development Height

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