Определение влажностных параметров внутриконтурного капельно увлажняемого почвенного пространства

ТЕХНИЧЕСКИЕ НАУКИ

удк 631.423.2:631.674.6 DOI 10.32786/2071-9485-2018-04-45

ОПРЕДЕЛЕНИЕ ВЛАЖНОСТНЫХ ПАРАМЕТРОВ ВНУТРИКОНТУРНОГО КАПЕЛЬНО УВЛАЖНЯЕМОГО ПОЧВЕННОГО ПРОСТРАНСТВА

DETERMINATION OF MOISTURE PARAMETERS OF INCONTOUR DRIP MOISTURED SOIL SPACE

С. М. Васильев, доктор технических наук, доцент В. Н. Шкура, кандидат технических наук, профессор А. С. Штанько, кандидат технических наук

S. M. Vasil'ev, V. N. Shkura, A. S. Shtanko

ФГБНУ «Российский научно-исследовательский институт проблем мелиорации» FSBS «Russian Scientific-Research Institute of Land Improvement Problems»

Целью исследования является разработка методики прогнозирования распределения уровней и величины влажности почвы в зонах промачивания (увлажнения) подкапельного почвенного пространства при капельном поливе растений. При капельном орошении в почвогрунтовой толще формируется контур увлажнения, геометрические и влажностные параметры которого являются критериальными показателями, определяющими качество капельного полива и качество условий для роста и развития корневой системы сельскохозяйственных культур. Несмотря на наличие ряда исследований в этой области, общепринятая методика определения влажностных параметров контуров капельного увлажнения почвы отсутствует, что предопределяет актуальность научной работы. В результате исследований предложена методика графоаналитического определения геометрических размеров и влажностных параметров контуров и построения очертания контуров капельного увлажнения почв для широкого диапазона почвенных, фитопочвенных и технологических условий капельного полива растений, позволяющая спрогнозировать и оценить геометрические размеры контуров и характер распределения уровней влажности почвы во внутриконтурном почвенном пространстве с приемлемой для практического использования точностью. Методика базируется на материалах авторских экспериментальных исследований и полученных на их основе экспериментальных зависимостях по определению очертания, глубин и диаметров контуров увлажнения и координат линий влажности (изоплет) почвы в пределах увлажняемого внутриконтурного пространства для широкого спектра условий капельного полива.

The object of the research is to develop a methodology for predicting the distribu-tion of levels and the value of soil moisture in moistening zones of the subsurface soil space during drip irrigation of plants. The humidification contour where geometric and moisture parameters are criteria indicators determining the quality of drip irrigation and the conditions for the growth and development of agricultural crops root system is formed in the soil layer with drip irrigation. Despite the presence of a number of studies in this area, there is no generally accepted method for determining the moisture parameters of the contours of drip soil moistening, which predetermines the relevance of scientific work. As a result of the research a method for graphoanalytical determination of geometric dimen-sions and moisture parameters of contours and for constructing contours of drip irrigation for a wide range of soil, phytosoil and technological conditions for drip irrigation of plants is proposed. It allows predicting and estimating the geometrical sizes of contours and the distribution of soil moisture levels in the incontour soil space with acceptable accuracy for practical use. The methodology is based on the materials of the author's experimental stud-ies and the obtained experimental dependences on determining the shapes, depths and di-ameters of moisture contours and the coordinates of moisture lines (iso-pleths) of soil with-in the moistured incontour space for a wide range of drip irrigation conditions.

Ключевые слова: капельное орошение, капельный полив, контур капельного увлажнения, геометрия контура, влажность почвы, линии влажности, зона увлажнения почвы.

Key words: drip irrigation, drip watering, drip moistening contour, contour geometry, soil moisture, humidity lines, soil moistening zone.

