DETERMINING THE SIZE OF LINEAR STOCKPILE FOR THE AFŞİN

Journal of Science and Technology
1 (2), 2007, 167-183
©BEYKENT UNIVERSITY
DETERMINING THE SIZE OF LINEAR
STOCKPILE FOR THE AFŞİN - ELBİSTAN (B)
POWER STATION (TURKEY)
Suphi URAL
Çukurova University, Department of Mining Engineering, Adana, Türkiye,
[email protected]
ABSTRACT
The aim of this study is to show how to estimate optimum size of a linear coal
stockpile for the Afşin - Elbistan (B) power station using geostatistical
modeling techniques. After a description of the geology of the area a
geostatistical study based on 2968 core samples of 193 drill holes is presented.
The coal samples present a symmetric frequency distribution. Experimental
variograms representing three main directions - strike plunge, cross-dip and
downdip directions - are calculated and a simple spherical-type model is fitted
to the variograms. The model parameters are validated. Storages capacities for
buffering and homogenisation are estimated taking into account of longer
operation standstill of the Çöllolar mine and daily fluctuations of quality
parameters. In order to estimate the storage capacity of stockpile and number
of stacking layers, block variance of varies size blocks are calculated utilising
F chart. Finally, numbers and dimensions of prismatic type linear stockpile
bunkers are determined.
Key Words: Afşin-Elbistan, geostatistics, coal, stock pile, block variance
ÖZET
Bu çalışmanın amacı, jeoistatistiksel modelleme teknikleri kullanılarak bir
kömür stok sahasının optimum boyutlarının nasıl kestirilebileceğini
göstermektir. Çalışmada, Afşin-Elbistan (B) termik santralinde, lineer bir
kömür stok sahasının stoklama kapasitesi ve boyutlandırılmasına ilişkin olarak
yapılan çalışmalar anlatılmıştır. Bölgenin jeolojik yapısı tanıtıldıktan sonra,
193 adet sondaj kuyusundan sağlanan 2968 adet karot numunesine ilişkin
veriler esas alınarak sahanın jeoistatistik modeli oluşturulmuştur. Üç ana
yönde de deneysel variogramlar ve bu variogramlara uygun modeller
seçilmiştir. Çöllolar maden işletmesinde ortaya çıkabilecek uzun süreli
duruşlar, kömür harmanlama gerekleri ve kalite kontrolüne ilişkin kısıtlarda
167
Determining The Size Of Linear Stockpile For The Afşin - Elbistan (B) Power Station
(Turkey)
dikkate alınarak, stok sahasının kapasitesi kestirilmiştir. Bu hesaplamalarda,
farklı büyüklüklere sahip blokların varyansları F kartları kullanılarak
hesaplanmıştır. Son olarak stok sahasındaki bunker sayıları hesaplanmış ve
bunlar boyutlandırılmıştır.
Anahtar kelimeler: Afşin-Elbistan,
varyansı
jeoistatistik,
kömür ,stok sahası , blok
1. INTRODUCTION
The targets of the mining-technical planning are to guarantee that the power
station operations will have to be stopped at no time due to lack of fuel,
fluctuations with regard to volume and quality of the product should be
balanced and the total calorific value of mining field should be used as fully as
possible by blending the coal. Stockpile with its capacity for long-term low
frequency quality variation filtering has been widely accepted as a key method
of raw coal homogenization. It has dual role within the raw coal preparation
system, being a buffer storage unit and capability for homogenization of raw
coal if properly designed and operated. The purpose of the present paper is to
estimate the size of linear stockpile for the Afşin - Elbistan (B) power station.
The Afşin-Elbistan lignite deposit with its 3.4 billion metric tons of reserve is
the biggest lignite basin and one of the most important resources for electrical
energy production in Turkey [10]. The Elbistan lignite is best suited to use in
thermal power generation and are capable of supplying 16 x 340 MW power
stations. The lignite deposit is located north of Elbistan, Kahramanmaras
Province, in southeast Turkey (Figure 1). The Çöllolar mining field is the
second open cast mine to feed Afsin-Elbistan (B) power station with the
installed capacity of 1376 MW. Afsin-Elbistan (A) power station has been
under operation since 1984 [11].
Figure 1. Geographic location of study area
168
Suphi Ural
The Elbistan basin encompasses an area of 900km , and its centre lies about
1150m above mean sea level. The sedimentation took place mainly in the
upper Pliocene period. The Collolar sub-basin of the Afsin-Elbistan basin is a
multi-layer lignite deposit with associated clays and gyttja partings. The
thickness of the coal bed varies from 60m to 80m respectively. The basal part
of the coal beds consists of clay, marl and sandy clay. These formations are
overlain by gyttja, which may also occur in the form of calcareous ooze,
enveloping coal beds. The thickness of the gyttja is up to 80m. Alluvium,
occurring very widespread along Hurman Stream, consists of gravels, sands
and silt [5]. Lithological profile of the Collolar field is presented in Figure 2.
Figure 2. Lithological profile of the Collolar coal -field
2. PRESENTATION OF THE DATA SET
The information available comes from 193 holes drilled vertically on about a
200m to 250m grid (Figure 3). The coordinates and the number of the hole,
and the results of the analyses for raw coal are supplied. The raw coal had been
analysed for moisture content, ash content and calorific value. Core recoveries
in the lignite are between 92% and 96% with an average of 94%.
169
Determining The Size Of Linear Stockpile For The Afşin - Elbistan (B) Power Station
(Turkey)
47000^
46000
N
45000^
t
44000^
43000^
42000^
41000^
21000
22000
23000
24000
25000
26000
27000
28000
29000
Figure 3. Drill hole location map of Çöllolar coal-field
The deposit is made up of a great many layers of coal and intercalation. In this
case the question of defining what is to be considered as "coal" and what is the
"waste" is very real. This choice is governed by the technical constraints
imposed by the mining method used. To be mineable, a coal seam must be at
least 50cm thick. Similarly any intercalation less than 25cm thick is mined
with the coal. These two constraints are applied to each drill hole in turn to
distinguish mineable coal from waste. Figure 4 illustrates this procedure for a
typical drill hole. In the first step the thicknesses of all contiguous coal seams
are summed. Any coal seams that are too thin to be mined are then redefined
as waste. Lastly any narrow bands of waste sandwiched in between exploitable
coal seams are included with the coal.
Since the coal is destined for power generation the variables of interest include
calorific value and ash content. The core samples are composited into a
constant length of 5m from sample lengths that varied from 0.5m to 5m. The
basic statistics of the variables to be studied is calculated in order to check
outliers (extremely large or small values) or nonhomogeneous population. The
histograms of the variables show that they present nearly normal distribution.
2968 composite samples have a mean of 1167 kcal/kg calorific value and
37.2% ash content in dry coal (Figure 5).
170
Suphi Ural
ORIGINAL
DRILL HOLE
MINING
DRILL HOLE
Waste
Coal
Waste< 0 5
m
Coal
1
1
1
1
Coal
Waste
Waste
Coal< 0.5 m
Waste
Figure 4. Procedure for calculating the quantities of coal and waste in a drill hole
Total number of samples: 2968
Mean: 1167; Variance: 99225
Total number of samples: 2968
Mean: 37.2; Variance: 42.25
14
12 n
12
10 8 6 4
2
0
İlli JUL
»il
CN
CN
„ojlil
E
CO CO
Ash content in dry basis (%)
Calorific value in raw coal (kcal/kg)
(a)
(b)
Figure 5. Histograms of variables: (a) ash content; (b) calorific value
171
Determining The Size Of Linear Stockpile For The Afşin - Elbistan (B) Power Station
(Turkey)
3. VARIOGRAM ANALYSIS AND MODELING
Geostatistical theory is based on the observation that the variability of such
quantities as grade and thickness has a particular spatial structure. The grades
z(x) and z(x+h) at point x and x+h are correlated; this spatial correlation
depends on the vector h, which separates the two points and decreases as the
distance I h I increases. Beyond a particular distance the two quantities are no
longer correlated. The spatial variability can also differ from one direction to
another [7]. A characteristic behaviour or structure of the spatial variability can
be discerned behind a locally erratic aspect. Geostatistics consider these
variables as regionalized variables and both the random and structured aspect
of the regionalized variables are expressed in the probabilistic language of
random functions, every point of the deposit defines a random variable and the
value at that point is interpreted as a realization of that random variable. A
random function is seen as a set of random variables defined at each point in
the deposit and a regionalized variable is a realization of the random function.
In geostatistics the spatial variability of a regionalized variable is characterized
by the variogram function, 2y(h) [6].
2y(h) = Var[Z(x) - Z(x+h)];
where h is the separation vector, y(h) is the semi-variogram and Z(x) is a
random variable defined at point x. The variogram defines a few parameters
like the zone of influence in all directions (which is named range, a), sampling
error (which is named nugget effect, C0) and the variance of the samples
(which is named sill, C) and is used in estimating an unknown value of the
regionalized variable under study. Experimental variograms and model
variograms of the variables of calorific value and ash content are seen in
Figure 6 and Figure 7. The first experimental variogram of calorific value is
computed for the downdip direction (Figure 6 (a)); this variogram shows a
range of around 40m with a nugget effect value 2400 kcal/kg. The
experimental variograms cross dip (Figure 6 (b)) and in strike-plunge (Figure 6
(c)) indicate the existence of geometric anisotropy, and the ranges in these
directions are 750m and 500m, respectively. The following three-dimensional,
simple, spherical type of variogram model is adopted
y(h) = 22080
y(h) = 0
h > a())
h = a(9)
172
Suphi Ural
The experimental variogram of ash content for the downdip direction (Fig.
7(a)); shows a range of around 40m. The experimental variograms cross dip
(Fig. 7(b)) and in strike-plunge (Fig. 7(c)) indicate the existence of geometric
anisotropy, and the ranges in these directions are 682m and 455m,
respectively.
(a)
Experimental variogram
Model
E
™ 30000
<0 d) ^ 20000
S
~
10000
0
0
500
1000
Distance (metres)
(b)
•Experimental variogram
Model
Jw
C
m
a re 40000
r c
o JC 20000
e
0
> rs
0
>
r
a
200 400 600 800 1000 1200
Distance (metres)
(c)
Figure 6. Experimental variograms and models for calorific value: (a) downdip;
(b) cross-dip; (c) strike plunge
173
Determining The Size Of Linear Stockpile For The Afşin - Elbistan (B) Power Station
(Turkey)
-Experimental variogram
( 15
ğ
10-
İ
•g
0
2
10
20
30
40
50
Distance (metres)
(a)
•Experimental variogram
3
e 15
10
>
m
a 5
0
>
0
200
400
600
800
Distance (metres)
(b)
"Experimental variogram
e 15
>
Ü
10
5
0
0
500
1000
1500
Distance (metres)
(c)
Figure 7. Experimental variograms and models for ash content: (a) downdip; (b) crossdip; (c) strike plunge
The following three-dimensional, simple, spherical type of variogram model is
adopted:
y(h) = 10 1.5-
h
a(9)
•-0.5
(
h
I a(G)
Y(h) = 12
^3
+2
h <a(G)
h > a(G)
174
Suphi Ural
y(h) = 0
h = a(d)
The parameters of the variogram model are validated by use of the backestimation technique, in which each sample is removed from the data set in
turn and then estimated from the remaining data and the variogram model.
For conditionally unbiased estimates, the linear regression of actual values on
estimated values should be close to a 45° line-that is, an intercept of zero and a
slope of 1 [4].
The linear regressions of actual values on estimated values are presented in
Figure 8 and Figure 9, respectively. Figures 8 and 9 can be taken as a
verification of the model chosen.
~
t
1600
!
1400 -
t
100
1
Act. Value = 0.889Est. Value +
85.88
R
^
=
W
Ü 1000 -
=
o
<
•,
800 J — ksC?*
'
800
1000
Î
P
^
'
*
1200
,
,
1400
1600
Estimated Value (kcal/kg)
Figure 8. Linear regression of actual values on estimated values for the calorific value
of coal
Act. Value = 0.898Est. Value +
20
30
40
50
Estimated Value (%)
Figure 9. Linear regression of actual values on estimated values for ash content of coal
175
Determining The Size Of Linear Stockpile For The Afşin - Elbistan (B) Power Station
(Turkey)
4. STORAGES CAPACITY FOR BUFFERING
The following points favour the use of stockpile as a buffer storages unit: (1)
activities both at the mine and during production need not interact
dynamically; (2) a steady supply of homogeneous raw coal to the power
station is ensured; and (3) reduces demand for complex automation of
subsequent processing stages.
The minimum buffering storages capacity of a power station coal stockpile can
be determined
Qb = T P (SE /LCV) /1000
Vsb = Qb / p
where Qb is stocking capacity of stockpile, tonne; T is duration of
maximum stoppage time for coal production, hour; P is installed capacity
of power station, kW; SE is specific energy consumption of power station,
kcal/kWh; LCV is lower calorific value of raw lignite, kcal/kg; Vsb is
stocking volume of stockpile, m3; and p is average bulk density of the coal,
t/m3.
It is reported that, duration of the longer operating standstill is 264 hours
[9]; the installed capacity of Çöllolar power station is 1376000 kW,
specific energy consumption of power station is 2300 kcal/kWh, average
bulk density of the Elbistan coal is 0.85 t/m 3 and average calorific value is
1167 kcal/kg. The stocking volume of stockpile can be calculated as
Qb = (264 x 1376000 x (2300 /1167))
/1000
Qb = 7160001
Vsb = 716000 / 0.85
Vsb =842353 m3
5. STORAGES CAPACITY FOR HOMOGENISATION
Raw coal arriving from the mine has certain unwanted characteristics
fluctuations in their quality composition that must be filtered out if a consistent
good quality coal feed is to be achieved. The following considerations favour
the use of stockpile as a raw coal homogenisation unit: (1) sufficient overall
input/output variability reduction ratio can be achieved, (2) raw coal from
poorer deposits with large scale heterogeneity can be utilised, (3) mixing of
various component raw coals can be achieved and (4) power station demands
for improved product quality can be met.
176
Suphi Ural
Several attempts have been made to determine the size of stockpile utilising
the theory of regionalised variables [1, 3] and also using unrealistic black-box
macroscopic design techniques based on standard Gaussian ideal statistical
formulae for assumed approximately normal distributed, statistically
independent coal quality variation [8].
Geostatistical approach has been applied to determining the size of stockpile.
In many cases a regionalized variable is defined as the average over a certain
volume or a surface rather than a point. The basic volume on which a
regionalized variable is defined is called its support. If we change the support
we obtain a new regionalized variable, which is related to the preceding one
but which has different structural characteristics. In case a stockpile is being
considered to reduce stockpile input quality standard deviations, one can
compute which size it should be so that stockpile output quality standard
deviations are reduced to a certain level. Geostatistics can be used to calculate
the variability of the average value of blocks of certain sizes, in other words
the coal is homogenised over the size of the stockpile [1]. For that reason, it is
essential to know the variance of blocks of a given size in the deposit. For the
sake of generality, the supports will be referred to as volume v (i.e. daily
production) and volume V (i.e. deposit). After Kriging's relationship [2] it is
known that the variance of a block within a deposit is equal to the variance of a
point within a deposit minus the variance of a point in the block. Thus all that
is needed is the within block variance for a point as the first one is
experimentally available. The variance of a block v within a large block V will
be obtained by (David, 1977)
a2 (v/V)
= a2(0/V)-a2(0/v)
= F(V)-F(v)
;
where v is volume of daily production, V is volume of stockpile, a2(0/V) is
variance of a point within stockpile and a2(0/v) is variance of a point within
daily production.
To determine the size of the stockpile, one has to find the size of block, which
has the value of a2(0/V) (David, 1988)
a2(0/V)
= a2 (v/V)
+
a2(0/v)
Assume that the effect of stockpiling and subsequent reclaiming is perfect,
in other words the coal is homogenized over the size of stockpile. Regarding
the quality parameters as calorific value and ash content, the necessary size of
stockpile for perfect homogenization can be calculated as follows:
Calorific value: The variogram of calorific value is spherical with the
horizontal ranges across-dip and a stnke-plmge are 750m and 500m, the vertical range
177
Determining The Size Of Linear Stockpile For The Afşin - Elbistan (B) Power Station
(Turkey)
is 40m with the nugget effect value of 2400 kcal/kg and the sill of 19680
kcal/kg. Each day 75m x 50m x 15m block, oriented as the main variation
axes, is mined out. Using the chart of Fig. 10 and taking the anisotropies into
account to calculate a2(0/v)
F [(75 / 750), (50 / 500), (15 / 40)] = 0.21
a2(0/v) = F(v) C = 19680 x 0.21 = 4133
It is desired to keep daily fluctuation of calorific value to 30 kcal/kg. To
determine the size of the stockpile, one has to find the size of block, which has
a variance (a2(v/V)) of 900 kcal/kg, as stockpiling and reclaiming to send the
power station will show the variability of these larger blocks. The block must
be 105m and 15m in two of its dimensions. To obtain a variance of 900
kcal/kg, we need a block within which the variance of a point will be
a2(0/V) = 900 + 4133= 5033
F [(h/750),(105/500),(15/40)]= 5033 /19680
F[(h/750), (0.210), (0.375)] = 0.26
Figure 10. Variance aa2(0/V) of a point within a parallelepipedic block (h, h, l) function
(1/C) F(h/a, h/a, l/a) [2]
178
Suphi Ural
Using the chart of Figure 10, we can determine the third dimension of the
block that should be 0.21 x 750m = 157.5m long. So, volume of the block is
V =157.5m x 105m x 15m = 248063 (bank) m3
and the volume of the stockpile regarding the calorific value is
Vc = V x
Vc = 248063 (bank) m3 x 1.18
Vc = 292714 (loose) m3,
where Vc is volume of the stockpile regarding the calorific value, m3; and y is
average swell factor of Elbistan coal.
Ash content: The variogram of ash content is spherical with the horizontal
ranges across_dip and astrike_plunge are 682m and 455m, the vertical range is 40m
with the nugget effect value of 2% and the sill of 10%. Using the chart of Fig.
10 and taking the anisotropies into account to calculate a2(0/v)
F[(75 / 682), (50 / 455), (15 / 40)] = 0.22
a2(0/v) = F(v) C = 10 x 0.22 = 2.2
It is desired to keep daily fluctuation of ash content to 1.0%. To determine the
size of the stockpile, one has to find the size of block, which has a variance
(a2(v/V)) of 1.0%. The block must be 105m and 15m in two of its dimensions.
To obtain a variance of 1.0%, we need a block within which the variance of a
point will be
a2(0/V) = 1.0% + 2.2% = 3.2%
F [(h/682), (105/455), (15/40)] = 3.2 /10
F[(h/682),
(0.230), (0.375)] = 0.32
Using the chart of Fig. 10, we can determine the third dimension of the block
that should be 0.23 x 682m = 157m long. So, volume of the block is
157m x 105m x 15m = 247275(bank) m3
and the volume of the stockpile regarding the ash content is:
Va = V y
Va = 247275 (bank) m3 x 1.18
179
Determining The Size Of Linear Stockpile For The Afşin - Elbistan (B) Power Station
(Turkey)
Va = 291785 (loose) m3
If it is desired to keep daily standard deviation of calorific value and ash
content to 30 kcal/kg and 1% respectively, the total volume of stockpile is to
be at least 292714m 3 .
6. STACKING AND RECLAIMING
Most of stockpile equipment manufacturers utilize design methods that are
essentially based on standard statistical macroscopic relationships [3]. These
are used primarily to specify the number of layers of coal (Figure 11) that
would be theoretically required in the stockpile to achieve a certain variability
reduction ratio within the pile. The variability reduction ratio is also shows the
homogenisation effect.
n = c /c
where n is homogenisation effect, c o is stockpile output sample-to-sample
quality standard deviation and c is stockpile input sample-to-sample quality
standard deviation.
Figure 11. Representing of stacking in layered linear stockpile
The assumptions made for the specification of the number of layers of coal are:
(1) that the variability of the raw coal to be homogenised are of a total random
character with no underlying trends and present a near normal or Gaussian
probability density function, (2) that adjacent quality sample values for any
sample quantity/size are statistically independent. The sampling distribution of
means is related to the original distribution for the data in having the same
mean and a variance proportional to the original variance divided by the
number of elements in each sample [8]
a = k a2 j / n
where n is number of elements that can be interpreted as being equal to the
number of layers (N) and k is the proportionality constant the value of which is
x
180
Suphi Ural
dependent on the amount of sample to sample correlation. If it can be assumed
that the variation within the individual reclaimed slices is so insignificant to be
of no interest and sample-to-sample correlation is zero then the a x can be
expressed as:
a2x = a2o
a o = a2, / N
N = a2, / a2o
Relating the Kriging's relationship [2], number of layers for perfect blending is
equal the variance of a point within the bunker divide to the variance of a point
within a layer.
N
_ a2
(0/Vb)
2
(0/Vl)
a
2
where a (0/Vb) is variance of a point within the bunker, Vb is volume of a
bunker, m 3 ; a2(0/Vl) is variance of a point within a layer and Vl is volume of a
layer, m3.
If it is desired to keep stockpile daily output standard deviation of calorific
value to 30 kcal/kg and a2(0/Vb) equal to a2(0/V) and a2(0/Vl) equal to a2(0/v);
number of layers for the variable of calorific value can be calculated
ir
N _
5033
900
N =6
7. DIMENSIONS OF STOCKPILE
Total volume of the stockpile and number of bunkers can be given by
Vsb > Vc and Vc = Va
Vb = Vc
Vb = 292714 m3
Nb = Vsb / Va
Nb = 842353 m3 /292714 m3
Nb =3
where Nb is number of bunkers. Total volume of the stockpile is equal to
buffering storages capacity and stockpile can consist of three identical bunkers.
The maximum storages volume of a prismatic stockpile bunker is given by
n
H3
Vb = — x
— +
3 tan2 a
181
H2L
tana
Determining The Size Of Linear Stockpile For The Afşin - Elbistan (B) Power Station
(Turkey)
where L is the length of a bunker, m; H is total height of pile, m; a is the angle
of repose of the coal; and D is the base with of the pile, m. Base with of the
pile can be expressed as
D = J H .
tana
Length of bunker can be expressed as:
2(Vb
D2H
6
L =
)
DH
It is reported that [9], total height of pile is 14m and the angle of repose of the
coal is 32°. The base width of the pile is
D =
2xI4
D = 46.6m
tan SI
Length of bunker is
2x(2927I4
L=
n
6
46.6 xI4
L = SSSm
x(46.6)2xI4)
Proposed dimensions of the stockpile bunkers are presented in Figure 12.
Figure 12. Modeling of prismatic type linear stockpile
182
Suphi Ural
8. CONCLUSIONS
The main findings of this study are the following:
The deposit is made up of a great many layers of coal and intercalation. To be
mineable, a coal seam must be at least 50 cm thick. Similarly any intercalation
less than 25 cm thick is mined with the coal. These two constraints are applied
to each drill hole in turn to distinguish mineable coal from waste.
Geostatistical modeling techniques are good tool for the coal quality
controlling purposes. The study presented here indicates that three-dimensional
variograms of the quality parameters as calorific value and ash content of the
Çollolar coal deposit can be modeled by spherical models. The models are then
used to estimate fluctuation of daily production and storage capacity of the
stockpile for perfect blending. The results of this study show that; the
necessary storage volume for blending is 292714 m 3 regarding to the coal
quality parameters. As a result, total volume of the stockpile is to be at least
842353m 3 and consists of 3 identical bunkers. The length and base width of
each bunker is to be 888.0m and 46.6m respectively.
REFERENCES:
1. Armstrong, M. and Carignan J.; Géostatististique linéaire, Application au domaine
minier, Paris: Les press de l'Ecole des Mines, (1997), 113 p.
2. David, M. Geostatistical ore reserve estimation, New York: Elsevier, (1977), 364p.
3. David, M.; Handbook of applied advanced geostatistical ore reserve estimation, New
York: Elsevier, (1988), 215 p.
4. Dowd, P.A.; Algorithms for three-dimensional interpolation between planar slices.
In Earnshaw R.A. ed. Fundamental algorithms for computer graphics, NATO ASI
Series F17, Berlin: Springer-Verlag, (1985), 531-54.
5. Gold, O. Lignite Deposit Afsin Elbistan Feasibility Report, Volume I, Cologne,
(1969), 354 p.
6. Journel, A.G. and Huijbregts Ch. J.; Mining geostatistics, London: Academic Press,
(1978), 600 p.
7. Matheron, G.; Principal of geostatistics. Economical Geol., (1963), 58, 1246-66.
8. Schofield, C. G.; Homogenisation/Blending Systems Design and Control for Mineral
Processing, Clausthal: Trans Tech Publication, (1980), 315 p.
9. Turkish Electricity Enterprises; Project Report for Afsin Elbistan Çollolar (B),
Volume I, Ankara, (1978), 123p.
10. Unver, O. and Kalafatcioglu, A.; Utilization of low calorific value Afsin-Elbistan
Lignites, Braunkohle - Tagebautechnik, Sondernummer EII, December (1979), 57-64.
11. Ural, S. and Onur A.H.; Application of continuous surface mining technology in
Afsin-Elbistan lignite basin, Transactions of the institution of mining and metallurgy
(Section A), (1998), Vol. 107, A48-51.
183