CN104346512A

CN104346512A - GIS (geographic information system)-based grid method coal resource/reserves estimation method

Info

Publication number: CN104346512A
Application number: CN201310341399.5A
Authority: CN
Inventors: 胡宝林; 郑凯歌; 程乔; 刘兵昌; 杨华泽
Original assignee: Anhui University of Science and Technology
Current assignee: Anhui University of Science and Technology
Priority date: 2013-08-01
Filing date: 2013-08-01
Publication date: 2015-02-11

Abstract

The invention discloses a method for reserves computation, which is realized through the following steps: based on the grid method, and through the adoption of the Kriging method and the surfer gridded data, contour line data is generated through the data conversion function, computation parameters are assigned to all grid cells through the mapgis software, and the spatial overlay technology is adopted for the overlapping with the computed grid range. The method overcomes the problem that the selection result of reserves computation parameter control data, such as coal-seam thickness and rock stratum dip angle, is affected by human factors.

Description

Based on the gridding method coal resources/estimation method of reserve of GIS

Art

Patent of the present invention relates to and space overlapping technology being combined with GIS grid method and the method that calculates of a kind of coal resources reserves carried out.

Background technology

Mineral resources Resource/Reserves computing method refer to the estimation carried out underground mine according to the geological information obtained in geology of mineral deposit exploration (comprising exploitation) process ^[1].It is an element task in coal exploitation process, and its result reliability how, is directly connected to the rationality of colliery investment risk size, coal production decision-making, and the validity of routine work.At present in inclination and gently inclined seam, coal resources reserves calculate the computing method that to be the methods such as block section method, gridding method be all based on paper drawing adopted.Although the utilization of computing machine, add the precision of Resource/Reserves reference area, but because boring point is always in a discrete distribution in plane, to in its piecemeal manufacturing process, need to carry out interpolation to raw data, this has a lot of human factor with regard to the selection caused for the Parameter of reserve volume calculation such as thickness of coal seam and formation dip control data, and the problem that result of calculation is larger with real data phase ratio error is more and more outstanding ^[2-3].

Citing document

[1] Wang Zhengshuai, Gu Hehe. Digital Mine resource reserve computing method [J]. coal field geology and exploration, 2010,38 (3): 6-7.

[2] Ma Hongbin, Xiong Junnan. based on the reserve estimate system [J] of geostatistics. coal journal, 2007 (3); 18-22.

[3] Li Cuiping, Li Zhongxue, Hao Xiuqiang, Hao Jinhui. in Geological Engineering visual simulating, the calculating of grade and reserves realizes [J]. University of Science & Technology, Beijing's journal, 2007 (9).

Summary of the invention

The selection of getting it right in the Parameter of reserve volume calculation such as thickness of coal seam and formation dip control data to overcome existing coal resources finimeter has a lot of human factor, the problem that error calculated is larger.Patent of the present invention provides a kind of based on gridding method, utilizes and can solve based on GIS space overlapping technology the error that human factor causes preferably, improve computational accuracy.

