автореферат диссертации по энергетике, 05.14.09, диссертация на тему:Русловые водные балансы за короткие интервалы времени

РГБ ОА

« «.к! '.ЧПЯ

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

ЗАЖИМЛРОВ Сергей Анатольевич

РУСЛОВЫЕ ВОДНЫЕ БАЛАНСЫ ЗА КОРОТКИЕ ИНТЕРВАЛЫ ВРЕМЕНИ

Специальность 05.14.09 - Гидравлика и инженерная гидрология

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

САНКТ-ПЕТЕРБУРГ

1995

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

Научный руководитель - доктор технических наук,

профессор И. Ф. КАРАСЕВ

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

доктор технических наук,

профессор Коваленко К В. ,

кандидат технических наук Лобанов Е А.

Ведущая организация - МосЦГКС

Защита диссертации состоится "30 " мач 1996 г. в II час на заседании диссертационного совета Государственного гидрологического института по адресу: 199053, Санкт-Петербург, 2-я линия, д.23 Сакс: (812) 213-10-28

С диссертацией можно ознакомиться в библиотеке Государственного гидрологического института ®

Автореферат разослан " 28 " апреля 1995 года

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

совета, кандидат технических наук (^огучюм^и^/^алкшIый И. Л

- а -

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

Цель работы - методическое обоснование наиболее важных аспектов расчета РВБ за короткие интервалы времени. К ним относятся:

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

2. Проверка надежности расчета баланса в целом и его отдель-пых элементов.

3. Качественная и количественная оценка генезиса невязок баланса.

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

Практическая ценность. Применение предложенных методов позволяет существенно повысить достоверность расчета суточных рус-

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

Апробация работы. Основные положения диссертации докладывались на V Всесоюзном гидрологическом съезде (1988), научных семинарах (1987-1994), конференции молодых ученых и специалистов (1992), итоговой сессии Ученого совета ГГИ (1994).

Публикации. Основные положения работы нашли отражение в трех опубликованных статьях и двух находящихся в печати.

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

Структура и объем работы. Диссертация состоит иэ введения, шести глав и заключения на 157 страницах машинописного текста, иллюстрирована 9 рисунками и 26 таблицами. Список литературы включает 117 наименований.

Содержание работы

В первой главе дается краткий очерк развития метода РВВ. современное состояние вопроса, постановка задачи исследований. Большой вклад в совершенствование метода внесли Вольфцун И. Б.. Устюжанин Б. С., Карасев И.Ф.. Савельева А. В. В частности, благодаря их усилиям подготовлены и изданы Штодические указания, положенные в основу расчетов РВВ за декадные и месячные интервалы времени, выполняемых в УГШ.

В 80х-90х годах ГГИ проведены гидрометрические исследования рек Куры, Москвы, Нижнего Дона. Амударьи, Маныч с участием автора. Существенной частью этих исследований являлась увязка стока по длине реки и метрологическая оценка гидрометрических данных.

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

чета русловых водных балансов за короткие интерваш времени.

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

Qb - Qh + Q6n + Qck + Qcb - Obs - Qm? + Qoo + Qcp + Qhii +

± Qpp ± Q6p t Оф ± Qa ± Oo - 0, ( 1 )

где Qb и Qh - расходы воды в верхнем и нижнем створах; Q6n и Qck - боковой и склоновый приток; Qcb и Qbs - сбросы и водозаборы воды; Qht - потери на испарение и транспирацию; Qoo и Qcp - поступление от осадков и таяния снегозапаеов в русле; Онп - приток напорных подзекшых вод; Qpp и Q6p - русловое и береговое регулирования; 0ф - фильтрационный отток или приток подземных вод в русло при установившемся режиме; Ол - потери на ледообразование или поступление от таяния льда; Qo - остаточный члзн уравнения, компенсирующий неучтенные элементы РВВ и совокупную погрешность учтенных.

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

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

т.е. тем или иным методом выполнять расчет неустановившегося движения воды. Для короткого интервала времени уравнение РВБ можно записать в следующем виде:

Кг ОвЛ-Т - <}нЛ + £ + <5оЛ - 0. ( 2 )

¿•1

где Ь - расчетный момент времени; т - время добегания расхода воды на участке: - сумма неосновных элементов баланса

(основными будем считать данные учета стока и русловое регулирование); Кт - коэффициент трансформации волны паводка,

