CN107886185A - 一种基于空间分析的矿井静涌水量计算方法及系统 - Google Patents
一种基于空间分析的矿井静涌水量计算方法及系统 Download PDFInfo
- Publication number
- CN107886185A CN107886185A CN201710948510.5A CN201710948510A CN107886185A CN 107886185 A CN107886185 A CN 107886185A CN 201710948510 A CN201710948510 A CN 201710948510A CN 107886185 A CN107886185 A CN 107886185A
- Authority
- CN
- China
- Prior art keywords
- water
- bearing strata
- quiet
- prediction
- model
- 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
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 442
- 238000012732 spatial analysis Methods 0.000 title claims abstract description 30
- 238000000205 computational method Methods 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 6
- 239000003673 groundwater Substances 0.000 claims description 35
- 238000009826 distribution Methods 0.000 claims description 25
- 238000004458 analytical method Methods 0.000 claims description 7
- 239000004744 fabric Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 abstract description 17
- 238000005065 mining Methods 0.000 abstract description 16
- 238000004321 preservation Methods 0.000 abstract description 8
- 239000003245 coal Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical group O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 description 5
- 238000010276 construction Methods 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000004575 stone Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000005641 tunneling Effects 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 235000013312 flour Nutrition 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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"
-
- 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
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Animal Husbandry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Mining & Mineral Resources (AREA)
- General Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于空间分析的矿井静涌水量计算方法及系统,该方法包括以下步骤:根据预先录入的充水含水层的特征参数,建立充水含水层的分布模型;根据预先设定的预测域,建立预测域模型;根据所述充水含水层的分布模型和所述预测域模型,建立预测域静涌水量模型;根据所述预测域静涌水量模型,计算矿井静涌水量。本发明通过对工作面开采影响范围的含水层所赋存地下水的静储量的计算评价,提高了对矿井涌水量的预测结果的准确性。
Description
技术领域
本发明涉及矿井水文地质领域,具体涉及一种基于空间分析的矿井静涌水量计算方法及系统。
背景技术
矿井涌水量是煤矿建井、防治水中长期规划、煤矿水文地质勘查的主要要素。准确的涌水量预测对于合理确定防治水措施,优选施工方案是极其重要的。井筒涌水量预测的准确与否,直接关系到建井工期、施工安全和经济效益。
我国西部的煤系地层上覆着的含水层的富水性较差,但由于其充水含水层厚度较大,因而也赋存较大量地下水,在工作面回采过程中,开采影响范围内的含水层所赋存的地下水也会较短时间内溃入矿井,从而威胁着煤矿安全生产。
然而,目前涌水量的计算方法主要是针对动涌水量,即导水裂隙带所影响的含水层的地下水单位时间流经预测范围横断面的体积,而对于工作面开采影响范围的含水层所赋存地下水的静储量没有进行计算评价,从而导致对矿井涌水量的预测结果不准确。
发明内容
本发明针对现有的涌水量计算方法存在的上述技术问题,提供一种基于空间分析的矿井静涌水量计算方法及系统。
本发明解决上述技术问题的技术方案如下:一种基于空间分析的矿井静涌水量计算方法,该方法包括以下步骤:
步骤1、根据预先录入的充水含水层的特征参数,建立充水含水层的分布模型;
步骤2、根据预先设定的预测域,建立预测域模型;
步骤3、根据所述充水含水层的分布模型和所述预测域模型,建立预测域静涌水量模型;
步骤4、根据所述预测域静涌水量模型,计算矿井静涌水量。
本发明的有益效果是:通过对工作面开采影响范围的含水层所赋存地下水的静储量的计算评价,提高了对矿井涌水量的预测结果的准确性。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,所述特征参数包括:给水度、地下水水位标高和充水含水层底板标高,所述充水含水层的分布模型包括给水度分布模型、预测前地下水水位标高分布模型和充水含水层底板标高分布模型;
所述根据充水含水层的特征参数,建立充水含水层的分布模型具体包括:
根据所述给水度建立所述给水度分布模型;
根据所述地下水水位标高建立所述预测前地下水水位标高分布模型;
根据所述充水含水层底板标高建立所述充水含水层底板标高分布模型。
进一步,所述充水含水层的分布模型还包括充水含水层厚度分布模型,所述根据所述充水含水层的分布模型和所述预测域模型,建立预测域静涌水量模型具体包括:
根据所述地下水水位标高和所述充水含水层底板标高,建立充水含水层的充水含水层厚度分布模型;根据所述充水含水层的充水含水层厚度分布模型、给水度分布模型和预测域模型,采用空间分析方法建立预测域静涌水量模型。
进一步,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层中任一层的单位面积静涌水总量,所述根据所述预测域静涌水量模型,计算矿井静涌水量具体包括:
步骤4.1、从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度和充水含水层厚度;
步骤4.2、分别根据所述至少一个预测域单元中每一个预测域单元的给水度和充水含水层厚度,计算每一个预测域单元的单位面积静涌水量;
步骤4.3、根据每一个预测域单元的单位面积静涌水量,统计所述预测域中第一充水含水层的单位面积静涌水总量,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层。
进一步,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层的静涌水总量,所述根据所述预测域静涌水量模型,计算矿井静涌水量具体包括:
步骤4.1、从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度、充水含水层厚度和面积;
步骤4.2、分别根据所述至少一个预测域单元中每一个预测域单元的给水度、充水含水层厚度和面积,计算每一个预测域单元的静涌水量;
步骤4.3、根据每一个预测域单元的静涌水量,统计所述预测域中第一充水含水层的静涌水总量;
步骤4.4、统计所述至少一层充水含水层中的每一层充水含水层的静涌水总量之和,作为所述矿井静涌水量,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层。
为实现上述发明目的,本发明还提供一种基于空间分析的矿井静涌水量计算系统,该系统包括:
第一建立模块,用于预先录入的充水含水层的特征参数,建立充水含水层的分布模型;
第二建立模块,用于根据预先设定的预测域,建立预测域模型;
第三建立模块,用于根据所述充水含水层的分布模型和所述预测域模型,建立预测域静涌水量模型;
计算模块,用于根据所述预测域静涌水量模型,计算矿井静涌水量。
本发明的有益效果是:通过对工作面开采影响范围的含水层所赋存地下水的静储量的计算评价,提高了对矿井涌水量的预测结果的准确性。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,所述特征参数包括:给水度、地下水水位标高和充水含水层底板标高,所述充水含水层的分布模型包括给水度分布模型、预测前地下水水位标高分布模型和充水含水层底板标高分布模型;
所述第一建立模块具体包括:
第一建立单元,用于根据所述给水度建立所述给水度分布模型;
第二建立单元,用于根据所述地下水水位标高建立所述预测前地下水水位标高分布模型;
第三建立单元,用于根据所述充水含水层底板标高建立所述充水含水层底板标高分布模型。
进一步,所述充水含水层的分布模型还包括充水含水层厚度分布模型,所述第三建立模块具体包括:
第四建立单元,用于根据所述地下水水位标高和所述充水含水层底板标高,建立充水含水层的充水含水层厚度分布模型;
第五建立单元,用于根据所述充水含水层的充水含水层厚度分布模型、给水度分布模型和预测域模型,采用空间分析方法建立预测域静涌水量模型。
进一步,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层中任一层的单位面积静涌水总量,所述计算模块具体包括:
第一获取单元,用于从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度和充水含水层厚度;
第一计算单元,用于分别根据所述至少一个预测域单元中每一个预测域单元的给水度和充水含水层厚度,计算每一个预测域单元的单位面积静涌水量;
第二计算单元,用于根据每一个预测域单元的单位面积静涌水量,统计所述预测域中第一充水含水层的单位面积静涌水总量,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层。
进一步,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层的静涌水总量,所述计算模块具体包括:
第二获取单元,用于从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度、充水含水层厚度和面积;
第三计算单元,用于分别根据所述至少一个预测域单元中每一个预测域单元的给水度、充水含水层厚度和面积,计算每一个预测域单元的静涌水量;
第一统计单元,用于根据每一个预测域单元的静涌水量,统计所述预测域中第一充水含水层的静涌水总量;
第二统计单元,用于统计所述至少一层充水含水层中的每一层充水含水层的静涌水总量之和,作为所述矿井静涌水量,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层。
附图说明
图1为本发明实施例提供的一种基于空间分析的矿井静涌水量计算方法的流程图;
图2为本发明实施例提供的一种基于空间分析的矿井静涌水量计算系统的结构框图;
图3为某矿含水层概化剖面示意图;
图4为预测前地下水水位标高分布模型的等值线计算图;
图5为厚度分布模型的等值线计算图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
图1为本发明实施例提供的一种基于空间分析的矿井静涌水量计算方法的流程图,如图1所示,该方法包括以下步骤:
S1、根据预先录入的充水含水层的特征参数,建立充水含水层的分布模型;
具体的,各个充水含水层的相关特征参数是通过对开采区域进行钻孔揭露及抽水试验等取得的,并且通过在如Mapgis,Arcgis等空间分析软件中输入特征参数来完成模型的建立。
可选地,在该实施例中,所述特征参数包括:给水度、地下水水位标高和充水含水层底板标高,所述充水含水层的分布模型包括给水度分布模型、预测前地下水水位标高分布模型和充水含水层底板标高分布模型;
所述根据充水含水层的特征参数,建立充水含水层的分布模型具体包括:
根据所述给水度建立所述给水度分布模型;
根据所述地下水水位标高建立所述预测前地下水水位标高分布模型;
根据所述充水含水层底板标高建立所述充水含水层底板标高分布模型。
S2、根据预先设定的预测域,建立预测域模型;
具体的,预测域是根据开采掘进方案进行设定的,用户按照开采方案在空间分析软件中绘制预测域的边界线来完成预测域模型的建立。
S3、根据所述充水含水层的分布模型和所述预测域模型,建立预测域静涌水量模型;
可选地,在该实施例中,所述充水含水层的分布模型还包括充水含水层厚度分布模型,所述根据所述充水含水层的分布模型和所述预测域模型,建立预测域静涌水量模型具体包括:
根据所述地下水水位标高和所述充水含水层底板标高,建立充水含水层的充水含水层厚度分布模型;
具体的,充水含水层厚度为水位标高减去底板标高,通过空间分析,将预测前地下水水位标高分布模型中的水位标高属性值减去底板标高分布模型中的底板标高属性值,即可得到充水含水层的厚度分布模型。
根据所述充水含水层的充水含水层厚度分布模型、给水度分布模型和预测域模型,采用空间分析方法建立预测域静涌水量模型。
具体的,经过空间分析后所建立的预测域模型由多个预测域单元组成。
S4、根据所述预测域静涌水量模型,计算矿井静涌水量。
具体的,预测域静涌水量模型的属性值中包含了各个预测域单元的给水度、充水含水层厚度等参数,根据预测域静涌水量模型的属性值即可计算矿井静涌水量。
可选地,在该实施例中,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层中任一层的单位面积静涌水总量,所述根据所述预测域静涌水量模型,计算矿井静涌水量具体包括:
S4.1、从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度和充水含水层厚度;
S4.2、分别根据所述至少一个预测域单元中每一个预测域单元的给水度和充水含水层厚度,计算每一个预测域单元的单位面积静涌水量Q单静ij;
具体的,Q单静ij的计算公式为:
Q单静ij=μijMij (1);
其中,m为预测域静涌水量模型中预测域单元的个数,μij为第i个充水含水层的第j个预测域单元的给水度,Mij为第i个充水含水层的第j个预测域单元的厚度;
S4.3、根据每一个预测域单元的单位面积静涌水量,统计所述预测域中第一充水含水层的单位面积静涌水总量Q单静j,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层;
具体的,Q单静j的计算公式为:
其中,n为开采影响范围内充水含水层的数量。
可选地,在该实施例中,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层的静涌水总量,所述根据所述预测域静涌水量模型,计算矿井静涌水量具体包括:
S4.1、从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度、充水含水层厚度和面积;
S4.2、分别根据所述至少一个预测域单元中每一个预测域单元的给水度、充水含水层厚度和面积,计算每一个预测域单元的静涌水量Q静ij;
具体的,Q静ij的计算公式为:
Q静ij=μijMijSij (3);
其中,μij为第i个充水含水层的第j个预测域单元的给水度,Mij为第i个充水含水层的第j个预测域单元的厚度,Sij为第i个充水含水层的第j个预测域单元的面积;
S4.3、根据每一个预测域单元的静涌水量,统计所述预测域中第一充水含水层的静涌水总量Q静i;
具体的,Q静i的计算公式为:
其中,m为预测域静涌水量模型中预测域单元的个数;
S4.4、统计所述至少一层充水含水层中的每一层充水含水层的静涌水总量之和,作为所述矿井静涌水量Q静,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层;
具体的,Q静的计算公式为:
其中,n为开采影响范围内充水含水层的数量。
本发明提供的一种基于空间分析的矿井静涌水量计算方法,通过对工作面开采影响范围的含水层所赋存地下水的静储量的计算评价,提高了对矿井涌水量的预测结果的准确性。
图2为本发明提供一种基于空间分析的矿井静涌水量计算系统的结构框图,如图2所示,该系统包括:
第一建立模块,用于预先录入的充水含水层的特征参数,建立充水含水层的分布模型;
具体的,各个充水含水层的相关特征参数是通过对开采区域进行钻孔揭露及抽水试验等取得的,并且通过在如Mapgis,Arcgis等空间分析软件中输入特征参数来完成模型的建立。
进一步,所述特征参数包括:给水度、地下水水位标高和充水含水层底板标高,所述充水含水层的分布模型包括给水度分布模型、预测前地下水水位标高分布模型和充水含水层底板标高分布模型;
所述第一建立模块具体包括:
第一建立单元,用于根据所述给水度建立所述给水度分布模型;
第二建立单元,用于根据所述地下水水位标高建立所述预测前地下水水位标高分布模型;
第三建立单元,用于根据所述充水含水层底板标高建立所述充水含水层底板标高分布模型。
第二建立模块,用于根据预先设定的预测域,建立预测域模型;
具体的,预测域是根据开采掘进方案进行设定的,用户按照开采方案在空间分析软件中绘制预测域的边界线来完成预测域模型的建立。
第三建立模块,用于根据所述充水含水层的分布模型和所述预测域模型,建立预测域静涌水量模型;
可选地,在该实施例中,所述充水含水层的分布模型还包括充水含水层厚度分布模型,所述第三建立模块具体包括:
第四建立单元,用于根据所述地下水水位标高和所述充水含水层底板标高,建立充水含水层的充水含水层厚度分布模型;
具体的,充水含水层厚度为水位标高减去底板标高,通过空间分析,将预测前地下水水位标高分布模型中的水位标高属性值减去底板标高分布模型中的底板标高属性值,即可得到充水含水层的厚度分布模型。
第五建立单元,用于根据所述充水含水层的充水含水层厚度分布模型、给水度分布模型和预测域模型,采用空间分析方法建立预测域静涌水量模型。
具体的,经过空间分析后所建立的预测域模型由多个预测域单元组成。
计算模块,用于根据所述预测域静涌水量模型,计算矿井静涌水量。
可选地,在该实施例中,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层中任一层的单位面积静涌水总量,所述计算模块具体包括:
第一获取单元,用于从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度和充水含水层厚度;
第一计算单元,用于按照下列公式(1),分别根据所述至少一个预测域单元中每一个预测域单元的给水度和充水含水层厚度,计算每一个预测域单元的单位面积静涌水量Q单静ij;
具体的,Q单静ij的计算公式为:
Q单静ij=μijMij (1);
其中,m为预测域静涌水量模型中预测域单元的个数,μij为第i个充水含水层的第j个预测域单元的给水度,Mij为第i个充水含水层的第j个预测域单元的厚度;
第二计算单元,用于根据每一个预测域单元的单位面积静涌水量,统计所述预测域中第一充水含水层的单位面积静涌水总量Q单静j,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层。
具体的,Q单静j的计算公式为:
其中,n为开采影响范围内充水含水层的数量。
可选地,在该实施例中,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层的静涌水总量,所述计算模块具体包括:
第二获取单元,用于从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度、充水含水层厚度和面积;
第三计算单元,用于分别根据所述至少一个预测域单元中每一个预测域单元的给水度、充水含水层厚度和面积,计算每一个预测域单元的静涌水量Q静ij;
具体的,Q静ij的计算公式为:
Q静ij=μijMijSij (3);
其中,μij为第i个充水含水层的第j个预测域单元的给水度,Mij为第i个充水含水层的第j个预测域单元的厚度,Sij为第i个充水含水层的第j个预测域单元的面积;
第一统计单元,用于根据每一个预测域单元的静涌水量,统计所述预测域中第一充水含水层的静涌水总量Q静i;
具体的,Q静i的计算公式为:
其中,m为预测域静涌水量模型中预测域单元的个数;
第二统计单元,用于统计所述至少一层充水含水层中的每一层充水含水层的静涌水总量之和,作为所述矿井静涌水量Q静,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层。
具体的,Q静的计算公式为:
其中,n为开采影响范围内充水含水层的数量。
本发明提供的一种基于空间分析的矿井静涌水量计算系统,通过对工作面开采影响范围的含水层所赋存地下水的静储量的计算评价,提高了对矿井涌水量的预测结果的准确性。
实例:某煤矿规划采区静涌水量预测
第一步,某煤矿位于鹤岗矿区中部,矿井主要充水含水层为煤层顶板白垩系下统石头庙子组砾岩裂隙含水层,石头庙子组含水层与开采煤层之间有一层相对隔水的粉砂岩、细砂岩隔水层,倾角15~35°(见图3)。矿井开采方式为分水平开采,主采11#煤层,规划主要开采水平为四水平,开采标高为-600~-900m。
第二步,根据钻孔揭露及抽水试验等取得的各个充水含水层的相关特征参数,在MAPGIS平台建立石头庙子组含水层规划采区的预测前地下水水位标高分布模型,如图4所示。
第三步,因规划采区西边为正在开采矿井,因此规划采区西边含水层厚度为0,直接赋值到预测前地下水水位标高分布模型中。
第四步,对规划采区的预测前地下水水位标高分布模型和充水含水层底板标高分布模型进行空间分析,建立充水含水层的厚度分布模型,如图5所示。
第五步,考虑给水度取经验值(取值0.1),直接在规划采区开采前水位模型中赋值。对规划采区的预测前地下水水位标高分布模型和充水含水层厚度分布模型进行空间分析,建立静涌水量模型。
第六步,根据式(4)可统计规划采区的总静涌水量6618901m3。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于空间分析的矿井静涌水量计算方法,其特征在于,该方法包括以下步骤:
步骤1、根据预先录入的充水含水层的特征参数,建立充水含水层的分布模型;
步骤2、根据预先设定的预测域,建立预测域模型;
步骤3、根据所述充水含水层的分布模型和所述预测域模型,建立预测域静涌水量模型;
步骤4、根据所述预测域静涌水量模型,计算矿井静涌水量。
2.根据权利要求1所述的一种基于空间分析的矿井静涌水量计算方法,其特征在于,所述特征参数包括:给水度、地下水水位标高和充水含水层底板标高,所述充水含水层的分布模型包括给水度分布模型、预测前地下水水位标高分布模型和充水含水层底板标高分布模型;
所述根据充水含水层的特征参数,建立充水含水层的分布模型具体包括:
根据所述给水度建立所述给水度分布模型;
根据所述地下水水位标高建立所述预测前地下水水位标高分布模型;
根据所述充水含水层底板标高建立所述充水含水层底板标高分布模型。
3.根据权利要求2所述的一种基于空间分析的矿井静涌水量计算方法,其特征在于,所述充水含水层的分布模型还包括充水含水层厚度分布模型,所述根据所述充水含水层的分布模型和所述预测域模型,建立预测域静涌水量模型具体包括:
根据所述地下水水位标高和所述充水含水层底板标高,建立充水含水层的充水含水层厚度分布模型;根据所述充水含水层的充水含水层厚度分布模型、给水度分布模型和预测域模型,采用空间分析方法建立预测域静涌水量模型。
4.根据权利要求3所述的一种基于空间分析的矿井静涌水量计算方法,其特征在于,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层中任一层的单位面积静涌水总量,所述根据所述预测域静涌水量模型,计算矿井静涌水量具体包括:
步骤4.1、从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度和充水含水层厚度;
步骤4.2、分别根据所述至少一个预测域单元中每一个预测域单元的给水度和充水含水层厚度,计算每一个预测域单元的单位面积静涌水量;
步骤4.3、根据每一个预测域单元的单位面积静涌水量,统计所述预测域中第一充水含水层的单位面积静涌水总量,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层。
5.根据权利要求3所述的一种基于空间分析的矿井静涌水量计算方法,其特征在于,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层的静涌水总量,所述根据所述预测域静涌水量模型,计算矿井静涌水量具体包括:
步骤4.1、从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度、充水含水层厚度和面积;
步骤4.2、分别根据所述至少一个预测域单元中每一个预测域单元的给水度、充水含水层厚度和面积,计算每一个预测域单元的静涌水量;
步骤4.3、根据每一个预测域单元的静涌水量,统计所述预测域中第一充水含水层的静涌水总量;
步骤4.4、统计所述至少一层充水含水层中的每一层充水含水层的静涌水总量之和,作为所述矿井静涌水量,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层。
6.一种基于空间分析的矿井静涌水量计算系统,其特征在于,该系统包括:
第一建立模块,用于预先录入的充水含水层的特征参数,建立充水含水层的分布模型;
第二建立模块,用于根据预先设定的预测域,建立预测域模型;
第三建立模块,用于根据所述充水含水层的分布模型和所述预测域模型,建立预测域静涌水量模型;
计算模块,用于根据所述预测域静涌水量模型,计算矿井静涌水量。
7.根据权利要求6所述的一种基于空间分析的矿井静涌水量计算系统,其特征在于,所述特征参数包括:给水度、地下水水位标高和充水含水层底板标高,所述充水含水层的分布模型包括给水度分布模型、预测前地下水水位标高分布模型和充水含水层底板标高分布模型;
所述第一建立模块具体包括:
第一建立单元,用于根据所述给水度建立所述给水度分布模型;
第二建立单元,用于根据所述地下水水位标高建立所述预测前地下水水位标高分布模型;
第三建立单元,用于根据所述充水含水层底板标高建立所述充水含水层底板标高分布模型。
8.根据权利要求7所述的一种基于空间分析的矿井静涌水量计算系统,其特征在于,所述充水含水层的分布模型还包括充水含水层厚度分布模型,所述第三建立模块具体包括:
第四建立单元,用于根据所述地下水水位标高和所述充水含水层底板标高,建立充水含水层的充水含水层厚度分布模型;
第五建立单元,用于根据所述充水含水层的充水含水层厚度分布模型、给水度分布模型和预测域模型,采用空间分析方法建立预测域静涌水量模型。
9.根据权利要求8所述的一种基于空间分析的矿井静涌水量计算系统,其特征在于,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层中任一层的单位面积静涌水总量,所述计算模块具体包括:
第一获取单元,用于从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度和充水含水层厚度;
第一计算单元,用于分别根据所述至少一个预测域单元中每一个预测域单元的给水度和充水含水层厚度,计算每一个预测域单元的单位面积静涌水量;
第二计算单元,用于根据每一个预测域单元的单位面积静涌水量,统计所述预测域中第一充水含水层的单位面积静涌水总量,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层。
10.根据权利要求8所述的一种基于空间分析的矿井静涌水量计算系统,其特征在于,所述预测域包含至少一层充水含水层;所述矿井静涌水量为所述所述预测域中所有充水含水层的静涌水总量,所述计算模块具体包括:
第二获取单元,用于从所述预测域静涌水量模型中获取与第一充水含水层对应的至少一个预测域单元的给水度、充水含水层厚度和面积;
第三计算单元,用于分别根据所述至少一个预测域单元中每一个预测域单元的给水度、充水含水层厚度和面积,计算每一个预测域单元的静涌水量;
第一统计单元,用于根据每一个预测域单元的静涌水量,统计所述预测域中第一充水含水层的静涌水总量;
第二统计单元,用于统计所述至少一层充水含水层中的每一层充水含水层的静涌水总量之和,作为所述矿井静涌水量,其中所述第一充水含水层为所述至少一层充水含水层中的任意充水含水层。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710948510.5A CN107886185A (zh) | 2017-10-12 | 2017-10-12 | 一种基于空间分析的矿井静涌水量计算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710948510.5A CN107886185A (zh) | 2017-10-12 | 2017-10-12 | 一种基于空间分析的矿井静涌水量计算方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107886185A true CN107886185A (zh) | 2018-04-06 |
Family
ID=61781456
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710948510.5A Pending CN107886185A (zh) | 2017-10-12 | 2017-10-12 | 一种基于空间分析的矿井静涌水量计算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107886185A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108915667A (zh) * | 2018-06-27 | 2018-11-30 | 中煤科工集团西安研究院有限公司 | 基于渗透系数随分层开采动态变化的矿井涌水量预测方法 |
CN108955649A (zh) * | 2018-05-21 | 2018-12-07 | 中国矿业大学 | 一种煤矿工作面顶板过程涌水量的预测方法 |
CN109166633A (zh) * | 2018-08-28 | 2019-01-08 | 中国煤炭地质总局勘查研究总院 | 矿井涌水量计算方法及装置 |
CN109377042A (zh) * | 2018-10-17 | 2019-02-22 | 中国煤炭地质总局水文地质局 | 一种煤层顶板突水危险性评价方法 |
CN110580655A (zh) * | 2019-09-17 | 2019-12-17 | 山东浪潮通软信息科技有限公司 | 一种基于lstm算法的煤矿涌水量预测的方法和系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105260537A (zh) * | 2015-10-12 | 2016-01-20 | 内蒙古上海庙矿业有限责任公司 | 裂隙型含水层富水性分析方法及工作面涌水量预计方法 |
-
2017
- 2017-10-12 CN CN201710948510.5A patent/CN107886185A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105260537A (zh) * | 2015-10-12 | 2016-01-20 | 内蒙古上海庙矿业有限责任公司 | 裂隙型含水层富水性分析方法及工作面涌水量预计方法 |
Non-Patent Citations (2)
Title |
---|
吴铁卫,张冀: "层次分析模型在"三图-双预测"法中的应用--以鹤岗矿区新陆煤矿为例", 《中国煤炭地质》 * |
王春,李春忠: "《深基坑工程降水技术研究与实践》", 1 September 2016 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108955649A (zh) * | 2018-05-21 | 2018-12-07 | 中国矿业大学 | 一种煤矿工作面顶板过程涌水量的预测方法 |
CN108955649B (zh) * | 2018-05-21 | 2020-06-16 | 中国矿业大学 | 一种煤矿工作面顶板过程涌水量的预测方法 |
CN108915667A (zh) * | 2018-06-27 | 2018-11-30 | 中煤科工集团西安研究院有限公司 | 基于渗透系数随分层开采动态变化的矿井涌水量预测方法 |
CN108915667B (zh) * | 2018-06-27 | 2019-09-24 | 中煤科工集团西安研究院有限公司 | 基于渗透系数随分层开采动态变化的矿井涌水量预测方法 |
CN109166633A (zh) * | 2018-08-28 | 2019-01-08 | 中国煤炭地质总局勘查研究总院 | 矿井涌水量计算方法及装置 |
CN109166633B (zh) * | 2018-08-28 | 2020-10-16 | 中国煤炭地质总局勘查研究总院 | 矿井涌水量计算方法及装置 |
CN109377042A (zh) * | 2018-10-17 | 2019-02-22 | 中国煤炭地质总局水文地质局 | 一种煤层顶板突水危险性评价方法 |
CN110580655A (zh) * | 2019-09-17 | 2019-12-17 | 山东浪潮通软信息科技有限公司 | 一种基于lstm算法的煤矿涌水量预测的方法和系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107886185A (zh) | 一种基于空间分析的矿井静涌水量计算方法及系统 | |
CN106703883B (zh) | 一种个性化确定采煤工作面底板突水危险等级的方法 | |
CN109577982A (zh) | 壁式连采连充保水采煤及水资源运移监测、水害预警方法 | |
CN103323887B (zh) | 一种煤矿采动稳定区煤层气储量的评估方法及系统 | |
CN103226732B (zh) | 一种基于gms的矿区不同开采中段的地下水渗流场预测方法 | |
CN106032754B (zh) | 一种基于地下水流速流向测定的煤矿开采防治水方法 | |
CN110162922A (zh) | 一种水驱油藏优势渗流通道的综合识别方法 | |
CN111691872B (zh) | 一种煤层覆岩裂隙发育范围及程度的动态精准探测方法 | |
CN107748928A (zh) | 一种基于空间分析的矿井动涌水量计算方法及系统 | |
CN104899358B (zh) | 奥灰岩溶裂隙水网络横向分布的预测方法 | |
CN111794740A (zh) | 一种适用在缝洞型碳酸盐岩油藏动态储量计算的方法 | |
CN103778480A (zh) | 一种基于敏感度分析的裂隙带高度预测方法 | |
CN111379562B (zh) | 一种复合水体下的控水采煤方法及装置 | |
CN114004178A (zh) | 一种单元流量分析方法及系统 | |
CN113375753B (zh) | 一种煤矿综采工作面采动对地下水影响的监测与分析方法 | |
CN110567854A (zh) | 一种基于矿井闭井后相邻矿井水文地质模型的涌水量预计和计算方法 | |
CN102155213B (zh) | 矿井采动裂隙的动态检测方法 | |
CN111119992B (zh) | 一种煤层顶板疏放水钻孔参数确定方法 | |
Light et al. | Mine-water flow between contiguous flooded underground coal mines with hydraulically compromised barriers | |
CN112948924A (zh) | 基于覆岩结构的近松散层开采导水裂隙带高度确定方法 | |
CN107542454A (zh) | 一种采场围岩裂隙发育程度评价方法 | |
CN116756807A (zh) | 耦合施工过程的线性地下工程涌水量数值预测方法 | |
CN110700884B (zh) | 一种基于解析法比拟预测综放工作面涌水量的方法 | |
Martinez et al. | Use of numerical groundwater modelling for mine dewatering assessment | |
CN114444266A (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 |