автореферат диссертации по строительству, 05.23.07, диссертация на тему:Фильтрационно-температурный режим системы "плотина-основание"



На правах рукописи

0 3 СЕН 2009

АНИСКИН Николай Алексеевич

ФИЛЬТРАЦИОННО-ТЕМПЕРАТУРНЬШ РЕЖИМ СИСТЕМЫ «ПЛОТИНА-ОСНОВАНИЕ»

Специальность 05.23.07 - Гидротехническое строительство

Автореферат диссертации на соискание ученой степени доктора технических наук

Москва 2009

003475883

Работа выполнена в Государственном образовательном учреждении высшего профессионального образования Московском государственном строительном университете (ГОУ ВПО МГСУ)

Официальные оппоненты:

доктор технических наук Беликов Виталий Васильевич

' доктор технических наук, профессор Анахаев Кошкинбай Назирович

доктор технических наук Покровский Геннадий Иванович

Ведущая организация -

ОАО «Институт Гидропроект»

Защита диссертации состоится « 29 » сентября 2009 года в 15 час. 30 мин. на заседании диссертационного совета Д 212.138.03 при ГОУ ВПО Московском государственном строительном университете по адресу: г. Москва, ул. Спартаковская, д.2/1, ауд. 212.

С диссертацией можно ознакомиться в библиотеке ГОУ ВПО Московского государственного строительного университета.

Автореферат разослан « » ¿"¿¿^2009 года

Ученый секретарь диссертационного совета

Общая характеристика работы.

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

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

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

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

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

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

Целью диссертационной работы явились:

• разработка методики решения фильтрационных, температурных и температурно-фильтрационных задач применительно к системе «плотина-основание» с учетом пространственности, неоднородности и нелинейности свойств материалов плотины и основания, возможной анизотропии и нестационарности процесса;

• создание численных фильтрационных и температурно-фильтрационных моделей гидротехнических сооружений и их оснований;

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

Основные задачи исследований:

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

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

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

Научная новизна работы состоит:

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

■ в полученных результатах численных исследований систем «плотина-основание» применительно к реальным гидротехническим объектам с выводами о работе сооружения и рекомендациями по его конструктивным элементам.

Достоверность научных результатов подтверждена:

• результатами решения ряда тестовых задач, имеющих точное аналитическое решение;

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

• сравнением результатов численных исследований Юмагузинского, Сангтудинского, Курейского гидроузлов с данными натурных наблюдений.

Практическое значение работы и полученных результатов заключается:

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

■ в разработке комплексов вычислительных программ «Ш/ГЯ» и «ТЕ11М1С» на базе метода конечных элементов для решения фильтрационных, температурных и совместных задач;

■ в создании численных моделей фильтрационных режимов плотин совместно с основаниями применительно к следующим гидроузлам: Бурейскому, Камбаратинскому, Юмагузинскому, Сангтудинскому, Загорской ГАЭС-2, температурно-фильтрационной модели и прогнозу работы для плотины Курейской ГЭС, в решении задач по прогнозу нестационарных фильтрационных режимов ограждающих дамб Северной и Мезенской ПЭС;

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

Апробация работы Основные положения методики и результаты исследований докладывались на научно-технических конференциях:

- на научно-технических конференциях МИСИ-МГСУ;

на Всесоюзном координационном совещании «Прочность и температурная трещиностойкость бетонных гидротехнических сооружений при температурных воздействиях, г. Нарва, 1989;

- на второй, третьей и четвертой научно-технических конференциях «Гидроэнергетика, новые разработки и технологии» в 2006, 2007 и 2008 годах в г. Санкт-Петербурге.

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

исследования по Бурейскому, Камбаратинскому, Юмагузинскому, Сангтудинскому, Загорской ГАЭС-2, Курейской ГЭС, решены задачи по прогнозу нестационарных фильтрационных режимов ограждающих дамб Северной и Мезенской ПЭС. Сделан анализ влияния технологических факторов на температурный разогрев бетонного массива при его возведении. Проведены численные расчеты температурных режимов гравитационных плотин из укатанного бетона Хлонг Та Дат (Тайвань) и Шон Ла (Вьетнам) в строительный и эксплуатационный периоды, позволившие дать оценку эффективности принятых параметров возведения.

Публикации. Основные положения диссертации опубликованы в 16 научных работах, в том числе 9 работ в ведущих рецензируемых научных журналах и изданиях, рекомендованных ВАК: в журнале «Гидротехническое строительство» -6, сборнике «Вестник МГСУ» -2, в журнале «Строительные материалы, оборудование, технологии XXI века» -1. Основные положения методики и некоторые результаты исследований помещены в учебник «Гидротехнические сооружения» (издательство АСВ, 2008 год).

На защиту выносятся:

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

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

■ численные математические модели фильтрации и на их основе анализ работы сооружений и протавофильтрационных элементов, выполненные для Бурейского, Камбаратинского, Юмагузинского, Сангтудинского гидроулов, Загорской ГАЭС-2, температурно-фильтрационная модель и прогноз работы на перспективу плотины Курейской ГЭС, результаты решения задач по прогнозу нестационарных фильтрационных режимов ограждающих дамб Северной и Мезенской ПЭС.

Структура и объем диссертации. Работа состоит из введения, пяти глав, заключения и списка литературы. Объем диссертации составляет 416 страниц, из них 244 страницы текста, список литературы состоит из 272 наименований.

Основное содержание работы.

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

Теория фильтрации, началом которой можно считать открытие в 1856 году французским инженером Г.Дарси (Darcy, 1856) закона движения фильтрационного потока, нашла свое развитие в работах Буссинеска (1904 )Дюпюи Ж. (1857)Жуковского Н.Е. (1889), Павловского Н.Щ1922, 1956), Маскета М. (Muskat 1937, 1949), Лейбензона Л.С.(1947, 1953), Форхгеймера

Ф. (1901), Щелкачева В.Н. и Лапука Б.Б.(1949), Полубариновой-Кочиной П.Я. (1952), Шейдеггера А.Э.( 1960), Чарного И.А. и т.д.

Огромные масштабы гидроэнергетического и гидромелиоративного строительства в СССР в период с 1920 до 1980 годы потребовали от ученых решения многих задач напорной и безнапорной фильтрации. В начале этого периода внимание исследователей было сконцентрировано на решениях так называемой плоской задачи теории установившейся фильтрации. Решений фильтрационных задач в пространственной постановке в те времена практически не существовало, и в тех случаях, когда фильтрация носила пространственный характер, решение сводилось к плоской задаче с использованием разумной схематизации и упрощений задачи.

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

Во всех вышеперечисленных работах пористая среда рассматривалась как однородная и изотропная. Однако в природе достаточно часто встречаются грунты с анизотропными свойствами. В работах Дахлера Р. (1933), Шаффернака Ф.(1933), Ризенкампфа Б.К. (1938), а затем Аравина В.И. (1937) и др. рассматривались задачи фильтрации в анизотропной пористой среде, решение которых сводится к решению в фиктивной изотропной пористой среде с последующим возвращением в исходную область. Задачи безнапорной фильтрации в анизотропной однородной пористой среде рассмотрены в работах Полубариновой-Кочиной П.Я.(1940), а затем Михайлова Г.К. (1951).

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

Лишь частные результаты были получены теоретическими методами при решении задач нелинейной фильтрации, когда скорость фильтрации и градиент напора связаны нелинейной зависимостью (например, по зависимости Смрекера). Такие задачи могут быть решены методом последовательных приближений линейных задач, или, по предложению Христиановича С.А.(1940) методом ЭГДА.

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

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

фрагментов является изучение напорной и полунапорной фильтрации в основаниях бетонных плотин. Этому вопросу посвящены работы Павловского Н.Щ1936), Давидовича В.Щ1937), Аравина В.И. и Нумерова С.Н. (1947,1955), Недриги В.Щ1956) и др.

Метод фрагментов подготовил почву для разработки широко применяемого и сегодня для приближенных оценок фильтрации метода фильтрационных сопротивлений, который был впервые сформулирован в 1951 году Нумеровым С.Н. и, почти одновременно с ним Борисовым Ю.П. (1951). Метод фильтрационных сопротивлений использован для решения задач фильтрации в работах Нумерова С.Н.(1953,1954), Аравина В.И.. Нумерова С.Н. (1953,1955), Павловской Л.Н. (1964), Шестакова В.М. (1964), Чугаева Р.Р. (1955) и др.

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

Поэтому помимо теоретических методов, параллельно шло развитие экспериментальных или модельных исследований. Одно из направлений экспериментальных исследований фильтрации - моделирование процесса путем пропускания вязкой жидкости в щели между двумя параллельными пластинками (так называемые щелевые лотки). Эта методика была использована в работах Аравина В.И., Замарина Е.А., Каранфилова Т.С. Следует отметить, что моделирование фильтрации с помощью вязкостной аналогии достаточно широко использовалось не только в СССР, но и в других странах: США, ГДР, Румынии и др.

Получили очень широкое развитие аналоговые методы. Наиболее широкое применение нашел метод электро-гидродинамических аналогий (метод ЭГДА), основанный на математической аналогии течения электрического тока в электропроводном материале и фильтрации жидкости в пористой среде.

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

Использование одного из численных методов - метода конечных разностей (МКР) при фильтрационных расчетах было впервые предложено в 1924 году Николаи Е.Л. (задолго до появления ЭВМ). Большинство последовавших за этим работ с использованием МКР были посвящены численному решению уравнения Буссинеска.

Метод математического моделирования получил в последние десятилетия широкое распространение, что связано с развитием численных методов и вычислительной техники. Если на первых этапах развития и применения численных методов превалировал МКР, то позднее наибольшее распространение получил метод конечных элементов (МКЭ). Преимущество

последнего метода над МКР объясняется его вариационной сущностью, большими возможностями при моделировании областей сложной геометрии и неоднородностью свойств материалов, возможностью использовать элементы различных типов. Основной недостаток МКЭ - большой порядок разрешающей системы уравнений, что вызывало определенные сложности при расчетах сложных систем в 70-80 -х годах прошлого столетия. В связи с этим появился метод «локальнных вариаций» (Черноусько Ф.Л., Банничук Н.В.), нашедший применение при решении статических и сейсмических задач для бетонных и грунтовых плотин.

Дальнейшим развитием численных методов стало появление метода суперэлементов (МСЭ), в котором расчетная область разбивается на отдельные подструктуры, каждая в свою очередь разделяется на отдельные части и т.д., пока не образуются достаточно малые и простые по форме «базисные элементы». Другой вариант совершенствования численных методов, также уменьшающий порядок разрешающей системы - метод «граничных элементов» (МГЭ).

В настоящее время появилось множество комплексов, позволяющих решать набор задач (статических, динамических, фильтрационных, температурных и т.д.) в плоской и пространственной постановках. В качестве примеров использования таких комплексов можно назвать комплекс «DRENA» (ВНИИГ им. Б.Е.Веденеева), «CRISP-PATH FEM» (Чехия, д-р М.Долежалова), расчетный комплекс «ANSYS». Достаточно широкое применение нашел программный комплекс «Visual MODFLOW», разработанный канадской фирмой Waterloo Yelrogeologic, использующийся для создания пространственных геофильтрационных моделей, комплекс вычислительных программ «Земля» (институт «Гидропроект») и др..

Тем не менее, решение пространственных фильтрационных задач для системы «плотина-основание», имеющей, как правило, сложную геометрию, с учетом множества факторов (неоднородность, анизотропность свойств материалов, нелинейность и нестационарность процесса, наличие в системе трещин, разломов, зон повышенной проницаемости и т.д.) является весьма трудоемкой задачей.

Аналогичное развитие претерпела современная теория теплопроводности и методы решения температурных задач. Эволюция методов решения температурных задач аналогична развитию фильтрационных методов. В первой половине XX века широко используются аналитические и графоаналитические методы. С середины века интенсивное развитие гидротехники и необходимость решения практических задач в СССР вызывает «бум» модельных исследований и аналоговых методов. Появление вычислительной техники и ее развитие во 2-ой половине XX столетия вызывают разработку численныых методов применительно к решению температурных задач.

Вопросы, связанные с определением температурного режима бетонных конструкций с применением различных методов рассматриваются в работах Александровского C.B., Белова A.B., Васильева П.И., Гутмана С.Г., Дзюбы К.И., Ламкина М.С., Маслова Г.Н., Орехова В.Г., Фрадкиной Н.И., Фрида С.А., Чилингаришвили Г.И., Цыбина A.M., Fanelli M. и Giuseppetti G. и многих других авторов.

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

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

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

Основное дифференциальное уравнение в частных производных неустановившейся фильтрации в пространственной постановке записывается в виде (уравнение Пуассона):

дН

■Л—- = 0, (1)

о t

где H=f(x,y,z,t) - искомая напорная функция в расчетной области, изменяющаяся во времени t\ Кх, Ку, Кг - коэффициенты фильтрации по направлениям координатных осей X, Y,Z; ц- коэффициент водоотдачи грунта.

Прямое численное решение уравнения (1) является весьма трудоемкой задачей. В настоящей работе отыскание напорной функции происходит в ходе минимизации некоторого, специальным образом подобранного функционала Ф, для чего и создан специальный расчетный аппарат, основанный на синтезе МКЭ и МЛВ (метод локальных вариаций).

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

д ( д Н] д г dlî д Г д Н]

— - +- К + — Kz-

д х s Ч д х ^ ду V У ду, dz 1

На основании теоремы Эйлера функционал Ф (для случая трёхмерной задачи) имеет вид:

Ф = 1II

V

К.

дН

+ К,

дН

¥

+ К,

дН

аН

+ /л—Н \dxdydz, (2) &