Кт _ а«,* - . ( з )

О-а.г-т

Такая интерпретация РВБ дает возможность:

- оценивать трансформацию волны паводк^ на участке реки;

- контролировать данные учета стока или элементы, входящие в 2 (?1( О, в том числе не поддающиеся непосредственным измерениям;

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

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

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

решения конкретной задачи.

В данном исследовании на натурных материалах специальных наблюдений ГГИ в нижнем бъейЬе ГЭС на р. Тверце (1959) сделана попытка оценить область применения и погрешность расчета наиболее распространенных методов, основанных на упрощенных уравнениях движения. Проверялись следующие методы: Маскингам, метод кривых добегания (интеграл Дюамеля). методы, основанные на уравнениях диффузионной и кинематической волны.

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

Он - <3в + Орр + Опр , ( 4 )

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

дрр „ --- I у ( б )

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

<6 >

зависит, по крайней мере, от двух факторов:

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

- параметра нестационарности (локальная характеристика)

'«** " 773Г Й • < 8 >

где са

волновая скорость, м/с; 2 - отметка уровня воды. м.

Наилучшей формой связи Qnp - f( о(.н), оказалось уравнение регрессии

Qnp - аД^-)г+ аг/оСн/ . ( 9 )

Параметры уравнения ( 9 ) рассчитаны по результатам экспедиционных наблюдений ГГИ на реках Кура (1984), Москва (1986), Дон (1988): а^ 23,2, аг- 2,1.

Проверка методов по материалам наблюдений на р. Тверце выполнялась как для участков, ограниченных соседними створами, так и для общего участка (I-VI) протяженностью 27.3 км. Первый участок представляет собой канализированное русло, следующие - естественное русло с поймой, причем, если на втором располагается высокая пойма, то, начиная с третьего, - пойма низкая. Боковая приточность на всех участках практически.отсутствует. Нестационарность водного режима оценивалась с помощью коэффициента н. Численные эксперименты относились к достаточно разнообразным условиям: отсутствие и наличие поймы, резкая нестационарность потока, наблюдаемая на первых трех участках (o¿h > 0,25) и слабо выраженная - на четвертом и пятом (<* н < 0,25).

Расчет производился в два этапа: 1) по измеренным данным трех попусков (8-9.08, 10-15.08, 15-16.08) подбирались оптимальные значения параметров каждой модели (критерий - минимальное значение среднего квадратического отклонения рассчитанных гидрографов от наблюденных); 2) для контрольных попусков 7- 8.08 и 16-17.08 по гидрографам входных створов и принятым параметрам рассчитывались расходы воды в замыкающих створах каждого участка. Результаты приведены в табл. 1 и на рис.

Апробация балансово-динамического метода (БДМ) по данным наблюдений на р. Тверцы позволила в результате численных экспериментов установить оптимальные значения коэффициентов а( и а* для разных условий нестационарности (табл.2).' По мере уменьшения нестационарности потока, параметры а« и a¿ закономерно стремятся к значениям, полученным при разработке метода (23.2 и 2.1). При выраженной нестационарности потока форма зависимости Qnp -<¿ н) изменяется, в этом случае параметры а< и аг следует рассматривать, как функции от <t н.

Таблица 1

Средние квадратические отклонения рассчитанных гидрографов р. Тверцы от измеренных (%)

метод расчета

м/и Нас кингам Кривых добег. да КВ БДМ РВВ

ств.

1 2 3 4 5 6

ГЬпуск 7-8.08

1-2 14.9 71.1 1.0 13.0 . 24.8 31.2

2-3 25.7 44.8 9.8 25.2 11.9 28.8

3-4 25. 4 24.6 15.5 48.5 19.1 23.0

4-5 6.6 10.2 13.5 13.8 9.2 7.0

5-6 8.5 8.4 2.5 16.3 8.4 6.7

1-г 15.2 49.6 25.1 33.8 19.3 31.2

Попуск 16-17.08

1-2 6.4 43.1 0.7 5.7 7.0 5.0

2-3 8.9 30.5 6.0 15.4 9.6 га 8

3-4 17. 4 19.1 9.3 32. 8 12.5 16.4

4-5 5.9 10.4 10.5 15.4 5.6 5.7

5-6 4.5 1.6 7.4 5.5 4.0 1

1-6 9.4 37.4 • 19.5 35.1 20.6 30.1

Весь период (6.08-17.08) 1

1-6 7.2 35.4 20.1 ..... 36. 3 1 21.3 " 32.2

0 п*/с

/го

100 ■ егн

80 ■ \

60 . \

АО .

20 ■ Л \

р.Таерца. Гидрографы попуска 7.08-8.08 <3м3/с и

а)

