(11)

Для определения ущерба урожая от недополива использовалась зависимость

(12)

где У - проектная урожайность в конкретные годы или в среднем за многолетний период; У_max - максимальная урожайность, получаемая в передовых хозяйствах с аналогичными почвами или на сортоиспытательных станциях; К_у - коэффициент снижения урожайности за счет неблагоприятного водного режима мелиорируемых земель; К_о, К_wi K_i, K_i - коэффициенты, учитывающие снижение урожайности в результате соответственно: запаздывания со сроками сева, отклонений от оптимальных значений влажности почвы W, уровня грунтовых вод и периодических колебаний уровня грунтовых вод с амплитудой ; а_i - вес фазы, отражающий чувствительность растений в каждой фазе развития с порядковым номером i к отклонениям водно-воздушного режима почвы от оптимальных условий.

Урожайность люцерны при полной водообеспеченности принималась по данным хозяйства У = 10 т/га за 4 укоса, К_wi и а_i принимаем соответственно по данным В.В. Шабанова. K_i и K_i принимаем равными 1 при глубоких уровнях грунтовых вод.

Величина оросительной нормы М при расчете по методу min-max У составила 4200 м³/га при средней влажности расчетного слоя почвогрунта 0,63 НВ в течение вегетационного периода люцерны.

При расчетах поливного режима для сухого 1981 года с использованием среднемноголетних данных норма орошения М была 4400 м³/га, влажность расчетного слоя - 0,65 НВ.

Расчет поливного режима люцерны с использованием прогноза, выдававшегося гидрометслужбой в 1981 году показал. Что оросительная норма должна была составить 4800 м³/га, при средней влажности - 0,68 НВ.

Сравнение поливного режима, рассчитанного различными методами, позволяет сделать вывод о возможности применения при оптимальном планировании метода «min-max У», и метода с использованием среднемноголетних погодных данных.

В таблице 4 приведен пример выбора величины управляющего воздействия при 3-м подходе. Все расчеты выполнялись по формулам (9), (10), (12). Так, при величине поливной нормы m = 600 м³/га, инфильтрация воды за пределы расчетного слоя будет g = 170 м³/га, стоимость вымыва гумуса составит в денежном выражении 0,25 х 170 = 42,5 руб/га. Стоимость вымыва ионов Са принята 8 руб/га. Затраты на воду принимаются в зависимости от величины поливной нормы и замыкающих затрат на воду в источнике , поличаем 432 руб/га.

Расходы на электроэнергию рассчитываются по формуле

руб/га (13)

где m = 600 м³/га - величина поливной нормы; Н = 100 м - средневзвешенный годовой напор; S - тариф на электроэнергию; 0,004 - удельные затраты на электроэнергию для подъема 1 м³ воды на высоту 1 м (кВт.ч).

Таблица 4

Расчет ущерба при орошении различными нормами, III декада июня, min-max У, руб/га

Ущерб от недобора урожая можно получить при использовании графиков и зависимостей Ю.Н. Никольского и В.В. Шабанова. В данном случае при поливной норме m = 600 м³/га во влажный год он будет составлять 72,1 руб/га.

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

M = 0 м³/га; У_max = 490 руб/га

m = 400 м³/га; У_max = 64,2 руб/га

m = 600 м³/га; У_max = 71,2 руб/га

а затем выбирается минимальный ущерб У_min = 64,2 руб/га при m = 400 м³/га. Управляющее воздействие при m = 400 м³/га считается оптимальным и применяется в дальнейших расчетах для текущего года.

Выбор того или иного метода принятия решения определяется целями хозяйства-водопользователя.

Оперативное определение сроков и норм поливов, основанное на предварительных расчетах с использованием ЭВМ сезонных планов для различных сочетаний погодных условий

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

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

Попытки решить проблему оперативного определения сроков и норм поливов делались неоднократно. Однако при применении разработанных методов требуется использование ЭВМ непосредственно в процессе оперативного планирования. В перспективе это не будет вызывать затруднений, но в настоящее время в практике эти методы не применяются.

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

