ISPRS Workshop on Service and Application of Spatial Data Infrastructure, XXXVI(4/W6), Oct. 14-16, Hangzhou, China 
VISUALIZED APPLICATION OF GIS IN MINING SUBSIDENCE ANALYSIS 
Huayang Dai a , Zuorua, Yin b , Zhiyi Deng a,c 
‘'China University of Mining and Technology, Beijing campus, Beijing, 100083, China - daihuayang@263.net 
b Kailuan (Group) Ltd. Liability Corporation , Tangshan, 063000, China - yzru@kailuan.com.cn 
c Kailuan (Group) Ltd. Liability Corporation , Tangshan, 063000, China - dengzy@kailuan.com.cn 
KEY WORDS: Visualized Application with GIS, Mining Subsidence Analysis System, Maplnfo/ MapBasic, Seam Inclination- 
based Subsidence Model 
ABSTRACT: 
To determine the degree of mining affection on surface or in strata, the GIS technique is used, with Mapinfo as the platform. 
Through programming of seam inclination-based subsidence mathematical models and mining disturbance analysis visualization, 
Mapbasic/Vc++/Vfortran programming design is realized under Windows. The program system provides not only an effective 
approach for model analysis, but also an efficient visual tool and interface for subsidence model application, laying therefore a solid 
foundation for widespread use of the model. 
1. INTRODUCTION ON MINING SUBSIDENCE 
MODEL 
1.1 General Instructions 
A new subsidence model is put forward to calculate surface 
movement and deformation induced by coal seam mining[l,2]. 
The approach features that: the seam inclination angle, as the 
most prominent influential factor, among others, governing the 
pattern of ground movement, is reckoned as a principal 
parameter while the individual elements in working seams with 
any dip angles are considered to be surface micro-element 
vectors which are then resolved according to the parallelogram 
law; introduction of mining influential functions respectively 
for vertical and horizontal components of surface micro-element, 
based on random medium theory; introduction of rock 
movement parameters irrespective of the occurrence pattern 
and size of workings; determination of the unit movement 
trough through integrating both vertical and horizontal 
components according to equivalent influence principle; 
development of a continuous mining subsidence mathematical 
model through surface integral of influential functions over the 
entire area above the workings, according to the principle of 
superimposition. A totally new theory and method, so-called 
vector prediction method that is applicable for seams with dip 
angles ranging 0-90° is thus developed through combining the 
seam dip angle based maximum subsidence calculation model, 
mining influence propagation angle calculation model as well as 
the design routine of these models' programmed algorithms. 
W eh (x,y) 
x 2 + (,y-Hctg0 o ) 2 
cos a * r i 
— 
(1) 
dSv is the vertical mining unit with an area of l*sina and a 
thickness of 1. It can also be reckoned as a horizontal mining 
unit with an area of 1><1 and a thickness of sina. Similarly, the 
expression for the subsidence basin for the vertical component 
dSv can be obtained as follows: 
x 2 +(y-Hctg0 o ) 2 
X , 
(2) 
sin Cl 
W ev (x,y) = — r e 
In Equations (1) and (2): 
Weh(x, y), Wev (x, y) - Mining subsidence basins of dSh and 
dSv respectively 
rh, rv - Main influence radii of mining of flat and vertical 
seams respectively 
H - Mining depth of dS at point P 
00 - Mining influence propagation angle 
By superimposition of the subsidence basins of the horizontal 
component and vertical component, we obtain the expression 
for the subsidence basin for dS: 
W e {x,y) = W eh {x,y)+W ev {x,y) (3) 
Considering the orthogonality of the micro-unit vector’s two 
basic components of dSh and dSv, we have: 
f 
I \[w eh {x,y)dS v =0 
and hence 
V' (w) = W mn iH s WJx,y)dS h + fjW tv (x,y)dSv (4) 
As can be seen from Equation (4), the rock and surface 
movements caused by mining of a pitching seam may be 
attributed to two mining influences. One is the influence 
produced by the horizontal component and the other by the 
vertical component. The superimposed area of the horizontal 
mining component is equal to the projected area of the working 
face on the horizontal plane; while the superimposed area of the 
vertical component is equal to the projected area of the working 
face on the vertical plane. Therefore, the horizontal and vertical 
projections can respectively be taken as the field of integration 
in the integration operation. The integrand is the function in the 
domain of area of S with varying Z coordinates. The 
significance of Equation (4) lies in that it enables the integration 
over the domain of area S to be converted to integration over its 
vertical and horizontal projected areas. This lays the foundation 
for the design of the algorithm of the vector method. 
According to the basic theory of ground movement, there 
exists a functional relation between ground movement and the