10 0 « .о * 1» Л £ я VI»

80-

о/

60 - 4

40.

Ц «V/

го. чо

80 70 60 50 АО 30 20 10

г)

/ А' -А '

/У' V *

( 7 / А' 'А

Ф

'Р > О

•

Рхл чь* V

£ *л

о чг>к

, л

у

« Л

I >

12

18

24

12

18

* V)

- измеренные данные

у«** балансово-динамический метод о о о о Маскингам — метод кривых добегания

--- кинематическая волна

ллдА диффузионная волна

12 1в 4 час

Расчетные участки: 0.4-1.4 км (а), 1.4-4.9 км (б), 4.9-П.8 кн (6) П.8-20. 6 км (г), 20.6-27.8 кн (д) от НовотвереЦкон ГЭС.

Таблица 2

Оптимальные значения коэффициентов а< и для р. Тверды.

Дата попуска

ограничивающие ствош 1-П II-III III-IV 1У-У V- VI

а< а? а^ а< а^ а< а^ а< аг

8-9.08 1.1 0.9 1.2 1.0 3.3 1.1 10.3 2.2 23.5 2.2 10-15.08 1.2 0.8 5.1 1.2 20.2 1.9 23.0 2.3 25.0 2.5 15-10. СО 0.7 1.0 1.0 1.0 3.1 1.1 21.0 1.0 22.0 2.2

Для сравнения выполнялся расчет по уравнению руслового водного баланса для бесприточного участка в его традиционном виде:

Изменение объема возы на участке ЛУ определялось о дсот:;/ точностью - по материалам аооо^/отое'.ч'мки. Установлено, ето л ли пойменных участков пои резко вкгагдкиом нестационарно;,: рг;та»о (участок £ и 3) замена уравнения ¡«разрывности на уравнение суслового водного баланса ( 10 ) приводит к существенной лоп^шно'.-ти расчета (см. табл. 1). В то же время для этих ме условии балан-сово-динамический метод (иначе, уравнение руслового годного баланса с учетом пойменного регулирования) дает значительно лучине результаты, что свидетельствует о его эффективности. С уменьшением нестацлонарности движения (по мере удаления от ГЭС) результат расчетов этими методами становятся близкими. Для участков 4 и 5 относительное расхождение значений Опр вычисленных по данным аэрофотосъемки и соотношению ( 4 ) составило 3-10 т. е. в этом случае новый метод позволяет непосредственно оценивать пойменное регулирование, для определения которого обычно требуется значительный объем детальной морфометрической информации.

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

Он - Ов + #

■

{ 10 )

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

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

1) выраженная нестационарность потока (ан > 0.25)

- уравнение ДВ (5-16 %) и БДМ (10-20 %);

2) умеренная нестационарность потока (обн < 0.25)

призматического русла

- уравнение ДВ (2-5 %), а при н < 0.11 - допустимо применение интеграла Дюамеля (5-10 7.) и уравнения КВ (5-15 %);

непризматтесного русла

- при наличии детальной морфометрической информации выполнять расчет по уравнению РВБ (5-10 %):

- при отсутствии морфометрических данных наиболее эффективны уравнение ДВ (2-5 X), метод Маекингам (5-10 X) и БДМ (5-10 X).

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

Проверка' методики, выполненная по данным наблюдений на реках Тверца, Москва' и Дон, позволила получить оценки точности определения уровня и обгема воды, а также руслового регулирования при различной дискретизации русла системой поперечников и уровенных постов. Для сравнения использовались результаты линейной интерполяции. При этом установлено, что:. - при отсутствии данных о профиле водной поверхности сплайн-интерполяция отметок уровня воды точнее линейной, так как она использует данные по участкам, сопряженным с расчетным;