Сущность способа оперативного расчета сроков и норм полива

В основу расчетов, как и во всех методиках оперативного планирования поливного режима, положено уравнение водного баланса орошаемого поля.

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

(14)

где i - номер декады вегетационного периода рассматриваемой сельскохозяйственной культуры; п - число декад в вегетационном периоде; W_нi ,W_кi - доступные растениям запасы воды в расчетном слое почвы соответственно в начале и конце декады; E_i - водопотребление сельскохозяйственной культуры; Р_i - усвояемые атмосферные осадки; q_i - поступление воды в расчетный слой почвы за счет капиллярного поднятия грунтовых вод; m_i -поливная норма.

Запасы влаги в расчетном слое почвы, соответствующие наименьшей влагоемкости (W_НВ), вычисляются по формуле

мм, (15)

где Н_i - глубина расчетного слоя почвы в i-ю декаду, м; - средняя по глубине плотность почвы, т/м³; щ_нв - средняя по глубине влажность в долях от массы сухой почвы, соответствующая наименьшей влагоемкости.

Запасы влаги в расчетном слое почвы, соответствующие предполивной влажности (W_п), составляют

мм, (16)

где щ_п - предполивная влажность в долях от веса сухой почвы; r - предполивная влажность в долях от влажности, соответствующей наименьшей влагоемкости.

Значение r может быть разным для различных сельскохозяйственных культур и изменяться в течение вегетационного периода. Однако во многих случаях принимают r = 0,7.

В течение вегетации происходит постепенное углубление расчетного слоя почвы, что вызвано ростом корневой системы, поэтому W_нвi и W_пi должны постепенно возрастать.

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

Доступные растениям влагозапасы в расчетном слое почвы в i-ю декаду могут изменяться от 0 до максимальной величины

мм. (17)

Для расчета водопотребления используется формула Г.К. Льглва:

мм. (18)

где t_i - сумма среднесуточных температур воздуха в i-ю декаду;

K_i - биоклиматический коэффициент.

Биоклиматический коэффициент принимается на основании рекомендаций научных учреждений для каждой культуры и региона. Он изменяется в течение вегетационного периода.

Формула (18) принята для расчета потому, что в ее основе лежит температура воздуха. Это в настоящее время единственная метеорологическая характеристика, на которую дается декадный прогноз.

Каждая таблица составляется для всего вегетационного периода с разбивкой по декадам. Число таблиц для каждого поля зависит от диапазона изменения температур воздуха и атмосферных осадков, наблюдавшихся в районе расположения за многолетний период. Учитывая точность прогноза, расчеты проводятся для температур с шагом 5⁰С. Например, если из имеющихся исходных данных на 25...30 лет видно, что температуры в период вегетации могут изменяться на 3-8⁰С, то таблицы нужно составлять для температур t = 10, 15, 20, 25, 30⁰С. Осадки в настоящее время гидрометеослужбой не прогнозируются, поэтому можно принимать их среднемноголетние значения для каждой декады соответствующего месяца при заданном температурном диапазоне. Например, из 25-летнего ряда наблюдений в 7 годах температура во II декаде июня была 202,5⁰С. Все годы, в которые для данной декады среднедекадная температура находилась в таком диапазоне, выписываются в отдельную строку и им сопоставляется сумма декадных осадков, наблюдавшаяся в этой декаде за каждый год выборки. Затем вычисляется средняя декадная сумма осадков за эти годы. Если наблюдающийся диапазон изменения осадков (О_imax - O_imin) велик (более 15...20 мм), то его можно разбить на ряд интервалов, например, 0_iк 5, мм, к - номер интервала в i-й декаде, к =1,…, п_i. Тогда число интервалов

(19)

