Современные методологические подходы к расчету параметров гидротехнических отстойников Текст научной статьи по специальности «Математика»

УДК 502/504:626.823.4

Н.Т. Кавешников, доктор техн. наук, профессор

Федеральное государственное образовательное учреждение высшего профессионального образования «Московский государственный университет природообустройства»

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

Показано, что традиционные подходы к разделению расчетов отстойников на энергетические и ирригационные несостоятельны. Мировой опыт свидетельствует о другом: на деривационной ГЭС в Непале, например, наносы очень мелкие (в основном 0,01...0,1 мм, из них мельче 0,05 мм составляет 80 %), но концентрация достигает 6,0 кг/м3. Если крупность наносов на подходе к гидроузлу соответствует условиям ирригационных, а отстойник энергетического назначения, то расчет следует выполнять по методике, служащей для ирригационного. Автор предлагает простую и многократно проверенную методику расчета осаждения мелких наносов в отстойнике, приводит конкретный пример.

It is shown that the traditional approaches to the division of calculations of settlers into power and irrigation ones are not consistent as the world experience shows. In particular there is given an example on the derivational power plant in Nepal where drifts are very small (mostly 0,01...0,1 mm, 80 % of them being smaller than 0,05 mm) but their concentration amounts to 6,0 kg/m3. The size of drifts in the approach zone to the hydro-electric scheme corresponds to the conditions which are characteristic of irrigation settlers though the settler is a power one. Hence, despite the purpose of the considered settler, it is necessary to carry out its calculation using the method worked out for the irrigation settler. In this article the author offers a simple and repeatedly checked by him method of calculation of settlement of small drifts in the settler illustrated by a particular example.

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

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

Для гидроэлектростанций, особенно деривационного типа, это задача, связанная с повреждением или разрушением направ-

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

более двух лет — при интенсивном истирании;

3-4 года — при умеренном истирании; менее 4 лет — при слабом истирании агрегатов наносами [1].

В соответствии с техническими условиями и нормами проектирования гидротехнических отстойников для гидроэлектростанций (ТУиН 24-110-48), действующими с 1948 г., по условиям допустимого абразивного воздействия на турбины через отстойник можно пропускать наносы крупностью до 0,25 мм с концентрацией в потоке не более 0,2 кг/м3.

Для оросительных систем, расположенных, например, на реках в Приморском крае, на Урале, Северном Кавказе, в Краснодарском и Ставропольском краях, основной проблемой является заиление наноса-

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

Инженерным решением такой задачи является устройство отстойников для осаждения наносов — отстойников ГЭС, ирригационных отстойников. Отечественный и зарубежный опыт защиты ирригационных систем от заиления наносами свидетельствует о том, что, как правило, в отстойниках необходимо осаждать наносы со средним диаметром фракции крупнее 0,01...0,0 5 мм (крайне редко менее 0,01 мм).

По допускаемой крупности наносов в агрегат ГЭС (0,25 мм) и в оросительную систему (0,05 мм) разработаны разные методологические подходы и методы определения рациональных параметров и условий эксплуатации отстойников. Это сложнейшие по трудоемкости и инженерному восприятию методы расчета, в том числе ранее известный метод А.Г. Хачатряна. В последнее десятилетие автором обнародован более простой и надежный метод расчета отстойников, строящихся на водных объектах, несущих большое количество мелких наносов, о чем будет сказано ниже [2].

Специалисты ОАО НИИЭС (В.В. Беликов, А.М. Прудовский и В.П. Середавкин) в своей статье отмечают, что мировой опыт не подтверждает рекомендации, заложенные в отечественном документе ТУиН 24-110-48, который более 55 лет не уточнялся [1]. Они констатируют: «Так, на деривационной ГЭС Marsyngdi в Непале, мощностью 69 МВт (напор — 95 м, расход — 100 м3/с), несмотря на наличие отстойника длиной 400 м и шириной 75 м, требуется проведение ежегодных ремонтно-восстановительных работ. Отстойник был запроектирован на осаждение частиц наносов крупнее 0,1 мм. По данным натурных измерений через турбинный тракт проходят частицы наносов с концентрацией в потоке до 6 кг/м3. При этом, обратим внимание читателя, частиц мельче 0,05 мм содержится более 80 % от общего количества наносов, и только размеры отдельных зерен песка превосходят 0,125 мм».

Отсюда можно сделать вывод, что отечественный метод расчета отстойников ГЭС, изложенный в ТУиН 24-110-48, к наносам