- наличие переломов водной поверхности, не учитываемых при ее аппроксимации (а такие переломы г(х) обнаружены практически на всех расчетных участках), приводит к существенной погрешности расчета V и суточных Орр, заметно возрастающей в паводочный период. Наибольшие значения погрешности определения емкости (^У) и суточного руслового регулирования (<У(}рр) имеют место на участках с относительно большим уклоном. Так, для верхнего течения р. Юсквы (I - 0,2 - 0,3 %..) они в отдельных случаях достигали соответственно 40 X и 39 м3/с:

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

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

Проверка предложенных соотношений проводилась по материалам исследований на реках Тверце (1959), Москве (1986), Лмударье (1989, 1990). Для сравнения использовались результаты расчета по измеренным морфометрическим характеристикам участков.

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

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

1. Результирующая погрешность

б0 -УГ57

/ до / < Бо -

( и )

где N - число элементов РВБ: ] - порядковый номер элемента; 6а, 6¡ - средние квадратические погрешности определения соответственно остаточного члена баланса и его отдельных состав ля ющх; So - прег^льная (допустимая) погрешность расчета: tp - квантиль-ный множитель, при доверительной вероятности р - 0,95 и нормальном распределении совокупной погрешности tp - 1,96.

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

2. Величина Ро - фактическая обеспеченность допуг- ямой погрешности расчета, иначе - отношение числа расчетных интервалов времени (п - п0), в течение которых невязки не превышали So, к общему их числу (п)

п ~ „ t л о \

Ро - —— . Ро > р , ( 12 )

где р - заданная доверительная вероятность расчета допустимой погрешности расчета (So).

3. Центр распределения выборки Qo, оцениваемый средней арифметической величиной

п

rtiflo - Qb - . /«е./ * sm - tp ^ . ( 13 )

где í - порядковый номер расчетного интервала времени: Sm - предельная погрешность определения т^. Несоблвдение ( 13 ) свидетельствует о значимости систематической составляющей погрешности расчета за рассматриваемый период.

4. Рассеяние невязок баланса, оцениваемое величиной среднего квадратического отклонения

jíQo.i - Q,)2 '

n-i- • б* • < 14 >

где ás - среднее квадратическое отклонение предельных погрешностей расчета < * —

При выполнении условия ( 15 ) будем считать, что учтены все зна-

- 1Б -

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

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

При сравнении двух вариантов РВБ в качестве целевых используем следующие функции:

п -(ро-р^)-б* . Гг-/т<!а/ Гз-бо.-й! (16-18)

Индексы 1 и 2 означают номер варианта расчета

Значения А* - назначаем в зависимости от поставленной задачи. Более точным признаем расчет максимизирующий взвешенную сумму. Эту схему можно применять для решения следующих задач.

1. Оценка надежности и детальности расчета балансз.

В этом случае характеристики р„, тс,„. 60о искомого варианта сравниваются с их предельными значениями, вычисленными пои проверке критериальных условий ( 12 ), ( 13 ), ( 14 ). При этом весовым коэффициентам присваиваются следующие значения: - 0.40, А* - 0.15, Аз — 0.45. полученные по результатам расчетов РВБ на участках реки Амударьи (1981-90). Различие значений А« вызвано тем обстоятельством, что при неполном балансе уточнение одного из элементов иногда приводит к увеличению шо„ при уменьшении <$ао и увеличении р0. В связи с этим вес Тг - наименьший. В случае выполнения условия ( 13 ) коэффициент полагается равным нулю.

Составление РВБ считаем надежным, а его расчетное уравнение достаточно детальным, если сумма 01 -2«".. АКГ* >0.

2. Выбор эффективных методик расчета элементов РВБ.

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

Методику определения рассматриваемого элемента признаем более эффективной, по сравнению с другими вариантами расчета, если ее использование приводит к максимизированию суммы 02 -3. Проверка значимости неучтенной составляющей. • Неучитываемый ранее элемент РВБ признаем значимым, если при включение его в расчетное уравнение соблюдается условие:

Ю- 1.966«,) > 0, (19)

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

б11 ( 20-22 )

При выявлении источников систематических погрешностей расчета элементов РВБ применен корреляционно-регрессионный анализ связей остаточного члена баланса с его элементами 0о-Г((^), для которых характерна взаимная коррелированность, нестационарность и нелинейность.

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

Щ - Г' - Г*""'« ==; , ( 23 )

Vа-ггю - пи'

где - частный коэффициенты корреляции I - го элемента о Оо при устранении влияния Ок; г; , г* , г;„ - парные коэффициенты корреляции.

Значимость 1?} устанавливается.посредством известных соотношений с учетом длины выборки. Считается, что элементы имевшие значимые частные коэффициенты корреляции Я] рассчитываются с систематической погрешностью. На основе.анализа значимых Я] элементы ранжируются по вкладу в остаточной член баланса

