_ ПРОБЛЕМЫ ЯДЕРНОЙ, РАДИАЦИОННОЙ _
И ЭКОЛОГИЧЕСКОЙ БЕЗОПАСНОСТИ -
УДК 631.6
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ТЕХНОЛОГИИ ОБЕСПЕЧЕНИЯ РАДИАЦИОННОЙ БЕЗОПАСНОСТИ ПРИ ПОЛИВЕ ПОЧВЫ МАЛЫМИ ДИФФЕРЕНЦИРОВАННЫМИ НОРМАМИ В САНИТАРНО-ЗАЩИТНОЙ ЗОНЕ АЭС
© 2015 г. В.И. Веприков*, Ю.В. Веприков**, Ю.В. Веприкова**
* Волгодонский инженерно-технический институт - филиал Национального исследовательского ядерного университета «МИФИ», Волгодонск, Ростовская обл.
** Институт сферы обслуживания и предпринимательства филиала «Донской государственный технический университет» в г. Шахты, Ростовская обл.
Моделируется инновационная технология полива дифференцированными малыми нормами, основанная на точном контроле влажности почвы и исключающая глубинную фильтрацию загрязнённой воды. Технология позволяет исключить деградацию почвы, отличается простотой и экономичностью. В работе исследовалась зависимость между скоростью перемещения фронта увлажнения и уменьшением плотности почвы с учётом величины межполивных интервалов. Приведённые экспериментальные результаты показывают принципиальную возможность улучшения плотности почвы применением технологии полива дифференцированными малыми нормами.
Ключевые слова: математическое моделирование, фильтрация воды, полив малыми нормами, плотность почвы, деградация почвы.
Поступила в редакцию 15.09.2015 г.
ПЛАНИРОВАНИЕ ЭКСПЕРИМЕНТА
Математическая модель технологии полива дифференцированными малыми нормами решает задачу прогнозирования глубины фильтрации радиационно загрязнённой воды для разных объёмов полива, что имеет важнейшее экологическое значение для безопасности человека. В работе исследовалась зависимость между скоростью перемещения фронта увлажнения и уменьшением плотности почвы с учётом величины межполивных интервалов.
Большинство исследований проводят для установления с помощью эксперимента функциональных или статистических связей между несколькими величинами или для решения экстремальных задач [1]. Классический метод постановки эксперимента предусматривает фиксирование на принятых уровнях всех переменных факторов, кроме одного, значения которого определенным образом изменяют в области его существования [2]. Этот метод составляет основу однофакторного эксперимента. При однофакторном эксперименте, варьируя один фактор и стабилизируя все прочие на выбранных уровнях, находят зависимость исследуемой величины только от одного фактора. Производя большое число однофакторных экспериментов при изучении многофакторной системы, получают зависимости, представленные графиками, имеющие иллюстративный характер [3]. Найденные таким образом частные зависимости невозможно объединить в одну общую.
©Издательство Национального исследовательского ядерного университета «МИФИ», 2015
Использование однофакторного эксперимента для всестороннего исследования многофакторного процесса требует постановки очень большого числа опытов. Для их выполнения в ряде случаев необходимо значительное время, в течение которого влияние неконтролируемых факторов на результаты опытов может существенно измениться. По этой причине данные большого числа опытов оказываются несопоставимыми. Отсюда следует, что результаты однофакторных экспериментов, полученные при исследовании многофакторных систем, часто малопригодны для практического использования [4].
Для решения интерполяционной задачи состоящей в построении формулы для предсказаний значений изучаемого параметра необходимо иметь математическую модель исследуемого объекта. Модель объекта получают, используя результаты опытов. При исследовании многофакторного процесса постановка всех возможных опытов для получения математической модели связана с огромной трудоемкостью эксперимента, так как число всех возможных опытов очень велико. Задача планирования эксперимента состоит в установлении минимально необходимого числа опытов и условий их проведения, в выборе методов математической обработки результатов опытов и в принятии решений [5]. Планирование эксперимента значительно сокращает число опытов, необходимых для получения модели процесса.
Применяемый в работе метод Бокса—Уилсона предусматривает проведение опытов небольшими сериями. В каждой серии одновременно варьируют все факторы по определенным правилам. Опыты проводят так, чтобы после математической обработки результатов предыдущей серии можно было спланировать следующую серию опытов. При планировании эксперимента цель исследования должна быть четко сформулирована и должна иметь количественную оценку. Характеристику цели, заданную количественно, называют параметром оптимизации. Параметр оптимизации является реакцией, или откликом, на воздействие факторов, определяющих поведение процесса. Результаты эксперимента используют для получения математической модели исследуемого процесса. Математическая модель — система математических соотношений, описывающих изучаемый процесс или явление. При планировании эксперимента под математической моделью часто понимают уравнение, связывающее параметр оптимизации с факторами. Такое уравнение называют также функцией отклика [6].
В общем случае объект исследования можно представить в виде структурной схемы, показанной на рисунке 1.
w
1
w
2
w„
Х1 X
2
X
к
► У1 ►У2
►У
р
Рис. 1. - Схема представления объекта исследования
Представление объекта в виде такой схемы основано на принципе «черного ящика». Имеем следующие группы параметров:
- управляющие (входные) Хь которые называются факторами;
- выходные параметры у!, которые называются параметрами состояния;
- "" - возмущающие воздействия.
Предполагается, что возмущающие воздействия не поддаются контролю и либо являются случайными, либо меняются во времени.
Управляющие параметры xi представляют собой независимые переменные, которые можно изменять для управления выходными параметрами. При поливе почвы малыми дифференцированными нормами управляющим параметром являлась норма полива. Выходным параметром являлась скорость движения фронта увлажнения.
1 ОБЩЕЕ РЕШЕНИЕ УРАВНЕНИЯ ГОЛОВАНОВА, ОПРЕДЕЛЯЮЩЕГО ЗАВИСИМОСТЬ НАПОРА ПОЧВЕННОЙ ВЛАГИ ОТ ВРЕМЕНИ И ГЛУБИНЫ ЗАЛЕГАНИЯ ПЛАСТА
Зависимость напора почвенной влаги от глубины и с течением времени описывается уравнением Голованова:
дЯ __д_ дН
дt дг В дг
Са- — _ — ^в— ) - Е (ю, г, t), (1)
где Сю - коэффициент влагоёмкости;
кв - коэффициент влагопроводности почвы, равен коэффициенту пропорциональности между скоростью потока воды и градиентом сил (изменяется от серии к серии); Н - напор почвенной влаги - изменяется - увеличивается со временем,
уменьшается по вертикальной координате; кв - коэффициент влагопроводности почвы,
кв _ £ф • ((ю-юм )/(Р -юм ))5,
где £ф - коэффициент фильтрации рассматриваемого слоя; ю - объёмная влажность почвы; юм - максимальная гигроскопичность почвы; Р - объёмная пористость почвы; е - интенсивность отбора влаги корнями растений. Коэффициент влагопроводности почвы кв зависит от коэффициента фильтрации £ф , который зависит от влажности почвы.
В эксперименте скорость фильтрации изменяется от слоя к слою, а так как кв - коэффициент влагопроводности почвы связан с коэффициентом фильтрации
£
ф зависимостью:
5 д£ф д£В £В _ £ф ■ ((Ю-Юм)/(Р-Юм))5, то следует ^
Уравнение (1) является линейным дифференциальным уравнением в частных производных второго порядка (неоднородным). Для упрощения модели величины
коэффициентов Сю,£В, Е(ю, г, t) будем полагать условно-постоянными. Тогда
уравнение (1) можно переписать в виде:
_ ЭН , Э 2Н _
Со"ЭН = кв (эЕ(0(2)
Применим к уравнению (2) метод, подобный методу разделения переменных. Сделаем в уравнении (2) следующую замену:
Н (г, г ) = А( г) - Б(г) + F (г); (3) После замены (3) уравнение (2) можно записать в форме: ^ л( \ 4Б . оМй2А , й2^
с»- ) —=к в ■ б( ^+к в ш -Е. <4)
Функцию F = ^ (г) будем определять из условия:
кв - Е = 0 (5)
й )
Общее решения уравнения (5) согласно [1] имеет вид:
Е
F (г ) = — г2 + С г + С2, (6)
2к в
где с1, с2 - произвольные постоянные.
Теперь функцию Н = Н (г, г) из (3) можем представить в форме:
Е
Н (г, г) = А( г) - Б (г) + — г2 + С г + С2. (7)
2к в
Подставим зависимость в виде (7) в уравнение (4), после упрощений получим: _ Л йБ , „м й2А
С»'А(г) 'йБ = к в'Б(г • (8)
Разделим обе части уравнения (8) на произведение А(г) Б(г), после упрощений уравнение (8) преобразуется к виду:
_ йБ / к йА
С°~йй/ кв (йг)\ /Б (г )= /А(г)
(9)
Левая часть уравнения (9) зависит только от переменной г, а правая - только от г , поэтому левая, и правая части уравнения (9) являются постоянной величиной,
которую обозначим через Я2 > 0 , то есть:
с йБ / кв
0 йг/, л = в (йг) V. , = Я2 (10)
/Б(г) /А(г) Я . (10)
Тогда уравнение (9) распадается на пару обыкновенных дифференциальных уравнений:
, й2А о2 л
квТ^2 = ЯА(г). (11) (йг)
С» йБ = Я2Б(г). (12) йг
Уравнение (11) является обыкновенным линейным дифференциальным
уравнением второго порядка с постоянными коэффициентами (однородным), соответственно его общее решение в соответствии с [1] имеет вид:
—j=z --т= z
A(z) = Axe^B + A2e ^
А2 е 1 (13)
где А1, А2 - произвольные постоянные.
Уравнение (12) является обыкновенным линейным дифференциальным уравнением первого порядка с постоянными коэффициентами (однородным), соответственно его общее решение в соответствии с [1] имеет вид:
B(t) = B1eCa , (14)
где B1 — произвольная постоянная.
Пусть теперь значение константы в уравнении (10) будет отрицательным, то есть:
d2 A
C
k
= 'В (dz)2
= —A2.
(15)
B(t )= /A(z)
В этом случае уравнение (9) распадается на следующую пару обыкновенных дифференциальных уравнений:
d2 A
" (16)
k
= —Л2 A(z).
5 (dz )2
с. f = ~Л2 B(t).
dt
(17)
Также как и ранее, уравнение (16) является обыкновенным линейным дифференциальным уравнением второго порядка с постоянными коэффициентами (однородным), соответственно его общее решение в соответствии с [1] для мнимых корней характеристического уравнения примет вид:
A(z) = A sin
Л
\M kв J
+ A2 cos
Л
VA/ k В J
(18)
где А1, А2 - произвольные постоянные.
Общее решение уравнения (17) в соответствии с [1] имеет вид:
Л2 —t
B(t) = Bxe C* ,
(19)
Решения в виде (13), (14) и (17), (18) являются базовыми для построения общего решения уравнения (2).
Таким образом, из (7) с учетом формул (13), (14) и (17), (18) следует, что решения уравнения (2) могут быть представлены в следующих видах:
H (z, t ) =
AekB + Ae Jk
л \
z
■ e
H 2 (z, t ) =
Aj sin
Л
■yf^B
+ A2 cos
Л
E
2k В
W
z2 + Cj z + C2.
(20)
-J^B
Л
--1
У J
E
■ e C + 2Tz2 + C1 z + C2. (21)
2k В
Общее решение уравнения (2) содержит функции вида (13) и (17) при различных
z
z
Л
t
a +
z
значениях параметра Я , например, при целых значениях. Так как уравнение (2) линейное, то его общее решение содержит и линейные комбинации решений вида:
Я N
Де^ + А2 е
еСш
а также решений вида: с
Л1 вШ
' я л
г
V
лД~В
+ Л2 008
У
Я
\\
уу
--г
- ст
Из проведенных рассуждений следует, что для случая положительного значения
параметра из уравнения (10), т.е. при значениях Я > 0 общее решение уравнения (2) можно записать в форме:
„2 Г
Н (г, г ) = X е
п=0
п
Л1,„е
п \ г —г
кв + Л е Чкв + Л2,пе
Е 2 +-г2 + С г + С2.
2к
(22)
В
V у
В случае отрицательного значения параметра из уравнения (10), т.е. при
значениях — Я < 0 общее решение уравнения (2) можно записать в виде линейной комбинации тригонометрических функций по переменной г :
/ л / \Л
Н 2 (г, г ) = X е
п2
~--г
С,
-ш
п=0
Л1,п вш
п
■у[кв
+ Л2,п оов
п
у
Е2
+ -—- г + С1г + С 2.(23) 2к в
Общие решения уравнения (2) в формах (22) и (23) являются базовыми для решения граничных задач по определению зависимости напора почвенной влаги от времени и глубины залегания пласта. Они будут использованы в последующем изложении для решения граничных задач, возникающих при определении давления почвенной влаги на разных почвенных горизонтах при поливе участков малыми нормами.
2 ГРАНИЧНЫЕ ЗАДАЧИ ПО ОПРЕДЕЛЕНИЮ НАПОРА ПОЧВЕННОЙ ВЛАГИ ОТ ВРЕМЕНИ И ГЛУБИНЫ ЗАЛЕГАНИЯ ПЛАСТА
Для нахождения решения граничной задачи по определению напора почвенной влаги необходимо описать все ее начальные и граничные условия.
Описание начальных и граничных условий задачи проведем в соответствии с физическим процессом фильтрации водного потока через грунт при поливе участка малыми нормами.
Как показали эксперименты, на каждом уровне залегания пласта напор почвенной влаги убывает с течением времени, то есть является монотонно убывающей функцией по времени. Из этого условия следует, что функции вида (22) не удовлетворяют
физическим условиям задачи фильтрации воды через почву, так как функции Н1 (г, г )
содержат множители е ш , которые являются монотонно возрастающими по времени. Условию убывания напора почвенной влаги с течением времени удовлетворяет вторая группа решений в форме (23), соответственно в дальнейшем решения граничной задачи
Я
г
г
к
в
ш
2
п
по определению напора будем искать в классе функций вида (23).
Рассмотрим граничную задачу по определению напора почвенной влаги с течением времени в первом поливе. Экспериментальные данные по распределению давлений на глубинах залегания пласта 15 см, 30 см, 45 см, 60 см представлены в таблице 1.
Таблица 1. - Временная зависимость показаний тензиометров для почвенных _горизонтов в первом поливе_
Глубина, см Первый полив нормой 20 мм
15 30 45 60
Т, мин. Р, сентибар
0 68 37 28 21
50 53 37 28 20
78 43 37 27 20
108 35 37,1 27 20
142 27 37 27 20
208 12,6 37,2 27 20
238 10 37,1 27 20
298 10 37,2 27 19
358 10 37,2 27 19
При фиксированном значении переменной г , т.е. для фиксированного уровня залегания пласта решения граничной задачи в форме (23) допускают функции в виде:
Н (г о, г ) = в
е
■—(г-а)
са — 1
-—(г—а) е Са +1
+ 7,
(24)
Л2
Параметры а, определялись из условия минимизации рассогласования
С а
теоретической кривой и экспериментальной по методу наименьших квадратов. Таким образом, на уровне залегания пласта г = 15 см напор почвенной влаги в первом поливе определяется по формуле:
Н1 (15, г ) = 57
е
-0,014 (г—10)
1
— 0,014 (г—10)
+ 65.
(25)
е ----------+1
Сравнение теоретической кривой и экспериментальных данных напора почвенной влаги из таблицы 1 на уровне залегания пласта г = 15 см приводится на рисунке 2. Р,сентибар
Н, см
Рис. 2. - Сравнение теоретической кривой напора почвенной влаги и точек эксперимента на уровне г = 15 см в первом поливе
Все расчеты проводились в среде прикладных математических программ Maple 9.5. Расчеты показывают, что рассогласование между теоретической кривой и экспериментальными данными не превышает 5%.
Рассмотрим граничную задачу по определению напора почвенной влаги с течением времени во втором поливе. Экспериментальные данные по распределению давлений на глубинах залегания пласта 15 см, 30 см, 45 см, 60 см представлены в таблице 2.
Таблица 2. - Временная зависимость показаний тензиометров для почвенных
горизонтов во втором поливе
Глубина, см Второй полив нормой 20 мм
15 30 45 60
Т, мин. Р, сентибар
0 46 42 46 30
25 44 42 46 30
72 4 42 46 30
138 4 29 46 30
199 4 17 46 30
234 4 11 46 30
299 4 10 45,8 30
422 4 11 46 29,8
Зависимость напора почвенной влаги с течением времени при фиксированном уровне залегания пласта будем снова искать в форме (24). При этом параметры
л ^
а будем определять из условия минимизации рассогласования
теоретической кривой и точек эксперимента по методу наименьших квадратов.
На уровне залегания пласта г = 15 см напор почвенной влаги во втором поливе предлагается определять по формуле:
е-0,14(г-50) _ 1
Н2 (15, * )= 21 ^ - 0,14(*-50) + 1 + 24'7- (26)
Номер индекса у функции напора в формуле (26) означает номер полива. Сравнение теоретической кривой и экспериментальных данных напора почвенной влаги из таблицы 1 на уровне залегания пласта г =
15 см приводится на рисунке 3.
Р, сентибар
30 25 20 15
0 50 100 150 200
х
Н, см
Рис. 3. - Сравнение теоретической кривой напора почвенной влаги и точек эксперимента на уровне г = 15 см во втором поливе
Как показывают расчеты в среде Maple 9.5, рассогласование между теоретической кривой и экспериментальными данными на уровне z = 15 см во втором поливе не превышает 2%.
На уровне залегания пласта z = 30 см напор почвенной влаги во втором поливе
предлагается определять по формуле:
- 0,037 (t-155 ) _ I
H2(30,') = 16,1 g_ 0,037 (.-,55 ) + 1 + 26,4- (27)
R
Параметры а, также определялись из условия минимизации
рассогласования теоретической кривой и точек эксперимента по методу наименьших квадратов.
Сравнение теоретической кривой и экспериментальных данных напора почвенной влаги из таблицы 2 на уровне залегания пласта z =
30 см во втором поливе
приводится на рисунке 4. Р, сентибар
Н, см
Рис. 4. - Сравнение теоретической кривой напора почвенной влаги и точек эксперимента на уровне Z = 30 см во втором поливе
Как показывают расчеты в среде Maple 9.5, рассогласование между теоретической кривой и экспериментальными данными на уровне Z = 30 см во втором поливе не превышает 3%.
Рассмотрим граничную задачу по определению напора почвенной влаги с течением времени в третьем поливе. Экспериментальные данные по распределению давлений на глубинах залегания пласта 15 см, 30 см, 45 см, 60 см представлены в таблице 3.
Зависимость напора почвенной влаги с течением времени при фиксированном уровне залегания пласта будем снова искать в форме (24). При этом параметры
л ^
а, P, будем, как и ранее, определять из условия минимизации рассогласования
теоретических кривых и соответствующих им точек эксперимента по методу наименьших квадратов.
На уровне залегания пласта Z = 15 см напор почвенной влаги в третьем поливе предлагается определять по формуле:
e - 0,7 (t-21) _ 1
H3(15't)= 21'3 e_0,7(t-21)+ 1 + 23 (28)
Таблица 3. - Временная зависимость показаний тензиометров для почвенных горизонтов в третьем поливе
Глубина, см Третий полив нормой 43 мм
15 30 45 60
Т, мин. Р, сентибар
0 44 33 35 36
19 41 33 35 36
23 5 33 35 36
38 2 33 35 36
53 2 31 35 36
56 2 6 35 36
68 2 2 32 36
78 1 2 2 35
97 2 2 2 22
108 2 2 2 12
109 2 2 2 12
110 2 2 2 12
111 2 2 2 10,2
112 2 2 2 8,4
113 2 2 2 8,8
114 2 2 2 7,5
115 2 2 2 6,4
116 2 2 2 6
117 2 2 2 6
118 2 2 2 6
119 2 2 2 6
204 2 2 2 6
215 2 2 2 6
216 2 2 2 6
Сравнение теоретической кривой и экспериментальных данных напора почвенной влаги из таблицы 3 на уровне залегания пласта z = 15 см приводится на рисунке 5.
Как показывают расчеты в среде Maple 9.5, рассогласование между теоретической кривой и экспериментальными данными на уровне z = 15 см во втором поливе не превышает 4%.
Р, сентибар -------,
2520 15-
0 50 100 150 200
х
га,
Рис. 5. - Сравнение теоретической кривой напора почвенной влаги и точек эксперимента на уровне Z = 15 См в третьем поливе
На уровне залегания пласта г = 30 см напор почвенной влаги в третьем поливе предлагается определять по формуле:
е"(?"55) -1
Нз (30, г) = 15,5 -М5) | 1 +17,5. (29)
+1
Сравнение теоретической кривой и экспериментальных данных напора почвенной влаги из таблицы 3 на уровне залегания пласта г = 30 см в третьем поливе приводится на рисунке 6.
Р, сентибар
Н, см
Рис. 6. - Сравнение теоретической кривой напора почвенной влаги и точек эксперимента на уровне г = 30 см в третьем поливе
Как показывают расчеты в среде Maple 9.5, рассогласование между теоретической кривой и экспериментальными данными на уровне z = 30 см в третьем поливе не превышает 3%.
На уровне залегания пласта z = 45 см напор почвенной влаги в третьем поливе предлагается определять по формуле:
-(t-70,5) _ 1
H3 (45, t) = 16,5 _(t-70,5) . 1 +18,5. (30)
е
+1
Сравнение теоретической кривой и экспериментальных данных напора почвенной влаги из таблицы 3 на уровне залегания пласта г = 45 см в третьем поливе приводится на рисунке 7.
Р,сентибар
Н, см
Рис. 7. - Сравнение теоретической кривой напора почвенной влаги и точек эксперимента на уровне г = 45 см в третьем поливе
Как показывают расчеты в среде Maple 9.5, рассогласование между теоретической кривой и экспериментальными данными на уровне z = 45 см в третьем поливе не превышает 2%.
На уровне залегания пласта z = 60 см напор почвенной влаги в третьем поливе предлагается определять по формуле:
-0,18 (t-98) _ i
Н3 (45, t) = 15 ^-0,18(t-98) + 1 + 21 (31)
Сравнение теоретической кривой и экспериментальных данных напора почвенной влаги из таблицы 3 на уровне залегания пласта z = 60 см в третьем поливе приводится на рисунке 8.
Р, сентибар
35 30 25 20 15 10
0 ' ' ' ' 50 ' ' ' 1□□ ' ' 150 ' ' 200
х Н, СМ
Рис. 8. - Сравнение теоретической кривой напора почвенной влаги и точек эксперимента на уровне Z = 60 см в третьем поливе
Как показывают расчеты в среде Maple 9.5, рассогласование между теоретической кривой и экспериментальными данными на уровне z = 60 см в третьем поливе не превышает 4%.
Рассмотрим граничную задачу по определению напора почвенной влаги с течением времени в четвертом поливе. Экспериментальные данные по распределению давлений на глубинах залегания пласта 15 см, 30 см, 45 см, 60 см представлены в таблице 4.
Зависимость напора почвенной влаги с течением времени при фиксированном
Л *
уровне залегания пласта будем определять в форме (24). Параметры СХ-,
будем находить из условия минимизации рассогласования теоретических кривых и соответствующих им точек эксперимента по методу наименьших квадратов.
На уровне залегания пласта z = 15 см напор почвенной влаги в четвертом поливе предлагается определять по формуле:
_0,5(t-34) _ 1
Н4 (15, t) = 24,5 ^_0,5(t-34) + 1 + 26,5 (32)
Сравнение теоретической кривой и экспериментальных данных напора почвенной влаги из таблицы 4 на уровне залегания пласта z =
15 см приводится на рисунке 9.
Таблица 4. - Временная зависимость показаний тензиометров для почвенных горизонтов в четвёртом поливе
Глубина, см Четвёртый полив нормой 20 мм
15 30 45 60
Т, мин. Р, сентибар
0 51 49 49 14
17 51 49 25 14
26 50 49 22 14
29 48 49 22 14
30 44 49 22 14
33 33 49 22 14
35 21 49 22 14
37 10 49 22 14
39 6 49 22 14
41 4 49 22 14
51 2 45 22 14
53 2 8 22 14
55 2 6 22 14
92 2 6 22 14
107 2 6 22 14
119 2 6 22 14
137 2 6 22 14
172 2 6 22 14
202 2 6 22 14
242 2 6 22 14
277 2 6 22 14
Р, сентибар
Н, см
Рис. 9. - Сравнение теоретической кривой напора почвенной влаги и точек эксперимента на уровне Z = 15 См в четвертом поливе
Как показывают расчеты в среде Maple 9.5, рассогласование между теоретической кривой и экспериментальными данными на уровне z = 15 см в четвертом поливе не превышает 2%.
На уровне залегания пласта z = 30 см напор почвенной влаги в четвертом поливе предлагается определять по формуле:
H 4 (30, t ) = 21,5
e
-2 (t-52)
1
e
-2 (t-52)
+ 1
+ 27,5.
(33)
Сравнение теоретической кривой и экспериментальных данных напора почвенной влаги из таблицы 4 на уровне залегания пласта г = 30 см в третьем поливе приводится на рисунке 10.
Р, сентибар
100
150
200
га, СМ
Рис. 10. - Сравнение теоретической кривой напора почвенной влаги и точек эксперимента на уровне z = 30 см в четвертом поливе
Как показывают расчеты в среде Maple 9.5, рассогласование между теоретической кривой и экспериментальными данными на уровне z = 30 см в третьем поливе не превышает 2%.
На уровне залегания пласта z = 45 см напор почвенной влаги в четвертом поливе предлагается определять по формуле:
,-0,4(t-11,5) - 1
+35,7. (34)
H 3 (45, t ) = 13,7
e
-0,4 (t-11,5)
+ 1
Сравнение теоретической кривой и экспериментальных данных напора почвенной влаги из таблицы 4 на уровне залегания пласта г = 45 см в четвёртом поливе приводится на рисунке 11.
Р, сентибар
Н, см
Рис. 11. - Сравнение теоретической кривой напора почвенной влаги и точек эксперимента на уровне Z = 45 см в четвёртом поливе
Как показывают расчеты в среде Maple 9.5, рассогласование между теоретической кривой и экспериментальными данными на уровне Z = 45 см в четвёртом поливе не превышает 2%.
ВЫВОДЫ
1) Расчёты в среде Maple 9.5, дают рассогласование между теоретической кривой и экспериментальными данными в пределах 2-4%, в зависимости от уровня почвенного горизонта, т.е. модель в пределах допустимых погрешностей адекватно отображает процесс движения горизонта увлажнения и может быть использована для прогнозирования результатов реальных экспериментов.
2) Установленная функциональная зависимость между напором почвенной влаги и глубиной фронта увлажнения полностью решает задачу прогнозирования глубинны фильтрации радиационно загрязнённой воды для разных объёмов полива.
СПИСОК ЛИТЕРАТУРЫ
1. Новик, Ф.С. Математические методы планирования экспериментов в металловедении. Раздел 3. Выбор параметров оптимизации и факторов [Текст] / Ф.С. Новик. - М.: МИСИС, 1971. - 106 с.
2. Адлер, Ю.П. и др. Планирование эксперимента при поиске оптимальных условий [Текст] / Ю.П. Адлер, Е.В. Маркова, Ю.В. Грановский. - М.: Наука, 1976. - 279 с.
3. Львовский, Е.Н. Статистические методы построения эмпирических формул [Текст] / Е.Н. Львовский. - М.: Высшая школа, 1992. - 224 с.
4. Веприков, В.И. Математическое моделирование процессов электрокристаллизации из водного раствора электролита, под действием постоянного магнитного поля. Юбилейный сборник научных трудов ЮРГУЭС [Текст] / В.И. Веприков // Известия высших учебных заведений. Северо-Кавказский регион. Технические науки. Приложение №6. - 2004. - С. 153-157.
5. Мелик-Гайказаян, И.В. и др. Методология моделирования нелинейной динамики сложных систем [Текст] / И.В. Мелик-Гайказаян, М.В. Мелик-Гайказян, В.Ф. Тарасенко. - М.: Физматлит, 2001. -272 с.
6. Чуличков, А.И. Математические модели нелинейной динамики [Текст] / А.И. Чуличков. - М.: Физматлит, 2000. - 296 с.
REFERENCES
[1] Novik F.S. Matematicheskie metody planirovaniya eksperimentov v metallovedenii [Mathematical methods of planning of experiments in metallurgical science]. Razdel 3 [Part 3]. Vybor parametrov optimizacii i faktorov [Choice of parameters of optimization and factors]. М. Pub. MISIS [MISiS], 1971, 106 p. (in Russian)
[2] Adler Yu.P., Markova E.V., Granovskij Yu.V. Planirovanie eksperimenta pri poiske optimalnyh uslovij [Planning of experiment by search of optimum conditions]. M. Pub. Nauka [Science], 1976, 279 p. (in Russian)
[3] Lvovskij E.N. Statisticheskie metody postroeniya empiricheskih formul [Statistical methods of creation of empirical formulas]. M. Pub. Vysshaya shkola [Higher Education], 1992, ISBN 5-06001264-6, 224 p. (in Russian)
[4] Veprikov V.I. Matematicheskoe modelirovanie processov еlektrokristallizacii iz vodnogo rastvora elektrolita, pod dejstviem postoyannogo magnitnogo polya. Yubilejnyj sbornik nauchnyx trudov YuRGUES [Mathematical modeling of processes of electrocrystallization of water solution of electrolyte, under the influence of a constant magnetic field. Anniversary collection of scientific works of SRSSTU]. Izvestiya vysshix uchebnyx zavedenij [News of higher educational institutions]. Severo-Kavkazskij region [North Caucasus region] Texnicheskie nauki [Technical science] Prilozhenie №6 [Appendix № 6]. 2004, ISSN 0321-2653, p. 153-157 (in Russian)
[5] Melik-Gajkazayan I.V., Melik-Gajkazyan M.V., Tarasenko V.F. Metodologiya modelirovaniya nelinejnoj dinamiki slozhnyh system [Methodology of modeling of nonlinear dynamics of difficult systems]. M. Pub. Fizmatlit [Phismatlit], 2001, ISBN 5-9221-0087-4, 272 p. (in Russian)
[6] Chulichkov A.I. Matematicheskie modeli nelinejnoj dinamiki [Mathematical models of nonlinear dynamics]. M. Pub. Fizmatlit [Phismatlit]. 2000, ISBN 5-922100-83-1, 296 p. (in Russian)
Mathematical Model of Radiation Safety Technology when Watering the Soil by Small Differentiated Standards in the NPP Sanitary Protection Zone
V.I. Veprikov*1, Yu.V. Veprikov**2, Yu.V. Veprikova***2
1 Volgodonsk Engineering Technical Institute the branch of National Research Nuclear University «MEPhI»,
73/94 Lenin St., Volgodonsk, Rostov region, Russia 347360
2 Institute of services industry and business branch of "Don State Technical University", 147, Shevchenko St.,
Shakhty, Rostov region 346500 * e-mail: vv_@list.ru ; ** e-mail: veprikow@mail.ru ; *** e-mail: sw_@list.ru
Abstract - The innovative technology of watering by the differentiated small standards based on exact control of humidity of the soil and excluding a deep filtration of the polluted water is modelled. The technology allows excludes degradation of the soil, differs in simplicity and profitability. In work dependence between the speed of movement of the front of moistening and reduction of density of the soil taking into account the size of interirrigation intervals was investigated. The given experimental results show basic possibility of improvement of density of the soil by application of technology of watering by the differentiated small standards.
Keywords: mathematical modeling, water filtration, watering with small standards, the soil density, soil degradation.