п_i - округляется до целого. Интервал осадков, равный 10 мм, вполне приемлем с точки зрения точности выполнения расчетов, ошибка в величине составляющих водного баланса, равная 5 мм, допустима, так как она приведет к ошибке в дате полива не более, чем на 1...2 суток. Если выбирается несколько расчетных значений осадков при данной температуре, то для каждого сочетания t_si, O_ki составляется отдельная таблица.

Расчет проводится для каждой декады вегетационного периода в следующем порядке. Выписываются сумма осадков и значения биоклиматических коэффициентов. Вычисляется величина водопотребления по формуле (18) и поступление воды в расчетную зону за счет капиллярного подпитывания. Затем определяется суточный дефицит водопотребления

мм/сут. (20)

где n_i - число дней в i-й декаде, n_i = {10.11}. Теперь, задаваясь различными значениями начальных влагозапасов (W_н), можно определить нужен ли полив в текущую декаду, какого числа его провести и какой поливной нормой. Учитывая, что все составляющие водного баланса вычисляются с не очень большой точностью, начальные влагозапасы следует принимать в расчетах с разницей между ближайшими значениями 10...20 мм. Это приведет при выборе поливной нормы к ее изменению на 5...10 мм.

Оценки показывают, что для этого достаточно изменять W_н с шагом (0,05...0,10) W_н . Сначала принимается W_н = 0,1 W_нв и проводятся дальнейшие расчеты.

Определяются доступные влагозапасы в расчетном слое почвы на начало i-й декады.

мм, (21)

и свободная влагоемкость в этом слое на тот же момент времени

мм. (22)

Следующим шагом является определение даты полива при разных поливных нормах. Номер дня i- й декады (N_mi), в который нужно проводить полив поливной нормой m_l, можно вычислить по соотношению

сут. (23)

Здесь D_ci - cуточный дефицит водопотребления, рассчитываемый по выражению (20).

В практике орошаемого земледелия в засушливой зоне России при использовании дождевальной техники обычно поливные нормы назначают кратными 10 мм и принимают равными m_l = {0, 40, 50, 60} мм. Для этих значений вычисляются возможные даты поливов в каждой декаде.

После этого определяются влагозапасы в конце i-й декады для каждой поливной нормы

мм (24)

Если в i-ю декаду не выполняется условие m_li - W_cвi 0, то в эту декаду полив не предусматривается и принимается m_l = 0.

Такие расчеты осуществляются для всех сельскохозяйственных культур при каждой t_si в диапазоне {t_min, t_max} с интервалом 5⁰С, при каждом значении 0_кi, каждых W_нi и m_l.

В результате всего комплекса расчетов можно получить альбом таблиц для орошаемых полей севооборота.

Основные выводы и предложения

Впервые для условий Саратовского Заволжья на основе системного подхода и теоретических обобщений выполнен анализ материалов, характеризующих агроклиматические и почвенные условия, где:

с использованием закона географической зональности почв Докучаева-Григорьева-Будыко был рассчитан и построен график внутризонального изменения агрохимических свойств почв в зависимости от гидротермического показателя;

определена почвенно-экологическая и биологическая продуктивность почв в зависимости от гидротермического показателя и количественных значений необходимых расчетных параметров, их оценки направленности и интенсивности биологического и геологического круговоротов;

на основе сравнительных расчетов биологической продуктивности почв и урожайности сельскохозяйственных культур в естественных условиях и после мелиорации, определена «граница оросительных мелиораций» в пределах агроклиматических зон с индексом сухости J = 1,0...1,8.

2. Практикой современного использования орошаемых земель и исследованиями установлено, что оросительные системы Саратовского Заволжья имеют большие резервы по дальнейшему совершенствованию методов оперативного управлениями водными ресурсами при изменяющихся стохастических условиях погоды, что является одной из важных задач повышения эффективности орошения.

На основании исследований на Энгельсской оросительной системе разработаны методики планирования и контроля водопользования на оросительных системах засушливой зоны. Основные положения методик внедрены в практику работы Энгельсской оросительной системы.

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