Решение фильтрационной задачи сводится к отысканию напорной функции Н(х, у, 2, которая была бы непрерывной внутри замкнутой области V, удовлетворяла бы начальным и граничным условиям, и доставляла бы минимум функционалу Ф (2).

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

ч2" д Н

+ /л--Н^хёуёгч-

N 1 [ Ггн^ 2 { д Н^ 2

Ф= ПН к* + к2

п=1 V 2 У \дх)

дХ

к М (дУр 2 (дН)

+ ЕЛяНШП + 111 { дг})

к-Щ ш=1 5

(3)

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

направлению осей локальной системы координат.

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

У=-Кф-Г, (4)

где п - показатель степени, определяемый по экспериментальным данным (изменяется в интервале от 1,0 до 0,5); при п = 1,0 выполняется закон Дарси; при п = 0,5 имеем турбулентный закон фильтрации; если 1,0< п< 0,5 имеет место переходный режим от ламинарного к турбулентному.

В общем случае составляющие скорости фильтрации по координатным осям ¡=х, у, г могут быть представлены в следующем виде (формула Смрекера): V, = =- К, У,""7 У, = - К" У, (5)

где К" - нелинейные коэффициенты фильтрации по осям координат ;'=дг, у, г.

Задача в нелинейной постановке решается методом последовательных приближений. Первоначально решается линейная фильтрационная задача, из которой определяются величины градиентов. Затем, в соответствии с зависимостью (5) определяются нелинейные коэффициенты фильтрации К"= К, У/1"' и процесс минимизации возобновляется. Итерационный процесс заканчивается в случае неизменности величин К" на текущем и предыдущем шагах (с заданной точностью).

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

Эх

а

Эх,

а

а г а'

+ — \аг — дг V Эг

Ц ЭЭ ср Эт

а

дх'

эу1у эу;

где Г = Г (х, у, г,г) - искомая температурная функция; ах, ау, а2 - температуропроводность материала по направлению координатных осей X, Г,т - время; с - удельная теплоемкость бетона; р - плотность бетона; Э - удельное тепловыделение цемента; Ц - расход цемента в 1 м3 бетона.

Решение дифференциального уравнения (6) аналогично решению основного дифференциального уравнения теории фильтрации (1) сводится к минимизации следующего объемного интеграла при заданных начальных и граничных условиях:

Эх

а

1

а ц ээ

Эт ср Эт

I \ ах • ау • Аг +

+ ячкю+ я-р(1-0 «ю, (7)

п, п2 2 '

где Пь Пг - поверхности расчетной области, на которой выполняются соответственно

граничные условия 2 (задание на поверхности удельного теплового потока я) и 3

(теплообмен с окружающей средой, имеющей температуру родов.

Решение совместной фильтрационно-температурной задачи при

равенстве температуры твердого скелета грунта и температуры

фильтрующей жидкости в любой точке области фильтрации сводится к

решению известного уравнения Фурье-Кирхгофа:

а

Эт:

д1г

--1- а

:Эх2 '

дЧ

э2/

ду ог

СгУг

св ув

& а

. —ну,,--1- у,

дх Уду 2

(8)

1(х,у,г,т) - температура твердого тела (материала плотины или основания); т - время; ах, ау, а. - коэффициенты температуропроводности твердого тела, насыщенного водой, по осям х, у, г; V*, V,,, - компоненты скорости фильтрации по осям х ,у ,г; С„ Ст - удельная объемная теплоемкость воды и твердого тела; у„ уг - плотности воды и твердого тела.

Решение дифференциального уравнения (8) равносильно минимизации функционала:

\ 2

+ 110X1+¡£/3(1

о г, ^

— + я„ —

+—/+(?, раЛУ \dxdydz + 9т

(9)

где (у,£гаЛ)= 1

а

.Уча,

а

ду

а

дг

скалярное произведение векторов

фильтрационной скорости и температурного градиента.

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

В качестве основной элементной базы используются объемный 8-узловой и плоский 4-узловой элементы произвольной формы, позволяющий достаточно точно аппроксимировать исследуемую область. Для

моделирования трещин используются 4-узловой элемент произвольной формы (пространственная задача) и 2-х узловой линейный элемент (плоская задача).

На основе изложенной методики были составлены программные комплексы для ЭВМ расчета фильтрации "ИЬТ11" и расчета температурного режима «ТЕИ-МГС». В состав программных комплексов входят несколько подпрограмм: подпрограмма подготовки исходной информации; подпрограмма проверки исходной информации; подпрограмма минимизации функционала в расчетной области; подпрограмма обработки к печати полученных результатов; подпрограмма визуализации полученных результатов.

Сетки разбивки на конечные элементы при современном уровне развития вычислительной техники могут включать до 800000 узлов для объемных задач и до 100000 узлов для плоских задач (при необходимости количество узлов может быть увеличено.

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

Грунтовая плотина Юмагузинского гидроузла- с центральным ядром из суглинка и боковыми призмами из гравийно-галечникового грунта имеет максимальную высоту 70 метров. Плотина возведена в сложных с точки зрения фильтрации инженерно-геологических условиях. Основание плотины в русловой части и, особенно, на правобережной пойме характеризуется крайней неоднородностью фильтрационных свойств: коэффициенты фильтрации различных слоев грунта изменяются от 0,3 до 90 м/сут. Аллювиальные отложения в пределах правобережной поймы представляют собой напластование грунтов с различными фильтрационными характеристиками суммарной глубиной до 40 метров. Целью проведенных исследований было сопоставление плоских и пространственных решений, проведение анализа влияния таких факторов, как фильтрационная неоднородность, размеры противофильтрационных элементов. Исследование фильтрационного режима плотины и основания Юмагузинского гидроузла проводились в несколько этапов с 1999 по 2006 год.

Первый этап фильтрационных исследований (1999-2000 год) проходил в период начальной стадии проектирования грунтовой плотины Юмагузинского гидроузла и начала ее строительства. Для оценки эффективности принятых в проектном варианте противофильтрационных элементов (ПФЭ) и их необходимости были проведены расчеты в плоской и пространственной постановках для двух конструктивных вариантов: с устройством ПФЭ в основании и без них. Проектом предусматривалось устройство 2-х рядной инъекционной противофильтрационной завесы под плотиной и в бортовых примыканиях, дополнительной 5-ти рядной завесы в основании ядра в русловой части створа и «стены в грунте» на участке залегания аллювиальных отложений на правобережной пойме.

На рис Л, а представлено плановое сечение расчетной области фильтрации в основании под плотиной со следами эквипотенциальных

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

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

Была составлена упрощенная модель основания плотины (с объединением слоев грунта с близкими значениями коэффициента фильтрации). Выделено шесть слоев грунтов, для каждого из которых на основании инженерно-геологических данных определялись максимальные, минимальные и средние коэффициенты фильтрации. Так, для галечника эти значения соответственно равны 100, 10 и 55 м/сут, для суглинков и глин -0,4, 0,01 и 0,205 м/сут.

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

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

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

<5г = 1,5 + 0,051X1+0,164 Х2 + 0,61Х5 + 0,052 Х2Х5 (10)

При этом величина суммарного расхода изменяется от СЬ =0,72 м3/сек (при минимальных коэффициентах фильтрации) до Ох =2,5 м3/сек (при

Рис. 1. Фильтрация в плотине в основании Юязгузянского гидроузла, з- распределение линий равного напора б горизонтальном оеченви по отаетке 130.0; б-фзльтрашя в левобережном примыкании иосж>Б£шш плотины;

максимальных коэффициентах фильтрации). Как видно из уравнения (10), наибольшее влияние на величину суммарного расхода оказывает фактор Х5 (коэффициент фильтрации известняка в пределах Кф = 5,0+30,0 м/сут) и фактор Х2 (коэффициент фильтрации суглинков и глин в пределах Кф = 0,014-0,4 м/сут).

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

материалы одновременно рассматриваются или с максимальными или с минимальными коэффициентами фильтрации. Рассматривались следующие факторы: X] - уровень воды в верхнем бьефе (Х| =-1 соответствуют отметке НПУ 260,0; X] =+1 - отметке ФПУ 270,0); Х2 - значения коэффициентов фильтрации грунтов, залегающих в основании плотины (Х2 =-1 - все значения коэффициентов фильтрации минимальны; Х2 =+1 - все значения коэффициентов максимальны). Уровень нижнего бьефа принимался на отметке 210,0. В качестве откликов, как и в предыдущих исследованиях, рассматривались величины фильтрационного расхода и максимального градиента напора по оси противофильтрационного устройства. Были получены следующие функции откликов:

-для расхода: <} = 1,87 + 0,2Х, + 1,125Х2 + 0,115 X, Х2 (11)

- для градиента: = 18,375 + 1',625Х,+ 9,125Х2 + 0,875 X, Х2 (12)

Можно отметить следующие различия в результатах решения задач в плоской и пространственной постановках: величина суммарного фильтрационного расхода, полученного из пространственной задачи несколько выше величины, полученной из плоского решения. При максимальных значениях коэффициентов фильтрации грунтов основания величина суммарного расхода в объемной задаче составляет 3,31 м3/сек, в плоской - 2,5 м3/сек. При минимальных значениях коэффициентов фильтрации эти значения соответственно равны 0,83 м3/сек и 0,72 м3/сек. Увеличение суммарного расхода в пространственной задаче получено за счет обходной фильтрации в бортовых примыканиях.

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

Последующие исследования (2001 -2006 г.г.) были вызваны уточнением данных по инженерно-геологическому строению основания, свойствам грунтов тела плотины и конструкции противофильтрационных элементов в основании и бортах сооружения. Была предпринята попытка создания более совершенной и достоверной численной фильтрационной модели, учитывающей уточненные инженерно-геологические условия, реальные параметры и конструкции как самой плотины и других сооружений, входящих в состав гидроузла (эксплуатационного водосброса, ГЭС) и противофильтрационных устройств в основании. Совершенствование программных и вычислительных возможностей также позволило значительно увеличить частоту аппроксимационной сетки метода конечных элементов (в 2001 году сетка МКЭ пространственной задачи состояла всего из 9729 узлов и 8398 элементов. Количество узлов и элементов новой сетки МКЭ в каждом

поперечном оси плотины сечении равно соответственно 2721 и 2595. Общее количество узлов и элементов аппроксимационной сетки составило соответственно 179586 и 168675.

Результаты решения фильтрационной задачи сравнивались с натурными данными. Это сравнение показало, что математическая геофильтрационная модель достаточно адекватно отражает реальную ситуацию: величины напоров, полученные в результате пьезометрических измерений достаточно близки (за небольшими исключениями) к расчетным в большинстве пьезометров: в 10 из 15 относительная разница не превышает 6,8%, в 6 пьезометрах относительная разница находится в пределах 5%.

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

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

Согласно рассмотренному первоначальному варианту в русловой части створа и на левом берегу в месте расположения современных аллювиальных пород (галечник с включением валунов с гравийно-песчаным заполнителем) с большим коэффициентом фильтрации-Кф=85 м/с предусмотрено устройство «стены в грунте». Помимо этого по всему створу устраивается площадная цементация под ядром и глубинная противофильтрационная завеса.

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

На основании проведенных на первом этапе численных исследований были сделаны следующие выводы. «Стена в грунте» в русловой части плотины (рис. 2,а) и левобережном примыкании в рассмотренном варианте работает неэффективно. Также малоэффективна глубинная цемзавеса в правобережном примыкании. Был сделан вывод о целесообразности

изменения конструкций противофильтрационных элементов в русловой части и бортовых примыканиях Фильтрационный расход через выделенную область в основании грунтовой плотины и ее береговые примыкания составляет ~ 5.7 м3/сек.

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

По результатам натурных наблюдений при уровне верхнего бьефа 554,6 фильтрационный расход на 05.05.08 составлял 4,2 м3/с, а на 30.06.08 - 4,0 м3/с, что несколько выше расхода, полученного по численной геофильтрационной модели «Института Гидропроект». Данная модель дала соответственно 3,4 м3/с и 3,45 м3/с, причем на долю обходной фильтрации в бортах приходится 1,13 м3/с и 1,03 м3/с.

На основании модели МГСУ 2006 года при уровне верхнего бьефа 554,6 величина расхода и его распределение по створу может быть оценена следующим образом: суммарный фильтрационный расход составляет примерно 4,39 м3/с из них на долю обходной фильтрации приходится ~ 1,0 м /с. Как видно, прогнозные значения по фильтрационным расходам, полученные в МГСУ в 2006 году несмотря на ограниченный объем информации и произошедшие изменения в конструкции и оценке фильтрационных свойств грунтов плотины и основания, достаточно хорошо корреспондируются с величинами натурных и модельных исследований «Института Гидропроект». На основании исследований, проведенных на участке создания «стены в грунте» (бурение скважин с отбором керна, сейсмическое просвечивание между скважинами, поверхностного сейсмопрофилирование и резистивиметрия на выборочных участках) был установлен факт отсутствия сплошности в противофильтрационном элементе. Были установлены «окна» в едином фронте противофильтрационного элемента на русловом участке. Для имитации реальных условий была составлена численная модель руслового фрагмента с учетом имеющихся данных по качеству противофильтрационных элементов. Моделируемая часть руслового фрагмента рассматривалась от сечения ПК 2+10 до сечения ПК 4+05. Была составлена сетка разбивки расчетной фильтрационной области на конечные элементы. При этом моделировалась зона повышенной проницаемости или зона разрыва между скважинами струйной и глубокой цементации (см. рис. 4).

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

ЯНПУ 531,50-

ПКЗ+10Й

V 51-3,00

300 до

9НПУ 571.50

уз?«,»

? 517.00

■ ЗООЙО

РнсЗ Распределение линий равных капере® по сеч«яих>9-9 по резупь таг мл рейхе пространственной фильтрационной задачи для плотккь: СангтудинсхойГЭ

ПК4+5Я

ПК2+10Л

в стенке

РисС «к а конечных мементов русл ос «го фрагмента (сечение по сн «стекые грунте»)

Рис 2 Распределение линий равного напора в основании грунтовой плотины Санпудинского гидроузла (пространственная задача) I - площади« цементация. 2-гпубиним ц««м«са. З-деггреесиоктя псеерзенот, 4- «стен« е грунт«», водепдоедезо; 6- контур ппоткны

галечник обладает большей проницаемостью, заложенной в фильтрационные модели МГСУ и «Института Гидропроект», что может быть вызвано и процессами фильтрационных деформаций, произошедших в подземном контуре сооружения.

Учитывая, что к моменту исследований фильтрационный расход упал до ~ 2 м3/с вероятно, был сделан вывод, что грунт «в окне» имеет коэффициент фильтрации ориентйровочно 10-^-20 м/сут, так как обходную фильтрацию можно оценить в 1 м3/с.

Дальнейшие фильтрационные исследования были проведены с целью создания более совершенной пространственной геофильтрационной модели системы «плотина- основание». Пространственная модель включала в себя численную модель руслового фрагмента (см. выше), в правобережном примыкании моделировалось устройство противофильтрационной глубокой завесы в соответствии с проектом, более детально воспроизведен рельеф в нижнем бьефе с учетом устройства отводящик каналов ГЭС и водосбросов. Была существенно сгущена сетка МКЭ: количество узлов и элементов сетки было увеличено соответственно до 451584 и 433675. Расчеты проводились для двух уровней верхнего бьефа: отметки 554,6 и отметки НПУ 571,5.

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

На рис.3 дана картина изменения фильтрационного напора в горизонтальном сечении расчетной области на отметке 500,0. Как видно, под плотиной основное падение напора происходит на протифильтрационной завесе. В берегах основным сопротивление является пласт глинистого грунта (Кф=0,003 м /с).

Полученные на численной модели результаты сравнивались с результатами натурных данных (показания установленных пьезометров) и численной модели института «Гидропроект». Сравнение показало в целом удовлетворительное совпадение с натурными данными по состоянию на 30.06.2008 и данными модели «Гидропроекта». В большинстве пьезометров (в 13 из 16) расхождения с натурой при величине действующего напора 46 м находятся в пределах от 0,03 м (пьезометр 117) до 4,3 м (пьезометр ПО 1-3.7). Несколько больше расхождения в пьезометрах ПО 1-4.10 (5,16 метров) и 123 (7,1 метра). Самое большое расхождение в 22,2 м, как и в модели института «Гидропроект» (23,5 м), получено в пьезометре НП 2.2, расположенном в русловом аллювии в основании ядра плотины. Высокие значения уровня в пьезометре НП 2.2 могут быть вызваны заниженной водопроницаемостью руслового аллювия или завышенной проницаемостью инъекционной завесы. Возможно, это связано с возникновением в аллювиальном слое закольматированной зоны с меньшими значениями коэффициента фильтрации.

Параллельно исследованиям в 2009 году проводились работы по снижению фильтрационного расхода в русловой части створа. Были

продолжены работы по устройству скважин трехрядной противофильтрационной завесы со стороны верхнего бьефа от «стены в грунте». Кроме того, было проведено «запесочивание» проницаемого слоя, располагающегося в русловой части под ядром плотины. В результате этих работ величина фильтрационного расхода к апрелю 2009 года была снижена с 4 м3/с до 1,2 м3/с, причем по оценкам ~0,8 м3/с приходится на русловую часть и -0,4 м3/с - на оставшуюся часть створа. Таким образом, можно утверждать, что в проведенных численных исследованиях была завышена проницаемость «бухарских» известняков (пачка № 23-25 с коэффициентом фильтрации Кф= 50 м/сут выше отметки 460,0 и Кф= 5 м/сут ниже этой отметки по всему створу) и, возможно, глубокой противофильтрационной завесы. Была проведена очередная калибровка численной модели и проведены расчеты при значении коэффициента фильтрации «бухарских» известняков Кф= 15 м/сут выше отметки 460,0. Это привело к снижению фильтрационного расхода до ~1,2 м3/с, что соответствует наблюдаемому. Дальнейшее совершенствование численной модели может быть связано с локальными уточнениями характеристик пород на основании натурных данных. Конечно, только численные исследования могли позволить провести настолько подробные исследования в столь сложных геологических условиях.

Часто негативные последствия фильтрации сказываются на устойчивости грунтовых склонов. В частности, весьма остро стоит вопрос об устойчивости склонов, на которых располагаются водоводы ГАЭС. Исследование фильтрационного режима склона и ограждающей дамбы ГАЭС проведено на примере сооружений Загорской ГАЭС-2. Пристальное внимание к фильтрационному режиму склона вызвано проблемами, возникшими при строительстве Загорской ГАЭС-1.

Грунты, складывающие верхнюю часть склона характеризуются достаточно низкими коэффициентами фильтрации (порядка 10'2 м/сут и ниже). Влияние колебаний уровней бассейнов, происходящие за достаточно короткие промежутки времени, будет сказываться незначительно и только вблизи бассейнов, не сказываясь на общем фильтрационном режиме всего склона. Вследствие этого решалась задача установившейся фильтрации для наихудшего случая: в верхнем бассейне максимальный уровень НПУ 266,5 метров, в нижнем бассейне УМО 152,5 метра (величина расчетного фильтрационного напора составляет 114,0 м).

По результатам фильтрационных расчетов склона водоводов можно отметить: движение фильтрационного потока из верхнего аккумулирующего бассейна имеет ярко выраженный пространственный характер. Это положительно сказывается на фильтрационном режиме склона, что приводит к значительному снижению депрессионной поверхности по сравнению с результатами плоской фильтрационной задачи. Четко просматривается влияние оврага №2 и других понижений рельефа, играющих роль естественного дренажа. Учет пространственности также значительно снизил фильтрационное давление на водоприемник ГАЭС.

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

проведены исследования фильтрационного режима дамбы верхнего аккумулирующего бассейна Загорской ГАЭС-2 в плоской постановке с учетом возможной отсыпки боковых призм дамбы слоями различной проницаемости, что вызвано стремлением к всепогодной технологии укладки грунта. При этом рассматривалось влияние таких факторов, как толщины отсыпаемых слоев, их протяженность, расположение водосборного коллектора в низовой призме дамбы. Решалась задача установившейся фильтрации при уровне воды в верхнем бассейне ГАЭС на уровне НПУ 266,5 метров для 2-х сечений дамбы: ПК 26+00 и ПК+11.00. Следует отметить, что рассматривался наихудший с точки зрения фильтрации случай: уровень воды верхнего бассейна ГАЭС рассматривался на отметке НПУ.

Конструкция дамбы в сечении ПК 26+00 представляет собой следующее. Вдоль верхового откоса дамбы выполнен экран из суглинистого грунта толщиной 5,0 метров с коэффициентом фильтрации Кф=0,005 м/сут. По оси дамбы устроен вертикальный дренаж толщиной 2,0 м из грунта с коэффициентом фильтрации Кф=1,0 м/сут. Вдоль основания под низовой призмой и частично под верховой уложен горизонтальный дренаж толщиной 2,0 м из грунта с коэффициентом фильтрации Кф=1,0 м/сут. Горизонтальный дренаж сопрягается с вертикальным и заканчивается водосборным коллектором. В проведенных расчетах варьировалось удаленность коллектора от оси дамбы. Верховая и низовая призмы дамбы представляют собой слоистую конструкцию, уложенную чередующимися слоями одинаковой толщины, выполненными из грунтов с различной проницаемостью: с коэффициентами фильтрации Кф=0,5 м/сут (укладываемый в любых погодных условиях) и Кф=0,01 м/сут. Варьировалась толщина укладываемых слоев (1 и 2м).

Аналогичные исследования были проведены для сечения дамбы в ПК 11+00, имеющей аналогичную конструкцию. Исследования

фильтрационного режима дамбы ВАБ ГАЭС позволили сделать вывод о возможности послойной укладки грунта различной проницаемости при заданных характеристиках грунтов и нормальной работе противофильтрационных элементов.

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

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

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

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

Камбаратинская плотина возводится направленным взрывом в сложных топографических условиях. Каньон плотины имеет сложную форму - в борту правого берега имеется естественная выемка, в пределах которой длина плотины по гребню достигает 400 метров. Затем по течению реки каньон сужается до 200 метров (на отметке НПУ) и после сужения наблюдается вновь плавное расширение каньона в сторону нижнего бьефа (рис. 5, а). На основании прогноза изменения геотехнических свойств грунта по высоте в рассмотренной плотине было выделено пять зон с различными коэффициентами фильтрации и показателями степени от Кф = 2.5 см/сек ( а = 0.6 ) в верхних слоях и в бортовых примыканиях до Кф = 0.3 см/сек ( а = 0.95 ) в нижней части плотины. Некоторы результаты пространственного расчета фильтрации в теле плотины даны на рис.5. На рис.5,а показаны в плане следы поверхностей равных напоров и поверхностей тока на депрессионной поверхности и низовом откосе плотины. Линии равных напоров (эквипотенциали) проведены с шагом 0.1 в долях относительной величины статического напора. Следы поверхностей тока построены по касательным к проекциям вектора скорости на плоскость чертежа. Можно отметить четко выраженный пространственный характер движения фильтрационного потока с отклонением линий тока к правому борту каньона. К нижнему бьефу линии тока концентрируются к русловой части створа. Направление движения потока на депрессионной поверхности иллюстрируется на рис.5, б, где по продольным сечениям построены эпюры составляющей скорости Vx. Как видно из данного рисунка, со стороны верхнего бьефа поток в плане направлен в сторону выемки в правобережном борту, где максимальная составляющая скорости Vx достигает величины 0.48 см/сек. Затем, в сечении В-В на выходе из выемки, поток меняет направление на противоположное, а по мере приближения к нижнему бьефу вновь изменяет направление, повторяя изгибы правого борта каньона. В результате, на выходе фильтрационного потока в нижний бьеф ниже кривой высачивания депрессионной поверхности на низовой откос образуются два встречных потока, смыкающихся в районе продольного сечения VI-VI. Максимальная выходная скорость здесь равна 1.87 см/сек.

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

эквипотенциален на границах зон с различными коэффициентами фильтрации.

Для поперечного сечения А-А, параллельного плоскости XOZ, построены положение депрессионной поверхности, следы поверхностей тока и направления проекций векторов скорости фильтрации на плоскость XOZ (рис. 5, г). Можно отметить достаточно плавный характер следов поверхностей тока в верховом сечении А-А, отклонение фильтрационного потока в сторону правобережной выемки. По мере продвижения потока в сторону нижнего бьефа и, особенно, в низовом сечении С-С, наблюдается интенсивная циркуляция воды в вертикальном направлении, которая образуется в результате взаимодействия встречных потоков.

Сформировавшаяся в плотине депрессионная поверхность имеет сложную пространственную форму. Как видно из рис. 5, г, положение депрессионной поверхности изменяется от борта к борту: ее средний уклон в направлении от верхнего к нижнему бьефу уменьшается от правого борта к левому борту, и в продольном сечении IX-IX она становиться горизонтальной. Одновременно наблюдается понижение уровня депрессионной поверхности от центра потока к бортам (сечения А-А, В-В, С-С). Высота выклинивания на низовой откос также непостоянна - в русловом сечении V-V она достигает максимальной отметки 21.0 м, уменьшаясь к бортам каньона.

Величина фильтрационного расхода для контроля определялась по двум сечениям - входному (вдоль верхового откоса) и выходному (вдоль низового откоса). Невязка по величинам расхода составила 5% при его средней величине 17 м3/сек. При этом получено неравномерное распределение расхода по высоте плотины. Основная часть потока фильтрует через две верхние зоны с максимальными коэффициентами фильтрации, а в нижних зонах значения скорости резко уменьшаются. В связи с этим была дана рекомендация о возможном снижении проницаемости верхних слоев за счет их замыва мелкозернистым грунтом или формирования более высокой плотины, что позволит исключить высокопроницаемый материал из фильтрационной области. При условии снижения проницаемости верхних слоев до уровня нижних слоев следует ожидать снижение фильтрационного расхода в 2.5 - 3.0 раза.

Нелинейность закона фильтрации значительно сказывается на проницаемости материала в теле плотины, которая в зонах, близких к верхнему бьефу, оказывается в 30-50 раз выше, чем со стороны нижнего бьефа. Это приводит к «запиранию» фильтрационного потока со стороны нижнего бьефа и подъему депрессионной поверхности по сравнению с режимом линейной фильтрации, где величина коэффициента фильтрации не зависит от градиента.

Для сравнения результатов этим же методом была решена плоская задача фильтрации через русловое сечение плотины. Определена величина удельного фильтрационного расхода - 0.085 м3/сек. При пересчете на объемную задачу в самом узком сечении (В-В) получен полный расход порядка 15 м3/сек. При этом в плоской задаче точка высачивания получена на 14 метров ниже, чем в русловом сечении V-V пространственной задачи (рис.

5, в). Это объясняется сильным сжатием фильтрационного потока в узкой части каньона.

Получены также другие отличия результатов пространственной и плоской задач. Так, в плоской задаче получено более равномерное распределение компоненты скорости Уу по глубине потока, что объясняется перераспределением плановых скоростей в пространственной задаче. Максимальные скорости на выходе фильтрационного потока в нижний бьеф в плоской задаче оказались в 1.2 раза ниже, чем в пространственной.

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

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

геофильтрационной модели ЦСГНЭО ОАО "Институт Гидропроект" каждая "тектоническая зона" представляет собой трехслойное образование: в центре располагается так называемая «зона сместителя» (зона дробления), заполненная материалами с определенной водопроницаемостью. По бокам от «зоны сместителя» располагаются «зоны влияния», представляющие собой области повышенной трещиноватости и водопроницаемости. В конечно-элементной аппроксимации в плоской и в пространственной постановках «зоны сместителей» моделировались трещинными элементами, а «зоны влияния» - обычными конечными элементами с толщинами, равными толщинам зон.

Фильтрационная задача в плоской постановке решалась для 3-х сечений по осям секций №№ 16, 21 и 29 бетонной плотины. Учет «тектонических зон» изменил распределение линий равных напоров и увеличил удельные расходы ~на 10-И9 % под секциями № 16,21 и ~ на 50% в основании секции № 29. Общий фильтрационный расход с учетом "тектонических зон" в русловой части, полученный из решения плоских задач, ориентировочно равен примерно 0,16 м3/с.

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

Плотння Бурейской ГЭС

| »»о

Рис,б. Пространствашая фильтрация в основании плотины Б;урейжой ГЭС. а) сетка МКЭв сечениит оси плотины; б) рез5шьтаты реикния в лею бережной чети створа пшгины

1 — 'Ьлекы" Экггат7генцкапей;2 - '£леды"депргсснокн:й тгиерхнссш; 3 — га&зеет граница 1-господ осюЕанпя (х =23 м/сут); 4 —нижняягражад 3-рэ агая сскса&шл (к-=2 м/суг); 5 —гектойк^спе зоны; 6 - хекешздюнная - дренаа: основания

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

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

Без учета «тектонических зон» величина фильтрационного расхода через расчетную область составляет величину ~0,313 м3/сек, большая часть которого приходится на основание под русловой частью створа (-0,160 м3/сек). Учет "тектонических зон" с большой водопроницаемостью значительно увеличил значения расходов под русловой часть и под правобережной глухой частью плотины, где сосредоточены эти нарушения. Расход под русловой частью плотины увеличился с 0,160 м3/сек до 0,709 м3/сек, под правобережной глухой частью плотины с 0,059 до 0,501 м3/сек. Величина суммарного расхода в створе увеличилась с 0,313 до 1,365 м3/сек. Полученная величина фильтрационного расхода С)ф и 1,4 м3/сек при учете "тектонических нарушений" с максимальной водопроницаемостью, по видимому, является близкой к верхнему пределу, так как в расчеты вводились максимальные значения коэффициентов фильтрации и толщины зон с повышенной водопроницаемостью. Эта величина может быть использована при выборе контрольного параметра при натурных измерениях фильтрационного расхода. Проницаемость «тектонических зон» значительно влияет на величину фильтрационного расхода. Так, при уменьшении величин коэффициентов фильтрации «зон влияния» в 10 раз по сравнению с максимальными принятыми, расход снижается до ~ 0,607 м3/с (из них 0,215 м3/с - под русловой частью).

В целом по проведенным исследованиям фильтрационного режима основания плотины Бурейской ГЭС можно сделать следующие выводы. Решение пространственной фильтрационной задачи позволило получить картину распределения фильтрационного потока с учетом расположения «тектонических зон» в объеме расчетной области. Основная часть фильтрационного расхода проходит в основании русловой части и правобережной глухой плотины, где и сосредоточены крупные трещины. Приведенное исследование позволило рассмотреть в новой постановке движение фильтрационного потока в скальном основании. Движение воды рассмотрено с использованием закона Дарси для цельных пористых блоков и гидравлических решений Г.М.Ломизе при движении потока по крупным трещинам, которые в то же время могут играть роль дренажа. Такой прием привел к увеличению фильтрационного расхода. Учет трещиноватости повысил фильтрационный расход до ~1,4 м3/с.

Нестационарный фильтрационный режим грунтовых плотин, создающих бассейны приливных электростанций (ПЭС), или так называемых «отсечных» дамб вызывает особое внимание при проектировании. Особенности такой конструкции в фильтрационном смысле заключаются в следующем. Дамба работает в двустороннем нестационарном режиме: уровень воды больше то со стороны моря, то со стороны бассейна. При этом период колебания достаточно короткий (для Северной ПЭС он составляет 12 часов). Водонепроницаемость дамбы должна обеспечивать фильтрационный расход через ее створ несоизмеримо меньший по сравнению с полезным объемом бассейна. В противном случае возможно снижение напора на агрегаты ПЭС и, как следствие, снижение выработки электроэнергии. В качестве основного материала для ограждающих бассейн дамб чаще всего

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

Исследования фильтрационного режима дамбы Северной ПЭС в губе Долгая проводились с целью оптимизации ее возможной конструкции. Дамба имеет максимальную высоту -42,55 метра (отметка гребня 13,55, нижняя отметка поверхности основания со стороны бассейна -29,0, со стороны моря -21,0). В основании залегают слои различной мощности и проницаемости.

Графики изменения во времени уровней моря и бассейна в случае малофильтрующей дамбы (когда фильтрационные потери через ограждающую дамбу незначительны и практически не сказываются на уровне воды в бассейне ПЭС) приведены на рис.7, а. Рассматривалось четыре варианта конструкции грунтовой ограждающей дамбы: однородная дамба из горной массы, каменно-набросная дамба с толстым трапецеидальным ядром толщиной 10 м по верху и 60 м по основанию, каменно-набросная дамба с прямоугольным ядром толщиной 20 м и каменно-набросная дамба с трапецеидальным ядром толщиной по верху 5 м и по низу 20 м.

Решалась нестационарная фильтрационная задача с колебаниями уровней со стороны моря и бассейна с заданным графиком изменения уровней (см. рис.7). Временной интервал разбивался на шаги по времени, равные 1 часу. При решении фильтрационной задачи учитывалось, что в горной массе имеет место турбулентная фильтрация, и использовалась известная формула (6) зависимости скорости фильтрации и фильтрационного градиента. Для гравийно-галечниковых и песчано-гравийных грунтов рассматривалась линейный закон фильтрации, и использовался, соответственно, закон Дарси.

Для первого варианта дамбы при коэффициенте фильтрации горной массы Кф= 65000 м/сут или 75 см/с фильтрационные потери настолько велики, что пришлось параллельно с решением фильтрационной задачи также определять ориентировочные положения уровней воды в бассейне. Для этого на каждом временном шаге определялся удельный фильтрационный расход. В том случае, если движение фильтрационного потока было направлено из бассейна в море, оценивался полный фильтрационный расход в течение часа через всю дамбу (примерная длина эквивалентной дамбы с ее высотой 42,55 м в рассматриваемом сечении была принята 500 м). Полученные фильтрационные потери вычитались из объема бассейна. По батиграфической характеристике бассейна уточнялся уровень воды в бассейне, закладываемый в расчет на следующем временном шаге. Аналогично корректировался уровень бассейна и в случаях, когда фильтрационный расход направлен из моря в бассейн. В результате был получен график изменения уровней в бассейне с учетом фильтрационных потерь, представленный на рис.7, б.

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

а) Уровень еоды, м

- - * -/

■ч у ч ■у

ч * г у \ Ч ч и у Х-

О 2 4 6 8 10 12 14 16 18 X 22 2»

Времясга1вла суток, "не

о) Уровень юдл.м

суток; щс

в) Уровень во до, м