Устранение нестационарности связей достигается вы-

делением периодов, однородных по генезису источников систематических погрешностей. С этой целью предлагается использовать критериальное условие (?о < Бо. а уточнение границ выделенных периодов выполнять, используя результаты корреляционного анализа выборок Оо фиксированной длины (20-30 суток) с последовательным их смещением на одни сутки.

Количественная оценка систематических погрешностей выполняется методом регрессионного анализа связей 0о-Г(0]). выделенных на основе сравнения коэффициентов корреляции Из . На значения ¡коэффициентов регрессионных уравнений накладываются ограничения. исходя из знания физически реальных пределов изменения систематической погрешности данного элемента. При линейных связях Оо-Г(оз) общее уравнение регрессии имеет вид:

00 > I (33 + в , ( 24 )

где Л} - коэффициенты уравнения; - элементы баланса, имеющие

значимые частные коэффициенты Р); е - случайная ошибка. В этом случае

т, - Л < с; . ■ ( Ж )

где т, - относительные систематические погрешности: с; - предельно возможные относительные систематические погрешности. При нелинейности искомых зависимостей

<Эо -¿ЛГШ + а . ( 26 )

- —— < с] . { ¿г )

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

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

баланса Пролетарского водохранилища (р.Маныч). Б данных случаях применение предложенных методических' разработок позволило:

- оценить эффективность выполненных уточнений расчета стока:

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

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

- существенно повысить надежность расчета РВЕ и ВБ за короткие интервалы времени (табл.3).

Таблица 3

Статистические характеристики суточных РВБ

Река Участок Период Вариант ОН Шв. Ро 01

РВБ м3/с м*/с и3/ с М3/с

Москва Можайский ГУ IV. 86 цгмс 165 1.8 25,0 24.4 0.63 -3.4

-Нестерово ГГИ 151 0.0 11,0 23.7 0.87 4.9

Нестерово- ЦГМС 180 -4,8 19,1 24,9 0.80 1,1

Звенигошд гги СО к* 0.2 12,9 24,5 0,87 4.4

Звенигород- цгмс 241 -15,5 21.1 28; 6 0,77 0.1

Петрово-Дал. гги 259 -0.5 10,9 28.9 0,97 8.7

Ыаныч а Егорлык- Ш-ХП ДБВО 40 0,6 12,1 8.7 0,58 -2.8

Пролет. ГУ 1993 гги 44 2.0 8.2 8.8 0,76 -1.1

В заключении приведены основные результаты и выводы:

1. Установлены области применения упрощенных методов расчета неустановившегося движения и количественные оценки'погрешности расчета. •

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

3. Для случая отсутствия детальной морфометрической инйюрма-

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

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

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

6. Обосновано применение корреляционно-регрессионного анализа пои выявлении источников систематических погрешностей расчета элементов РВВ и их количественной оценки.

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

По теме диссертации опубликованы следующие работы:

1. Зажимаров С. А. Приближенный расчет трасЮрмации волн паводка па участках русловых водньг/ балансов. - Труди V Веоео.тан. гидрол. съезда, Т.З, 1988, с. 253-258.

2. Яковлева Т. 'Л. . Зягамаров С. А. Оперативно-прогностические мод ••ли учета стока на гидрологическом посту р.Дон - ст. Раздорс-кая. - С-П.: ГГИ. Деп. в ВИНИТИ М 2696 - Е04, 1094. - 13 с.

3. Зажимаров С. А. Критерии и условия надежности расчета русловых водных балансов. - С-П.: ГГИ. Деп. в ВИНИТИ N 2695- В94. 1994. - 7 с.

4. Зажимаров С. А. Сравнительный анализ и область применения упрощенных методов расчета неустановившегося движения воды. - СП.: РГГМИ (в печати).

5. Зажимаров С. А. Применение одномерных сплайнов при расчете русловых емкостей и суточного руслового регулирования. - С-II: РГГМИ (в печати).

Ротп.ААНИИ.Зак.22-100 Уч.ИЗД*.Л.0.8 2J.0-1.35