Введение. При капельном поливе растений в подкапельном почвенном пространстве формируются локальные (единичные) контуры увлажнения (промачивания), характеризуемые геометрическими (линейными, площадными и объемными) размерами и влажностными параметрами в пределах увлажняемой почвогрунтовой зоны. Исследованию геометрических и влажностных параметров контуров капельного увлажнения почвы посвящены работы А. Д. Ахмедова, В. В. Бородычева, В. С. Бочарникова, Е. В. Мелиховой, А. С. Овчинникова, Е. А. Ходякова, М. Ю. Храброва, О. Е. Ясониди, М. Л. Ромащенко, Ю. С. Уржумовой и других специалистов в области капельного орошения [1-10]. Несмотря на значительный объем исследований, общепринятая методика определения геометрических и влажностных параметров контуров капельного увлажнения к настоящему времени отсутствует, а одна из попыток ее разработки и есть цель настоящей статьи.

Материалы и методы. Распределение влажности почвы во внутриконтурном почвенном пространстве зависит от ряда факторов капельного полива, в качестве которых рассматриваются нижеследующие.

1. Почвенные характеристики: содержание в промачиваемом почвенном пространстве глинистых частиц - Wr/4 , измеряемое в процентах от массы сухой почвы (% МСП); средняя плотность сложения почвы - /об, т/м и наименьшая влагоемкость почвы - WHB, % МСП.

2. Фитопочвенные параметры: заданная глубина увлажняемого слоя почвы -Л (глубина контура - Лкон ), м; дополивная влажность почвы - Дп , % МСП; постполивная влажность почвы - Д^п , % МСП.

3. Технологические параметры: поливная норма - ^пол, м3/кап. и продолжительность формирования контура увлажнения ?пол, ч.

Для достижения цели необходимо решить следующие задачи.

1. Установить наличие функциональной связи вида - Лкон = f (/об,

W/ч, WHB, , Д^п, ^пол ) и получить описывающую ее зависимость.

2. Для определенных значений уровней влажности почвы Д = (0,6; 0,65; 0,7; 0,75; 0,8; 0,85; 0,9; 0,95; 1,0; 1,05; 1,1; 1,15 и 1,2) • Дш (где Дш - влажность почвы, соответствующая наименьшей влагоемкости) во внутриконтурном почвенном пространстве установить глубинные и плановые координаты внутриконтурных изоплет и геометрические параметры разноизоплетных («разновлажностных») контуров и в частности значений (Лиз/п)д , (^из/п)д , где (Лиз/п)д - максимальная величина заглубленности

изоплеты с уровнем влажности Д, отмеряемая по оси капания от поверхности земли, м; (^из/п)д - максимальный диаметр контура, ограниченного изоплетой с уровнем влажности Д .

3. Определить вертикальные и горизонтальные координаты точек (глубину (Лиз/П)д и радиус изоплетного контура (гиз/п)д = 0,5 • (^из/п)д ) изоплеты с влажностью

почвы, соответствующей дополивному ее значению Д /п, и изоплет с заданными уров-

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА: НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

нями влажности почвы от Д = Дд/п до Д = 1,2 • Дш, ограничивающих «изоплетные» контуры влажности в границах общеконтурного почвенного пространства, и построить систему внутриконтурных линий влажности, характеризующих распределение влажности почвы в увлажняемом пространстве.

4. Построить формирующийся в заданных условиях капельного полива контур увлажнения почвы с внутриконтурными изоплетами, позволяющий определить его геометрические размеры и влажностные параметры и оценить условия формирования и функционирования корневой системы растений во внутриконтурном почвенном пространстве.

Результаты и обсуждение. Методика определения распределения влажности почвы в пределах контуров увлажнения, формирующихся в подкапельном почвенном пространстве при капельном поливе, заключается в проведении нижеследующих рас-четно-графических операций.

1. Определение расчетных значений факторов влияния (Жг/ч и ЖНВ в % МСП;

/об в т/м3; Дд/п и Дп/п в % МСП; ^ м и ^ м3/кап).

д/п /~п/п ' " — --у