20 22 24 Врешсш'вла суток, *вс

Рис. 7. Графики кз лебания уровней моря и бассейна и фильтрационных напорсе в ядое

даьйы во вр вини (вт иенш суток), а-дтшлэфтгрунвдйсграждаощй даи6ы(сяд)омвлиезамывом горней массы песком).

6 - дня оувя о дноро дзой огрвждаюцрй дамбы из горной мас(ы (= 6500 0 м' суг); в - дгншшфшшгрующей ограждающей дамбы страпецодальнымядрсм

йлзкчые. обозте>ш\- уровшь моря;

_ —,, давни бассйна;

......фильтрационные напоры в ядое дамбы (по оси я фа соответственно

на отметках 2 3,0; -1 ОД -24,5)

перепад между уровнями воды в бассейне и море. Если для варианта с малофильтрующей дамбой (рис. 7, а) этот перепад достигает почти 2 метров момент времени 9-10 часов с начала суток), то для рассмотренного варианта ограждающей дамбы этот перепад снижается до 0,8-1,0 метра (момент времени 14 часов на рис.7, б). Остальные рассмотренные варианты с

противофильтрационными ядрами различной конструкции также допустимы с точки зрения фильтрационных расходов и могут рекомендоваться к возведению при экономическом обосновании. Вариант однородной дамбы из (горной массы возможен только за счет увеличения мелких фракций в материале наброски с <¿<5 мм.

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

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

Рассматривалось послойное возведение бетонного столба на массиве основания. Расчеты проведены для трех случаев температурного воздействия среды: температура воздуха принималась постоянной и равной 30°С, 20°С (летний период) и 5 °С (зимний период). Температура укладываемой бетонной смеси принималась постоянной и равной 1б с= 16 °С.

Для анализа влияния факторов и создания имитационной модели использовалась методика факторного анализа. Рассматривался полно факторный эксперимент, для которого функция откликов принимает следующий вид:

= Ь0 + Ь,Х, + Ь2Х2 + Ь3Х3 + Ь4Х4 + Ь,2 X, Х2 + Ь13 X, Х3 + Ь,4 X, Х4 + Ь23 Х2 Х3 + Ь24 Х2 X) + Ь34 Х3 Х4 + Ь[234 X, Х2 Х3 Х4 (13)

В качестве факторов в проведенных исследованиях рассматривались следующие величины. X] - расход цемента (принято, что Х]= -1 - нижний уровень фактора, соответствующий минимальному значению расход цемента - 50 кг/м3, Х]= +1 - верхний уровень фактора, соответствующий максимальному значению расход цемента 170 кг/м3); Х2 - толщина укладываемого слоя бетона ( в первом интервале от 0.3м до 1,0 м, во втором интервале от 1,0 до 5,0 м); Х3 - интенсивность (скорость) бетонирования массива по высоте (рассмотрено два интервала изменения этого фактора:

первый интервал - от 0,15 до 0,6 м/сут, второй интервал - от 0,6 до 2,0 м/сут); Х4 - полное тепловыделение цемента (от 120 КДж/кг до 350 КДж/кг). В качестве откликов рассматривались следующие величины: максимальная температура tmax, возникающая в возводимом массиве или максимальное приращение температуры Atmax по отношению к начальной температуре бетонной смеси (разница между максимальной температурой и температурой укладки бетонной смеси). Были получены функции откликов в виде полиномов для определения максимальных температур внутри возводимого бетонного массива при интенсивности бетонирования в пределах от 0,15 до 0,6 м/сут и толщинах слоев от 0,3 до 1,0 метра (после исключения малозначимых членов):