деривационной ГЭС Магвупд& в Непале, т.е. к наносам мелким, применять нельзя, ибо этот метод базируется на том условии, что наносы будут осаждаться сравнительно крупного размера (более 0,25 мм). В то же время происходит интенсивное абразивное воздействие на агрегаты ГЭС (инженеры и ученые Непала занимаются этой проблемой).

Итак, для отстойников ГЭС и ирригационных отстойников критическая максимальная фракция, на которую их следует рассчитывать, зависит от многих факторов [3]. При наносах, образующихся в результате прохождения интенсивных дождей и ливней, расчетной является фракция размером 0,01...0,05 мм. В этом случае обосновывающие рекомендации, служащие для проектирования и эксплуатации ирригационного либо гидроэнергетического отстойника, практически не отличаются.

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

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

Термины «транспортирующая способность» и «критическая мутность потока» зачастую трактуют одинаково, однако следует четко разделять эти понятия: транспортирующая способность — это потенциальные возможности потока нести предельное количество взвеси, если бы он был насыщен всей гаммой возможных фракций наносов, что на практике может быть крайне редко; критическая мутность — это предельное насыщение потока наносами данного фракционного состава, как пра-

вило, не всегда имеющего полную гамму разнородных частиц. Далее автор рассматривает вопросы именно критического состояния насыщения потока наносами в открытых руслах.

Существует свыше двух десятков различных эмпирических и полуэмпирических формул для определения величины критической мутности потока, авторами которых являются: В.В. Пославский, Г.С. Чекулаев, А.Г. Хачатрян, Г.О. Хорст, Е.А. Замарин, В.А. Шаумян, А.Н. Гостунс-кий, М.А. Великанов, П.В. Михеев, Г.В. Лопатин, С.Х. Абальянц, И.И. Леви, И.С. Яковлев, А.В. Ефремов, С.А. Гиршкан, И.И. Горошков, А.В. Караушев, А.А. Черкасов, К.И. Россинский, Ф.С. Салахов и др. Из числа зарубежных авторов, получивших также расчетные формулы для определения критической мутности, можно назвать следующие имена: Р. Кеннеди, Ша-юй цзинь, Фан Цзахуа, Чжан Жуэй-цзинь, Доу-Гуо-жень, Л. Богарди и др. Следует также отметить теоретические и экспериментальные работы по данной проблеме Н.Е. Жуковского, В.Д. Журина, В.Н. Гончарова, С.Т. Алтунина, Г.В. Железнякова, Ц.Е. Мирцхулавы, В.С. Алтунина, А.Б. Век-слера, О.Г. Натишвили, Н.А. Ржаницина, Е.К. Рабковой, Х.Ш. Шапиро, Т.А. Алиева, Г. А. Тер-Абрамянца, Н.А. Михайловой, В.И. Астраханцева, Г.И. Баренблатта, Ю.М. Кузьминова, Ю.А. Ибад-Заде и др.

Большинство из установленных зависимостей неоднократно анализировались в работах многих ученых, что позволяло им получать свои эмпирические зависимости для определения критической мутности потока в конкретных условиях. Общим недостатком этих зависимостей является разнородность их структур и показателей степеней влияния различных факторов (от 0,25 до 2,0). Многих авторов это привело к необходимости введения эмпирических коэффициентов, значения которых также изменялись от нескольких сотых долей единицы (0,018 — С.Х. Абальянц) до нескольких единиц (11,0 — Е.А. Замарин), т.е. отличающихся в 500 и более раз. Этот факт не мог не побудить автора статьи к работе над совершенствованием расчетных зависимостей для определения критической мутности потока в отстойниках и каналах. В результате обработки свыше 500 натурных опытов по транспортирующей способности потоков

двадцати оросительных каналов Средней Азии и других водотоков, проходящих по равнинным территориям различных регионов, а также благодаря результатам ученых КНР и США была получена универсальная эмпирическая формула для расчета коэффициента критической мутности потока, несущего мелкие и средние наносы, которая зависит от величины насыщения потока наносами ро [2]:

Кт = 0,25ро07. ° (1)

На основе обработки указанных, в основном натурных, данных разработана новая методика определения мутности потока и гидравлической крупности на выходе из отстойников. Получена кривая и универсальная зависимость для определения критической мутности потока в каналах, несущих большое количество мелких и средних наносов [3]:

Ркр = К т

v

(2)