2. Расчет глубины контура увлажнения (¿из/пД , м при Д = Д

д/п '

(Лиз/п )Дд/п

N

0,00196 • /0б •

Полученное по зависимости (1) значение (¿из/п)^ соответствует глубине контура увлажнения почвы и может быть записано как Лкон. В случае когда (¿из/п)д ^ ¿увл, изменяют предварительно принятое значение поливной нормы, добиваясь равенства значений (¿Из/п)дд/п = ¿Увл.

3. Определение величины максимальной заглубленности изоплеты (¿из/п)0 7 дв, м, с уровнем влажности, соответствующим Диз/п = 0,7 • Дш:

(Аиз/п ) Дд/п

из/п/ 0,7Днв

0,5 •

1,0 -1,1

0,1^гА

г в л Дд/п- -1,0

Д

0,7Днв

+1,0 - 1,20,1^НВ

Д

Дд/п- -1,0

Д0

0,7Днв

4. Определение значений заглубленности изоплет (¿из/п )д, м, с уровнем влажно

сти от Д = Дд/п до Д = 1,2 • Дш с шагом изменения Д - АД = (0,05 - 0,1) • Дш (т. е. Д = (0,6; 0,65; 0,7; 0,75; 0,8; 0,85; 0,9; 0,95; 1,0; 1,05; 1,1; 1,15 и 1,2) • Дш по нижеследу-

ющей зависимости:

(¿из/п)д = 0,5 •

1,0 -1,1'

0,шгЬ

Д

--1,0

+1,0 -1,2'

0Жн,

Д

-1,0

'

0,7Днв .

5. Для установленных значений по п. 4 заглубленности изоплет (¿из/п)д определяются соответствующие им значения максимальных диаметров (^из/п)д , м, и радиусов контуров (г^/п)^ , м, по зависимостям:

(<з/п)Д = 0,5• [(0,51 + 0,009-Жг/Ч) + (0,073 + 0,038)]• (^Л ,

(ги/п )д = 0,5 • (^п )Д .

318

0,333

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА: НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

6. Для каждого принятого к рассмотрению уровня влажности почвы Д выполняются расчеты координат точек, соответствующих изоплетам, ограничивающим изо-плетные контуры, с соблюдением нижеприведенной последовательности расчетно-графических операций.

6.1. По каждому изоплетному контуру с установленной по п. 4 глубиной (Ииз/п)д намечается 11 вертикальных сечений (вертикальных координат ограничивающей контур изоплеты (Иу , где ] - номер сечения) начиная с (\ = 0 для сечения 1 с шагом по

глубине ДИд = 0,1 • (Ииз/п)д в соответствии со схемой, приведенной на рисунке 1.

Данные разбивки профиля по изоплетным контурам с влажностью почвы Д на 11 горизонтальных сечений заносятся в таблицу 1.

Таблица 1 - Матрица вертикальных (глубинных) координат граничной линии изоплетных контуров влажности почвы

Показатель Значение показателя

№ сечения 1 2 3 4 5 6 7 8 9 10 11

И )Д /(Ииз/п)д 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0

% )Д , м (Мд. (И)Д (Из)д (Юд (И7)Д. (И8)Д. (И9 ) Д. (И10)Д. (И11)Д.

_4 4 ч 1 ч £ -« л о  и 3

/ Ч 2 ч \ " _\| V / 1 (

/ Ч 3 ч гл -л с -т я

ч 4 ч \ 1 о 4) II _ |П О II * I Г й ^ с * -г

ч 5 ч о II С» X . »-

ч 6 ч II л 1 ц

7 ч , в 1

ч 3 ч

\ ч

10 ч» у * 2 1

II. _

1 3

Рисунок 1 - Расчетная схема изоплетного контура капельного увлажнения почвы: 1 - капельница; 2 - оконтуривающая линия; 3 - ось капания; 4 - поверхность земли; (И^ (ИОд - расстояния от поверхности земли до характерных сечений, разделяющих контур по вертикали на 10 расчетных слоев, м; тн + тн - радиусы контура, соответствующие характерным сечениям; 1-11 - точки на оси капания, соответствующие характерным сечениям контура увлажнения почвы

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА: НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

6.2. Значения горизонтальных координат гй по каждому из сечений изоплетных контуров глубиной (Иу. )д определяются по зависимости:

Я = (Гиз/п)Д  [к0 + к1  [(И] )д /(Ииз/п)Д ] - к2  [(И] )д /(Ииз/п)Д ] + + К [(к] )д /(¿из/п )Д ] - [(к] )Д /(^из/п )Д}°]

где Г - радиус ограничивающей изоплетный контур линии, соответствующий сечению изо-плетного контура с глубиной И-, м; (гиз/п)д - максимальный радиус изоплетного контура с влажностью ограничивающей линии Д , определяемый по соотношению (5); к0 - экспериментально установленный свободный член полинома:

К = 0,5  (2 - 0,005  Жг/ч  г.об - 0,001  ЖоВ84 - 0,002  Ж^45);

К - коэффициент при члене полинома [(И])д /(Ииз/п)д ]:

К = 1,03/ К01,2;

к2 - коэффициент при члене полинома [(И )д /(Ииз/п)д. ] :

— к ^ к\ ;

к, к4 - эмпирические коэффициенты при членах полинома со значениями

[(И] )Д /(Ииз/п)Д]3 и [(И] )Д /(Ииз/п)Д ]0 соответственно:

к — к — 0,255  к00,01.

6.3. Расчетные значения координат (И и Г заносят в таблицу 2.

Таблица 2 - Матрица координат изоплеты, ограничивающей изоплетный контур влажности почвы

Показатель Значение показателя

№ сечения 1 2 3 4 5 6 7 8 9 10 11

и )д (И2)Д (И3)Д (И4)Д (и5)д (ил (И7)Д (ил (И9)Д ^Д ^Д

% ГЮ ГИ2 ГИ3 ГИ5 % ГИ9 ГИ00 ГИ11

6.4. По приведенным в таблице 2 данным и в соответствии со схемой по рисунку 1 строятся изоплетные («влажностные») контуры.

Пример расчета по методике определения влажностных параметров контуров капельного увлажнения почвы выполнен для условий капельного полива, проведенного Ю. С. Уржумовой [6]. Приведенная в [6] матрица влажности почвы в подкапельном пространстве получена при капельном поливе южного тяжелосуглинистого средне-мощного чернозема для нижеследующих почвенных, фитопочвенных и технологических условий: Жг/ч = 58,0 % МСП; ЖНВ = 25,2 % МСП; уоб = 1,4 т/м3; Дп/п = 0,95  ДНВ

= 0,9525,2 = 23,9 % МСП; Дд/п = 0,6  Днв = 0,6  25,2 = = 15,1 % МСП; Иувл = 0,83 м; Nпол = 37,3 л/кап.; дкаП = 3,2 л/ч; I_ = 11,7 ч.

***** ИЗВЕСТИЯ ***** № 4 (52), 2018

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА: НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

1. Расчетное значение поливной нормы для условий опыта:

^ = 0,00196 • ^ • С • (0,0765 • Ж™ + 0,0292 • Жнв)2 • (Д^ - Д*) =

= 0,00196 • 1,4 • 0,833 • (0,0765 • 58,00,6 + 0,02 92 • 25,2)2 х х (0,95 • 25,2 - 0,6 • 25,2) = 35,9 л/кап. Полученная расчетная величина Жпол меньше фактической на Д^ = 37,3 - 35,9 = 1,4 л/кап или на ам = 3,8 %, что допустимо.

2. По приведенной в [6] матрице влажности почвы в подкапельном пространстве построен контур увлажнения, приведенный на рисунке 2.

г, см 40 30 20 10 0 10 20 30 40 г, см

Рисунок 2 - Контур капельного увлажнения почвы, построенный по первичным опытным данным Ю. С. Уржумовой. Подписи на изоплетах - влажность в % от наименьшей влагоемкости почвы

3. Расчетное значение соотношения диаметра «из/п и глубины йш/п контуров, формирующихся в вышеуказанных условиях, составит:

«з/п /¿ш/п)расч = 0,5 [(0,51 + 0,009 • Ж/ч) + (0,073 + 0,038 • Жт)] = 0,5 х х [(0,51 + 0,009 • 58,0) + (0,073 + 0,038 • 25,2)] = 1,03.

Опытное (среднее по всем изоплетным контурам) значение соотношения

«з/п / \з/п )опыт = 1,05 и превышает расчетное на 2 %.

4. Соответствующие изоплетным контурам значения глубин (^из/п)расч и диаметров («из/п) ч контуров влажности в сопоставлении их с фактическими значениями этих параметров приведены в таблице 3.