- при температуре наружного воздуха 20°С

W< = 24,65+3,21Хр 0,02Х2+1,77Хз+2,55Х4+ 1,43X1X3+1,88X1X4 +

0,32Х2Х3+0,9ЗХ3Х4+0,74Х,Х3Х4; (14)

- при температуре наружного воздуха 5°С

tmax=16,15+3,32X1+1,01X2+3,08X3+3, IOX4+1,14Х,Х3+1,91Х,Х4+

0,34Х2Х3+0,94Х3Х4+0,43Х,Х3Х4. (15)

Полученные зависимости удовлетворяют условию адекватности.

В рассмотренных примерах интенсивность бетонирования задавалась в интервале от 0,15 до 0,6 м/сут, что, в основном, и имеет место при современных методах возведения плотин из укатанного бетона. Однако, в некоторых случаях возможно возведение и при более высоких темпах возведения. Например, при бетонировании на высоких отметках, когда объемы бетонирования снижаются ввиду уменьшения ширины бетонируемого профиля, а производительность бетонного хозяйства достаточно высока. Кроме того, постоянное совершенствование технологии не исключает в недалеком будущем повышения интенсивности бетонных работ. Поэтому исследования были расширены и рассматривались примеры с бетонированием бетонного массива при интенсивности бетонирования в интервале от 0,6 до 2,0 м/сут и толщинах слоев от 0,3 до 1,0 метра. Были получены функции откликов аналогичные выражениям (14, 15)

Для получения математической модели при возведении сооружения из вибрированного бетона были проведены исследования также с применением факторного анализа. Рассматривались те же факторы, что и в предыдущих исследованиях, но с другими интервалами изменения: X] - расход цемента от 170 до 400 кг/м3; Х2 - толщина укладываемого слоя бетона от 1,0 до 5,0 м; Хз - интенсивность возведения массива по высоте для первого факторного пространства от 0,15 до 0,6 м/сут, для второго факторного пространства от 0,6 до 2,0 м/сут; полное тепловыделение цемента от 120 до 350 КДж/кг. В качестве откликов также рассматривались максимальная температура или максимальный разогрев внутри бетонного массива. Были получены зависимости, аналогичные уравнениям (14-15).

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

Номограмма для определения максимальной температуры бетона после укладки (/в0=5°С, расход цемента 50-170 кг/м3, толщина слоя 0.3-1.0 м, Этах 120-350 КДж/кг, интенсивность бетонирования

0.6-2.0 м/сут)

Номограмма для определения максимальной температуры бетона после укладки 0„о -20°С, расход цемента 50-170 кг/м3, толщина слоя 0.3-1.0 м, Этах 120 - 350 КДж/кг, интенсивность бетонирования

0.6-2.0 м/сут)

Толщина (м);

Рис. 8. Номограммы для определения максимальной температуры в бетонном массиве.

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

Полученная имитационная модель сравнивалась с аналогичной моделью возведения бетонного столба плотины Бурейской ГЭС, полученной во ВНИИГе им Б.Е.Веденеева. Значения температурного разогрева вследствии экзотермии, полученные по двум разным математическим моделям, практически равны, что подтверждает достоверность полученной модели.

По разработанной методике были исследованы температурные режимы гравитационных плотин из укатанного бетона: Хлонг Та Дат (Тайвань) и Шон Ла (Вьетнам). Получены температурные поля для различных моментов времени периодов строительства и эксплуатации, необходимые для дальнейшего решения задачи термонапряженного состояния.

Исследования термического режима грунтовой плотины Богучанской ГЭС проводились для определения температурного состояния асфальтобетонной диафрагмы. Прогноз температурного режима асфальтобетонной диафрагмы необходим как для определения временных шагов технологического процесса укладки асфальтобетона и грунта, так и для определения свойств асфальтобетона, влияющих на напряженно-деформированное состояние системы «плотина-диафрагма»

Рассматривался процесс остывания асфальто-бетонной диафрагмы при ее послойной укладке при различных температурах воздуха: зимний вариант -24,1°С° и летний вариант +20,0°С. Возведение диафрагмы производится слоями по 0,6 метра, каждый из которых, в свою очередь, укладывается в триприема по 0,2 метра с интервалом укладки ~ 6 часов. Толщина асфальтобетонной диафрагмы -1,2 метра. Укладка и остывание асфальтобетона от температуры укладки 160°С до 40°С происходит в теплоизолирующей опалубке (фенолпласт толщиной 6 см) со стороны боковых поверхностей. В резульате численного решения получены минимально допустимое время охлаждения асфальтобетона: в зимнее время -48 часов и в летнее время ~72 часа, что хорошо корреспондируется с результатами физического эксперимента опытного фрагмента в г. Братске (ВНИИГ, 1984 г.). Для прогноза температурного режима конструкции в строительный и эксплуатационный периоды решалась 2-х мерная нестационарная задача теории теплопроводности.

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

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

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

Глава 5 посвящена совместной фильтрационно-температурной задаче, решение которой необходимо для сооружениий, возводимых в условиях вечной мерзлоты. Вопросы обеспечения устойчивости и надежности таких сооружений невозможно решать без четкого представления о фильтрационном и температурном режимах плотин и их оснований. Исследования, связанные с совместным решением фильтрационно-температурной задачи, были проведены на примере плотины Курейской ГЭС. Рассматривался участок грунтовой плотины во втором правобережном понижении. Данный участок плотины представляет собой интерес в связи с поведением этой конструкции в период эксплуатации. В 1989 году в пределах данного участка (у ПК 10) произошла аварийная ситуация, связанная с выходом фильтрационного потока на низовой откос грунтовой плотины. Осадки гребня плотины за 10-летний период эксплуатации (с 1990 по 2000 гг.) достигли 64 см и продолжали постепенно увеличиваться. Дальнейшее продолжение роста осадок могло вызвать понижение гребня противофильтрационного элемента плотины до отметки ниже уровня водохранилища и необходимость его снижения.

Расчеты температурного режима грунтовой плотины проведены в соответствии со схемой поэтапного возведения и наполнения водохранилища. К моменту возведения плотины до ее полного профиля (конец 1987 года), все ее поперечное сечение за исключением ядра выше отметки ~ V90,0 м находилось в промороженном состоянии. К моменту наполнения водохранилища это состояние сохранялось. В средней части профиля отрицательные температуры составляли величины - 0,4-ь -0,6 °С.

На начальных этапах наполнения водохранилища воздействие фильтрационного потока не ощущалось. В дальнейшем тепломассоперенос за счет фильтрационного потока начинает сказываться интенсивнее и к декабрю 1989 года тепловой поток проникает через каменную наброску и температура на верховой грани экрана достигала положительных температур ~1-г2°С.

Температурное поле грунтовой плотины в сечении ПК 12+50 в 1993 году (рис.10) характеризовалось следующим. В талом состоянии находится довольно значительная часть экрана. По подошве экрана положительная температура достигла низовой грани экрана. Низовая призма плотины - в мерзлом состоянии. В средней ее части температура в пределах -0,2 4- -1,0°С. Продолжается тепломассоперенос на границе экрана и цокольной части плотины. Показательна динамика изменения температуры в точке

пересечения низовой грани экрана с поверхностью цокольной части (точка А, рис. 10): - 0,092 °С в январе, - 0,008 °С в июле и + 0,05 °С в сентябре месяце.

К 2002 году процесс проникновения тепла фильтрационным потоком продолжался. В январе 2002 года по сравнению с январем 1993 года граница оттаявшего грунта сместилась еще дальше в сторону нижнего бьефа. По контакту плотины с цокольной частью оттаял грунт экрана и незначительная часть грунта низовой призмы. В течение 2002 года наблюдалось увеличение положительных температур в области фильтрации (область талых грунтов). Для сравнения с 1993 годом в точке А значения температуры в течение 2002 года изменяются следующим образом: январь - +0,74°С, апрель - +0,76°С, июль - +0,78°С, сентябрь - +0,81°С.

К 2015 году тенденция проникновения теплового потока по контакту экрана и цокольной части плотины сохранялась, хотя скорость этого процесса достаточно низка. Значения температуры в точке А в 2015 году следующие: в январе +1,14°С, в апреле +1,15°С, в июле+1,16°С и в сентябре +1,18°С. Расстояние от точки А до границы талого и мерзлого грунтов по горизонтали в сторону нижнего бьефа в сентябре 2015 года составляет примерно 17,0 м. Для сравнения в сентябре 1993 года это расстояние равнялось ~ 0, а в сентябре 2002 года - 12,3 м. Таким образом, средняя горизонтальная скорость продвижения нулевой изотермы в сторону нижнего бьефа составляет ~1,4 м/год за период 1993 н- 2002 годы и ~0,36 м/год за период 2002ч-2015 годы. Если судить по данным расчёта, то будет иметь место некоторая стабилизация температурного режима плотины к 2015 году.

Расчеты в пространственной постановке подтвердили качественную картину фильтрационно-температурного режима, полученную из решения плоской задачи. Несколько меньшее проникновение нулевой изотермы внутрь конструкции получено у бортов плотины. Были примерно оценены скорости продвижения нулевой изотермы в нижней части экрана и цоколя. Продвижение границы мерзлых и талых грунтов за период 19934-2002 годов характеризовалось средней скоростью по горизонтали ~ 0,9 м/год. В вертикальном направлении средняя скорость продвижения нулевой изотермы составляла примерно 0,5 м/год. За период 20024-2015 гг. эти величины соответственно равны ~ 0,3 м/год по горизонтали и ~0,1 м/год в вертикальном направлении. Таким образом, скорость продвижения тепловой волны в пространственной задаче меньше, чем в плоской. Сохраняется тенденция ее уменьшения со временем.

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

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

Общие выводы.

1. На основе метода конечных элементов в локально-вариационной постановке разработана методика расчетов фильтрационных, температурных и совместных фильтрационно-температурных задач в плоской и пространственной постановках, позволяющая решать задачи с достаточно полным учетом влияющих на исследуемые процессы факторов. Разработанная методика реализована в алгоритмах и программных комплесах ШЛИ и ТЕКМЮ. Достоверность получаемых с использованием методики и программ расчета подтверждена решением ряда тестовых задач и сравнением результатов численных исследований с апробированными методами и натурными данными.

2. Были проведены сопоставления результатов плоских и пространственных задач в различных условиях: при линейных и нелинейных законах фильтрации, в стационарной и нестационарной постановках, для нескальных оснований и скальных грунтов с трещинами, при различной густоте сеток МКЭ, при двусторонних и односторонних колебаниях уровня воды, при одновременном действии фильтрационного и температурного потоков. Эти сопоставления показали, что практически во всех случаях надо стремиться к решению пространственных задач. Даже в относительно «простом» створе плотины Юмагузинского гидроузла с коэффициентом створа порядка 10 фильтрационный поток имеет ярко выраженный пространственный характер.

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

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

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

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

4. Решения нестационарных и нелинейных фильтрационных задач всегда вызывали достаточно большие сложности. Сейчас получена возможность решения таких задач в плоской и пространственной постановках. Решение задач применительно к взрывонабросным плотинам снова приобретет большое значение, когда мы вплотную приблизимся к освоению таких створов как Камбаратинский №1(р. Нурек) и Достимжумский (р. Пяндж). К этому уже надо быть готовыми.

5. Фильтрационно-температурный класс задач очень важен для освоения Севера. Решения этих задач диктовались проблемами строительства в условиях крайнего Севера. И Курейская ГЭС тому пример. Успешное решение задач тепломассопереноса позволило прогнозировать движение нулевой изотермы в плотине в пространстве и времени. Полученное в результате расчетов распределение температуры в плотине достаточно приемлемо корреспондируется с данными натурных наблюдений.