где юо — средняя гидравлическая крупность наносов в канале, м/с; V и И — соответственно средняя скорость, м/с, и гидравлический радиус потока в канале, м.

Критическую мутность потока в канале ркр. рассчитывают по характерным фракциям I. Зная фактическую мутность отдельной фракции на входе в отстойник рвх., вычисляют мутность каждой фракции, осаждаемой в отстойнике (рБ. = рвх. - ркр.), и мутность всех фракций, подлежащих осаждению в нем, по зависимости р = Е(р . - р .). (3)

1 Б м вх.1 1 кр.17 '

Зная интервал времени Т между промывками, определяют предполагаемый объем отложившихся наносов в отстойнике:

W = 86,4£Грв/ув, (4)

где — расход отстойника, м3/с; Т — время заиления отстойника между промывками, ч или сут; рБ — мутность фракций, осаждаемых в отстойнике, кг/м3; уБ — средняя плотность наносов, осаждаемых в отстойнике, равная 1,2...1,4 т/м3.

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

рости воды в нем, наличия волнения поверхности потока и т.д.

Условия осаждения наносов в отстойниках ухудшаются в ветреную погоду, особенно в период, когда направление ветра примерно совпадает с направлением продольной оси отстойника. В таком случае возникают волны, вызывающие повышенную турбулентность потока и значительно ухудшающие условия осаждения наносов. Чтобы обеспечить допускаемую мутность на выходе из отстойника, для расчетного уровня воды в отстойнике необходимо вводить поправку. По данным А.И. Егорова, уменьшение процента осаждения наносов при ветре можно оценить по зависимости, справедливой при скорости ветра до 10 м/с: р . = р - 1,5у .1Д6, (5)

1 w.l 1 О 7 w.l 7 '

где р — процент осаждения наносов при отсутствии ветра; V . — скорость ветра на высоте 2,0 м от свободной поверхности воды в отстойнике, м/с.

Для получения уравнения динамики осаждения взвеси в потоке примем следующие допущения:

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

крупность и неоднородность наносов характеризуется среднестатистической гидравлической крупностью.

В прямоугольном русле выделим элемент потока длиной X, шириной В = 1,0 и глубиной Н, который движется со скоростью V и имеет на входе мутность рх. Рассмотрим условия осаждения взвеси на длине рассматриваемого участка X. Если в начале участка мутность равна рх, то в конце его она изменится на величину dр и составит рт. Тогда уравнение баланса выпадения наносов за время £ будет иметь следующий вид [2]:

^рги = (рх - р^и^г, (6)

где dр — изменение мутности потока за время dt; & — площадь основания рассматриваемого элемента потока; иВ — скорость осаждения наносов в турбулентном потоке; dрFИ — количество наносов, выпадающих за время dt; (рх - рт)&и^£ — изменившееся количество наносов в потоке.

Из выделенного объема воды за время

dt выпадет часть взвеси, а величина мутности

потока может быть выражена зависимостью

¿Р _ Рх ~Рт

ж

Н

-и в.

Разделив переменные и интегрируя уравнение (7), получим:

---—(8)

Рх-Рт Н к >

_1п(рх -Рт)=~~ + С-

(9)

Рассмотрим граничные условия: при

£ = 0 рх = Po, тогда С = - 1п (ро - рт)-Примем при X = Ь: рх = рВ

_ X

Принимая ¿ — —> получим:

-1п(рв -рт) = ^--1п(р0 -рт). НУ

(10)

Выполнив некоторые преобразования, несложно получить уравнение динамики осаждения наносов:

Ро-Рт

с

Ну

(11)

Рв =Рт+(Ро-Рт)еХР

Отсюда запишем уравнение изменения мутности потока по длине отстойника:

(12)

При использовании уравнения (11) А.Н. Гостунский предлагал принимать скорость осаждения наносов избыточной мутности иВ, равной их среднестатистической гидравлической крупности (и^. По-видимому, это была одна из главных причин значительных отклонений расчетных величин мутности, определенных по зависимости (12), от измеренных. Чтобы учесть этот недостаток, П.В. Михеев принял концепцию взаимосвязи работ турбулентного потока и наносов при их выпадении. Он считал, что взвесь при выпадении производит работу, равную разности нагрузок (роЩо - рТЮо) и выпадает с несколько меньшей скоростью движения избыточной мутности в турбулентном потоке ив, производя работу роив. Это дало ему возможность записать уравнение:

(ро - ртй= роив. (13)

По мнению П.В. Михеева, нагрузка р является условной работой, а фактическая работа должна быть равна:

Уг ~ У о

Р®, (14)

Уг

Уг-Уо

(7)

где

Уг

— коэффициент уменьшения веса

частиц наносов в воде.

Зависимость (13) позволяет выразить значение скорости выпадения наносов в турбулентном потоке:

t/ =PiZPr-Po

(15)

Lcoo (PQ-Pt) Hw po

с —

PB =Рт+(Ро-рт)еХР

L&o (Ро-Рт)

_p0tÛo-pTQ>r откуда uв _

Ри

Избыточная мутность

Рв = Ро - Рт-

В итоге уравнение (19) принимает следующий вид:

Подставив зависимость (15) в (11), получим уравнение динамики осаждения наносов, которое может быть применено для предварительного определения мутности потока в конце отстойника:

Ро-Рт

и* =

(20)

Рв'Рт

Ро-Рт

ро(0о -^рт0>г

Ро-Рт

Подставив уравнение (20) в (11), получим: ' L (роЮо-^рт0)в)Л

= ехр

Ну

(21)

(16)

V "У

Для определения мутности потока на выходе из отстойника эту зависимость можно записать так:

Ро-Рт

Критическую мутность определим по формуле (2) с использованием формулы (1): v

(22)

рт = 0,25р|

\/ЯСОо

(17)

ч Ро .

Это уравнение может быть применено только на предварительной стадии расчетного обоснования динамики мутности потока по следующим соображениям. Во-первых, критическая мутность в рассматриваемом сечении или в конце отстойника не может иметь гидравлическую крупность наносов, равную о^, так как крупность критической мутности при осаждении частиц существенно меньше, чем начальная. Во-вторых, избыточная мутность потока всегда меньше, чем начальная мутность (о,,, т.е. в правой части уравнения (13) некорректно включены также наносы критической мутности, что отмечалось и в других работах. Величину скорости выпадения наносов в турбулентном потоке целесообразно выразить с учетом избыточной нагрузки рИ^В как разности между общей роах, и критической нагрузками рт0>г:

РоСОо -рто>г = Ри*/в, (18)

Для определения коэффициента К вычислим величину гидравлической крупности наносов в выходном сечении по зависимости, полученной в работе [4]:

- _(«Ч>.Г

(0в =

( Л10

Рв

+ 0,04

Рв

Ро

(23)

С этой же целью вычислим гидравлическую крупность критической мутности потока по формуле (6) [2]:

- _(<а>Ро)°,?

СОг =

Рт

0,14

Г Л10 Рт

Ро

+ 0,04— Ро

(24)

Формула (22) позволяет использовать ее для определения окончательного значения мутности по длине потока:

Рв=Рт+(Ро-Рт)х

хехр

— л

L роЮо-^ртС0в Ну Ро-Рт

(25)

(19)

В конечном сечении потока наносы имеют гидравлическую крупность большую или равную гидравлической крупности наносов критической мутности о^. Поэтому выразим через гидравлическую крупность наносов в концевом сечении о^:

СОг =^(0в-

Зная требующуюся мутность рВ, несложно определить длину отстойника:

, ^У(Ро-Рт) ,, Р0-Рт Ь = —=- - 1п-. (26)

РоС0о -/<Грт03в Рв-Рт Полученные уравнения (23), (25) и (26) позволяют проектировщику вычислить характеристики наносов (мутность и гидравлическую крупность), которые поступают в каналы, по длине и в конце отстойника. В связи с тем что мутность на выходе из отстойника отличается от мутности на входе, уравнения дают возможность правильно определить критическую мутность каналов и

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

Расчет мутности и гидравлической крупности наносов в конце отстойника на одном из реальных примеров Каракумского канала-отстойника.

Исходные данные: расход отстойника Q = 30,2 м3/с; длина отстойника Ь = 500 м; среднестатистическая глубина потока в отстойнике и = 1,26 м; среднестатистическая скорость потока в отстойнике V = 0,53 м/с; средневзвешенная мутность наносов на входе в отстойник ро = 3,112 кг/м3; средневзвешенная гидравлическая крупность наносов на входе в отстойник юо = 0,00109 м/с.

Расчеты рекомендуется выполнить в следующей последовательности [2, 5].

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

рТо=0,25р° = 0,25-3,112

о,7 0,53

^1,26-0,00109

= 2,639 кг/м.3

2. Определить приближенное значение мутности потока на выходе из отстойника:

Рв =Рт.о+(Ро-Рт.о)ехР

ЬЮо (Ро~Рт.о)

Нм

Рс

= 2,639 + (3,112-2,639) х

500 0,00109 (3,112-2,639)

хехр

1,26-0,53-3,112

3. Вычислить гидравлическую крупность потока на выходе из отстойника по обобщенной формуле (23):

СОв =

(®>Ро)°

Рв

0,14

Г 1 Л10

Рв

Ро

_ (0,00109-3,112)

0,7

3,057 = 0,000955 м/с.

0,14

1

+ 0,04— Ро

V0

3,057 V 3,112,

+ 0,04

3,057 3,112

4. Определить критическую мутность в конце отстойника для параметров наносов и потока, соответствующих концевому сечению:

/ ! \0,7 V

Рт.в = 0,251 рв ) = 0,25-3,057

0,7 V

НСОв 0,7 0,53

^/1,26-0,000955

= 2,723 кг/м3

5. Вычислить гидравлическую крупность критической мутности потока:

- (®>Ро)°

Юг =

Рт.в

0,14

(0,00109-3,112)' 2,723

= 0,000493 м/с.

0,7

( V0

Рт.в

V Ро ,

0,14

+ 0,04^-Ро

'2,723^°

ч 3,112,

+ 0,04

2,723 3,112

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

^^0000493 Юв 0,000955

7. Вычислить мутность потока в выходном сечении отстойника:

Рв=Рт.В+(Ро-Рт.о)еХР

г Ь РоС0о -^рТ оС0оЛ

НУ Рв-Рт.,

° У

= 2,723+ (3,112-2,723) ехр

500

1,26-0,53

-х

= 3,057 кг/м3.

X

3,112 ■ 0,00109 - 0,516 ■ 2,723 • 0,000955

3,112-2,723

= 2,731 кг/м3.

Выводы

1. Результаты обработки свыше 500, в основном натурных данных, полученных на каналах, в ирригационных отстойниках и других открытых водотоках, несущих мелкие наносы, продемонстрировали, что существующие формулы (около 20), полученные для определения критической мутности потока, носят частный характер, так как стоящие перед этими формулами коэффициенты имеют непостоянные значения. Исследования автора показали, что их значения зависят от величины абсолютной степени насыщения потока наносами.

2. Впервые на основе обобщения сотен натурных опытов, проведенных различными авторами, получена универсальная формула для вычисления величин критической мутности потока, апробированная на 19 отечественных и зарубежных водотоках в следующих пределах: Я = 0,1...2,6 м; V = 0,1...1,3 м/с, I = 0,00001...0,001; 0,0003...0,006 м/с; ро = 0,5...30,0 кг/м3 и й < 0,25 мм.

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

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

Список литературы

1. Беликов, В. В. К определению оптимальных параметров отстойников ГЭС [Текст] /

B. В. Беликов, А. М. Прудовский, В. П. Сере-давкин // Безопасность энергетических сооружений. - Вып. 12. - М.: ОАО «НИИЭС», 2003. -

C. 166-183.

2. Кавешников, Н. Т. Разработка методов расчетного обоснования, конструирования и эксплуатации сооружений комплексных низконапорных гидроузлов [Текст] / Кавешников Николай Трофимович : дис. ... д-ра техн. наук ; защищена 02.02.1998 ; утв. 05.05.1998. - 529 с.

3. Кавешников, Н. Т. Критическое состояние насыщения потока наносами в открытых руслах [Текст] / Н. Т. Кавешников // Гидротехническое строительство. - 1997. - Вып. 1. -С. 20-23.

4. Кавешников, Н. Т. Расчетное обоснование гидравлической крупности наносов на выходе из отстойника [Текст] / Н. Т. Кавеш-ников // Гидротехническое строительство. -1997. - Вып. 3. - С. 53-56.

5. Кавешников, Н. Т. Методы расчета осаждения наносов и параметров мутности потока в отстойниках [Текст] / Н. Т. Кавешников // Гидротехническое строительство. - 1997. -Вып. 10. - С. 26-30.

6. Kaveshnikov, N. T. Falluja Barrage Water intake Structure Settling basin Hidraulic studes on the Euphrates River in Republic of Irag. [Text] / N. T. Kaveshnikov, N. P. Rozanov // Report and recommendation contract 21305 v/o Selkhospromexport. - Baghdad, 1981. - 94 р.