4. На основании анализа системных планов водопользования и годовых отчетов управления эксплуатации Энгельсской оросительной системы 1966-2000 гг. выявлено, что объемы заборов воды на орошение в годы различной засушливости значительно отличаются от плановых и расчетных значений. При стохастических условиях погоды недостаточно достоверных нормативных данных по эксплуатационным режимам орошения культур для планового водопользования.

5. Оптимальное суммарное водопоступление на территорию с развитым орошением составляет для условий Саратовского Заволжья (Энгельсская система) 7200 м³/га. Годовая его величина может достигать 11000м³/га.

6. Расчетами установлено значительное влияние погодных условий на режимы водопользования. За период 1951 по 1995 годы оросительные нормы по недостаткам водопотребления изменялись от 958 до 216 мм. Отмечались критические периоды в отдельные месяцы, когда не обеспечивались поливы, или были избытки воды в системе, которыми можно поливать дополнительные площади. Необходимо определять вероятностные режимы водопользования на основе анализа климатических данных за последние 25 лет и намечать мероприятия по устранению пиковых потребностей воды на орошение и по использованию оросительной воды на поливы дополнительных площадей.

7. Рекомендуется плановое водопользование на оросительных системах засушливой зоны проводить на основе дифференцированных режимов орошения. Определять для характерных по засушливости лет расчетные ординаты графиков гидромодуля по декадам, расчеты водоподачи проводить для года 75% обеспеченности. При проведении водопользования в начале декады корректировать расчеты заборов и подачи воды на орошение по прогнозу недостатков водопотребления.

8. Разработка методов оперативного и достаточно точного управления водными ресурсами оросительной системы засушливой зоны является важной задачей повышения эффективности орошения в этой зоне. На протяжении последних 30 лет ХХ века проводились исследования на Энгельсской оросительной системе, на основании которых разработаны методы планирования водопользования, основные положения методик внедрены в производство.

9. Экономическая эффективность рекомендуемых мероприятий от внедрения обоснованного в данной диссертации мелиоративного режима орошаемых земель и совершенствования водопользования в Саратовском Заволжье определялась по методическим рекомендациям Министерства сельского хозяйства РФ. . Анализ расчетов показывает, что простой срок окупаемости составил семь лет, а чистый дисконтированный доход - 17611 руб/га. Сравнивая результаты можно сделать вывод, что внедрение научных разработок улучшает экономические показатели. Так простой срок окупаемости сокращается на один год, а чистый дисконтированный доход увеличивается в 1,12 раза.

СПИСОК ОПУБЛИКОВАННЫХ РАБОТ

В изданиях, рекомендуемых ВАК РФ:

1. Обобщенный метод определения потребности оросительной системы в водных ресурсах //Мелиорация и водное хозяйство: научно-практический журнал. - М., 1998, двухмесячник, № 4. с. 4/2. Бегляров Д.С., Рыбкин В.Н. ISBN 0235-2524

2. Управление оросительными системами: состояние, проблемы и пути их решения. Монография. - М.: МГУП, 2006. - 217 с., 500 экз. Рыбкин В.Н. ISBN 5-89231-187-2

3. Анализ стохастического характера функционирования оросительных систем (на примере Энгельсской ОС) //Мелиорация и водное хозяйство: научно-практический журнал. - М., 2006, двухмесячник, № 4. с. 6/3. Буркова Ю.В., Рыбкин В.Н. ISSN 0235-2524.

4. Методика прогнозирования водопотребности на оросительных системах. //Мелиорация и водное хозяйство: научно-практический журнал. - М., 2007, двухмесячник, № 5. 6/3. Перминов А.В., Рыбкин. ISSN 0235-2524.

5. Анализ эксплуатационных данных по водопользованию //Вестник Российской академии сельскохозяйственных наук.- М., 2008, двухмесячник, № 6, с. 33-34. ISSN 0869-3730.