6. XX век был веком решения фильтрационных и температурных задач аналитическими методами, в основном, плоских, в стационарной и линейной постановках. Во второй половине века появились численные методы, но возможности вычислителной техники были еще весьма ограничены. К концу XX столетия, с появлением мощных быстродействующих ЭВМ происходит «бум» численных решений. Это был, в определенном смысле, прорыв. XXI век будет, вероятно, веком решения нестационарных нелинейных задач в пространственной постановке. Современная вычислительная техника уже сейчас позволяет решать эти задачи, но она быстро совершенствуется, и ее возможности будут необъятны.

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

Основные положения диссертации изложены в следующих работах:

1. Анискин H.A. Расчет температурного режима и напряженно-деформированного состояния арочной плотины в «объемной» постановке. Труды Всесоюзного координациооного совещания «Прочность и температурная трещиностойкость бетонных гидротехнических сооружений при температурных воздействиях, -Л.,Энергоиздат, 1989.

2. Анискин H.A., То Ван Тхань. Прогноз фильтрационного режима грунтовой плотины Юмагузинского гидроузла и ее основания. «Гидротехническое строительство», 2005, № б, М., Энергопрогресс, с. 19-25.

3. Анискин H.A. Температурный режим гравитационной плотины из укатанного бетона. «Гидротехническое строительство», 2005, № 12, М., Энергопрогресс, с. 13-17.

4. Анискин H.A. Фильтрация в основании и бортах бетонной гравитационной плотины Бурейского гидроузла. «Вестник МГСУ», 2006, №2, с. 87-96.

5. Анискин H.A. Температурно-фильтрационный режим основания и плотины Курейской ГЭС во втором правобережном понижении. «Вестник МГСУ», 2006, №2,43-52.

6. Анискин H.A. Исследования фильтрационного и температурно-фильтрационного режимов плотин и оснований. Материалы Второй научно-технической конференции «Гидроэнергетика. Новые разработки и технологии», ВНИИГ им. Б.Е.Веденеева, г. Санкт-Петербург, 2006 г.

7. Анискин Н.А.,Нгуен Данг Жанг. Численное моделирование температурного режима гравитационной плотины из укатанного бетона. Строительные материалы, оборудование, технологии XXI века, 2006, №10, МГСУ, Москва, с. 30-32.

8. Анискин H.A. Численное моделирование фильтрации в трещиноватых скальных основаниях плотин. «Инженерные изыскания», декабрь 2007 года.

9. Анискин H.A. Неустановившаяся фильтрация в грунтовых плотинах и основаниях. Сборник «Вестник МГСУ», 2009, №2, с.70-79.

10. Анискин H.A., Нгуен Данг Жанг Прогноз температурного режима бетонных гравитационных плотин из укатанного бетона. «Гидротехническое строительство», 2007, № 12, М., Энергопрогресс, с. 8-14.

11. Рассказов Л.Н., Анискин H.A. и др. Фильтрация в грунтовых плотинах в плоской и пространственной постановке. «Гидротехническое строительство», 1989, №11.

12. Рассказов JI.H., Анискин H.A. Фильтрационные расчеты гидросооружений и оснований. «Гидротехническое строительство», 2000, №11, с. 2-7.

13. Рассказов JI.H., Анискин H.A. и др. Сангтудинский гидроузел: напряженно-деформированное состояние и фильтрация в основании плотины и в обход гидроузла. «Гидротехническое строительство», 2008, №5.

14. Рассказов Л.Н., Орехов В.Г., Анискин H.A. и др. Гидротехнические сооружения. Учебник для ВУЗов, АСВ, М., 2008.

15. Рассказов JI.H., Анискин H.A. и др. Конструктивные решения морских грунтовых плотин при строительстве ПЭС, «Гидротехническое строительство», 2009, №7.

КОПИ-ЦЕНТР св. 7:07:10429 Тираж 100 экз. г. Москва, ул. Енисейская, д.36 тел.: 8-499-185-7954, 8-906-787-7086

Введение.г.

Глава 1. Развитие теории и методов решения задач фильтрации и теплопроводности в гидротехнике.

1.1. Развитие теории фильтрации.

1.2. Решение фильтрационных задач в области гидротехники и гидромелиорации.

1.3. Методы решения температурных задач.

1.4. Выводы по главе 1.

Глава 2. Численное решение фильтрационных и температурных задач на основе метода конечных элементов.

2.1. Математические основы фильтрационно-температурных задач.

2.2. Основные теоретические положения решения фильтрационных задач.

2.3. Теоретические основы решения задач теории теплопроводности.

2.4. Совместная температурно-фильтрационная задача.

2.5. Метод конечных элементов в локально-вариационной постановке.

2.6. Выводы по главе 2.

Глава 3. Фильтрационный режим плотин и оснований.

3.1. Фильтрационный режим системы «плотина-основание».

3.1.1. Плоская и пространственная фильтрационные задачи на примере плотины Юмагузинского гидроузла.

3.1.2. Математическая модель фильтрационного режима системы «плотина-основание» на примере Сангтудинского гидроузла.

3.1.3. Фильтрация склона напорных водоводов и ограждающей дамбы Загорской ГАЭС-2.

3.2. Фильтрационные задачи в нелинейной постановке на примере взрывонабросной плотины Камбаратинского гидроузла).

3.3. Фильтрация в трещиноватых скальных основаниях на примере плотины Бурейской ГЭС.

3.4. Неустановившаяся фильтрация в грунтовых плотинах и основаниях.

3.4.1. Постановка задачи и обзор развития методов решения нестационарных фильтрационных задач.

3.4.2. Фильтрационный режим земляной плотины.

3.4.3. Решение нестационарной фильтрационной задачи на примере отсечной дамбы Северной ПЭС.

3.5. Выводы по главе 3.

Глава 4. Температурный режим гидросооружений.

4.1. Температурный режим бетонных плотин.

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

4.1.2. Температурный режим гравитационных плотин из укатанного бетона.

4.2. Температурный режим каменнонабросной плотины с асфальтобетонной диафрагмой (на примере плотины Богучанского гидроузла).

4.3. Выводы по главе 4.

Глава 5. Фильтрационно-температурный режим плотин и их оснований.

5.1. Исследования фильтрационно-температурного режима грунтовых плотин и оснований.

5.2. Температурно-фильтрационный режим грунтовой плотины и ее основания на примере плотины Курейской ГЭС.

5.3. Выводы по главе 5.

Представленная работа выполнена на кафедре Гидротехничеких сооружений Московского государственного строительного университета. Разработки и исследования, приведенные в данной работе, в основном, были выполнены в период с 1999 по 2009 год.

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

Большую признательность автор выражает также профессорам кафедры Орехову Валентину Гавриловичу, Малаханову Вячеславу Васильевичу, доцентам Толстикову Виктору Васильевичу, Бестужевой Александре Станиславовне, Саинову Михаилу Петровичу за консультации и помощь в работе. Автор также благодарен секретарям кафедры Овчинниковой Ирине Леонидовне и Сусариной Татьяне Валентиновне за помощь в оформлении диссертационной работы.

ВВЕДЕНИЕ

Настоящая работа выполнялась в период с 1990 г. по 2009 г. и является обобщением научных исследований, выполненных автором на кафедре Гидротехнических сооружений Московского государственного строительного университета.

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

Гидротехнические сооружения и, в частности, плотины, работают в чрезвычайно сложных условиях, подвергаясь в процессе строительства и эксплуатации действию множества нагрузок и воздействий. Одним из основных воздействий является фильтрация как непосредственно в плотинах (особенно в грунтовых), так и в основаниях и их бортовых примыканиях. Помимо фильтрационных потерь из водохранилища, это явление вызывает дополнительные фильтрационные нагрузки и может привести к возникновению фильтрационных деформаций. Чрезвычайно важной задачей при проектировании плотины является прогноз фильтрационного режима системы «плотина-основание». Как правило, каждая плотина является по-своему уникальным сооружением, включающим в свой состав различные конструктивные элементы (например, противофильтрационные экраны, ядра из глинистых грунтов, противофильтрационные элементы из негрунтовых материалов - экраны и диафрагмы из бетона, железобетона, асфальтобетона, дренажи, переходные зоны и т.д.). Часто плотина возводится в сложных инженерно-геологических условиях, и ее основания и борта представляют собой весьма разнородную по фильтрационным свойствам структуру, часто с явно выраженной анизотропией свойств, с присутсвием в ней трещин, разломов , зон повышенной проницаемости и т.д. Топография створа строительтва, посадка в нем сооружения, система противофильтрационных мероприятий в основании и бортах сооружения часто вызывают сложный пространственный характер фильтрационного потока. Учет множества факторов, влияющих на фильтрационный режим системы «плотина-основание» при решении фильтрационной задачи является весьма сложной проблемой.

О важности правильного учета воздействия фильтрационного потока на сооружение и его основание свидетельствуют многочисленные исследования причин аварий грунтовых плотин. По опубликованным данным /40/, в качестве основных причин аварий грутовых плотин США в первой половине XX века называются: фильтрационный прорыв через тело плотины или через основание, фильтрационные деформации грунтов плотины и основания, оползни в зоне усиленной фильтрации. По оценке Дж. Шерарда (США)/26б/, сделанной в 1963 году на основе анализа 214 случаев аварий грунтовых плотин, 41% из них связан с фильтрацией. По другим оценкам /27/ отмечается, что около 60% повреждений и разрушений произошло, в конечном счете, в результате неблагоприятного фильтрационного воздействия. Из них примерно 45% связаны с фильтрацией непосредственно и произошли из-за фильтрационных деформаций (прежде всего механической и химической суффозий, фильтрационного выпора). В остальных 15% фильтрация является косвенной причиной разрушений, проявившихся в обрушении откосов, разжижении грунта и т.д. В качестве примера можно привести разрушение высокой плотины Титон (96 м, 1976, США), произошедшее вследствие нарушения фильтрационной прочности тела плотны у основания на контакте с правобережным примыканием /40/. Негативное воздействие фильтрации на грунтовую плотину оказывается как при наполнении водохранилища и первые годы эксплуатации (примерно около половины аварий произошло именно в этот период /40/), так и после многолетней эксплуатации (примерно пятая часть разрушений произошла после 50 лет эксплуатации). Известны и многочисленные аварии бетонных плотин и их оснований, вызванные фильтрационным воздействием: например, авария высокой арочной плотины Мальпассе (66,5 м, 1959 г., Франция).

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

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

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

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

Отепляющее воздействие фильтрационного потока может вызвать нежелательное чрезмерное оттаивание грунтов, что в свою очередь вызовет увеличение проницаемости плотины и основания и утечку воды из водохранилища, появление зон локальной усиленной фильтрации, увеличение пористости и осадок сооружения. Примером такого проявления может служить ситуация на плотине Курейской ГЭС (см. главу 5), хотя ее первопричиной являются скорее всего технологические факторы. Фильтрационный и температурный режимы в этом сложном нестационарном процессе взаимно влияют друг на друга. Очевидно, что прогноз фильтрационно-температурного режима должен постоянно вписываться в процесс проектирования, так как при принятии проектных решений необходимо учитывать результаты таких расчетов и периодически вносить коррективы в проект.

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

Целью диссертационной работы явились:

• разработка методики решения фильтрационных, температурных и температурно-фильтрационных задач применительно к системе «плотина-основание» с учетом пространственности, неоднородности и нелинейности свойств материалов плотины и основания, возможной анизотропии и нестационарности процесса;

• создание численных фильтрационных и температурно-фильтрационных моделей гидротехнических сооружений и их оснований;

• анализ влияния неонородности и нелинейности свойств материалов, пространственности, нестационарности процесса и других факторов на формирование фильтрационного и температурного режимов плотин и их оснований.

Основные задачи исследований:

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

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

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

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

Достоверность научных результатов подтверждена:

• результатами решения ряда тестовых задач, имеющих точное аналитическое решение;

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

• сравнением результатов численных исследований Юмагузинского, Сангтудинского, Курейского гидроузлов с данными натурных наблюдений.

Практическое значение работы и полученных результатов заключается: в разработке и внедрении метода расчета фильтрационных, температурных и совместных температурно-фильтрационных задач, позволившего создать пространственные математические модели системы «плотина-основание» с учетом природных и инженерно-геологических условий, конструктивных особенностей сооружений и технологии их возведения; в разработке комплексов вычислительных программ «БИ/ГЯ» и «ТЕИМТС» на базе метода конечных элементов для решения фильтрационных, температурных и совместных задач; в создании численных моделей фильтрационных режимов плотин совместно с основаниями применительно к следующим гидроузлам: Бурейскому, Камбаратинскому, Юмагузинскому, Сангтудинскому, Загорской ГАЭС-2, температурно-фильтрационной модели и прогнозу работы для плотины Курейской ГЭС, в решении задач по прогнозу нестационарных фильтрационных режимов ограждающих дамб Северной и Мезенской ПЭС; в использовании результатов численных исследований фильтрационных и температурных режимов гидроузлов в практике проектирования и строительства.

Апробация работы Основные положения методики и результаты исследований докладывались на научно-технических конференциях:

- на научно-технических конференциях МИСИ-МГСУ; на Всесоюзном координациооном совещании «Прочность и температурная трещиностойкость бетонных гидротехнических сооружений при температурных воздействиях, г.Нарва, 1989;

- на второй, третьей и четвертой научно-технических конференциях «Гидроэнергетика, новые разработки и технологии» в 2006, 2007 и 2008 годах в г. Санкт-Петербурге.

Личный вклад автора Диссертантом непосредственно разработана методика решения фильтрационных, температурных и совместных задач в плоской и пространственной постановках, методика реализована в программных вычислительных комплексах. Проведены численные исследования по Бурейскому, Камбаратинскому, Юмагузинскому, Сангтудинскому, Истринскому гидроузлам, Загорской ГАЭС-2, Курейской ГЭС, решены задачи по прогнозу нестационарных фильтрационных режимов ограждающих дамб Северной и Мезенской ПЭС. Сделан анализ влияния технологических факторов на температурный разогрев бетонного массива при его возведении. Проведены численные расчеты температурных режимов гравитационных плотин из укатанного бетона Хлонг Та Дат (Тайвань) и Шон Ла (Вьетнам) в строительный и эксплуатационный периоды, позволившие дать оценку эффективности принятых параметров возведения.

Публикации. Основные положения диссертации опубликованы в 16 научных работах, в том числе 9 работ в ведущих рецензируемых научных журналах и изданиях, рекомендованных ВАК: в журнале «Гидротехническое строительство» -6, сборнике «Вестник МГСУ» -2, в журнале «Строительные материалы, оборудование, технологии XXI века» -1. Основные положения методики и некоторые результаты исследований помещены в учебник «Гидротехнические сооружения» (издательство АСВ, 2008 год). На защиту выносятся: методика численного решения пространственных фильтрационных, температурных и температурно-фильтрационных задач с учетом множества действующих факторов; разработанные алгоритм и комплекс вычислительных программ по решению фильтрационных, температурных и совместных температурно-фильтрационных задач с учетом неоднородности, анизотропии и нелинейности свойств, пространственного характера и нестационарности процессов; ■ численные математические модели фильтрации и анализ влияния пространственности, нелинейности, нестационарности и других факторов на работу сооружений и противофильтрационных элементов, выполненные на примере Бурейского, Камбаратинского, Юмагузинского, Сангтудинского гидроулов, Загорской ГАЭС-2, температурно-фильтрационная модель и прогноз работы на перспективу на примере плотины Курейской ГЭС, результаты решения задач по прогнозу нестационарных фильтрационных режимов ограждающих дамб Северной и Мезенской ПЭС. Структура и объем диссертации Работа состоит из предисловия, введения, пяти глав, заключения и списка литературы. Объем диссертации составляет 416 страниц, из них 265 страниц текста, список литературы состоит из 272 наименований.

