# DETERMINATION OF THE HORIZONTAL DEFORMATION FACTOR FOR MINERAL AND FLUIDIZED DEPOSITS EXPLOITATION.

INTRODUCTION

Avershin (1948) published the relationship between horizontal displacement and tilt, known in the literature as the Avershin's hypothesis. It assumes a linear proportionality between the horizontal displacement vector and the subsidence trough tilt vector. The analysis of the results of horizontal deformations and movement of points due to the convergence of salt caverns used for storing gaseous and liquid fuels together with the analysis of the results of measurements made during the extraction of fluid deposits show clearly that the value of the Avershin's proportionality index, known in the literature as the horizontal deformation factor B, for the deformations over the cavern fields in a salt rock mass or surface deformations resulting from the extraction of fluid deposits are vastly different from the standard values of this index as provided by Budryk (1953). It was then determined based on measurements of deformations caused solely by the hard coal mining operations.

The solution of the proper determination of the horizontal deformation factor B is an extremely important aspect in the realistic estimation of the expected deformation, in particular of the horizontal deformation, as the most accurately describing the expected damage (Knothe, 1953b; Budryk and Knothe, 1956). This is a necessary action, especially in the increasingly difficult general situation of mining in Europe, where this industry has to deal with increasingly frequent restrictions and formal and political demands.

THEORETICAL FOUNDATIONS

Based on the solutions provided by Salustowicz (1953) and Litwiniszyn (1953), Budryk (1953) proved that the maximum value for horizontal deformation [[epsilon].sub.max] for the flat issue (so-called: infinite semi-plane) can be described with the following formula:

[[epsilon].sub.max] = -c * [T.sub.max], (1)

where:

[T.sub.max] - maximum tilt of the subsidence trough and

C - constant.

Based on the analysis of measurement data for coal seams with the tilt up to 35[degrees], Budryk calculated the c parameter value within the limit of 0.52 / 0.60.

In his magnum opus describing the concepts for calculating horizontal deformations and movements, Budryk (1953) writes that: "in actuality, the largest horizontal deformations usually do not exceed:

[[epsilon].sub.max]= 0.60 * [T.sub.max]". (2)

This shows that the solution provided by Budryk is applicable to the ultimate values of [[epsilon].sub.max], and therefore should be treated as the worst case scenario.

After accepting the Averhin's assumption and Knothe's theory formulas, Budryk uses the determination of (2) to specify the formula for the standard value of the B factor as follows:

[formula not reproducible] (3)

where:

R - main influences range radius acc. to Knothe's definition (1953a).

The value of the influence range radius R is determined by following formula (4):

R = H * cot[beta] (4)

where:

H - depth of mining operation,

[beta] - angle of main influences acc. to Knothe's definition.

According to Budryk, the standard value of the B factor corresponds with the range of maximum mining influence on the surface.

Both Knothe and Litwiniszyn in their theories assumed (following Avershin) that between the U horizontal displacement vector and the T tilt vector, there is proportionality:

U = -B * T. (5)

Formula (5) shows that the relation between the horizontal deformation tensor [epsilon] and the vertical curvature tensor K is as follows:

[epsilon] = -B * K. (6)

The most important indicator of deformation, from the point of view of threat to civil structures, is horizontal deformation (i.a. Knothe, 1953b; Tajdus and Tajdus, 2015; Sroka et al., 2015; Malinowska and Hejmanowski, 2016).

For today's conditions of mining operation and various application of the Knothe's theory, formula (6) is more general than (1).

Mining operations at very significant depths cause the mined out fields to be vastly geometrically different from the infinite semi-plane (i.a. Blachowski and Ellefmo, 2012). This prevents the application of formula (1) to forecasting calculations.

The results of the analyses of subsidence and horizontal displacement measurements for the points on the surface above the fields of cavern in a salt rockmass used to store liquid and gaseous fuels and for the fluidized bed extraction demonstrate that the value of the horizontal deformation factor B in these cases is significantly higher than the standard value provided by Budryk (1953). For instance, the analysis of vertical and horizontal displacements over the Etzel caverns made by Quasnitza (1988) in his doctoral dissertation shows that the value of this factor resulting from identification is:

B = 1.0 * R [+ or -] 0.1 * R. (7)

This proves that the most probable value of the B factor in this case is equal to the main effects range radius acc. to the Knothe's definition--i.e. it is some 2.5 times higher than the standard value provided by Budryk.

A comprehensive description of the formation of the B factor value studies to date can be found in the works of (Knothe and Sroka, 2010) as well as (Tajdus and Tajdus, 2015) and (Sroka et al., 2015). They are mostly applicable to the dependencies described by the following formulas:

B = k * H, or (8)

B = [lambda] * R. (9)

It should be noted that formula (9) is preferred in the literature. For hard coal and copper ore mining, the forecast values of [lambda] are lower than the standard value provided by Budryk ([lambda] = 0.4) with an average of 0.32 (Popiolek and Ostrowski, 1978).

SOLVED ISSUES

The studies of Dzegniuk (1970) and Sroka (1973) show that the sum of the values of horizontal deformations of the main components of the deformation tensor is 5/10 times larger than the value of vertical deformation.

[[epsilon].sub.xx] + [[epsilon].sub.yy] + [mu] * [[epsilon].sub.zz] = 0 (10)

where:

[[epsilon].sub.xx], [[epsilon].sub.yy]- deformation in the horizontal plane,

[[epsilon].sub.yy]- deformation in the vertical plane,

5.0 < [mu] < 10.0.

According to Sroka (1973), the value of [mu] for the surface horizon, when assuming the soil medium defined as a linear-elastic medium, amounts to:

[mu] = [1 - [nu]]/[nu] (11)

where:

[nu] - Poisson's ratio.

According to Forster (1996) for the small and slowly increasing deformations of the soil medium, the value of the Poisson's ratio is within the range of 0.10 / 0.20 while in the case of secondary deformations with high values and large increase rate --from 0.20 to 0.40.

The comparison of the aforementioned deformation values shows that the deformations caused by mining operations should be categorised as small. The observations of Dzegniuk (1970) and Sroka (1973) regarding the [mu] factor value we can infer that the corresponding range of values of the Poisson's ratio [nu] is:

0.09 < [nu] < 0.17

and corresponds rather well with the scope provided by Forster.

According to Drzezla (1989), for the rock mass with the tan[beta] parameter value from 1.5 to 2.5, the value of the Poisson's ratio is within the range from 0.096 to 0.150. This range is fully compliant with the findings of Dzegniuk (1970) and Sroka (1973).

By analysing the current condition of the Budryk-Knothe theory, Pielok and Sroka (1980/81), under the assumption of a standard horizontal deformation factor value, obtained the following dependency for the surface horizon:

[formula not reproducible] (12)

where:

n - so-called influence area factor in the rock mass.

The n factor is a parameter of an equation describing the changes in the value of the main mining operation influence range radius in the rock mass (formula 13).

[formula not reproducible] (13)

where:

R(z) - main influences range radius in the rock mass,

R - main influences range radius on the surface,

H - depth of mining operation,

Z - vertical distance from the roof of the working.

The values of the n parameter published by various authors are provided in Table 1.

It should be noted that the values of this parameter determined by in situ measurements are practically within the limits:

0.47 < n < 0.55.

This is applicable to the studies of: Sroka and Bartosik-Sroka (1974), Drzezla (1979, 1989) as well as Preusse (1990), who made the analysis of vertical movements in the rock mass for the mines in the Ruhr Valley.

By transforming formula (12), for any value of the B factor, we receive:

[formula not reproducible] (14)

For the assumption B = [lambda] * R and by assuming an equality between formulas (11) and (14), we receive:

[formula not reproducible] (15)

and, ultimately:

[formula not reproducible] (16)

The quality of formula (16) is also confirmed by previous partial results included in the studies by Drzezla (e.g. 1975, 1978, 1979, 1989). Drzezla (1975) provided approximate solutions of the elasticity theory as applied to the calculations of the area surface and rock mass deformation factors. In this study, he assumed that the vertical movement formula is in line with the solution provided by Knothe.

The result shows that the value of the n factor describing the changes in the range of main influences in the rocks mass is constant and equal to n = 0.665.

Drzezla (1989) asserts that the ratio S between the value of horizontal deformations/movement calculated in accordance with his solution and the values calculated according to the Budryk-Knothe theory is:

[formula not reproducible] (17)

which leads to the value of [lambda] described with formula (18).

[formula not reproducible] (18)

Comparing formula (16) to formula (18), it should be noted that they are identical in terms of quality. The difference between them is only that the formula based on Drzezla's assertions is based on a constant value of n = 0.665.

VERIFICATION OF THE ANALYTICAL MODEL ON THE BASIS OF HITHERTO EXPERIENCE

These considerations show that the value of [lambda] factor depends on the following:

* average value of the Poisson's ratio [nu] for the near-surface layer of soil,

* main mining influence area factor in the rock mass n, and

* angle of main influences [beta].

It seems obvious that the use of v and n will permit better estimation of [lambda] factor. On the other hand, however, it enables the calculation of [lambda] values, which are close to reality.

