Определение водохозяйственных балансов

СОДЕРЖАНИЕ

ВВЕДЕНИЕ

1 ОБЩАЯ ХАРАКТЕРИСТИКА СТВОРА 3

1.1 Краткая характеристика природных условий района проектируемого водохранилища 3

1.2 Гидрологическая изученность и особенности гидрологического режима р.Иртыш 4

1.3 Составление батиграфических характеристик водохранилищ 5

2 ГИДРОЛОГИЧЕСКИЕ РАСЧЕТЫ 8

2.1 Подбор теоретической кривой обеспеченности 11

2.2 Внутригодовое распределение годового стока 13

3 ВОДОПОТРЕБИТЕЛИ И ВОДОПОЛЬЗОВАТЕЛИ 15

3.1 Расчет годовых объемов водопотребления 15

3.2 Расчет объема возвратных вод и вод, необходимых для разбавления сточных вод 18

3.3 Расчет объемов попусков 18

4 ВОДОХОЗЯЙСТВЕННЫЕ БАЛАНСЫ (ВХБ) 21

4.1 Годовой ВХБ и способы его увязки 22

4.2 Составление месячного ВХБ 24

4.3 Аппроксимация кривой Q = f (H) 24

4.4 Расчет полезной емкости водохранилища 26

4.5 Варианты возможных ВХБ 27

5 ВОДНОЭНЕРГЕТИЧЕСКИЙ РАСЧЕТ 36

6 ПРОГНОЗ ПОДПОРА ГРУНТОВЫХ ВОД 39

ЗАКЛЮЧЕНИЕ………………………………………………………………………..43

СПИСОК ИСПОЛЬЗУЕМОЙ ЛИТЕРАТУРЫ……………………..……………….44