1. На основе метода конечных элементов в локально-вариационной постановке разработана методика расчетов фильтрационных, температурных и совместных фильтрационно-температурных задач в плоской и пространственной постановках, позволяющая решать задачи с достаточно полным учетом влияющих на исследуемые процессы факторов.Разработанная методика реализована в алгоритмах и программных комплесах FILTR и TERMIC. Достоверность получаемых с использованием методики и программ расчета подтверждена решением ряда тестовых задач и сравнением результатов численных исследований с апробированными методами и натурными данными.2. Были проведены сопоставления результатов плоских и пространственных задач в различных условиях: при линейных и нелинейных законах фильтрации, в стационарной и нестационарной постановках, для нескальных оснований и скальных грунтов с трещинами, при различной густоте сеток МКЭ, при двусторонних и односторонних колебаниях уровня воды, при одновременном действии фильтрационного и температурного потоков. Эти сопоставления показали, что практически во всех случаях надо стремиться к решению пространственных задач. Даже в относительно «простом» створе плотины Юмагузинского гидроузла с коэффициентом створа порядка 10 фильтрационный поток имеет ярко выраженный пространственный характер.3. Создание математической фильтрационной, температурной или совместной модели является в современных условиях процессом поэтапным, развивающимся параллельно с проектированием, строительством и мониторингом сооружения. На начальном этапе исследований дается общая оценка фильтрационного режима плотины и основания. Делаются выводы о необходимости проектирования противофильтрационных элементов в основании под плотиной и их размерах, оценивается фильтрационный расход. Сопоставляются результаты плоских и пространственных задач.На следующем этапе уточняются инженерно-геологические условия в основании и бортах створа, и корректируется пространственная модель.Даются рекомендации по совершенствованию и оптимизации противофильтрационных элементов. Возможно использование методики факторного анализа, что позволяет связать математическую модель с возможными комбинациями характеристик грунтов, входящих в расчетную область.Создание математической модели, в конечном счете, приводит к методике непрерывных в контакте с проектировщиками уточняющих исследований.Такая откалиброванная математическая модель позволяет осуществить прогноз поведения конструкции, что необходимо для мониторинга и декларации безопасности сооружения.4. Решения нестационарных и нелинейных фильтрационных задач всегда вызывали достаточно большие сложности. Сейчас получена возможность решения таких задач в плоской и пространственной постановках. Решение задач применительно к взрывонабросным плотинам снова приобретет большое значение, когда мы вплотную приблизимся к освоению таких створов как Камбаратинский №1(р. Нурек) и Достимжумский (р. Пяндж). К этому уже надо быть готовыми.5. Фильтрационно-температурный класс задач очень важен для освоения Севера. Решения этих задач диктовались проблемами строительства в условиях крайнего Севера. И Курейская ГЭС тому пример. Успешное решение задач тепломассопереноса позволило прогнозировать движение нулевой изотермы в плотине в пространстве и времени. Полученное в результате расчетов распределение температуры в плотине достаточно приемлемо корреспондируется с данными натурных наблюдений.6. XX век был веком решения фильтрационных и температурных задач аналитическими методами, в основном, плоских, в стационарной и линейной постановках. Во второй половине века появились численные методы, но возможности вычислителной техники были еще весьма ограничены. К концу XX столетия, с появлением мощных быстродействующих ЭВМ происходит «бум» численных решений. Это был, в определенном смысле, прорыв. XXI век будет, вероятно, веком решения нестационарных нелинейных задач в пространственной постановке.Современная вычислительная техника уже сейчас позволяет решать эти задачи, но она быстро совершенствуется, и ее возможности будут необъятны.7. Большинство из приведенных в данной работе исследований базировались на достаточно скудных исходных данных (по крайней мере, на начальных этапах). При решении фильтрационных задач главное - иметь достаточно полную и точную исходную информацию, но получить ее современными методами, которые практически себя исчерпали, невозможно.Нужны новые, видимо, геофизические методы изысканий. И, в первую очередь, это относится к скальным трещиноватым основаниям плотин.Методы эти должны быть не только качественными, но и количественными, что особенно важно.

1. Адлер Ю. П., Маркова Е. В. И Грановский Ю. В. Планированиеэксперимента при поиске оптимальных условий. Изд-во "наука" 1976, 280 с.

2. Александровский С В . Расчет бетонных и железобетонных конструкцийна изменения температуры и влажности с учетом ползучести. -М.: Стройиздат, 1973, 432с.

3. Алексеева Л. Г., Кузнецова О. М., Фрид А. Программа для расчетатемпературного режима массивных бетонных гидротехнических сооружений. - Материалы симпозиума "ЭВМ-ГЭС 73". -Л.: Энергия, 1973, 7с.

4. Алексеева Л. Г., Фрид А. Опыт применения ЭВМ для исследованиянестационарного температурного поля бетонных гидротехнических сооружений. Труды координационных совещаний по гидротехнике. -Л.: Энергия, вып. 103. 1975, с. 76-79.

5. Аликин В.Г. Фильтрационный расчет взрывонабросных плотин напроницаемых основаниях при нелинейном законе движения. Тр. ВОДГЕО, 1980.

6. Анахаев К.Н. Выбор экрана для каменно-земляных плотин,«Мелиорация и водное хозяйство», 1991, №3.

7. Анахаев К.Н., Пространственная фильтрация в ядрах каменно-земляныхплотин, «Гидротехническое строительство», 1994, №7, М., Энергопрогресс.

8. Анахаев К.Н. Совершенствование конструкций, методов расчетногообоснования и проектирования противофильтрационных устройств грунтовых плотин: Дис. доктора техн. наук, М.,МГУП, 1997.

9. Анахаев К.Н., Амшоков Б.Х., Ищенко А.В. О фильтрационном расчетеземляных плотин с ядром, «Гидротехническое строительство», 2006, № 5, М., Энергопрогресс.

10. Анахаев К.Н., Ляхевич Р.А. Фильтрация в анизотропных грунтовыхплотинах, «Гидротехническое строительство», 2005, № 4, М., Энергопрогресс.

11. Анахаев К.Н., Ляхевич Р.А. К фильтрационному расчету анизотропныхземляных плотин с наслонным дренажем, «Гидротехническое строительство», 2006, № 9, М., Энергопрогресс.

12. Анискин Н.А., То Ван Тхань. Прогноз фильтрационного режимагрунтовой плотины Юмагузинского гидроузла и ее основания. «Гидротехническое строительство», 2005, № 6, М., Энергопрогресс.

13. Анискин Н.А. Температурный режим гравитационной плотины изукатанного бетона. «Гидротехническое строительство», 2005, № 12, М., Энергопрогресс.

14. Анискин Н.А. Фильтрация в основании и бортах бетоннойгравитационной плотины Бурейского гидроузла. «Вестник МГСУ», 2006, №2.

15. Анискин Н.А. Температурно-фильтрационный режим основания иплотины Курейской ГЭС во втором правобережном понижении. «Вестник МГСУ», 2006, №2.

16. Анискин Н.А. Исследования фильтрационного и температурнофильтрационного режимов плотин и оснований. Материалы Второй научно-технической конференции «Гидроэнергетика. Новые разработки и технологии», ВНИИГ им. Б.Е.Веденеева, г. Санкт-Петербург, 2006 г.

17. Анискин Н.А.,Нгуен Данг Жанг. Численное моделированиетемпературного режима гравитационной плотины из укатанного бетона. Строительные материалы, оборудование, технологии XXI века, 2006, №10, МГСУ, Москва.

18. Анискин Н.А. Численное моделирование фильтрации в трещиноватыхскальных основаниях плотин. «Инженерные изыскания», декабрь 2007 года.

19. Анискин Н.А. Неустановившаяся фильтрация в грунтовых плотинах иоснованиях. Сборник «Вестник МГСУ», 2009, №2.

20. Анискин Н.А., Нгуен Данг Жанг Прогноз температурного режимабетонных гравитационных плотин из укатанного бетона. «Гидротехническое строительство», 2007, № 12, М., Энергопрогресс.

21. Антоне Э.Р. Электроинтегратор для моделирования нелинейногоуравнения Буссинеска в электролитической ванне. Труды координационных совещаний по гидротехнике. Выпуск 93. ВНИИГ им. Б.Е.Веденеева, «Энергия», Л., 1974.

22. Аравин В.И. К вопросу о фильтрации в анизотропно-водопроницаемыхгрунтах. Тр. ЛИИ, 1937, № 9.

23. Аравин В.И. Расчет фильтрации в обход гидротехнических сооружений,Изв. НИИГ, 1940, т.27.

24. Аравин В.И., Мошкова М.А. Исследование на щелевом лотке влиянияформы водоупора при неустановившейся фильтрации. Изв. ВНИИГ, 1964, т.76.

25. Аравин В.Н., Нумеров Н. Теория движения жидкостей и газов впористой недеформируемой среде. М., Гостехиздат, 1953.

26. Аравин В.Н., Нумеров Н. Фильтрационные расчеты гидротехническихсооружений-М.,-Л., Госстройиздат, 1955

27. Аравин В.И., Носова О.Н. Натурные исследования фильтрации(теоретические основы). Л.: Энергия. 1969.

28. Аргал Э.С., Тужихин Г.Г. Противофильтрационные устройства восновании плотины Бурейского гидроузла (опыт проектирования и первые результаты производства работ), Гидротехническое строительство № 8, Энергопрогресс, М., 2001.

29. Арсенин В.Я., Математическая физика. Основные уравнения испециальные функции. Издательство «Наука», М., 1966.

30. Арсеньева А. Н., Вайнер М. М. Опыт использования электроинтеграторадля исследования температурного режима бетонных и земляных гидротехнических сооружений - Труды координационных совещаний по гидротехнике. -Л.: Энергия, вып. 103. 1975, с. 79-81.

31. Бабаян А.Г. Конечно-элементная методика для совместных расчетовфильтрационного режима и статической работы системы «бетонная плотина-скальное основание», «Гидротехническое строительство», № 4.Д992

32. Бабушкин И.Д., Плотников Н.И., Чуйко В.М. Методы изученияфильтрационных свойств неоднородных пород. М., «Недра», 1974

33. Банничук Н. В., Петров В. М., Черноусько Ф. Л. Численное решениевариационных и краевых задач методом локальных вариаций. Теория вычислительной техники и математической физики. 1966, т. 6, с .6.

34. Банничук Н. В., Черноусько Ф. Л. Вариационные задачи механики иуправления. Численные методы. -М .: Наука, 1973, 283 с.

35. Банничук И. В. Алгоритм метода локальных вариаций для задач счастными производными. - М.: Наука, 1971

36. Баренблатт Г.И., Желтов Ю.П. Об основных уравнениях фильтрацииоднородных жидкостей в трещиноватых породах. - «Докл. АН СССР», 1960, т. 132, №3

37. Белан В.И. Влияние фильтрации на температурный режим мерзлыхбереговых примыканий грунтовых плотин талого и мерзлого типов. Изв. ВНИИГ, сб. научных трудов, том 151, Л., Энергоатомиздат, 1981.

38. Белан В.И. Методика расчета пространственной фильтрации в талойчасти примыкания плотины к мерзлому борту. Изв. ВНИИГ, сб. научных трудов, том 158, Л., Энергоатомиздат, 1984.

39. Белаш П.М., Чен-Син Э., Сенюков Р.В. Метод направленногостатистического поиска при. решении задач разработки нефтяных и газовых пластов (метод Монте-Карло). Тез. Докл. Научн.-техню конф. МИНХиГП, 1964

40. Беллендир Е.Н., Ивашинцов Д.А., и др. Вероятностные методы оценкинадежности грунтовых гидротехнических сооружений. СПб.: Изд-во ОАО «ВНИИГ им. Б.Е.Веденеева», 2003.

41. Белов А.В. Графический способ определения температурныхнапряжений в бетонной стенке. -Труды ЛПИ, 1948, вып.5.

42. Белов А.В. Температурные напряжения в бетонных плитах пригармонических колебаниях температуры. -Л.: Известия ВНИИГ, 1951, т. 45.

43. Белов А.В. К определению температурных напряжении в бетоннойлог плите с учетом экзотермии и теплоизоляции при переменной температуре окружающей среды. Известия ВНИИГ, 1952, т.47.

44. Белов А.В. Температурные напряжения в круглой плите, заделанной пооснованию.-Известия ВНИИГ, т.66,1960.

45. Берг В. А., Решение уравнения теплопроводности методом конечныхразностей. - Информационный сборник Ленгидэпа, Л., No21, 1961.

46. Бильдюг Е.И. Расчет на фильтрацию земляной дамбы с экраном ишпунтом. Изв. НИИГ, 1940, т.28.

47. Биндеман Н.Н. Гидрогеологические расчеты подпора грунтовых вод ифильтрации из водохранилищ., м., Углетехиздат, 1951.

48. Био М. Вариационные принципы в теории теплообмена. -М.: Энергия,1975, 209 с.

49. Битюрин А.К., Соболь С В . Методическое обеспечениеинженерныхрасчетов температурно-фильтрационного режима гидроузлов в условиях вечной мерзлоты на ЭВМ. Гидротехническое строительство, 1993.- №11.

50. Битюрин А.К., Соболь С В . Численное моделирование термокарста.Известия ВУЗов, Энергетика, 1991,- №12.

51. Битюрин А.К., Соболь СВ., Февралев А.В. Температурнофильтрационный режим грунтов долины реки в створе гидроузла. Известия ВУЗов. Строительство, 1997, №7.

52. Биянов Г.Ф. Плотины на вечной мерзлоте.- М.: Энергоатомиздат, 1983.

53. Бетонные плотины (на скальных основаниях). Гришин М.М., РозановН.Н., Белый Л.Д. и др. - М.: Стройиздат, 1975, 352 с.

54. Богословский П.А. Расчет многолетних изменений температурыземляных плотин, основанных на толще мерзлых грунтов.- Труды ГИСИ им. В.П.Чкалова, 1957, вып. 27, с. 123- 178.

55. Богословский П.А. Моделирование температурного режима грунта прифильтрации.- Известия ВУЗов, раздел «Строительство и архитектура», №5, 1959,с.81-89.

56. Богословский П.А. Перспективы исследования предельноготемпературного состояния примыкания земляной нефильтрующей плотины к мерзлому борту долины. -Труды координационных совещаний по гидротехнике /ВНИИГ им.Б.Е.Веденеева, 1977, вып.117, с.67-71.

57. Богословский П.А., Горохов Е.Н., Жданов В.А., Соболь СВ., ФевралевА.В. Прогноз температурно-влажностного режима каменно-земляной плотины. Проблемы инженерного мерзлотоведения в гидротехническом строительстве.- М.: Наука, 1986.

58. Богословский П.А., Соболь СВ., Февралев А.В. Температурный режимгрунтового ядра каменнонабросной плотины Усть-Хантайской гидростанции. Научные исследования по гидротехнике в 1975 году.- Л.: ВНИИГ им. Б.Е.Веденеева, 1976.- Том 1.

59. Борисов А.В. Приближенный метод вычисления напоров в основанииводонапорных сооружений. Тр. объед. семинара по гидротехн. и водн. хоз-ву, вып. 2, Харьков,»Строительство», 1959.

60. Бочевер Ф.М. Расчет подпора подземных вод в берегах на участкахплотин (при наличии оврага в верхнем или нижнем бьефах). Вопросы фильтрационных расчетов гидротехнических сооружений. 1959, №3.

61. Бребиа К., Уокер Применение метода граничных элементов втехнике. Пер. С англ., М.: Мир. 1982.

62. Васильев П.И. Температурный режим массивных плотин. Вопросыпроектирования высоких плотин. - Труды ЛИИ, No251, 1965, с. 19-37.