Таблица 3 - Расчетные и опытные значения

Показатель Значение показателя

Дз/п , % Рнв 60 65 70 75 80 85 90 95 100

(^из/п)опыт , м 0,83 0,76 0,68 0,64 0,53 0,46 0,44 0,24 0,14

(«из/п )опыт, м 0,94 0,86 0,74 0,60 0,52 0,48 0,43 0,26 0,14

(^из/п)расч , м 0,83 0,75 0,67 0,60 0,52 0,44 0,36 0,28 0,20

(«из/п)расч , м 0,86 0,78 0,69 0,62 0,53 0,45 0,37 0,29 0,21

5. Значения коэффициентов к0, к , к2, к и к4 в зависимости (6): к0 — 0,5  [(1 - 0,005  Жг/Ч  Гоб) + (1 - 0,001  Ж0в84 - 0,002  ЖНВ45)] = 0,5

№ 4 (52) 2018

НВ 1,84

[(1 - 0,005  58,0  1,4) + (1 - 0,001  25,21,84 - 0,002  25,20,45)] = 0,603

к —1,03/к!'2 = 1,03/0,6031,2 = 1,889; к — к + к = 0,603 + 1,889 =

= 2,492; к — К — 0,255  к00'01 = 0,254.

6. Для заданных уровней влажности оконтуривающих линий Д = (0,6 - 1,0) Дш,

соответствующих им значений

и данных, приведенных в таблице 3, получены расчетные значения вертикальных и горизонтальных координат точек оконтуривающих линий по изоплетным контурам, позволяющие построить прогнозируемый для соответствующих опыту почвенных, фито-почвенных и технологических условий полива контур увлажнения подкапельного почвенного пространства, который приведен на рисунке 3.

Рисунок 3 - Расчетный контур капельного увлажнения подкапельного почвенного пространства

Заключение. Формы очертаний опытных и прогнозных изоплетных контуров увлажнения почвы являются аффинными. Параметры опытного и прогнозного контуров, ограниченных изоплетами с уровнями влажности, составляющими Диз/п = (0,6; 0,7;

0,8; 0,9; 1,0) ДНВ отличаются по глубине и радиусам в среднем не более чем на 12 %, что позволяет заключить о приемлемости методики и расчетных зависимостей для использования в практике проектирования и эксплуатации систем капельного орошения.

Библиографический список

1. Ахмедов, А. Д. Контуры увлажнения почвы при капельном орошении [Текст]/ А. Д. Ахмедов, Е. Ю. Галиуллина // Известия Нижневолжского агроуниверситетского комплекса: наука и высшее профессиональное образование. - 2012. - № 3. - С. 183-188.

2. Бородычев, В. В. Современные технологии капельного орошения овощных культур [Текст]/ В. В. Бородычев. - Волгоград, 2010. - 241 с.

X

X

3. Бочарников, В. С. Новые приемы возделывания овощных культур в системе водосбере-гающего орошения [Текст]/ В. С. Бочарников, М. П. Мещеряков // Роль науки в инновационном развитии АПК: материалы Всероссийской научно-практической конференции, посвященной 80-летию со дня рождения известного ученого, профессора А. П. Иофинова /Башкирский ГАУ. -Уфа, 2012. - 204 с.

4. Мелихова, Е. В. Моделирование и обоснование ресурсосберегающих параметров капельного орошения при возделывании корнеплодов [Текст]: монография / Е. В. Мелихова. -Волгоград: ФГБОУ ВО Волгоградский ГАУ, 2017. - 112 с.

5. Овчинников, А. С. Применение ресурсосберегающих способов полива при возделывании сельскохозяйственных культур [Текст]/ А. С. Овчинников, М. П. Мещеряков // Известия Нижневолжского агроуниверситетского комплекса: наука и высшее профессиональное образование. - 2007. - № 1. - С. 46-49.

6. Уржумова, Ю. С. Технологические и конструктивные элементы локального низконапорного орошения садов для условий южных черноземов Ростовской области [Текст]: дис. ... канд. техн. наук / Уржумова Юлия Сергеевна. - Новочеркасск, 2004. - 147 с.