ПРИЛОЖЕНИЯ 

 

     ВВЕДЕНИЕ

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

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

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

 

  1.   ОБЩАЯ ХАРАКТЕРИСТИКА СТВОРА
    1. Краткая характеристика природных  условий района проектируемого водохранилища
    
  • Река Васюган  расположена на юге Западно-Сибирской равнины, левый приток Оби. Река Васюган протекает на территории Западной Сибири и Алтая. По характеру рельефа большая часть территории представляет пониженную, слабо дренированную равнину, на которой чередуются самостоятельные низменности и возвышенности. К юго – востоку  Западно – Сибирская низменность сменяется зоной предгорных расчлененных равнин, переходящих в горные хребты Салаирского кряжа, Кузнецкого Алатау и Горной Шории, окаймляющие Кузнецкую котловину.
  •     
  • Река Васюган (с. Средний Васюган) является притоком реки Оби, впадает в неё на 2169-м км от устья. Длина реки 820км, площадь водосбора 31700км2. Васюган берет начало в болотах Обско-Иртышского водораздела, имеет равнинный характер, очень извилиста, бассейн сильно заболочен.
  •     
  • В верховьях  река протекает среди низких заболоченных берегов. Пойма изобилует множеством озер и стариц. В среднем течении  высота берегов увеличивается, местами  до 50м. В нижнем течении река становится многоводной, ширина её увеличивается  до 600м. Пойма на этом участке  широкая, с большим количеством отмелей, перекатов и проток.
  •           В пределах Васюганского плато неогеновые и палеогеновые породы выходят на дневную поверхность или перекрываются обводненными четвертичными песками и галькой мощностью до 30-50 м. Питание реки происходит за счет снеготаяния  и дренирования грунтовых вод четвертичных и более глубоких водоносных горизонтов. Величина подземного стока в пределах Васюганского плато составляет 20-40%  общего речного, а коэффициент подземного стока - 5-10%. Модуль подземного стока - 1,5-2,0л/сек.км2.

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

        

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

  • Самый холодный месяц – январь (средняя температура  от –18,0 ºС до –22,0 ºС), а самый тёплый – июль (средняя температура которого колеблется в пределах 17,0–19,0 ºС).
  •     

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

      1. Гидрологическая изученность и  особенности гидрологического режима р. Васюган
     

         В данной курсовой работе рассматривается р. Васюган–с. Средний Васюган.

         Наблюдения  за уровнями и расходами в створе водохранилища продолжались с 1936 по 1962 гг.

         Размер  водозаборной площади составляет 31700 км2.

         Средняя высота водосбора – 110 м.

         Расстояние  от истока реки до створа водохранилища  – 812км.

         Средний уклон реки – 0,09%.

         Средневзвешенный  уклон реки – 0,07%.

         Лесистость  – 58%.

         Озёрность – 2%.

         Заболоченность  –40%.

         Вариант назначения водохранилища -  I . 

      1. Составление батиграфических характеристик водохранилищ
     

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

       Объем воды в водохранилище находят  суммированием частичных объемов 

    ΔVi, заключенных между смежными горизонталями.  Объёмы слоёв воды в проектируемом водохранилище между двумя соседними горизонталями рассчитываются либо по формуле:

                                                                                  (1.1)

       либо  при  Ωi+1/ Ωi > 1,5 по зависимости

        ,                          (1.2)

       где Ωi и Ωi+1 – площади водной поверхности, соответствующие уровням Hi и Hi+1, км2;

             ΔH – вертикальное расстояние между горизонталями, м.

             Объём воды под первой горизонталью вычисляется по формуле:

        ,                                               (1.3)

       где Ω1 – площадь зеркала водохранилища при отметке уровня воды на первой горизонтали, км2.

       Интегральные  суммы объёмов в млн м3, соответствующие всем уровням воды Hi, находят по формуле:

        .                                                 (1.4)

       

     м3.

       Важными характеристиками водохранилищ являются также средняя глубина воды в  водохранилище:

                                                               (1.5)

       и критерий площади литорали (мелководья):

        ,                                                      (1.6)

       где Vi и Ωi – объём и площадь при уровне Hi;

             ΩLi – площадь литорали (мелководья), км2, вычисляемая по формуле:

        ,                                                (1.7)

       где ΩH-2 – площадь водной поверхности, соответствующая уровню воды (Hi–2 м).

             Результаты расчётов заносятся в таблицу 1.1 и затем  строятся соответствующие кривые зависимостей hср = f(H), Li = f(H), Ω = f(H), V = f(H)   (Рисунок А1).     
     
     
     
     
     

       Таблица 1.1 - Батиграфические характеристики водохранилища 

    Hi, м Ωi, км² Ωср, км² ∆Vi, млн.м³ Vi, млн.м³ hср, м Ω(Hi-2), км² ΩLi, км² Li
    1 2 3 4 5 6 7 8 9
    0       0 0 0 0 0
    5 0,5 0,25 1,67 1,67 3,33 0,2 0,3 0,6
    10 2 1,25 12,5 14,17 7,08 1,1 0,9 0,45
    15 8,5 5,25 78,75 92,92 10,93 5,3 3,2 0,38
    20 22 15,25 305 397,92 18,09 15,2 6,8 0,31
    25 45 33,5 837,5 1235,42 27,45 35,2 9,8 0,22
    30 62 53,5 1605 2840,42 45,81 51,6 10,4 0,17
    35 180 121 4235 7075,42 39,31 128 52 0,29
    40 280 230 9200 16275,42 58,13 220 60 0,21
    45 350 315 14175 30450,42   312 38 0,11
     
     
     
     
     
     
     
     
     
     
     

       
     
     
     
     
     
     
     
     
     
     
     
     
     

    1. ГИДРОЛОГИЧЕСКИЕ РАСЧЕТЫ

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

          Параметрами кривой обеспеченности являются: среднее  арифметическое, коэффициент вариации и коэффициент асимметрии, рассчитанные по рядам наблюдений.

          Норма стока – среднее значение величины стока за многолетний период такой  продолжительности, при увеличении которой полученное значение существенно  не меняется.

          При наличии данных наблюдений за n лет и расчётах методом моментов норма рассчитывается как среднее арифметическое:

    ,                                                       (2.1)

    где Qi – текущие значения расхода;

          n – число лет наблюдений.

          Коэффициент вариации Cv определяется по формуле

    ,                                (2.2)

    где  – модульный коэффициент.

          Коэффициент асимметрии рассчитывается в соответствии со СНиП 2.01.14–83 по формуле

                                             (2.3)

          или по формуле:                         .                                            (2.4)

         Средние квадратические ошибки Q, CV и CS при их определении методом моментов оцениваются по зависимостям:

    %;                                             (2.5)

    %;                                      (2.6)

    %.                                 (2.7)

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

         Определение и CV удобно выполнять в табличной форме (таблица 2.1).

         Рассчитаем  основные параметры кривой обеспеченности:

    163,94 м3/с.

    0,27.

    0,64.

    Определим средние квадратические ошибки методом моментов:

    5,4%<10%.

    4,09%<10%.

    16,36% >10%. 

    Таблица 2.1 - Расчеты нормы и коэффициента вариации. Среднегодовой сток. 

    годы Qi Qi-Qср (Qi-Qср)² K Ki - 1 (Ki - 1)² m p, % (Ki - 1)³ lg Ki K*lg Ki
    1936 288 124,06 15391,88 1,76 0,76 0,57 1 3,85 0,433 0,245 0,430
    1937 227 63,06 3977,07 1,38 0,38 0,15 2 7,69 0,057 0,141 0,196
    1938 224 60,06 3607,68 1,37 0,37 0,13 3 11,54 0,049 0,136 0,185
    1939 199 35,06 1229,48 1,21 0,21 0,05 4 15,38 0,010 0,084 0,102
    1940 194 30,06 903,84 1,18 0,18 0,03 5 19,23 0,006 0,073 0,087
    1941 189 25,06 628,20 1,15 0,15 0,02 6 23,08 0,004 0,062 0,071
    1942 187 23,06 531,95 1,14 0,14 0,02 7 26,92 0,003 0,057 0,065
    1943 183 19,06 363,44 1,12 0,12 0,01 8 30,77 0,002 0,048 0,053
    1944 183 19,06 363,44 1,12 0,12 0,01 9 34,62 0,002 0,048 0,053
    1945 176 12,06 145,54 1,07 0,07 0,01 10 38,46 0,000 0,031 0,033
    1946 172 8,06 65,03 1,05 0,05 0,00 11 42,31 0,000 0,021 0,022
    1948 168 4,06 16,52 1,02 0,02 0,00 12 46,15 0,000 0,011 0,011
    1949 165 1,06 1,13 1,01 0,01 0,00 13 50,00 0,000 0,003 0,003
    1950 160 -3,94 15,49 0,98 -0,02 0,00 14 53,85 0,000 -0,011 -0,010
    1952 155 -8,94 79,85 0,95 -0,05 0,00 15 57,69 0,000 -0,024 -0,023
    1953 151 -12,94 167,34 0,92 -0,08 0,01 16 61,54 0,000 -0,036 -0,033
    1954 140 -23,94 572,93 0,85 -0,15 0,02 17 65,38 -0,003 -0,069 -0,059
    1955 138 -25,94 672,68 0,84 -0,16 0,03 18 69,23 -0,004 -0,075 -0,063
    1956 132 -31,94 1019,91 0,81 -0,19 0,04 19 73,08 -0,007 -0,094 -0,076
    1957 127 -36,94 1364,27 0,77 -0,23 0,05 20 76,92 -0,011 -0,111 -0,086
    1958 125 -38,94 1516,01 0,76 -0,24 0,06 21 80,77 -0,013 -0,118 -0,090
    1959 112 -51,94 2697,35 0,68 -0,32 0,10 22 84,62 -0,032 -0,165 -0,113
    1960 110 -53,94 2909,09 0,67 -0,33 0,11 23 88,46 -0,036 -0,173 -0,116
    1961 105 -58,94 3473,45 0,64 -0,36 0,13 24 92,31 -0,046 -0,193 -0,124
    1962 88,4 -75,54 5705,69 0,54 -0,46 0,21 25 96,15 -0,098 -0,268 -0,145
      1. Подбор  теоретической кривой обеспеченности  

         Эмпирическая  обеспеченность в % каждого члена  статистического ряда вычисляется  по формуле Крицкого-Менкеля:

    %,                                              (2.8)

    где m – порядковый номер члена ранжированного (убывающего) ряда;

         n – число членов ряда.

         Поскольку при расчётах обеспеченности гидрологических  характеристик оперируют с относительно непродолжительными рядами наблюдений, всегда возникает необходимость  экстраполяции рассчитываемой величины за пределы наблюдений, осуществляемой с помощью кривой обеспеченности. Для сглаживания и экстраполяции  эмпирических кривых распределения, как  правило, применяют трехпараметрическое  распределение Крицкого-Менкеля при любом отношении Сs/Cv, распределение Пирсона III типа при Сs/Cv³2 и другие.

         При использовании распределения Пирсона  III типа по значениям Q, Сs и Cv, найденным по формулам, приведенным в п.2, находят нормированные отклонения теоретической кривой.

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

    KP%=1+ФP%·CV;                                                (2.9)

    QP%=KP%· ;                                                   (2.10)

    где KP% – модульный коэффициент обеспеченностью P%;

          QP% – расход обеспеченностью P%.

         Значения  KP% рассчитываются в диапазоне P=0,1…99%. Результаты расчётов сводятся в таблицу 2.2. Затем эти значения (KP%) выносятся на поле эмпирических точек (Рисунок А2, Рисунок А3) и производится построение теоретической кривой обеспеченности QP%=f(P). 
     

         Таблица 2.2 – Ординаты теоретических кривых обеспеченности для среднегодового стока

    P,% Фр,% Кр,% Q,%
    1 2 3 4
    0,1 3,66 2,0614 337,9377
    1 2,61 1,7569 288,0192
    3 2,04 1,5916 260,9205
    5 1,75 1,5075 247,1335
    10 1,32 1,3828 226,6907
    20 0,82 1,2378 202,92
    30 0,47 1,1363 186,2805
    40 0,19 1,0551 172,9689
    50 -0,07 0,9797 160,6081
    60 -0,31 0,9101 149,1982
    70 -0,57 0,8347 136,8374
    75 -0,71 0,7941 130,1816
    80 -0,85 0,7535 123,5258
    90 -1,23 0,6433 105,46
    95 -1,52 0,5592 91,67301
    99 -2,03 0,4113 67,42688
     

         Определение параметров кривой обеспеченности графоаналитическим методом Г. А. Алексеева выполняют  в следующей последовательности. Выносятся и осредняются эмпирические точки QP%=f(P%) плавной «эмпирической» кривой обеспеченности (Рисунок А 2), снимаются значения расходов Q5, Q50 и Q95 и рассчитывается коэффициент скошенности S по формуле:

    ,                                             (2.11)

    по значениям которого находится значение коэффициента асимметрии CS и чисел Фостера Ф5, Ф50 и Ф95.

         Среднее квадратическое отклонение σQ, среднее арифметическое и коэффициент вариации рассчитывают по зависимостям

    ;                                                 (2.12)

      =Q50 – σQ·Ф50;                                               (2.13)

    .                                                    (2.14)

    Q5=250 м3/с; Q50=163 м3/с; Q95=94 м3/с.

    0,11.

         Тогда, интерполируя, получим следующие  числа Фостера 5-, 50- и 95-про-центной обеспеченности:

    для S=0,11 

    Ф5 = 1,75; Ф50 = - 0,07; Ф95 = - 1,52;

    47,71%;

    =Q50 – σQ·Ф50=163–47,71·(–0,07)=166,34 м3/с;

    Определение водохозяйственных балансов