6 Методика расчета режима орошения с учетом вероятностного характера погодных условий //Мелиорация и водное хозяйство. - М., 2008, двухмесячник, № 6, с. 39-41. Буркова Ю.Г., Рыбкин В.Н. ISSN 0235-2524.

7. Методика корректировки планового водопользования /Доклады Российской академии сельскохозяйственных наук. - М., 2009, № 1, с. 40-41. ISSN 0869-6128.

8. Орошение и водопользование в Саратовском Заволжье //Мелиорация и водное хозяйство. - М., 2009, двухмесячник, № 3, с. 22-24. Гостищев Д.П., Рыбкин В.Н. ISSN 0235-2524.

В других изданиях:

9. Учет воды на оросительных системах / Вопросы мелиорации. 1988, № 5-6. 5/2,5. Бегляров Д.С., Рыбкин В.Н.

10. Методы расчета поливного режима сельскохозяйственных культур с учетом вероятностного характера погодных условий //Сельскохозяйственные гидротехнические мелиорации: Сб. науч. трудов. - М.: МГМИ, 1989, 92 с.,8/4. Рыбкин В.Н., Галямин Е.П. ISBN 0-000-00000.

11. Совершенствование техники полива виноградника при реконструкции внутрихозяйственной оросительной сети в совхозе «Коммунизм» Ура-Тюбинского района Таджикской ССР [Текст]: отчет о НИР (заключит.). М.: МГМИ, 1990. Деп. № ГР 0186.0117141. - Инв. № 029.10041211. 63/20. Рыбкин В.Н., Сурин В.А.

12. Водосбрегающая технология поверхностного полива сельскохозяйственных культур при реконструкции внутрихозяйственной оросительной сети на больших уклонах Ленинабадской области Таджикской ССР [Текст]: отчет о НИР (заключит.). - М.: МГМИ, 1991. Деп. № ГР 0186.0117141. - Инв. № 029.10041212. 81/20. Рыбкин В.Н., Сурин В.А.

13. Обоснование реконструкции водохозяйственной оросительной сети // Эксплуатация мелиоративных и водохозяйственных сооружений: Информационный сборник. - М.: ЦБНТИ, 1991, № 53. - 3 с.

14. Внедрение хозрасчета на оросительных системах засушливой зоны при платном водопользовании // Сельскохозяйственные гидротехнические мелиорации: Докл. научно-техн. конф. МГМИ. - М.: МГМИ, 1991. - 114/2. ISBN 0-000-00000

15. Проблемы развития мелиораций в засушливой зоне //Вопросы мелиорации: научно-практический журнал. - М., 1998. - двухмесячник. - № 4. 5/2,5. Бегляров Д.С., Рыбкин В.Н.

16. Обобщенный метод определения потребности оросительной системы в водных ресурсах //Мелиорация и водное хозяйство: научно-практический журнал. - М., 1998. - двухмесячник. - № 4. 4/2. Бегляров Д.С., Рыбкин В.Н. ISBN 0235-2524.

17. Учет воды на оросительных системах // Вопросы мелиорации: научно-практический журнал. - М., 1998. - двухмесячник. - № 5-6. 5/2,5. Бегляров Д.С., Рыбкин В.Н.

18. Методы проектирования открытых каналов оросительных систем, обеспечивающие водо- и энергосбережение. Пособие к СНиП 2.06.03-85. М.: Минсельхоз РФ Депмелиоводхоз, ГП СНЦ «Госэкомелиовод». 2000. 10/6. Рыбкин В.Н., Казаков Н.Т.

19. Влияние погодных условий на режим работы насосных станций оросительных систем // Природоохранное обустройство территорий. Материалы Международной научно-технич. конф. - М.: МГУП. 2002. 2/1. Буркова Ю.Г., Рыбкин В.Н. ISBN 978-5-89231-243-1.