The obtained equation in the form of the general formula for determining the value of parameter B (formulas 9 and 16) was confronted with the results of the identification of this parameter value based on in situ measurements for two cases. The first concerns surface deformation in case of Etzel cavern and the second deformation will be by the exploitation of hard coal deposits in the Upper Silesian Basin.

For the characteristic data for the Etzel caverns field:

* Poisson's ratio: [nu] = 0.1,

* main mining influence area factor in the rock mass: n = 0.5,

* angle of main influences: [beta] = 34.8[degrees],

when using formula (16), we receive the [lambda] factor value of 1.03.

This result is consistent with the calculations made by Quasnitza (1988) and more recently by Sodmann (verbal information 2016--[lambda] = 1.06).

When assuming the following values for the Upper Silesian Coal Basin:

* [lambda] = 0.11,

* n = 0.5,

* tan[beta] = 2.0,

the resulting value of [lambda]= 0.32 corresponds well with the observed actual conditions (i.a. Popiolek and Ostrowski, 1978; Knothe and Sroka, 2010).

CONCLUSION

The considerations of this article enable us to provide a general formula for determining the value of the horizontal deformation factor B for various applications of the Knothe's theory. This applies to the calculations of surface deformations not only in the case of hard coal deposits mining, but also for borehole salt mining, deformations caused by the convergence of caverns used to store liquid and gaseous fuels or the deformations related to the fluidized bed extraction. According to formula (16) formulated based on the Knothe's theory, the value of the [lambda] factor (for B = [lambda] * R) depends on the following values:

* angle of main influences [beta],

* Poisson's ratio for the near-surface layer of soil [nu],

* main influence area factor in the rock mass n.

REFERENCES

Avershin, S.G.: 1948, Mining subsidence engineering. Ugletechizdat, Moskva, (in Russian).

Blachowski, J. and Ellefmo, S.: 2012, Numerical modeling of rock mass deformation in sublevel caving system. Acta Geodyn. Geomater., 9, 3 (167), 379-388.

Budryk, W.: 1953, Determining of the horizontal deformations values of the terrain. Archiwum Gornictwa i Hutnictwa, Warszawa, PWN, 1, 1, 63-74, (in Polish).

Budryk, W. and Knothe, S.: 1956, Rules for the classification of the Upper Silesian Industrial District areas due to the possibility of their development. Warszawa, 1-23, (in Polish).

Drzezla, B.: 1975, Approximate solution of elasticity theory equations for use in rock mechanics. Archiwum Gornictwa, 20, 2, Krakow, (in Polish).

Drzezla, B.: 1978, Solution of a certain spatial linear elasticity theory for use in predicting rock deformation due to mining exploitation, together with software. Politechnika Slaska, Zeszyty Naukowe, 588, Gliwice, 113 pp., (in Polish).

Drzezla, B.: 1979, Variability of main influences range in rock mass. Przeglad Gorniczy, 10, (in Polish).

Drzezla, B.: 1989, Description of software used in rock deformation prediction under mining exploitation, Current state of the software. Zeszyty Naukowe Politechniki Slaskiej, Gornictwo, 165, Gliwice, (in Polish).

Dzegniuk, B.: 1970, Experimental determination of the relationship between horizontal and vertical deformations in the rock mass. Zeszyty Naukowe Akademii Gorniczo-Hutniczej im. Stanislawa Staszica, Geodezja, 263, 17, Krakow, (in Polish).

Dzegniuk, B., Niedojadlo, Z. and Sroka, A.: 2003, Basics of dimensioning and exploitation in the safety pillars. Materialy Szkoly Eksploatacji Podziemnej 2003. Szczyrk, 17-21 lutego 2003. Wydawnictwo IGSMiE PAN, Nr. 58, Krakow, (in Polish).

Forster, W.: 1996, Mechanical properties of the loose rocks. Teubner Studienbucher. Bauwesen. Verlag B.G. Teubner Stuttgart, Leipzig, (in German).

Knothe, S.: 1953a, The profile equation of the finally formed subsidence trough. Archiwum Gornictwa i Hutnictwa, I, 1, PWN, Warszawa, 22-38, (in Polish).

Knothe, S.: 1953b, Project of mining areas classification from the point of usefulness for construction purposes view. Przeglad Geodezyjny, 9, 4, 97-99, (in Polish).

Knothe, S. and Sroka, A.: 2010, Stochastic assessment of extraction impact on building objects in the process of mining exploitation planning. III Konferencja Naukowo-Szkoleniowa "Bezpieczenstwo i ochrona obiektow budowlanych na terenach gorniczych", Ustron--Zawodzie, 4-10 pazdziernika 2010. Prace Naukowe GIG, Gornictwo i Srodowisko, 4, 1, 155-173, (in Polish).