63. Васильев П.И. Выбор температур замыкания арочных плотин,Гидротехническое строительство, No6, 1965, с.26-39.

64. Васильев П.И. Пути уменьшения температурных напряжений вбетонных плотин. - Труды ЛИИ, No292, 1968, с. 14-29.

65. Ведерников В.В. Методы решения некоторых задач по фильтрации сосвободной поверхностью. Гидротехническое строительство, 1935, №9.

66. Веригин Н.Н. О течениях грунтовых вод при местной усиленнойинфильтрации. Докл. АН СССР, 1950, т.70, №5.

67. Вовкушевский А.В. Расчет фильтрации в основании сослабопроницаемым слоем методом конечных элементов. - Известия ВНИИГ, сб. научных трудов, том 241, 2002.

68. Волынчиков А.Н., Мгалобелов Ю.Б. Обоснование надежности основныхсооружений гидроузла Шон Ла во Вьетнаме., Гидротехническое строительство, 2007, №12, НТФ «Энергопрогресс».

69. Гераськин Н.Н. Расчет пространственного температурного режимаземляной плотины и фильтрующего берегового примыкания. Изв. ВНИИГ, сб. научных трудов, том 151, Л., Энергоатомиздат, 1981.

70. Гинзбург СМ. , Рукавишникова Т.Н. , Шейнкер Н.Я. , Имитационныемодели для оценки температурного режима бетонной плотины на примере Бурейской ГЭС. - Известия ВНИИГ 2002, т. 241, с. 173-178.

71. Гинзбург СМ. , Рукавишникова Т.Н. , Шейнкер Н.Я. , Применениеимитационных моделей при идентификации параметров температурного режима бетонных массивов в стадии возведения. - Известия ВНИИГ 2002, т. 241, с. 178-187.

72. Гинзбург СМ. , Шейнкер Н. Я. , Оценка температурного режимабетонных массивов в процессе их возведения с учетом случайных факторов. - Известия ВНИИГ 2002, т. 241, с. 188-192.

73. Гиринский Н.К. Грунтовый поток земляных плотин. Научн. Зап. МГМИ,1937, вып. 4.

74. Гоголев Е.С Температурное состояние фильтрующей плотины с ядромна мерзлом основании. -Труды V Всесоюзного совещания -семинара по обмену опытом строительства в суровых климатических условиях.Красноярск, 1968, вып. 1, с. 120-135.

75. Гольдин А.Л., Рассказов Л.Н. Проектирование грунтовых плотин. М.;Изд. АСВ, 2001.

76. Горохов Е.Н. Метод расчета температурного режима каменно-землянойплотины с учетом сублимационного ледонакопления в наброске. Известия ВНИИГ им. Б.Е.Веденеева, 1986.- том 188.

77. Горохов Е.Н., Логинов В.И. Программный комплекс «NORD»трехмерного моделирования температурного режима каменно-земляных плотин. Материалы международной конференции «ИГК-2000», ВНИИГ, СПб.,2000

78. Горохов Е.Н. Температурный режим грунтов левобережногопримыкания Вилюйской ГЭС-3. «Гидротехническое строительство», 2003, №2.

79. Гутман Г. Определение тепловых напряжений при гармоническихколебаниях температуры.-Л.:Известия ВНИИГ, 1952, т.47.

80. Давидович В.И. К вопросу о фрагментном способе фильтрационныхрасчетов. Гидротехническое строительство, 1937, №8.

81. Данилова Г.И., Бучко Н.А. Приближенный способ определениятемператур в бетонной кладке строящихся плотин. -Гидротехническое строительство, 1963, №4.

82. Девисон Б.Б. Движение грунтовых вод. В кн.: Христианович А.,Михлин Г. Девисон Б.Б. Некоторые вопросы механики сплошной среды, ч.З, М.-Л.,Изд-во АН СССР, 1938.

83. Девисон Б.Б. Некоторые точные решения задачи о движении грунтовыхвод, получаемые упрощенным приемом. Изв. АН СССР, ОТН, 1938.

84. Дзюба К. И. Исследования лгермонапряженного состояния элементовгидротехнических сооружений, выполненных в НИСе Гидропроекта. труды координационных совещаний по гидротехнике, вып. IY, Л. : Госэнергоиздат, 1962, с. 168-191.

85. Дзюба К. И. Применение численных методов для исследованиятермонапряженного состояния гидросооружений. - труды координационных совещаний по гидротехнике, вып. 103, Л. : Энергия, 1975, с. 9-13.

86. Дзюба К. И., Гунн Я. Расчет толстых прямоугольных плит методамитрехмерной теории упругости. - Труды гидропроекта, сб. No28, М. : 1973, с. 40-45.

87. Дзюба К. И., Деткова М. И. , Прелова Л. Е. , Фрадкина Н. И.Теоретические исследования термонапряженного состояния блоков бетонирования, укладываемых Токтогульским методом. - Труды гидропроекта, сб. No28, М. : 1973, с. 224-232.

88. Иоффе А.Я. К вопросу о фильтрации в неоднородной среде. Ж. тех.физики, 1935, т.5, вып. 1.

89. Исякаев В.А. Решение одной задачи пространственной фильтрацииметодом статистических испытаний. ПМТФ, 1967, №2.

90. Калинин Н.К. Некоторые приближенные приемы решения задачи офильтрации в двухслойной среде. Докл. АН СССР, 1941, т.ЗО.

91. Каранфилов Т.С. О точности фильтационных расчетов на щелевомлотке. Гидротехническое строительство, 1960, №2.

92. Карлслоу Г.С. Теория теплопроводности.- М.- Л.: ОГИЗ, 1947- 288с.Карслоу Г., Егер Д. Теплопроводность твердых тел. - М. :Наука, 1964, с. 488.

93. Клейн И.С. Метод расчета температурного режима каменно-земляныхплотин. Труды ВОДГЕО, М., 1981.

94. Клейн И.С. Численное моделирование оттаивания мерзлых основанийгрунтовых плотин под , действием фильтрации. Инженерное мерзлотоведение в гидротехнике. Материалы конференций и совещаний по гидротехнике.- Л.: Энергоатомиздат, 1989.- с.97-100.

95. Коган Е.А. Плотины из укатанного бетона. Анализ зарубежных данныхо трещинообразовании и рекомендации по обеспечению термической трещиностойкости. Научно технический сборник «Безопасность энергетических сооружений», выпуск 6, АО НИИЭС, М., 2000.

96. Корн Г., Корн Т. Справочник-по математике (для научных работников иинженеров). -М.: Наука, 1977, 831 с.

97. Коротков Б.И. Исследование фильтрации двухшпунтового фрагментапри наличии щели между грунтом основания и флютбетом плотин. Тр. Новосиб. Ин-та инж. жел-дор. Трансп., 1963, вып.36.

98. Кроник Я.А., Кадкина Э.Л., Лосева Г. Анализ состояния грунтоврусловой плотины Усть-Хантайской ГЭС. Гидротехническое строительство, 1978.- №11.

99. Кроник Я. А., Кадкина Э.Л., Скоблин Г. А. Температурный ифильтрационный режим русловой плотины Усть-Хантайской ГЭС в первые годы эксплуатации. Энергетическое строительство, 1980.- №6.

100. Кузнецов Г.И., Кочубиевская Р.Л. Методика и результаты расчетатемпературного режима бурта зимнего хранения суглинка. Изв. ВНИИГ, сб. научных трудов, том 198, Энергоатомиздат, Л.: 1987

101. Кузнецова А.К. Исследования термонапряженного состояния элементовна специализированном интеграторе ЭМБУ-6. -Труды координационных совещаний по гидротехнике, вып. 4, Л., Госэнерго, 1962.

102. Ламкин М.С. Некоторые вопросы практических методов расчетатемператур и температурных напряжений,- Труды координационных совещаний по гидротехнике. Л.: Госэнергоиздат, 1962, вып. 4.

103. Лейбензон Л.С. Движение природных жидкостей и газов в пористойсреде. М.-Л., Гостехиздат, 1947.

104. Ломизе Г.М. Фильтрация в трещиноватых породах. -М.- Л.:Госэнергоиздат, 1951.- 127 с.

105. Лыков А.В. Теория теплопроводности. - М.: Высшая школа, 1976, 599 с.

106. Лукьянов B.C. Гидравлические приборы для технических расчетов,«Известия АН СССР», отд. техн. наук, №2, 1939.

107. Максимов И. А. Расчет промерзания фильтрующего основаниягрунтовой плотины в районах Крайнего севера. Изв. ВНИИГ, сб. научных трудов, том 146, Л., Энергоатомиздат, 1981.

108. Маскет М. Течение однородных жидкостей в пористой среде. М.-Л.,Гостоптехиздат, 1949.

109. Маслов Г.Н. Опыт теоретического анализа температурных измененийкладки плотин и подпорных стенок треугольного профиля.- Известия ВНИИГ, 1936.

110. Маслов Г.Н. Термическое напряженное состояние бетонных массивовпри учете ползучести бетона. - Известия ВНИИГ, 1941, т.28.

111. Микеладзе Ш.Е. Численные методы интегрирования дифференциальныхуравнений с частными производными.- Издательство АН СССР, 1936.

112. Минский Е.М. Статистическое обоснование уравненийфильтрационного движения. Докл. АН СССР, 1958, т.118, №2.

113. Михайлов Г.К. К задаче о фильтрации в анизотропных земляныхплотинах трапецеидального профиля на горизонтальном водоупоре. Доклады АН СССР, 1951, т.80, №4.

114. Михайлов Г.К. К классификации задач теории установившегосядвижения грунтовых вод в вертикальной плоскости. Изв. Ясского политехи, ин-та, 1959, т.5.

115. Можевитинов А.Л., Сулейманов И.А. Критерий необходимости учетанеустановившейся фильтрации в грунтовых плотинах при оценке устойчивости откоса. Известия ВНИИГ им. Б.Е.Веденеева. 1979. т. 135.

116. Мухетдинов Н.А. Термический режим низовой призмыкаменнонабросной плотины. Изв. ВНИИГ, сб. научных трудов, том 90, Л., Энергоатомиздат, 1969.

117. Мухетдинов Н.А. Влияние нелинейной фильтрации воздуха натермический режим каменнонабросных плотин. Изв. ВНИИГ, сб. научных трудов, том 96, Л., Энергоатомиздат, 1971.

118. Недрига В.П. Расчет фильтрации в обход гидротехническихсооружений. Гидротехническое строительство, 1947, №5.

119. Недрига В.П. Расчет фильтрации в основании плотин с учетомводопроницаемости шпунтов. Сб.: Вопросы фильтрационных расчетов гидротехнических сооружений, №2, М., Госстройиздат, 1956.

120. Недрига В.П. Расчет фильтрации в обход гидротехнических сооруженийпри сопряжении их с берегами. Вопросы фильтрационных расчетов гидротехнических сооружений, М., Госстройиздат, 1959.

121. Недрига В.П., Покровский Г.И. Фильтрационный расчет однородныхвзрывонабросных плотин. Тр. Гидропроекта, 1978, №59.

122. Нумеров Н. Приближенные способы расчета фильтрации черезземляные • плотины на водопроницаемых основаниях. Тр. ЛПИ, 1947, №4.

123. Нумеров Н. Обобщение формулы Дюпюи. Изв. ВНИИГ, 1953, т.50.

124. Нумеров Н. Приближенный способ расчета напорной фильтрации восновании гидротехнических сооружений. Изв. ВНИИГ, 1953, т.50.

125. Нумеров Н. Методы исследования плоской установившейсяфильтрации (тяжелой несжимаемой жидкости) в однородной среде. Автореф. докторской дисс. Ин-т механики АН СССР, 1954.

126. Нумеров Н., Барсегян P.M. Об оценке основных допущений методикирасчета фильтрации жидкости в горизонтальных гидравлически связанных пластах. Изв. ВНИИГ, 1965, т.78.

127. Орехов В.Г. Исследование термонапряженного состояния бетонных ижелезобетонных конструкций гидротехнических сооружений. - Труды координационных совещаний по гидротехнике, вып. IY, 1962, Л. : Госэнергоиздат, с. 76-97.

128. Орехов В.В., Хохотва Н. Объемная математическая модельгеофильтрации скального массива, вмещающего подземные сооружения ГЭС Яли во Вьетнаме. «Гидротехническое строительство», 2004, №12.

129. Павловская Л.Н. Моделирование ирригационной сети при исследованииплановой фильтрации методом ЭГДА. Изв. ВНИИГ, 1964, т.75.

130. Павловский Н.Н. Теория движения грунтовых вод подгидротехническими сооружениями и ее основные приложения. Петроград, 1922, Научн-мелиорац. Ин-т; Собр. соч., т.2, М.-Л., Изд-во АН СССР, 1956.

131. Павловский Н.Н. Неравномерное движение грунтовых вод. М.,Гостехиздат, 1930.

132. Павловский Н.Н. О фильтрации воды через земляные плотины. Л., Издво Сектора гидротехн. сооружен., 1931, вып.24.

133. Павловский Н.Н. Гидравлическая теория турбулентного движениягрунтовых вод. Тр. Ноябрьск. сессии АН СССР, Проблемы ВолгоКаспия, 1934.

134. Павловский Н.Н. Гидромеханический расчет плотин системы Сенкова.М. -Л., ОНТИ ГРСЛ, 1937.

135. Пехович А.И. Основы гидроледотермики. - Л.: Энергоатомиздат, 1983,200с.

136. Плят Ш.Н. Расчеты температурных полей бетонных гидросооружений.Издательство "Энергия", 1974.

137. Плят Ш.Н. Метод расчета теплового режима бетонныхгидротехнических сооружений. Известия ВНИИГ, 1964, т. 84, с. 298-342.

138. Плят Ш. Н., Овчиникова Т. Т. Метод расчета термонапряженногосостояния бетонных массивов гидротехнических сооружений в процессе их возведения. - Труды координационных совещаний по гидротехнике, вып. 103, Л. : Энергия, 1975, с. 180-184.

139. Покровский Г.И. Рациональные методы возведения и расчетапротивофильтрационных устройств взрывонабросных плотин. «Гидротехническое строительство», 1996, № 10, М., Энергопрогресс.

140. Полубаринова-Кочина П.Я. О фильтрации в анизотропном грунте.ПММ, 1940, т.4.

141. Полубаринова-Кочина П.Я. О фильтрации в анизотропном грунте.ПММ, 1940, т.4.

142. Полубаринова-Кочина П.Я. О неустановившихся движениях грунтовыхвод. Докл. АН СССР, 1950, т.75, №3.

143. Полубаринова-Кочина П.Я. Теория движения грунтовых вод. Гос.Издательство технико-теоретической литературы, М., 1952.

144. Проектирование и строительство больших плотин. Аварии иповреждения больших плотин / Н.С.Розанов, А.И.Царев, Л.П.Михайлов и др./ Под редакцией А.А.Борового. М.: Энергоатомиздат. 1986.

145. Развитие исследований по теории фильтрации в СССР. Институтпроблем механики АН СССР, Институт гидродинамики СО АН СССР, ВНИИГ им. Б.Е.Веденеева, ВНИИ Природных газов, М., «Наука».

146. Рассказов Л.Н. Схема возведения и напряженно-деформированногосостояния грунтовой плотины с центральным ядром. - Энергетическое строительство, 1977, No2, с.65-75.

147. Рассказов Л.Н., Беляков А.А. Расчет пространственного напряженнодеформированного состояния каменно-земляной плотин. Гидротехниическое строительство, 1982, No2, с. 16-22.

148. Рассказов Л. Н. и др., Гидротехнические сооружения. Учебник дляВУЗов. Под редакцией, М., Стройиздат, 1996.