Patent of the present invention solves the technical scheme that its technical matters adopts: be first optimized Parameter of reserve volume calculation, 1, thickness of coal seam: survey coal dot thickness data according to each boring and down-hole in bore database, generate thickness of coal seam data file [.dat], and in surfer software based on Kriging method to data file gridding process, save as Grd file, in gridding process, mesh spacing is chosen with the rectangular node calculated onesize, with the grid cell one_to_one corresponding calculated, the data transformations function of surfer software is utilized to be [.dat] data by Grd file translations, and process data into excel form, in MapGis software, utilize No. ID to mount to each grid cell by thick for coal parameter, 2, seam inclination: in GIS software, isocontour for every bar line attribute is composed as actual level line elevation.Utilize the DTM analysis module line drawing point function of spatial analysis to extract point height, save as TmpGrid.Grd file (see Fig. 2).Utilizing the function of surfur infinitesimal analysis module mesorelief modeling to obtain inclination angle Grd file (see Fig. 3), obtaining seam inclination [.dat] data by transforming.The parameter acquiring thick with coal is similar, uses the gridding of suffer software and transformation function to generate [.dat] data file, and processes data into excel form, utilize No. No. ID to mount to each grid cell by dip angle parameter by mapgis software; 3, area: the graticule mesh of the CGIS form of foundation is converted into mapgis form, accurately obtains the area of each graticule mesh.Secondly, Resource/Reserves calculates, and 1, utilize the boundary lines such as boundary line, border, field with "nine squares", coal class boundary line, reserves type boundary line to split, and first left and then right according to block segment number on the basis of original piece of segment number, principle divided block segment number from shallow to deep, the block segment number in each coal seam is independent separately; 2, when calculating runs into trap-up, coil, in space plane projection, there is additivity because the special nature of trap-up determines tomography two, this just makes the superposition will considering upper and lower dish reserves in the calculating of reserves, the reserves double counting namely between the upper and lower intersection line of broken coal seam of tomography; 3, original coal seam in the reserves figure of coal seam is utilized to adopt boundary line, to carry out calculating contrast to distinct methods.4, the block segment limit that reserves block section method calculates is carried out space overlapping with the rectangular node comprising coal seam computer capacity by MapInfo software, initial parameter and Optimal Parameters are assigned to each grid cell and utilize reserves computing formula to carry out calculating and contrasting.

The beneficial effect of patent of the present invention is, this method improves reserves computational accuracy compared with the methods such as geological block, solves the problem that reserves computational accuracy is on the low side, decreases the unnecessary error that artificial origin causes; Inherit the advantage of block section method arithmetic mean, have certain originality, contrast with production real data, result of calculation is objective vivid, may be used for spot-check the accuracy that reserves examine the result of calculation verifying report with reserves.

Accompanying drawing explanation

Below in conjunction with accompanying drawing table and embodiment, patent of the present invention is further illustrated.

Fig. 1 is raster data superposition schematic diagram

Fig. 2 is coal thick data gridding process schematic diagram

Fig. 3 is coal thick Grd file translations data function schematic diagram

Fig. 4 is the instruction of the generative process to coal seam spot elevation.

Fig. 5 is the instruction to seam inclination data generating procedure.

Fig. 6 is grid unit attribute schematic diagram.

Fig. 7 can adopt boundary block section schematic diagram.

Fig. 8 is reserves summarized results contrast table.

Fig. 9 is that contrast table calculated by partial block section finimeter.

embodiment

Fig. 1 be to GIS space overlapping function simple displaying.

Fig. 2 and Fig. 3 carries out process of refinement to the data gridding function and grid conversion function that utilize suffer software to thickness of coal seam, and in conjunction with the attribute connecting function of MapGis software, utilizes No. ID by thick for the coal signal being assigned to grid cell.

Fig. 4 and Fig. 5 is to utilizing GIS software that actual level line elevation is assigned to each level line.Utilize the DTM analysis module line drawing point function of MapGis software space analytic function to extract point height, save as TmpGrid.Grd file.Utilize the function of surfur infinitesimal analysis module mesorelief modeling to obtain inclination angle Grd file, data gridding function and the grid conversion function of recycling suffer software carry out process of refinement to thickness of coal seam, utilize No. ID by thick for the coal signal being assigned to grid cell.

Fig. 6 is the schematic diagram processed computer capacity border and tomography etc.

Fig. 7 calculates at parameters and utilizes MapGis software to be assigned to attribute function parameters to be assigned to each unit grid and the result obtained, final derivation attribute form, treated as Fig. 9, utilize Fig. 9 data to carry out calculating contrast to each computing method, obtain Fig. 8 data.

Claims

1. one kind based on GIS gridding method, surfer gridded data is utilized by means of Kriging method, and maintenance data transformation function generates isoline data, by mapgis software, calculating parameter is assigned to each grid cell, and superposes the method for carrying out reserves calculating in conjunction with space overlapping technology with the graticule mesh scope of calculating.

2. Reserve Estimation Method according to claim 1, it is characterized in that: computer capacity is divided into small square net, utilize GIS spatial overlay analysis function the layer of vector data structure to be carried out the space overlapping in territory and territory and the layer space overlapping of raster data structure, the Reserve Estimation Method of function to routine such as data gridding and data conversion in conjunction with suffer software improves.