Litwiniszyn, J. 1953: Differential equation of rock displacement. Archiwum Gornictwa i Hutnictwa, 1, 1, 9-21, (in Polish).

Malinowska, A. A. and Hejmanowski, R.: 2016, The impact of deep underground coal mining on Earth fissure occurrence. Acta Geodyn. Geomater., 13, No. 4 (184), 321-330. DOI: 10.13168/AGG.2016.0014

Pielok, J. and Sroka, A.: 1980-81, New concepts for surface and rock mass protection. Komisja Ochrony Terenow Gorniczych PAN, Oddzial Katowice, Wydzial Gorniczy Politechniki Slaskiej w Gliwicach, Metody wyznaczania i eksploatacji filarow ochronnych, Katowice- Gliwice, (in Polish).

Popiolek, E. and Ostrowski, J.: 1978, Dependence between slopes and horizontal displacements of terrain in finally formed subsidence trough. Ochrona Terenow Gorniczych, 46, (in Polish).

Preusse, A.: 1990, Mine surveyor analysis and prognosis of the vertical stress of shaft in the area of underground coal mining influence. TU Clausthal, Dissertation, (in German).

Quasnitza, H.: 1988, Calibration strategy of rock mass movements and rock elements prediction. Dissertation, Technische Universitat Clausthal, (in German).

Salustowicz, A.: 1953, Profile of subsidence trough as deflection of the layer on elastic subsoil. Archiwum Gornictwa i Hutnictwa, Warszawa, PWN, 1, 1, 39-50, (in Polish).

Sroka, A.: 1973, Relationships between the components of the deformation state on the surface. Rudy i Metale Niezelazne, 18, 12, Katowice, (in Polish).

Sroka, A. and Bartosik-Sroka, T.: 1974, Varibility of T. Kochmanski and T. Knothe theory in the rock mass. Rudy i Metale Niezelazne, 19, 7, Katowice, (in Polish).

Sroka, A., Preusse, A. and Kateloe, H.-J.: 2005, Basics on the dimensioning and the extraction of shaft safety zones. 24th International Conference on Ground Control in Mining, August 2-4, 2005, Morgantown, WV, USA.

Sroka, A., Knothe, S., Tajdus, K. and Misa, R.: 2015, Point movement trace vs. the range of mining exploitation effects in the rock mass. Archives of Mining Sciences, 60, 4, 921-929. DOI: 10.1515/amsc-2015-0060

Tajdus, K. and Tajdus, A.: 2015, Some considerations on horizontal displacement and horizontal displacement coefficient B. Studia Geotechnica et Mechanica, 37, 4, 75-82. DOI: 10.1515/sgem-2015-0047

Anton SROKA (*), Rafal MISA and Krzysztof TAJDUS

Strata Mechanics Research Institute of the Polish Academy of Sciences, ul. Reymonta 27, 30-059, Krakow, Poland

(*) Corresponding author's e-mail: sroka@img-pan.krakow.pl

ARTICLE INFO

Article history:

Accepted 29 November 2017

Available online 21 December 2017

Keywords:

Horizontal deformation

Horizontal displacement

Mining operation

Avershin's hypothesis

DOI: 10.13168/AGG.2017.0030
```Table 1 Value of the n parameter according to different authors
(Dzegniuk et al., 2003; Sroka et al., 2005).

Author          Year  Value

Budryk          1953  n = [square root of (2[pi])] * tan[beta]
Mohr            1958  n = 0.65
Krzyszton       1965  n = 1.0
Drzezla         1972  n = 0.525
Sroka,
Bartosik-Sroka  1974  n = 0.50
Drzezla         1975  n = 0.665
Gromysz         1977  n = 0.61
Drzezla         1979  0.47 [less than or equal to]
n [less than or equal to] 0.49
Kowalski        1984  0.48 [less than or equal to]
n [less than or equal to] 0.66
Zych            1985  n = 0.55
Drzezla         1989  0.50 [less than or equal to]
n [less than or equal to] 0.70
Preusse         1990  n = 0.54
```
Title Annotation: Printer friendly Cite/link Email Feedback ORIGINAL PAPER Sroka, Anton; Misa, Rafal; Tajdus, Krzysztof Acta Geodynamica et Geromaterialia Report Jan 1, 2018 2843 A METHOD OF FORECASTING EARTHQUAKE BY OPTICAL ASTRONOMIC TIMELATITUDE RESIDUALS. CLUSTERING OF 2008 WEST BOHEMIAN SEISMIC SWARM MECHANISMS AND ITS POSSIBLE INTERPRETATION. Coal industry Coal mining Geodynamics Mathematical models