149. Рассказов Л.Н., Анискин Н.А. и др. Фильтрация в грунтовых плотинах вплоской и пространственной постановке. «Гидротехническое строительство», 1989, №11.

150. Рассказов Л.Н., Анискин Н.А. Фильтрационные расчетыгидросооружений и оснований. «Гидротехническое строительство», 2000, №11.

151. Рассказов Л.Н., Анискин Н.А. и др. Сангтудинский гидроузел:напряженно-деформированное состояние и фильтрация в основании плотины и в обход гидроузла. «Гидротехническое строительство», 2008, №5.

152. Рассказов Л.Н., Анискин Н.А. и др. Конструктивные решения морскихгрунтовых плотин при строительстве ПЭС, Малая энергетика, 2008, №4, ОАО НИИЭС, Москва.

153. Рассказов Л.Н., Орехов В.Г., Анискин Н.А. и др. Гидротехническиесооружения. Учебник для ВУЗов, АСВ, М., 2008.

154. Рассказов Л.Н., Анискин Н.А. и др. Конструктивные решения морскихгрунтовых плотин при строительстве ПЭС, «Гидротехническое строительство», 2009, №7.

155. Ризенкампф Б.К. Гидравлика грунтовых вод. Ч. 1. Уч. зап. Саратовск.Ун-та, 1938, т. 14, вып.1.

156. Ромм Е.С. Исследование фильтрации в трещиноватых породах в связи сих коллекторскими свойствами. - «Тр. ВНИГРИ», 1958, вып. 121.

157. Ромм Е.С. Фильтрационные свойства трещиноватых горных пород. М.,«Недра», 1966.

158. Саваренский Е.Ф. Обобщенная формула Кене и три случая ееприменения. Гидрогеол. и инж. геол., 1936, №2.

159. Савельева Ю.Ю. Использование численного моделированияпространственной фильтрации при обосновании инженернотехнических решений по безопасному обустройству городских территорий. «Гидротехническое строительство», 2003, №9.

160. Сегерлинд Л. Применение метода конечных элементов.-М.: Мир, 1979.

161. Скворцов В.В. Вычисление статистических характеристик дебитовметодами Монте-Карло и малых возмущений. Сб.: Тезисы докладов Юбилейной научн. техн. Конференции, посвящ. 20-летию Казанск. Физ.техн. Инст-та АН СССР, Казань, 1966.

162. Соболь В. Стационарное температурное состояние примыкания кборту долины, имеющей профиль в виде двугранного угла.- Труды координационных совещаний по гидротехнике, ВНИИГ им. Б.Е.Веденеева, 1977, вып. 117, с.74-79.

163. Соболь В. Прогнозирование тепловых и механических процессов восновании и берегах водохранилищ криолитозоны. Известия ВУЗов. Строительство, 1997.- №8, с.36-39.

164. Соболь СВ., Февралев А.В. Температурный режим фильтрующихталиков в основании гидроузла и берегах водохранилища. Известия ВУЗов. Строительство, 1992.- №5-6, с. 106-110.

165. Сольский СВ., Гусакова И.Н. Применение численного моделированиядля расчета фильтрационных полей в основании энергетических объектов. Изв.ВНИИГ: юбилейный сборник, 1996, том 231.

166. Ставровский А.П. Температурное состояние оттаявшей части бортаречной долины в местах примыкания к нему фильтрующей плотины.Известия ВУЗов, раздел «строительство и архитектура».- Новосибирск, 1974, №9, с.96-100.

167. Суровцев Б.П. Напорно-безнапорная фильтрация в основании плоскогофлютбета." Тр. Ташкентск. ин-та инж. ирригации и механизации с. Х.,1961,вып.18.

168. Тихонов А.Н., Самарский А.А. Уравнения математической физики.- М.:Наука, 1977.-736 с.

169. Ухов С Б . Расчет сооружений и оснований методом конечныхэлементов. - М.: МИСИ, 1973, 118 с.

170. Угинчус А.А. Расчет фильтрации через земляные плотины. М.-Л.,Стройиздат, 1940.

171. Угинчус А.А. Расчет фильтрации через земляные плотины. М.-Л.,Госэнергоиздат, 1960.

172. Февралев А.В. Аналитический расчет температурного режимафильтрующего основания бетонной плотины. Известия ВУЗов. Энергетика, 1983, №4.- с.109-112.

173. Февралев А.В. Температурный режим тела и фильтрующего основаниябетонной плотины в системе криволинейных ортогональных координат. Известия ВУЗов.Строительство и архитектура, 1981.- №1.-с. 105-110.

174. Февралев А.В., Янченко А.В. Исследование пространственногонестационарного температурного поля под дном водоемов криолитозоны. Береговые процессы в криолитозоне.- Новосибирск: Наука, 1984.-с. 111-115.

175. Фильчаков П.Ф. Электромоделирование задач фильтрации вразнородном грунте. Докл. АН СССР, 1949, т.66, №4.

176. Фрадкина Н.И. К решению плоской нестационарной задачитеплопроводности на произвольных нерегулярных сетках. Труды координационных совещаний по гидротехнике, вып. 103, л.: Энергия, 1975.

177. Фрид А. Температурные напряжения в бетонных и железобетонныхконструкциях гидросооружений. - М: Госэнергоиздат, 1959, 79с.

178. Фрид А. Расчет изменения температуры бетонных массивов подвлиянием экзотермии цемента. - Известия ВНИИГ, т.41, 1949, с. 67-76.

179. Хованский Г. Основы номографии. Издат-во "наука" 1976, 348 с.

180. Храпков А.А., Цыбин A.M. Некоторые варианты метода дробных шаговдля решения плоской задачи теплопроводности. -Известия ВНИИГ, т.129, 1979.

181. Христианович А. Движение грунтовых вод, не следующее законуДарси. Прикладная математика и механика, вып. 1, т.4, 1940.

182. Цыбин A.M. К созданию метода расчета температурного режимабетонной плотины при наличии фильтрационного потока в основании — Л., Энергия, Труды координационных совещаний по гидротехнике, вып. 103, 1975.

183. Цыбин А. М., Вариационно - разностное решение температурной задачио послойном наращивании системы, состоящей из длинных бетонных блоков. Известия ВНИИГ, 1978, т. 124, с. 42-48

184. Цыбин A.M. Одномерная нестационарная задача фильтрации сосвободной поверхностью. Известия ВНИИГ, 1981, т. 146, с. 26-30.

185. Цыбин A.M. Некоторые вопросы расчета температурных полей,связанные со строительством и эксплуатацией гидросооружений, работающих в районах крайнего севера и вечной мерзлоты.- Петербург: ВНИИГ им. Б.Е.Веденеева, 1995.- 344 с.

186. Цыбин A.M. Температурный режим грунта под руслом реки иопределение границы талика. Гидротехническое строительство в районах вечной мерзлоты и сурового климата. Материалы конференций и совещаний по гидротехнике.- Л.: Энергия, 1979. -с. 36-39.

187. Чарный И.А. Подземная гидромеханика. М.-Л., Гостехиздат, 1948.

188. Чарный И.А. Строгое доказательство формулы Дюпюи для безнапорнойфильтрации с промежутком высачивания. Докл. АН СССР, 1951, т.79, №6.

189. Черноусько Ф. Л. Метод локальных вариантов для численного решениявариационных задач. Теория вычисленной математики и математической физики. 1965, т.5, No4, с.749-754.

190. Чернышев Н. Движение воды по сетям трещин. М.: Недра, 1979

191. Чилингаришвили Г.И. Определение температур в твердеющихпоследовательно уложенных бетонных блоках строящихся плотин. Гидротехническое строительство, 1965, No5, с. 18-23.

192. Чилингаришвили Г.И., Гогоберидзе В.А. Экпериментальнолабораторные исследования термики бетонных плотин. —Труды координационных совещаний по гидротехнике, вып. 103, Л., Энергия, 1975.

193. Чилингаришвили Г.И., Джингарадзе К.В. Определение температур вбетонных плотинах в период строительства с применением ЭЦВМ. Известия ТНИИСГЭИ, 1967, т. 17.

194. Чилингаришвили Г. И. Аппроксимация граничных условий при решениизадач теплопроводности методом сеток для бетонных гидросооружений. - Известия ВНИИГ, 1975, т. 109, с. 109-206.

195. Чугаев P.P. Проектирование подземного контура плотин,расположенных на нескальных грунтах. Изв. ВНИИГ, 1955, т.53.

196. Шаманский В.Е. Численное решение профильной задачи фильтрациигрунтовых вод. Тр. Координац. Совещаний по гидротехнике, 1967, вып.35.

197. Шаманский В.Е. Численное решение задач фильтрации грунтовых водна ЭВМ.- Киев: Наукова думка, 1969.- 374 с.

198. Шейдеггер А.Э. Физика течения жидкостей через пористые среды.Гостоптехиздат, 1960.

199. Шестаков В.М. Фильтрационный расчет земляных плотин и перемычекпри колебании бьефов. Гидротехническое строительство, 1953, №7, с.36-39.

200. Шестаков В.М. Расчет кривых депрессии в земляных плотинах и дамбахпри понижении горизонта водохранилища. Гидротехническое строительство, 1954, №4, с. 32-36.

201. Шестаков В.М. Определение гидродинамических сил в земляныхсооружениях и откосах при падении уровней в бьефах. Сб. «Вопросы фильтрационных расчетов гидротехнических сооружений», ВОДГЕО, 1956, №2.

202. Шестаков В.М. Некоторые вопросы моделирования неустановившейсяфильтрации. Сб. «Вопросы фильтрационных расчетов гидротехнических сооружений», ВОДГЕО, 1956, №2.

203. Шестаков В.М. Расчет фильтрации через земляные плотины напроницаемом прослое. Гидротехническое строительство, 1957, №1.

204. Шестаков В.М. Оценка сопротивления ложа водоемов пригидрогеологических расчетах. Разведка и охрана недр, 1964, №5.

205. Шорин Н. Теплопередача. - Гос. Издательство по строительству иархитектуре. - М. - Л., 1952, 335 с.

206. Штейнберг В.М., ПрокоповичИ.Е., Гольдфарб И.В. Практический методисследования температурного поля в теле массивных бетонных гидротехнических сооружений. - Сб. трудов Одесского гидротехнического института, Одесса, 1953, вып.4.

207. Шугаева Р.Т. Расчет термического режима земляных плотин с учетомпослойного их возведения. Изв. ВНИИГ, Сб.научных трудов, т. 96, Энергоиздат, 1971

208. Шугаева Р.Т. Учет фильтрации в расчетах на ЭЦВМ термическогорежима неоднородных грунтовых плотин. Изв. ВНИИГ, Сб.научных трудов, т. 151, Энергоиздат, 1981

209. Шугаева Р.Т. Прогноз термического режима грунтовой плотиныВилюйской ГЭС-Ш. Изв. ВНИИГ, Сб.научных трудов, т. 158, Энергоиздат, 1984

210. Шугаева Р.Т., Распопова Р.Х. Прогноз термического режимафильтрующего основания грунтовой плотины, промораживаемой парожидкостными охлаждающими устройствами. ВНИИГ, Сб.научных трудов, т. 178, Энергоиздат, 1985

211. Шугаева Р.Т., Шугалей Р.Т. Изменение термического состояниямноголетнемерзлого массива под влиянием экзотермии хладостойкого инъекционного раствора. Изв. ВНИИГ, Сб.научных трудов, т. 198, Энергоиздат, 1987

212. Шульц Г.Л. Расчет температуры трехмерной области методом конечныхэлементов. - Материалы конференций и совещаний по гидротехнике.Л.: Энергия, 1979, с. 70-73.

213. Щелкачев В.Н., Лапук Б.Б. Подземная гидравлика. М.-Л.,Гостоптехиздат, 1949.

214. Эйдельман Я. Натурные исследования плотины Братской ГЭС. - Л. :Энергия, 1968, 253 с.

215. Юшков П.П. Приближенное решение задач нестационарнойтеплопроводности методом конечных разностей. Труды Института энергетики АН БССР, 1958, вып.6.

216. Яковлева М.В. Исследование температурного поля в массивном бетоне.Инженерно-физический журнал, 1960, т.З, №11.

217. Янченко А.В. Метод расчета пространственного нестационарноготемпературного режима грунтов ложа водохранилища и основания плотины. Известия ВУЗов. Строительство и архитектура, 1983.- №10.с.93-96.

218. Выполнение работ по инженерным изысканиям для обоснованиярабочей документации строительства Сангтудинской ГЭС-1 на р.Вахш в Республике Таджикистан. 1. Инженерно-геологические изыскания. Этап

219. Информационный отчет за 2-ой квартал 2008 г. Филиал ОАО«Инженерный центр ЕЭС» -«Институт Гидропроект», М., 2008 г.

220. Заключение по результатам- производственных работ по созданиюпротивофильтрационного элемента методом струйной цементации в основании руслового участка ядра плотины. Филиал ОАО «Инженерный центр ЕЭС» -«Институт Гидропроект», М., 2007 г.

221. СНиП 2.06.06-85. Плотины бетонные и железобетонные. М:,ЦИТПГосстроя СССР, 1986.

222. Brebbia А. Finite Element Systems, Springer-Verlag, ComputationalMechanics Centre, Soutgempton, 1982.

223. Carrere A., Hamon M., Poniet P. Three-dimensional finite element analysesof the Laparan Dam.- Water Power and Dam Constraction, 1983, August.

224. Cedergreen H.R. Investigation of drainage rates affecting stability of earthdam. Discussion. Proceeding of ASCE, February, 1948.

225. Cedergreen H.R. Seepag, drainage and flow nets. New-York - LondonSidney-Toronto, 1977

226. Darcy H. Les fontaines publiques de la ville de Dijon. Общественныеколодцы в городе Дижон, Paris, 1856

227. Dachler R. Uber Stekerwasserstromungen in geschichtetem Material.Фильтрация воды в слоистом материале. Wasserwirtschaft, 1933, N 2.

228. Dolegalova М., Shvanchar J. Analysis of unusual behavior and stability ofZermanice dam. Cheh National Committee on Large Dams, 2006.

229. Fanelli M., Giuseppetti G., Numerical analysis of the thermal state of a dam,Num. Meth.Theem. Probl., Pt.2, Proc. 4-th Int. Conf. Swonsea, 1985.

230. Irmay S. Flow of liquid through cracked media. Bull. Rec. Council of Israel,1955, v. 5a, №1. a/

231. Hamel G/ Uber Grundwasserstromung. О течении грунтовых вод. ZAMM,1934, Bd. 14

232. Musket M. The seepage of water throught dams with vertical faces.Фильтрация воды через дамбы с вертикальными стенками. Phisics, 1935, v.6

233. Sherard J.J., Woodward R.J., Earth and earth-rock dams. J.Wiley&sons. Inc.1967.

234. Snow D.T. Anisotropic permeability of fractured media. Water Resour. Res.1969,5(6).

235. Reinius E., The stability of the upstream slope of earth dams, Stockholm,1948.

236. Visual MODFLOW v.4.0. User's Manual. For Professional Applications inThree-Dimensional Groundwater Flow and Contaminant Transport Modeling. 2004, Waterloo Hidrogeologik Inc.

237. Wittke W. Three-Dimensional Percolation of Fissured Rock. Planning openpit mines. Proc. of the Symposium on the theoretical Background to the planing of Open Pit Mines with special reference to slope Stability. 1.hannesburg. 1970.

238. Wittke W., Louis C. Zur Berechnung des Einflusses der Bergwasserstromungauf die Standsicherheit von Boschungen und Bauwerken in zerkluftetem Fels. Proc. I. Congresses Int. Soc. Of Rock Mechanics. Vol. II, Lisboa, 1966.

239. Wittke W., Louis C. Modelversuche zur Durchstromung kluftiger medien.Felsmechanik und Ingenieurgeologie, 1968, Suppl. 4.