CN113792909A - 矿井突涌水预警分析方法 - Google Patents

矿井突涌水预警分析方法 Download PDF

Info

Publication number
CN113792909A
CN113792909A CN202110936634.8A CN202110936634A CN113792909A CN 113792909 A CN113792909 A CN 113792909A CN 202110936634 A CN202110936634 A CN 202110936634A CN 113792909 A CN113792909 A CN 113792909A
Authority
CN
China
Prior art keywords
water
aquifer
mine
characteristic information
water level
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.)
Pending
Application number
CN202110936634.8A
Other languages
English (en)
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.)
Tiandi Science and Technology Co Ltd
CCTEG Coal Mining Research Institute
Original Assignee
Tiandi Science and Technology Co Ltd
CCTEG Coal Mining Research Institute
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 Tiandi Science and Technology Co Ltd, CCTEG Coal Mining Research Institute filed Critical Tiandi Science and Technology Co Ltd
Priority to CN202110936634.8A priority Critical patent/CN113792909A/zh
Publication of CN113792909A publication Critical patent/CN113792909A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
    • E21F17/18—Special adaptations of signalling or alarm devices
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Agronomy & Crop Science (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (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)

Abstract

本发明提供了一种矿井突涌水预警分析方法。包括获取矿井的多个层结构的特征信息，以根据多个特征信息筛选出重点含水层和监测含水层，当含水层的水位下降至预设水位值时，对含水层的涌水量进行分析的步骤；根据多个特征信息判定存在充水危险性的层结构的步骤；根据多个特征信息的改变动态调整充水危险性的步骤。应用本实施例中的矿井突涌水预警分析方法时，可以获取矿井的多个层结构的特征信息，并进行比对分析来筛选具有充水危险性的层结构，并且可以根据获取的特征信息进行动态调整，相较于传统的单一参数判定方法，本方不仅可以提高预警分析的准确度，还可以进行实时监测，从而保证作业的安全性，使用效果好。

Description

矿井突涌水预警分析方法
技术领域
本发明涉及煤炭开采技术领域，尤其涉及一种矿井突涌水预警分析方法。
背景技术
目前，煤矿顶板含水层涌水危险性监测预警是当前煤矿水害防控的一项重要措施，但在目前煤矿监测预警中，主要是以单一阈值的改变作为依据进行预警分析，例如当含水层水位发生一定程度的改变即进行报警，而受丰水期、枯水期等因素影响，容易造成含水层中出现水位波动，这就导致上述监测方法容易出现水害误报的现象，监测准确度不佳。
此外，在对工作面进行突涌水评价时，往往以地面抽水试验所获取的参数为准，但实际生产过程中，相关数据往往与实际情况存在较大的差异，这就导致工作面充水危险性预警和工作面涌水量的预测均与实际作业需求存在较大的偏差，作业的安全性无法得到有效保障。
因此，需要针对上述问题进行改进，以改变现状。
发明内容
本发明提供一种矿井突涌水预警分析方法，用于解决传统矿井中监测突涌水存在监测精度较差，从而容易导致作业容易出现安全风险的问题。
本发明提供了一种矿井突涌水预警分析方法，包括：
获取矿井的多个层结构特征信息，以根据多个所述特征信息筛选出重点含水层和监测含水层，当含水层的水位下降至预设水位值时，对含水层的涌水量进行分析；
根据多个所述特征信息判定存在充水危险性的所述层结构；
根据多个所述特征信息的改变动态调整所述充水危险性。
根据本发明的一个实施例，所述获取矿井的多个层结构的特征信息，以根据多个所述特征信息筛选出重点含水层和监测含水层的步骤之后，所述方法还包括：
当含水层的富水性达到预设值之后定义为监测含水层，并对所述监测含水层设置监测钻孔，且所述监测钻孔设置在采动影响范围以外。
根据本发明的一个实施例，所述获取矿井的多个层结构的特征信息的步骤中，还包括获取顶板岩层的底板标高Hi；所述当含水层的水位下降至预设水位值时，对工作面的涌水量进行分析的步骤，还包括：
通过所述监测钻孔获取含水层的水位变化数值，当含水层的水位下降至预设水位值时，可对工作面的涌水量进行分析操作；
获取煤层开采深度H，计算得出采空区重新压实范围为0.3H。
根据本发明的一个实施例，所述获取矿井的多个层结构的特征信息，以根据多个所述特征信息筛选出重点含水层和监测含水层的步骤还包括，获取顶板含隔水层厚度Mi、煤层与顶板岩层的底板标高Hi；
所述当含水层的水位下降至预设水位值时，对工作面的预计涌水量进行分析的步骤包括：
采用稳定流计算公式计算工作面预计涌水量，所述稳定流计算公式为：
Figure BDA0003213463520000021
其中，Q预为预计涌水量，K为渗透系数，H0为自然水位至含水层底板的距离，h0′为水文监测钻孔获取的含水层水位标高与含水层底板标高的差值，R0为引用影响半径，r0为引用半径
根据本发明的一个实施例，所述引用半径的计算公式为：
Figure BDA0003213463520000031
其中，其中H为煤层埋深，L为工作面斜长。
根据本发明的一个实施例，所述引用影响半径的计算公式为：
Figure BDA0003213463520000032
其中，L0为水位降深。
根据本发明的一个实施例，所述通过所述监测钻孔获取含水层的水位变化数值，当含水层的水位下降至预设水位值时，可对工作面的涌水量进行分析操作的步骤之后，还包括：
获取工作面的涌水量数值，并与预计涌水量比对；当工作面仅第一含水层受采动扰动影响时，渗透系数的计算公式为：
Figure BDA0003213463520000033
Figure BDA0003213463520000034
其中，K01′为修正后的渗透系数，Q测均1为所获取工作面的涌水量，K01为初始渗透系数，q01′为修正后的工作面涌水量数值，q01为初始工作面涌水量数值。
根据本发明的一个实施例，所述通过所述监测钻孔获取含水层的水位变化数值，当含水层的水位下降至预设水位值时，可对工作面的涌水量进行分析操作的步骤之后，还包括：
获取工作面的涌水量数值，并与预计涌水量比对；当工作面同时受第一、第二含水层受采动扰动影响时，渗透系数的计算公式为：
Figure BDA0003213463520000041
Figure BDA0003213463520000042
其中，K02′为修正后的渗透系数，Q测均为所获取工作面的涌水量，K02为初始渗透系数，q02′为修正后的工作面涌水量数值，q02为初始工作面涌水量数值。
根据本发明的一个实施例，所述获取矿井工作面的多个层结构的特征信息，以根据筛选出重点含水层和监测含水层的步骤，还包括：
获取覆岩破坏高度H_H并计算裂采比，所述裂采比的计算公式为：
Figure BDA0003213463520000043
其中，Hlie为裂采比数值，m为煤层开采厚度。
根据本发明的一个实施例，所述根据多个所述特征信息判定存在充水危险性的所述层结构的步骤，还包括：
获取导水裂缝带的预计高度L1和隔水层的厚度L2并进行判断，若L1>L2，则判断充水危险性大于预设危险性；若L1<L2则判断充水危险性小于预设危险性。
根据本发明的一个实施例，所述根据多个所述特征信息筛选出重点含水层和监测含水层的步骤，还包括：
多个所述特征信息包括含水层厚度、单位涌水量、渗透系数。
根据本发明的一个实施例，所述获取矿井的多个层结构的特征信息的步骤，还包括：
获取矿井的覆岩结构中含水层和隔水层的特征信息。
根据本发明的一个实施例，所述根据多个所述特征信息判定存在充水危险性的所述层结构的步骤，还包括：
多个所述特征信息包括含水层富水性参数、含水层与煤层之间的间距。
根据本发明的一个实施例，所述根据多个所述特征信息的改变动态调整所述充水危险性的步骤，还包括：
多个所述特征信息包括距离煤层顶板最近的含水层水位变化、煤层顶板最近含水层的相邻含水层水位变化、各含水层之间的位置关系。
实施本发明实施例，具有如下有益效果：
应用本实施例中的矿井突涌水预警分析方法时，可以获取矿井的例如覆岩结构中含隔水层等特征信息，并进行比对分析来筛选具有充水危险性的层结构，并且可以根据获取的例如含水层水位变化的特征信息进行动态调整，同时针对性监测相邻含水层的水位变化，并调整相邻含水层的监测频率等，相较于传统的单一参数判定方法，本方法不仅可以提高预警分析的准确度，还可以进行实时监测，从而保证作业的安全性，使用效果好。
附图说明
为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
其中：
图1是本发明的实施例中矿井突涌水预警分析方法的流程示意图；
图2是本发明的实施例中矿井突涌水预警分析方法的层结构分布示意图；
图3是本发明的实施例中矿井突涌水预警分析方法的监测钻孔的布置示意图；
图4是本发明的实施例中工作面采动影响范围的结构示意图；
附图标记：
10、煤层；11、第一隔水层；12、第一含水层；13、第二隔水层；14、第二含水层；
100、导水裂缝；200、监测钻孔；300、岩层移动边界角；400、工作面；500、采动影响范围；600、采空区。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在本发明实施例的描述中，需要说明的是，术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明实施例和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明实施例的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
在本发明实施例的描述中，需要说明的是，除非另有明确的规定和限定，术语“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明实施例中的具体含义。
在本发明实施例中，除非另有明确的规定和限定，第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触，或第一和第二特征通过中间媒介间接接触。而且，第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方，或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方，或仅仅表示第一特征水平高度小于第二特征。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
参阅图1所示，本发明实施例提供了一种矿井突涌水预警分析方法，其包括：
步骤S100、获取矿井的覆岩结构中含隔水层的特征信息，以根据含水层厚度、单位涌水量、渗透系数等特征信息筛选出重点含水层和监测含水层，当含水层的水位下降至预设水位值时，对工作面的涌水量进行分析；
步骤S200、根据多个特征信息判定存在充水危险性的层结构；
步骤S300、根据多个特征信息的改变动态调整充水危险性。
应用本实施例中的矿井突涌水预警分析方法时，可以获取矿井的覆岩结构中含隔水层等特征信息，并进行比对分析来筛选具有充水危险性的含水层，并且可以根据获取的含水层水位变化特征信息进行动态调整，同时针对性监测相邻含水层的水位变化，并调整相邻含水层的监测频率等，相较于传统的单一参数判定方法，本方不仅可以提高预警分析的准确度，还可以进行实时监测，从而保证作业的安全性，使用效果好。
进一步地，步骤S100之后，参阅图3和图4所示，本矿井突涌水预警分析方法还包括：步骤S400、当含水层的富水性满足预设条件时定义为监测含水层，并对监测含水层设置监测钻孔200，且监测钻孔200设置在采动影响范围500以外。
具体地，步骤S100还包括：获取覆岩破坏高度H_H并计算裂采比，裂采比的计算公式为：
Figure BDA0003213463520000081
其中，Hlie为裂采比数值，m为煤层10开采厚度。
应用本实施例的矿井突涌水预警分析方法的第一阶段时，可构建工作面三维参数模型：
第一步、统计工作面400的煤层10厚度、含水层和隔水层的厚度：具体为，统计研究区域钻孔柱状图，提取各钻孔煤层10厚度m0、顶板含隔水层厚度Mi，煤层10与顶板岩层的底板标高Hi，并构建工作面400三维模型；
第二步、开展工作面400含水层富水性分析：具体为，开展工作面400顶板含水层抽水试验，获取含水层渗透系数K、单位涌水量q、含水层水位标高h等参数；
第三步、预计工作面400覆岩破坏范围：具体为，统计工作面400及周边钻孔柱状图的岩性，分析上覆岩层力学参数，结合《建筑物、水体下、铁路及主要井巷煤柱留设与压煤开采规范》给出的相关参数，以确定覆岩的类型，以工作面开采厚度为基数，预计覆岩破坏高度H_H，分析工作面开采过程中可能波及含水层。
参阅图3所示，应用本实施例的矿井突涌水预警分析方法的第二阶段时，可建立工作面上覆含水层监测：
第一步、确定工作面400需要监测的含水层层位：具体为，结合覆岩破坏高度预计范围和含水层富水特性，以富水性、导水裂缝100带波及范围为依据，筛选出监测的重点含水层，并针对富水性中等及以上、可能采动影响范围500波及到的含水层来布置监测钻孔200；
第二步、确定水文长观孔布置的位置：具体为，根据工作面400开采影响边界以外，如图3所示，角300为岩层移动边界角，将监测钻孔200布置在采动影响范围500以外，可将测钻孔200布置在开切眼附近、工作面400见方附近、以及工作面停采线附近。
具体地，步骤S100中，还包括获取顶板岩层的底板标高H；当含水层的水位下降至预设水位值时，对工作面的涌水量进行分析的步骤，还包括：通过监测钻孔200获取含水层的水位变化数值，当含水层的水位下降至预设水位值时，可对含水层的涌水量进行分析操作；获取工作面400的煤层埋藏深度，计算采空区600重新压实范围为0.3H。
进一步地，步骤S100还包括获取顶板含隔水层厚度Mi、煤层10与顶板岩层的底板标高Hi；当含水层的水位下降至预设水位值时，对含水层的预计涌水量进行分析的步骤包括：采用稳定流计算公式计算矿井预计涌水量，稳定流计算公式为：
Figure BDA0003213463520000091
其中，Q预为预计涌水量，K为渗透系数，H0为自然水位至含水层底板的距离，h0′为水文监测钻孔获取的含水层水位标高与含水层底板标高的差值，R0为引用影响半径，r0为引用半径
具体地，在上述实施例中，引用半径r0的计算公式为：
Figure BDA0003213463520000101
其中，其中H为工作面400煤层的埋深，L为工作面400斜长。
具体地，引用影响半径R0的计算公式为：
Figure BDA0003213463520000102
其中，L0为水位降深。
具体地，步骤S100还包括：获取导水裂缝100带的预计高度L1和隔水层的厚度L2并进行判断，若L1>L2，则判断充水危险性大于预设危险性；若L1<L2则判断充水危险性小于预设危险性。应用本实施例的矿井突涌水预警分析方法的第三阶段时，可开展工作面的充水危险性分析：
第一步、计算开采对上覆含水层的影响：具体为，按照导水裂缝带100的预计高度、煤层顶板与含水层之间的隔水层厚度，计算两者的差值，当导水裂缝带高度大于煤层顶板与含水层之间的隔水层厚度，则具有充水危险性，反之充水危险性较小，计算并绘制出工作面的开采影响程度等值线图；具体地，在本实施例中，可以将各个钻孔的两者差值输入到ArcGIS软件中，绘制出工作面的开采影响程度等值线图；
第二步、分析含水层富水特性：具体为，根据含水层单位涌水量q，绘制含水层富水性等值线图；具体地，可以采用ArcGIS软件，绘制含水层富水性等值线图；
第三步、分析工作面充水危险性：具体为，将开采影响程度等值线图和富水性等值线图进行叠加，获取工作面充水危险性图；具体地，可以采用ArcGIS软件中的叠加功能，以将开采影响程度等值线图和富水性等值线图进行叠加，获取工作面充水危险性图。
应用本实施例的矿井突涌水预警分析方法的第四阶段时，可开展工作面400的充水危险性分析：
第一步、监测含水层水位变化：具体为，通过地面水位监测钻孔200监测含水层水位变化，当含水层水位下降为预设水位值时，可以开展工作面涌水量分析工作；
第二步、确定工作面采动影响范围：具体为，统计首采工作面煤层埋深，确定首采工作面回采过程中采空区重新压实范围。根据矿山压力岩层移动经验，如图4所示，确定埋深的0.3倍为采空区重新压实范围；
第三步、预计工作面涌水量：具体为，可按照“大井法”，结合周边钻孔获取的含水层水位监测数据，针对出现水位变化的含水层，分析前后水位变化情况，采用承压转无压完整井的稳定流计算公式，预计工作面回采过程中矿井涌水量，具体地，稳定流计算公式为：
Figure BDA0003213463520000111
第四步、对比分析工作面涌水量、修正渗透系数和单位涌水量：具体为，通过地面水文长观孔进行监测含水层的变化，结合水位变化情况，分别制定相应的调整方案。
进一步地，通过监测钻孔200获取含水层的水位变化数值，当含水层的水位下降至预设水位值时，可对工作面的涌水量进行分析操作的步骤之后，还包括步骤：获取工作面的涌水量数值，并与预计涌水量比对；
当工作面仅第一含水层受采动扰动影响时，通过工作面监测获取的涌水量Q测均1与Q预之间的关系，对比分析两者之间的差值。观测工作面获取的最近10组观测数据的平均值Q测均1，结合两者对比，以R01(第一含水层的引用影响半径)为影响半径，调整影响半径范围以内的水文孔参数，以K01为初始渗透系数为基数进行调整渗透系数，渗透系数的计算公式为：
Figure BDA0003213463520000121
K01＝K01′
Figure BDA0003213463520000122
其中，K01′为修正后的渗透系数，Q测均1为所获取工作面的涌水量，K01为初始渗透系数，q01′为修正后的工作面涌水量数值，q01为初始工作面涌水量数值。
具体地，当工作面同时受第一含水层、第二含水层受采动扰动影响时，工作面导水裂缝带波及到第二含水层时，修正工作面裂采比，将原有裂采比Hlie修正为Hlie′，且Hlie′的计算公式为：
Figure BDA0003213463520000123
需要说明的是，在本实施例中，MM为煤层顶板距第一层含水层之间的距离，Mi为上覆含隔水层厚度。M2x-1(x为正整数，且x>1)为第2x-1隔水层，M2x(x为正整数，且x>1)为第x含水层。
通过工作面监测获取的涌水量Q测均与Q预(Q预1+Q预2)之间的关系，对比分析两者之间的差值，观测工作面获取的最近10组观测数据的平均值Q测均，结合两者对比，以K01和K02为初始渗透系数为基数进行调整渗透系数。考虑到K01在前期回采过程中已结合实测数据进行了修改、调整，本次调整主要对K02进行调整，因此主要对比Q测均-Q预1和Q预2之间的关系，以R02(第二含水层的引用影响半径)为影响半径，调整影响半径范围以内的水文孔参数，并修正K02，渗透系数的计算公式为：
Figure BDA0003213463520000131
K02＝K02′
Figure BDA0003213463520000132
其中，K02′为修正后的渗透系数，Q测均为所获取工作面的涌水量，K02为初始渗透系数，q02′为修正后的工作面涌水量数值，q02为初始工作面涌水量数值。
应用本实施例的矿井突涌水预警分析方法的第五阶段时，可实现动态调整工作面突涌水危险性：
具体为，结合修正后的工作面涌水量q01′、q02′和工作面导水裂缝带高度
Figure BDA0003213463520000133
运用ArcGIS软件，修正工作面上覆含水层富水性等值线图和工作面开采影响程度等值线图，运用ArcGIS的叠加功能，将开采影响等值线图和富水性等值线图进行叠加，获取修正后的工作面突涌水危险性等值线图。
具体地，参阅图2所示，在本实施例中：
HH＝MM+M1+M2+M3+M4
在本实施例的矿井突涌水预警分析方法的步骤中，构建了工作面涌水量动-静态监测预警算法，并且实现对工作面涌水量进行动态监测、对相关参数进行动态调整的功能，另外，通过对工作面含水层单位涌水量参数、导水裂缝带高度的调整，实现工作面充水危险性评价的动态调整，使用效果好。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (14)

1.一种矿井突涌水预警分析方法，其特征在于，包括：
获取矿井的多个层结构特征信息，以根据多个所述特征信息筛选出重点含水层和监测含水层，当含水层的水位下降至预设水位值时，对含水层的涌水量进行分析；
根据多个所述特征信息判定存在充水危险性的所述层结构；
根据多个所述特征信息的改变动态调整所述充水危险性。
2.根据权利要求1所述的矿井突涌水预警分析方法，其特征在于，所述获取矿井的多个层结构的特征信息，以根据多个所述特征信息筛选出重点含水层和监测含水层的步骤之后，所述方法还包括：
当含水层的富水性达到预设值之后定义为监测含水层，并对所述监测含水层设置监测钻孔，且所述监测钻孔设置在采动影响范围以外。
3.根据权利要求2所述的矿井突涌水预警分析方法，其特征在于，所述获取矿井的多个层结构的特征信息的步骤中，还包括获取顶板岩层的底板标高Hi；所述当含水层的水位下降至预设水位值时，对工作面的涌水量进行分析的步骤，还包括：
通过所述监测钻孔获取含水层的水位变化数值，当含水层的水位下降至预设水位值时，可对工作面的涌水量进行分析操作；
获取煤层开采深度H，计算得出采空区重新压实范围为0.3H。
4.根据权利要求3所述的矿井突涌水预警分析方法，其特征在于，所述获取矿井的多个层结构的特征信息，以根据多个所述特征信息筛选出重点含水层和监测含水层的步骤还包括，获取顶板含隔水层厚度Mi、煤层与顶板岩层的底板标高Hi；
所述当含水层的水位下降至预设水位值时，对工作面的预计涌水量进行分析的步骤包括：
采用稳定流计算公式计算工作面预计涌水量，所述稳定流计算公式为：
Figure FDA0003213463510000021
其中，Q预为预计涌水量，K为渗透系数，H0为自然水位至含水层底板的距离，h0′为水文监测钻孔获取的含水层水位标高与含水层底板标高的差值，R0为引用影响半径，r0为引用半径。
5.根据权利要求4所述的矿井突涌水预警分析方法，其特征在于，所述引用半径的计算公式为：
Figure FDA0003213463510000022
其中，其中H为煤层埋深，L为工作面斜长。
6.根据权利要求5所述的矿井突涌水预警分析方法，其特征在于，所述引用影响半径的计算公式为：
Figure FDA0003213463510000023
其中，L0为水位降深。
7.根据权利要求4-6任意一项所述的矿井突涌水预警分析方法，其特征在于，所述通过所述监测钻孔获取含水层的水位变化数值，当含水层的水位下降至预设水位值时，可对工作面的涌水量进行分析操作的步骤之后，还包括：
获取工作面的涌水量数值，并与预计涌水量比对；当工作面仅第一含水层受采动扰动影响时，渗透系数的计算公式为：
Figure FDA0003213463510000024
Figure FDA0003213463510000025
其中，K01′为修正后的渗透系数，Q测均1为所获取工作面的涌水量，K01为初始渗透系数，q01′为修正后的工作面涌水量数值，q01为初始工作面涌水量数值。
8.根据权利要求4-6任意一项所述的矿井突涌水预警分析方法，其特征在于，所述通过所述监测钻孔获取含水层的水位变化数值，当含水层的水位下降至预设水位值时，可对工作面的涌水量进行分析操作的步骤之后，还包括：
获取工作面的涌水量数值，并与预计涌水量比对；当工作面同时受第一、第二含水层受采动扰动影响时，渗透系数的计算公式为：
Figure FDA0003213463510000031
Figure FDA0003213463510000032
其中，K02′为修正后的渗透系数，Q测均为所获取工作面的涌水量，K02为初始渗透系数，q02′为修正后的工作面涌水量数值，q02为初始工作面涌水量数值。
9.根据权利要求1所述的矿井突涌水预警分析方法，其特征在于，所述获取矿井工作面的多个层结构的特征信息，以根据筛选出重点含水层和监测含水层的步骤，还包括：
获取覆岩破坏高度H_H并计算裂采比，所述裂采比的计算公式为：
Figure FDA0003213463510000033
其中，Hlie为裂采比数值，m为煤层开采厚度。
10.根据权利要求1所述的矿井突涌水预警分析方法，其特征在于，所述根据多个所述特征信息判定存在充水危险性的所述层结构的步骤，还包括：
获取导水裂缝带的预计高度L1和隔水层的厚度L2并进行判断，若L1>L2，则判断充水危险性大于预设危险性；若L1<L2则判断充水危险性小于预设危险性。
11.根据权利要求1所述的矿井突涌水预警分析方法，其特征在于，所述根据多个所述特征信息筛选出重点含水层和监测含水层的步骤，还包括：
多个所述特征信息包括含水层厚度、单位涌水量、渗透系数。
12.根据权利要求1所述的矿井突涌水预警分析方法，其特征在于，所述获取矿井的多个层结构的特征信息的步骤，还包括：
获取矿井的覆岩结构中含水层和隔水层的特征信息。
13.根据权利要求1所述的矿井突涌水预警分析方法，其特征在于，所述根据多个所述特征信息判定存在充水危险性的所述层结构的步骤，还包括：
多个所述特征信息包括含水层富水性参数、含水层与煤层之间的间距。
14.根据权利要求1所述的矿井突涌水预警分析方法，其特征在于，所述根据多个所述特征信息的改变动态调整所述充水危险性的步骤，还包括：
多个所述特征信息包括距离煤层顶板最近的含水层水位变化、煤层顶板最近含水层的相邻含水层水位变化、各含水层之间的位置关系。
CN202110936634.8A 2021-08-16 2021-08-16 矿井突涌水预警分析方法 Pending CN113792909A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110936634.8A CN113792909A (zh) 2021-08-16 2021-08-16 矿井突涌水预警分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110936634.8A CN113792909A (zh) 2021-08-16 2021-08-16 矿井突涌水预警分析方法

Publications (1)

Publication Number Publication Date
CN113792909A true CN113792909A (zh) 2021-12-14

Family

ID=79181676

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110936634.8A Pending CN113792909A (zh) 2021-08-16 2021-08-16 矿井突涌水预警分析方法

Country Status (1)

Country Link
CN (1) CN113792909A (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109377042A (zh) * 2018-10-17 2019-02-22 中国煤炭地质总局水文地质局 一种煤层顶板突水危险性评价方法
CN109767136A (zh) * 2019-01-17 2019-05-17 内蒙古上海庙矿业有限责任公司 煤层顶板含水层涌突水危险性评价与预测四双工作法
CN111191849A (zh) * 2020-01-03 2020-05-22 中煤科工集团西安研究院有限公司 一种西部矿区深埋工作面涌水量预测方法
CN111859256A (zh) * 2020-07-17 2020-10-30 中煤科工集团西安研究院有限公司 采煤工作面顶板水害时空控制疏水防治方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109377042A (zh) * 2018-10-17 2019-02-22 中国煤炭地质总局水文地质局 一种煤层顶板突水危险性评价方法
CN109767136A (zh) * 2019-01-17 2019-05-17 内蒙古上海庙矿业有限责任公司 煤层顶板含水层涌突水危险性评价与预测四双工作法
CN111191849A (zh) * 2020-01-03 2020-05-22 中煤科工集团西安研究院有限公司 一种西部矿区深埋工作面涌水量预测方法
CN111859256A (zh) * 2020-07-17 2020-10-30 中煤科工集团西安研究院有限公司 采煤工作面顶板水害时空控制疏水防治方法

Similar Documents

Publication Publication Date Title
CN102705010B (zh) 一种工作面突出危险评价指标体系分析方法
CN111691872B (zh) 一种煤层覆岩裂隙发育范围及程度的动态精准探测方法
CN109767136A (zh) 煤层顶板含水层涌突水危险性评价与预测四双工作法
CN104847379B (zh) 一种奥灰顶部注浆加固工作面底板的方法
CN106337426A (zh) 一种增加承压水地层基坑抗突涌稳定性的降水方法
CN111612643B (zh) 一种瓦斯抽采对象与抽采措施优选匹配方法
CN102691522A (zh) 一种工作面突出危险动态预测图形成方法
CN111932128B (zh) 一种煤层底板承压含水层定向钻注浆层位确定方法
CN108536891B (zh) 地下水库储水量的计算方法、存储介质和装置
CN110359841A (zh) 水平井的轨迹确定方法、水平井的钻井控制方法和装置
CN101936154B (zh) 一种煤层气水平井钻井架桥导向的方法
CN113792909A (zh) 矿井突涌水预警分析方法
CN116070416A (zh) 一种确定采场覆岩竖向三带发育高度的综合测试方法
CN109635357A (zh) 一种考虑采动岩体碎胀性的覆岩离层动态位置预测方法
CN112502775B (zh) 一种半固结砂岩含水层疏水降压方法
CN113324506A (zh) 双模态盾构隧道下穿水库引起导水裂隙带高度的预测方法
CN108952569A (zh) 一种底水层状油气藏水平井钻探轨迹控制方法
CN116792071B (zh) 一种适用于弱胶结储层的流砂塑岩固砂防砂方法
CN116611567B (zh) 一种矿区采动覆岩顶板突水复合灾害危险性分区预测方法
CN114718558B (zh) 一种侵入岩潜山随钻快速判别储层方法
CN113833432B (zh) 水力压裂弱化卸压施工方法及施工系统
CN114184153B (zh) 基于光纤及渗压计的采场覆岩及土层复合导高监测方法
CN114294059B (zh) 一种坚硬岩层水力压裂层位高精度控制方法
CN113982574B (zh) 地面水害区域治理过程中的地质补勘方法
CN114021489A (zh) 一种基于抗覆岩渗透破坏评价的防砂煤柱留设方法

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