CC BY
11
Geology • Geography
Dnipro university bulletin
Journal home page: geology-dnu-dp.ua
ISSN 2313-2159 (print) ISSN 2409-9864(online)
Dniprop. Univer.bulletin. Geology, geography., 25(2), 123-128.
Tatyina P. Mokritskaya,
Ksenia O. Samoylych, Andrii O. Stanchev
doi: 10.15421/111727
Dniprop. Univer. bulletin, Geology, geography., 25(2), 123-128.
Towards the issue of the methods of calculating fractal dimension of the disperse pul-verescent-clayey soil structure
T. P. Mokritska, K. A. Samoylych, A. O. Stanchev
Oles Honchar Dnipro National University, Dnipro, Ukraine, e-mail:[email protected]
Received 04 October 2017
Received in revised form 30 October 2017
Accepted 19 November 2017
Abstract.Our study on the structure of pulverescent-clayey loess deposits was conducted at different levels of organization with the purpose of clarifying the characteristics of their deformational behaviour.In our description of the structure and microstructure, the theory of fractals is used in the study ofmacroporosity, inter-aggregate porosity and interparticle porosity of loess loams and sand loams. The calculations of indicators of fractal dimension structure can be made using analysis of electronic images of microstructure, granulometry and dispersity of soils, experimental calculations of penetrability and porosity. Comparison of the results of quantitative fractal dimension indicators obtained using different methods is of interest at the stage of choosing a more efficient and meaningful method. The article analyses the results of calculating the values of the fractal dimension of particle distribution function, which were obtained using different methods. We conducted analysis of discretion of differently formed stratigraphic-genetic horizons. The images were analyzed in the "JmageL" programme, the calculation of the Hurst exponent and fractal dimension of the particle distribution function was made using "Fraktan" software and EXELL electronic tables. The study analyzedimages on different scales. The results of the calculations of the particle sizes were considered as a series of fractal dimension series sufficient for calculation, Hurst exponent, correlational dimension and entropy of length. The results of granulometric examination (pipette method, three methods for preparing the samples) were analyzed using classified data on the fractions. The following fractions were considered: the fractions with particles less than 0.001 mm, 0.002 mm, 0.005 mm, 0.01 mm, 0.05 mm, 0.1 mm, 0.25 mm, 0.5 mm, 1 mm, 2 mm, 5 mm, 10 mm. The values of fractal dimension of particle distribution at the level of inter-aggregate and interpar-ticle porosity were different. When the disperse method for preparing the sample was used, the fractal dimension of particle (fraction) distribution obtains higher values compared to the aggregate and semi-disperse methods. We proved the relation between the values of fractal dimension and the genesis, discretion and method for preparing the soil samples. The dimension values of loess sediments were lower compared to paleosolregardless of the method for the preparation. Comparability of the assessments of fractal dimension of distribution of the solid phase of different topological dimension (flat sections and volumetric particles) was observed during aggregate and semi-disperse methods of preparation.
Keywords:loess, porosity, fractal dimension,particle distribution function.
До питання про методи оцшки фрактальноТ p03MipH0CTi структури дисперсного пилу-вато-глинистого грунту.
Т.П. Мокрицкая, К.О. Самойлич, А. O. Станчев
Днiпровський нацюнальний утверситет шет Олеся Гончара, Днiпро, Украша, e-mail: mokritska@i. ua.
Резюме. Вивчення структури пилувате-глинистих лесовидних BiAKTa^B виконуеться на pi3Hnx рiвнях оргашзаци з метою уточнения особливостей !х деформацшно! поведшки. Застосування теорп фракталгв до опису структури i мжроструктури виконуеться при вивченш макропориста, мiжaгpегaшоI i мiжчaстковоI пористост лесовидних суглинюв i сутсюв. Юльюст оцшки значень фрактально! pозмipностi структури можна отримати при ана^ електронних зтмюв мжроструктури, грану-лометрп i дисперсном! грунпв, експериментальних оцшок проникност i пористом! Зютавлення результата кшьюсних оцшок фрактально! pозмipностi, отриманих piзними методами, становить штерес на етaпi вибору бшьш результативного i обгрунтованого методу. У данш pоботi пpоaнaлiзовaнi результати кшьюсно! оцiнки значень фрактально! pозмipностi функ-цiI pозподiлу часток, яю отpимaнi piзними методами. Виконано aнaлiз дискретност piзних за умовами походженнястратиг-рафо-генетичних гоpизонтiв.Обpобкa знiмкiв виконувалася за допомогою програми "JmageL", розрахунок показника Херсту i фрактально! pозмipностi функцп pозподiлу часток виконувався за допомогою ПО "Fraktan" i електронних таблиць EXELL. Дослiджено знiмки в piзних масштабах. Поодинокш запис результата тдрахунку pозмipiв частинок розглядався як ряд дос-татньою для розрахунку фрактально! pозмipностi ряду, показника Херсту, кореляцшно! pозмipностi i
ешропй довжини. Обробка результатiв гранулометричного аналiзу (пiпеточнiй метод, три способи тдготовки зразкгв) ви-конувалася по згрупованим даними - фракцiям. враховувалося змiст частинок розмiрами менше 0,001мм, 0,002мм, 0,005мм, 0,01мм, 0,05мм, 0,1мм, 0,25 мм, 0,5 мм, 1 мм, 2 мм, 5 мм, 10 мм. Значення фрактально! розмiрностi розподшу часток на рiвнi мiжагрегатноl i мiжчастковоl пористостi вiдрiзняються. Фрактальна розмiрнiсть частинок (фракцiй) при дисперсному способi тдготовщ зразка приймае бшьш високi значення, нiж за агрегатно! або толдисперсношдготовкою. Шдт-верджено залежтсть значень фрактально! розмiрностi ввд генезису, дискретност i способу тдготовки зразка Грунту. Значення розмiрностi приймають меншi значення в лесових вiдкладах, нiж в палеогрунтах, незалежно вiд способу пiдготовки. Порiвняннiсть оцшок фрактально! розмiрностi розподiлу твердо! фази рiзноl тополопчно! розмiрностi (плоских перетитв i об'емних часток) спостерiгаеться при агрегатному i полiдисперсному способах тдготовки.
Ключовi слова: лес, пористкть, фрактальнарозмiрнiсть, функцiя розподыу часток.
The problem statement and state of study.The
usage of the theory of fractals in studying natural objects, phenomena and processes has rapidly increased in recent yearsPuetz, S.J., Borchardt G. (2015); Cheng Wen-Chiehet al, (2017); X. Du • et al, (2017).
The study ofpulverescent-clayey soils as a natural multifractal system is a matter of great of interest (Zhan Ce Wang et al, 2009). The description of soil structure from the perspective of the theory of fractals, search for relation between the genesis and peculiarities of the structure, attempts of defining patterns in deformation and decrease in the strength of soils as natural fractals is studied by scientists of various countries (Russell AR, 2014; Munoz-Castelblanco JA and all, 2012; Li P.et al, 2016).
The classical methods of predicting the deformation of soil environments are based on using the model of continuous medium, and less often on the discrete medium model; the processes are considered reversible (Jiang M. and all, 2017). A principal difference of the non-classic method is defining patterns of deformation in soil as a natural fractal, a dynamic evolving system (Ikeda, Kiyohiroet al, 2008), capable of chaotic behaviour. Methods.We studied the peculiarities of the structure of loess sediments in Prychornomorsko-Dofinovski, Vitachevski, Udaiski, Dniprovski deposits from the perspective of the theory of fractals. The samples were collected from natural outcrops around the city of Dnipro: in the slopes of the Yevpatoriiska, Tonnelna, and Sazhenka ravines and from the RybalskyQuarry.
The basin of the Yevpatoriiska ravine includes Aeolian-diluvialloam sands of the Tiasmyn-sky horizon (vd PIIts), Dniprovski Aeolian-diluvial loess loam sands (vd PIIdn) and Zavadovskieluvial loams (e PIIzv). In the basin of the Tonnelna ravine, we collected samples of the Prychornomorsko-Dofinovsky Aeolian horizon (vd PIII pc + df), Buhskieluvial-diluvial loam sands (vd P IIIbg), DniprovskiAeolian-diluvial loam sands (vd P IIdn) (Fig. 4.1 a). The samples of Buhski Aeolian-diluvial loam sands (vd P IIIbg), eluvial loam sands of the Vitachevsky horizon (e PIIIvt), Aeolian-
diluvial loam sands of the Dniprovsky horizon (vd PIIdn) were also collected in the basin of the Sa-zhevka ravine.
We calculated the full complex of indicators of physical properties, focusing mainly on the characteristics of dispersity. The calculation of fractal dimension of the surfaces, images, ramifieldsys-tems is performed in materials science, medicine, geology using the results of analyzing a series of data of sufficient length. The digital model of an electronic image of the soil structure can be converted into a series of values of area, diameter, perimeter and contour (of particles and pores), using the corresponding programme, in particular, "Jma-geL"
The developed sequence is a random series of the contours'occurrence on the surface, i.e. a characteristic of microstructure. The analysis of the images made on different scales allows analysis of different elements of the soil structure - aggregates, particles and pores. The sufficient length of the series allows the calculation of the fractal dimension of a row, Hurst exponent, correlational dimension and entropy ["Fraktan" software].
Due to anisotropy and complex form of particles, aggregates and pores, the calculations of fractal dimension of random section do not allow accurate characterization of the soil structure, but allow making selective estimations of distribution of the number of particles of particular size.
Russell A R (2011)has suggested a method of calculating fractal dimension of function of particles distribution by using the results of standard granulometric analysis. This case does not include the calculation of plainly visible particles.Loess pulverescent-clayey soils are characterized by predisposition of the microaggregatesto the break-up under the change of external conditions, which is the manifestation of their degradation during the process of anthropization. This method (Riashenko T H, 2010) allows one to obtain more accurate calculations of content of free particles and particles of different fractions assembled in aggregates than when using standard granulometric analysis. The main results of the study.We performed the calculations of fractal dimension of function of par-
ticle distribution during the analysis of plane images and on the basis of results of analysis of samples' microstructure with calculation of composition of aggregates and particles of loess from the Dniprovsky horizon. Electronic images were made using an electronic microscope at 80, 200 and 500 x
magnification (Fig. 1). The analysis of images was made using the programme "JmageL". The calculations of the Hurst exponent, correlational dimension, fractal dimension of function of distribution of particles and pores were performed using "Fractan" programme.
During the preparation of the initial image for the analysis of particles, we used operations of calibration, scaling, conversion of the image. We performed the calculations of radius and diameters of particles, number of particles included at a certain interval. The values of intervals were set in correspondence with current classification of the dispersed rocks; the content of particles with sizes less than 0.001 mm, 0.002 mm, 0.005 mm, 0.01 mm, 0.05 mm, 0.1 mm, 0.25 mm, 0.5 mm, 1 mm, 2 mm, 5 mm, 10 mm were considered. If the particles of these sizes were not visible in the image, the fraction was omitted.
The second stage was calculation of the fractal dimension of the random ordered sequence of the image contours (form 1) and the random un-
-8
ordered sequence (form. 2). In the first case, the data was grouped according to the fractions; the second case used no grouping methods. Charts of dependencies between the total number of particles in the fraction or the number of particles and the diameter (fig. 2, 3) were developed. The fractal dimension of the grouped data Ds was calculated using the formula:
Ds = 2-k
(form. 1)
WhereDs - fractal dimension of the distribution of the particles' contours; 2 - topological dimension of the image; K - the slope of the lineLnNR = f{LnDR).
6 5 4 3 2 1
y = 1,008x + 6,46 R2 = 0,957
0 -1 -2
Fig. 2.The chart of relation between the total number of particles NR and the diameter DR In the form of LnNR = f(LnDR) (Udaiski,
vdPIIud).
Notes to Fig. 2:
1. Nr - the total number of particles: DR - the largest diameter
2. R2 - determination coefficient
3. x- experimental values;--linear trend.
The fractal dimension of random sequence was calculated using the Hurst exponent H (form. 2).
D = 2 - H (form. 2).
The values of fractal dimension D and Ds are significantly different, which indicates the difference in
the calculations of fractal dimension of functions of fraction and certain particle distribution (Table 1). The values of dimension change in relation to the increase, i.e. the type of porosity (macroporosity or inter-particle porosity).
Selected results of calculations of values of fractal dimension of the structure of the Dniprovsky horizon Table 1.
Magnification The fractal dimension of the function of particle (fraction) size distribution The fractal dimension of the particle contour sequence
k Ds D H
80 0.045 1.955 1.119 0.881
200 0.034 1.966 1.143 0.857
800 0.037 1.963 1.145 0.855
Notes to Table 1.
K - slope coefficient of the straight line; Ds- fractal dimension of the sequence of grouped data; D - fractal dimension of contour sequence; H - the Hurst exponent
The calculations of the dimension of the function of particle distribution using the results of granulometric analysis with three methods of sample preparation indicate the presence of a dependency between the dispersity (aggregativity) and the dimension.
The sample preparation method also affects the dispersity of soil (Fig.3). The disperse method of preparation achieves the highest decomposition of the aggregates and higher number of fine fractions compared to the samples prepared using aggregate or semi-disperse method.
edPIIIdf
0,001 0,01 0,1 1
vdPIIIud
0,001 0,01 0,1 1
Fig. 3. Integral curves of granulometric composition of DofinovskyedP///d/, UdaiskyvdPnlud and Dniprovskyvd PH dn horizons. Notes to Fig. 3, 4.
,3, - disperse method of sample preparation; A- aggregate method of sample preparation; n- semi-disperse method of sample preparation.
The fractal dimension of fraction distribution function was calculated using the method described above (see form. 1, Fig. 2), the value of topological dimension was considered as equal to 3. The analysis of the results showed that using the disperse method of sample preparations, the values
of the fractal dimension of the particles (fractions) are higher compared to the aggregate and semiaggregate methods (Fig. 4).
This indicates the natural relation between the dispersity and the fractal dimension of the particle distribution function of loess soil.
3 2,5 2 1,5 1 0,5 0
edPIIIdf
ed PIIIvt
vPIIIud
Fig. 4. The fractal dimension of DofinovskyedPIIIdf, UdaiskyvdPIIIud and Dniprovskyvd PI1dn horizons
Also we calculated the fractal dimension of the distribution of particles in the samples with altered condition: compressed and filtrated (Table 2). The fractal dimension of the distribution of particle sizes is related not only to the dispersity, but also to the hydrophility. In the horizons of fossil soils (water-saturated condition) the values of dimension were higher, and in the subaerial loess sediments were lower compared to the dimension typical for natural conditions. A sharp rise in the content of free clayey particles as a result of compression in
the water-saturated condition is typical for the horizons of fossil soils. The loess sediments were observed to have change in the content of pul-verescent and larger particles and aggregates.
As a result, the value of the fractal dimension were lower after compression and moisturizing of the soil. Usage of magnetic stirrer in the preparation of the sample for the granulometric analysis leads to increase in the value of dimension compared to the values after using mechanical stirring
Table 2. The fractal dimension of the soil particles distribution function (change in the condition).
Horizon S Experimental conditions
V Z M F
Prychornomorsko-Dofinovsky(vd,ePIUpc+df) 2.7 2.67 2.63 2.41 -
Udaisky horizon (vd Pm ud) 2.25 - - 2.61 -
Tiasmynsky horizon (vdPnts) 2.45 - - - 2.37
Dniprovsky horizon (vd Pu dn) 2.57 - - 2.65 -
2.59 2.73 - 2.67 -
2.76 - - - 2.45
Kaidaksky horizon (e Pukd), 2.78 2.39 2.32 - -
Zavadovskyhorizon (ePuzv) 2.62 - - - 2.5
Notes to table 2:
S - standard preparation of the sample forgranulometric analysis; V - the condition after the compression in water-saturated condition; Z - the sample tested for the subsidence using the method of a single curve; M - aggregate preparation (magnetic stirrer); F - the sample affected by the filtration.
Results and discussion.The method of sample preparation forgranulometric analysis affects the measurement of the content of fine fractions, fractal dimension of particles distribution function. Disregarding the methods of preparation, the values of dimension are lower in loess sediments compared to the paleosols.
Comparability of the values of fractal dimension of distribution of solid phase of different topological dimension (flat sections and volumetric particles) was observed during the aggregate and semi-disperse methods of preparation.
The external effect on the loess soil causes a change in the microaggregate structure, which was observed in the changes of fractal dimension of the soil particle distribution function.
References
Puetz, S. J., Borchardt G. 2015. Quasi-periodic fractal patterns in geomagnetic reversals, geological activity, and astronomical events. Chaos, Solitons and Fractals 81: 246-270. Jiang Mingjing, Sima Jun, Cui Yujun, Hu Haijun, Zhou Chuangbing, Lei Huayang 2017. Experimental
127
Investigation of the Deformation Characteristics of Natural Loess under the Stress Paths in Shield Tunnel ExcavationInt. J. Geomech., 17(9): 04017079.
Cheng Wen-Chieh, Cui Qing-Long, Shui-Long Shen Jack, Arulrajah Arul and Yuan Da-Jun 2017. Fractal Prediction of Grouting Volume for Treating Karst Caverns along a Shield Tunneling Alignment. Appl. Sci. 7, 652; doi: 10.3390/app7070652
X. Du P. Zhang, L. Jin R. Zhang 2017. A mesoscopic model of sand-gravel stratum in Beijing based on fractal theory. Chinese Journal of Rock Mechanics and Engineering. Vol. 36 Issue (2): 437-445.
Zhan Ce Wang, Zhang Yang Liu, Shi-min Fan 2017. Geometric and fractal analysis of dynamic cracking patterns subjected to wetting-drying cycles. -Soil and Tillage Research, 1-13.
Ikeda, Kiyohiro and Yamakawa, Yuki and Desrues, Jacques and Murota, Kazuo. Bifurcations to Diversify 2008. Geometrical Patterns of Shear Bands on Granular Material. Phys. Rev. Lett.
100(19): 198001-198005.
DOI: 10.1103/PhysRevLett. 100.198001.
Russell A. R. 2011. A compression line for soils with evolving particle and pore size distributions due to particle crushing. Geot Let 1:5-9.
Russell A. R. 2014. How water retention in fractal soils depends on particle and pore sizes, shapes, volumes and surface areas. Ge'ot 64 (5):379-390.
Russell A. R., Buzzi O. 2012. A fractal basis for soil-water characteristics curves with hydraulic hysteresis. Ge'ot 62 (3):269-274.
Riashenko T. H. 2010. Regional soil science (Eastern Siberia). SD RAS, Irkutsk.
Munoz-Castelblanco JA, Pereira JM, Delage P and Cui YJ 2012. The water retention properties of natural unsaturated loess from Northern France. Geot 62 (2): 95-106.
Li P, Vanapalli S, Li T 2016. Review of collapse triggering mechanism of collapsible soils due to wetting. J Rock Mech and GeotEng 8:256 - 274 https://doi.org/10.1016/jjrmge.2015.12.002.