20. Эксплуатационная надежность закрытого горизонтального дренажа //Проблемы экологии и безопасности жизнедеятельности в ХХ1 веке: Материалы Международной научно-практич. конф.: Норма, 2002, вып. 3. 6/2. Ильина Н.С., Ильин С.П., Рыбкин В.Н., Сильченков И.С., Гурин В.А. ISBN 5-85302-352-7.

21. Погодные условия и режим работы насосных станций // Вопросы мелиорации: научно-практический журнал. - М., 2002, двухмесячник. - № 5-6. 5/3. Буркова Ю.Г., Рыбкин В.Н.

22. Природоохранные мероприятия при освоении и рекультивации площадей //Проблемы экологии и безопасности жизнедеятельности: материалы Международной научно-практич. конф. - М.: Норма, 2003, вып. 4. - 213 с., 3/2. Ильин С.П., Рыбкин В.Н. ISBN 5-85302-352-7.

23. Эксплуатационные мероприятия по защите земель от эрозии почвы на мелиоративных системах // Проблемы экологии и безопасности жизнедеятельности: материалы Международной научно-практич. конф. - М.: Норма, 2003, вып. 4. - 213 с. 6/2. Ильина Н.С., Ильин С.П., Рыбкин В.Н., Сильченков И.С. ISBN 5-85302-352-7.

24. Учет влияния климатических факторов на требуемые подачи насосных станций при имитационном моделировании их работы //Экологическая устойчивость природных систем и роль природообустройства в ее обеспечении: Сб. материалов Всерос. научно-технич. конф. - М.: МГУП, 2003. - 198 с. 4/3. Буркова Ю.Г., Рыбкин В.Н. ISBN 5-89231-116-3.

25. Охрана атмосферного воздуха //Проблемы экологической безопасности и природопользования: Материалы Междунар. практ. конф. - М.: Норма, 2004, вып. 5. - 320 с. 5/2. Ильина Н.С., Ильин С.П., Рыбкин В.Н., Сильченков И.С. ISBN 5-85302-352-7.

26. Охрана растительности //Проблемы экологической безопасности и природопользования: Материалы Междунар. практ. конф. - М.: Норма, 2004, вып. 5. - 320 с. 6/2. Ильина Н.С., Ильин С.П., Рыбкин В.Н., Сильченков И.С. ISBN 5-85302-352-7.

27. Учет режимов эксплуатации крупных насосных станции при имитационном моделировании их работы // Проблемы научного обеспечения развития эколого-мелиоративного потенциала России: Сб. матер. научно-техн. конф. - М.: МГУП, 2004. 6/2. Буркова Ю.Г., Рыбкин В.Н. ISBN 5-89231-140-6.

28. Полевые исследования эксплуатационных режимов поливов при дождевании. Опытные участки. Методика исследований. //Проблемы экологической безопасности и природопользования. Матер. Международ. научно-практической конф. - М.: Норма, МАЭБП, 2006. Вып 7. 5/3. ISBN 5-85302-352-7. Рыбкин В.Н., Кузнецова Е.В.

29. Эксплуатация и мониторинг мелиоративных систем. Учебник для вузов. - М.: ООО «Инлайт», 2006. 390/130. Ольгаренко В.И., Ольгаренко Г.В. ISВN 5-91040-002-9.

30. Экологизация водопользования на оросительных системах (на примере Энгельсской оросительной системы) // Проблемы экологической безопасности и природопользования: Материалы Междунар. научно-практ. конф. - М.: Норма, 2006, вып. 7. 6/3. Кузнецова Е.В., Рыбкин В.Н. ISBN 5-85302-352-7.

31. Оптимизация эксплуатационного режима орошения сельскохозяйственных культур //Роль мелиорации и водного хозяйства в реализации национальных проектов. Материалы Мждународной научно-практич. конф. Часть 1. - М.: МГУП, 2008. - 6/3. Рыбкин В.Н., Буркова Ю.Г. ISBN 978-5-89231-243-1.

Размещено на Allbest.ru