基于GIS的网格法煤炭资源/储量估算方法技术

一种以网格法为基础，借助于克里格法利用surfer网格化数据，并运用数据转化功能生成等值线数据，通过mapgis软件将计算参数赋给各个网格单元，并结合空间叠加技术与计算的格网范围叠加进行储量计算的方法。以克服对于煤层厚度及岩层倾角等储量计算参数控制数据的选择具有很多人为因素造成的问题。

【技术实现步骤摘要】
基于GIS的网格法煤炭资源/储量估算方法 所属
本专利技术专利涉及将空间叠加技术与GIS格网法相结合而进行的一种煤炭资源储量计算的方法。
技术介绍
矿产资源资源/储量计算方法是指依据矿床地质勘探(包括开采)过程中所获得的地质信息对地下矿藏进行的估算[1]。它是煤炭开发过程中的一项基础工作，其结果可靠性如何，直接关系到煤矿投资风险大小、煤矿生产决策的合理性，以及日常工作的有效性。目前在倾斜及缓倾斜煤层中，煤炭资源储量计算采用的是块段法、网格法等方法均是基于纸质图纸的计算方法。虽然计算机的运用，增加了资源/储量计算面积的精度，但由于钻孔点总在平面呈离散分布，在对其分块制作过程中，需要对原始数据进行插值，这就造成了对于煤层厚度及岩层倾角等储量计算参数控制数据的选择具有很多人为因素，计算结果与实际数据相比误差较大的问题越来越突出[2_3]。 引用文献 [I]王正帅，顾和和.数字化矿山资源储量计算方法[J].煤田地质与勘探，2010，38(3):6-7. [2]马洪滨，熊俊楠.基于地质统计学的储量估算系统[J].煤炭学报，2007(3)；18-22. [3]李翠平，李仲学，郝秀强，郝晋会.地质工程可视化仿真中品位与储量的计算实现[J].北京科技大学学报，2007 (9).
技术实现思路
为了克服现有煤炭资源储量计算对于煤层厚度及岩层倾角等储量计算参数控制数据的选择具有很多人为因素，计算结果误差较大的问题。本专利技术专利提供一种以网格法为基础，利用基于GIS空间叠加技术可以较好的解决人为因素造成的误差，提高计算精度。 本专利技术专利解决其技术问题所采用的技术方案是:首先对储量计算参数进行优化，1、煤层厚度:根据钻孔数据库中各个钻孔及井下测煤点厚度资料，生成煤层厚度数据文件[.dat]，并在surfer软件中基于克里格法对数据文件网格化处理，保存为Grd文件，网格化过程中格网间距选取与计算的矩形网格同样大小，与计算的网格单元一一对应，利用surfer软件的数据转化功能将Grd文件转化为[.dat]数据，并将数据处理为excel格式，在MapGis软件中利用ID号将煤厚参数挂接给各个网格单元；2、煤层倾角:在GIS软件中将每条等高线的线属性赋为实际等高线高程。利用空间分析的DTM分析模块线提取点功能提取点高程，保存为TmpGrid.Grd文件(见图2)。利用surfur微积分模块中地形建模的功能得到倾角Grd文件(见图3)，通过转化获得煤层倾角[.dat]数据。同煤厚的参数获取相似，运用suffer软件的网格化和转化功能生成[.dat]数据文件，并将数据处理为excel格式，利用ID号号将倾角参数通过mapgis软件挂接给各个网格单元；3、面积:将建立的CGIS格式的格网转化为mapgis格式，精确求出各个格网的面积。其次，资源/储量计算，1、利用井田边界界线、煤类界线、储量类型界线等界线进行分割，并在原有块段编号的基础上按照块段编号先左后右，由浅到深的原则划分块段编号，各煤层的块段编号各自独立；2、当计算遇到逆断层时，由于逆断层的特殊性质决定了断层两盘在空间平面投影上具有叠加性，这就使得储量的计算中要考虑上、下盘储量的叠加，即断层上、下断煤交线之间的储量重复计算；3、利用煤层储量图中原有煤层可采边界线，以便对不同方法进行计算对比。4、将储量块段法计算的块段范围与包含煤层计算范围的矩形网格借助MapInfo软件进行空间叠加，将原始参数和优化参数赋给各个网格单元利用储量计算公式进行计算和对比。 本专利技术专利的有益效果是，该法与地质块段等方法相比提高了储量计算精度，解决了储量计算精度偏低的问题，减少了人为原因造成的不必要误差；继承了块段法算术平均的优点，具有一定的独创性，与生产实际数据对比，计算结果客观真切，可以用于抽查储量核实和储量核查报告的计算结果的准确性。 【附图说明】 下面结合附图表和实施例对本专利技术专利进一步说明。 图1是栅格数据叠加示意图 图2是煤厚数据网格化处理示意图 图3是煤厚Grd文件转化数据功能示意图 图4是对煤层高程点的生成过程的指示。 图5是对煤层倾角数据生成过程的指示。 图6是格网单元属性示意图。 图7是可采边界块段示意图。 图8是储量汇总结果对比表。 图9是部分块段储量计算对比表。 【具体实施方式】 图1是对GIS空间叠加功能的的简单展示。 图2和图3是对利用suffer软件的数据网格化功能和网格转换功能对煤层厚度进行精细化处理，并结合MapGis软件的属性挂接功能，利用ID号将煤厚赋给网格单元的示 O 图4和图5是对利用GIS软件将实际等高线高程赋给各个等高线。利用MapGis软件空间分析功能的DTM分析模块线提取点功能提取点高程，保存为TmpGrid.Grd文件。利用surfur微积分模块中地形建模的功能得到倾角Grd文件,再利用suffer软件的数据网格化功能和网格转换功能对煤层厚度进行精细化处理，利用ID号将煤厚赋给网格单元的示意。 图6是对计算范围边界以及断层等处理的示意图。 图7是在各个参数计算并利用MapGis软件赋给属性功能将各个参数赋给每个单元网格而得到的结果，最终导出属性表格，经处理如图9，利用图9数据对各个计算方法进行计算对比，得到图8数据。本文档来自技高网...

【技术保护点】
一种以GIS网格法为基础，借助于克里格法利用surfer网格化数据，并运用数据转化功能生成等值线数据，通过mapgis软件将计算参数赋给各个网格单元，并结合空间叠加技术与计算的格网范围叠加进行储量计算的方法。

【技术特征摘要】
1.一种以GIS网格法为基础，借助于克里格法利用surfer网格化数据，并运用数据转化功能生成等值线数据，通过mapgis软件将计算参数赋给各个网格单元，并结合空间叠加技术与计算的格网范围叠加进行储量计算的方法。2.根据权利要求1所述...

【专利技术属性】
技术研发人员：胡宝林，郑凯歌，程乔，刘兵昌，杨华泽，
申请(专利权)人：安徽理工大学，
类型：发明
国别省市：安徽;34