CC BY
37
Том 150, кн. 4
Естественные науки
2008
УДК 631.43
ЭНЕРГЕТИЧЕСКАЯ ОЦЕНКА СОСТОЯНИЯ ПОЧВЫ И СТЕПЕНИ ВЛИЯНИЯ НА НЕЕ ПОЧВООБРАБАТЫВАЮЩЕЙ ТЕХНИКИ
В.В. Сироткин, В.М. Сироткин
Аннотация
Разработан метод оценки энергетического состояния почвы и степени механического воздействия на нее почвообрабатывающих машин на основе единого интегрального параметра, определяемого экспериментально-расчетно на основе измерения общепринятых величин: удельной поверхности твердой фазы, удельной поверхности конденсированной фазы, коэффициента влагопроводности, пористости, влажности, потенциала почвенной влаги.
Ключевые слова: удельная поверхность твердой фазы, удельная поверхность конденсированной фазы, коэффициент влагопроводности, пористость, влажность, потенциал почвенной влаги.
Известные способы оценки состояния почвы и степени воздействия на нее по измерениям твердости, по изменению объемной массы и основной гидрофизической характеристики (ОГХ) имеют недостатки связанные, в основном с длительностью эксперимента построения зависимости между потенциалом влаги и влажностью почвы (ОГХ) и с погрешностями при определении общепринятых характеристик почвы.
С целью сокращения времени эксперимента для определения большинства общепринятых характеристик почвы предлагается использовать аэродинамический метод, предусматривающий последующий перерасчет аэродинамических параметров в гидродинамические на основе теории подобия.
Рассматривая образец правильной геометрической формы, например цилиндрической, можно утверждать, что суммарный объем пор интерпретируется как объем одной цилиндрической поры эквивалентного радиуса, площадь сечения которой равна суммарной площади всех пор образца:
Введение
Метод
,2
>2
е!
>2 '
р _ _ еЧ __<Щ_
У, п %-Ах Я2'
Р_
(1)
где Ур - объем пар, м3; У - объем образца, м3; Ах - длина образца, м; ^ - радиус образца, м; Яед - эквивалентный радиус, м (см. рис. 1).
Я
Я
Рис. 1. Идеализированная модель образца почвогрунта
Боковая поверхность эквивалентной трубки тока равна поверхности конденсированной фазы образца. Отношение площади этой поверхности к объему конденсированной фазы V/, м3 есть ее удельная поверхность 0,с/, м2/м3:
п _ 2Яед п-М
"V
(2)
с/
где длина А/ трубки тока эквивалентного радиуса определяется экспериментально, исходя из того, что при протекании газа через образец за счет разности давлений между его торцами поток теряет часть своей кинетической энергии на преодоление трения о поверхность конденсированной фазы. Из теории подобия следует, что в нашем случае потери энергии на трение сравнимы с кинетической энергией самого потока, то есть поток описывается уравнением Пуа-зейля [1]:
п _п МР Я 4
У _о---7/4,
8ц М/ у
(3)
3 Мр 2 2
где У - объемный расход газа, м /с;--градиент давления, кг/с -м ; щ - диМ/
намическая вязкость газа, кг/с-м; Я/ - эффективный радиус, м (см. рис. 1). Выражая Яед из уравнения (1), запишем для конденсированной фазы:
2ЯпМ/ л/Р
Я/ _-
(4)
V (1 - Р)
Удельная поверхность твердой фазы, то есть сухого образца, определяется
выражением
Я _ 2ЯП о
V,
(1 - Ро).
(5)
С величиной удельной поверхности тесно связан потенциал почвенной влаги, определяемой как энергия, необходимая для переноса единицы массы
жидкости из образца в свободную жидкость и выражаемая обычно через эквивалентное давление. Для почвы или другого пористого материала можно записать:
¥ = + (6)
где ¥ - полный потенциал, Дж/кг, или Па (для воды); V' - адсорбционный потенциал, отражающий влияние твердой фазы (матрицы) на понижение свободной энергии воды в образце; V* - капиллярный потенциал, определяемый вкладом поверхностной энергии на границе раздела жидкость - газ.
Эквивалентное давление, отражающее адсорбционный потенциал, можно найти по формуле Б.В. Дерягина [1, с. 147]:
Ар = л/к3, (7)
где А - постоянная, меняющаяся в пределах (5 • 10-21 - 5 • 10-19) Дж. Определяя к по уравнению Б.Н. Мичурина [2]:
и = цу/а, (8)
где Жу - объемная влажность образца, м3/м3; О - удельная поверхность, м2/м3; к - толщина водной плёнки, м, можно записать:
А ' ла0 ^ Ар ла0
Ар =—3р, или V = — =-, (9)
Ж р „ рХ
где рк - плотность жидкости, кг/м3.
Величина V" определяется с помощью выражения
а гО,
V'' = _JLJL, (10)
где ох - удельная свободная поверхностная энергия на границе раздела жидкость - газ, Дж/м2.
Выделяя из (4) и (5) постоянной сомножитель и приравнивая правые части, можно записать
аг(1 -Р) ар(1 -Рр)
А!1У[Рг А!о4Р '
(11)
где индекс i относится к любому промежуточному значению величины при объемной влажности в пределах от 0 до Х = P0; с учетом (1) уравнение (11) может быть записано в виде
а0 =
= аг (1 - р) р24р = аг (1 - р)
(Р ^ 1 г
2.5
0
чРу
(12)
(1 -Р0)Р^,/Р (1 -р0)
где Рг - измеренное или заданное значение пористости в долях единицы.
С учетом того, что Р0 = Ж + Р есть расчетная формула для полного потенциала, выражающая его зависимость от объемной влажности, основную гидрофизическую характеристику можно записать в виде:
- Ро) (р0 - К,)2" + (13)
РР25 1 - (Р - Ки)
Первый член уравнения (13) отражает влияние капиллярной составляющей потенциала, второй - пленочной. Переход от капиллярно-пленочной формы влаги к капиллярной определяется уравнением А. Д. Воронина [3]:
^ 4-2.18 + . (14)
Поэтому при влажности, большой, чем определяемая по уравнению (14), второй член в уравнении (13) равен нулю,
Потенциал влаги обычно выражается через эквивалентное давление
р = р' + * ^ (1 - Ро) (Ро - К )25 + М
Ро25 1 - (Ро - К) жу У '
Интегрируя это выражение по объемной влажности, получим
\pciW - (р2 - а2)+ ЪЩ^ х
3 2 р а2 1 ' Р25
а г
х {(р^ + (р^ + (р - а)о.5 - агс1Ь(р - а)0-5 4
(16)
Нижний предел интегрирования (Ж = а) отвечает значению влажности в момент взятия образца.
Для сравнения образцов, взятых при различной влажности, значения а выбирается соответствующим влажности «устойчивого завядания». По своей физической сущности физическая величина интеграла выражает удельную объемную энергию Гиббса конденсированной фазы Е, Дж/м3.
Почвогрунты представляют собой дисперсную систему с сильно развитой поверхностью. Их подвижные отдельности удерживаются между собой благодаря молекулярным силам взаимодействия между твердыми частицами и жидкостью. Поэтому физико-механические свойства, такие, как твердость, пластичность и т. д., зависят от количества содержащейся в нем влаги для каждого конкретного почвогрунта. При разрыхлении или уплотнении почвогрунта мы, с точки зрения термодинамики, совершаем над ним работу, которая идет на изменение энергии связей между подвижными отдельностями почвогрунта. Измерив разность энергии связи, после механического воздействия и до воздействия, мы получаем значения работы, совершенной над системой. Эта работа характеризует степень воздействия на почвогрунт и позволяет оценить как энергетическое состояние почвогрунта, так и эффективность механического воздействия на него.
Если до воздействия величина интеграла составляла Ео, а после него - Е, то разность ДЕ = Е - Ео покажет изменение удельной свободной энергии Гиббса. А отношение ДЕ/Е можно использовать как показатель степени воздействия на почвогрунт.
6 А з 2 I 5 2—
Таким образом, мы получили возможность количественно оценить механические воздействия того или иного обрабатывающего орудия.
Величина Е0 может быть использована для объективной оценки качества почвы сельскохозяйственного назначения. Принимая значение Е0 для наиболее плодородной почвы (чернозема) за эталон и измеряя значение той же величины для данной почвы, можно оценить ее качество относительно эталона объективно и однозначно.
Первичные измеряемые величины определяются экспериментально при помощи разработанного нами устройства. На рис. 2 показана принципиальная схема устройства.
Устройство состоит из кассеты 1, содержащей образец почвогрунта 2 между крупнопористыми мембранами 3. Кассета 1 соединена с торцов с входной 4 и выходной 5 предкамерами, к которым присоединены патрубки 6 и 7 одинакового поперечного сечения, к которым, в свою очередь, подсоединены колена и-образного манометра 8. Патрубок 7 через уплотнительную пробку 9 введен в заполненный водой аспиратор 10 с нижним выходным отверстием 11, снабженным герметизирующей прокладкой 12 и сливным краном 13. Для измерения объема вытекающей жидкости под краном 13 установлен мерный цилиндр 14.
Способ и устройство реализованы следующим образом. Собранное устройство приводится в действие открытием крана 13. При этом расход жидкости остается постоянным, пока уровень жидкости в аспираторе 10 выше нижнего среза патрубка 7. Скорость вытекания жидкости регулируется разностью уровней между нижними срезами патрубка 7 и крана 13, а также степенью открытости крана 13.
После установления постоянной разности уровней в коленах манометра 8 под кран 13 устанавливают мерный цилиндр 14 и измеряют время вытекания определенного объема жидкости, равного объему воздуха, протекшего через образец. Из эксперимента определяют разность уровней в коленах манометра,
Методика измерений и оборудование
объем газа, протекшего через образец, время протекания. Предварительно измеряют пористость влажного образца, например, способом по а. с. СССР № 1260763 А1, М.кл. G 01 N 15/08, и влажность образца методом сушки.
Полученные данные подставляют в формулу для определения удельной поверхности твердой фазы
Q = 1 П. (pw + W)2 • R -УPw + wvApAt 0 4 пэ AxAV '
Для определения удельной поверхности конденсированной фазы используется следующая формула:
1 п R3 • P2 'л/Pw •ApAt
bV --т—Г77-, (18)
^ 4 п AxAV
где п" - удельная поверхность твердой фазы, м2/м3; Pw - пористость влажного образца в долях единицы; Wv - объемная влажность образца в долях единицы; Ap - перепад давления между торцами образца, Па; R - радиус образца, м; At - время протекания измеренного объема, с; AV - объем газа, прошедшего через образец, м3; Ax - длина образца, м; Qcf -удельная поверхность конденсированной фазы, м2/м3.
Результаты
Примеры реализации способа приведены в табл. 1-4.
Табл. 1
Гидрофизические параметры для исследуемых типов почв
а)
Материалы Wv P 1 w ApAt, Па с Ar, мм AV , м3
00 045
песок аллювиальный 01 035
02 025 90.23 100 О О
03 015 о
045 00
00 052
темно-серая 01 042
02 032 1183.35 100 о
лесная песчаная 03 022 о о
04 012
052 00
00 055
01 045
чернозем 02 035 1316.52 100 0.0001
выщелоченный 03 025
04 015
055 00
б)
Материалы П-10-5, кг/м3 и л /\—3 V -10 3, кг/м3 Цу -106, м/м3 ^, Дж/кг
песок аллювиальный 1.72 1.05 3.32 -
1.77 111.23
0.76 41.11
0.21 10.12
0 -
темно-серая лесная песчаная 1.72 1.05 62.5 -
36.64 2235.0
18.56 972.0
7.27 331.4
1.59 65.57
0 -
чернозем выщелоченный 1.72 1.05 80 -
48.44 3149.2
25.84 1338.1
11.14 497.56
3.10 124.1
0 -
Табл. 2
Вспашка под однолетние травы после однолетних трав
Орудие Средняя удельная поверхность, м2/см3 Объемная масса, кг/м3 Пористость Удельная энергия, кДж/м3 АЕ/Е0
до после до после до после
Т-25Н+ экспер.плуг 117 1.24 1.14 0.581 0.62 287 315 0.09
ДТ-75Н+ ПЛН-4-35 1.21 0.59 293 0.02
Т-25Н+э кспер.плуг 127 1.30 1.16 0.562 0.619 271 316 0.17
ДТ-75Н+ ПЛН-4-35 1.19 0.607 307 0.13
Т-25Н+ экспер.плуг 142 1.34 1.15 0.573 0.641 338 397 0.17
ДТ-75Н+ ПЛН-4-35 1.17 0.621 381 0.125
На первом поле: объемная влажность на поле 0.23, плотность тв фазы 2.41 кг/м3. На втором поле: объемная влажность на поле 0.22, плотность тв фазы 2.47 кг/м3. На третьем поле: объемная влажность на поле 0.22, плотность тв фазы 2.52 кг/м3.
Из данных таблиц следует, что на всех полях агрегат Т-25Н+экспер.плуг более эффективен. Преимущество этого агрегата более выражено для почв с меньшей удельной поверхностью.
Табл. 3
Боронование под пшеницу и ячмень
Орудие Средняя удельная поверхность, м2/см3 Объемная масса, кг/м3 Пористость Удельная энергия, кДж/м3 АЕ/Е0
до после до после до после
Т25+эксп. аф- ронтальн. игольч.борона 105 1.46 1.33 0.55 0.60 219 254 0.16
Т25+эксп. фронтальн. игольч.борона 1.23 0.63 276 0.26
ДТ-75Н+ БИГ-3 1.38 0.58 243 0.11
ДТ-75Н+ БЗТС-1,0 1.39 0.58 241 0.10
Т25+эксп. афронтальн. игольч.борона 131 1.49 1.28 0.58 0.63 310 350 0.13
Т25+эксп. фронтальн. игольч.борона 1.26 0.70 392 0.26
ДТ-75Н+ БИГ-3 1.35 0.58 311 0.003
ДТ-75Н+ БЗТС-1,0 1.36 0.58 315 0.016
На первом поле: объемная влажность на поле 0.29, плотность тв. фазы 2.63 кг/м3. На втором поле: объемная влажность на поле 0.31, плотность тв. фазы 2.51 кг/м3.
Из табл. 3 следует, что из экспериментальных борон наиболее эффективна фронтальная игольчатая борона на всех полях. Из серийных агрегатов более эффективен агрегат ДТ-75Н+БЗТС-1,0.
Табл. 4
Прикатывание после сева гороха (МТЗ-82+КЗК-10), Культивация с выравниванием и прикатыванием (ДТ-75Н+РВК-3,6)
Орудие Средняя удельная поверхность, м2/см3 Объемная масса, кг/м3 Пористость Удельная энергия, кДж/м3 АЕ/Е0
до после до после до после
МТЗ-82+КЗК-10 106 1.226 1.397 0.55 0.48 215 165 -0.23
ДТ-75Н+РВК-3,6 112 1.42 1.40 0.49 0.49 144 150 0.04
На первом поле: объемная влажность на поле 0.11, плотность тв. фазы 2.46 кг/м3. На втором поле: объемная влажность на поле 0.16, плотность тв. фазы 2.45 кг/м3.
Из данных таблицы следует, что агрегат МТЗ-82+КЗК-10 уплотняет почву, а агрегат ДТ-75Н+РВК-3,6 прикатывает почву без уплотнения.
Summary
V.V. Sirotkin, V.M. Sirotkin. Energetic Estimation of Soil Condition and Degree of Mechanical Influence on It by Soil-Cultivating Machines.
Method of estimating the energy condition of soil and degree of mechanical influence on it by soil-cultivating machines is developed on the basis of a uniform integrated parameter. The latter is determined by experimental and calculate means on the basis of measuring standard values, such as firm phase specific surface, condensed phase specific surface, moisture conductivity coefficient, porosity, humidity, soil moisture potential.
Key words: firm phase specific surface, condensed phase specific surface, moisture conductivity coefficient, porosity, humidity, soil moisture potential.
Литература
1. Сироткин В.В., Сироткин В.М. Прикладная гидрофизика почв. - Чебоксары: Изд-во Чуваш. ун-та, 2001. - 250 с.
2. Мичурин Б.Н. Энергетика почвенной влаги. - Л.: Гидрометоиздат., 1975. - 137 с.
3. Воронин А.Д. Энергетическая концепция физического состояния почв // Почвоведение. - 1990. - № 5. - С. 7-19.
Поступила в редакцию 22.04.08
Сироткин Вячеслав Владимирович - доктор географических наук, профессор кафедры ландшафтной экологии Казанского государственного университета.
E-mail: [email protected]
Сироткин Владимир Михайлович - доктор технических наук, профессор, заведующий кафедрой физики Чувашской государственной сельскохозяйственной академии, г. Чебоксары.
