SOBZYCK - New Technological Solutions in Underground Mining
SOBZYCK - New Technological Solutions in Underground Mining
SOBZYCK - New Technological Solutions in Underground Mining
=
tg l
x 2
y
;
tg l
x 2
y
3
2
2
2
4
2
1
1
(1)
According to (1), we can calculate load P
1
on first crosshole pillar:
( )( )
(
(
=
3
l
l l 4 l l l
tg l 2
P
3
4
2
1 4 3 1 4
3
1
(2)
where average volumetric weight of rocks inside the arch of natural balance.
The maximum load P
2
on crosshole pillar is defined similarly. When it is situated in the middle
of distance l
3
, in account that
=
tg 2
l
h
3
3
, it will be:
( )
(
(
=
3
l
l l l
tg 2
P
2
4
4 1 3 2
(3)
Formulae (2) and (3) are the basis for further choice of pillar widths.
Load P
3
on the basic pillar is defined by the formula:
|
.
|
\
|
+
=
3
l
l
tg 2
l
P
3
2
3
3
(4)
Obvious direction towards the reduction in losses of coal in pillars is the reduction of their
width (l
1
and l
2
) till the minimum permissible sizes. As to crosshole pillars parameters I proposed
the following effective solution. The pillar nearest to drilling hole should keep stability during the
whole period of the drilling works. Apart from that it is in a relaxing zone of the arch of natural
balance of the seam roof. As soon as the production face will depart from pillar by two three
drilling holes, the necessity to keep it steady disappears. Thus, crosshole pillar carries out the func-
tions of protection for drilling hole from jamming. Here is the zone of smallest load.
As to the basic pillar, we can get such distance l
3
between basic pillars, when losses of coal V
will be the smallest. Taking into account that coal losses in pillars depend on widths l
1
and l
2
,
which are calculated from geotechnical characteristics, we can get the function of the purpose:
( )
( ) min , D , m , l , , , r , l , l
l l m
V
4 6 5
3 2
o
+
(5)
where D agency diameter of an auger installation, o coal strength, r radius of the preparatory
working, l
5
and l
6
the length of a drill hole and the length of a basic pillar.
Expressing width of the pillar through symbols l
1
and l
2
, it is possible to simplify the function
of the purpose (5):
( )
( )
min
m l l
1 n V V nV
3 2
3 2 1
+
+ + +
(6)
where V
1
, V
2
and V
3
losses of coal in crosshole pillars, basic pillars, and in the drill hole (usually
a pack of coal is left in the roof of the seam), n number of crosshole pillars.
Copyright 2006 Taylor & Francis Group plc, London, UK
30
For definition of parameters V
1
, V
2
, V
3
and n the following expressions are received:
( )
( )
=
t
=
t
+ =
t
+ =
D 3 l
D 3 l
n
; D
2
Dm 2 V
;
4
D
mD D m 2 l V
;
4
D
mD D m 2 l V
1
3
2
3
2
2 2
2
1 1
(7)
Solving in common the equations (7) and (6), we can come to the criterion of the minimum coal
loss:
min
m
D
4
3
D 3
m
D
2 l
m
D
4
3
D 3
m
D
2 l
D 3 l
D 3 l
l l
1
2
2
2
1
1
3
3 2
t + |
.
|
\
|
+
(
(
t + |
.
|
\
|
+
(8)
To determine the rational value l
3
of the distances between basic pillars according to the formu-
la (8), it is necessary to define widths l
1
and l
2
, and then to search (8) for the minimum. Therefore
the following steps are the calculation of parameters of pillars according to their functional purpose.
Among the existing methods of such calculations there are two basic ways: according to the
admitted stress equal to the strength on compression o
cx
of the pillar material; and according to the
admitted load P
nec
, which is determined from the stress condition of the pillar.
In the first way the estimation of pillar stability is done on the basis of comparison of some ave-
rage (in square of pillar) tensions o with durability of the pillar material in compression o
cx
along
one axis. This parameter is determined from laboratory samples of the pillar material. Such me-
thodology is chosen because defining mechanical properties of pillars directly in a mine is rather
inconvenient. It is necessary to take into account the scale factor, which mentioned various geo-
logical and technological discontinuities. It is necessary to enter correction factors into the results
of laboratory tests. The most important of them is the form factor K
]
. It shows the relationship be-
tween the strength of a pillar in a cubic and other than cubic (prismatic) form. For the crosshole pil-
lar with the proportion of pillar height, equal to the thickness m of the seam, to its width l
1
in the
most typical range:
, 3
l
m
1
1
s s
it is possible with sufficient accuracy to recommend Erofeeva's formula:
1
]
l
m
15 , 0 15 , 1 K = (9)
The specified range allows estimating pillar parameters from the first method, i.e. to estimate
its stability from the strength of coal in uniaxial compression o
cx.y
. It is necessary to take into ac-
count that undermining breaks coal seam, that's why it is necessary to do the calculations from resi-
dual uniaxial compression strength of coal o
o
cx.y
. Its value is determined in laboratory conditions.
Basic pillars, bearing the load from most part of rock in the arch of natural balance, have much
more width l
2
, which frequently exceeds pillar height. So, for basic pillars the ratio m/l
2
< 1 is fair.
A complex stress condition exists in the pillar [Ruppenejt et al. 1960] so it is inadmissible to esti-
mate it from the first method. Stability of the pillars is expediently assessed by the second method.
Actual load on a pillar is compared with its bearing capacity P
nec
. Bearing capacity of a pillar is de-
fined on the basis of the theory of limited balance. It means that pillar material passes into a con-
dition of limited balance [3].
Copyright 2006 Taylor & Francis Group plc, London, UK
31
Thus pillar-bearing capacity is determined (in our case basic) in view of rectilinear bending around
the Mohre circles:
(
|
.
|
\
|
+ =
y
2
2 nec
,
m
l
F 2 cl 2 P ,
where c and
y
cohesion and angle of internal friction of coal (cohesion is defined in view of the
partial crushing of the seam), F(l
2
/m,
y
) the tabulated function, reflecting a high stress condition
in pillar; it is defined in [3].
According to the above-stated, it is possible to write down two conditions for the choice of cross-
hole pillar width l
1
. They determine the minimum value needed for pillar stability l
1min
near drilling
hole; and the maximum value of pillar width l
1max
, when it is in the middle of the distance l
3
.
Then it is possible to write down a system of inequalities:
o >
o s
max 1
o
y . cx ] 2
min 1
o
y . cx ] 1
l K P
; l K P
(10)
The first inequality (10) of the system is the determining. Having solved it together with the
formula (2), we will get
( )( )
(
(
+
o
>
3
l
l l 4 l l l
tg l K 2
l
3
4
2
min 1 4 3 min 1 4
o
y . cx 3 ]
min 1
(11)
The second inequality is for verifying, whether finding out the minimum pillar width l
1min
exceeds such value, in which crosshole pillar remains stable being in the middle of the distance l
3
.
It is practically impossible to obtain from the analysis of the settlement circuit such variant (when
l
1max
< l
1min
), as load P
1
will always be smaller than load P
2
on the crosshole pillar in the middle. Ne-
vertheless it is interesting to see what safety factor can be chosen near a crosshole pillar so that in
the middle of the distance l
3
all the same occurred it pressing. For this purpose we shall determine
the maximum pillar width l
1max
by solving the second inequality of system (10) and expression (3):
( )
3
o
y . cx ]
o
y . cx ]
4
3 4
max 1
l K tg 2 K tg 2
3
l
l l
l
+ o o
|
.
|
\
|
< (12)
Both equations (11) and (12) are transcendental and their solutions are found only approximate-
ly by numerical methods. For convenience, the calculation of parameters l
1min
and l
1max
is conduc-
ted on the basis of the numerical analysis of the equations (11) and (12). Corresponding monogra-
ms were created [Bondarenko et al. 1999].
Smallest allowed width of a basic pillar l
2
in condition of its stability is determined by the equa-
tion (13):
|
.
|
\
|
+
|
|
.
|
\
|
>
tg 2
l
,
m
l
F 2 c 2
tg 3
l
h l
l
3
y
2
3
4 3
2
(13)
where h
4
height of the column of undermined rock above the basic pillar.
However equation (13) is transcendental, in an obvious kind is not solved, and width l
2
we can
determined by monogram [4].
On Figure 2 the schedule of dependences of smallest allowed width of crosshole pillar l
1min
on
the distance between basic pillars l
3
is given for different mechanical characteristics of rock.
Copyright 2006 Taylor & Francis Group plc, London, UK
32
1
MHa 6 tg
o
y . cx
= o
; 2
MHa 5 tg
o
y . cx
= o
; 3
MHa 4 tg
o
y . cx
= o
;
4
MHa 3 tg
o
y . cx
= o
; 5
MHa 2 tg
o
y . cx
= o
; 6
MHa 1 tg
o
y . cx
= o
Figure 2. Dependence of minimum width l
1min
of crosshole pillar from distance between basic pillars l
3
Pillar width l
1min
changes depending on the distance between basic pillars l
3
and for l
3
s 20 m it
changes from 0,55 m to 1,1 m at parameter o
o
cx.y
tg s 4 MPa (Fig. 2). An increase in the distance
between basic pillars l
3
leads to essential losses of coal in crosshole pillars. From this point of view
it is expedient to suppose that the distance between basic pillars l
3
will be 1820 m (that corresponds
to drilling 910 holes between them). The analysis of dependence l
1min
from l
3
confirms the expe-
diency of cutting crosshole pillars of variable width as a multiplication of coal drilling between
neighboring basic pillars. It allows reducing losses of coal together with the preservation of the
functions stipulated by the technology.
On Figure 3 dependence of smallest allowed width of basic pillar l
2
on the distance between
them l
3
is shown in condition of variable seam thickness m with other parameters fixed. Connec-
tion l
2
with l
3
has non-linear character. So, at 5 m s l
3
s 15 m there is a substantial growth of the re-
quired width of the basic pillar l
2
(from 0,95 to 1,3 m at l
3
= 5 m to 1,31,85 m at l
3
= 15 m at seam
thickness variation from 0,7 m to 1,2 m). At further increase of the distance l
3
(l
3
> 15 m) the cur-
ves occur more flat l
2
(l
3
) and increase of width of basic pillar l
2
occurs less intensely, and after l
3
=
= 20 m within the limits of several percent.
The described law says that with an increase in width l
2
of basic pillar its bearing capacity shar-
ply grows. Intensity of width growth l
2
at increase in distance l
3
depends on mechanical properties
of coal and roof rocks.
Thus, it is established, that rational value of smallest allowed width of crosshole pillars is change-
able l
3
between basic pillars. Crosshole pillars of variable widths allow reducing losses of coal at
auger mining. Also, the dependence of losses of coal from distance l
3
between basic pillars was re-
vealed. Thus the presence of minimum coal losses was established at certain values l
3
. Depending
on geometrical parameters of auger technology and mechanical characteristics of coal and roof rocks
rational value of the distance between basic pillars changes within limits from 10 to 20 m.
Copyright 2006 Taylor & Francis Group plc, London, UK
3
3
Figure 3. Dependence of smallest allowed width l
2
basic pillar from distance between them
1 m = 0,7 m, 2 m = 0,8 m, 3 m = 0,9 m, 4 m = 1,0 m, 5 m = 1,2 m
Copyright 2006 Taylor & Francis Group plc, London, UK
34
REFERENCES
Zborshchik M.P., Nazimko V.V.: Deep Mines Workings Protection in Zones of Unloading. KYIV, Technicka
1991, p. 247.
Zborshchik M.P.: Deep Mines Workings Protection in Goaf. KYIV, Technicka 1978, p. 176.
Ruppenejt K.V., Liberman Yu.M.: Introduction in Mechanics of Rocks. M.: Gosgortehizdat 1960, p. 356.
Bondarenko V.I., Kovalevska I.A., Bugajev I.V., Simanovich G.A.: Reduction in Losses of Coal at Augering
Mining of Thin Coal Seams. The Scientific Bulletin of National Mining Academy of Ukraine 1999, X 4,
p. 3639
Copyright 2006 Taylor & Francis Group plc, London, UK
35
Geomechanics: History, Modern State
and Prospects of Development
Olexandr Shashenko
National Mining University. Dnipropetrovsk, Ukraine
Tadeusz Majcherczyk
AGH University of Science and Technology. Cracow, Poland
ABSTRACT: The history and the main steps of development of geomechanics is considered. The
modern state of modelling of the mechanical processes taking place in rock mass is given.
KEYWORDS: Geomechanics, mathematical modeling, mechanical processes, longwall method,
rope saw mining, gas-dynamics manifestation
The history of development of geomechanics covers a relatively short period of time from the be-
ginning of the last century till our days. From 1904 M.M. Protodiakonov publishes a series of arti-
cles dedicated to the problems of mining stress, which later become the base of his monograph
Mining rock stress to mining support published in 1907. That was the first great theoretical work
in the field of rock mechanics. It considered the hypothesis of natural balanced arch over mining
works and described the calculation method of supports for shallow mines. So the development of
geomechanics in Russia began.
W. Budryk's works in the field of mathematical modelling of mechanical processes taking place
in the rock mass during the extraction of layers of coal with the longwall method became the be-
ginning of geomechanics development in Poland. The ideas of W. Budryk were later advanced in
works of his disciples: A. Salustowicz, S. Knothe, J. Litwiniszyn, etc.
In 1911 T. Karmman's works on research of mining rocks in conditions of three-dimensional
compression appeared. He also checked the correctness of the Moor's theory.
In the years 19201950 O.F. Graf, F.A. Belaenko, A.A. Borisov, W. Budryk, F... Van-
Iterson, R. Kvapil, M.I. Koyfman, G.N. Kuznecov, V.D. Slesarev, A. Salustowicz, D.S. Rostovcev,
A.I. Celigorov, P.M. Cimbarevich and other authors published their works, mostly dedicated to the
problems of mining stress and its influence on support, investigation of physical and mechanical
rock properties, rockmass stability. Research projects of this period are mostly based on simple
models of mediums and explored objects.
Later, in the years 19501980, fundamental works in the field of mining stress based on limi-
ting condition theories (V.V. Sokolovsky, N.N. Maslov, G.L. Fisenko, V.N. Zemisev, etc.), on
theories of elasticity, plasticity and creep (F.A. Belaenko, V.V. Vinogradov, J.S. Erjanov, L.V. Er-
shov, B.A. Kartozia, A. Labass, U.V. Liberman, G.G. Litvinsky, L.Y. Parchevsky, A.G. Protosen-
ya, K.V. Ruppeneyt, M.M. Rozovsky, R. Fenner, A.N. Shashenko, E.I. Shemyakin and others)
were published. It was in this period when papers, in which the natural rock mass was considered
on the basis of statistic models (I.V. Baklashov, M.A. Dolgih, V.V. Matvienko, L.Y. Parchevsky,
K.V. Ruppeneyt, A.N. Shashenko, V.I. Sheynin and others) appeared.
Copyright 2006 Taylor & Francis Group plc, London, UK
36
In the same time deep investigations in the field of rock destruction, including controlled load-
ing mode, were conducted in such well-known geomechanic scientific research centres as VNIMI,
Mining Institute named after Skochinsky, IGTM NAS of Ukraine, Dnipropetrovsk Mining Institu-
te, Doneck Polytechnic Institute, Leningrad Mining Institute, Moscow Mining Institute and others
leading scientific research centres in SAR, Japan, England and Poland.
Reviews of these researches were done by V.V. Vinogradov, Z. Bieniawski, A. Cook, A.N. Sta-
vrogin, L.I. Baron, E.I. Ilnitskaja, J.M. Kartashov, Z. Kleczek, G.N. Kuznetsov, B.V. Matveev,
T. Majcherczyk, A.G. Protasenja, V.V. Rzhevsky, V.S. Jamshchikov, G.T. Kirinichansky and other
authors. They allowed creating new, qualitative rockmass models, putting forward and solving
elastic-plastic tasks in the new statement.
By that time an independent line of rock mechanics was formed. It was connected to the re-
search into the gas and dynamic manifestations of the mining stress: sudden burst of coal, rock,
gas, seismicity (S.G. Avershin, B. Drzezla, Z. Kleczek, J. Dubinski, M. Chudek, L.N. Bikov, A.N.
Zorin, I.M. Petuhov, V.V. Hodot, E.I. Shemyakin and others).
Thorough studies of surface movement under the influence of mining were carried out, which
are represented by the works of, for instance, S.G. Avershin, V.F. Galahov, V.N. Zemisev, M.A.
Iofis, S. Knothe, J. Litwiniszyn, R.A. Muller, V.I. Mjakenkij and other authors.
In that period, such frequent in mining workings a phenomenon as floor heave was explored
(P.M. Cimbarevich, V.D. Slesarev, A.P. Maximov, Y.Z. Zaslavskiy, I.L. Chernyak, V.A. Litkin,
V.T. Glushko, A.N. Shashenko and others).
A large piece of work on studying mining stress appearing in natural conditions was carried out
by scientific collectives under the supervision of V.T. Glushko, Y.Z. Zaslavsky, A.N. Zorin, A.G.
Krupennikov, K.V. Koshelev, A.P. Maximov, T. Majcherczyk, I.D. Nasonov, G.I. Pokrovskij, I.L.
Chernyak and others.
In the same time, wide laboratory research on modelling geomechanical processes (G.N. Kuz-
necov, G.I. Pokrovskij, V.F. Trumbachev and others) was undertaken.
In the 70s and 80s, determinative models in geomechanics were explored to a satisfactory de-
gree. As a rule, they were carried out in order to analyse the stability of a single long gallery, which
was the main underground construction structure unit. The determinative model required idealiza-
tion of the research object (mining working, rock mass), and boundary conditions. This circumsta-
nce allowed using solid mechanics techniques and obtaining distinct solutions of the problems. But
comparing actual measurements with analytical calculations showed that high idealization of the
object in geomechanics led to such a situation: a range of results in actual conditions fluctuated
around the forecasted result, corresponding to a single distinct result of a theoretical calculation.
The range of factors that were not included into the physical model due to their unimportance (ac-
cording to the initial viewpoint) was the main reason for these differences.
This circumstance caused a distrust among the practising mine designers and engineers, com-
pelling them to correct analytical calculations and to accept technological and technical parameters
intuitively, first of all basing them on the experience gained in similar underground mining condi-
tions. This approach very rarely results in making optimal engineering decisions, which leads to
the increase of expenditure on maintaining the mining workings or superfluous support stability
and higher costs of support. It is especially seen in the mines, which were built in difficult mining-
geological conditions in the Eastern and Western Donbass, Krasnoarmeysk coal region, Lviv-Vo-
lyn coal region etc.
Obviously it is the unconformity of determinative models with real objects that is the reason for
such a situation. It is impossible to forecast behaviour of such difficult objects like rock mass
containing a mining excavation using only models where everything is defined and there is no pla-
ce for variability. Any highly organized system is formed and exists effectively when it has some
internal instability. Then the strict connection of one structural unit, existing by action of some
Copyright 2006 Taylor & Francis Group plc, London, UK
37
common physical laws, would not make any strong restrictions to other elements, which should
have an opportunity to be changed in a random way.
This shows that geomechanical research should be based on stochastic models, which are more
complicated and reflect fluctuation of parameters around determinative values. Understanding this
problem led to the occurrence of some research projects based on probability models during the
last 1520 years (K.V. Rupeneyt, L.Y. Parchevsky, A.N. Shashenko, S.B. Tulub, V.P. Pustovoy-
tenko, G.T. Rubec, E.A. Sdvijkova and others). The spectrum of these works is considerable:
strength-strain state of rock, floor heave, expense optimisation for long-life underground construc-
tions, deformation of support, quick analysis of management decisions etc.
Now we can surely say that statistical geomechanics became an independent branch of science
dealing with the mechanical processes taking place in rocks.
At the end of the 80s geomechanical models became so difficult that using analytical methods
(as they are) ceased to be effective. At that time personal computers became more powerful and ac-
cessible for carrying out research projects. With the development of personal computers numerical
methods became very effective and popular, as it may be judged from the research projects by B.Z.
Amusin, J.S. Erjanov, T.D. Karimbaev, S. Krauch, L.V. Novikova, E.A. Sdvijkova, A. Starfield,
B.A. Fadeev and others.
This direction in the development of geomechanics is considerable and perspective.
The combination of continuum and discrete medium methods with new geoinformation techno-
logies, geophysics, theory of mathematical statistics constitutes the basis for creating a system of
global geomechanical computer monitoring. Then the exploration of underground space can be-
come more effective and safer.
Copyright 2006 Taylor & Francis Group plc, London, UK
39
Speed of Roof Rock Separation and a Type of Workings Support
Tadeusz Majcherczyk, Piotr Makowski, Zbigniew Niedbalski
AGH University of Science and Technology, Cracow. Poland
ABSTRACT: The paper presents measurement results of speed of roof rock separation for three
selected headings protected with various types of support, i.e. individual roof bolting, combined
system of prop support and roof bolting, and prop support. The measurements were carried out for
a period of dozen to several dozens of months. Separations were recorded with the use of low and
high multilevel telltales as well as extensometers. Speed of separations was determined in relation
to the time of measurement. Recording separations for a long period of time allowed for an esti-
mation of stability of the analysed workings as well as for relating movement of roof layers to va-
rious types of support applied in their construction.
1. INTRODUCTION
Speed of roof rock separation is one of the basic factors affecting stability of headings and it is
most often understood as a quotient of separation value and an assumed unit of time [Butra &
Orzepowski 2001]. The measurement of separation speed can be carried out with the use of special
measuring instruments or standard telltales and extensometers if the changes around the heading
are going on for a long period of time [Butra 1998, Fabianczyk 2001, Majcherczyk et al. 2005].
A considerable intensiveness of the phenomenon requires the use of a proper method of limiting
roof rock separation, thanks to which it is possible to preserve the workings stability [Gale &
Blackwood 1987]. Apart from geological conditions and mining situation, a type of support applied
in the working is one of the crucial factors affecting changes of separation in time [Majcherczyk &
Niedbalski 2002a, 2002b, Szymiczek et al. 2001, Fabich et al. 2000]. Separation speed will not
influence the stability of heading with a prop support or a combined system of prop support and
roof bolting, nevertheless in the case of heading protected with an individual roof bolting, high
speed of separation can cause hazard to the safety of miners.
The paper analyses changes of separation speed in roof rocks, where various types of headings
support were applied. In the analysed workings, i.e. in the heading 2-E1 in the seam 703/1, in the
inclined drift Izn in the seam 358/1 and in the inclined drift B-1 in the seam 404/1, the following
types of support were applied respectively: the individual roof bolting, the combined system of
prop support and roof bolting, and the traditional P-type arch support. Although the workings
were located in three different hard coal mines, their location was selected in such a way that the
roof conditions were to be comparable in all cases. Due to the fact that it was expected that the
rock mass movements would be happening in a long period of time, the workings were drilled in
the areas where mining exploitation had not been carried out before.
Copyright 2006 Taylor & Francis Group plc, London, UK
40
2. CHARACTERISTICS OF MINING-GEOLOGICAL CONDITIONS
AND SUPPORT SCHEMES
2.1. Rise gallery I-E1 in the seam 703/1 individual roof bolting
The area of the seam 703/1, where the rise gallery I-E1 was drilled, is located at the depth of appro-
ximately 970 m. The thickness of the seam is 1.92.2 m with the inclination of up to 5. The roof
and the floor of the working consist of shale clays and arenaceous shales. The rise gallery I-E1 is in
the neighbourhood of the walls caving goaf. It may be also assumed that the working is not sub-
jectted to the interference of exploitation edges as their vertical distance considerably exceeds 200 m.
The rise gallery I-E1 in the seam 703/1 has a rectangular cross section with the dimensions:
width 5.0 m and height 2.5 m. In the headings roof, seven steel bolts with the total length of 2.35 m
were fixed in a row (Fig. 2.1). The distance between the bolts in the row is 0.8 m and the distance
between the rows is 0.9 m. The wall of the heading from the side of unmined coal is protected with
three rows of steel bolts with the total length of 2.35 m, but from the side of the front, the wall was
protected with two rows of workable bolts with the total length of 2.0 m. The length of the section
of the heading I-E1, in which individual roof bolting was applied, is 182 m.
5,0 m
2
,
0
m
0,8 m
0,8 m 0,8 m 0,8 m
0
,
8
m
steel bolts
lenght 2,35 m
1
,
4
m
0,8 m 0,8 m
workable bolts
lenght 2,20 m
Figure 2.1. Scheme of bolt distribution in the rise gallery I-E1 in the seam 703/1
2.2. Inclined drift Izn 358/1 combined system of prop support and roof bolting
The seam 358/1 in the area of the inclined drift Izn is located at the depth of approximately 900 m.
The thickness of the seam is approximately 2.5 m with the inclination of up to 5. The 1.5-meter
thick layer of shale clay constitutes the immediate roof, above which there are sandstones and are-
naceous shale. The neighbouring seams were not exploited and the inclined drift Izn was drilled
parallel to the caving goaf located in the distance of 5070 m.
The inclined drift Izn in the seam 358/1 was constructed in the P9/V29-type support, so the
width of the working is 5.0 m with the height of 3.5 m (Fig. 2.2). The spacing of the sets in the
analysed working is 1.5 m. The support is reinforced with two pairs of bolts fixed to the roof-bar
arch at the distance of approximately 1.0 m from the workings axis. The applied bolts have a total
length of 2.5 m and were fixed on all the length. The section of the inclined drift Izn subjected to
measurements had a length of approximately 100 m.
Copyright 2006 Taylor & Francis Group plc, London, UK
41
steel bolts
lenght 2,50 m
5,0 m
P9/V29
Figure 2.2. Scheme of bolt distribution in the inclined drift Izn in the seam 358/1
2.3. The inclined drift B-1in the seam 404/1 prop support
The seam 404/1 in the area of the inclined drift B-1 has the thickness of 2.052.15 m and the
inclination of approximately 5. In the seam there is an interlayer of shale with coal laminas with
the thickness of 0.350.40 m. There are no exploitation edges in the premises of the working. How-
ever, in its neighbourhood there used to be active wall fields in the seam 403/1 and 404/2. The
distance between the seam 404/1 and the seam 403/1 in the area of the inclined drift B-1 is between
approximately 10 m and 20 m.
In the immediate roof of the working there is shale clay with occasional interlayers of arenaceous
shale and sandstone. In the floor, to approximately 25 m, shale clay with occasional interlayers of
arenaceous shale can be found. Immediately below the working, carbonaceous shale or shale clay
with coal laminas may occur.
No tectonic disturbance occurs in the runway of the inclined drift B-1.
The working was built in the P-V29/9-type support, so it has a length of 5 m and the height of
3.5 m, with the spacing of 1 m.
3. ANALYSIS OF SEPARATION SPEED
The measurements of separation were carried out using telltales of high and low separation, multi-
level telltales and extensometric probes.
In the rise gallery of the wall I-E1, 11 stations for measuring low and high separations as well as
multilevel telltales were fixed. The distance between measurement stations is approximately 20 m.
The high telltales were fixed at the height of 4.5 m and the low telltales 2.3 m. The critical value
of separation in the rise gallery I-E1 was determined as follows: for low separation 25 mm, for
high separation 40 mm. The frequency of measurements was between 2 and 4 times per week du-
ring the drilling of the rise gallery, and between 5 and 8 times per month after terminating the drill-
ling process.
Five stations for measuring low and high separation were fixed in the inclined drift Izn. High
telltales were fixed at the height of 5.3 m, and low telltales 2.5 m. The distance between the tell-
tales was approximately 20 m. Due to the application of a combined system of prop support and
Copyright 2006 Taylor & Francis Group plc, London, UK
42
roof bolting, the critical value of separation was not assumed. The frequency of measurements was
once per month in the first year and later it was once per three months.
In the inclined drift B-1 the measurements of separation were carried out with the use of exten-
sometric probes. They were fixed in the roof of the working in the two measurement stations, ap-
proximately 67 m distant from one another. Roof separations were controlled using the above-
mentioned instruments on 18 levels of the immediate and the proper roof to the depth of approxi-
mately 7.5 m. The precision of the record, based on the electromagnetic signal, was 0.025 mm. The
measurements were initially carried out every week and later once per month.
3.1. Rise gallery of the wall I-E1
The measurements of the rise gallery I-E1 were carried out for a period of 22 months. For the sake
of the analysis of separation speed, the low and high telltales No 7 and No 8 as well as the
multilevel telltales No 3 were used. Fig. 3.1 presents the separation speed recorded with the use the
telltales No 7. The RN symbol refers to low separation, whereas the RW symbol refers to high
separation.
0
0,05
0,1
0,15
0,2
0,25
0,3
m
m
/
d
a
y
50 100 150 200 250 300 350 400 450 500 550 600
day
RN7
RW7
Figure 3.1. Speed of separation in time in the telltales No 7 in the rise gallery I-E1
It may be assumed from the graphs presented that the highest separation speed occurs in the
period of approximately 100 days and it reaches nearly 0.3 mm/day. It should be pointed out here
that in the period of the initial 50 days of measurement, separation speed was similar for low and
high separation. In the consecutive measurements it was observed that separation speed was de-
creasing systematically, reaching the value of 0.05 mm/day in the final period. The decrease of the
analysed parameter could be also observed in relation to high telltales, although in the period be-
tween 250th and 300th day of measurement separation speed increased from 0.112 mm/day to 0.16
mm/day. However, this phenomenon occurred only in the high telltales and after this a systematic
decrease of separation speed could be observed, and its final value was 0.08 mm/day.
A slightly different character of changes of separation speed in particular periods of time could
be observed with the use of the telltales No 8 (Fig. 3.2).
Copyright 2006 Taylor & Francis Group plc, London, UK
43
0
0,1
0,2
0,3
0,4
0,5
0,6
m
m
/
d
a
y
50 100 150 200 250 300 350 400 450 500 550 600
day
RN8
RW8
Figure 3.2. Speed of separation in time in the telltales No 8 in the rise gallery I-E1
On the basis of the above graph it may be assumed that separation speed was decreasing sys-
tematically, initially having its peak value of 0.44 mm/day for the low telltales and 0.60 mm/day
for the high telltales. Increased value of the analysed parameter can be also observed for high se-
paration, which suggests more considerable increase of separation. The final values of separation
speed were slightly lower than the ones recorded in the station No 7 and were 0.036 mm/day and
0.053 mm/day respectively.
0
0,1
0,2
0,3
0,4
0,5
0,6
m
m
/
d
a
y
50 100 150 200 250 300 350 400 450 500 550 600
day
0-1m
0-2m
0-3m
0-4m
Figure 3.3. Speed of separation in time in the multilevel telltales No 3 in the rise gallery I-E1
Copyright 2006 Taylor & Francis Group plc, London, UK
44
The graph representing separation speed for the multilevel telltales seems to offer some interes-
ting possibilities for interpretation (Fig. 3.3). After studying it carefully it may be assumed that the
dynamics of separation on particular levels increases significantly from 0.1 mm/day for the level
01 m to 0.44 mm/day for the level 04 m. Even more significant differences, more than quintuple,
occur in the period between the 50th and 100th day of measurement. Then the separation speed
reaches 0.6 mm/day for the separation in the packet 04 m.
After this period of time the dynamics of changes decreases systematically for all levels, and
the differences in values for particular levels, which had occurred in the initial phase of measure-
ments, remained unchanged till the very end of the measurements.
In the case of the multilevel telltales it may be assumed that separation speed depends on an
analysed rock thickness. Usually, the thicker the analysed packet of roof rocks, the bigger the total
separations. Similarly, using a similar way of monitoring, separation speed of particular sections of
the workings roof, which is a difference of separations speed on the analysed lengths, can be de-
termined. For instance, for the rise gallery I-E1 in the initial 50 days of measurement, the separa-
tion speed between the second and the third metre of the roof is 0.1 mm/day.
To sum up, it may be claimed that the dynamics of roof rock separation in the rise gallery I-E1
in particular telltales has a very similar course. The biggest values of separation were recorded in
the initial period of measurements and, with time, the values were systematically decreasing. As far
as the maximal absolute values are concerned, separation speed changed from 0.3 mm/day to 0.6
mm/day depending on the location of the measuring instrument.
3.2. Inclined drift Izn
The results from the inclined drift Izn were obtained during a 37-month period of measurements
and the low and high telltales No 6 and No 9 were used in order to analyse the separation speed.
Figure 3.4 presents the separation speed recorded by the telltales No 6. The symbol RN refers to
low separation and RW to high separation.
0,00
0,05
0,10
0,15
0,20
0,25
0,30
0,35
0,40
0,45
0,50
m
m
/
d
a
y
0 42 99 155 311 371 417 617 837 1102
day
RN6
RW6
Figure 3.4. Speed of separation in time in the telltales No 6 in the inclined drift Izn
It may be assumed from the presented graphs that the highest sped of separation occurs in the
period of the initial 10 days and it is 0.5 mm/day. Until the 155th measurement day, a significantly
Copyright 2006 Taylor & Francis Group plc, London, UK
45
high value of separation speed occurs and it is 0.2 mm/day. In the later period of time the dynamics
of separation speed decreases systematically to the point in which the values of low and high sepa-
rations are similar. It reflects the changes going on in the packet of the immediate roof at the depth
of approximately 2.5 m.
Significantly higher values of separation speed were observed in the case of the telltales No 9
(Fig. 3.5), as in this case separation speed in the period of the initial 10 days was 1.2 mm/day, both
for low and high separation. Almost two times lower value of the analysed parameter was obtained
in the following 30 measurement days. Also in the case of the telltales No 9, a regular decrease of
separation speed occurs and the final values were 0.044 mm/day for high separation and 0.043
mm/day for low separation.
0,0
0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1,0
1,1
1,2
m
m
/
d
a
y
0 42 99 155 311 371 417 617 837 1102
day
RN9
RW9
Figure 3.5. Speed of separation in time in the telltales No 9 in the inclined drift Izn
On the basis of the comparison of the obtained results from the ones obtained in the rise gallery
I-E1 built in an individual roof bolting, it may be concluded that in the case of the working with the
combined system of prop support and roof bolting, separation occurs much faster. The time, in
which extreme values are reached, is another factor differentiating the changes of separation speed;
for the working with a combined system of prop support and roof bolting they occur in a shorter
period of time. However, some similarities can be observed in the character of changes of the ana-
lysed parameter as separation speed decreases with time regardless of the type of support system.
3.3. Inclined drift B-1
The measurements of separation in the inclined drift B-1 were carried out for the period of 28
months. This paper presents the results recorded with the use of the probe No 2. Out of the set of
18 results, only the results recorded in 8 selected sensors fixed at the heights of 0.59 m, 1.26 m,
2.30 m, 3.43 m, 4.79 m, 5.82 m, 6.92 m and 7.52 m were used for the sake of the presentation in
this paper (Fig. 3.6).
On the basis of these results it may be assumed that separation speed, similarly to the above-
mentioned cases, decreases with time, but its character of changes is different. In the period of the
initial 7 days, separation speed has a peak value of between 0.26 mm/day and 0.68 mm/day. In the
following periods of time separation speed decreases. The graph representing separation speed has
Copyright 2006 Taylor & Francis Group plc, London, UK
46
a form of a wavy curve, as its values increase and decrease alternately, however with the tendency
towards zero. After 171 days of measurement, separation speed of roof rocks does not exceed 0.07
mm/day, and after 444 days 0.003 mm/day. Despite practically stopping the process of separation
until the 594th measurement day, after the 600th day the increase of fissure opening occurs again.
As the strata control in the working was maintained by means of a prop support, whose work cha-
racteristics is based on cyclic slides alongside with the progressing pressure of the rock mass on the
support, the oscillatory movement of roof layers can be caused by the work of the support.
It should be also pointed out that separation speed of roof rocks might be negative, which means
the compaction of rock layers and closing fissures. Such a situation may occur not only in the first
metres of the roof (sensor at the depth of 0.59 m) but also in the further distance from the work-
ings contour (sensor at the depth of 4.79 m). Due to a very high precision of the measurement and
a high density of distribution of measurement sensors, such a result could be obtained only in the
case of controlling rock mass movement with the help of the extensometric probe.
-0,050
0,000
0,050
0,100
0,150
0,200
0,250
0,300
0,350
0,400
0,450
0,500
0,550
0,600
0,650
0,700
m
m
/
d
a
y
7 21 39 73 112 140 171 200 231 266 322 357 382 417 444 475 504 531 566 594 630 656 685 747
day
0,588 m
1,266 m
2,298 m
3,437 m
4,790 m
5,821 m
6,920 m
7,520 m
Figure 3.6. Roof separation speed in the inclined drift B-1
4. CONCLUSIONS
On the basis of the above considerations and evaluations of separation speed in underground head-
ings in relation to a type of support system applied, the following conclusions may be drawn:
1. Separation speed always asymptotically decreases to zero, having a value of between 0.3 mm/day
and 1.2 mm/day in the initial period of time.
2. The process of disappearance of separation in the case of applying an individual roof bolting is
the longest. A rock beam created by fixing the bolts subsides for a much longer period of time
and even after 600 days of measurements the speed of dislocations may still have the value of
0.15 mm/day.
3. Roof separation has a changeable character, i.e. it can increase or decrease in the consequent
periods of time.
Copyright 2006 Taylor & Francis Group plc, London, UK
47
4. Separation speed of roof rocks can be negative, which suggests compaction of rock layers and
closing fissures.
5. Only extensometric probes allow for a precise control of roof strata above a working. In this
case, the precision of measurement is 0.001 inch. Controlling roof layers with the use of tell-
tales does not provide such possibilities.
REFERENCES
Butra J. & Orzepowski S. 2001: Ocena statecznoci stropu wyrobiska na podstawie prdkoci jego rozwars-
twienia. Miesicznik WUG, nr 5, pp. 2230.
Butra J. 1998: Kierunki rozwoju metod monitoringu stropu w warunkach wzmoonych cinie i przemiesz-
cze grotworu. Cuprum, nr 3, pp. 7796.
Fabich S., Pytel W. & Skopiak L. 2000: Badania zachowania si obudowy kotwowej stropu wyrobisk na pod-
stawie pomiarw wielopunktowymi sondami ekstensometrycznymi i kotwami oprzyrzdowanymi. Cuprum,
nr 17, pp. 7394.
Fabjanczyk M. 2001: The Use of Monitoring for the Design and Management of Bolting System. Nowoczes-
ne technologie grnicze 2001 Kotwienie, pp. 135148.
Gale W.J. & Blackwood R.L. 1987: Stress Distributions and Rock Failure around Coal Mine Roadways. Int. Jour-
nal Rock Mech. Min. Scien. & Geomech., Vol. 24, nr 3, pp. 165173.
Majcherczyk T. & Niedbalski Z. 2002a: Badania zasigu strefy spka wok wyrobiska korytarzowego. Prze-
gld Grniczy, nr 2, pp. 914.
Majcherczyk T. & Niedbalski Z. 2002b: Ocena obudowy podporowo-kotwiowej na podstawie wybranych ba-
da in situ. Przegld Grniczy, nr 12, pp. 17.
Majcherczyk T., Makowski P. & Niedbalski Z. 2005: Separation of Roof Rock Observed in Headings under
Development. Proceedings of the 6th Int. Mining Forum 2005, 2327 February, Cracow, Szczyrk, Wie-
liczka. Rotterdam, Balkema, pp. 5362.
Szymiczek W., Ficek J. & Sobik M. 2001: Nowoczesne metody monitoringu pracy obudowy kotwiowej. No-
woczesne technologie grnicze 2001 Kotwienie. Gliwice, Ustro, pp. 453468.
Copyright 2006 Taylor & Francis Group plc, London, UK
49
The Use of Pre-stressed Tendons
in the Polish Coal Mining Industry
Zbigniew Rak, Jerzy Stasica, Micha Stopyra
AGH University of Science and Technology, Cracow. Poland
ABSTRACT: The paper outlines the practice of roof support with the use of long string wire cable-
bolts in Polish coal mines from its beginning till now. Shown are the results of the observations
conducted to determine the interaction between the cables and the rockmass in various geological
and mining conditions. The documented cases of cablebolt support installation and performance
monitoring have been done on such mines as: Bogdanka, Sobieski Jaworzno III, Murcki,
Wieczorek, Bielszowice, Piast. The authors also present the results of laboratory research
conducted at a suitably prepared test stand as well as new cable construction designs together with
their technical characteristics. Further on described are the most effective methods to install cable-
bolts in various rocks, characterised by diverse structural and strength properties.
KEYWORDS: Pre-stressed tendons, coal mining, mine tunnels, support reinforcement
1. INTRODUCTION
The beginnings of the use of rockbolts for rock reinforcement dates back to the end of the 19th
century. On a wider scale, however, this support method was employed only half a century later.
As from the 1950s rockbolting became the basic means of support used in Europe in underground
mining and civil engineering excavations. The current world usage or rock anchors is estimated at
some tens of millions of units installed per year.
Rope tendons, an important type of rock anchors, have been in use in the mining industry for
more than thirty years. They were first installed during the field trials conducted in underground
metal ore mines at the late 1960s. Grouted rope strands 15.2 mm in diameter were used in the
1970s as rock reinforcements. At that time they did not yet come with auxiliary fittings such as fa-
ceplates and locking nuts. By the end of the 1980s rope bolts were already in common use in un-
derground coal mines of South Africa, Australia and the USA.
Thanks to their construction rope bolts possess several advantages making them particularly
useful in the most difficult mining and geological conditions, particularly conditions of high static
stress and large de-stressed zones around excavations. In high seismic hazard areas they perform
the function of so called barrier support.
Cablebolts have been in wider use in the Polish coal mining industry for the past 10 years only.
Extensive development has taken place during this period of time, both in the design of cable con-
struction and in the technology of longhole bolting itself. Tendon construction, technology and me-
chanization of installation and support performance monitoring development took place in all
these areas. The paper presents the practical experience gained in Polish mines in the use of long-
hole rope and string wire cablebolt support so far. The authors also present the results of laboratory
Copyright 2006 Taylor & Francis Group plc, London, UK
50
research conducted on this support type and present new cable anchor designs together with the
most important technical characteristics of the described units. Further on described are the most
effective methods to install wire cablebolts in rocks characterised by diverse structural and strength
properties.
2. HIGH BOLTING WITH THE USE OF STRING WIRE AND ROPE CABLEBOLTS
It is often necessary in underground mines to install, either as primary or secondary support, rock-
bolts whose lengths exceed (sometimes manifold) the height of the excavation. Traditionally the
Polish nomenclature would call rockbolting high whenever the length of support units exceeded
the height of the supported excavation. However, in the case of large-sized excavations this de-
finition does not convey the true meaning of high bolting. Many authors formulated their own de-
finitions, stipulating different rockbolt lengths to differentiate between high and low bolting. It
is largely accepted in the Polish coal mining industry to refer to rockbolting as high in the cases
when anchor length exceeds 3,0 metres.
It is relatively seldom that stiff steel rods are used for high bolting. The reason is the necessity
to join rod sections (usually with threaded joins), making the installation process troublesome. Lar-
ger diameter holes are normally needed to accommodate them as well. The notch effect of the joint
weakens the whole unit. For these reasons cables rather than rods are utilized. The most popular
definition of a cablebolt describes its construction as comprising several individual or twisted steel
wires (anchors suitable for pre-tensioned installations) or twisted strand ropes [Stilborg 1994]. In
the Polish nomenclature these two types are called respectively string wire and rope cablebolts.
The most important feature of a cablebolt is its flexibility, rendering it possible to install cable
of practically any length regardless the dimensions of the excavation. This makes it possible to re-
inforce a wider zone of rock around an excavation than can be done with low bolting. The im-
portant advantage of cable bolts over plain rods is their notably greater tensile strength for the same
diameter tendons. Subjected to the same loads cables normally display smaller deformations. How-
ever, it is possible to increase yielding capacity of a cable without changes to its construction.
Thanks to the named advantages, cablebolts are used increasingly more frequently for low rock-
bolting as well [Rataj 1995, Rataj 2000].
Figure 2.1 shows types of string wire cablebolts most commonly used in the world mining in-
dustry. Of them, only the types shown as a, d and f are used on a larger scale in Polish coal mines.
In the 1990s several field trials with 18 and 20 mm in diameter rope cablebolts were conducted. An
example of such construction is a unit manufactured and marketed as IR-3 by INTERRAMM. It is
manufactured from standard 18 mm diameter twisted strand ropes type 67+A
0
(18,0-67+A
0
-Z/s-
n-II-g-140) (Photo 2.1a, b and c). The tendon of this rope comprises 6 strands, each of them con-
sisting of 7 twisted steel wires. Constructed in this way rope comprises 42 steel wires and a soft co-
re, usually made from jute fibres. The cross-sectional area of steel wires in a 18 mm diameter rope
is approximately 131 mm
2
. Its breaking force is at least 160 kN. The tendon of this cablebolt is mo-
dified by untwisting the original rope over the length corresponding to the length of bond. The rope
prepared in this way and grouted in a hole forms a kind of disorderly steel reinforcement.
The maximum diameter of the bulbed part of a cable called cage is up to 38 mm. In its thin-
nest part, where the wires remain twisted the diameter of the tendon equals to approximately 24 mm.
The advantage of this construction is that the cable remains securely fastened even if the grout re-
veals certain contraction. Moreover, contrary to smooth cable, after grouting the bulbed strands do
not have the tendency to break loose from the grout under load. The locking mechanism of this bolt
comprises the 3-part Gifford lock.
Figure 2.2 shows an example of a load-deformation graph of IR-3 cable bolt. The units used for
the test were left ungrouted over a length of 1,0 metre. The measured load-bearing capacity was not
Copyright 2006 Taylor & Francis Group plc, London, UK
51
less than 175 kN and the total elongation of the tendon at the moment of failure was approx. 5%. Ca-
blebolts of this type were frequently used as reinforcement in face gates supported solely by steel
dowels at, among others, KWK Ziemowit, KWK Wesoa and KWK Czeczott mines. They
were also used at road crossings at KWK Piast and KWK Wesoa and to reinforce footwall at
KWK Wesoa mine. Despite promising laboratory test results their use have been discontinued
due to inability to overcome technological problems, such as difficulties with correct installation of
the cablebolts locking mechanism (Gifford lock) in practical underground working conditions and
problems with forcing the cables into boreholes pre-filled with thick cement grout. It was also in
the 1990s when the use of string wire cablebolts designed for pre-stressed installations began in
Poland on a larger scale. They were 15,5 mm in diameter, which made them easier to insert into
boreholes. Cablebolts of this type are manufactured by ARNALL POLAND and are known under
market name KL-15,5 or KL-15,5K (Photo 2.2a, b and c).
Figure 2.1. Most common cable bolt designs [Windsor 1992]
a b
c
Photo 2.1. IR-3 rope cablebolt: a lock, b cable, c untwisted length of tendon (birdcage)
Copyright 2006 Taylor & Francis Group plc, London, UK
52
0
20
40
60
80
100
120
140
160
180
200
0 10 20 30 40 50 60 70 80 90 100 110 120
yield [mm]
l
o
a
d
[
k
N
]
Figure 2.2. Working characteristics of IR-4 cablebolt
The cable of this bolt, manufactured according to PN-71/M-80236 standard, comprises seven
axially twisted steel wires. The total cross-sectional area of the wires is 143 mm
2
. The breaking
force of such tendon is approximately 250 kN. The cable may be either smooth (KL-15,5) or geo-
metrically modified (KL-15,5K) by the formation of small cages, about 100 mm long and 26 to 28 mm
in diameter. The cages are usually spaced every 200 mm. Cablebolts with greater birdcage spacings
are also manufactured. The locking mechanism designed for this bolt is also the 3-part Gifford lock.
a)
b)
c)
Photo 2.2. KL-15,5 (KL-15,5K) wire cablebolt: a cable, b birdcage, c lock
0
50
100
150
200
250
0 10 20 30 40 50 60 70 80 90
yield [mm]
l
o
a
d
[
k
N
]
Figure 2.3. Working characteristics of KL 15,5 cablebolt
Copyright 2006 Taylor & Francis Group plc, London, UK
53
The stress-strain characteristics of KL 15,5 cablebolt is shown on picture 2.3. The total elonga-
tion of the bolt at the moment of failure was approx. 4%.
Another wire cablebolt commonly used in Polish coal mines is IR-4, manufactured by INTER-
RAMM BIS (Photo 2.3a, b and c). The cable of this bolt comprises eight 6 mm diameter bound in
parallel wires.
a) b) c)
Photo 2.3. IR-4 wire cablebolt: a cable, b birdcage, c lock
The total cross-sectional area of the wires is approximately 226.2 mm
2
. The breaking force of
such cable is approximately 320 kN, with the manufacturer-guaranteed load-bearing capacity of
280 kN. One meter of the cable weighs 1.75 kg. The wires comprising the cable are made of BII
spring wire as per PN-71/M-80057 standard. The cable of the IR-4A cablebolt may be either smooth
or modified to form birdcages (Fig. 2.3b). The diameter of the smooth cable is 2224 mm. In the
case of the modified one, its maximum diameter is 38 mm. The lock of IR-4 cablebolt comprises
a threaded sleeve with a nut. The cable is secured in the sleeve by means of a steel wedge. At load,
individual wires of the cable move and gradually wedge themselves tight in the sleeve until their
further movement is stopped. Such design provides for proper wedging of the wires in the lock and
allows eliminating their movement relative to each other already at installation, when initial load is
applied. Wedging of the wires at an early stage of pre-tensioning is important for the cable to work
properly. The strain-stress characteristics of IR-4 wire cablebolt is shown on Figure 2.4. The total
movement of the cable at failure, which is the sum of wire movement on the lock and the cable elon-
gation, does not exceed 5%. Based on the IR-4 design new types of cablebolts similar in construc-
tion were created IR-4C and IR-4W. IR-4C cablebolt was designed jointly by the scientists from
the Underground Mining Department of AGH (the Technical University of Mining and Metallurgy,
Krakw) and INTERRAMM BIS. The cable of this bolt comprises nine twisted 6 mm diameter
each wires. The total cross-sectional area of the wires is approximately 254.34 mm
2
. The breaking
force of the cable is not less than 350 kN, and the elongation at failure does not exceed 5% (Figure
2.5). The diameter of IR-4C cable is between 22 and 24 mm.
The cable of IR-4W bolt comprises seven 8 mm diameter wires. The cablebolt is about 26 mm
in diameter.
Copyright 2006 Taylor & Francis Group plc, London, UK
54
0
50
100
150
200
250
300
350
0 20 40 60 80 100
przemieszczenie [mm]
o
b
c
i
e
n
i
e
[
k
N
]
Fig. 2.4. The strain-stress characteristics of IR-4 wire cablebolt
0
50
100
150
200
250
300
350
400
0 20 40 60 80 100
yield [mm]
l
o
a
d
[
k
N
]
0
50
100
150
200
250
300
350
400
450
500
0 10 20 30 40 50 60 70 80 90
yield [mm]
l
o
a
d
[
k
N
]
Figure 2.5 (left). Working characteristics of IR-4C cablebolt
Figure 2.6 (right). Working characteristics of IR-4W cablebolt
The total cross-sectional area of the wires is approximately 351.68 mm
2
. The outside diameter
of the sleeve forming the lock has in this case been increased to 39 mm. It is made from the same
material, i.e. 40 HM steel (Cr-Mo steel, Polish classification). The maximum load at failure for this
bolt is approximately 440 kN, with the elongation not exceeding 5% Figure 2.6. Because of its
relatively good strength characteristics (highest load-bearing capacity of all the cablebolts used in
the Polish coal mining industry), IR-4W cablebolts were subjected to dynamic loading tests. The
tests were conducted in the specialised laboratory of the Central Mining Institute in Katowice. Bolts
were subjected to dynamic loads created by the impact of a falling weight. Kinetic energy of each
impact was 25 kJ. Some of the results of the research, presented as load-time graphs, are shown on
Figures 2.7 and 2.8.
yield [mm]
l
o
a
d
[
k
N
]
Copyright 2006 Taylor & Francis Group plc, London, UK
55
0
100
200
300
400
500
600
1,8 1,9 2 2,1 2,2
t ime [s]
d
y
n
a
m
i
c
l
o
a
d
[
k
N
]
0
100
200
300
400
500
600
1,8 1,9 2 2,1 2,2
t ime [s]
d
y
n
a
m
i
c
l
o
a
d
[
k
N
]
Figure 2.7 (left). Load-time graph for IR-4C cablebolt obtained by impact-loading it
with 25 kJ of dynamic energy
Figure 2.8 (right). Load-time graph for IR-4W cablebolt obtained by impact-loading it
with 25 kJ of dynamic energy
The obtained results allowed classifying IR-4W string wire cablebolt as suitable for installation
in dynamic load conditions.
3. WIRE CABLEBOLT INSTALLATION METHODOLOGIES
The choice of installation technique is determined by several factors, such as type of grout, the
length and diameter of the support hole, the grouted length of the hole and the available equipment.
In general, three cablebolt installation methods can be named:
1. With the use of thin cement grouts and liquid resins.
2. With the use of thick cement grouts.
3. With the use of cement or resin capsules.
Until recently the first method was the only one utilized to install long cablebolts in roof and
sidewalls. The principle of the method is schematically explained on Figure 3.1. Used in this case
particular are thin cementitious grouts or two-component resins. After the hole has been drilled
a cablebolt equipped with two tubes is inserted a breather tube, with length equal to that of the
cable, and a shorter tube for injecting the grout. The mouth of the hole is then sealed with a 30 cm
long polyurethane foam plug. To the prepared in this way hole liquid grout is injected. The hole is
filled starting from the bottom, with the air escaping from the hole via the breather tube. Pumping
of the grout is stopped when it starts to come out the breather tube. The major disadvantage of this
method is its high labour intensity. Additional problem in the cases where cement grouts are used
is their long setting time and the possibility of the grout to penetrate from the hole into fractures.
Additional disadvantage are the relatively poor strength parameters of the cement stone formed
from the setting grout. The use of two-component resins considerably increases the cost of grouting
and does not eliminate the risk of the grout escaping into the cracks. In the second method thick-
ened cement grouts are used (Figure 3.2). In this case filling of the hole starts at its further end. The
bolt is inserted by hand (or mechanically) into the hole, which is completely or partially filled with
grout. The method allows achieving higher bond strengths, mainly thanks to the relatively high
Copyright 2006 Taylor & Francis Group plc, London, UK
56
strength parameters of grout (even above 50 MPa). With the lower water-cement ratio the setting
time is also shortened (below 1 h). Troublesome with this method is the necessity to prepare fresh
grout directly before application and the labour-consuming process of thorough cleaning of the
grouting equipment.
Breathing
tube
Anchor head
Injection tube
Tendon
Plug
hangingwall
1 2 3
Figure 3.1. Cablebolt installation method with the use
of thin cementitious grouts and liquid two-component resins
1 2 3
Injection tube Rope tendon
Figure 3.2. Cablebolt installation
with the use of thick cementitious grouts
At the moment, the most often used in Polish mines for cablebolt installation is the method with
cement or resin capsules. Thanks to its simplicity and susceptibility to mechanisation the method
allows to substantially shortening the time of the cablebolt installation process. This in turn affords
installation directly at the face of the driven heading or, in the worst case, just behind the road-
header. Bolting of rock surface as soon as possible after it is exposed is beneficial to the stability of
the excavation. The use of fast-setting resin capsules allows to pre-tension the bolt directly after the
installation, resulting in its almost instantaneous interaction with the rockmass. The disadvantage
of the method is impossibility of grouting the cable over length greater than 1,5 metre. Usually two
60 cm long 3-minute cartridges per bolt are used due to the flexibility of the cable and the resistan-
ce encountered while mixing the capsules. It must be noted that both laboratory and in situ tests
show that cable failure loads can be achieved for bond lengths as small as 50 cm.
One of the most important parameters determining the proper installation of cablebolts with re-
sin capsules is the torque available on the installation equipment. The most common type of equip-
ment used with this method is a pneumatic roofbolter (e.g. Gopher). However if pressure of the
compressed air in the supply line is too low, it does not guarantee correct bolt installation. Relati-
vely recently introduced hydraulic roofbolters VSP-1 powered by their own hydraulic aggregate
were introduced. High efficiency of this equipment resulted in substantial increase in the quality of
bolt installation. And awakened interest in high rockbolting as well. In average rock conditions
the units allow to install up to ten 6-metre-long wire cablebolts per shift.
4. APPLICATIONS OF STRING WIRE CABLEBOLTS IN POLISH COAL MINES
The use of long bolting is relatively widespread in the Polish coal mining industry.
Some of the possible applications for string wire cablebolts are named below:
As reinforcement for tunnel steel set support, to eliminate the need for additional props. This is
particularly important at face-gate intersections, where no space for additional props is available.
Copyright 2006 Taylor & Francis Group plc, London, UK
57
To reinforce tunnel support in areas of increased stress related to geological disturbances and
mining.
To reinforce set support of a face gate maintained along the bog-side of a longwall.
To reinforce tunnel support at tunnel crossings.
To reinforce cross-gate support to facilitate equipping and starting of a new longwall.
To reinforce support of large-size excavations.
To reinforce support of capital excavations with expected long service life.
To reinforce tunnel support to facilitate the installation of overhead rail transport ways.
As a preventive measure in excavations located in areas with high seismic hazard.
As reinforcement in tunnels supported solely by rockbolt support.
To reinforce rockmass around mining excavations.
Concentration of extraction, high face advance rates and increasingly difficult geological con-
ditions in the mined coalbeds force engineers to look for new ways to support face-gate intersec-
tions. Until recently the role of support reinforcement was played by short steel rods used to anchor
the caps of the installed steel sets. It happened frequently, however, that the fracture zone around
the excavation was developed to such an extent that the length of the used rods was insufficient.
For several years now some mines, among them ZG Sobieski, LW Bogdanka and KWK Mur-
cki, use string wire cablebolts to support longwall-gate crossings and do away with prop support
in these areas. Anchoring or fastening of arch segments is usually done with the use of additional
V-section horizontal members of various lengths installed along the axis of the tunnel. Anchors are
installed through the especially for the purpose burned holes (Figure 4.1). The length of the bracing
and the spacing between the anchors depend on local mining and geological conditions as well as
the distance of the place of the installation from the face. Installation directly at the face requires
using short braces, extending only between two sets. This limitation does not apply where instal-
lation is done further away from the face. Bracing with lengths up to 6 metres may then be used.
0
.
9
m
~6.0 m
Mudstone
Mudstone
V-36 roof member anchored with
a single string wire bolt
Fig. 4.1. Anchoring the sets at face-gate
Copyright 2006 Taylor & Francis Group plc, London, UK
58
The described method allows maintaining longwall-gate intersections without resorting to the
use of additional prop support. Combined with a solid backfill rib it dramatically improves condi-
tions in gates located on the bog side of the longwalls (ZG Sobieski and KWK Murcki mines),
(Figure 4.2). Increased interest in high bolting to support cross-gates is recently observed (among
others LW Bogdanka, KWK Murcki, KWK Piast, KWK Polska-Wirek, KWK Bielszo-
wice mines). These excavations are characterised by extreme widths. No roof set members are
able to resist the loads exerted on them by roof rock without additional reinforcement.
The surface of the exposed roof of such excavations is frequently flat (for example where Pr,
Pro set types are used) and impossible to hold without additional support. The traditional method,
utilizing a line of props installed along the middle axis of the tunnel, hinders the subsequent long-
wall equipping process.
Two steel dowels
p
a
s
p
o
d
s
a
d
z
k
o
w
y
B
a
c
k
f
i
l
l
r
i
b
Fig. 4.2. The face-gate maintained behind the face
0
.
9
m
Imudstone
Mudstone
Coal
6,4 m
3
,
5
m
Figure 4.3. Cross-gate supported by a Pro set-anchor system
Copyright 2006 Taylor & Francis Group plc, London, UK
59
Replacing props with cablebolts is effective and improves comfort and safety of work of the
equipping crew. Anchoring roof members allows also reclaiming the used props and side members
of the installed arches.
Examples of anchor-set support of cross-gates are shown on Figures 4.3 and 4.4.
1
.
2
m
Mudstone
Mudstone
Coal
Steel prop
String wire cablebolt
Steel dowel
Steel set head member
5.7 m
6.0 m
4
-
6
m
3
,
5
m
Figure 4.4. Cross-gate supported by a rectangular set-anchor system
Figure 4.5 shows schematically the moment of the newly-installed mechanized longwall suppo-
rts stepping out of the cross-gate supported by rectangular sets and cablebolts. In this case the set
members are stripped successively after the loading of each mechanized support unit. Additional
advantage of the method is the flat roof of the excavation, which eliminates the need to use crib-
bing to fill the void between the supports roof shield and the roof.
Figure 4.5 (left). Mechanized support units stepping out from cross-gate
supported by rectangular sets and cablebolts
Figure 4.6 (right). Longwall traversing a tunnel
Copyright 2006 Taylor & Francis Group plc, London, UK
60
String wire cablebolts
Figure 4.7. The method of tunnel crossing support used at ZG Sobieski mine
String wire cablebolt
Two steel rod bolts
Gate B-419
Haulage gate
Figure 4.8. The method of tunnel crossing support used at KWK Murcki mine
High rockbolting serves similar purpose when longwall face traverses a tunnel or enters the
cross-gate developed at the furthermost limit where the longwall equipment is to be stripped (Figu-
re 4.6). In this case the extent of the fractured zone around the excavation is so wide that only the
use of long bolts can guarantee their stability.
For the same reasons string wire cablebolts may be used to reinforce support installed at tunnel
crossings and break-ways. Examples of bolting patterns at tunnel crossings are schematically
shown on Figures 4.7 and 4.8 (KWK Murcki and ZG Sobieski mines). Figure 4.9 shows a me-
thod to support a two-level crossing developed at KWK Murcki mine. In all the presented ex-
amples cablebolts were used mainly to eliminate or substantially limit the number of prop supports,
which increased the working space for the crew and transport and haulage equipment. Cablebolts
installed in the roof substantially limit the extend of the fracture zone and increase the life of the
excavations reinforced in this way.
Copyright 2006 Taylor & Francis Group plc, London, UK
61
String wire cablebolts
String wire cablebolts
Figure 4.9. The method of support of a two-level crossing at KWK Murcki mine
5. FINAL REMARKS
The practical examples of utilization of long bolting presented in this paper certainly do not cover
all the aspects related to this topic. We are aware that the mining industry creates opportunity for
many more applications for string wire and rope cablebolts. One of such applications, at the mo-
ment not yet experimented with in Poland, seems to be the use of cablebolts as barrier support in
conditions of high seismic hazard. The research currently under way seems to show that high bol-
ting reduces the intensity of seismicity, hence increasing work safety of the crew. The authors hope
that this paper, by summarising the achievements in and applications of high bolting, will serve to
promote this support system and stem research into new designs and applications.
REFERENCES
Rataj M. 1995: Poszukiwanie nowych rozwiza obudowy kotwiowej na przykadzie firmy ANI Arnall, naj-
wikszego producenta obudw kotwiowych w Australii. Materiay Szkoy Eksploatacji Podziemnej 95,
str. 165172.
Rataj M. 1998: Zastosowanie kotwi linowych w australijskich kopalniach wgla. Materiay konferencyjne
II Konferencja pt. Obudowa kotwiowa jako skuteczny sposb zabezpieczania wyrobisk grniczych,
wieradw, str. 175.
Stilborg B. 1994: Professionals Users Handbook for Rock Bolting. Series on Rock and Soil Mechanics, Trans-
Tech Publications.
Windsor C.R. 1992: Cable Bolting for Underground and Surface Excavations. Rock Support, (Kaiser P.K. and
McCreath D. eds), Rotterdam: A.A. Balkema, pp. 349376.
Copyright 2006 Taylor & Francis Group plc, London, UK
63
Contribution to the Prediction of Ground Surface Movements
Caused by a Rising Water Level in a Flooded Mine
Anton Sroka
Technische Universitt Bergakademie Freiberg. Freiberg, Germany
ABSTRACT: Survey measurements indicate that ceasing to pump mine water out of closed coal
mines results in surface heave. In the practice of the German mine surveying the mathematical
basis for the projection of such ground movements are the formulae derived by Pttgens (1985)
and Fenk (2000).
The paper critically evaluates these methods and proposes modifications from the rock enginee-
ring and the science of mining damages points of view.
1. INTRODUCTION
Ceasing to pump water from closed underground coal mines usually results in the mine water level
rising considerably over vast expanse. The research conducted in many countries showed that such
planned flooding of mines leads to noticeable heave of the ground surface over large areas e.g.
[Oberste-Brink 1940, Pttgens 1985, Fenk 2000, Goerke-Mallet 2000, Baglikow 2003, Heitfeld et
al. 2004, Mhlenbeck 2005, Sroka 2005]. The question that needs to be asked is whether the up-
ward movement of the previously sinking ground surface may have a negative impact on the sur-
face buildings, i.e. whether it may lead to mine damages.
Heitfeld et al. (2004) noticed that during flooding of Sophia-Jakoba mine in the Aachen region
a discontinuity formed near the outcrops of Meinweg and Ruhrrand faults. The resultant surface
deformation was 89 km long and 67 cm high. About 110 damages to structures, 9 of them be-
yond repair, were observed along this line.
The question regarding possible mine damages is closely connected to the question whether it is
possible to predict ground movements caused by mine flooding. This presents researchers with a task
to describe and analyse the current state of knowledge relating to the aforementioned problem. The
obtained results must be then generalised and extrapolated to apply in the specific rock-mass and
mining conditions characteristic to the examined coal mining region.
2. THE PTTGENS METHOD (1985, 1998)
Pttgens says that the zone of broken rock formed in the roof of a mined-out coal-seam is chara-
cterised by increased porosity and permeability. A rise of the water level in a flooded mine leads to
the increase of pressure in the caving zone and results in its expansion due to a decreased effective
vertical stress. Subsequently the hydrostatic force acting on the solid strata above the caving zone
heaves the rock right up to the surface. Pttgens claims that the extent of the movements depends
on the increase of mine water level, mining depth, seam thickness and the extend of extraction.
Copyright 2006 Taylor & Francis Group plc, London, UK
64
The mathematical foundation of the Pttgens method is the Geerstma solution (1973), dealing
with surface subsidence related to natural gas exploitation.
The solution was based on the poroelastic theory and two assumptions:
Linear relationship between stress and strain, and
Homogeneity of the whole rock-mass.
The convergence or compaction of the porous gas-bearing strata i.e. decrease in its thickness
due to gas extraction resulting in decreasing its pressure, can be calculated after Geertsma from the
following formula:
p M c M
m
A = A (1)
where AM deposit compaction [m]; M deposit thickness [m]; Ap drop in the pore pressure
[Bar], and c
m
compaction coefficient [Bar
-1
].
Compaction coefficient c
m
is equal to the relative change of deposit width at a pressure drop of
1 Bar. Generally speaking its value depends on rock-mass characteristics, porosity, initial pore pre-
ssure and deposit depth.
According to Geertsma a unit subsidence of a point on the surface can be defined for an infini-
tely small unit volume of deposit dV from the following equation:
( )
dV
H r
H ) 1 ( c
) H , r ( ds
2 / 3
2 2
m
+
t
u
= (2)
where ds (r, H) unit subsidence of a point on the ground surface; u Poisson ratio; H depth of
the deposit unit volume, and r horizontal distance between the point and the deposit unit volume.
The solution is known in literature under the name of the nucleus of strain concept.
For a circular gas reservoir with radius R the formula to calculate the maximum surface subsi-
dence is as follows:
p M c
C 1
C
1 ) 1 ( 2 M
C 1
C
1 ) 1 ( 2 ) H , R ( s
m
2 2
max
A
(
(
+
u = A
(
(
+
u = (3)
where C = H/R.
The point of maximum subsidence is located above the centre of the circular gas deposit. The
equation (3) was used by Pttgens as the basis for his method of calculating the ground surface
heave caused by a rising level of water flooding a mine.
As already mentioned Pttgens assumed that an increase of pressure in the zone of broken roof
rock leads to its elastic expansion. Such expansion was defined as the opposite to compaction:
p h d h
m
A = A
where Ah broken roof rock zone expansion; h thickness of the broken roof rock zone; Ap
pressure increase caused in the zone of broken roof rock by the rising level of water flooding the
mine; d
m
expansion coefficient of the broken roof rock.
The vertical and horizontal positions and the extent of the zone of broken roof rock were de-
fined by the following assumption:
The zone is situated on the exploitation level.
The height of the zone is equal to four widths of the mined-out coal seam, and
The horizontal extent of the zone is equal to the extent of the extraction.
The solution describing the maximum heave of a point on the ground surface situated above
a circular reservoir (i.e. a circular mining working filled with water) proposed by Pttgens is ana-
logical to that proposed by Geertsma (3).
Copyright 2006 Taylor & Francis Group plc, London, UK
65
It assumes the following form:
p h d
C 1
C
1 ) 1 ( 2 h
C 1
C
1 ) 1 ( 2 ) H , R ( h
m
2 2
max
A
(
(
+
u = A
(
(
+
u = (4)
It is possible to assess the heave of any point on the surface from a graphical representation of
the function:
( )
2
C 1
C
1 C y
+
= (5)
In the case of extraction of a larger number of coal seams the principle of linear superposition
with a constant value of the broken roof rock expansion coefficient d
m
is applied [Wings, Misere,
Pttgens 2004].
Pttgens estimated the value of the broken roof rock expansion coefficient for the Limburg (Hol-
land) coal mining region from in-situ subsidence and heave measurements.
The value is equal to d
m
.
. MN / m 10 35 , 0 Bar 10 35 , 0 d
2 2 3
m
= =
Some of the measurements are shown in Figure 1.
Figure 1. Heave and subsidence of the surface along the Aalbeek-Hoensbroek-Schinveld profile
The above figure shows ground surface heave measurements done in the years 19742003 and
the mine subsidence calculated for the years 19151974.
Pttgens points out the seemingly evident similarity of both graphs. Moreover, an analysis of
the results of the measurements combined with the geological and mining profile showed that not
only discontinuous subsidence but also discontinuous heave of the surface is possible in the area of
the faults surface outcrops.
Copyright 2006 Taylor & Francis Group plc, London, UK
66
According to Pttgens this is possible only if extraction at each side of the fault was conducted
at different rate and even so only in places where discontinuous subsidence previously occurred.
The relationship between the maximum* subsidence and the maximum* heave may for all prac-
tical purposes be presented by the constructed by Pttgens (1985, 1998) nomograph (Fig. 2) (*the
authors interpretation).
The nomograph allows calculating the expected maximum ground surface heave as a percent-
age of the previous maximum subsidence if the following are known: extraction geometry (GA),
mine water geometry (GG) and the mining depth.
According to Wings, Misere and Pttgens (2004) it is particularly difficult to properly assess
GA and GG.
The example shown in Figure 2 pertains to mining depth of 800 m. The extraction geometry is
theoretically equivalent to a circle with a radius equal to the radius of influence of the mining (GA =
= 100%) and the water horizon will rise right to the surface.
Figure 2. Heave-subsidence relationship
Calculated from the nomograph maximum heave is approximately equal to 5% of the maxi-
mum subsidence. According to Pttgens (1998) the maximum ground surface heave may reach 25%
of the maximum subsidence if the mine water horizon reaches the ground surface.
Upward movements of this magnitude do not normally cause any damage to surface structures.
Pttgens sets the limit value of deformations, which may result in damages at 0,25 mm/m in com-
pression or tension. This obviously does not pertain to outcrops of faults on the surface and the
areas of previous discontinuous deformations.
3. ANALYTICAL METHOD OF FENK (2000)
Similarly to Pttgens, Fenk based his method on the assumption that a rise in the water level in a
flooded mine leads, due to hydrostatic upward pressure, to an elastic expansion of the zone of
broken rock situated in the roof of the mined-out coal seam.
Copyright 2006 Taylor & Francis Group plc, London, UK
67
The formula for calculating the upward movement of the rock-mass he came up with following
the assumption is as follows:
e
(
= s W 1
H
E
d h
G
s
w m
(6)
where: h maximum heave of the ground surface, at the moment when the level of water flooding
the mine stabilizes [m]; d
m
expansion coefficient of the broken roof rock [m/MN];
w
specific
gravity of water [0,01 MN/m];
G
specific gravity of the rock-mass [MN/m]; H mining depth;
E
S
modulus of rigidity of the broken roof rock [MN/m]; W increase of the water level in the
rockmass [m], and s
e
maximum surface subsidence caused by the extraction [m].
In the case when W = H, i.e. the water level rises up to the ground surface, the formula (6) as-
sumes the following form:
e
= s H
E
d h
G
s
w m
(7)
The analytical solution presented by Fenk was derived for a so-called one seam model. According
to the author it only allows calculating the maximum possible heave of the ground surface and con-
trary to the Pttgens solution does not require determining the height of the caving zone. Because of
that it is impossible to compare d
m
coefficients appearing in the Pttgenss and Fenks formulae.
Formula (7) has been verified by Fenk (2000) for the Zwickau coal mining region by in-situ
measurements of the surface heave.
The verification was based on a regression function, which is the theoretical solution of (7) and
has the form:
( )
e
= s H b a h .
The unknowns a and b were established from the measured values of heave and maximum sub-
sidence and the depth of mining under the 26 measuring stations. Due to the fact that more than
one seam was mined in the Zwickau mining region, each time Fenk accepted the mining depth as
equal to the depth of the deepest coal seam (!).
After applying the regression analysis Fenk derived the following formula:
( )
e
= s H 2024 024 , 0 h (h w mm; H , s
e
w m) (8)
Hence the values of the geo-technical parameters d
m
and E
s
are equal to:
, MN / m 10 24 , 0 d
2 2
m
=
. m / MN 023 , 0 dla m / MN 6 , 46 E
3
G
2
s
= =
The measuring stations were individually assigned corresponding values of mining depths as equal
to that directly underneath them. For this reason the measuring stations located outside the mining area
were excluded from the assessment, as assigning them any value of mining depth was impossible.
In order to calculate the value of upward movement of the points situated outside the mined-out
area Fenk proposed the following empirical formula:
e
= s p h (9)
Specifically for the Zwickau mining region Fenk proposed the following equation based on the
data collected there in 1997:
e
= s 1 , 58 h (h w mm;
e
s w m) (10)
According to Fenk (2000) the advantage of his solution over Pttgenss is that the analyst does
not have to possess the knowledge of the height of the caving zone, the geometry of the conducted
mining and the interaction between individual excavations to be able to calculate the extent of hea-
ve of the ground surface.
Copyright 2006 Taylor & Francis Group plc, London, UK
68
4. CRITICAL REMARKS REGARDING PTTGENSS AND FENKS METHODS
The methods proposed by Pttgens and Fenk are the first analytical methods used to predict the
amount of ground surface heave caused by water flooding a mine.
Both methods are based on the relationship between the amount of heave and the amount of the
previously occurred mining-related surface subsidence. However, in the case of the Pttgenss me-
thod, combining the two is not necessary (see formula 4) and the author alleges that the connection
was created purely and somewhat artificially to afford possibility to estimate the expected heave
from the earlier subsidence.
It may appear at the first glance that a relationship between the two graphs shown in Figure 1
really exists. A careful analysis shows, however, that the places of maximum subsidence and heave
do not correspond. The same applies to other characteristic parameters of the graphs, e.g. the points
of their maximum slope.
Creating a ground heave distribution graph from subsidence distribution data will lead to ob-
taining substantially erroneous, both quantitatively and qualitatively, result.
The recent works of Goerke-Mallet (2000) and Mhlenbeck (2005) contradict the claim that
any relationship between the amounts of ground heave and previous subsidence exists.
According to Pttgens (1985) heave is primarily calculated from formula (4). The graphical in-
tegration grid produced by solving the formula allows establishing the value of heave at any point
on the surface provided that the input data i.e. the extraction geometry, the coal seam width, the
mining depth and the mine water level are known. It is also necessary to know the values of the
expansion coefficient d
m
and Poissons ratio u.
For the cases of multiple coal seam extraction Pttgens assumed the applicability of the linear
superposition principle, which in the form presented by him seems to be very complicated and time
consuming and from the theoretical point of view at least doubtful.
Comparative calculations of the maximum heave and the maximum subsidence reveal signifi-
cant differences between the subsidence influence (Pttgenss based on Grond (1957)) and the hea-
ve influence functions.
For example, for a circular extraction area with dimensions R = H and the angle of draw = 50
gon, the maximum subsidence, located directly above the middle of the circular excavation, is eq-
ual to the so-called full subsidence (GA = 100%). However, the ground heave calculated for this
point is only equal to 29,3% of the maximum heave, obtainable theoretically for an infinitely large
area extracted. 90% of the maximum heave is obtained for R = 10 H and 95% for R = 20 H.
This means that the graphical integration grid for calculating the heave is much larger than that
for calculating subsidence, and the surface heave extends over a much greater area than the mine
subsidence. The Gronds method for calculating the relative subsidence over circular mined-out
areas (GG) and the Geertsmas method for calculating the relative heave (GG) served Pttgens to
create a nomograph shown in Figure 2.
The nomograph was constructed for the following data:
Calculating mine subsidence Gronds method for = 50 gon.
Calculating heave d
m
= 0,3510
-2
m
2
/MN, u = 0,25, h = 4 M.
Pttgenss statements and assumptions are in a number of points absolutely logical. His quanti-
tative assumptions are, however, questionable.
One of them is the assumption that the height of the caving zone in the roof is always four ti-
mes the width of the extraction. According to the results of comprehensive research presented by
Borecki and Chudek (1972), the width of the zone falls in the range 320 times the mining width.
Also the value of the Poissons ratio is generally accepted as constant at 0.25.
For these reasons the value of the expansion coefficient of d
m
= 0,3510
-2
m
2
/MN must be per-
ceived only as applicable to the assumed values of the caving zone height and the Poissons ratio.
Copyright 2006 Taylor & Francis Group plc, London, UK
69
An analysis of the Geertsma solution (1973) shows that for a very large gas deposit the maxi-
mum value of surface subsidence is greater than compaction of the deposit itself. The ratio of the
two values equals 2 (1 u) meaning that for u = 0.25 the maximum subsidence is approximately
50% greater than the compaction of the gas deposit.
In his elaborations on the problem Pttgens does nor correct this error.
Simulations showed that applying Pttgenss method to multiple seam mining conditions re-
sults, due to extensive differences between the subsidence influence function and the heave influ-
ence function, in producing heave distributions substantially different from the calculated subsi-
dence distributions. The discrepancies were seen when comparing such characteristic values as e.g.
the maxima.
In the authors opinion the Fenks method (2000) is useful only in the case of single coal seam
extraction and then for estimating the maximum value of ground heave.
The expansion coefficient d
m
used in the Fenks formula can in no way be compared to the sym-
bol expansion coefficient represented by the same in the Pttgenss equation.
5. MATHEMATICAL MODEL FOR PREDICTING GROUND SURFACE MOVEMENTS
CAUSED BY RISING WATER LEVEL IN A FLOODED MINE
The formulae presented below are derived for a finite caving zone element with a base of Ax Ax and
a height h.
Due to rising of the mine water level the height of the zone increases by Ah. This leads to in-
finitesimal heave of the ground surface (Figure 3).
Figure 3. Graphical representations of mathematical equations
Assuming equal cause and result volumes the following solution of Geertsmas influence func-
tion is derived:
( )
2 / 3
2 2
m
H r
H
2
V p d
) r ( h
+
t
A A
= A (11)
where AV = Ax
2
h; Ap = (z
w
z
F1
)
w
, and . p
H
V
2
d
h
2
m
max
=
Copyright 2006 Taylor & Francis Group plc, London, UK
70
Division of the extracted areas into elements allows calculating heave values for any geometry
of mining voids by means of linear superposition. Analyses done for multiple seam mining condi-
tions revealed the existence of strong dependence of the expansion coefficient d
m
on mining depth
and width.
Sroka (1990) and Hejmanowski (1993) applied Knothes theory to predict subsidence caused
by gas and oil extraction. The results of the calculations showed a high level of correlation with the
results of in-situ measurements [Sroka 1990]. The subsequent analysis of the parameters of Kno-
thes theory permited to establish the value of the draw angle |. For Grningen gas deposit this
value equals to approx. 7 gon, from which follows the value of the subsidence influence radius as
R
K
= 9H. Recalculating these values into the values characteristic for a so called Ruhrkohle method
results in obtaining the so called radius of full subsidence through equal to 11H, i.e. 11 times grea-
ter than the mining depth. The results confirm the validity of the Geertsma solution.
Applying Knothes theory results in obtaining the formula (12):
|
|
.
|
\
|
t
A A
= A
2
K
2
2
K
m
R
r
exp
R
V p d
) r ( h (12)
| = cot H R
K
,
p tan
H
V
d h
2
2
m max
A |
A
= A .
Formulae (11) and (12) can be used alternatively to calculate the surface heave for known ex-
traction geometries. Simulations ran for multiple seam mining conditions showed disparity, in so-
me cases substantial, between the obtained predicted subsidence and heave profiles. The presented
method allows taking into account numerous geometrical variables of the extracted coal seam di-
mensions as well as to introduce a mining-depth-dependent expansion coefficient. Further develop-
ment of the analytical method is determined, however, by extensive research into and in-situ mea-
surements of surface heave caused by mine flooding.
REFERENCES
Baglikow V. 2003: Bergschden nach Beendigung der Grubenwasserhaltung im tiefen Bergbau. Markscheide-
wesen 110, 2003, Heft 2, s. 4549.
Bals R. 1931/32: Beitrag zur Frage der Vorausberechnung bergbaulicher Senkungen. Mitteilungen aus dem
Markscheidewesen, 42/43, s. 98111.
Borecki M., Chudek M. 1972: Mechanika Grotworu. Wydawnictwo lsk, Katowice.
Fenk J. 2000: Eine analytische Lsung zur Berechnung von Hebungen der Tagesoberflche bei Flutung unter-
irdischer Bergwerksanlagen. Das Markscheidewesen 107, 2000, Heft 2, s. 420422.
Flaschentrger H. 1938: Die Kostenverteilung bei gemeinsam verursachten Bergschden im Ruhrgebiet. Mit-
teilungen aus dem Markscheidewesen, 49, s. 95137.
Geertsma J. 1973: A Basic Theory of Subsidence Due to Reservoir Compaction: the Homogeneous Case. The
Analysis of Surface Subsidence Resulting from Gas Production in the Groningen Area, the Nederlands.
Editor Nederlandse Aardolie Maatschappinj B.V., Deel 28, 1973, s. 4361.
Geertsma J., van Opstal G. 1973: A Numerical Technique for Predicting Subsidence above Compacting Re-
servoirs, Based on the Nucleus of Strain Concept. The Analysis of Surface Subsidence Resulting from
Gas Production in the Groningen Area, the Nederlands. Editor Nederlandse Aardolie Maatschappinj B.V.,
Deel 28, 1973, s. 6378.
Goerke-Mallet P. 2000: Untersuchungen zu raumbedeutsamen Entwicklungen im Steinkohlenrevier Ibben-
bren unter besonderer Bercksichtigung der Wechselwirkungen von Bergbau und Hydrogeologie. Praca
doktorska RWTH Aachen, Wydawnictwo Mainz.
Goerke-Mallet P., Preue A. und Coldewey W.G. 2001: Hebungen der Tagesoberflche ber betriebenen und
gefluteten Bergwerken. 43. Wissenschaftliche Tagung des Deutschen Markscheider Vereins, Trier 2001,
s. 125139.
Copyright 2006 Taylor & Francis Group plc, London, UK
71
Grond G.J.A. 1957: Ground Movements Due to Mining. Colliery Engineering 34, s. 157197.
Heitfeld M., Rosner P., Mhlenkamp M. und Sahl H. 2004: Bergschden im Erkelenzer Steinkohlenrevier. 4.
Altbergbaukolloquium, Leoben 2004, s. 281295.
Hejmanowski R. 1993: Zur Vorausberechnung frderbedingter Bodensenkungen ber Erdl- und Erdgasla-
gersttten. Praca doktorska, TU Clausthal, 07.07.1993.
Hejmanowski R., Sroka A. 2000: Time-Space Ground Subsidence Prediction Determined by Volume Extrac-
tion from the Rock Mass. Proceedings of the Sixth International Symposium on Land Subsidence, Ra-
venna, Italia, 2429 September 2000, Vol. II, s. 367375.
Mhlenbeck H. 2005: Hebungen nach Einstellung von Wasserhaltungen im Ruhrgebiet. Markscheidewesen 112,
2005, Zeszyt 3, s. 97102.
Oberste-Brink K. 1940: Die Frage der Hebungen bei Bodenbewegungen infolge des Bergbaus. Glckauf 76,
s. 249256.
Pttgens J.J.E. 1985: Bodenhebung durch ansteigendes Grubenwasser. 6. Internationaler Kongress fr Mark-
scheidewesen, Harrogate 1985, s. 928938.
Pttgens J.J.E. 1998: Bodenhebung und Grundwasseranstieg aus geotechnischer und markscheiderisch-geo-
dtischer Sicht im Aachen Limburger Kohlenrevier. Freiberger Forschungshefte, Bergbau und Geotech-
nik A 847, s. 193207.
Sroka A. 1990: Studie zur Analyse und Vorhersage der Bodensenkungen und des Kompaktionsverhaltens des
Erdgasfeldes Groningen/Emsmndung. Ekspertyza, Goslar 1990.
Sroka A. und Fenk J. 2003: Studie zu Auswirkungen des Grubenwasseranstiegs auf Bewegungen der Tages-
oberflche in stillgelegten Bergbaubereichen des Saarlandes. F+E Bericht der TU Bergakademie Freiberg.
Sroka A. 2005: Ein Beitrag zur Vorausberechnung der durch den Grubberwasseranstieg bedingten Hebungen.
5. Altbergbaukolloquium, Clausthal-Zellerfeld 2005, s. 453464.
Wings R.W.M.G., Misere W.M.H. und Pttgens J.J.E. 2004: Bodensenkung Bodenhebung Bergschden?
44. Wissenschaftliche Tagung des Deutschen Markscheider Vereins, Bochum 2004, s. 258269.
Copyright 2006 Taylor & Francis Group plc, London, UK
73
A Geophysical Model for the Analysis of Seismic Emissions
in the Area of Rock Beam Splitting Caused by Mining Operations
Henryk Marcak
AGH University of Mining and Metallurgy, Institute of Geophysics. Krakw, Poland
ABSTRACT: In Poland the rock-burst risk in underground mines is an important issue in the mi-
ning industry. The predictive rock-burst information is extracted from different kind of information
based on measurements. One of them is seismic data, but the predictive information included in it
is limited, due to small number of strong seismic shocks and the chaotic development of rock de-
struction before the strong seismic events occur. In the paper, a geomechanical model of rock de-
struction prior strong seismic shocks is proposed. It can help in an optimal organization of the mea-
surement system and its interpretation for better assessment of the seismic risk in underground mines.
INTRODUCTION
There are two mining regions in Poland subject to rock-burst risk; one is the coal basin in Upper
Silesia, the second, the copper region in Lower Silesia. In both regions strong sedimentary layers
appear in the roof of exploited parts of the rock masses. They consist of sandstone in the coalmine
basin and limestone and dolomite in the copper mining areas.
The significant influence of strong sedimentary layers on building up stresses in exploited rock
masses has been considered and accepted many years ago [Saustowicz 1965]. Calculation of elastic
deformations, results of rock-beam bending in roof of exploited seam, has been used for modelling
stresses in exploited areas. In this paper the stresses produced by rock-beam in the roof of an ex-
ploited area are also considered, but its plastic behaviour is analysed.
In the paper it is shown, that both, seismic data, and the results of measurements of changes in
ratio of vertical and horizontal deformations of roof layers over the exploited area, can be used in
assessment of strong seismic hazard. In the simplest attitude the deformation of roof layer can be
considered on the base of elastic theory. The exploited area is here described in a two-dimensional
system, the x-axis being perpendicular to the long wall and the z-axis being the depth.
The stress-distribution may thus be described by the equation:
0
z
z) F(x,
z x
z) F(x,
2
x
z) F(x,
4
4
2 2
4
4
4
=
c
c
+
c c
c
+
c
c
(1)
where F
j
(x, z) is a function, which defines the elements of the stress tensor.
= o ) z , x (
x
2
2
z
) z , x ( F
c
c
= o ) z , x (
z
2
2
x
) z , x ( F
c
c
(2)
t
xz
(x, z) =
x z
) z , x ( F
2
c c
c
,
Copyright 2006 Taylor & Francis Group plc, London, UK
74
and
)
`
o t
t o
= t
z xz
zx x
j , i
(3)
The components of the displacement vector (u, w) fulfill the relation:
2
2
2
2
z
) z , x ( F
x
) z , x ( F
) 1 (
x
) z , x ( u
G 2
c
c
v
c
c
v =
c
c
(4)
2
2
2
2
z
) z , x ( F
z
) z , x ( F
) 1 (
z
) z , x ( w
G 2
c
c
v
c
c
v =
c
c
(5)
where G, v are the stiffness coefficient and Poissons ratio respectively.
Let the roof layer in the vicinity of the exploited area have a thickness h and elastic coefficients
G
1
, v1. Harder rocks overlying the seam are treated as a half-plane with elastic coefficients G
2
, v
2
.
This model allows obtaining an elastic model of energy concentration before rock burst. More so-
phisticated models can be found in literature e.g. [Salomon and Hunro 1967, Bieniawski 1976,
Brady and Brown 1993]. The general idea in this approach is to model the stress field and the
associated field of deformation by partial differential equations. The equations solved with a finite
element method, provide the stress distribution for complex geological structures. There are how-
ever, limits in the assessment of real rock deformations by such calculations. Firstly, it is difficult
to obtain true initial and boundary conditions. Furthermore, an unstable behaviour of rocks after crack-
ing and fracturing of its masses cannot be properly modelled with equations that assume continuity
of the medium.
In polish underground rock bursting mines, recently, a relation between deformation of bore-
holes located in the roof of the productive layer and emission of strong seismic events has been ob-
served. Sudden jumps of deformations were recorded 5 to 20 days before strong tremors occur.
[Orzepowski 1998, Matwiejszyn and Ptak 2000]. Measurement of roof deformation by Orzepowski
in coppers mine, with sensors measuring deformations of a borehole diameter in fixed horizontal
directions and at fixed depth has been carried out over several years. For example, [Orzepowski
1998] obtained results with sensors installed in four boreholes 30 m above the roof of the ore seam
in a copper mine. The results of measurements prior to a strong seismic event are shown in Figure 1.
The distance between the shock epicentre and the epicentre azimuth has been estimated and is
shown in Figure 1 caption. It can be seen that up to 25.01.97 changes in the borehole diameter were
slight. Later deformations of all boreholes exhibited a similar, time varying pattern, until 11.02.97,
when a strong mining shock occurred. There is no deterministic relation between anomalous defor-
mations in roof beam and strong seismic occurrence. Only regional deformations (those which are
registered with larger number of sensors) are used for prediction. Even, using such criteria only the
hazard of strong seismic shocks markedly increases after registration of roof deformations It is si-
gnificant, that after registration of anomalous deformation, failure of prediction have form of strong
events lack after roof deformation not lack of roof deformation before strong seismic shock. For
explanation of this results the model of brittle deformation of rocks should be used.
These measurements indicate, that the inelastic deformations have appeared in roof beam. It can
be expected, that it have form of plastic horizontal deformations, which can turn into splitting and
bending downward of roof layer. The horizontal stresses appear in specific time, when the inelastic
deformation appears.
The deformations in strong rocks overlying splitting are the best explained by the use of a mo-
del of inelastic deformation. For this purpose Rudnickis dilatancy-plastic model is suitable [Rud-
nicki and Rice 1975]. The plastic deformations
p
and c
p
are associated with shear stress t and com-
pressive pressure p, the elastic deformations
e
and c
e
are linked with the shear stress t and the
compressive pressure p through the coefficients G and K.
u
Copyright 2006 Taylor & Francis Group plc, London, UK
75
A
70
60
50
40
30
d
e
f
o
r
m
a
t
i
o
n
s
,
m
65
55
45
35
B
70
60
50
30
d
e
f
o
r
m
a
t
i
o
n
s
,
m
C
60
50
40
30
d
e
f
o
r
m
a
t
i
o
n
s
,
m
55
45
35
20
10
d
e
f
o
r
m
a
t
i
o
n
s
,
m
40
30
D
1
1
.
1
2
.
9
6
1
1
.
0
1
.
9
7
1
0
.
0
1
.
9
7
1
5
.
0
1
.
9
7
2
0
.
0
1
.
9
7
2
5
.
0
1
.
9
7
3
0
.
0
1
.
9
7
4
.
0
2
.
9
7
9
.
0
2
.
9
7
1
4
.
0
2
.
9
7
1
9
.
0
2
.
9
7
4
.
0
2
.
9
7
1
.
0
3
.
9
7
6
.
0
3
.
9
7
1
1
.
0
3
.
9
7
1
6
.
0
3
.
9
7
2
1
.
0
3
.
9
7
2
6
.
0
3
.
9
7
3
1
.
0
3
.
9
7
(date)
Figure 1. Diameter changes of four boreholes before a strong mining event,
which occurred on 11.02.97. in a copper mine.
a) borehole located 26 m from epicenter, azimuth 83, b) borehole located 95 m from epicenter, azimuth 15,
c) borehole located 288 m from epicenter azimuth 78, d) borehole located 301 m from epicentre azimuth 60
(from [Orzepowski 1998])
This model assumes that:
) p (
h
1
p
- -
t = (6)
-
t =
G
1
e
(7)
p
K
1
e
= c (8)
p p
| = c (9)
{ }
ij
tr p o = (10)
where
t c
t c
= t , G, K, h are the stiffness, elasticity and plasticity modules respectively, | is the dila-
tancy coefficient, is the coefficient of internal friction.
Copyright 2006 Taylor & Francis Group plc, London, UK
76
If the energy release is dissipated in the form of heat, and is, therefore, ignored and a loss of
cohesion is assumed the following equation may be written.
Y
p
+ qp
p
= t
p
- pc
p
(11)
where Y is the yield plastic coefficient, q is the coefficient of friction, resulting from the movement
of the rock-blocks.
Taking into account the relationship:
c
p
= |
p
(12)
The result of transformation formula (11) and (12) the following formula is obtained:
t = Y + (| + q)p (13)
and plastic yielding is described by
Y = t - p (14)
where is the coefficient of internal friction.
The coefficient q is related to the petrography structure of the rock and the quality of its cemen-
tation. The second coefficient is related to dilatancy. This term means the increase of the vo-
lume of the rock resulting from inelastic deformation. The dilatancy may introduce a negative effect
to the coefficient. It distorts the relationship between coefficient q and the total stress q
p
. The fo-
llowing relationship now emerges:
h = p
c
g
Y
g
t
p p p
c
c
+
c
c
=
c
c
(15)
The last equation can be used to describe the bending layer. Movement between the upper and
lower planar surfaces of the roof layers is assumed to be a consequence of bending.
0
h
k
h
) ( ) (
2 2
s c
o
+ + | c c
The relationship in this beam, based on the dilatancy-plastic model suggested by [Nikitin and
Ryzak et al. 1977, 1984] is:
0 c
o
h
k
h
) (c ) (c
2 2
s + + (16)
Applying Leganders method for maximalization left side of equation (16) gives a value for h
o
.
It can be shown [Nikitin and Ryzak et al. 1997, 1984] that from h
o
value the parameters of inclu-
sion zones can be estimated.
a
b
c
Figure 2. Parameters describing the inclusion zone in the roof rock beam.
a) perpendicular fracturing, b) parallel fracturing, c) parameters of inclusion zone b = x/l
Copyright 2006 Taylor & Francis Group plc, London, UK
77
Those zones have the following properties:
Constrain a number of fractures whose length distribution has fractal properties parallel to inc-
lusion zone boundary or parallel to them.
They are repeated at regular intervals.
The fracture zones have the following parameters (Fig. 2):
( -|)
2
=
k
h 4
0
(17)
o = |
.
|
\
|
+
|
G
h
1
2
0
(18)
k 3 G
k 2
2
1
+
o = + (19)
|
|
|
|
.
|
\
|
+
+
+
+
=
K
4
3
G
K
16
3
G
4
3
1 1
K
16
3
G
K
4
3
G
l
b
(20)
The direct application of formulas 1720 to seismic emission in underground mines [Marcak
1985, 1993], proved the relevance of this model for seismic emission in underground mines. The
intensity of seismic activity during exploitation has extremes located periodically and the distances
between them agree with values described by formulas 1720.
The inclusion zones are built up by rather complicated fracturing. However, appearance of very
strong shocks is rare, and the special mechanism describing their creation should be considered.
FRACTURING IN THE ROOF OF ROCK-BEAM
The rock-burst is a result of rock-mass movements in the fracture surface. The fracturing of the
rock overlying exploited seams starts when the difference between the greatest principal stresses o
I
and the smallest one o
III
have reached a maximum, called ultimate shear strength. In the case when
two dimensional stress distributions describe critical state of stresses, the limit condition has the form:
) ( f
2
| |
III I
max
o =
o o
= t (21)
The direction of fracture plane can be presented with use of two envelopes in Coulomb-Mohr
stress circles as is shown on Figure 3.
Envelopes can be described with equation:
u o + t = t
tan | |
0
(22)
where u angle of internal friction, t
0
the cohesive shear strength (depending on the lithological
properties of the rock), o
t
.
According the properties of formula (2) o
o + o
= o (24)
where t
v,H
shear stresses on the horizontal border.
H , v
2 2 H v
III I
)
2
( 2 t +
o o
= o o (25)
In result, taking into account, that o
I
+ o
III
is invariant in transformation of stresses, the stresses
o
.
where N
2
nitrogen concentration [%].
For the calculation of the indicator level, the regular measurements of the air composition in
different points of the longwall were taken while the fly ash suspensions with carbon dioxide were
being supplied to the gob.
The second part of the tests concerned the evaluation of the gob inertization influence on the
carbon dioxide concentration in the outlet stream of air. It is especially important in view of mining
regulations, which allow for a 1% concentration of CO
2
.
3. TECHNOLOGY OF FLY ASH SUSPENSION WITH CARBON DIOXIDE PRODUCTION
The suspensions are produced in a special installation (scheme on Figure 1). The energy wastes are
delivered from the Rybnik Power Plant to the mine by the rail transport and pneumatically
unloaded into the cisterns. The preparation of the suspensions is realised in the closed mixer, where
the ashes and water are added. While the ready suspensions flow to the shafts backfill pipeline, the
gaseous carbon dioxide is added. Now the hydraulic transport of mixtures into the gob begins.
The suspensions were delivered from the backfill pipeline, end parts of which were left into the
gob in the liquidated upper gateway. The distance between the end of the pipeline and the long-
wall head varied from a few to dozens meters, it was caused by local conditions. After some time
the end of the pipeline got plugged and the end parts of the pipeline were disconnected in an ac-
cessible place in the upper gateway.
Carbon dioxide was delivered to the mine in a liquid phase and stored in the tank of 30 tonnes
capacity. The gasification of CO
2
was conducted by an evaporating dish and an electric heating
unit. The efficiency of gasification installation reached 4,5 m
3
/min.
Copyright 2006 Taylor & Francis Group plc, London, UK
109
F
l
y
a
s
h
t
a
n
k
C
e
m
e
n
t
t
a
n
k
M
i
n
i
n
g
w
a
s
t
e
s
t
a
n
k
Water
Shatf
Mixer
Indirect
tank
CO cistern
2
Processing
plant
Figure 1. Scheme of the fly ash suspension with carbon dioxide production installation.
(Source: Borynia Coalmine)
4. THE PLACE OF THE TESTS
The in-situ tests have been conducted in the longwall D-31 in the Borynia Coalmine, which be-
longs to the Jastrzebie Coal Company in Southern Poland. The tests started in May 2004, and fi-
nished in April 2005.
The fly ash suspensions with carbon dioxide were supplied to the gob of the D-31 longwall,
which led in the upper layer of the 405/1 stratum. The length of the longwall [1] was from 140 m in
the beginning to 210 m in the end, the planned advance 600 m and the average height 3 m. The
lowest point of the longwall was located at the beginning of the lower gateway D-31 (first part of
the gob), the highest points at the end of the upper gateway D-31a. Both the gateways were li-
quidated behind the head of the longwall.
The longwall geometry helped to the carbon dioxide in migration into deeper parts of the gob.
The analyses of the air composition (especially concentration of oxygen, carbon dioxide, car-
bon monoxide, nitrogen and methane) were carried out in measurement stations, which location is
showed in Figure 2. On this basis the Graham indicator and the carbon monoxide indicator were
calculated.
Denotation and location of measurement station:
IN inlet airflow;
S10, S20 and S70 station located behind 10, 20 and 70 section of the hydraulic support;
FS first section of the hydraulic support (crossing of the head of the longwall and the upper
gateway);
GL D-31 upper gateway liquidation;
OU outlet of the airflow.
Copyright 2006 Taylor & Francis Group plc, London, UK
110
1
.
0
5
.
0
4
r
.
1
.
0
7
.
0
4
r
.
1
.
0
6
.
0
4
r
.
1
1
.
0
5
.
0
4
r
.
2
1
.
0
5
.
0
4
r
.
1
1
.
0
6
.
0
4
r
.
2
1
.
0
6
.
0
4
r
.
-633,31 m
-636,41 m
-650,03 m
Longwall D-31
Strata 405/1
m
a
x
.
6
m
Reserve backfilling
pipeline 150 mm
Backfilling pipeline
150 mm
Inlet airflow station
GL - gallery liquidation
LS - last section
S70 station
S20 station
S10 station
D
-
3
1
g
a
lle
r
y
D
-3
1
a
g
a
lle
ry
R
i
s
e
g
a
l
l
e
r
y
D
-
3
1
Outlet airflow station
Figure 2. Scheme of location of measurement station (flags) in the D-31 longwall.
(Source: Borynia Coalmine)
5. RESULTS THE EVALUATION OF THE FIRE PREVENTION METHOD
The frequency of application of the fly ash suspensions with carbon dioxide to the gob was very
irregular. It depended on spontaneous heating development, and was in accordance with the value
of fire indicators. During the tests there were some periods, when the heating process in the gob
advanced. Then the fire prevention intensified the quantity and frequency of the fly ash suspen-
sion with carbon dioxide supplied to the gob increased.
In Figure 3 the data concerning the quantity of supplied CO
2
is shown, and in Figure 4 the
quantity of supplied fly ash as a part of suspensions.
0
1
2
3
4
5
6
7
01-
may
21-
may
10-
jun
30-
jun
20-jul 09-
aug
29-
aug
18-
sep
08-
oct
28-
oct
17-
nov
07-
dec
27-
dec
16-
jan
05-
feb
25-
feb
17-
mar
C
O
2
a
p
p
l
i
c
a
t
i
o
n
i
n
t
o
t
h
e
g
o
b
[
t
]
Figure 3. Application of the carbon dioxide into the D-31 longwall gob
Copyright 2006 Taylor & Francis Group plc, London, UK
111
0
100
200
300
400
500
600
700
800
900
01-
may
21-
may
10-
jun
30-
jun
20-jul 09-
aug
29-
aug
18-
sep
08-
oct
28-
oct
17-
nov
07-
dec
27-
dec
16-
jan
05-
feb
25-
feb
17-
mar
F
l
y
a
s
h
a
p
p
l
i
c
a
t
i
o
n
(
i
n
s
u
s
p
e
n
s
i
o
n
s
)
i
n
t
o
t
h
e
g
o
b
[
t
]
Figure 4. Application of the fly ash, as a part of suspensions, into the D-31 longwall gob
There was a period between September and November 2005, when heating process in the gob
advanced, and the situation was serious. The evaluation of the use of the fly ash suspensions with
carbon dioxide as a fire prevention method focused on this period. The values of fire indicators
(CO and G) and the amount of carbon dioxide (transported with suspension) applicated into the
gob during this period, measured in the FS station, are presented in Figures 5 and 6. In the begin-
ning of September the value of the Graham indicator exceeded 40 [10
-4
]. After the supply of 9 ton-
nes of CO
2
and about 1800 tonnes of fly ashes (into the suspensions) within six days the Graham
indicator value was decreased to 20 [10
-4
]. The situation was stable, but for a few days only. Be-
tween 21st September and 6th October 24 tons of the carbon dioxide and 730 tons of fly ash were
supplied to the D-31 longwall gob. Thanks to that the decrease of the value of the fire indicators
could have been observed again. And again the situation was stable for a few days. This scenario
repeated twice, and the value of the Graham indicator reached more than 120 [10
-4
], and after a few
days of application of the fly ash suspension with CO
2
(2224 Mg CO
2
was supplied) was decrea-
sing to 2030 [10
-4
] level. The situation was under full control before the end of November 2004.
0
200
400
600
800
1000
1200
1400
1600
1800
3
0
-
0
8
9
-
0
9
1
9
-
0
9
2
9
-
0
9
9
-
1
0
1
9
-
1
0
2
9
-
1
0
8
-
1
1
1
8
-
1
1
2
8
-
1
1
8
-
1
2
1
8
-
1
2
2
8
-
1
2
7
-
0
1
1
7
-
0
1
Date
C
O
c
o
n
c
e
n
t
r
a
t
i
o
n
[
p
p
m
]
.
0
1000
2000
3000
4000
5000
6000
7000
8000
C
O
2
a
p
p
l
i
c
a
t
i
o
n
[
k
g
]
CO concentration
CO2 application
Figure 5. CO concentration measured in FS station and the amount of CO
2
supplied to the D-31 longwall gob, in the described period
Copyright 2006 Taylor & Francis Group plc, London, UK
112
0
50
100
150
200
250
300
350
400
3
0
-
0
8
9
-
0
9
1
9
-
0
9
2
9
-
0
9
9
-
1
0
1
9
-
1
0
2
9
-
1
0
8
-
1
1
1
8
-
1
1
2
8
-
1
1
8
-
1
2
1
8
-
1
2
2
8
-
1
2
7
-
0
1
1
7
-
0
1
Date
G
r
a
h
a
m
i
n
d
i
c
a
t
o
r
G
[
1
0
-
4
]
0
1000
2000
3000
4000
5000
6000
7000
8000
C
O
2
a
p
p
l
i
c
a
t
i
o
n
[
k
g
]
Graham indicator G
CO2 application
Figure 6. The value of the Graham indicator (computed on the base of measurements in FS station)
and the amount of CO
2
supplied to the D-31 longwall gob, in the described period
At this time the second part of the tests began. It was a single application of carbon dioxide
with maximum efficiency of the installation. It caused a decrease of value of Graham indicator.
This action and later regular of the fly ash application with carbon dioxide supplying helped to
stabilise the situation (the level of the G value no more than 40 [10
-4
]) until April 2005. When the
value of the Graham indicator reached 3040 [10
-4
] it was a signal to start suspensions with CO
2
start application.
6. INFLUENCE OF CARBON DIOXIDE APPLICATION UPON THE MINE ATMOSPHERE
During the test in the D-31 longwall, the permanent control of the mine air composition, in front of
and behind the head, was conducted. In few cases an increasing value of CO
2
was noted, especially
in FS and GL stations, when heating processes in the gob went on. The concentration of carbon di-
oxide in the outlet airflow never reached the permissible level (1%).
-6 000
-5 000
-4 000
-3 000
-2 000
-1 000
0
1 000
2 000
3 000
4 000
5 000
6 000
7 000
8 000
9 000
10 000
11 000
0 2 4 6 8 10 12 14 16 18 20 22 24 26
Hour of the tests [h]
C
O
2
a
p
p
l
i
a
c
t
i
o
n
[
k
g
]
0
0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1
1,1
1,2
1,3
1,4
1,5
1,6
1,7
C
O
2
c
o
n
c
e
n
t
r
a
t
i
o
n
i
n
o
u
t
l
e
t
a
i
r
f
l
o
w
[
%
]
D-31a - CO2 application
D-31 - CO2 application
D-31a - CO2 concentration in outlet airflow
D-31 - CO2 concentration in outlet airflow
Figure 7. The quantities of carbon dioxide applicated into the gob during the tests (upper graphs, left scale)
and CO2 concentration measured in the outlet airflow OU station (lower graphs, right scale)
Copyright 2006 Taylor & Francis Group plc, London, UK
113
Special tests were carried out in order to precisely evaluate the CO
2
application to the gob and
on the carbon dioxide concentration in the outlet airflow. The first of them, mentioned above, was
realized on 21st November 2004. About 6.9 tonnes of CO
2
was added into the gob within 21 hours
by the pipeline in the liquidated part of the lower gateway D-31. It is significant that the direction
of CO
2
applications was in accordance with the airflow direction through the gob (from the lower
gateway D-31 to the upper gateway D-31a). The concentration of the CO
2
was measured in the OU
station in the upper gateway D-31a. The results of the tests are presented in Figure 7. The increase
of the CO
2
concentration was not noted.
The second test was conducted in April 2005, when 10.2 tonnes of CO
2
was added to the gob
within 20 hours from the pipeline in the upper gateway D-31a. The pipeline in the liquidated part
of the upper gateway was 55 m long. In that case the direction of the CO
2
application was opposite
to the direction of the airflow through the gob. The measurements taken in the OU station showed
an increase in the values of CO
2
concentration in comparison with the results of the previous test.
The permissible level of CO
2
1% concentration was not reached.
7. SUMMARY
The example of a new method of the fire prevention, the application of fly ash suspension with
carbon dioxide during regular longwall exploitation in an underground coalmine was described.
Several times higher values of fire indicators (CO concentration and the Graham indicator) were
observed. It pointed at spontaneous heating in the gob. The new method allowed for retarding and
controlling all the processes without disturbing the exploitation.
During the research the limitation of heating process was observed (abating of heating symp-
toms) especially after the periods of intensified application of fly ash suspension with carbon di-
oxide. The best results were obtained when stable amounts of fly ash suspension with carbon dioxi-
de were regularly applicated. Longer breaks, even several days, caused reappearance of spontane-
ous heating process.
A good example of this is the action undertaken from September 2004 to November 2004. Af-
ter the application of fly ash suspension with carbon dioxide within this period, there was a long
period without spontaneous heating symptoms.
The conducted research indicates that it is possible to avoid exceedance of permissible level of
CO
2
concentration in the airflow even if a maximum dose of CO
2
was applicated. It may be stated
that, according to assumptions, CO
2
migrates into deeper parts of the gob. At the same time the in-
fluence of the place of CO
2
application on its flow in the gob was observed.
REFERENCES
Budryk W.: Grnictwo tom XI: Poary i wybuchy w kopalniach. Wyd. Grniczo-Hutnicze, Stalinogrd 1956.
Mazurkiewicz M., Piotrowski Z., Tajdu A.: Lokowanie odpadw w kopalniach podziemnych. Biblioteka SEP,
Krakw 1997.
Mazurkiewicz M., Piotrowski Z., Pomykaa R.: Zawiesina popioowo-wodna jako rodek transportu CO
2
do
zrobw kopalni podziemnej. Materiay Konferencyjne SGO, Krakw 2004.
Mazurkiewicz M., Piotrowski Z., Pomykaa R.: Zawiesina popioowo-wodna z dodatkiem CO
2
jako rodek w pro-
filaktyce przeciwpoarowej. Materiay Konferencyjne SEP, Szczyrk 2005.
Pomykaa R.: Waciwoci zawiesin popioowo-wodnych z ditlenkiem wgla stosowanych w profilaktyce po-
arowej, Praca doktorska, niepublikowana, Krakw 2006.
Strumiski A.: Zwalczanie poarw w kopalniach gbinowych. Wyd. lsk, Katowice 1996.
Tor A. et al.: Dowiadczenia w inertyzacji z zastosowaniem dwutlenku wgla do zwalczania zagroenia po-
arowego w kopalniach Jastrzbskiej Spki Wglowej S.A. Materiay 3 SAG, Zakopane 2004.
Copyright 2006 Taylor & Francis Group plc, London, UK
115
A New Method for Assessing the Status
of Sealed-off Coal Mine Fires
Yuan Shujie
Anhui University of Science and Technology. Huainan, Anhui, China
Nikodem Szlazak, Dariusz Obracaj
AGH University of Science and Technology. Cracow, Poland
Si Chunfeng
Wanbei Coal Electricity Company. Suzhou, Anhui, China
ABSTRACT: The paper presents a method to assess the status of sealed-off underground coal mine
fires. The assessment is based on the analysis of gas samples taken from the fire areas. Considering
the mining an geological factors influencing the results, the authors suggested their own method of
analysing and assessing the condition of a sealed-off fire with the use of so-called fire character-
ristic, which graphically describes the variations of five fire gas components oxygen, nitrogen,
carbon dioxide, carbon monoxide and hydrocarbons in time. An analysis of the fire gas compo-
nents tendencies and their correlations allows drawing proper conclusions about the fire status. In
order to mark fire indicators tendencies in time the authors applied the time series analysis. The
cases studied confirmed that the suggested method is beneficial for analysing and assessing the sta-
tes of fires located in sealed-off areas.
KEYWORDS: Coal mines, mine fires, spontaneous combustion, sealed-off fire areas, fire status
assessment, fire indicators
1. INTRODUCTION
Underground fires are one of the serious hazards in coal mines. Fire fighting is often done by seal-
ing off the area containing the fire. After the fire area is sealed, the necessity arises to assess the
fire status, in order to decide whether to open the sealed area. Heterogeneity of the medium makes
it difficult to determinate the physical phenomena occurring in the sealed space. The methods of
fire status assessment depend mainly on gas samples taken from the sealed area and analyses of
their composition. The various methods of fire status assessment are based on fire indicators,
which are calculated from the concentrations of fire gases in the samples. Known fire indicators are
not always correctly used and not always describe the status of the fire contained in a sealed area
properly. The changes of fire indicators values in time have more significance in the proper inter-
pretation of fire state. Up to now no reliable, practical indicators signalising extinction of a sealed-
off fire has been worked out. A fire may exist in a sealed area from several weeks to several de-
cades. Therefore, proper interpretation of fire status parameters is very important in the safe open-
ing of sealed areas as soon as possible and recovering the locked up coal reserves. Although a ran-
ge scientific investigations were done, there are no unambiguous criteria for assessing the fire sta-
tus and the problem is still not solved.
Copyright 2006 Taylor & Francis Group plc, London, UK
116
2. FACTORS INFLUENCING SEALED-OFF FIRE STATUS ASSESSMENT
Fires in coal mines can be divided into endogenous and exogenous, according to their genesis. En-
dogenous fires are caused by spontaneous combustion, and exogenous fires by open flame. In the
case of tight sealing of fire areas both of the two kinds of fires go through the same process leading
to their extinction. In the process of fire extinction, if oxygen is present in a sealed area and tem-
perature higher than normal atmospheric temperature, coal spontaneous heating or combustion may
occur.
The necessary condition for spontaneous coal combustion to occur is concurrent presence of the
following [Struminski 1996, Maciejasz 1977]:
Fragmentised coal prone to low-temperature oxidization.
Appropriate airflow.
Heat accumulation during coal oxidization.
The first of the conditions is connected with the natural character of coal, and the other two de-
pend on mining-technological conditions, namely the mining system, ventilation method and so on.
Except for its natural tendency to spontaneous combustion, the following physical characteristics
of coal also affect its spontaneous combustion [Cygankiewicz 2000, Struminski 1996]:
Degree of metamorphism.
Petrography composition.
Moisture content.
Initial temperature.
Fragmentation.
Thickness of coalbed.
Coalbeds dip angle.
For certain coal types the main factor influencing spontaneous combustion is physical parame-
ters of the air flowing through the fragmentised coal pile, for example flow speed, temperature and
humidity, which are affected by mining-geological conditions: exploitation system, ventilation sys-
tem and strength parameters of rock. A substantial number of theories were formulated on sponta-
neous combustion of coal [Maciejasz 1977, Struminski 1996, Strzeminski 1984].
Determination of sealed fire extinction criteria is very important for safe opening of sealed are-
as and resuming ventilation there. At the moment of opening of sealed areas spontaneous heating
or combustion often re-starts, and sometimes explosion of fire gas occurs. Fire extinction depends
on oxygen concentration in the atmosphere and coal temperature. The two parameters dont change
identically with time. Minimum concentration of oxygen supporting flaming combustion is shown
in Tables 2.1 and 2.2.
Table 2.1. Minimum concentration of oxygen supporting flaming combustion
of some combustible gases [Greuer 1973]
Other non-combustible gases in mixture
Combustible gas
Air and nitrogen Air and carbon dioxide
CH
4
12,1 14,6
H
2
4,9 5,9
CO 5,6 5,9
Table 2.2. Minimum concentration of oxygen supporting flaming combustion of some fuels [Mitcher 1996]
Fuels Timber Low rank coal Middle rank coal High rank coal Methane
Ignition 2 5 10 18 12
Flaming combustion 2 5 5 5 12
Smouldering 0 0 0 0 -
Copyright 2006 Taylor & Francis Group plc, London, UK
117
The minimum concentration of oxygen supporting combustion is far smaller for solid fuels than
for combustible gases. Research shows that the minimum concentration of oxygen supporting com-
bustion of coal and timber is 2% [Mason 1993, Stark 1965]. In general, its accepted that combu-
stion of solid fuels stops when concentration of oxygen falls below 2%.
The main reason of wrong decisions to open sealed fire areas is incorrect analysis and assess-
ment of fire status.
Wrong assessment of fire status may be caused by following factors:
The seal doesnt completely stop air from flowing through the fire area, the air carries oxygen
to the fire source and allows spontaneous heating to progress.
Because of carbon monoxide adsorption by coke it cant be unambiguously said that fire is ex-
tinguished only because the concentration of carbon monoxide decreased to zero.
Pressure equalization on stoppings disables airflow through them. In sealed areas small amount
of carbon monoxide may remain for a long time. It may be produced by low-temperature oxi-
dation of other organic substances, not only coal combustion.
In case of untight stoppings the concentration of oxygen may increase with the concentrations
of CO, CH
4
, H
2
, C
x
H
y
and CO
2
decreasing. If the sealed area is situated near a big fault and clo-
se to the surface, the concentration of oxygen in the atmosphere in the sealed area may stay high.
Methane inflow from strata or nitrogen supply to the fire source may change the composition of
fire gases and influence proper assessment of fire status on the basis of gas samples taken from
the sealed area.
Because of unknown temperature distribution in a sealed area its difficult to estimate the fire
source temperature. If the temperature remains high, smouldering may occur even at very small
concentration of oxygen, which may cause fire to re-start when the sealed area is opened.
The above shows that we cant properly assess fire status on the basis of a single indicator nei-
ther by its spot value, nor its changes in time. To avoid mining and geological factors influencing
the results, fire status assessment can be done by analysing the fire gas composition, the correlation
between its component elements and the time trends.
3. METHOD OF PROCESSING DATA FROM SEALED-OFF FIRE AREAS
Spontaneous heating or combustion in sealed-off areas is a slow, sequential process of coal oxi-
dation. Fire extinction is also a sequential process. Sudden changes of gas composition cant be re-
garded as a basis for assessment of a sealed-off fire status. We can analyse fire state by studying
changes of gas concentration in some period of time by taking a series of samples. In order to avoid
drawing incorrect conclusions regarding the changes of fire gas composition caused by mining-
geological reasons or inaccurate analyses of gas samples, three period moving average method is
used to analyse the tendency of fire gas composition in time.
As a result of serial (at certain time intervals) surveying of concentration of gases behind stop-
pings we can get some number of measurement values. When taking a large number of measure-
ments certain statistical dispersion of results always appears. The results can be regarded as ran-
dom variables. The changes of gas concentration may result in difficulties in interpreting the phe-
nomena occurring in the sealed area. As an example broken line in Figure 3.1 shows oxygen con-
centration changes in a sealed area. On its basis its difficult to estimate tendency of gas concentra-
tion, decreasing or increasing, in time. In this case, time series analysis can be applied to describe
the course of gas concentration changes.
Time series analysis can be applied in the case of unknown functional relationship between
random variables y and time intervals t [Brandt 1999]. Time series can be defined as series of ran-
dom variables with connected distribution {y
i
; t = 1 n}. Results of gas composition measurement
and fire indicators calculated on their basis are the variables. Random variables are the sum of true
Copyright 2006 Taylor & Francis Group plc, London, UK
118
measurement values
i
and measurement errors
i
. If random variable function is subject to certain
trend, it can be described by a curve passing the measurement points. In order to obtain possibly
smooth function of variable t, every value can be written as:
k i
k i j
i i
y
1 k 2
1
u (3.1)
Its an average (unweighted) value of measurements taken in moments t
i-k
, t
i-k+1
t
i
, t
i+1
t
i+k
,
and its named as moving average of variable y
i
. Variable k presents the time interval, on which
values of variable y
i
are averaged. When k = 1 the formula (3.1) yields three period averages, when
k = 2 five period averages, and etc. For the examined case of oxygen concentration shown in
Figure 3.1, their three period averages are calculated and presented in Figure 3.1 by solid line. On
the basis of changes of gas components concentrations and fire indicators values by three period
average, we can easily interpret their trends in time. In the following analysis the three period ave-
rage method will be used.
0
5
10
15
20
25
0 5 10 15 20 25 30 35 40 45 50 55 60
Time (day)
O
x
y
g
e
n
c
o
n
c
e
n
t
r
a
t
i
o
n
(
%
)
Figure 3.1. Oxygen concentration changes in a sealed area.
------ trend on basis of measured value of oxygen concentration,
trend on basis of three period average value of oxygen concentration
4. APPLICATION OF FIRE CHARACTERISTIC TO ANALYSE SEALED-OFF FIRE STATE
4.1. Idea of fire characteristic
A general analysis of fire tendency (development or extinction) can be supported by fire charac-
terristic, which in the same figure graphically describes the concentration tendencies of the five
components of fire gas oxygen, nitrogen, carbon dioxide, carbon monoxide and hydrocarbons
in time. Development trends of particular components of fire gas in an examined period can show
useful information.
Change tendency of oxygen concentration in sealed areas. If coal oxidation in sealed areas oc-
curs, oxygen concentration decreases. The rate of oxygen concentration decrease depends on oxi-
dation intensity. The bigger the oxidation intensity, the greater the rate of oxygen concentration de-
Copyright 2006 Taylor & Francis Group plc, London, UK
119
crease. In case of spontaneous heating or combustion of coal, oxygen concentration fluently de-
creases, and its course can be described by exponential function in dependence on the time from
the moment of sealing the fire area. When explosion of gas or coal dust occurred, oxygen concen-
tration in the sealed area intensively decreases. If spontaneous heating, combustion or explosion doe-
snt occur, course of oxygen concentration changes is monotonous. If fire is subject to extinction,
oxygen concentration in the atmosphere increases or stabilizes around a certain value. When gases,
such as methane from strata, inflow to the sealed area, oxygen concentration decreases.
Change tendency of carbon monoxide concentration. Carbon monoxide is the main fire indica-
tor gas. The rate of carbon monoxide concentration increase in a sealed area is considerably greater
than that of other gaseous products of coal oxidation. Carbon monoxide concentration increase in-
dicates combustion or low-temperature oxidation of coal. As opposed to carbon dioxide, the origin
of carbon monoxide is only combustion or low-temperature oxidation of coal and other organic
substances. The lower limit of carbon monoxide concentration in the atmosphere below which it
does not burn in fire is 12.5%. It is almost not changed by other reasons except fire. Generally after
fire extinction, carbon monoxide concentration in atmosphere slowly gradually decreases. In some
cases, carbon monoxide concentration in sealed areas may remain high for a long time after its
extinction, which may be the reason for erroneous assessment of the fire status.
Change tendency of carbon dioxide concentration. Carbon dioxide is the main gas produced du-
ring coal combustion. Carbon dioxide production is higher than that of other gaseous combustion
products. During a fire carbon dioxide concentration increases with oxygen concentration decrease
at the same or higher rate. Inversely, if carbon dioxide concentration in the sealed area atmosphere
decreases, oxygen concentration will increase. If oxygen concentration decreases slowly, the inten-
sity of fire development is not big. The fact of small rate of carbon dioxide concentration decrease
cant be the basis to state that the intensity of fire development is small. In coal mines carbon
dioxide may come from strata or be produced in a reaction of carbonates with acid water. Carbon
dioxide can be adsorbed by coke or soot produced in the fire and be dissolved in water.
Change tendency of nitrogen concentration. If the process of spontaneous heating in a sealed
area doesnt occur, the tendency of nitrogen concentration change is similar to that of oxygen.
When the process of spontaneous heating of coal occurs, the tendency of nitrogen concentration
changes is contrary to that of oxygen-nitrogen concentration increases and oxygen concentration
decreases. When open fire exists, nitrogen concentration remains stable or slowly decreases, while
oxygen concentration decreases fast. When explosion occurred, nitrogen concentration remains
around its normal concentration in the air, and oxygen concentration intensively decreases.
Change tendencies of hydrogen and hydrocarbons concentrations. Hydrogen (H
2
), ethane (C
2
H
4
),
propylene (C
3
H
6
) and sometimes acetylene (C
2
H
2
) or butylene (C
4
H
8
) release in certain sequence
with the increase of coal temperature. Detection of hydrogen, ethylene, propylene and acetylene in
gas samples in certain sequence is the sign of spontaneous heating development, while during
reduction in spontaneous heating the gases disappear in reverse sequence.
Change tendencies of fire indicators in time supply better information about the processes ta-
king place in sealed areas than their spot values.
On the basis of results from general analysis of fire development according to fire characterris-
tic, we can farther more analyse the course of complex fire indicators changes.
4.2. An example of assessing a sealed-off fire status with the use of fire characteristic
In X-Q coal mine longwall 1123(3) area was sealed off from the ventilation system by upper and
lower stoppings because of spontaneous heating of coal. After sealing of the fire area gaseous ni-
trogen was injected into the sealed area from the side of the lower stopping to create inert at-
mosphere. Gas samples were taken from behind the stoppings. On the basis of gas composition
analyses fire characteristic were obtained (Figures 4.1 and 4.2).
Copyright 2006 Taylor & Francis Group plc, London, UK
120
0
20
40
60
80
100
1
9
9
9
-
0
1
-
1
2
1
9
9
9
-
0
1
-
2
2
1
9
9
9
-
0
2
-
0
1
1
9
9
9
-
0
2
-
1
1
1
9
9
9
-
0
2
-
2
1
1
9
9
9
-
0
3
-
0
3
1
9
9
9
-
0
3
-
1
3
1
9
9
9
-
0
3
-
2
3
1
9
9
9
-
0
4
-
0
2
1
9
9
9
-
0
4
-
1
2
1
9
9
9
-
0
4
-
2
2
1
9
9
9
-
0
5
-
0
2
1
9
9
9
-
0
5
-
1
2
Date
O
2
,
N
2
(
%
)
0
2
4
6
8
10
12
14
C
H
4
,
C
O
2
(
%
)
,
C
O
(
p
p
m
)
CH4 CO2 CO N2 O2
Figure 4.1. Fire characteristic obtained on the basis of gas samples taken from behind the upper stopping
0
20
40
60
80
100
120
1
9
9
9
-
0
1
-
1
2
1
9
9
9
-
0
1
-
2
2
1
9
9
9
-
0
2
-
0
1
1
9
9
9
-
0
2
-
1
1
1
9
9
9
-
0
2
-
2
1
1
9
9
9
-
0
3
-
0
3
1
9
9
9
-
0
3
-
1
3
1
9
9
9
-
0
3
-
2
3
1
9
9
9
-
0
4
-
0
2
1
9
9
9
-
0
4
-
1
2
1
9
9
9
-
0
4
-
2
2
1
9
9
9
-
0
5
-
0
2
1
9
9
9
-
0
5
-
1
2
Date
O
2
,
N
2
(
%
)
0
2
4
6
8
10
12
C
H
4
,
C
O
2
(
%
)
,
C
O
(
p
p
m
)
CO CH4 CO2 N2 O2
Figure 4.2. Fire characteristic obtained on the basis of gas samples taken from behind the lower stopping
General analysis of fire characteristics on the basis of gas samples from behind both of the
stoppings shows that the change tendencies of carbon monoxide concentrations behind the two
stoppings are identical. After sealing off the longwall 1121(3) carbon monoxide concentration de-
creased fast and on January 21, 1999 reached zero. The maximum value of carbon monoxide con-
centration was only 13 ppm, which probably is the result of dilution by nitrogen injected into the
sealed area. In this situation we cant estimate the status of spontaneous heating or combustion of
coal only on basis of carbon monoxide concentration value. A comparison of carbon monoxide
concentration trend with the trends exhibited by the other components of fire gas can be helpful in
fire status estimation. In Figure 4.2 can be seen that nitrogen concentration is about 90% and at the
end of the period it presents increasing trend. Meanwhile oxygen concentration is about 10%, and
then decreases to about 7%. Carbon dioxide concentration is below 1%, and methane concentration
below 2%. Methane and carbon dioxide concentrations have decreasing trends similar to that of
oxygen, which indicates that their decrease is a result of dilution by nitrogen injected into the sea-
led area, and oxygen decrease is not due to oxidation of coal. After January 21, carbon monoxide
Copyright 2006 Taylor & Francis Group plc, London, UK
121
disappeared, which can indicate, that from this day oxidation process stopped. Before January 21,
1999 even though carbon monoxide appeared oxygen concentration remained at a level of about
10% and didnt show decreasing trend. The course of oxygen, nitrogen, methane and carbon di-
oxide concentrations changes had similar (parallel) trends. Therefore it can be supposed, that after
the moment of sealing the area to January 21, 1999 spontaneous heating process didnt develop.
Nitrogen was injected from the side of the lower stopping, so it can be supposed, that methane and
carbon dioxide moved in the direction of the upper stopping. So behind the upper stopping with
a nitrogen concentration increase, methane and carbon dioxide concentrations increased, which
wasnt a result of spontaneous heating or combustion of coal. Oxygen concentration decrease cor-
responds to nitrogen concentration increase. Oxygen concentration decrease mainly results from
nitrogen concentration increase. Even though before January 21, 1999 carbon monoxide appeared,
no other hydrocarbons appeared, which can be acknowledged as low-temperature spontaneous hea-
ting of coal.
In order to confirm the general estimation of the assessed fire status some complex fire in-
dicators were analysed. Trends of complex fire indicators in time are presented in Figures 4.34.6.
In the case of nitrogen injection into the sealed area, fire indicators basing on nitrogen concen-
tration and oxygen consumption cant be used as a base to interpret the fire status. In this situation
only CO/CO
2
index can be used to assess the fire state. The maximum value of CO/CO
2
index
behind the lower stopping didnt attain 0.002, and behind the upper stopping it even was below
0.0005 and had decreasing trend in the monitored period. Too small a CO/CO
2
index value indi-
cates that spontaneous combustion didnt occur, only low-temperature heating occurred at the be-
ginning of the monitored period, which confirms the results of the analysis done on the basis of the
fire characteristics.
Figure 4.3. Trend of CO/CO
2
index
behind the upper stopping
Figure 4.4. Trend of CO/CO
2
index
behind the lower stopping
Figure 4.5. Trend of Graham index
behind the upper stopping
Figure 4.6. Trend of Graham index
behind the lower stopping
0,000
0,001
0,002
1 1 1 1 1 1 1 1 1 1 1 1 1
Date
C
O
/
C
O
2
0,000
0,001
0,002
0,003
0,004
0,005
1 1 1 1 1 1 1 1 1 1 1 1 1
Date
C
O
/
C
O
2
0,000
0,003
0,006
0,009
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
1
9
9
9
-
.
.
.
Date
1
0
0
G
0,000
0,003
0,006
0,009
1 1 1 1 1 1 1 1 1 1 1 1 1
Date
1
0
0
G
Copyright 2006 Taylor & Francis Group plc, London, UK
122
CONCLUSION
The suggested in the paper method of sealed-off fire status analysis is based on an analysis of fire
indicators time trends and their correlation. The analysis is conducted in two steps. In the first step
the fire characteristic is used, which allows drawing preliminary conclusions about the coal spon-
taneous combustion process in progress, moreover allows to eliminate the factors influencing the
results of fire state assessment. The final confirmation the fire status is done by an analysis of com-
pound fire indicators tendencies in time. The suggested method can appropriately indicate the sta-
tus of a fire in a sealed-off area, which has the vital significance in practice for taking decisions
about opening or abandoning the sealed-off fire areas. The cases studied confirmed that the sugges-
ted method is beneficial for analysing and assessing the states of fires located in sealed-off areas.
REFERENCES
Brandt S.: Analiza danych. Wydawnictwo Naukowe PWN, Warszawa 1999.
Cygankiewicz J.: Niektre wskazniki pozarowe w kopalniach wgla. Biblioteka Szkoly Eksploatacji Podziem-
nej, Seria z Perlikiem Nr 2, Krakw 2000, ISBN 83-87854-72-7, pp. 45.
Greuer R.E.: Influence of Mine Fires on the Ventilation of Underground. Contract, No. S0122095, United Sta-
tes Bureau of Mines, 1973, pp. 173.
Maciejasz Z., Kruk F.: Pozary podziemne w kopalniach, Czsc 1, Wydawnictwo Slask, Katowice 1977, pp. 11,
22, 2428, 23.
Mason T.N. i inni: Gob Fires, Part 2, The Revival of Heatings by Inleakage of Air. Paper of Safety in Mines
Research Board, No. 76, 1993.
Mitcher D.W.: Mine Fires, Prevention, Detection, Fighting. 3rd ed., Maclean Huster, USA, 1996, pp. 167.
Stark G.W.V.: Use of Jet Engines for Fires Fighting, Tests at the Fire Research Station. Institution of Fire En-
gineers, Quarterly, No. 60, 1965, pp. 345360.
Struminski A.: Zwalczanie pozarw w kopalniach glbinowych. Wydawnictwo Slask Sp z o.o., Katowice
1996, ISBN 83-7164-023-4.
Strzeminski J.: Wczesne wykrycie pozarw endogenicznych. Przeglad Grniczy, Nr 5, 1984.
Yuan Shujie: Selection and Estimation of Sealed-off Fire State Indication Methods. Doctorial Dissertation,
AGH University of Science and Technology, Cracow, Poland.
Copyright 2006 Taylor & Francis Group plc, London, UK
123
Development of Strategy for Revival
of Regional Coal-Mining Companies
A.A. Malysheva
Industrial-Ecological Laboratory, The Centre for Monitoring the Socio-Ecological Consequences
of Closing the Eastern Donbass Coal Mines Ltd
The established approach to the development of market relations in the TEC (technical and econo-
mic complex) of Russia resulted in essential change of fuel and energy balance of the country where
in 2004 the share of petroleum and gas made up more than 80%, whereas the share of coal is
reduced to a critically dangerous level 1214%. At the same time certain mining-geological and
investment problems are aggravating and the tendency of the rapidly exhausting natural gas and
petroleum prices to rise at the home market is steady. This will inevitably result in the necessity of
their partial replacement with coal.
The contemporary coal industry in Russia is very business-orientated, so its main attention of is
concentrated on the operation of highly profitable deposits yielding high investment return and lo-
cated far from the consumers. Besides, carried out without sufficient scientifically-grounded esti-
mation of the available potential, re-structuring of coal industry resulting in closure of more than
200 collieries, essentially aggravated regional social, economic and environmental problems, prob-
lems of employment, especially in the areas where coal enterprises were of major importance for
the citys economy, that is were city-forming. It should be noted that small coal mining companies
commercially unprofitable for objective reasons were much more likely to be doomed for closure
than big ones.
The reforms of economic relations in the fuel and energy sector of Russia and the restructuring
carried out in the recent years resulted in a number of negative consequences to the functioning of
the coal industry: a significant reduction of coal share in the fuel and energy balance (12% in
2003); a decrease in coal output (6065% of the level of 1990); lower asset amortization rates (up
to 15%); full elimination of small deposits ready for extraction from the stock exchange (due to
their objectively low adaptation to the market); under utilization and in some cases complete demo-
lition of the available basic production assets; retrenchment of qualified workers and scientific staff
that aggravated social and environmental problems in the regions traditionally focused on coal
mining.
A transition to the use of natural gas for power generation and temporarily low internal prices
of this commodity connected with an underestimation of the investment profit forming factors,
exhaustibility and rarity of its deposits, have considerably lowered the competitiveness of the coal
mining industry in Russia as a whole.
Not only individual enterprises but also whole regional coal mining complexes, those represen-
ting small and average deposits of local importance in particular, have undergone closure. While
the average decrease in coal extraction for the period of re-structuring was less than 50%, some
coal-mining regions reduced their coal production ten-fold and more.
Now the production share of coal enterprises working small deposits of local importance makes
up less than 1% of the total Russian coal output. Meanwhile, an analysis of the foreign experience
shows that countries with available coal resources increase coals share in their fuel and energy
Copyright 2006 Taylor & Francis Group plc, London, UK
124
balances (USA, China, Canada), and small deposit operations are encouraged by the state, as for
the most part working of small coal deposits promotes development of small independent power
engineering, does not do essential harm to the environment and allows to solve problems of unem-
ployment. On average, small coal mining enterprises abroad are responsible for 8 to 15% of the
total production.
An intention to increase the coal output and bring the absolute volume of extraction up to 350
million t. by 2020 begins to take shape now, making it necessary to work out a new strategy of coal
business development. However, just as it was during the period of re-structuring, insufficient at-
tention of the economists and administrative bodies is devoted today to developing coal-mining
companies of local importance, serving the regional market, capable of solving a number of fuel
and energy, social, ecological and economic problems by effective utilization of mineral resources
available at their territory.
Various methods and aspects of increasing the efficiency of operating coal enterprises in the
market conditions have been devised in sufficient detail. Among publications on this subjects the
works of A.S. Astahova, N.A. Arhipova, L.A. Belashova, M.M. Gurena, B.A. Davydova, V.T. Ko-
valja, S.S. Lihtermana, N.J. Lobanova, J.N. Malysheva, A.A. Petrosova, V.P. Ponomareva, V.A. Har-
chenko, A.B. Janovskogo, M.A. Jastrebinskogo and other scientists are widely known.
An analysis of these works showed that they often disregard modern specificity of development
of coal-mining base of local importance and their regional potential, and underestimate the oppor-
tunities of their effective utilization for increase of power safety and social and economic stability
in the coal mining regions.
From these facts transpires the existence of an actual scientific problem how to perfect the avail-
able and develop new methodical approaches to ecological and economic maintenance and effect-
tive running of regional coal-mining companies, working small and average-sized coal deposits, to
promote an integrated use of coal and development of diversification, and the need to solve it. The
urgency and the importance of this research for the coal sector and the economy of coal-mining
regions will grow in line with the growth of prices of energy resources.
Development of a methodical approach to the maintenance of ecological and economic effi-
ciency of running coal companies working small and average-sized deposits of local importance,
promoting growth of extraction, processing and consumption of locally available fuel and energy
resources, increase in output of commodity from a unit of extracted and processed resources and
employment growth is one of the primary goals of the modern mining economic science.
To attain this aim it is necessary to solve the following problems:
To analyse the role of local coal resources, associating minerals and waste products of coal pro-
duction in the development of a regional economy.
To investigate the condition of regional coal basins after re-structuring and to reveal the basic
ecological, social, economic and energy supply regional problems.
To systematize the internal (regional) and external factors promoting an increase of ecological
and economic efficiency of functioning of coal companies.
To establish the criteria of ecological and economic estimation of the overall performance of
coal-mining companies exploiting small and average-sized deposits and using the strategy of
diversifications by reclaiming associated useful minerals, waste products of mining and the te-
chnogenic deposits formed as a result of partial closure of coal enterprises.
To work out an algorithm allowing to reveal an ecologically and economically effective way of
developing coal companies of local importance promoting growth of coal output with consi-
deration for diversification of production and provision of social living standards of mining
workers.
To develop administrative measures for stimulating the development of small coal businesses
and reveal ways to maintain their effective work.
Copyright 2006 Taylor & Francis Group plc, London, UK
125
The methodological approach to overcome the above-stated problems utilizes the scientifically
proved method of managing the efficiency of working small coal deposits of local importance by
increasing competitive commodity output from coal and coal wastes having market value, on the
basis of their complex ecological and economic estimation and making use of the specificity of the
resource potential of coal-mining regions.
The existing technology of coal production is objectively connected with production of great
volumes of solid waste products. Technogenic deposits in places of closed coal enterprises and
solid waste products from coal production occupy extensive areas and negatively affect the con-
dition of water and air. The enterprises of the coal sector spend significant means on transportation
and disposal of solid waste products, pay environmental faxes for their disposal and the pollution
of environment arising from the harmful emissions in places where waste dumps are located.
Transition to chiefly opencast methods of coal mining and paying insufficient attention to an
integrated development of coal deposits result in further increase of specific parameters of forma-
tion and accumulation of solid waste products per 1 t. of coal, despite the reduction of volumes of
produced solid waste products, caused by closing coal-mining enterprises: in 1990 the yield of
waste products amounted to 4.2 t, and in 2004 it rose to 7.6 t.
Building industry
Possible ways of practical use of solid technogenic coal
Ceramics, paint and varnish industry Agriculture
Metallurgy The chemical industry Secondary fuel
1. Silting and filling material,
filling excavated mine workings
(sand, polymictic, aleurolites,
sandstones)
2. Grout, concrete (sand,
polymictic, sandstones,
aleurolites, argillites, clay
slates, ash, ceramsit,
agloporit)
3. Rock for building and
filling, crush stone for various
purposes (sandstones,
aleurolites, carbonaceous
rock, ash)
4. Construction of
artificial filtering areas
(aleurolites and sandy
slates)
5. Manufacture of fire-resistant,
ceramic, silicate brick (porcelain
clay, hydromicaceous rock,
argillites, clay slates,
carbonaceous rock)
6. Manufacture of
glass wool of
GOST 4640/81*
(aleurolites,
argillites)
1. Conditioning and deoxidizing of
soil, enhancing its biological
properties, granulating fertilizers
(glauconitic sandy rock, aleurolites and
sandy slates, clay and
sandy-argillaceous rock, carbonaceous
rock)
2. Enhancing agricultural
chemical properties of
recultivated soil layer
(ash)
3. Original nitric and
phosphoric fertilizers,
additives for manufacture
of complex fertilizers
(waste products of
concentration)
1. Manufacture silicate
alloys, carborundum,
silicon (quarzitic
sandstones)
2. Iron melting (sideritic
concretions)
3. Producing of metal by
a biotechnological method
(ferriferous-sandy
concretions)
1. Manufacture of pottery
(glauconitic rock, porcelain clay,
argillites, clay slates)
2. Pigment for manufacture of
building paints on various bases
(ash)
1. Sorbents, bleaches (porcelain
clay, clay and sandy-argillaceous
rock, siliceous rock)
1. Power and household fuel -
after secondary concentration
and bricketing (waste products of
concentration)
Figure 1. Possible ways of practical use of solid technogenic coal wastes
While the volume of waste products accumulated by the industry is considerable, a mere 3% of
them are used as a commodity. Meanwhile, technogenic deposits of coal wastes contain a great
deal of useful elements which get an increasing industrial value and social, ecological and eco-
nomical importance due to the exhaustion and limitation of primary natural resources, lack of fi-
nancial means for nature protection work and a growing problem of unemployment among miners
discharged in the result of closure of coal enterprises. Possible ways of practical use of solid tech-
nogenic coal wastes are presented in figure 1.
Copyright 2006 Taylor & Francis Group plc, London, UK
126
Thus, the integrated approach to the utilization of small coal deposits of local importance, often
developed and ready for extraction, will allow to increase the effectiveness and financial viability
of their mining using the existing production assets, qualified manpower and infrastructure and, in
conditions shortage of investments in the coal sector, solve the increasing problems of power sup-
ply in regions.
Effective and competitive operation of small coal deposits should be fulfilled on a basis of di-
versification and integrated development of coal resources with the use of market and administra-
tive measures and keeping in mind the need to establish long-term economic relations with the end
consumers.
Copyright 2006 Taylor & Francis Group plc, London, UK
127
Mining of Adjacent Coal Seams in the Mines of Western Donbass
Volodymyr Buzylo, Tatiana Morozova, Vitaliy Vasilev, Olexandr Koshka
National Mining University. Dnepropetrovsk, Ukraine
Its well known that mining of adjacent coal seams in the mines of Western Donbass became
difficult due to special mining and geological conditions. These stem from the fact of the adjacent
coal seams being separated by rock middling between 5 and 20 m, strong and tough coal, having
the resistance for cutting of up to 300 kN/m, combined with low mechanical strength of the sur-
rounding rock, rock tendency for separation and swelling, high vulnerability of rock to moisture
resulting in 50% up to complete loss of its carrying capacity in humid conditions, disjunctive and
frequent disturbances resulting in fractured and weak rock zones.
Mining and geological conditions in all the working mines of Western Donbass are the fol-
lowing: the sequence of preparation is level by level, mining system is the longwall method (from
700 to 1600 m) by coupled or single faces with individual lengths of 150200 m, mining of seams
is carried out in descending order and it takes 22,5 years to mine two mining blocks on neigh-
bouring seams with longwalls working in stepped configuration. The mining blocks are developed
by raises, all faces are equipped with mechanized longwall sets.
All analytical analyses and calculations of stress and strain conditions of the rock mass avail-
able in literature confirm that carrying out stoping in adjacent seams (with rock thickness between
the seams of 520 m) results in the increase in the intensity of parameters and rock stress in the
overmined and undermined zones. Underground investigations were carried out with the aim to
prove the truth of the basic theoretical calculations and to obtain actual values of the parameters by
which rock stress manifests itself in longwall faces in the situation when adjacent seams are mined.
All investigations were carried out in Stepnaya mine OAS MCK Pavlogradugol in 2001
2004 in coupled faces 213 and 215, located in the
6
/
seam and in coupled faces 117 and 119 in
the C
6
seam. The distance between the seams is representative to the area that is the rock thickness
between the seams within the perimeters of the examined production sites varied from 8 to 12 m.
Figure 1 shows the mutual positions of the longwall faces in question. The length of the mining
block is 11001200 m, the length of longwall faces 213, 215 and 119 is 200 m and the length of
longwall face 117 is 175 m.
All longwall faces are equipped with KMK 97M mechanized sets with MK-98 powered support
and MK 67M cutter loader.
The main goal of the on-site investigations was to establish functional characteristics and de-
pendence of rock stress manifestation parameters in column and in zones of mutual influence of
production and development workings.
Investigations were done according to special methodology, taking into account the basic points
and demands of universally accepted branch methodology. Actual values of rock joint parameters
in the production faces were established. The character of the roof rock and coal seam caving was
examined.
While observing stability conditions and rockwall displacement near the face area great atten-
tion was paid to fissures, faults, separation of roof rock and cave-ins, the moments of appearing of
roof disturbances were established; the distances from the face where they happened were recorded.
Copyright 2006 Taylor & Francis Group plc, London, UK
128
The crack width and rock displacements were measured. The character of deformation and rock ca-
ving (rock size, thickness of rock fall, caving sequence, and the influence of roof disturbances on
the character of roof caving near the face area) was considered. In each case the angle of intersec-
tion between the direction of the main jointing of roof rock at the production face and the cleavage
present in the coal seam were established.
Figure 1. Relative positioning of the examined longwall faces
The results of these analyses showed that the character and regularity of the rock stress mani-
festing itself in all the examined faces were practically the same. It should be pointed out that in
some cases there was up to 0,15 m of coal yield from the seam in faces 213 and 215 and no such
phenomenon was observed in the seam C
6
. Average collapse volume in the faces in seam
6
/
were
1520% smaller than in seam C
6
. Roof rock collapses most often took place in the middle of the
faces and at their perimeters.
The fracture formation process in all longwall faces was the same, that is, after passing of the
cutter loader vertical cracks appeared, which were opened up having the width of 1020 m after
extracting the next coal seam. After the third cycle the rock of the immediate roof was broken up
by chaotically arranged cracks into large pieces of irregular shape, but nevertheless the contour and
integrity of block between the technological cracks remained the same. The immediate roof rock
caved into the worked-out area after each movement of a powered support section. In some cases
a 1,63,2 m long cantilever of hanging rock formed above the worked-out area in the faces located
in the upper seam.
The parameters of the powered support were recorded with the help of pressure gauge and self-
writing manometer MP66A, connected to head ends of powered support hydraulic props. While
carrying out investigations concerning support power parameters the pressure of the medium in the
hydraulic system was changed in the range of 1822 MPa being 20 MPa at average. The pressure
in the front props was 515% higher than in the back ones, which is typical for mining and geo-
logical conditions of Western Donbass. As a rule, the support hydraulic props where operating in
the condition of increasing resistance. Safety valves of the hydraulic props were setting off in prac-
tically all the sections. In the longwall faces in the upper seam the working resistance of the sup-
port hydraulic props was in average 510% higher, than in faces 117 and 119.
The results of the experiments showed that the actual extraction ratio in all the examined faces
was on average 0,95. At the bottom end of the face rock with the thickness of 0,250,35 m was
Seam
6
Seam
6
/
6
/
6
200 m
170200 m
200 m
170200 m
200 m 175 m
longwall 215
longwall 213
8
1
2
m
longwall 117 longwall 119
w
Pokj
R. Polska
Polska
Wirek
Halemba
lsk
Katowice
Kleofas
R
.
K
a
to
w
ic
e
Wujek
Siemianowice
Saturn
Grodziec
Pary
Sosnowiec
Porbka
Klimontw
Kazimierz
Juliusz
Wieczorek
Mysowice
Niwka
Modrzejw
Jan
Kanty
Staszic
ZGE
Sobieski Jaworzno
Wesoa Murcki
M300
Z
a
b
r
z
e
B
ie
ls
z
o
w
ic
e
Makoszowy
Sonica
Bolesaw
miay
Budryk
Knurw
Szczygowice
Dbiesko
Pstrowski
Gliwice
Siersza
Ziemowit
Piast
Janina
Czeczott
Brzeszcze
Silesia
Rydutowy
R.
Rymer
Chwaowice
Jankowice
Marcel
Anna
ory
Krupiski
Borynia
R. 1 Maja
Jas-Mos
Pniwek
R. Mos
Zofiwka
Morcinek
Thorez
Victoria
Wabrzych
Wabrzych
Szczawno Zdrj
Gliwice
Tychy Rybnik
Chrzanw
Spytkowice
Owicim
Czechowice
Bielsko-Biaa
Jaworzno
Czech Republic
Czech Republic
POLAND
mines finishing production
in 2002
operating mines
closed mines
ZGE
Sp.z o.o.
ZG
Piekary
w
Pokj
R. Polska
Polska
Wirek
Halemba
lsk
Katowice
Kleofas
R
.
K
a
t
o
w
i
c
e
Wujek
Siemianowice
Saturn
Grodziec
Pary
Sosnowiec
Porbka
Klimontw
Kazimierz
Juliusz
Wieczorek
Mysowice
Niwka
Modrzejw
Jan
Kanty
Staszic
ZGE
Sobieski Jaworzno
Wesoa Murcki
M 300
Z
a
b
r
z
e
B
i
e
l
s
z
o
w
i
c
e
Makoszowy
Sonica
Bolesaw
miay
Budryk
Knurw
Szczygowice
Dbiesko
Pstrowski
Gliwice
Siersza
Ziemowit
Piast
Janina
Czeczott
Brzeszcze
Silesia
Rydutowy
R.
Rymer
Chwaowice
Jankowice
Marcel
Anna
ory
Krupiski
Borynia
R. 1 Maja
Jas-Mos
Pniwek
R. Mos
Zofiwka
Morcinek
Thorez
Victoria
Wabrzych
Wabrzych
Szczawno Zdrj
Gliwice
Tychy
Rybnik
Chrzanw
Spytkowice
Owicim
Czechowice
Bielsko-Biaa
Jaworzno
Czech Republic
Czech Republic
POLAND
mines finishing production
in 2002
operating mines
closed mines
ZGE
Sp.z o.o.
ZG
Piekary
Others
Figure 1. RIAM Conceptual Model Outline [Liming 2002]
The potential benefits stemming from the application of the integral AP-913 ERP are as follows:
Consistent, systematic and rigorous approach to major investment value-based decision-ma-
king.
Consistent treatment and integration of plant safety, reliability, efficiency, and cost factors in the
decision-making process.
Consistent framework for continuous monitoring, and projection of plant safety and production
performance.
Operating & Maintenance procedure improvement and procedure training prioritisation.
Copyright 2006 Taylor & Francis Group plc, London, UK
147
Trade-offs between on-line and off-line maintenance.
Capital spares procurement analysis and prioritisation.
Unit efficiency improvement, etc.
2.3. Preventive Maintenance Optimisation (PMO)
Following several references, a standard or classical RCM was not found to be an optimal process
with regard to the needs and expectations of the nuclear industry. It usually led to the creation of
a new preventive maintenance program rather than to enhance and revise the existing ones. Fur-
thermore, a considerable amount of time and effort required to perform an analysis and implement
PM tasks recommendations has been considered as another major weakness for standard RCM ap-
plications [Johnson 1998, Turner 2002, Harazim & Ferguson 2003, Messaoudi et al. 2003, Messa-
oudi 2005].
The experience has shown that approximately 60% of the time required to accomplishing a clas-
sical system RCM analysis is spent on the following two major steps:
1. Failure Modes and Effects Analysis (FMEA) identification of equipment critical for each sys-
tem function.
2. Logic Tree Analysis (LTA) identification of the most feasible and effective PM tasks to pre-
vent critical component failure modes and causes of concern.
This situation is due to the level of details required and excessive documentation produced as
a result of the rigid process steps. In this context, in mid-1990s the nuclear industry initiated rese-
arch projects to develop ways to perform RCM more cost effectively, based on its own experience.
Some concepts have been developed in order to demonstrate that the conventional RCM methodo-
logy [Moubray 1997] could be streamlined to reduce the cost of analysis while maintaining a high
quality product. More information about the results of these initiatives (works of EPRI, Fractal-so-
lutions, INPO, etc.) may be found, among other sources, in the work of Messaoudi (2005).
The PMO concept is based on a more streamlined RCM approach. With PMO, many nuclear
generating stations have achieved very high capacity factors (over 90%) at a lower cost, while en-
suring a safe production [Johnson 1998 and 2003, Turner 2002, Harazim & Ferguson 2003, Mes-
saoudi et al. 2003, Coppock 2004, Messaoudi 2005]. Table 1 (below) shows a comparison between
the classical RCM analysis and streamlined methods.
Table 1: Comparison of various RCM related concepts after [Messaoudi 2005]
Streamlined analyses
Classical RCM analysis phases
Streamlined
Classical Process
PMO
Process
Criticality
Checklist
Determine System Boundaries X X X
Determine Subsystems 6)
Data Collection and Plant History Review X X X
Identify system Functions
and Functional Failures
1) 1)
FMEA 2) 3) 4)
Instrument Matrix
7)
X
Non-critical Evaluation X X X
LTA PM Task Recommendations
(Critical/Non-critical)
5) 5) 5)
Task Comparison X X X
Copyright 2006 Taylor & Francis Group plc, London, UK
148
1) Identify all system functions, and classify into two major groups: (1) Important functions, and (2) Non-
Important functions.
2) Perform FMEA only on components and equipment that support important/critical functions.
3) Perform a streamlined FMEA, which combines dominant failure modes and plant effects into one compo-
nent record on components (equipment) that support important functions.
4) Component (Equipment) is critical if its failure will result in one of several pre-selected plant effects, which
one wants to prevent through a PM program.
5) Maintenance Templates added to formal LTA. LTA only formally documented in modified classical analysis.
6) Only for very large (complex) systems.
7) Replaces FMEA for Instrumentation & Control (I&C) equipment.
As shown in Table 1, the PMO employs many of the same analysis techniques as the classical
RCM, but it represents a more streamlined approach. The classical RCM analysis begins at the top
with a system, breaks it down into subsystems, identifies critical components, and recommends PM
tasks. The latter are subsequently compared to those already in place (existing PM tasks), which
leads to final task recommendations. Contrary to that, PMO begins at the opposite end. The PM
procedure is split into tasks that are reviewed to identify the failure they are intended to prevent.
Subsequently, related data are collected, critically reviewed and analysed in order to make final
task recommendation [Johnson 1998, Messaoudi 2005]. Another major difference in the PMO pro-
cess vs. the classical RMC analysis is that the standard FMEA and LTA have been streamlined and
combined into one record (see Table 1). The PMO process documents the most important failure
modes for a component/equipment (usually three dominant ones), and allows documenting the most
dominant plant effects (typically three) in one component/equipment record, used to determine com-
ponent/equipment degree of criticality. For those considered critical, PM tasks are recommended
using maintenance templates.
The most significant elements contributing to the successful completion of the system analysis
using this concept are as follows:
Combination of failure modes and plant effects into one record.
Combination of FMEA and LTA into one component record with less formal documentation.
Efficiency and savings related to performing LTA in the nuclear industry have been achieved
through a concept of Maintenance Template for a given component or equipment type. At the
stage of a template development, several component characteristics such as equipment frequency
of usage, operational and service environment as well as functional importance are taken into ac-
count. Each maintenance template recommends appropriate tasks for a given component or equip-
ment type in terms of condition monitoring tasks, time directed tasks, and surveillance tasks. They
provide major additional savings. In the nuclear industry there currently exist maintenance tem-
plates for nearly 80 major components and/or equipment considered critical with regard to both sa-
fety and power generation. They have been developed through a teamwork involving manufac-
turers, nuclear industry people, and EPRI experts who benefited from over 20 years of operating
experience, including data collection and analysis.
The experience gathered in the nuclear industry has shown that experienced analysts using the
PMO process have been able to achieve results similar to the modified classical RCM process, and
the concept itself has proven to be highly efficient in terms of achieved savings and overall results
[Johnson 1998 and 2003, Harazim & Ferguson 2003, Messaoudi 2005].
As an example of the PMO process advantages, some results concerning improvement of equip-
ment reliability and plant efficiency at Kewaunee Nuclear Power Plant (KNPP) will be presented
[Johnson 1998 and 2003]. The plant is situated close to the city of Green Bay (Wisconsin, USA).
The PM Optimisation was initiated following the need to shift the plants PM program from a 12-
to 18-month cycle. The technical basis (justification) of maintenance tasks and frequencies needed
to be re-established because the original records were not retrievable. Initial efforts to compile si-
milar information consisted of a conventional RCM program. In five years, RCM evaluations were
Copyright 2006 Taylor & Francis Group plc, London, UK
149
completed on only 8 of the plants 65 systems. Critical review of the RCM results and progress re-
vealed that although the process was rendering useful information, it was advancing too slowly gi-
ven KNPPs needs and expectations. In mid-1990s, the KNPP management decided to turn to PMO.
A plan, considering long-term plant operation, was defined with a following main objective:
Maximize plant generating capability and equipment reliability. The other objective was to es-
tablish a general plant maintenance program that would minimize power reductions and forced out-
ages resulting from equipment failures, procedural errors, technical inadequacies, and human error.
Also, specific objectives were formulated and PMO streamlined approach was used to achieve them.
One year after the PMO program had been initiated, the analysis of all plant systems was accomp-
lished. As a result, PM task change recommendations were drawn and, subsequently implemented
into existing KNPP programs [Johnson 1998].
An overall evaluation performed four years later [Johnson 2003] revealed that since the begin-
ning of KNPPs Living Program the following significant results have been achieved (only some of
them are listed here):
Savings resulting from PM interval extension and task elimination over a four year period rea-
ched 1.5 M$ (US).
Many PM tasks dropped and some new added.
Number of PM work orders dropped drastically due to grouping related tasks between plant or-
ganisations.
12% reduction in maintenance staffing.
Reliability Program recognized as strength by US Nuclear Regulatory Commission (USNRC)
and INPO.
Spin-off of other maintenance processes.
Acceptance of a new program was made through staff involvement; etc.
As shown above, the PMO process has proven its efficiency, and is now widely used in the nu-
clear industry.
In brief, its potential benefits are the following:
Elimination of unnecessary PM tasks.
Decrease of corrective maintenance in terms of time, man-labour and cost.
Reduction of forced outages.
Improvement of plant availability.
Focus of maintenance resources on critical equipment and components.
PMO balances corrective, preventive and predictive maintenance tasks along with equipment
availability and maintenance cost considerations. The ultimate goal consists in maximizing plant
value while remaining safe [Liming 2002]. RIAM provides a framework for employing those tools
and methods. It extends the scope from plant maintenance policies and practices to all aspects of
asset management. Rather than optimising on the basis of reliability, availability and cost, as is the
case of a conventional approach, RIAM does an optimisation on the net present value (NPV) as
a principal profitability parameter [Liming 2002].
Following encouraging results obtained through PMO applications in the nuclear industry, there
are the reasons to believe that its implementation could be rewarding in other industries, too. In this
context, the next section discusses some possible applications of the PMO process in mining.
3. POTENTIAL FOR A TRANSFER OF NUCLEAR INDUSTRY KNOW-HOW
AND PRACTICES TO MINING OPERATIONS
As mentioned before, conventional RCM is usually a time-consuming process. At the same time,
several applications of the ERP and PMO processes in the nuclear industry have delivered encou-
raging results in terms of saving time, effort and money while remaining safe.
Copyright 2006 Taylor & Francis Group plc, London, UK
150
One may expect that these processes may be transferable to the mining industry, particularly in
the following circumstances:
Mills and mineral treatment plants.
High-scale open-pit mines where capital and operational cost of the equipment as well as its
performance are significant contributors to the production cost.
Existing and future underground mines where production depends (or will depend) to a large
extent on highly automated production equipment.
In all cases, in the context of fierce international competition, the mines need the highly effect-
tive equipment. Concerning the case of mineral treatment plants, they operate a considerable num-
ber of stationary equipment, which makes them somehow similar to nuclear power plants. In this
respect, a transfer of know-how and experience could be probably the easiest. In the two remaining
cases there are some significant differences between them. Open-pit mines employ huge fleets of
capital-intensive equipment and their production systems have got a limited redundancy. Many of
them operate around the clock for 365 days a year often-in extreme atmospheric conditions. In the-
se circumstances, ensuring high equipment reliability and availability is a very challenging task
and a success depends to a large extent on adequate maintenance policies and strategies. Even if in
some cases classical RCM seems to deliver anticipated results, it is worth trying less time- and re-
source-intensive processes, such as those reviewed in preceding sections.
In what concerns the automated equipment in the underground mines, its operational and ser-
vice environment is very specific and substantially different from that of man-operated equipment.
Due to safety concerns and regulations, zones where automated equipment operates are usually off
the limits for human workers. In the case of a failure that requires human intervention, it is an
imperative that all man-less traffic and operations in the area concerned be suspended and/or shut
down. This may bring a whole section of the mine to a virtual halt and, thus, induce considerable
production losses. Also some minor (from a technical point of view) problems that would other-
wise be rapidly solved by a human operator may lead to production interruptions. Since in under-
ground mines there is usually even less redundancy (if any) in production systems than at open-
pits, automated drill jumbos and/or rigs, LHDs and trucks are really critical for a whole production
process.
In each of the above cases, the improvement of reliability and optimisation of proactive main-
tenance (preventive, predictive, condition based) is a key factor with regard to production efficien-
cy and cost. In this context there is apparently some potential for the application of the ERP and
PMO processes in the mining industry.
At the beginning, these processes need to be conceived and developed in the way to adapt them
to a specific character and context of mining. At the early stages, the work should focus on some
selected categories of equipment considered critical with regard to its role in the production system
(including a degree of redundancy), its working conditions and environment, intensity of their use,
maintainability and maintenance support. Cable (electric) shovels, hydraulic shovels and crushers
in open-pit mines, and underground automated jumbos, production drill rigs and LHD are good ex-
amples here. Databases and maintenance templates for those may be developed with an input from
OEM, mine operators and maintenance management people, consultants and research institutions.
They are founded reasons to believe that a successful development of specific ERP and PMO
and their implementation at mine sites will bring substantial benefits such as:
Less time and effort deployed, but with similar results compared to conventional RCM.
Substantial reduction of reactive maintenance, whose share at many mine sites is by far too
high.
Rationalization and optimisation of proactive maintenance tasks will definitely contribute to an
increase of availability and reduction of a maintenance cost.
Better utilization of maintenance resources due to focusing on the truly critical equipment and/
/or its components.
Copyright 2006 Taylor & Francis Group plc, London, UK
151
ERP may be considered almost a closed-loop system that once properly implemented will lay
ground for continuous reliability and PM improvement (see Fig. 2).
Increase of equipment effectiveness and better life-cycle management.
Figure 2. Top Level Diagram of AP-913 Equipment Reliability Process [INPO 2001]
It should not be forgotten however, that the development and implementation of these two pro-
cesses can never be done instantaneously and their success requires some time, effort, proper un-
derstanding, support and commitment as well as patience of all the involved parties (especially
from the management).
CONCLUSIONS
Facing tough international competition, mining companies in develop countries are constantly for-
ced to produce more at the lesser cost. One of the most important avenues in this respect is an in-
crease of overall equipment effectiveness. Adequate preventive and proactive maintenance policies
and tasks (particularly in what concerns their scope, content and frequency) are among the most
crucial factors. However, the mining industry has not yet developed its own original processes for
this purpose and it relies predominantly on the know-how and experience gathered in other, more
innovative industries. This paper intends to indicate that there is a potential for transferring some
important know-how from the nuclear power generation industry. This includes two processes, na-
mely: Equipment Reliability Process (ERP) and Preventive Maintenance Optimisation (PMO). Their
successful implementation requires that they be adequately adapted and tailored to fit within the
operational and business context of the mining industry. Another indispensable condition is a close
collaboration and open-minded exchange between OEM, mine operators and their maintenance
management personnel as well as other involved personnel. This approach will be undoubtedly be-
neficial for mining companies facing a tough competition at the international level.
Copyright 2006 Taylor & Francis Group plc, London, UK
152
REFERENCES
Campbell J.D. 1995: Uptime Strategies for Excellence in Maintenance Management. Productivity Press,
Portland, Oregon (USA).
Coppock S. 2004: Raising the Bar, Palo Verdes Secret for Good Equipment Reliability. Nuclear Professional,
Atlanta.
DOE-DP-STD-3023-98. 1998: Guidelines for Risk-Based Prioritisation of Doe Activities. (Doe Limited Stan-
dard), U.S. Department of Energy AREA MISC, Washington, DC 20585.
Ednie H. 2004: RCM Helps Operations Get the Most out of Equipment, CIM Bulletin, vol. 95, no 1065, pp.
2025.
Faitakis Y., Mackenzie C. and Powley G.J. 2004: Reducing Maintenance Cost through Predictive Fault De-
tection. CIM Bulletin, vol. 97, no 1076, pp. 6366.
Harazim M.L. & Ferguson J.B. 2003: North American Nuclear Maintenance Best Practices Compared to
Japanese Utility Maintenance Practices. 24th Annual Conference of the Canadian Nuclear Society, Toronto.
Institute of Nuclear Power Operations (INPO) 2001: AP-913, Revision 1, Equipment Reliability Process De-
scription. Atlanta.
Johnson L.P. 1998: Improving Equipment Reliability and Plant Efficiency through PM Optimisation. Fractal
Solutions, Suwanee.
Johnson L.P. 2003: Continuous Improvement with PM Optimisation. Lessons Learned, Fractal Solutions, Su-
wanee.
Komljenovic D., Vaillancourt R. & Messaoudi D. 2004: Approach Regarding the Implementation of S-98
Reliability Program at the Gentilly-2 Power Plant Considering AP-913 Equipment Reliability Process.
25th Annual Conference of the Canadian Nuclear Society, Toronto.
Knights P.F. 1999: Analysing Breakdowns. Mining Magazine, vol. 181, no 3 (September), pp. 165171.
Knights P. & Oyaneder P. 2005: Best-In-Class Maintenance Benchmarks in Chilean Open-Pit Mines, CIM
Bulletin, vol. 98, no 1088, p. 93 (full text available on the CIM webpage: www.cim.org).
Kumar U. 2003: Economic Aspects of Equipment Breakdowns in Mines and Control Measures. Lecture 5,
Division of Mining Equipment Engineering, Lulea University of Technology, Lulea, Sweden.
Kumar U. 2003: Reliability Technique a Powerful Tool for Mine Operators. Lecture 13; Application of Re-
liability Centred Maintenance (RCM) Strategy, Lecture 14, Division of Mining Equipment Engineering,
Lulea University of Technology, Lulea, Sweden.
Lewis M.W. 2000: Remote-Control Monitoring Cuts Maintenance Cost, Mining Engineering, vol. 52, no 8,
pp. 1519.
Liming K.J. 2002: Risk-Informed Asset Management for Electric Power Plants. Energy Pulse RIAM Paper.
Mercier E., Valle P., Tapp B., Mthot R. & Clouet J.M. 2004: Maintenance at Falconbridge Ltd.s Raglan
Mine. CIM Bulletin, vol. 97, no 1084, pp. 5358.
Messaoudi D, Abdul-Nour G., Komljenovic D. & Vaillancourt R. 2003: RCM in Nuclear Power Plants. Ca-
nadian Reliability and Maintainability Symposium, Ottawa.
Messaoudi D. 2005: Processus de fiabilit des quipements la centrale nuclaire Gentilly-2. MSc. Thesis,
Universit du Qubec Trois-Rivires.
Mitchell J.S. (editor-in-chief) 2002: Physical Asset Management Handbook. 3rd edition, Clarion Technical
Publishers.
Moubray J. 1997: Reliability-Centered Maintenance. 2nd edition, New York, Industrial Press Inc.
NUREG-0800. 2002: Chapter 19, Use of Probabilistic Risk Assessment in Plant-Specific, Risk-Informed De-
cision-making: General Guidance, Revision 1 of Standard Review Plan Chapter 19, U.S. Nuclear Regula-
tory Commission, Washington, DC 20555-001.
NUREG-1513. 2001: Integrated Safety Analysis Guidance Document, U.S. Nuclear Regulatory Commission,
Washington, DC 20555-001.
Paraszczak J., Vachon J. & Grammond L. 1997: Benefits of Studies on LHD Reliability and Availability for
Mines. In Strakos et al. (eds), Proc. 6th Int. Symposium on Mine Planning and Equipment Selection, Stra-
kos et al. (eds), Ostrava, Czech Republic, 36 September 1997, pp. 469475, A.A. Balkema.
Turner S. 2002: PM Optimization, Maintenance Analysis of the Future. OMCS Ph 0419397035, Webpage:
http://www.pmoptimisation.com.au
Copyright 2006 Taylor & Francis Group plc, London, UK
153
Problems of Underground Ore Mining at Great Depths
in Kryvyi Rih Basin
Yuriy Vilkul, Yuriy Kaplenko, Mykola Stupnik, Victor Sydorenko
Kryvyi Rih Technical University, Ukraine
ABSTRACT: The paper presents an analysis of mining methods applied at Kryvyi Rih Iron Ore
Complex. Considers the technical and economic indicators of sublevel caving. Offers directions for
further development of the technology of ore mining at great depths. Corresponding conclusions
are presented.
KEYWORDS: Underground ore mining, naturally rich ores, technical and economic indicators,
mining ratio useful mineral
Further development of underground iron ore mining technologies in Kryvyi Rih Basin the main
iron ore mining area in Ukraine must be directed first of all at increasing the quality of market-
able output, reducing production costs and improving the ecological situation in the region.
Mines that extract naturally rich ores carry out stoping at the depth of 10001250 m, develop-
ment at 12001300 m, opening new reserves at 13001500 m. Different mining and geological
conditions require application of different working methods and their variants. At the same time
a specific character of mining in Kryvyi Rih Basin great depths, high rock stress, relatively soft
unstable ores, fractured unstable surrounding rocks limit the range of methods possible to apply.
Mining and geological conditions do not allow applying highly efficient self-propelled machi-
nes and block-caving stoping methods. Most mines limit or even exclude the application of block-
and-room methods.
The summaries of mining methods applied at the rich ore mines of Open Joint Stock Company
Krivoy Rog Iron Ore Complex are given in Table 1.
Table 1. Mining methods applied at Kryvyi Rih Iron Ore Complex
Mines
Block-room
[%]
Sublevel-room
[%]
Sublevel-caving
[%]
Lenin 82.1 3.2 14.7
Gvardeyskaya 54.6 22.2 23.2
Oktyabrskaya 29.0 71.0
Rodina 97.8
Total in the Complex 45.0 24.9 29.7
As can be seen from Table 1, the most efficient block-room methods are responsible for less
than half of the output. The block-room mining methods provide high labour productivity and least
Copyright 2006 Taylor & Francis Group plc, London, UK
154
cost per 1 tonne of ore mined. However, under high rock stress, which often occurs in Kryvyi Rih
Basin mines their application, is being reduced. Nowadays in mines of the Basin sub-level room
systems are widely used though greater volumes of development workings characterize them. Thus,
for block-and-room mining A methods the specific quantity of development workings averages
4.56.5 m/1000 t, for sublevel-room methods it fluctuates within the limits of 6.08.0 m/1000 t, i.e.
1.31.5 times higher. Labour productivity differs correspondingly: 4060 t/cm for block-level
methods and 3545 t/cm for sublevel-room mining methods.
Sublevel caving is used to mine ores of medium and low hardness and stability under high rock
stress conditions. Among variants of the caving method applied in Kryvyi Rih Basin are caving
with sets of deep vertical holes into a free-face opening and onto pre-caved ore.
Technical and economic indicators of sublevel caving are given in Table 2.
Table 2. Technical and economic indicators of sublevel caving
Variants
Indicators
Caving into a vertical
free-face opening
(Kryvorizhstal Mine)
Caving into
a horizontal cut
(Sukhaya Balka Mine)
Caving onto
pre-caved ore
(Rodina Mine)
Panel area, m
2
500600 500600 500600
Total labour productivity, t/cm 3545 3040 2535
Ore losses, % 1517 1415 1415
Ore dilution, % 46 911 89
Quantity of development
workings, m/1000 t
6.26.5 6.57.5 7.58.5
Large specific quantity of development workings makes negative influence on technical and
economic indicators of ore mining. This results from the fact that the output per drift miner is 810
times smaller than that per driller employed at stoping. Besides, specific power consumption for
ore extraction at developing is 46 times higher than analogous index at stoping. Nowadays, the
cost of cartridge explosives widely used in mining is 34 times higher than granulated explosives.
Taking into account specific quantity of explosives and their cost, acquisition expenses for explo-
sives used in development are 1416 times higher than those for stoping. That is why the expenses
of preparing a block for stoping make up 3240% of all the operating costs of mining per block.
The process of formation of a vertical free-face space is technologically complicated and ex-
pensive. During the development of a raise the extraction of ore occurs in narrow rock exposure
conditions, which decreases hole efficiency parameters and increases the volume of drilling works
as well the quantity of explosives for breaking.
The cost involved in the formation of free-face spaces in the operating costs of the ore mining
per block varies from 10 to 20%. On the whole, the operating costs per block at Lenin and Ro-
dina mines vary from 0.8 to 2.00 USD per tonne correspondingly.
The extraction indicators depend to a great extent on the type of deposit and the elements of
blocks mined out during the current year.
In general, metal grade loss at the Complex improved by 0.33%.
Actual loss of ore at the Complex have amounted to 9.04 million tonnes over the past 11 years
with predicted loss of 9,8 million tonnes. This result was achieved by applying an ore concen-
tration principle enabling to produce ore from blocks with lower Fe content than is possible with
the standard method.
Copyright 2006 Taylor & Francis Group plc, London, UK
155
The extracted ore is a raw material for the production of marketable iron ore in preparation
plants. Kryvyi Rih deposits of hematite ores contain intrusions (up to 6% of volume) of rock with
low Fe content. Besides, in the process of mining, dilution of ore with surrounding waste rock ta-
kes place. The ores have hardness coefficient of 26 and rarely 810, whereas waste intrusions ha-
ve 1012 and more. Thats why the ore is crushed comparatively easily with a cascade of crushers
whereas substandard hard rocks can be crushed by means of repeated return to additional crushing
in a closed circuit.
The simple technology of ore concentration is based on this specific character of ores. It con-
sists in dividing the ore in the mines into two streams: pure ore and diluted ore. Pure ore is extrac-
ted and crushed completely in a closed circuit. Diluted ore is taken out in a separate stream. After
crushing according to a closed circuit, the product of screening is marketable sintering ore and over-
size fraction +10 mm with Fe content of 4445% makes the tailings of the concentration process.
Such a technology has enabled to increase the content of iron in sintering ore by 24% de-
pending on the volume of intrusions and the degree of dilution. All this has increased competiti-
veness of the ore. At the same time, the ores of the Complex are still inferior to the ores of other
suppliers at the international market.
The operating costs of ore mining depend on the mining and geological factors (thickness of
deposit, hardness and water content of ores, quantity of mining, drilling and support work) as well
as subjective factors: meeting the planned production volumes, labour management and efficient
use of resources. Thus the high operating costs of ore mining at Rodina mine are explained by
high water hazard (water inflow equals to those in the three remaining mines combined) and by
high quantity of development workings and support.
An analysis of expenditure indicates that 5 expenditure items: auxiliary materials; power; wa-
ges; maintenance of basic assets; mine general purpose expenditure form 77% of manufacturing
cost. So, they are to be the main subject of in-depth study for taking up appropriate steps for cost
saving, which is one of the main components in the struggle for the competitiveness of production.
As it was mentioned above the mines of the Complex use the simple concentration principle
based on a contrast of strength characteristics of ore and waste rocks. This enables to get quality-
sintering ore with increased Fe content due to crushing and sifting. It is necessary to take into con-
sideration the following fact. During the process of ore mining, the dilution rocks with an average
Fe content of 36.74% mix with it and the tailings with an average Fe content of 44.75% are dum-
ped during the process of ore concentration. Preliminary calculations show that 2499,4 thousand
tons of Fe were mined along with diluted rocks and 4208,6 thousand tons of Fe were lost along
with the tailings. This testifies to low efficiency of the concentrating process and calls for urgency
in its upgrading. It would be expedient to concentrate ore after each crushing stage. It would enable
to increase the grinding efficiency and to decrease power consumption, as power consumption du-
ring crushing is 1015 times higher than power consumption used for screening.
The applied technological methods of producing sintering ore focus on the raw material mined
directly in stope blocks.
For this reason, the major problem of underground mining technologies is the increasing of
qualitative and quantitative characteristics of extraction. Over the past 11 years the mines of the
Complex have lost 9042,8 thousand tonnes of ore with average Fe grade of 55.75% and at the same
time 6803,4 thousand tonnes of intrusive rocks with 36.74% Fe content were mined. The annual
tonnage of Complex averaged 5957.3 thousand tonnes. An average annual loss comes to 822.1
thousand tonnes. Only 1% in loss reduction enables to market commodity for the sum of U$76.5
thousand and to reduce costs for mining of intrusive rocks with low Fe content.
With regard to the level of 5% losses, the indicators of ore mining and concentration are cal-
culated with fixed crude ore productivity at the Complex (Fig. 1). The given results show that re-
ducing the losses from 13.8% to 5% results in increasing pure ore output and increasing Fe content
in it, increasing sintering ore output, decreasing the volume of diluting rock.
Copyright 2006 Taylor & Francis Group plc, London, UK
156
The share of individual processes in the mine and site operating costs are given in Figures 2 and 3
correspondingly. Using these data we can conclude that considerable improvement of technical and
economic indicators in mines is possible by upgrading the mining methods. It is determined by the
fact that the upgrading of underground mining technologies doesnt require capital investments,
and there are some reserves in other expense items of operating costs (for example, supplementary
materials). Specific quantity of expenditures for supplementary materials in mine operating costs
varies from 20.7% (Rodina mine) to 28.4% (Lenin mine). An aver-age for the Complex is equ-
al to 25.3%.
Development works is the main consumer of mine timber and metal support. Upgrading of mi-
ning methods in order to reduce the expenditure on development works will enable reducing the
expenditure for the abovementioned supplementary materials. Expenditures on explosives in the
Complex amounted to U$1378.2 thousand in 1999, for cartridge explosives U$655.05 thousand.
Reducing specific development works cost will allow to purchase expensive cartridge explosives.
As the development is the most laborious technological process of underground mining, reducing
the volume of development work will enable to increase labour productivity in general.
55,5
56,1
56,9
56,6
54
54,5
55
55,5
56
56,5
57
57,5
58
58,5
5 6 7 8 9 10 11 12 13
ore loss [%]
F
e
r
r
u
m
c
o
n
t
e
n
t
[
F
e
%
]
Ferrum content in sintering ore
Ferrum content in crude ore
4,7
5,0
5,5
5,3
5,9
0,4
0,6
0,9
1,2
0,0
1,0
2,0
3,0
4,0
5,0
6,0
7,0
5 6 7 8 9 10 11 12 13
ore loss [%]
O
r
e
m
i
n
n
i
n
g
[
m
l
n
t
o
n
n
e
s
]
Crude ore
Throw
Sintering ore
Figure 1. Change of work indicators of the Complexs mines
with the reduction of ore loss
Copyright 2006 Taylor & Francis Group plc, London, UK
157
breaking
16,7%
output & haulage
9,2%
crusher
6,6%
transport
5,3%
hoisting
9,0%
pumping
7,6%
ventilation
6,4%
amortization
4,6%
others
19,3%
development works
15,2%
Figure 2. Share of processes in mine operating cost
output & haulage
22,4%
breaking
40,6%
development works
37,0%
Figure 3. Share of processes in operating cost per block
The characteristic feature of underground mining is the complicated structure of mining me-
thods and a large volume of development workings. This results in large amount of manual work to
support the workings and large consumption of support materials. Consumption of mine timber is
2.14 m
3
per 1,000 tonnes of mined ore, consumption of metal support is 0.585 tonnes per 1,000
tonnes. Foreign enterprises use similar methods, but their distinctive feature is simplicity of method
and a sharp decrease in the amount of mining and support works.
The above examples indicate that considerable possibilities for improvement of technical and
economic indicators of the underground mining can be achieved by upgrading the mining methods.
The key element in increasing the efficiency of extracting minerals is upgrading the actual enginee-
ring methods and development of new systems and techniques of underground mining. Along with
development works, one of the basic and the most laborious processes of underground mining is
ore breaking, which is mostly (7480%) performed by wells. The specific quantity of labour input
in drilling and blasting is 40.6%.
The weak point in the technology of stoping is the process of formation of a free-face breaking
space. With a space amount of 2530%, labour inputs in working the free-face space are equal to
labour input in extracting the remaining 70% of the block.
So, the necessity of upgrading the drilling and blasting technique of stoping works is vivid.
Moreover, it is relevant not only to the technique of formation of a free-face, but also to the pre-
paration of deep wells, as its share in mining cost is 16.7%.
Copyright 2006 Taylor & Francis Group plc, London, UK
158
Blasting is a multi-factor process, in which the influence of various factors on the final result
isnt the same.
To establish the factors influencing the drilling and blasting efficiency, 9 years work results of
Kryvbass mines were processed. Analysis was done with the random balance method.
The following characteristic (further factors) were processed: the line of least resistance (W),
strength of broken ore (f), the specific quantity of explosives consumption for breaking (q), the
direction of breaking in respects to layering (), the mining depth (H), the number of drill-hole
rows (N), the volume of broken ore (Q), the type of explosives (Z), the area of free-face rooms (S),
and the orientation of the free-face space ().
Estimation of breaking efficiency was done by analysing specific consumption of explosives at
secondary breaking.
Diameter of the blastholes was not considered, as it is the same for all the Kryvyi Rih mines.
The values of linear and paired factors affecting the quality of ore breaking are given in Table 3.
The studies showed that the depth of pit development, the line of least resistance, physical and
mechanical properties of the rock mass, specific consumption of explosives and the value of rock
mass exposure are most important, i.e. exert most influence on the quality breaking of rock mass.
Thus not only physical and mechanical properties of rock, parameters of blasting and drilling but
also mining conditions where rock mass breaking and crushing occur influence the quality of ex-
traction.
Nowadays the share of ore haulage in mine manufacturing cost is 9,2% and is determined by
rock size distribution in the blasted rock and loading-hauling equipment type.
Two types of equipment: scraper and vibro-haulage installations are used at the mines of Kry-
vyi Rih. That is why the real factor providing the increase of haulage efficiency is improving the
granulation of the blasted rock mass.
The cost of transforming crude ore into sinter ore is 6.9% of the production cost. The above
mentioned data testify of low efficiency of crude ore processing resulting in considerable losses of
useful component to dumps, considerable content of 0.25 mm fractions in sinter ore and in increa-
sed energy consumption for crushing. This forces the necessity to perfect crude ore processing.
Table 3. Value and grade of basic line paired factors
Grade Factors Value, %
1 16,6
2 f 15,1
3 W 14,5
4 HW 13,4
5 f 11,6
6 Hq 10,7
7 HS 10,2
8 fq 9,4
9 q 8,9
10 Wf 8,8
11 7,4
12 fS 7,2
13 S 7,0
14 Wq 6,8
15 HN 6,5
16 WS 6,3
17 qS 5,9
Copyright 2006 Taylor & Francis Group plc, London, UK
159
Analysis of technology and technical and economical indicators of underground ore mining in
Kryvyi Rih Basin and abroad as well as of the state of iron ore product market in Europe allows to
come to the following conclusions:
Increase of iron ore product quality and technical and economical indicators of ore mining in
market conditions is the most important problem of national economy which must be solved in
the Complex by perfecting the technological processes of underground ore mining and the ore
processing technology at crushing and sorting plants.
Ore mining cost depends on the extraction ratio, as there exist several cost items whose share in
the total cost value is the function of extraction ratio. Prospecting, opening up and deposit deve-
lopment, block cutting and ore extraction costs are these items. That is why considerable cost
decrease can be reached by increasing the quality and quantity indicators of recovery and im-
provement of the named technological processes.
Useful mineral indicators are determined by the choice of development method, which speci-
fies the expenditure on block development and extraction. So, improving development method
with the aim to decrease specific volume of cut workings and the expenditure on creating free-
face breaking spaces together with the increase in the quality of recovery and quantity indica-
tors will ensure decreasing the mining cost.
Development works are the most labour consuming process of underground mining. This is
caused by high energy expenditures on rock extraction in conditions of narrow working space,
by high cost of manual labour required to support the workings. Thus, for example, expenditure
on explosives in development is 1416 times higher than in production operations taking into
consideration energy consumption and the cost of explosives. Besides, the specific expenditures
on timber and metal sets are today respectively 214 m
3
and 585 kg per 1000 t of mined ore.
That is why the improvement of production indicators can be reached not only by cutting the
production expenses but also by improving the drivage methods.
The bottleneck in the production operation technologies is the process of formation of the free-
face breaking space by widening of a cut raise. It requires increased number of shot-holes and
higher consumption of explosives for extraction. The expenditures on free-face space formation
per block cost change from 10 to 20%.
The share of expenditure on production of sintering ore in the total mining cost is 6.9% on
average at a mine. It is caused by imperfection of the crude ore processing technology in that
ore is run through a closed system of crushers leading to over-crushing of its richest part, the
increase of energy consumption and high iron losses in the tailings.
Copyright 2006 Taylor & Francis Group plc, London, UK
161
Satellite Remote Sensing as Means to Assess Primary Impact
of Open-Pit Mining Development
R. Latifovic
Natural Resources Canada/Earth Sciences Sector. Ottawa, ON, Canada
K. Fytas & J. Paraszczak
Dept. of Mining, Metallurgical and Materials Engineering, Universit Laval.
Qubec City, PQ, Canada
ABSTRACT: Assessment of the possible impact of existing and future mining projects located in
remote areas may be difficult and costly. This may have an adverse effect on the process of issuing
mine permits and development of new mines. A remote sensing based land cover change assess-
ment methodology presented in this paper is a promising avenue to render this process more mea-
ningful, precise and affordable. It has been applied to a case study of the Oil Sands Mining De-
velopment in the north-east of Canadian province of Alberta. So-called Athabasca Oil Sands Re-
gion (AOSR) is a place of unprecedented growth of open-pit mining operations. This study was
focused on the primary impact onto the landscape and vegetation, considered as all surface distur-
bances resulting from exploration, mining development, urban development and logging. This im-
pact was assessed using an information extraction method applied to two LANDSAT scenes. The
analysis based on derived land cover maps shows a decrease of natural vegetation in the study area
(715,094 ha) for 2001 approximately 8.64% relative to 1992. It has been demonstrated that the
methodology presented here provides reliable results and is fully applicable for the assessment of
the environmental impact of huge-scale mining operations.
1. INTRODUCTION
The mining industry does not have a good track record in what concerns its impact onto environ-
ment. Despite the impressive progress in reducing and mitigating its negative effects it is still per-
ceived by many as a harmful, big-scale polluter, leaving behind deep scars in the nature that are
difficult and long to heal. This is definitely one of the reasons that the expansion of existing ope-
rations or the attempts to open new ones usually encounter considerable reluctance (if not fierce
resistance) of the general public and the authorities involved. Understanding of all kinds of inter-
actions between mining activities and the environment is a key factor if a proper balance between
the benefits of resource exploitation and its possible adverse effects is sought. Therefore, there is
an urgent need to develop and implement tools and methods that would enable the objective assess-
ment of the effects exerted by open-pit mining activities onto surrounding ecosystems.
A successful monitoring approach for evaluating surface processes and their dynamics at a re-
gional scale requires observations with frequent temporal coverage over a long period of time in
order to differentiate natural changes from those associated with human activities. However, in so-
me countries such as Canada, mining activities move gradually to remote areas, for which long-
Copyright 2006 Taylor & Francis Group plc, London, UK
162
term field observations are usually unavailable and/or taking such measurements prior to mine de-
velopment may not be justified from the economical point of view. In most cases when historical
records are needed for studying the long-term vegetation cycles, satellite based remote sensing re-
presents virtually the only alternative to time-consuming and expensive data collection on the
ground. This modern technology has already gained widespread recognition as a powerful tool in
the studies over different aspects of ecosystems at local, regional and global scales. In addition to
being the only available data source in many areas, remote sensing has the added advantage of ac-
quiring data with sufficient area coverage and temporal frequency for studying and monitoring
primary impacts caused by surface mining at low cost. It can also be used for studying atmospheric
emissions and water pollution indirectly by monitoring green vegetation, an indicator of ecosystem
health and conditions. A number of published papers [Ress & Williams 1997, Schmidt & Glaesser
1998, Prakash & Gupta 1998, Wright & Stow 1999] suggest usefulness of such techniques for de-
tecting contamination, determining success in reclaiming open cast mined areas and for providing
other relevant spatial data for assessing mining impact on the environment.
Despite its immense potential with regard to environmental monitoring, it seems that satellite
based remote sensing has not yet been adequately used in the mining sector. In this context, the
principal objective of the study presented here was twofold:
to provide an initial assessment approach based on remote sensing data, and
to evaluate the use of such data in investigation of trends and intensity concerning land cover
changes in the Athabaska Oil Sands Region (AOSR).
This paper will focus on primary impacts of mining activities onto the landscape and vegeta-
tion, considered as all surface disturbances resulting from exploration, mining development, urban
development and logging.
2. MATERIALS AND METHODS
2.1. Study area
As the world crude oil reserves are steadily getting depleted and the demand for it is on the rise,
there has been a growing interest in mining oil-bearing tar sands, whose deposits are abundant in
north-eastern part of the Canadas province of Alberta. So-called Athabasca Oil Sands (AOS) de-
posit, situated north of the city Fort McMurray is one of the worlds largest known oil deposits.
According to estimates based on exploration results, there are at least 300 billions barrels of re-
coverable bitumen within 400 km of Fort McMurray (Fig. 1). Two other large oil deposits found
also in northeast Alberta are Cold Lake and Pearce River. The AOS deposit has become the site of
unprecedented development, with investments in excess of 20 billion US dollars. With several mi-
nes already in operation, many more being developed or awaiting statutory approvals, concerns are
often raised about the long-term impact they may exert on the surrounding area. The scale, intensi-
ty and future perspective of the AOS development have made the region a unique area for studying
different aspects of environmental impact. The satellite remote sensing technique offers an additio-
nal potential that can be combined with conventional approaches for analyzing spatial and temporal
extent of affected areas. Accordingly, this particular region has been chosen as a study area for
evaluating a remote sensing based approach for assessment of primary impacts on the surrounding
vegetation (Fig. 1).
The study area is situated in the Boreal Plains ecozone, and includes fractions of three ecore-
gions according to the Canadian Ecological Land Classification System (1995): Wabasca Lowland
(ecoregion number 142) in the southwest, Slave River Lowland (ecoregion number 136) in the
north and Boreal Uplands (ecoregion number 147) in the west. Cool summers and long cold win-
ters characterize the subhumid mid-boreal ecoclimate type with mean annual temperature of 0.5C
and mean annual precipitation ranging from 350 to 500 mm. Dominant vegetation types are me-
Copyright 2006 Taylor & Francis Group plc, London, UK
163
dium to closed canopy stands of aspen (Populus tremuloides) and balsam poplar (Populus balsa-
mifera) with white spruce (Picea glauca), black spruce (Picea mariana) and balsam fir (Abies balsa-
mea). Cold poorly drained fens and bogs are covered with tamarack (Larix laricina) and black
spruce. Organic soil is dominant covering 50% of the region. Land uses include forestry, water
oriented recreation, wildlife hunting and trapping, as well as oil and gas exploration. Apart from
mining and processing of oil sands, other industrial activities are present in the region that may also
contribute to stress on the surrounding environment. These include timber harvesting, urban gro-
wth of the municipalities of Fort McMurray and Fort McKay, as well as new linear disturbances
such as exploration seismic lines, roads and pipelines.
Figure 1. Location of the Athabasca Oil Sands Region (AOSR) study area in the province of Alberta, Canada.
The white boundary outlines area of recoverable oil sands reserves
2.2. Experimental design
In order to relate surface disturbance and trends in vegetation growth to the proximity of surface
mining operations, the study area was spatially delineated into three separate impact zones using
potential acid input (PAI) prediction. PAI is the preferred method for evaluating overall effects of
acid forming chemicals on the environment since it accounts for the acidifying effect of sulphur
and nitrogen species, as well as the neutralizing effect of available base cations. The PAI prediction
used in this study is adopted from the Suncor In Situ Oil Sands application. Annual concentration
and deposition were predicted using CALPUFF, a multi-layer, multi species non-steady state puff
dispersion model which can simulate the effects of time and space varying methodological condi-
tion on pollutant transport, transformation and removal (EPA 1995). Three impact zones were de-
fined using isopleths maps representing annual PAI and generic critical load values taken from the
Clean Air Strategic Alliance (CASA): low deposition (0.25 keq/ha/year), medium deposition (0.5
keq/ha/year) and high deposition (1.00 keq/ha/year) (Fig. 2). In the following text these areas are
referred as Zone-PAI 0.25, Zone-PAI 0.50 and Zone-PAI 0.75 respectively. Direct surface distur-
bances are assessed for each impact zone by comparing the baseline 1992 versus the current 2001
conditions as defined by land cover maps.
Copyright 2006 Taylor & Francis Group plc, London, UK
164
Figure 2. Spatial distribution of annual potential acid input predicted by CALPUFF model run
(adopted from in situ Oil Sands Application Firebag project)
3. SATELLITE DATA AND ITS PROCESSING
As mentioned before, in order to properly assess the impact of mining activities onto environment,
long-term observations are indispensable. So far, the only available source of long-term terrestrial
observations is data collected by several generations of NOAA AVHRR and LANDSAT satellites.
The AVHRR (advanced very high resolution radiometer) onboard a series of NOAA satellites have
provided daily measurements since 1979 with a coarse spatial resolution of 1.1 km. Despite nume-
rous limitations, the AVHRR historical time series is considered as a primary data set for terrestrial
monitoring at regional or global scales. LANDSAT MSS, TM and ETM+, with medium spatial re-
solutions (3080 m pixel size) have provided measurements since 1972 at a temporal resolution of
1618 days.
3.1. Scene selection and study area boundary
Remote sensing based land cover mapping was conducted over the AOS region within the extents
of a single LANDSAT scene.
The three images used in this study are:
LANDSAT 5 TM path: 43, row: 20 acquired on 11 June 1992;
LANDSAT 7 ETM+ path: 43, row: 20 acquired on 17 August 1999;
LANDSAT7 ETM+ path: 42, row 20 acquired on 23 August 2001.
The selected images span a 9-year period during the most intensive mining activity.
Copyright 2006 Taylor & Francis Group plc, London, UK
165
3.2. Scene geometric rectification
Geometric rectification of the selected images was performed using 20 geographically dispersed
ground control points (GCPs) with a first order polynomial and nearest neighbour resampling. An
average RMS error below 1.0 pixel was targeted as an acceptable rectification goal. For the 1992
LANDSAT scene an RMS error of 0.65 pixel was achieved. In order to correctly co-register the
1999 and 2001 scenes with the geo-rectified 1992 scene, pixels were identified in each scene corres-
ponding to geo-referenced pixels in the 1992 scene. Rectification of the 1999 and 2001 scenes was
performed using these pixels as GCPs. The RMS errors for each geo-rectification were 0.35 pixels
and 0.38 for the 1999 and 2001 scene, respectively.
3.3. Scene radiometric normalization
Variation in solar illumination condition, phenology and detector performance results in different-
ces in radiance values unrelated to the reflectance of the land surface. Radiometric normalization
was performed using the approach given by Du et al. (2001). The process consists of the following
steps: (1) selection of pixels pairs in overlap area; (2) principal component analysis for the selec-
tion of characteristic pixels; (3) calculation of gain and offset; (4) radiometric normalization.
3.4. Land cover classification
Selected scenes were classified into 15 land cover classes following the method described in Cihlar
et al. (1998) and Latifovic et al. (1999). The images were enhanced using a linear stretch as in Beau-
bien (1986) to equalize the dynamic range of each spectral band. All pixels in the full-resolution
enhanced image were classified into 150 spectral clusters using a K-means classifier. Further clus-
ter agglomeration was performed based on cluster spectral similarity and spatial proximity. The
procedure employs basic statistical parameters that determine between cluster spatial proximity and
within cluster spectral space data distributions. It uses the arithmetic mean of clusters to compute
the Euclidean distance between cluster pairs and the standard deviation to represent clusters ellip-
soid major axis and orientation relative to other clusters. Spectral clusters are first sorted according
to decreasing size and then examined for the most similar cluster pairs, starting with the smallest
clusters. The spatial criterion is based on the property that many spectral clusters represent gra-
dients within cover types and for these clusters there frequently exists a strong relationship bet-
ween their spatial and spectral characteristics. The agglomeration procedure utilizes this relation-
ship as a weighting function in computing cluster similarity. More detail on the agglomeration pro-
cedure is provided in Latifovic et al. (1999). The agglomeration yielded 55 significantly different
spectral clusters further grouped into 15 landcover ones (see Table 1). Post-classification refine-
ment included corrections in delineating forest cutblocks, urban and mining development leading
to the final land cover maps presented in Figure 3.
Two land cover maps were compared to each other to assess map consistency and then to refe-
rence data to assess accuracy. The reference data were project footprints obtained from Environ-
mental Impact Assessment (EIA) Studies of Millennium, Muskeg River and Fort Hills oil sands
projects. Other reference data used for the assessment included segments of AlPac Alberta vege-
tation inventory based on interpretation of 1:20,000 scale aerial photographs from 1991 and series
of reports on the Terrestrial Environmental Effect Monitoring (TEEM) Program of the Wood Bu-
ffalo Environmental Association (WBEA). Map consistency assessment was performed on a pixel
basis over the area that did not undergo change from direct impacts by comparison of pixel label
agreement. The assumption was that if both classifications map the same temporally stable pixel
with the same label, confidence in mapping accuracy increases. For example, 90% of all pixels ma-
pped as broadleaf forest in the 1992 land cover map were also classified with the same label in the
2001 land cover map. Overall average agreement of all classes weighted by class extent was 87%.
Copyright 2006 Taylor & Francis Group plc, London, UK
166
Higher consistency was achieved in mapping forest classes than in mapping non-forest classes
(92% and 75% respectively). An analysis of pixels in disagreement showed confusion between
spectrally similar or transitional land cover types such as broadleaf forest and broadleaf shrub,
treed and wetland or different coniferous sub-classes. Spectrally distinguishable classes, for exam-
ple, disturb area and forest classes were separated well with an average accuracy of 94%.
Figure 3. AOSR landscape classifications baseline condition in 1992 and current condition in 2001
derived from LANDSAT TM/ETM+ satellite data
Table 1. Description and hierarchical structure of landscape classes
Landscape class Description Label
VEGETATED NATURAL ORIGIN
DECIDUOUS FOREST
Deciduous forest Deciduous trees (aspen, birch, balsam poplar) greater than 60% canopy cover 1
MIXEDWOOD FOREST
Broadleaf
forest
Deciduous and conifer trees: aspen-pine, aspen-white spruce or balsam pop-
lar-white spruce where deciduous fraction is more than 60%. Canopy cover
more than 60%
4
Coniferous
mixed
Coniferous and deciduous trees where coniferous fraction is more than 60%.
Canopy cover more than 60%
5
CONIFEROUS FOREST
Black spruce
forest
Coniferous trees, mainly black spruce with fraction of white spruce and pine.
Mosses and shrub dominate in understory. Canopy cover greater than 60%
6
Black spruce
young forest
Generally young forest after old perturbation 9
Pine forest
Coniferous trees, mainly jack pine with lichen and some shrub in understory.
Canopy cover less than 60%
10
Pine black
spruce forest
Coniferous trees, co-dominant black spruce, jack pine. Mosses, forbs and shrub
in understory. Canopy cover more than 60%
13
Copyright 2006 Taylor & Francis Group plc, London, UK
167
Landscape class Description Label
NONFOREST
Treed
Sparse coniferous trees, less than 25% canopy cover. Commonly located on
poor sites associated with wetland. This class usually includes a proportion of
broadleaf vegetation, mostly shrub
8
Deciduous
shrub
Dominated by shrubland or grassland 11
Burn
Burn is usually characterized by sparse vegetation cover varying with the age
and intensity of burn
14
Open wetland Open fen, marsh or swamp 15
NON VEGETATION
Barren
land
Bare soil, rock space herbaceous and grass 20
Water 16
ANTHROPOGENIC ORIGIN
Disturbance Forestry cutblocks, build-up, roads and exploration lines 18
Industrial
development
19
NO DATA Clouds, shadows 17
4. METEOROLOGICAL DATA
The meteorological data set over the AOSR was extracted from the DS 472.0-TDL US and Canada
surface hourly observations generated by the Data Support Section (DSS) of the National Center of
Atmospheric Research (NCAR) in Colorado, USA. The complete data set includes measurements
of about 1000 meteorological stations across the USA and Canada for the period from 1976 to
2000. In this study, measurements of air temperature at the five stations closest to the AOS region
(Table 2) were extracted and analyzed. In addition to NCAR, the metrological data collected by Air
Quality Monitoring Stations of the Wood Buffalo Environmental Association (WBEA) for the pe-
riod from 1998 to 2002 were also included and analyzed.
Table 2. Geographical locations of meteorological stations
Stations Station ID
Longitude
[dd]
Latitude
[dd]
Fort McMurray YMM -111.22 56.65
High Level YOJ -117.17 58.62
Peace River YPE -117.43 56.23
Fort Chipewyan YPY -111.12 58.77
Slave Lake YZH -114.78 55.30
5. CUMULATIVE AND PROJECT SPECIFIC IMPACTS ON LAND COVER DISTRIBUTION.
RESULTS AND DISCUSSION
Cumulative and project specific impacts on the landscape are quantified using a post-classification
change detection method. This method assumes that reference and compared images are classified
to a common legend and that the classification method utilized for landscape mapping provides
a high accuracy for both images. Landscape changes are simply detected as differences between
Copyright 2006 Taylor & Francis Group plc, London, UK
168
pixels labels. Changes in the landscape distribution in the AOSR for the period 19922001 were
quantified as the difference between classified images from 1992 and 2001. The assessed areas
spatially coincide with the impact zones. The Zone-PAI 0.25 that incorporates the other two zones
encompasses 715,094 ha. In 1992 Zone-PAI 0.25 was dominated by native forest vegetation, with
natural forest landscape classes occupying 75.6% of the area (Fig. 4 and Tab. 3).
Figure 4. Area occupied by forest, non-forest and anthropogenic landscape classes in AOSR in 1992 and 2001
Table 3: AOSR landscape fraction and class distribution for baseline 1992 and current 2001 conditions
Land cover type 1992 2001
[ha] PAI > 0.25 PAI > 0.50 PAI > 0.75 PAI > 0.25 PAI > 0.50 PAI > 0.75
Broadleaf 95,417 39,589 12,639 104,410 32,819 6,934
Coniferous Mixed 35,860 14,521 5,325 39,199 15,333 6,122
Broadleaf Mixed 69,635 27,135 9,475 75,212 29,570 9,769
Coniferous HD 238,475 95,031 47,577 212,751 79,730 38,805
Coniferous Young 58,551 23,193 8,467 62,835 26,679 11,383
Coniferous JP 18,312 3,066 1,705 16,912 2,154 939
Coniferous JP + BS 24,529 6,916 4,156 16,295 5,158 2,564
Coniferous LD 84,855 29,111 14,444 104,587 38,225 17,752
Broadleaf Shrub 24,595 9,609 4,235 9,961 2,857 758
Burn 11,171 481 220 3,447 40 10
Open wetland 1,193 357 171 231 77 47
Water 16,571 8,062 5,034 12,507 5,892 3,882
Barren Land 8,698 2,495 1,225 5,819 2,887 1,294
Logging 6,838 2,645 431 12,920 5,097 709
Mining 20,393 17,482 14,423 38,008 33,177 28,559
SUM 715,094 279,693 129,526 715,094 279,693 129,526
Forest 540,779 209,451 89,343 527,613 191,442 76,516
Non-Forest 147,083 50,115 25,328 136,552 49,977 23,742
Anthropogenic 27,231 20,127 14,854 50,929 38,274 29,268
Forest (%) 75.6% 74.9% 69.0% 73.8% 68.4% 59.17%
Non-Forest (%) 20.6% 17.9% 19.5% 19.1% 17.9% 18.3%
Anthropogenic (%) 3.8% 7.2% 11.5% 7.1% 13.7% 22.6%
Copyright 2006 Taylor & Francis Group plc, London, UK
169
Coniferous forest covered approximately half of the study area (48%), whereas deciduous and
mixed-wood forest covered 13% and 15% respectively. The remaining 24% of the area consisted
of broadleaf shrub, non-forested and open wetland, water bodies, developed land (including active
surface mining, roads and built-up), barren land and forest of anthropogenic origin (cutblocks). Two
anthropogenic classes together covered 3.8% (or 27,231 ha) of the area. In 2001 Zone-PAI > 0.25
exhibited a similar landscape class distribution as in 1992 but with a slight decrease in the area co-
vered by natural vegetation. From 1992 to 2001, forest-covered area had decreased from 75.6 to
73.8%, while the area covered by natural non-forest vegetation decreased from 20.6 to 19.1%. Du-
ring the same period, the areas occupied by the two-anthropogenic classes have increased from 3.8
to 7.1% (Table 3). Mining area almost doubled from 20,393 to 38,008 ha and forest cutblocks in-
creased from 6,838 to 12,920 ha.
The spatial distributions of the areas that have undergone change are illustrated in Figs. 5 and 6.
Figure 5. Anthropogenic landscape classes spatial distribution in AOSR in 1992
Figure 6. Anthropogenic landscape classes spatial distribution in AOSR in 2001
Copyright 2006 Taylor & Francis Group plc, London, UK
170
These figures show an increase in the disturbed area resulting mainly from new developments
north and west of the mining development that existed in 1992. An increase in logging activities is
also evident, in most cases spatially adjacent to the existing mine site or planned mining develop-
ments. However, the amount of logging activity associated with the existence of mining develop-
ment cannot be evaluated from the data used in this assessment. An analysis of the landscape dis-
tribution as a function of proximity to the mining development offers additional insights into the
class area changes. Table 3 specifies the landscape class distributions in three surrounding rings at
different distances from the mine site. Between 1992 and 2001 the landscape distributions in the
Ring-PAI 0.25, the farthest from the mine sites, are more similar than in the Ring-PAI 0.75, which
is the closest to the mining development. This difference can be explained mainly by direct impacts
due to new surface mining developments in the Ring-PAI 0.75 i.e. the Muskeg River, Millennium,
Steepbank, Aurora and North mines, and secondly by logging activities. Areas occupied by the
broadleaf and coniferous high-density forest classes in 2001 compared to their respective shares in
1992 have shrunk by 4 and 7%, respectively. Figure 7 shows the changes in forest, water bodies,
mining and logging classes relative to 1992 in each ring.
Figure 7. Difference in class area relative to 1992: (a) forest; (b) water bodies; (c) mining; (d) logging
The graph (b) in Figure 7 merits special attention since it exhibits (~25%) decrease in the area
of open water bodies. Significant decreases are evident in the vicinity of the Mildred Lake and
North mines, suggesting impacts on local hydrology due to these long-term mining operations.
Examination of LANDSAT scenes from 1999 and 2002 as well as an ASTER scene from 2000 (15 m
spatial resolution) confirmed shrinking of water bodies in the area south west of the Mildred Lake
and North mines. These findings need to be confirmed by more detailed analysis of in situ hy-
drological measurements because the influence on wetland and surrounding vegetation due to a lo-
wering water table could cause significant impacts. The other three graphs in Figure 7 show an
overall decrease in natural forest in 2001, an increase in area affected by logging activities for all
three rings and a significant increase in disturbed area by surface mining. The results shown on the
a) b)
c) d)
Copyright 2006 Taylor & Francis Group plc, London, UK
171
graph (c) indicate that the area affected by mining in 2001 has doubled in size since 1992. Regard-
less of relatively high mapping consistency we consider the direct impact presented here as appro-
ximate due to the limited capability of medium resolution data to depict small changes less than
a pixel size (30 m), such as seismic lines and small clearings for exploration drill holes.
6. CONCLUSIONS
A remote sensing based approach for quantifying primary impacts of surface mining is presented
and applied for the assessment of the mining development in Athabasca Oil Sands Region (AOSR).
Affected areas were identified as the difference between land cover maps derived from LANDSAT
data (30 m resolution) acquired in 1992 and 2001. Maps produced from remote sensing data pro-
vide information for subsequent impact assessments from surface mining development on land
cover, as well as forming the basis for reclamation planning and monitoring.
Understanding the cause of variability in vegetation conditions allows for distinguishing natural
from human-induced perturbations of ecosystems; an important pre-requisite in quantifying environ-
mental impacts caused by mining and other developments. However, it is important to point out that
the procedure presented here requires a somewhat complicated procedure for correction and nor-
malization of remote sensing data, including across sensor calibration, atmospheric correction, screen-
ing for cloud and other contaminants often present on remote sensing observations. In addition to
high quality remote sensing data, a comprehensive knowledge of conditions that a study area was
subject to in the past and present is essential for appropriate interpretation leading to higher confi-
dence in obtained results.
REFERENCES
Beaubien J. 1986: Visual Interpretation of Vegetation through Digitally Enhanced Landsat Mss Images.
Remote Sens. Rev. 2, 1143.
Cihlar J., Latifovic R., Chen J.M. & Li Z. 1999: Testing Near-Real Time Detection of Contaminated Pixels in
AVHRR Composites. Can. J. Remote Sens. 25, 160170.
Cihlar J., Hung L. & Xiao Q. 1996: Land Cover Classification with AVHRR Multichannel Composite in
Northern Environments. Remote Sens. Environ. 58, 3651.
Cihlar J., Xiao Q., Beaubien J., Fung K. & Latifovic R. 1998: Classification by Progressive Generalization:
a New Automated Methodology for Remote Sensing Multichannel Data. Int. J. Remote Sens. 19, 26852704.
Du Y., Cihlar J., Beaubien J. & Latifovic R. 2001: Radiometric Normalization, Compositing and Quality Con-
trol for Satellite High Resolution Image Mosaics over Large Areas. IEEE Trans. Geosci. Remote Sens. 39 (3).
Environmental Protection Agency (EPA). 1995. A Users Guide for the CALPUFF Dispersion Model. EPA-
454/B-95-006.
Latifovic R., Cihlar J. & Beaubien J. 1999: Clustering Methods for Unsupervised Classification. In: Procee-
dings of the 21st Canadian Remote Sensing Symposium, II-509-II-515, Ottawa.
Prakash A. & Gupta R.P. 1998: Land-use Mapping and Change Detection in a Coal Mining Area a Case
Study in the Jharia Coalfield, India. Int. J. Remote Sens. 19 (3), 391410.
Ress W.G. & Williams M. 1997: Monitoring Changes in Land Cover Induced by Atmospheric Pollution in
Kola Peninsula, Russia, using Landsat MSS data. Int. J. Remote Sens. 18 (8), 17031723.
Schmidt H. & Glaesser C. 1998: Multitemporal Analysis of Satellite Data Their Use in the Monitoring of the
Environmental Impact of Open Cast Mining Areas in Eastern Germany. Int. J. Remote Sens. 19 (12),
22452260.
Wright P. & Stow R. 1999: Detecting Mining Subsidence From Space. Int. J. Remote Sens. 20 (6), 11831188.
Copyright 2006 Taylor & Francis Group plc, London, UK