7. Ходяков, Е. А. Режим орошения сельскохозяйственных культур при капельном и внутрипочвенном способах полива [Текст]: монография / Е. А. Ходяков. - Волгоград: ВГСХА, 2002. - 143 с.

8. Храбров, М. Ю. Оценка способов малообъемного орошения [Текст]/ М. Ю. Храбров // Вестник РАСХН. - 2007. - № 5. - С. 53-56.

9. Ясониди, О. Е. Капельное орошение [Текст]: монография / О. Е. Ясониди; НГМА. -Новочеркасск: Лик, 2011. - 322 с.

10. Системи краплинного зрошення: навчальний пошбник / М. Л. Ромащенко [и др.]; за ред. акад. УААН М. Л. Ромащенка. - Дшропетровськ: Оксамит-текст, 2007. - 175 с.

References

1. Ahmedov, A. D. Kontury uvlazhneniya pochvy pri kapel'nom oroshenii [Tekst]/ A. D. Ah-medov, E. Yu. Galiullina // Izvestiya Nizhnevolzhskogo agrouniversitetskogo kompleksa: nauka i vysshee professional'noe obrazovanie. - 2012. - № 3. - Р. 183-188.

2. Borodychev, V. V. Sovremennye tehnologii kapel'nogo orosheniya ovoschnyh kul'tur [Tekst]/ V. V. Borodychev. - Volgograd, 2010. - 241 р.

3. Bocharnikov, V. S. Novye priemy vozdelyvaniya ovoschnyh kul'tur v sisteme vodosbere-gayuschego orosheniya [Tekst]/ V. S. Bocharnikov, M. P. Mescheryakov // Rol' nauki v innovacionnom razvitii APK: materialy Vserossijskoj nauchno-prakticheskoj konferencii, posvyaschennoj 80-letiyu so dnya rozhdeniya izvestnogo uchenogo, professora A. P. Iofinova /Bashkirskij GAU. - Ufa, 2012. - 204 р.

4. Melihova, E. V. Modelirovanie i obosnovanie resursosberegayuschih parametrov kapel'no-go orosheniya pri vozdelyvanii korneplodov [Tekst]: monografiya / E. V. Melihova. - Volgograd: FGBOU VO Volgogradskij GAU, 2017. - 112 р.

5. Ovchinnikov, A. S. Primenenie resursosberegayuschih sposobov poliva pri vozdelyvanii sel'skohozyajstvennyh kul'tur [Tekst]/ A. S. Ovchinnikov, M. P. Mescheryakov // Izvestiya Nizhnevolzhskogo agrouniversitetskogo kompleksa: nauka i vysshee professional'noe obrazovanie. - 2007. -№ 1. - Р. 46-49.

6. Urzhumova, Yu. S. Tehnologicheskie i konstruktivnye ]lementy lokal'nogo nizkonapornogo orosheniya sadov dlya uslovij yuzhnyh chernozemov Rostovskoj oblasti [Tekst]: dis. ... kand. tehn. nauk / Urzhumova Yuliya Sergeevna. - Novocherkassk, 2004. - 147 р.

7. Hodyakov, E. A. Rezhim orosheniya sel'skohozyajstvennyh kul'tur pri kapel'nom i vnutripoch-vennom sposobah poliva [Tekst]: monografiya / E. A. Hodyakov. - Volgograd: VGSXA, 2002. - 143 р.

8. Hrabrov, M. Yu. Ocenka sposobov maloob'emnogo orosheniya [Tekst]/ M. Yu. Hrabrov // Vestnik RASXN. - 2007. - № 5. - Р. 53-56.

9. Yasonidi, O. E. Kapel'noe oroshenie [Tekst]: monografiya / O. E. Yasonidi; NGMA. - Novocherkassk: Lik, 2011. - 322 р.

10. Sistemi kraplinnogo zroshennya: navchal'nij posibnik / M. L. Romaschenko [i dr.]; za red. akad. UAAN M. L. Romaschenka. - Dniropetrovs'k: Oksamit-tekst, 2007. - 175 р.

E-mail: rosniipm@yandex.ru